The power of single molecule real-time sequencing technology in the de novo assembly of a eukaryotic genome

Chromosome-level reference genome and resequencing of 322 accessions reveal evolution, genomic imprint and key agronomic traits in adzuki bean

Adzuki bean (Vigna angularis) is an important legume crop cultivated in over 30 countries worldwide. We developed a high-quality chromosome-level reference genome of adzuki bean cultivar Jingnong6 by combining PacBio Sequel long-read sequencing with short-read and Hi-C technologies. The assembled genome covers 97.8% of the adzuki bean genome with a contig N50 of approximately 16 Mb and a total of 32 738 protein-coding genes. We also generated a comprehensive genome variation map of adzuki bean by whole-genome resequencing (WGRS) of 322 diverse adzuki beans accessions including both wild and cultivated. Furthermore, we have conducted comparative genomics and a genome-wide association study (GWAS) on key agricultural traits to investigate the evolution and domestication. GWAS identified several candidate genes, including VaCycA3;1, VaHB15, VaANR1 and VaBm, that exhibited significant associations with domestication traits. Furthermore, we conducted functional analyses on the roles of VaANR1 and VaBm in regulating seed coat colour. We provided evidence for the highest genetic diversity of wild adzuki (Vigna angularis var. nipponensis) in China with the presence of the most original wild adzuki bean, and the occurrence of domestication process facilitating transition from wild to cultigen. The present study elucidates the genetic basis of adzuki bean domestication traits and provides crucial genomic resources to support future breeding efforts in adzuki bean.

A view of the pan-genome of domesticated Cowpea (Vigna unguiculata [L.] Walp.)

Cowpea, Vigna unguiculata L. Walp., is a diploid warm-season legume of critical importance as both food and fodder in sub-Saharan Africa. This species is also grown in Northern Africa, Europe, Latin America, North America, and East to Southeast Asia. To capture the genomic diversity of domesticates of this important legume, de novo genome assemblies were produced for representatives of six subpopulations of cultivated cowpea identified previously from genotyping of several hundred diverse accessions. In the most complete assembly (IT97K-499-35), 26,026 core and 4963 noncore genes were identified, with 35,436 pan genes when considering all seven accessions. GO terms associated with response to stress and defense response were highly enriched among the noncore genes, while core genes were enriched in terms related to transcription factor activity, and transport and metabolic processes. Over 5 million single nucleotide polymorphisms (SNPs) relative to each assembly and over 40 structural variants >1 Mb in size were identified by comparing genomes. Vu10 was the chromosome with the highest frequency of SNPs, and Vu04 had the most structural variants. Noncore genes harbor a larger proportion of potentially disruptive variants than core genes, including missense, stop gain, and frameshift mutations; this suggests that noncore genes substantially contribute to diversity within domesticated cowpea.

Translocations and inversions: major chromosomal rearrangements during Vigna (Leguminosae) evolution

KEY MESSAGE

Inversions and translocations are the major chromosomal rearrangements involved in Vigna subgenera evolution, being Vigna vexillata the most divergent species. Centromeric repositioning seems to be frequent within the genus. Oligonucleotide-based fluorescence in situ hybridization (Oligo-FISH) provides a powerful chromosome identification system for inferring plant chromosomal evolution. Aiming to understand macrosynteny, chromosomal diversity, and the evolution of bean species from five Vigna subgenera, we constructed cytogenetic maps for eight taxa using oligo-FISH-based chromosome identification. We used oligopainting probes from chromosomes 2 and 3 of Phaseolus vulgaris L. and two barcode probes designed from V. unguiculata (L.) Walp. genome. Additionally, we analyzed genomic blocks among the Ancestral Phaseoleae Karyotype (APK), two V. unguiculata subspecies (V. subg. Vigna), and V. angularis (Willd.) Ohwi & Ohashi (V. subg. Ceratotropis). We observed macrosynteny for chromosomes 2, 3, 4, 6, 7, 8, 9, and 10 in all investigated taxa except for V. vexillata (L.) A. Rich (V. subg. Plectrotropis), in which only chromosomes 4, 7, and 9 were unambiguously identified. Collinearity breaks involved with chromosomes 2 and 3 were revealed. We identified minor differences in the painting pattern among the subgenera, in addition to multiple intra- and interblock inversions and intrachromosomal translocations. Other rearrangements included a pericentric inversion in chromosome 4 (V. subg. Vigna), a reciprocal translocation between chromosomes 1 and 5 (V. subg. Ceratotropis), a potential deletion in chromosome 11 of V. radiata (L.) Wilczek, as well as multiple intrablock inversions and centromere repositioning via genomic blocks. Our study allowed the visualization of karyotypic patterns in each subgenus, revealing important information for understanding intrageneric karyotypic evolution, and suggesting V. vexillata as the most karyotypically divergent species.

Unique Salt-Tolerance-Related QTLs, Evolved in Vigna riukiuensis (Na+ Includer) and V. nakashimae (Na+ Excluder), Shed Light on the Development of Super-Salt-Tolerant Azuki Bean (V. angularis) Cultivars

Wild relatives of crops have the potential to improve food crops, especially in terms of improving abiotic stress tolerance. Two closely related wild species of the traditional East Asian legume crops, Azuki bean (Vigna angularis), V. riukiuensis "Tojinbaka" and V. nakashimae "Ukushima" were shown to have much higher levels of salt tolerance than azuki beans. To identify the genomic regions responsible for salt tolerance in "Tojinbaka" and "Ukushima", three interspecific hybrids were developed: (A) azuki bean cultivar "Kyoto Dainagon" × "Tojinbaka", (B) "Kyoto Dainagon" × "Ukushima" and (C) "Ukushima" × "Tojinbaka". Linkage maps were developed using SSR or restriction-site-associated DNA markers. There were three QTLs for "percentage of wilt leaves" in populations A, B and C, while populations A and B had three QTLs and population C had two QTLs for "days to wilt". In population C, four QTLs were detected for Na⁺ concentration in the primary leaf. Among the F₂ individuals in population C, 24% showed higher salt tolerance than both wild parents, suggesting that the salt tolerance of azuki beans can be further improved by combining the QTL alleles of the two wild relatives. The marker information would facilitate the transfer of salt tolerance alleles from "Tojinbaka" and "Ukushima" to azuki beans.

Mining legume germplasm for genetic gains: An Indian perspective

Legumes play a significant role in food and nutritional security and contribute to environmental sustainability. Although legumes are highly beneficial crops, it has not yet been possible to enhance their yield and production to a satisfactory level. Amid a rising population and low yield levels, per capita average legume consumption in India has fallen by 71% over the last 50 years, and this has led to protein-related malnutrition in a large segment of the Indian population, especially women and children. Several factors have hindered attempts to achieve yield enhancement in grain legumes, including biotic and abiotic pressures, a lack of good ideotypes, less amenability to mechanization, poorer responsiveness to fertilizer input, and a poor genetic base. Therefore, there is a need to mine the approximately 0.4 million ex situ collections of legumes that are being conserved in gene banks globally for identification of ideal donors for various traits. The Indian National Gene Bank conserves over 63,000 accessions of legumes belonging to 61 species. Recent initiatives have been undertaken in consortia mode with the aim of unlocking the genetic potential of ex situ collections and conducting large-scale germplasm characterization and evaluation analyses. We assume that large-scale phenotyping integrated with omics-based science will aid the identification of target traits and their use to enhance genetic gains. Additionally, in cases where the genetic base of major legumes is narrow, wild relatives have been evaluated, and these are being exploited through pre-breeding. Thus far, >200 accessions of various legumes have been registered as unique donors for various traits of interest.

Advancing Grain Legumes Domestication and Evolution Studies with Genomics

Grain legumes were domesticated in parallel with cereals in several regions of the world and formed the economic basis of early farming cultures. Since then, legumes have played a vital role in human and animal diets and in fostering agrobiodiversity. Increasing grain legume cultivation will be crucial to safeguard nutritional security and the resilience of agricultural ecosystems across the globe. A better understanding of the molecular underpinnings of domestication and crop evolution of grain legumes may be translated into practical approaches in modern breeding programs to stabilize yield, which is threatened by evolving pathogens and changing climates. During recent decades, domestication research in all crops has greatly benefited from the fast progress in genomic technologies. Yet still, many questions surrounding the domestication and diversification of legumes remain unanswered. In this review, we assess the potential of genomic approaches in grain legume research. We describe the centers of origin and the crucial domestication traits of grain legumes. In addition, we survey the effect of domestication on both above-ground and below-ground traits that have economic importance. Finally, we discuss open questions in grain legume domestication and diversification and outline how to bridge the gap between the preservation of historic crop diversity and their utilization in modern plant breeding.

Thirty Years of Mungbean Genome Research: Where Do We Stand and What Have We Learned?

Mungbean is a socioeconomically important legume crop in Asia that is currently in high demand by consumers and industries both as dried beans and in plant-based protein foods. Marker-assisted and genomics-assisted breeding are promising approaches to efficiently and rapidly develop new cultivars with improved yield, quality, and resistance to biotic and abiotic stresses. Although mungbean was at the forefront of research at the dawn of the plant genomics era 30 years ago, the crop is a "slow runner" in genome research due to limited genomic resources, especially DNA markers. Significant progress in mungbean genome research was achieved only within the last 10 years, notably after the release of the VC1973A draft reference genome constructed using next-generation sequencing technology, which enabled fast and efficient DNA marker development, gene mapping, and identification of candidate genes for complex traits. Resistance to biotic stresses has dominated mungbean genome research to date; however, research is on the rise. In this study, we provide an overview of the past progress and current status of mungbean genomics research. We also discuss and evaluate some research results to provide a better understanding of mungbean genomics.

Cantaloupe melon genome reveals 3D chromatin features and structural relationship with the ancestral cucurbitaceae karyotype

Cucumis melo displays a large diversity of horticultural groups with cantaloupe melon the most cultivated type. Using a combination of single-molecule sequencing, 10X Genomics link-reads, high-density optical and genetic maps, and chromosome conformation capture (Hi-C), we assembled a chromosome scale C. melo var. cantalupensis Charentais mono genome. Integration of RNA-seq, MeDip-seq, ChIP-seq, and Hi-C data revealed a widespread compartmentalization of the melon genome, segregating constitutive heterochromatin and euchromatin. Genome-wide comparative and evolutionary analysis between melon botanical groups identified Charentais mono genome increasingly more divergent from Harukei-3 (reticulatus), Payzawat (inodorus), and HS (ssp. agrestis) genomes. To assess the paleohistory of the Cucurbitaceae, we reconstructed the ancestral Cucurbitaceae karyotype and compared it to sequenced cucurbit genomes. In contrast to other species that experienced massive chromosome shuffling, melon has retained the ancestral genome structure. We provide comprehensive genomic resources and new insights in the diversity of melon horticultural groups and evolution of cucurbits.

•

We provide a chromosome scale C. melo var. cantalupensis Charentais mono genome

•

Epigenomic analysis revealed a widespread compartmentalization of the melon genome

•

We reconstructed the ancestral Cucurbitaceae karyotype

•

Melon has retained the ancestral Cucurbitaceae genome structure

BAC- and oligo-FISH mapping reveals chromosome evolution among Vigna angularis, V. unguiculata, and Phaseolus vulgaris

Cytogenomic resources have accelerated synteny and chromosome evolution studies in plant species, including legumes. Here, we established the first cytogenetic map of V. angularis (Va, subgenus Ceratotropis) and compared this new map with those of V. unguiculata (Vu, subgenus Vigna) and P. vulgaris (Pv) by BAC-FISH and oligopainting approaches. We mapped 19 Vu BACs and 35S rDNA probes to the 11 chromosome pairs of Va, Vu, and Pv. Vigna angularis shared a high degree of macrosynteny with Vu and Pv, with five conserved syntenic chromosomes. Additionally, we developed two oligo probes (Pv2 and Pv3) used to paint Vigna orthologous chromosomes. We confirmed two reciprocal translocations (chromosomes 2 and 3 and 1 and 8) that have occurred after the Vigna and Phaseolus divergence (~9.7 Mya). Besides, two inversions (2 and 4) and one translocation (1 and 5) have occurred after Vigna and Ceratotropis subgenera separation (~3.6 Mya). We also observed distinct oligopainting patterns for chromosomes 2 and 3 of Vigna species. Both Vigna species shared similar major rearrangements compared to Pv: one translocation (2 and 3) and one inversion (chromosome 3). The sequence synteny identified additional inversions and/or intrachromosomal translocations involving pericentromeric regions of both orthologous chromosomes. We propose chromosomes 2 and 3 as hotspots for chromosomal rearrangements and de novo centromere formation within and between Vigna and Phaseolus. Our BAC- and oligo-FISH mapping contributed to physically trace the chromosome evolution of Vigna and Phaseolus and its application in further studies of both genera.

Unequal contribution of two paralogous CENH3 variants in cowpea centromere function

In most diploids the centromere-specific histone H3 (CENH3), the assembly site of active centromeres, is encoded by a single copy gene. Persistance of two CENH3 paralogs in diploids species raises the possibility of subfunctionalization. Here we analysed both CENH3 genes of the  diploid dryland crop cowpea. Phylogenetic analysis suggests that gene duplication of CENH3 occurred independently during the speciation of Vigna unguiculata. Both functional CENH3 variants are transcribed, and the corresponding proteins are intermingled in subdomains of different types of centromere sequences in a tissue-specific manner together with the kinetochore protein CENPC. CENH3.2 is removed from the generative cell of mature pollen, while CENH3.1 persists. CRISPR/Cas9-based inactivation of CENH3.1 resulted in delayed vegetative growth and sterility, indicating that this variant is needed for plant development and reproduction. By contrast, CENH3.2 knockout individuals did not show obvious defects during vegetative and reproductive development. Hence, CENH3.2 of cowpea is likely at an early stage of pseudogenization and less likely undergoing subfunctionalization.

The effect of zebularine on the heat-activated retrotransposon ONSEN in Arabidopsis thaliana and Vigna angularis

The Ty1/copia-like retrotransposon ONSEN is conserved among Brassica species, as well as in beans, including adzuki bean (Vigna angularis (Willd.) Ohwi & Ohashi), which is one of the economically important crops in Japan. ONSEN has acquired a heat-responsive element that is recognized by plant heat stress defense factors, resulting in its transcription and the production of full-length extrachromosomal DNA under conditions with elevated temperatures. DNA methylation plays an important role in regulating the activation of this transposon in plants. Therefore, chemical inhibition of DNA methyltransferases has been utilized to study the effect of DNA methylation on transposon activation. To understand the effect of DNA methylation on ONSEN activation, Arabidopsis thaliana and adzuki bean seedlings were treated with zebularine, which is known to be an effective chemical demethylation agent. The results showed that ONSEN transcription levels were upregulated in zebularine-treated plants. Extrachromosomal DNA of ONSEN also accumulated in the treated plants.

Comparative genomics of muskmelon reveals a potential role for retrotransposons in the modification of gene expression

Melon exhibits substantial natural variation especially in fruit ripening physiology, including both climacteric (ethylene-producing) and non-climacteric types. However, genomic mechanisms underlying such variation are not yet fully understood. Here, we report an Oxford Nanopore-based high-grade genome reference in the semi-climacteric cultivar Harukei-3 (378 Mb + 33,829 protein-coding genes), with an update of tissue-wide RNA-seq atlas in the Melonet-DB database. Comparison between Harukei-3 and DHL92, the first published melon genome, enabled identification of 24,758 one-to-one orthologue gene pairs, whereas others were candidates of copy number variation or presence/absence polymorphisms (PAPs). Further comparison based on 10 melon genome assemblies identified genome-wide PAPs of 415 retrotransposon Gag-like sequences. Of these, 160 showed fruit ripening-inducible expression, with 59.4% of the neighboring genes showing similar expression patterns (r > 0.8). Our results suggest that retrotransposons contributed to the modification of gene expression during diversification of melon genomes, and may affect fruit ripening-inducible gene expression.

Same Locus for Non-shattering Seed Pod in Two Independently Domesticated Legumes, Vigna angularis and Vigna unguiculata

Loss of pod shattering is one of the most important domestication-related traits in legume crops. The non-shattering phenotypes have been achieved either by disturbed formation of abscission layer between the valves, or by loss of helical tension in sclerenchyma of endocarp, that split open the pods to disperse the seeds. During domestication, azuki bean (Vigna angularis) and yard-long bean (Vigna unguiculata cv-gr. Sesquipedalis) have reduced or lost the sclerenchyma and thus the shattering behavior of seed pods. Here we performed fine-mapping with backcrossed populations and narrowed the candidate genomic region down to 4 kbp in azuki bean and 13 kbp in yard-long bean. Among the genes located in these regions, we found MYB26 genes encoded truncated proteins in azuki bean, yard-long bean, and even cowpea. As such, our findings indicate that independent domestication on the two legumes has selected the same locus for the same traits. We also argue that MYB26 could be a target gene for improving shattering phenotype in other legumes, such as soybean.

Breaks of macrosynteny and collinearity among moth bean (Vigna aconitifolia), cowpea (V. unguiculata), and common bean (Phaseolus vulgaris)

Comparative cytogenetic mapping is a powerful approach to gain insights into genome organization of orphan crops, lacking a whole sequenced genome. To investigate the cytogenomic evolution of important Vigna and Phaseolus beans, we built a BAC-FISH (fluorescent in situ hybridization of bacterial artificial chromosome) map of Vigna aconitifolia (Vac, subgenus Ceratotropis), species with no sequenced genome, and compared with V. unguiculata (Vu, subgenus Vigna) and Phaseolus vulgaris (Pv) maps. Seventeen Pv BACs, eight Vu BACs, and 5S and 35S rDNA probes were hybridized in situ on the 11 Vac chromosome pairs. Five Vac chromosomes (Vac6, Vac7, Vac9, Vac10, and Vac11) showed conserved macrosynteny and collinearity between V. unguiculata and P. vulgaris. On the other hand, we observed collinearity breaks, identified by pericentric inversions involving Vac2 (Vu2), Vac4 (Vu4), and Vac3 (Pv3). We also detected macrosynteny breaks of translocation type involving chromosomes 1 and 8 of V. aconitifolia and P. vulgaris; 2 and 3 of V. aconitifolia and P. vulgaris; and 1 and 5 of V. aconitifolia and V. unguiculata. Considering our data and previous BAC-FISH studies, six chromosomes (1, 2, 3, 4, 5, and 8) are involved in major karyotype divergences between genera and five (1, 2, 3, 4, and 5) between Vigna subgenera, including mechanisms such as duplications, inversions, and translocations. Macrosynteny breaks between Vigna and Phaseolus suggest that the major chromosomal rearrangements have occurred within the Vigna clade. Our cytogenomic comparisons bring new light on the degree of shared macrosynteny and mechanisms of karyotype diversification during Vigna and Phaseolus evolution.

Next-Generation Sequencing: An Emerging Tool for Drug Designing

BACKGROUND

With the outbreak of high throughput next-generation sequencing (NGS), the biological research of drug discovery has been directed towards the oncology and infectious disease therapeutic areas, with extensive use in biopharmaceutical development and vaccine production.

METHOD

In this review, an effort was made to address the basic background of NGS technologies, potential applications of NGS in drug designing. Our purpose is also to provide a brief introduction of various Nextgeneration sequencing techniques.

DISCUSSIONS

The high-throughput methods execute Large-scale Unbiased Sequencing (LUS) which comprises of Massively Parallel Sequencing (MPS) or NGS technologies. The Next geneinvolved necessarily executes Largescale Unbiased Sequencing (LUS) which comprises of MPS or NGS technologies. These are related terms that describe a DNA sequencing technology which has revolutionized genomic research. Using NGS, an entire human genome can be sequenced within a single day.

CONCLUSION

Analysis of NGS data unravels important clues in the quest for the treatment of various lifethreatening diseases and other related scientific problems related to human welfare.

Comprehensive evaluation of non-hybrid genome assembly tools for third-generation PacBio long-read sequence data

Long reads obtained from third-generation sequencing platforms can help overcome the long-standing challenge of the de novo assembly of sequences for the genomic analysis of non-model eukaryotic organisms. Numerous long-read-aided de novo assemblies have been published recently, which exhibited superior quality of the assembled genomes in comparison with those achieved using earlier second-generation sequencing technologies. Evaluating assemblies is important in guiding the appropriate choice for specific research needs. In this study, we evaluated 10 long-read assemblers using a variety of metrics on Pacific Biosciences (PacBio) data sets from different taxonomic categories with considerable differences in genome size. The results allowed us to narrow down the list to a few assemblers that can be effectively applied to eukaryotic assembly projects. Moreover, we highlight how best to use limited genomic resources for effectively evaluating the genome assemblies of non-model organisms.

Whole genome sequencing data of Vigna nakashimae var. Ukushima and G418

Vigna nakashimae is one of the closely related species of Vigna angularis (Adzuki bean). Two strain of ‘Ukushima’ and ‘G418’ were identified as salt tolerance strains in Vigna nakashimae from gene bank collection. F₂ populations from an inter- or intra-specific cross between the sensitive and tolerant strains are useful for the detection of salt tolerance QTL in Vigna nakashimae. Although Vigna angularis reference genome is available and useful for genetic analysis by genotyping-by-sequencing/RADseq in closely related species, it is not enough for isolation of responsible genes. To reveal sequence variation in Vigna nakashimae “Ukushima” and “G418”, the whole genome sequencing was performed using Illumina HiSeq X Ten system (411,174,986 and 478,116,282 read). NGS data was deposited in the DNA Data Bank of Japan (DDBJ) under accession number DRA009307.

Updates on Legume Genome Sequencing

Legumes (Fabaceae) are agronomically and economically one of the most important crops. Because legumes serve as a source of food, feed, and industrial materials, many studies in the field of legume genomics, including genome sequencing, have been conducted over the last decade. Here, we update the progress in genome sequencing of legume crops, including soybean (Glycine max [L.] Merr.), mung bean (V. radiata var. radiata), adzuki bean (V. angularis var. angularis), common bean (Phaseolus vulgaris L.), pigeon pea (Cajanus cajan), chickpea (Cicer arietinum), and peanut (Arachis hypogaea). Since the publication of the first reference genome sequence of each species, many accessions have been resequenced to study genetic diversity, speciation, and polyploidization in the legume lineage.

Whole genome sequencing of Entamoeba nuttalli reveals mammalian host-related molecular signatures and a novel octapeptide-repeat surface protein

The enteric protozoa Entamoeba histolytica is the causative agent of amebiasis, which is one of the most common parasitic diseases in developed and developing countries. Entamoeba nuttalli is the genetically closest species to E. histolytica in current phylogenetic analyses of Entamoeba species, and is prevalent in wild macaques. Therefore, E. nuttalli may be a key organism in which to investigate molecules required for infection of human or non-human primates. To explore the molecular signatures of host-parasite interactions, we conducted de novo assembly of the E. nuttalli genome, utilizing self-correction of PacBio long reads and polishing corrected reads using Illumina short reads, followed by comparative genomic analysis with two other mammalian and a reptilian Entamoeba species. The final draft assembly of E. nuttalli included 395 contigs with a total length of approximately 23 Mb, and 9,647 predicted genes, of which 6,940 were conserved with E. histolytica. In addition, we found an E. histolytica-specific repeat known as ERE2 in the E. nuttalli genome. GO-term enrichment analysis of mammalian host-related molecules indicated diversification of transmembrane proteins, including AIG1 family and BspA-like proteins that may be involved in the host-parasite interaction. Furthermore, we identified an E. nuttalli-specific protein that contained 42 repeats of an octapeptide ([G,E]KPTDTPS). This protein was shown to be localized on the cell surface using immunofluorescence. Since many repeat-containing proteins in parasites play important roles in interactions with host cells, this unique octapeptide repeat-containing protein may be involved in colonization of E. nuttalli in the intestine of macaques. Overall, our draft assembly provides a valuable resource for studying Entamoeba evolution and host-parasite selection.

The genome of cowpea (Vigna unguiculata [L.] Walp.)

Cowpea (Vigna unguiculata [L.] Walp.) is a major crop for worldwide food and nutritional security, especially in sub-Saharan Africa, that is resilient to hot and drought-prone environments. An assembly of the single-haplotype inbred genome of cowpea IT97K-499-35 was developed by exploiting the synergies between single-molecule real-time sequencing, optical and genetic mapping, and an assembly reconciliation algorithm. A total of 519 Mb is included in the assembled sequences. Nearly half of the assembled sequence is composed of repetitive elements, which are enriched within recombination-poor pericentromeric regions. A comparative analysis of these elements suggests that genome size differences between Vigna species are mainly attributable to changes in the amount of Gypsy retrotransposons. Conversely, genes are more abundant in more distal, high-recombination regions of the chromosomes; there appears to be more duplication of genes within the NBS-LRR and the SAUR-like auxin superfamilies compared with other warm-season legumes that have been sequenced. A surprising outcome is the identification of an inversion of 4.2 Mb among landraces and cultivars, which includes a gene that has been associated in other plants with interactions with the parasitic weed Striga gesnerioides. The genome sequence facilitated the identification of a putative syntelog for multiple organ gigantism in legumes. A revised numbering system has been adopted for cowpea chromosomes based on synteny with common bean (Phaseolus vulgaris). An estimate of nuclear genome size of 640.6 Mbp based on cytometry is presented.

Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants

BACKGROUND

The sharp increase of plant genome and transcriptome data provide valuable resources to investigate evolutionary consequences of gene duplication in a range of taxa, and unravel common principles underlying duplicate gene retention.

RESULTS

We survey 141 sequenced plant genomes to elucidate consequences of gene and genome duplication, processes central to the evolution of biodiversity. We develop a pipeline named DupGen_finder to identify different modes of gene duplication in plants. Genes derived from whole-genome, tandem, proximal, transposed, or dispersed duplication differ in abundance, selection pressure, expression divergence, and gene conversion rate among genomes. The number of WGD-derived duplicate genes decreases exponentially with increasing age of duplication events-transposed duplication- and dispersed duplication-derived genes declined in parallel. In contrast, the frequency of tandem and proximal duplications showed no significant decrease over time, providing a continuous supply of variants available for adaptation to continuously changing environments. Moreover, tandem and proximal duplicates experienced stronger selective pressure than genes formed by other modes and evolved toward biased functional roles involved in plant self-defense. The rate of gene conversion among WGD-derived gene pairs declined over time, peaking shortly after polyploidization. To provide a platform for accessing duplicated gene pairs in different plants, we constructed the Plant Duplicate Gene Database.

CONCLUSIONS

We identify a comprehensive landscape of different modes of gene duplication across the plant kingdom by comparing 141 genomes, which provides a solid foundation for further investigation of the dynamic evolution of duplicate genes.

Construction of genetic linkage map and genome dissection of domestication-related traits of moth bean (Vigna aconitifolia), a legume crop of arid areas

The moth bean (Vigna aconitifolia), possibly the most primitive crop of the genus Vigna, is a highly drought- and heat-resistant legume grown in arid areas. Moth bean domestication involved phenotypic changes, including reduction of seed dormancy and pod shattering, increased organ size, and earlier flowering and maturity. However, the genetics of the domestication process in moth bean is not known. In this study, we constructed a genetic linkage map for moth bean and used the map to identify quantitative trait loci (QTL) for domestication-related traits of an F₂ population of 188 individuals produced from a cross of wild moth bean (TN67) and cultivated moth bean (ICPMO056). The genetic linkage map comprised 11 linkage groups (LG) of 172 simple sequence repeat markers and spanned a total length of 1016.8 centiMorgan (cM), with an average marker distance of 7.34 cM. A comparative genome analysis showed high genome synteny between moth bean and mungbean (Vigna radiata), adzuki bean (Vigna angularis), rice bean (Vigna umbellata), and yardlong bean (Vigna unguiculata). In total, 50 QTLs and 3 genes associated with 20 domestication-related traits were identified. Most of the QTLs belonged to five LGs (1, 2, 4, 7, and 10). Key traits related to domestication such as seed dormancy and pod shattering were controlled by large-effect QTLs (PVE > 20%) with one or two minor QTLs, whereas all other traits were controlled by one-seven minor QTLs, apart from seed weight, which was controlled by one major and seven minor QTLs. These results suggest that a small number of mutations with large phenotypic effects have contributed to the domestication of the moth bean. Comparative analysis of QTLs with related Vigna crops revealed that there are several domestication-related large-effect QTLs that had not been used in moth bean domestication. This study provides a basic genetic map and identified genome regions associated with domestication-related traits, which will be useful for the genetic improvement of the moth bean and related Vigna species.

Domesticating Vigna Stipulacea: A Potential Legume Crop With Broad Resistance to Biotic Stresses

Though crossing wild relatives to modern cultivars is a usual means to introduce alleles of stress tolerance, an alternative is de novo domesticating wild species that are already tolerant to various kinds of stresses. As a test case, we chose Vigna stipulacea Kuntze, which has fast growth, short vegetative stage, and broad resistance to pests and diseases. We developed an ethyl methanesulfonate-mutagenized population and obtained three mutants with reduced seed dormancy and one with reduced pod shattering. We crossed one of the mutants of less seed dormancy to the wild type and confirmed that the phenotype was inherited in a Mendelian manner. De novo assembly of V. stipulacea genome, and the following resequencing of the F2 progenies successfully identified a Single Nucleotide Polymorphism (SNP) associated with seed dormancy. By crossing and pyramiding the mutant phenotypes, we will be able to turn V. stipulacea into a crop which is yet primitive but can be cultivated without pesticides.

Genetic Dissection of Azuki Bean Weevil (Callosobruchus chinensis L.) Resistance in Moth Bean (Vigna aconitifolia [Jaqc.] Maréchal)

The azuki bean weevil (Callosobruchus chinensis L.) is an insect pest responsible for serious postharvest seed loss in leguminous crops. In this study, we performed quantitative trait locus (QTL) mapping of seed resistance to C. chinensis in moth bean (Vigna aconitifolia [Jaqc.] Maréchal). An F₂ population of 188 plants developed by crossing resistant accession 'TN67' (wild type from India; male parent) and susceptible accession 'IPCMO056' (cultivated type from India; female parent) was used for mapping. Seeds of the F₂ population from 2014 and F2:₃ populations from 2016 and 2017 were bioassayed with C. chinensis, and the percentage of damaged seeds and progress of infestation severity were measured. Segregation analysis suggested that C. chinensis resistance in TN176 is controlled by a single dominant gene, designated as Rcc. QTL analysis revealed one principal and one modifying QTL for the resistance, named qVacBrc2.1 and qVacBrc5.1, respectively. qVacBrc2.1 was located on linkage group 2 between simple sequence repeat markers CEDG261 and DMB-SSR160 and accounted for 50.41% to 64.23% of resistance-related traits, depending on the trait and population. Comparative genomic analysis suggested that qVacBrc2.1 is the same as QTL Brc2.1 conferring C. chinensis resistance in wild azuki bean (V. nepalensis Tateishi and Maxted). Markers CEDG261 and DMB-SSR160 should be useful for marker-assisted selection for C. chinensis resistance in moth bean.

Chromosome-scale assemblies of plant genomes using nanopore long reads and optical maps

Plant genomes are often characterized by a high level of repetitiveness and polyploid nature. Consequently, creating genome assemblies for plant genomes is challenging. The introduction of short-read technologies 10 years ago substantially increased the number of available plant genomes. Generally, these assemblies are incomplete and fragmented, and only a few are at the chromosome scale. Recently, Pacific Biosciences and Oxford Nanopore sequencing technologies were commercialized that can sequence long DNA fragments (kilobases to megabase) and, using efficient algorithms, provide high-quality assemblies in terms of contiguity and completeness of repetitive regions^1-4. However, even though genome assemblies based on long reads exhibit high contig N50s (>1 Mb), these methods are still insufficient to decipher genome organization at the chromosome level. Here, we describe a strategy based on long reads (MinION or PromethION sequencers) and optical maps (Saphyr system) that can produce chromosome-level assemblies and demonstrate applicability by generating high-quality genome sequences for two new dicotyledon morphotypes, Brassica rapa Z1 (yellow sarson) and Brassica oleracea HDEM (broccoli), and one new monocotyledon, Musa schizocarpa (banana). All three assemblies show contig N50s of >5 Mb and contain scaffolds that represent entire chromosomes or chromosome arms.

Assembling the genome of the African wild rice Oryza longistaminata by exploiting synteny in closely related Oryza species

The African wild rice species Oryza longistaminata has several beneficial traits compared to cultivated rice species, such as resistance to biotic stresses, clonal propagation via rhizomes, and increased biomass production. To facilitate breeding efforts and functional genomics studies, we de-novo assembled a high-quality, haploid-phased genome. Here, we present our assembly, with a total length of 351 Mb, of which 92.2% was anchored onto 12 chromosomes. We detected 34,389 genes and 38.1% of the genome consisted of repetitive content. We validated our assembly by a comparative linkage analysis and by examining well-characterized gene families. This genome assembly will be a useful resource to exploit beneficial alleles found in O. longistaminata. Our results also show that it is possible to generate a high-quality, functionally complete rice genome assembly from moderate SMRT read coverage by exploiting synteny in a closely related Oryza species.

Stefan Reuscher et al. assembled the genome of an African wild rice species to facilitate breeding efforts and functional genomic studies. They used SMRT sequencing, chromosomal synteny between rice species, and a linkage map to assemble the 351 Mb genome into 12 chromosomes.

Diversity of Drought Tolerance in the Genus Vigna

Wild relatives of crop plants are thought as reservoir of prominent genetic resources for abiotic stress tolerance. However, insufficient information on genetic variation and phenotypic traits restricts their use for crop breeding. This study focused on wild species of genus Vigna (family Fabaceae) originated from highly humid to arid regions. To clarify the diversity of drought tolerance during the vegetative stage, 69 accessions, including 15 domesticated, and 54 wild accessions, were evaluated under two drought conditions of non-terminal and terminal stresses. In the non-terminal drought condition, the plants were grown in pipes of different heights where surface soil water content decreased faster in pipes with greater height. Relative shoot biomass was used for tolerance evaluation and we identified 19 drought tolerant accessions. Almost of them were wild accessions showing higher relative shoot biomass than that in the domesticated accessions. Domesticated species were mostly classified as drought susceptible but could be improved using tolerant conspecific wild ancestors with cross-compatibility. The tolerance was related with higher plant water status presumably due to small water consumption. However, the variation of drought tolerance could not be explained by simple tolerance factor alone, and other tolerance mechanisms such as deep rooting and increasing in root biomass were found in the tolerant accessions. In the terminal drought condition, the plants were grown in small pots, and the watering was stopped to expose them extreme and rapid soil water scarcity. The tolerance was evaluated as the number of days until wilting. However, the accessions found to be tolerant in the pot experiment were not the same as those in the pipe experiment. In this condition, plant water status was not related with the length of days to wilting. This indicates that different mechanisms are necessary for adaptation to each of the non-terminal and terminal drought conditions. Many accessions were tolerant to one of the conditions, although we identified that some accessions showed tolerance in both experiments. The great diversity in drought tolerance in the genus Vigna might serve to both improve crop drought tolerance and understand the mechanisms of adaptation in drought-prone environments.

Characterization of a heat-activated retrotransposon in Vigna angularis

In plants, several transposable elements are conserved across species. We found a homolog of ONSEN, which is a heat-activated retrotransposon originally isolated from Arabidopsis thaliana, in Vigna. The ONSEN-like elements (VaONS) were detected in all the analyzed Japanese accessions of Vigna angularis (adzuki bean) by Southern blot analysis. However, VaONS sequences were observed to be polymorphic in the different accessions. Interestingly, extrachromosomal DNA (ecDNA) was detected in some accessions of adzuki bean, indicating the conserved heat-activation of VaONS. Furthermore, we successfully induced retrotransposition of VaONS in adzuki plant regenerated through callus. Findings of our study should provide a new tool for molecular breeding of adzuki bean.

Advances in Sequencing and Resequencing in Crop Plants

DNA sequencing technologies have changed the face of biological research over the last 20 years. From reference genomes to population level resequencing studies, these technologies have made significant contributions to our understanding of plant biology and evolution. As the technologies have increased in power, the breadth and complexity of the questions that can be asked has increased. Along with this, the challenges of managing unprecedented quantities of sequence data are mounting. This chapter describes a few aspects of the journey so far and looks forward to what may lie ahead. Graphical Abstract.

Molecular cytogenetic characterization and comparison of the two cultivated Canavalia species (Fabaceae)

The two cultivated Canavalia (Adanson, 1763) species, Canavalia gladiata (N. J. von Jacquin, 1788) A. P. de Candolle, 1825 and Canavalia ensiformis (Linnaeus, 1753) A. P. de Candolle, 1825 are closely related based on morphological and molecular phylogenetic data. However, the similarities and differences in genome organization between them have not been evaluated at molecular cytogenetic level. Here, detailed karyotypes of both species were constructed using combined PI and DAPI (CPD) staining, rDNA-FISH and self-genomic in situ hybridization (sGISH). For further comparison, comparative genomic in situ hybridization (cGISH) and sequence analysis of 5S rDNA were applied. Their chromosomes were accurately identified by sGISH and rDNA-FISH signals. Both species had the karyotype formula 2n = 22 = 18m + 4m-SAT, but the karyotype of C. ensiformis was shorter and more asymmetric than that of C. gladiata. They displayed similar CPD bands at all 45S rDNA sites and centromeres. C. gladiata had ten centromeric 5S rDNA loci and two SC (secondary constriction)-associated 45S rDNA loci. C. ensiformis had nine centromeric and one interstitial 5S loci, two SC-associated and one proximal 45S loci. Their sGISH signal patterns displayed both basic similarities and distinct differences. Reciprocal cGISH generated prominent signals in all pericentromeric regions and 45S sites. There was lower level of sequence identity of the non-transcribed spacer between their 5S rDNA repeats. These data confirmed the evolutionary closeness between C. gladiata and C. ensiformis and demonstrated obvious differentiation between their genomes, and supported the opinion that C. ensiformis is more advanced in evolution than C. gladiata.

Gapless genome assembly of Colletotrichum higginsianum reveals chromosome structure and association of transposable elements with secondary metabolite gene clusters

Background

The ascomycete fungus Colletotrichum higginsianum causes anthracnose disease of brassica crops and the model plant Arabidopsis thaliana. Previous versions of the genome sequence were highly fragmented, causing errors in the prediction of protein-coding genes and preventing the analysis of repetitive sequences and genome architecture.

Results

Here, we re-sequenced the genome using single-molecule real-time (SMRT) sequencing technology and, in combination with optical map data, this provided a gapless assembly of all twelve chromosomes except for the ribosomal DNA repeat cluster on chromosome 7. The more accurate gene annotation made possible by this new assembly revealed a large repertoire of secondary metabolism (SM) key genes (89) and putative biosynthetic pathways (77 SM gene clusters). The two mini-chromosomes differed from the ten core chromosomes in being repeat- and AT-rich and gene-poor but were significantly enriched with genes encoding putative secreted effector proteins. Transposable elements (TEs) were found to occupy 7% of the genome by length. Certain TE families showed a statistically significant association with effector genes and SM cluster genes and were transcriptionally active at particular stages of fungal development. All 24 subtelomeres were found to contain one of three highly-conserved repeat elements which, by providing sites for homologous recombination, were probably instrumental in four segmental duplications.

Conclusion

The gapless genome of C. higginsianum provides access to repeat-rich regions that were previously poorly assembled, notably the mini-chromosomes and subtelomeres, and allowed prediction of the complete SM gene repertoire. It also provides insights into the potential role of TEs in gene and genome evolution and host adaptation in this asexual pathogen.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4083-x) contains supplementary material, which is available to authorized users.

Advantages of genome sequencing by long-read sequencer using SMRT technology in medical area

PacBio RS II is the first commercialized third-generation DNA sequencer able to sequence a single molecule DNA in real-time without amplification. PacBio RS II's sequencing technology is novel and unique, enabling the direct observation of DNA synthesis by DNA polymerase. PacBio RS II confers four major advantages compared to other sequencing technologies: long read lengths, high consensus accuracy, a low degree of bias, and simultaneous capability of epigenetic characterization. These advantages surmount the obstacle of sequencing genomic regions such as high/low G+C, tandem repeat, and interspersed repeat regions. Moreover, PacBio RS II is ideal for whole genome sequencing, targeted sequencing, complex population analysis, RNA sequencing, and epigenetics characterization. With PacBio RS II, we have sequenced and analyzed the genomes of many species, from viruses to humans. Herein, we summarize and review some of our key genome sequencing projects, including full-length viral sequencing, complete bacterial genome and almost-complete plant genome assemblies, and long amplicon sequencing of a disease-associated gene region. We believe that PacBio RS II is not only an effective tool for use in the basic biological sciences but also in the medical/clinical setting.

Organelle_PBA, a pipeline for assembling chloroplast and mitochondrial genomes from PacBio DNA sequencing data

Background

The development of long-read sequencing technologies, such as single-molecule real-time (SMRT) sequencing by PacBio, has produced a revolution in the sequencing of small genomes. Sequencing organelle genomes using PacBio long-read data is a cost effective, straightforward approach. Nevertheless, the availability of simple-to-use software to perform the assembly from raw reads is limited at present.

Results

We present Organelle-PBA, a Perl program designed specifically for the assembly of chloroplast and mitochondrial genomes. For chloroplast genomes, the program selects the chloroplast reads from a whole genome sequencing pool, maps the reads to a reference sequence from a closely related species, and then performs read correction and de novo assembly using Sprai. Organelle-PBA completes the assembly process with the additional step of scaffolding by SSPACE-LongRead. The program then detects the chloroplast inverted repeats and reassembles and re-orients the assembly based on the organelle origin of the reference. We have evaluated the performance of the software using PacBio reads from different species, read coverage, and reference genomes. Finally, we present the assembly of two novel chloroplast genomes from the species Picea glauca (Pinaceae) and Sinningia speciosa (Gesneriaceae).

Conclusion

Organelle-PBA is an easy-to-use Perl-based software pipeline that was written specifically to assemble mitochondrial and chloroplast genomes from whole genome PacBio reads. The program is available at https://github.com/aubombarely/Organelle_PBA.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3412-9) contains supplementary material, which is available to authorized users.

Genome sequence and analysis of the Japanese morning glory Ipomoea nil

Ipomoea is the largest genus in the family Convolvulaceae. Ipomoea nil (Japanese morning glory) has been utilized as a model plant to study the genetic basis of floricultural traits, with over 1,500 mutant lines. In the present study, we have utilized second- and third-generation-sequencing platforms, and have reported a draft genome of I. nil with a scaffold N50 of 2.88 Mb (contig N50 of 1.87 Mb), covering 98% of the 750 Mb genome. Scaffolds covering 91.42% of the assembly are anchored to 15 pseudo-chromosomes. The draft genome has enabled the identification and cataloguing of the Tpn1 family transposons, known as the major mutagen of I. nil, and analysing the dwarf gene, CONTRACTED, located on the genetic map published in 1956. Comparative genomics has suggested that a whole genome duplication in Convolvulaceae, distinct from the recent Solanaceae event, has occurred after the divergence of the two sister families.

Diversity and Evolution of Salt Tolerance in the Genus Vigna

Breeding salt tolerant plants is difficult without utilizing a diversity of wild crop relatives. Since the genus Vigna (family Fabaceae) is comprised of many wild relatives adapted to various environmental conditions, we evaluated the salt tolerance of 69 accessions of this genus, including that of wild and domesticated accessions originating from Asia, Africa, Oceania, and South America. We grew plants under 50 mM and 200 mM NaCl for two weeks and then measured the biomass, relative quantum yield of photosystem II, leaf Na+ concentrations, and leaf K+ concentrations. The accessions were clustered into four groups: the most tolerant, tolerant, moderately susceptible, and susceptible. From the most tolerant group, we selected six accessions, all of which were wild accessions adapted to coastal environments, as promising sources of salt tolerance because of their consistently high relative shoot biomass and relative quantum yield. Interestingly, variations in leaf Na+ concentration were observed between the accessions in the most tolerant group, suggesting different mechanisms were responsible for their salt tolerance. Phylogenetic analysis with nuclear DNA sequences revealed that salt tolerance had evolved independently at least four times in the genus Vigna, within a relatively short period. The findings suggested that simple genetic changes in a few genes might have greatly affected salt tolerances. The elucidation of genetic mechanisms of salt tolerances in the selected accessions may contribute to improving the poor salt tolerance in legume crops.

Towards integration of population and comparative genomics in forest trees

Contents 338 I. 338 II. 339 III. 340 IV. 342 343 References 343 SUMMARY: The past decade saw the initiation of an ongoing revolution in sequencing technologies that is transforming all fields of biology. This has been driven by the advent and widespread availability of high-throughput, massively parallel short-read sequencing (MPS) platforms. These technologies have enabled previously unimaginable studies, including draft assemblies of the massive genomes of coniferous species and population-scale resequencing. Transcriptomics studies have likewise been transformed, with RNA-sequencing enabling studies in nonmodel organisms, the discovery of previously unannotated genes (novel transcripts), entirely new classes of RNAs and previously unknown regulatory mechanisms. Here we touch upon current developments in the areas of genome assembly, comparative regulomics and population genetics as they relate to studies of forest tree species.

Undiscovered Physiology of Transcript and Protein Networks

The past two decades have witnessed a rapid evolution in our ability to measure RNA and protein from biological systems. As a result, new principles have arisen regarding how information is processed in cells, how decisions are made, and the role of networks in biology. This essay examines this technological evolution, reviewing (and critiquing) the conceptual framework that has emerged to explain how RNA and protein networks control cellular function. We identify how future investigations into transcriptomes, proteomes, and other cellular networks will enable development of more robust, quantitative models of cellular behavior whilst also providing new avenues to use knowledge of biological networks to improve human health. © 2016 American Physiological Society. Compr Physiol 6:1851-1872, 2016.

Draft genome sequence of an inbred line of Chenopodium quinoa, an allotetraploid crop with great environmental adaptability and outstanding nutritional properties

Chenopodium quinoa Willd. (quinoa) originated from the Andean region of South America, and is a pseudocereal crop of the Amaranthaceae family. Quinoa is emerging as an important crop with the potential to contribute to food security worldwide and is considered to be an optimal food source for astronauts, due to its outstanding nutritional profile and ability to tolerate stressful environments. Furthermore, plant pathologists use quinoa as a representative diagnostic host to identify virus species. However, molecular analysis of quinoa is limited by its genetic heterogeneity due to outcrossing and its genome complexity derived from allotetraploidy. To overcome these obstacles, we established the inbred and standard quinoa accession Kd that enables rigorous molecular analysis, and presented the draft genome sequence of Kd, using an optimized combination of high-throughput next generation sequencing on the Illumina Hiseq 2500 and PacBio RS II sequencers. The de novo genome assembly contained 25 k scaffolds consisting of 1 Gbp with N50 length of 86 kbp. Based on these data, we constructed the free-access Quinoa Genome DataBase (QGDB). Thus, these findings provide insights into the mechanisms underlying agronomically important traits of quinoa and the effect of allotetraploidy on genome evolution.

The Vigna Genome Server, 'VigGS': A Genomic Knowledge Base of the Genus Vigna Based on High-Quality, Annotated Genome Sequence of the Azuki Bean, Vigna angularis (Willd.) Ohwi & Ohashi

The genus Vigna includes legume crops such as cowpea, mungbean and azuki bean, as well as >100 wild species. A number of the wild species are highly tolerant to severe environmental conditions including high-salinity, acid or alkaline soil; drought; flooding; and pests and diseases. These features of the genus Vigna make it a good target for investigation of genetic diversity in adaptation to stressful environments; however, a lack of genomic information has hindered such research in this genus. Here, we present a genome database of the genus Vigna, Vigna Genome Server ('VigGS', http://viggs.dna.affrc.go.jp), based on the recently sequenced azuki bean genome, which incorporates annotated exon-intron structures, along with evidence for transcripts and proteins, visualized in GBrowse. VigGS also facilitates user construction of multiple alignments between azuki bean genes and those of six related dicot species. In addition, the database displays sequence polymorphisms between azuki bean and its wild relatives and enables users to design primer sequences targeting any variant site. VigGS offers a simple keyword search in addition to sequence similarity searches using BLAST and BLAT. To incorporate up to date genomic information, VigGS automatically receives newly deposited mRNA sequences of pre-set species from the public database once a week. Users can refer to not only gene structures mapped on the azuki bean genome on GBrowse but also relevant literature of the genes. VigGS will contribute to genomic research into plant biotic and abiotic stresses and to the future development of new stress-tolerant crops.