Differential genome evolution and speciation of Coix lacryma-jobi L. and Coix aquatica Roxb. hybrid guangxi revealed by repetitive sequence analysis and fine karyotyping

Omics-driven utilization of wild relatives for empowering pre-breeding in pearl millet

MAIN CONCLUSION

Pearl millet wild relatives harbour novel alleles which could be utilized to broaden genetic base of cultivated species. Genomics-informed pre-breeding is needed to speed up introgression from wild to cultivated gene pool in pearl millet. Rising episodes of intense biotic and abiotic stresses challenge pearl millet production globally. Wild relatives provide a wide spectrum of novel alleles which could address challenges posed by climate change. Pre-breeding holds potential to introgress novel diversity in genetically narrow cultivated Pennisetum glaucum from diverse gene pool. Practical utilization of gene pool diversity remained elusive due to genetic intricacies. Harnessing promising traits from wild pennisetum is limited by lack of information on underlying candidate genes/QTLs. Next-Generation Omics provide vast scope to speed up pre-breeding in pearl millet. Genomic resources generated out of draft genome sequence and improved genome assemblies can be employed to utilize gene bank accessions effectively. The article highlights genetic richness in pearl millet and its utilization with a focus on harnessing next-generation Omics to empower pre-breeding.

Construction of a high-density genetic map based on large-scale marker development in Coix lacryma-jobi L. using specific-locus amplified fragment sequencing (slaf-seq)

Coix lacryma-jobi L. is one of the most economically and medicinally important corns. This study constructed a high-density genetic linkage map of C. lacryma-jobi based on a cross between the parents 'Qianyi No. 2' × 'Wenyi No. 2' and their F2 progeny through high-throughput sequencing and the construction of a specific-locus amplified fragment (SLAF) library. After pre-processing, 325.49 GB of raw data containing 1628 M reads were obtained. A total of 22,944 high-quality SLAFs were identified, among which 3952 SLAFs and 3646 polymorphic markers met the requirements for the construction of a genetic linkage map. The integrated map contained 3605 high-quality SLAFs, which were grouped into ten genetic linkage groups. The total length of the map was 1620.39 cM, with an average distance of 0.45 cM and an average of 360.5 markers per linkage group. This report presents the first high-density genetic map of C. lacryma-jobi. This map was constructed using an F2 population and SLAF-seq approach, which allows the development of a large number of polymorphic markers in a short period. These results provide a platform for precise gene/quantitative trait locus (QTL) mapping, map-based gene separation, and molecular breeding in C. lacryma-jobi. They also help identify a target gene for tracking, splitting quantitative traits, and estimating the phenotypic effects of each QTL for QTL mapping. They are of great significance for improving the efficiency of discovering and utilizing excellent gene resources of C. lacryma-jobi.

Cyto-molecular characterization of rDNA and chromatin composition in the NOR-associated satellite in Chestnut (Castanea spp.)

The American chestnut (Castanea dentata, 2n = 2x = 24), once known as the "King of the Appalachian Forest", was decimated by chestnut blight during the first half of the twentieth century by an invasive fungus (Cryphonectria parasitica). The Chinese chestnut (C. mollissima, 2n = 2x = 24), in contrast to American chestnut, is resistant to this blight. Efforts are being made to transfer this resistance to American chestnut through backcross breeding and genetic engineering. Both chestnut genomes have been genetically mapped and recently sequenced to facilitate gene discovery efforts aimed at assisting molecular breeding and genetic engineering. To complement and extend this genomic work, we analyzed the distribution and organization of their ribosomal DNAs (35S and 5S rDNA), and the chromatin composition of the nucleolus organizing region (NOR)-associated satellites. Using fluorescent in situ hybridization (FISH), we have identified two 35S (one major and one minor) and one 5S rDNA sites. The major 35S rDNA sites are terminal and sub-terminal in American and Chinese chestnuts, respectively, originating at the end of the short arm of the chromosome, extending through the secondary constriction and into the satellites. An additional 5S locus was identified in certain Chinese chestnut accessions, and it was linked distally to the major 35S site. The NOR-associated satellite in Chinese chestnut was found to comprise a proximal region packed with 35S rDNA and a distinct distal heterochromatic region. In contrast, the American chestnut satellite was relatively small and devoid of the distal heterochromatic region.

The potentialities of omics resources for millet improvement

Millets are nutrient-rich (nutri-rich) cereals with climate resilience attributes. However, its full productive potential is not realized due to the lack of a focused yield improvement approach, as evidenced by the available literature. Also, the lack of well-characterized genomic resources significantly limits millet improvement. But the recent availability of genomic data and advancement in omics tools has shown its enormous potential to enhance the efficiency and precision faced by conventional breeding in millet improvement. The development of high throughput genotyping platforms based on next-generation sequencing (NGS) has provided a low-cost method for genomic information, specifically for neglected nutri-rich cereals with the availability of a limited number of reference genome sequences. NGS has created new avenues for millet biotechnological interventions such as mutation-based study, GWAS, GS, and other omics technologies. The simultaneous discovery of high-throughput markers and multiplexed genotyping platform has aggressively aided marker-assisted breeding for millet improvement. Therefore, omics technology offers excellent opportunities to explore and combine useful variations for targeted traits that could impart high nutritional value to high-yielding cultivars under changing climatic conditions. In millet improvement, an in-depth account of NGS, integrating genomics data with different biotechnology tools, is reviewed in this context.

Genome-wide identification and expression analysis of VQ gene family under abiotic stress in Coix lacryma-jobi L

BACKGROUND

Valine-glutamine (VQ) proteins are non-specific plant proteins that have a highly conserved motif: FxxhVQxhTG. These proteins are involved in the development of various plant organs such as seeds, hypocotyls, flowers, leaves and also play a role in response to salt, drought and cold stresses. Despite their importance, there is limited information available on the evolutionary and structural characteristics of VQ family genes in Coix lacryma-jobi.

RESULTS

In this study, a total of 31 VQ genes were identified from the coix genome and classified into seven subgroups (I-VII) based on phylogenetic analysis. These genes were found to be unevenly distributed on 10 chromosomes. Gene structure analysis revealed that these genes had a similar type of structure within each subfamily. Moreover, 27 of ClVQ genes were found to have no introns. Conserved domain and multiple sequence alignment analysis revealed the presence of a highly conserved sequences in the ClVQ protein. This research utilized quantitative real-time PCR (qRT-PCR) and promoter analysis to investigate the expression of ClVQ genes under different stress conditions. Results showed that most ClVQ genes responded to polyethylene glycol, heat treatment, salt, abscisic acid and methyl jasmonate treatment with varying degrees of expression. Furthermore, some ClVQ genes exhibited significant correlation in expression changes under abiotic stress, indicating that these genes may act synergistically in response to adversarial stress. Additionally, yeast dihybrid verification revealed an interaction between ClVQ4, ClVQ12, and ClVQ26.

CONCLUSIONS

This study conducted a genome-wide analysis of the VQ gene family in coix, including an examination of phylogenetic relationships, conserved domains, cis-elements and expression patterns. The goal of the study was to identify potential drought resistance candidate genes, providing a theoretical foundation for molecular resistance breeding.

The Differences of Nutrient Components in Edible and Feeding Coix Seed at Different Developmental Stages Based on a Combined Analysis of Metabolomics

Coix lachryma-jobi L. is an excellent plant resource that has a concomitant function for medicine, foodstuff and forage in China. At present, the commonly used cultivar for both medicine and foodstuff is Xiaobaike, and the cultivar for foraging is Daheishan. However, differences in the internal composition of plants lead to the expression of different phenotypic traits. In order to comprehensively elucidate the differences in nutrient composition changes in Coix seeds, a non-targeted metabolomics method based on ultra-performance liquid chromatography quadrupole time-of-flight mass spectrometry (LC-Q-TOF-MS) was used to analyze the metabolic changes in Coix seeds at different developmental stages. An edible Coix relative (Xiaobaike) and a feeding Coix relative (Daheishan) were selected as the research subjects. In the metabolome analysis of Coix seed, 314 metabolites were identified and detected, among which organic acids, carbohydrates, lipids, nucleotides and flavonoids were the main components. As an important standard for evaluating the quality of Coix seed, seven lipids were detected, among which fatty acids included not only even-chain fatty acids, but also odd-chain fatty acids, which was the first time detecting a variety of odd-chain fatty acids in Coix seed. The analysis of the compound contents in edible and feeding-type Coix lachryma-jobi L. and the lipid content at the mature stage showed that, among them, arachidic acid, behenic acid, heptadecanoic acid, heneicosanoic acid and pristanic acid may be the key compounds affecting the lipid content. In addition, in the whole process of semen coicis maturation, edible and feeding Coix show similar trends, and changes in the third period show clear compounds in the opposite situation, suggesting that edible and feeding Coix not only guarantee the relative stability of species but also provide raw materials for genetic breeding. This study provides valuable information on the formation of the edible and medicinal qualities of Coix.

Efficient Anchoring of Erianthus arundinaceus Chromatin Introgressed into Sugarcane by Specific Molecular Markers

Erianthus arundinaceus is a valuable gene reservoir for sugarcane improvement. However, insufficient molecular markers for high-accuracy identification and tracking of the introgression status of E. arundinaceus chromatin impede sugarcane breeding. Fortunately, suppression subtractive hybridization (SSH) technology provides an excellent opportunity for the development of high-throughput E. arundinaceus-specific molecular markers at a reasonable cost. In this study, we constructed a SSH library of E. arundinaceus. In total, 288 clones of E. arundinaceus-specific repetitive sequences were screened out and their distribution patterns on chromosomes were characterized by fluorescence in situ hybridization (FISH). A subtelomeric repetitive sequence Ea086 and a diffusive repetitive sequence Ea009, plus 45S rDNA-bearing E. arundinaceus chromosome repetitive sequence EaITS were developed as E. arundinaceus-specific molecular markers, namely, Ea086-128, Ea009-257, and EaITS-278, covering all the E. arundinaceus chromosomes for high-accuracy identification of putative progeny. Both Ea086-128 and Ea009-257 were successfully applied to identify the authenticity of F1, BC1, BC2, BC3, and BC4 progeny between sugarcane and E. arundinaceus. In addition, EaITS-278 was a 45S rDNA-bearing E. arundinaceus chromosome-specific molecular marker for rapid tracking of the inherited status of this chromosome in a sugarcane background. Three BC3 progeny had apparently lost the 45S rDNA-bearing E. arundinaceus chromosome. We reported herein a highly effective and reliable SSH-based technology for discovery of high-throughput E. arundinaceus-specific sequences bearing high potential as molecular markers. Given its reliability and savings in time and efforts, the method is also suitable for development of species-specific molecular markers for other important wild relatives to accelerate introgression of wild relatives into sugarcane.

Beyond a reference genome: pangenomes and population genomics of underutilized and orphan crops for future food and nutrition security

Underutilized crops are, by definition, under-researched compared to staple crops yet come with traits that may be especially important given climate change and the need to feed a globally increasing population. These crops are often stress-tolerant, and this combined with unique and beneficial nutritional profiles. Whilst progress is being made by generating reference genome sequences, in this Tansley Review, we show how this is only the very first step. We advocate that going 'beyond a reference genome' should be a priority, as it is only at this stage one can identify the specific genes and the adaptive alleles that underpin the valuable traits. We sum up how population genomic and pangenomic approaches have led to the identification of stress- and disease-tolerant alleles in staple crops and compare this to the small number of examples from underutilized crops. We also demonstrate how previously underutilized crops have benefitted from genomic advances and that many breeding targets in underutilized crops are often well studied in staple crops. This cross-crop population-level resequencing could lead to an understanding of the genetic basis of adaptive traits in underutilized crops. This level of investment may be crucial for fully understanding the value of these crops before they are lost.

Comparative Repeat Profiling of Two Closely Related Conifers (Larix decidua and Larix kaempferi) Reveals High Genome Similarity With Only Few Fast-Evolving Satellite DNAs

In eukaryotic genomes, cycles of repeat expansion and removal lead to large-scale genomic changes and propel organisms forward in evolution. However, in conifers, active repeat removal is thought to be limited, leading to expansions of their genomes, mostly exceeding 10 giga base pairs. As a result, conifer genomes are largely littered with fragmented and decayed repeats. Here, we aim to investigate how the repeat landscapes of two related conifers have diverged, given the conifers' accumulative genome evolution mode. For this, we applied low-coverage sequencing and read clustering to the genomes of European and Japanese larch, Larix decidua (Lamb.) Carrière and Larix kaempferi (Mill.), that arose from a common ancestor, but are now geographically isolated. We found that both Larix species harbored largely similar repeat landscapes, especially regarding the transposable element content. To pin down possible genomic changes, we focused on the repeat class with the fastest sequence turnover: satellite DNAs (satDNAs). Using comparative bioinformatics, Southern, and fluorescent in situ hybridization, we reveal the satDNAs' organizational patterns, their abundances, and chromosomal locations. Four out of the five identified satDNAs are widespread in the Larix genus, with two even present in the more distantly related Pseudotsuga and Abies genera. Unexpectedly, the EulaSat3 family was restricted to L. decidua and absent from L. kaempferi, indicating its evolutionarily young age. Taken together, our results exemplify how the accumulative genome evolution of conifers may limit the overall divergence of repeats after speciation, producing only few repeat-induced genomic novelties.

Complete chloroplast genome sequencing of Job's tears (Coix L.): genome structure, comparative analysis, and phylogenetic relationships

Job's tears, also known as adlay, is a valuable plant that has commonly been used in traditional Chinese medicine, as well as an edible food. Due to the lack of knowledge of its genetics and gaps in its evolutionary analysis, breeding of adlay has been hindered. Here, we report five complete chloroplast genomes of various species and varieties in the genus by Illumina sequencing, while their genome structure, comparative analysis, and phylogenetic relationships were conducted. Genome sizes ranged from 140,860 to 140,864 bp in length, GC contents were 38.43%, and genome architecture was of a typical quadripartite structure. We annotated 82～83 protein-coding genes and 46～47 non-coding RNA genes in each genome and they functionally associated with self-replication, photosynthesis, cytochrome synthesis and other unknown functions. Three codons that encoded tryptophan, arginine and leucine were used frequently at rates of 41.42, 37.98, and 32.28% respectively. The preferred codons consistently ended with A or T. A total of 146 simple sequence repeats (SSR), 9 insertions and deletions (InDels) and 143 single nucleotide polymorphisms (SNPs) were observed among genomes. The InDel and SNP variations were mostly distributed in intergenic regions. It confirmed that Coix, Sorghum, Saccharum, Zea, Tripsacum and Saccharum were closely genera and the genetic distance of Sorghum to Coix was closer than Zea to Coix. These results give us more insight into the evolution of Coix in a wide range of evolutionary studies.

Genome evolution during bread wheat formation unveiled by the distribution dynamics of SSR sequences on chromosomes using FISH

BACKGROUND

During the bread wheat speciation by polyploidization, a series of genome rearrangement and sequence recombination occurred. Simple sequence repeat (SSR) sequences, predominately located in heterochromatic regions of chromosomes, are the effective marker for tracing the genomic DNA sequence variations. However, to date the distribution dynamics of SSRs on chromosomes of bread wheat and its donors, including diploid and tetraploid Triticum urartu, Aegilops speltoides, Aegilops tauschii, Triticum turgidum ssp. dicocoides, reflecting the genome evolution events during bread wheat formation had not been comprehensively investigated.

RESULTS

The genome evolution was studied by comprehensively comparing the distribution patterns of (AAC)_n, (AAG)_n, (AGC)_n and (AG)_n in bread wheat Triticum aestivum var. Chinese Spring and its progenitors T. urartu, A. speltoides, Ae. tauschii, wild tetroploid emmer wheat T. dicocoides, and cultivated emmer wheat T. dicoccum. Results indicated that there are specific distribution patterns in different chromosomes from different species for each SSRs. They provided efficient visible markers for identification of some individual chromosomes and SSR sequence evolution tracing from the diploid progenitors to hexaploid wheat. During wheat speciation, the SSR sequence expansion occurred predominately in the centromeric and pericentromeric regions of B genome chromosomes accompanied by little expansion and elimination on other chromosomes. This result indicated that the B genome might be more sensitive to the "genome shock" and more changeable during wheat polyplodization.

CONCLUSIONS

During the bread wheat evolution, SSRs including (AAC)_n, (AAG)_n, (AGC)_n and (AG)_n in B genome displayed the greatest changes (sequence expansion) especially in centromeric and pericentromeric regions during the polyploidization from Ae. speltoides S genome, the most likely donor of B genome. This work would enable a better understanding of the wheat genome formation and evolution and reinforce the viewpoint that B genome was originated from S genome.

New evidence confirming the CD genomic constitutions of the tetraploid Avena species in the section Pachycarpa Baum

The tetraploid Avena species in the section Pachycarpa Baum, including A. insularis, A. maroccana, and A. murphyi, are thought to be involved in the evolution of hexaploid oats; however, their genome designations are still being debated. Repetitive DNA sequences play an important role in genome structuring and evolution, so understanding the chromosomal organization and distribution of these sequences in Avena species could provide valuable information concerning genome evolution in this genus. In this study, the chromosomal organizations and distributions of six repetitive DNA sequences (including three SSR motifs (TTC, AAC, CAG), one 5S rRNA gene fragment, and two oat A and C genome specific repeats) were investigated using non-denaturing fluorescence in situ hybridization (ND-FISH) in the three tetraploid species mentioned above and in two hexaploid oat species. Preferential distribution of the SSRs in centromeric regions was seen in the A and D genomes, whereas few signals were detected in the C genomes. Some intergenomic translocations were observed in the tetraploids; such translocations were also detected between the C and D genomes in the hexaploids. These results provide robust evidence for the presence of the D genome in all three tetraploids, strongly suggesting that the genomic constitution of these species is DC and not AC, as had been thought previously.

Satellite DNA landscapes after allotetraploidization of quinoa (Chenopodium quinoa) reveal unique A and B subgenomes

If two related plant species hybridize, their genomes may be combined and duplicated within a single nucleus, thereby forming an allotetraploid. How the emerging plant balances two co-evolved genomes is still a matter of ongoing research. Here, we focus on satellite DNA (satDNA), the fastest turn-over sequence class in eukaryotes, aiming to trace its emergence, amplification, and loss during plant speciation and allopolyploidization. As a model, we used Chenopodium quinoa Willd. (quinoa), an allopolyploid crop with 2n = 4x = 36 chromosomes. Quinoa originated by hybridization of an unknown female American Chenopodium diploid (AA genome) with an unknown male Old World diploid species (BB genome), dating back 3.3-6.3 million years. Applying short read clustering to quinoa (AABB), C. pallidicaule (AA), and C. suecicum (BB) whole genome shotgun sequences, we classified their repetitive fractions, and identified and characterized seven satDNA families, together with the 5S rDNA model repeat. We show unequal satDNA amplification (two families) and exclusive occurrence (four families) in the AA and BB diploids by read mapping as well as Southern, genomic, and fluorescent in situ hybridization. Whereas the satDNA distributions support C. suecicum as possible parental species, we were able to exclude C. pallidicaule as progenitor due to unique repeat profiles. Using quinoa long reads and scaffolds, we detected only limited evidence of intergenomic homogenization of satDNA after allopolyploidization, but were able to exclude dispersal of 5S rRNA genes between subgenomes. Our results exemplify the complex route of tandem repeat evolution through Chenopodium speciation and allopolyploidization, and may provide sequence targets for the identification of quinoa's progenitors.

The Coix Genome Provides Insights into Panicoideae Evolution and Papery Hull Domestication

Coix is a grass crop domesticated as early as the Neolithic era. It is still widely cultivated for both highly nutritional food and medicinal use. However, the genetic study and breeding of this crop are hindered by the lack of a sequenced genome. Here, we report de novo sequencing and assembly of the 1619-Mb genome of Coix, and annotation of 75.39% repeats and 39 629 protein-coding genes. Comparative genomics analysis showed that Coix is more closely related to sorghum than maize, but intriguingly only Coix and maize had a recent genome duplication event, which was not detected in sorghum. We further constructed a genetic map and mapped several important traits, especially the strength of hull. Selection of papery hull (thin: easy dehulling) from the stony hull (thick: difficult dehulling) in wild progenitors was a key step in Coix domestication. The papery hull makes seed easier to process and germinate. Anatomic and global transcriptome analysis revealed that the papery hull is a result of inhibition of cell division and wall biogenesis. We also successfully demonstrated that seed hull pressure resistance is controlled by two major quantitative trait loci (QTLs), which are associated with hull thickness and color, respectively. The two QTLs were further fine mapped within intervals of 250 kb and 146 kb, respectively. These resources provide a platform for evolutionary studies and will facilitate molecular breeding of this important crop.

Evolution and Domestication Footprints Uncovered from the Genomes of Coix

Coix lacryma-jobi, a plant species closely related to Zea and Sorghum, is an important food and medicinal crop in Asia. However, no reference genome of this species has been reported, and its exact phylogeny within the Andropogoneae remains unresolved. Here, we generated a high-quality genome assembly of coix comprising ∼1.73 Gb with 44 485 predicted protein-coding genes. We found coix to be a typical diploid plant with an overall 1-to-1 syntenic relationship with the Sorghum genome, despite its drastic genome expansion (∼2.3-fold) due mainly to the activity of transposable elements. Phylogenetic analysis revealed that coix diverged with sorghum ∼10.41 million years ago, which was ∼1.49 million years later than the divergence between sorghum and maize. Resequencing of 27 additional coix accessions revealed that they could be unambiguously separated into wild relatives and cultivars, and suggested that coix experienced a strong genetic bottleneck, resulting in the loss of about half of the genetic diversity during domestication, even though many traits have remained undomesticated. Our data not only provide novel comparative genomic and evolutionary insights into the Andropogoneae lineage, but also an important resource that will greatly benefit molecular breeding of this important crop.

Integrated Karyotypes of Diploid and Tetraploid Carrizo Citrange (Citrus sinensis L. Osbeck × Poncirus trifoliata L. Raf.) as Determined by Sequential Multicolor Fluorescence in situ Hybridization With Tandemly Repeated DNA Sequences

Carrizo citrange [Citrus sinensis (L.) Osbeck × Poncirus trifoliata (L.) Raf., CC] is one of the most widely used rootstocks in citriculture worldwide, but its cytogenetic study has been hampered by its inherent small size, morphological similarity to mitotic chromosomes, and lack of accessible cytological landmarks. In our previous study, a spontaneously occurring tetraploid CC seedling was discovered. The main goals of this study were to elucidate the chromosome constitution and construct the karyotypes of diploid CC rootstock and its corresponding spontaneously occurring tetraploid. To accomplish these, the chromosomal characteristics were investigated by sequential multicolor fluorescence in situ hybridization (FISH) with eight properly labeled repetitive DNA sequences, including a centromere-like repeat, four satellite repeats, two rDNAs, and an oligonucleotide of telomeric (TTTAGGG) _n repeat. The results nicely demonstrated that these repetitive DNAs are reliable cytogenetic markers that collectively facilitate simultaneous and unequivocal identification of homologous chromosome pairs. Based on chromosome size and morphology together with FISH patterns of repetitive DNAs, an integrated karyotype of CC rootstock was constructed, consisting of 2n = 2x = 12m (1sat) + 6sm with karyotype asymmetry degree being divided into 2B category. Cytogenetically speaking, the variable and asymmetric distribution patterns of these repetitive DNAs were fully confirmed the hybrid nature of CC rootstock. In addition, comparative distribution patterns and chromosomal localizations of these repetitive DNAs convincingly showed that this tetraploid CC material arose from somatic chromosome doubling of diploid CC rootstock. This study revealed, for the first time, the integrated karyotype and chromosomal characteristics of this important citrus rootstock as well as its spontaneously occurring tetraploid plant. Furthermore, this study is a good prospective model for study species with morphologically indistinguishable small chromosomes.

Genome Assembly and Annotation of Soft-Shelled Adlay (Coix lacryma-jobi Variety ma-yuen), a Cereal and Medicinal Crop in the Poaceae Family

Coix lacryma-jobi, also called adlay or Job's tears, is an annual herbal plant belonging to the Poaceae family that has been cultivated as a cereal and medicinal crop in Asia. Despite its importance, however, genomic resources for better understanding this plant species at the molecular level and informing improved breeding strategies remain limited. To address this, we generated a draft genome of the C. lacryma-jobi variety ma-yuen (soft-shelled adlay) Korean cultivar, Johyun, by de novo assembly, using PacBio and Illumina sequencing data. A total of 3,362 scaffold sequences, 1.28 Gb in length, were assembled, representing 82.1% of the estimated genome size (1.56 Gb). Genome completeness was confirmed by the presence of 91.4% of the BUSCO angiosperm genes and mapping ratio of 98.3% of Illumina paired-end reads. We found that approximately 77.0% of the genome is occupied by repeat sequences, most of which are Gypsy and Copia-type retrotransposons, and evidence-based genome annotation predicts 39,574 protein-coding genes, 85.5% of which were functionally annotated. We further predict that soft-shelled adlay diverged from a common ancestor with sorghum 9.0-11.2 MYA. Transcriptome profiling revealed 3,988 genes that are differentially expressed in seeds relative to other tissues, of which 1,470 genes were strongly up-regulated in seeds and the most enriched Gene Ontology terms were assigned to carbohydrate and protein metabolism. In addition, we identified 76 storage protein genes including 18 seed-specific coixin genes and 13 candidate genes involved in biosynthesis of benzoxazinoids (BXs) including coixol, a unique BX compound found in C. lacryma-jobi species. The characterization of those genes can further our understanding of unique traits of soft-shelled adlay, such as high seed protein content and medicinal compound biosynthesis. Taken together, our genome sequence data will provide a valuable resource for molecular breeding and pharmacological study of this plant species.

Coix lacryma-jobi chymotrypsin inhibitor displays antifungal activity

A novel chymotrypsin inhibitor, named ClCI, was purified from coix seed (Coix lacryma-jobi L.) by aqueous two-phase extraction, chymotrypsin-Sepharose 4B affinity chromatography and centrifugal ultrafiltration. ClCI was a 7.9 kDa competitive inhibitor with pI 6.54. The inhibition constants (K_i) for bovine pancreatic chymotrypsin and bacterial subtilisin were 1.27 × 10^-10 M and 1.57 × 10^-9 M respectively. ClCI had no inhibitory activity against bovine trypsin and porcine elastase. ClCI had wide pH stability and good heat resistance. It can maintain >90% inhibition activity against chymotrypsin at 20-80 °C for 1 h. The primary structure of ClCI was highly similar (57%-92%) to those of several inhibitors belonging to the Gramineae crop potato protease inhibitor- I superfamily and showed the typical sequence motif of the protease inhibitor of the seed storage protein group. ClCI (12.5 mg) inhibited mycelial growth of the phytopathogenic fungi Mycosphaerella melonis, Helminthosporium turcicum, Alternaria solani, Phytophthora capsici, Isariopsis griseola, and Colletotrichum gloeosporioides, and caused 89% inhibition of the proteases from spore germination of plant-pathogenic fungi. The results of the present study indicate that ClCI had biotechnological potential as an alternative agent to combat the important phytopathogenic fungi.

The repetitive DNA landscape in Avena (Poaceae): chromosome and genome evolution defined by major repeat classes in whole-genome sequence reads

BACKGROUND

Repetitive DNA motifs - not coding genetic information and repeated millions to hundreds of times - make up the majority of many genomes. Here, we identify the nature, abundance and organization of all the repetitive DNA families in oats (Avena sativa, 2n = 6x = 42, AACCDD), a recognized health-food, and its wild relatives.

RESULTS

Whole-genome sequencing followed by k-mer and RepeatExplorer graph-based clustering analyses enabled assessment of repetitive DNA composition in common oat and its wild relatives' genomes. Fluorescence in situ hybridization (FISH)-based karyotypes are developed to understand chromosome and repetitive sequence evolution of common oat. We show that some 200 repeated DNA motifs make up 70% of the Avena genome, with less than 20 families making up 20% of the total. Retroelements represent the major component, with Ty3/Gypsy elements representing more than 40% of all the DNA, nearly three times more abundant than Ty1/Copia elements. DNA transposons are about 5% of the total, while tandemly repeated, satellite DNA sequences fit into 55 families and represent about 2% of the genome. The Avena species are monophyletic, but both bioinformatic comparisons of repeats in the different genomes, and in situ hybridization to metaphase chromosomes from the hexaploid species, shows that some repeat families are specific to individual genomes, or the A and D genomes together. Notably, there are terminal regions of many chromosomes showing different repeat families from the rest of the chromosome, suggesting presence of translocations between the genomes.

CONCLUSIONS

The relatively small number of repeat families shows there are evolutionary constraints on their nature and amplification, with mechanisms leading to homogenization, while repeat characterization is useful in providing genome markers and to assist with future assemblies of this large genome (c. 4100 Mb in the diploid). The frequency of inter-genomic translocations suggests optimum strategies to exploit genetic variation from diploid oats for improvement of the hexaploid may differ from those used widely in bread wheat.

In silico mining and FISH mapping of a chromosome-specific satellite DNA in Capsicum annuum L

BACKGROUND

A large proportion of eukaryote nuclear genomes is composed of repetitive DNA. Tracing the dynamics of repetitive elements in the genomes of related taxa can reveal important information about their phylogenic relationships as well as traits that have become distinct to a lineage.

OBJECTIVE

Study the genomic abundance and chromosomal location of repetitive DNA in Capsicum annuum L. to understand the repeat dynamics.

METHOD

We quantified repeated DNA content in the C. annuum genome using the RepeatExplorer pipeline.

RESULTS

About 42% of the C. annuum genome dataset comprised repetitive elements. Of these, 0.011, 0.98, 3.09, and 0.024% represented high and low confidence satellite repeats, putative long-terminal repeats (LTRs), and rDNA sequences, respectively. One novel high confidence 167-bp satellite repeat with a genomic proportion of 0.011%, Ca167TR, was identified. Furthermore, FISH with Ca167TR on metaphase chromosomes of C. annuum revealed signals in the subtelomeric regions of the short and long arms of chromosome 3 and 4, respectively.

CONCLUSION

Further understanding of the origin and associated functions of Ca167TR and other repeats in Capsicum will give us insights into the genomic relationships and functions of the genome.

Karyotype Analysis of Diploid and Spontaneously Occurring Tetraploid Blood Orange [Citrus sinensis (L.) Osbeck] Using Multicolor FISH With Repetitive DNA Sequences as Probes

Blood orange [Citrus sinensis (L.) Osbeck] has been increasingly appreciated by consumers worldwide owing to its brilliant red color, abundant anthocyanin and other health-promoting compounds. However, there is still relatively little known about its cytogenetic characteristics, probably because of the small size and similar morphology of metaphase chromosomes and the paucity of chromosomal landmarks. In our previous study, a naturally occurring tetraploid blood orange plant was obtained via seedling screening. Before this tetraploid germplasm can be manipulated into a citrus triploid seedless breeding program, it is of great importance to determine its chromosome characterization and composition. In the present study, an integrated karyotype of blood orange was constructed using sequential multicolor fluorescence in situ hybridization (FISH) with four satellite repeats, two ribosomal DNAs (rDNAs), a centromere-like repeat and an oligonucleotide of telomere repeat (TTTAGGG)3 as probes. Satellite repeats were preferentially located at the terminal regions of the chromosomes of blood orange. Individual somatic chromosome pairs of blood orange were unambiguously identified by repetitive DNA-based multicolor FISH. These probes proved to be effective chromosomal landmarks. The karyotype was formulated as 2n = 2x = 18 = 16m+2sm (1sat) with the karyotype asymmetry degree belonging to 2B. The chromosomal distribution pattern of these repetitive DNAs in this spontaneously occurring tetraploid was identical to that of the diploid, but the tetraploid carried twice the number of hybridization sites as the diploid, indicating a possible pathway involving the spontaneous duplication of chromosome sets in nucellar cells. Our work may facilitate the molecular cytogenetic study of blood orange and provide chromosomal characterization for the future utilization of this tetraploid germplasm in the service of seedless breeding programs.

Comparative karyotype analysis among six species of Ipomoea based on two newly identified repetitive sequences

Sweet potato is one of the most important crops worldwide; however, basic research in this crop is limited. In this study, we aimed to construct a detailed karyotype of six species of Ipomoea (hexaploid Ipomoea batatas and five related species, namely, one tetraploid, I. tabascana and four diploids, I. splendor-sylvae, I. trifida, I. tenuissima, and I. × leucantha) and understand the relationship among these species. Two satellite repeats (viz., Itf_1 and Itf_2) were identified from the diploid I. trifida genome sequence using RepeatExplorer on Galaxy. Together with the ribosomal DNA (rDNA), although without distinguishable chromosomes, a detailed karyotype was constructed for the six species. Our results showed a similar karyotype between I. tenuissima and I. × leucantha, indicating their close relationship. The signal distribution pattern of Itf_1, 45S rDNA combination, detected only in I. trifida, I. tabascana, and I. batatas, implied their close relationships. The chromosomes carrying 5S rDNA could be conserved among the six species as they always carried the Itf_2 signals, which generated a similar signal distribution pattern. The results enabled a detailed comparative cytogenetic analysis, providing valuable information to understand the relationship among these species and help assemble the genome sequence of the six species of Ipomoea.

Molecular cytogenetic analysis of genome-specific repetitive elements in Citrus clementina Hort. Ex Tan. and its taxonomic implications

BACKGROUND

Clementine mandarin (Citrus clementina Hort. ex Tan.) is one of the most famous and widely grown citrus cultivars worldwide. Variations in relation to the composition and distribution of repetitive DNA sequences that dominate greatly in eukaryote genomes are considered to be species-, genome-, or even chromosome-specific. Repetitive DNA-based fluorescence in situ hybridization (FISH) is a powerful tool for molecular cytogenetic study. However, to date few studies have involved in the repetitive elements and cytogenetic karyotype of Clementine.

RESULTS

A graph-based similarity sequence read clustering methodology was performed to analyze the repetitive DNA families in the Clementine genome. The bioinformatics analysis showed that repetitive DNAs constitute 41.95% of the Clementine genome, and the majority of repetitive elements are retrotransposons and satellite DNAs. Sequential multicolor FISH using a probe mix that contained CL17, four satellite DNAs, two rDNAs and an oligonucleotide of (TTTAGGG)₃ was performed with Clementine somatic metaphase chromosomes. An integrated karyotype of Clementine was established based on unequivocal and reproducible chromosome discriminations. The distribution patterns of these probes in several Citrus, Poncirus and Fortunella species were summarized through extensive FISH analyses. Polymorphism and heterozygosity were commonly observed in the three genera. Some asymmetrical FISH loci in Clementine were in agreement with its hybrid origin.

CONCLUSIONS

The composition and abundance of repetitive elements in the Clementine genome were reanalyzed. Multicolor FISH-based karyotyping provided direct visual proof of the heterozygous nature of Clementine chromosomes with conspicuous asymmetrical FISH hybridization signals. We detected some similar and variable distribution patterns of repetitive DNAs in Citrus, Poncirus, and Fortunella, which revealed notable conservation among these genera, as well as obvious polymorphism and heterozygosity, indicating the potential utility of these repetitive element markers for the study of taxonomic, phylogenetic and evolutionary relationships in the future.

De novo transcriptome assembly of the Chinese pearl barley, adlay, by full-length isoform and short-read RNA sequencing

Adlay (Coix lacryma-jobi) is a tropical grass that has long been used in traditional Chinese medicine and is known for its nutritional benefits. Recent studies have shown that vitamin E compounds in adlay protect against chronic diseases such as cancer and heart disease. However, the molecular basis of adlay's health benefits remains unknown. Here, we generated adlay gene sets by de novo transcriptome assembly using long-read isoform sequencing (Iso-Seq) and short-read RNA-Sequencing (RNA-Seq). The gene sets obtained from Iso-seq and RNA-seq contained 31,177 genes and 57,901 genes, respectively. We confirmed the validity of the assembled gene sets by experimentally analyzing the levels of prolamin and vitamin E biosynthesis-associated proteins in adlay plant tissues and seeds. We compared the screened adlay genes with known gene families from closely related plant species, such as rice, sorghum and maize. We also identified tissue-specific genes from the adlay leaf, root, and young and mature seed, and experimentally validated the differential expression of 12 randomly-selected genes. Our study of the adlay transcriptome will provide a valuable resource for genetic studies that can enhance adlay breeding programs in the future.

Organization and evolution of four differentially amplified tandem repeats in the Cucumis hystrix genome

Three subtelomeric satellites and one interstitial 5S rDNA were characterized in Cucumis hystrix, and the pericentromeric signals of two C. hystrix subtelomeric satellites along C. sativus chromosomes supported the hypothesis of chromosome fusion in Cucumis. Tandem repeats are chromosome structural fractions consisting of highly repetitive sequences organized in large tandem arrays in most eukaryotes. Differentiation of tandem repeats directly affects the chromosome structure, which contributes to species formation and evolution. Cucumis hystrix (2n = 2x = 24) is the only wild Cucumis species grouped into the same subgenus with C. sativus (2n = 2x = 14), hence its phylogenetic position confers a vital role for C. hystrix to understand the chromosome evolution in Cucumis. However, our knowledge of C. hystrix tandem repeats is insufficient for a detailed understanding of the chromosome evolution in Cucumis. Based on de novo tandem repeat characterization using bioinformatics and in situ hybridization (ISH), we identified and characterized four differentially amplified tandem repeats, Cucumis hystrix satellite 1-3 (CuhySat1-CuhySat3) located at the subtelomeric regions of all chromosomes, and Cucumis hystrix 5S (Cuhy5S) located at the interstitial regions of one single chromosome pair. Comparative ISH mapping using CuhySat1-3 and Cuhy5S revealed high homology of tandem repeats between C. hystrix and C. sativus. Intriguingly, we found signal distribution variations of CuhySat2 and CuhySat3 on C. sativus chromosomes. In comparison to their subtelomeric signal distribution on C. hystrix chromosomes, CuhySat3 showed a pericentromeric signal distribution and CuhySat2 showed both subtelomeric and pericentromeric signal distributions on C. sativus chromosomes. This detailed characterization of four C. hystrix tandem repeats significantly widens our knowledge of the C. hystrix chromosome structure, and the observed signal distribution variations will be helpful for understanding the chromosome evolution of Cucumis.

TAREAN: a computational tool for identification and characterization of satellite DNA from unassembled short reads

Satellite DNA is one of the major classes of repetitive DNA, characterized by tandemly arranged repeat copies that form contiguous arrays up to megabases in length. This type of genomic organization makes satellite DNA difficult to assemble, which hampers characterization of satellite sequences by computational analysis of genomic contigs. Here, we present tandem repeat analyzer (TAREAN), a novel computational pipeline that circumvents this problem by detecting satellite repeats directly from unassembled short reads. The pipeline first employs graph-based sequence clustering to identify groups of reads that represent repetitive elements. Putative satellite repeats are subsequently detected by the presence of circular structures in their cluster graphs. Consensus sequences of repeat monomers are then reconstructed from the most frequent k-mers obtained by decomposing read sequences from corresponding clusters. The pipeline performance was successfully validated by analyzing low-pass genome sequencing data from five plant species where satellite DNA was previously experimentally characterized. Moreover, novel satellite repeats were predicted for the genome of Vicia faba and three of these repeats were verified by detecting their sequences on metaphase chromosomes using fluorescence in situ hybridization.

TAREAN: a computational tool for identification and characterization of satellite DNA from unassembled short reads

Satellite DNA is one of the major classes of repetitive DNA, characterized by tandemly arranged repeat copies that form contiguous arrays up to megabases in length. This type of genomic organization makes satellite DNA difficult to assemble, which hampers characterization of satellite sequences by computational analysis of genomic contigs. Here, we present tandem repeat analyzer (TAREAN), a novel computational pipeline that circumvents this problem by detecting satellite repeats directly from unassembled short reads. The pipeline first employs graph-based sequence clustering to identify groups of reads that represent repetitive elements. Putative satellite repeats are subsequently detected by the presence of circular structures in their cluster graphs. Consensus sequences of repeat monomers are then reconstructed from the most frequent k-mers obtained by decomposing read sequences from corresponding clusters. The pipeline performance was successfully validated by analyzing low-pass genome sequencing data from five plant species where satellite DNA was previously experimentally characterized. Moreover, novel satellite repeats were predicted for the genome of Vicia faba and three of these repeats were verified by detecting their sequences on metaphase chromosomes using fluorescence in situ hybridization.

Advancing Eucalyptus Genomics: Cytogenomics Reveals Conservation of Eucalyptus Genomes

The genus Eucalyptus encloses several species with high ecological and economic value, being the subgenus Symphyomyrtus one of the most important. Species such as E. grandis and E. globulus are well characterized at the molecular level but knowledge regarding genome and chromosome organization is very scarce. Here we characterized and compared the karyotypes of three economically important species, E. grandis, E. globulus, and E. calmadulensis, and three with ecological relevance, E. pulverulenta, E. cornuta, and E. occidentalis, through an integrative approach including genome size estimation, fluorochrome banding, rDNA FISH, and BAC landing comprising genes involved in lignin biosynthesis. All karyotypes show a high degree of conservation with pericentromeric 35S and 5S rDNA loci in the first and third pairs, respectively. GC-rich heterochromatin was restricted to the 35S rDNA locus while the AT-rich heterochromatin pattern was species-specific. The slight differences in karyotype formulas and distribution of AT-rich heterochromatin, along with genome sizes estimations, support the idea of Eucalyptus genome evolution by local expansions of heterochromatin clusters. The unusual co-localization of both rDNA with AT-rich heterochromatin was attributed mainly to the presence of silent transposable elements in those loci. The cinnamoyl CoA reductase gene (CCR1) previously assessed to linkage group 10 (LG10) was clearly localized distally at the long arm of chromosome 9 establishing an unexpected correlation between the cytogenetic chromosome 9 and the LG10. Our work is novel and contributes to the understanding of Eucalyptus genome organization which is essential to develop successful advanced breeding strategies for this genus.

B chromosome contains active genes and impacts the transcription of A chromosomes in maize (Zea mays L.)

BACKGROUND

The dispensable maize (Zea mays L.) B chromosome is highly heterochromatic and widely believed to be devoid of functional genes. Although low-copy B chromosome causes no obvious phenotype variation, its existence might influence A genome gene expression. Previous studies suggested that B chromosomes are evolved from standard chromosomes; therefore, they might contain genic regions showing homology with A chromosome sequences.

RESULTS

Our data suggested that maize B chromosome influences the A-genome transcription with stronger effect associated with an increase in copy number of B chromosome. In total 130 differently expressed genes were detected in comparison between with and without B chromosome lines. These differentially expressed genes are mainly involved in cell metabolism and nucleotide binding. Using Starter + B, we amplified ten B chromosome loci with high sequence similarity to A-genome genes. Fluorescence in situ hybridization (FISH) confirmed that at least four ~5 kb-sized genes are located on the B chromosome. In addition, through de novo assembly of the reads not unmapped to maize B73 reference genome together with PCR validation, we found three B-located LTR; in particular, one of them, the 3.2 kb comp75688, is expressed in a B-dosage dependent manner.

CONCLUSION

We found that in the presence of maize B chromosome, the transcription of A genome genes was altered, with more impact by the increase of the B chromosome number. The B-located transcriptionally active genes showed high similarity to their A-genome homologues, and retrotransposons on B chromosome also have partial homologous to A genome sequences. Our data shed more lights on the genome structure and evolution of the maize B chromosome.

In Depth Characterization of Repetitive DNA in 23 Plant Genomes Reveals Sources of Genome Size Variation in the Legume Tribe Fabeae

The differential accumulation and elimination of repetitive DNA are key drivers of genome size variation in flowering plants, yet there have been few studies which have analysed how different types of repeats in related species contribute to genome size evolution within a phylogenetic context. This question is addressed here by conducting large-scale comparative analysis of repeats in 23 species from four genera of the monophyletic legume tribe Fabeae, representing a 7.6-fold variation in genome size. Phylogenetic analysis and genome size reconstruction revealed that this diversity arose from genome size expansions and contractions in different lineages during the evolution of Fabeae. Employing a combination of low-pass genome sequencing with novel bioinformatic approaches resulted in identification and quantification of repeats making up 55–83% of the investigated genomes. In turn, this enabled an analysis of how each major repeat type contributed to the genome size variation encountered. Differential accumulation of repetitive DNA was found to account for 85% of the genome size differences between the species, and most (57%) of this variation was found to be driven by a single lineage of Ty3/gypsy LTR-retrotransposons, the Ogre elements. Although the amounts of several other lineages of LTR-retrotransposons and the total amount of satellite DNA were also positively correlated with genome size, their contributions to genome size variation were much smaller (up to 6%). Repeat analysis within a phylogenetic framework also revealed profound differences in the extent of sequence conservation between different repeat types across Fabeae. In addition to these findings, the study has provided a proof of concept for the approach combining recent developments in sequencing and bioinformatics to perform comparative analyses of repetitive DNAs in a large number of non-model species without the need to assemble their genomes.

Genomics and Evolution in Traditional Medicinal Plants: Road to a Healthier Life

Medicinal plants have long been utilized in traditional medicine and ethnomedicine worldwide. This review presents a glimpse of the current status of and future trends in medicinal plant genomics, evolution, and phylogeny. These dynamic fields are at the intersection of phytochemistry and plant biology and are concerned with the evolution mechanisms and systematics of medicinal plant genomes, origin and evolution of the plant genotype and metabolic phenotype, interaction between medicinal plant genomes and their environment, the correlation between genomic diversity and metabolite diversity, and so on. Use of the emerging high-end genomic technologies can be expanded from crop plants to traditional medicinal plants, in order to expedite medicinal plant breeding and transform them into living factories of medicinal compounds. The utility of molecular phylogeny and phylogenomics in predicting chemodiversity and bioprospecting is also highlighted within the context of natural-product-based drug discovery and development. Representative case studies of medicinal plant genome, phylogeny, and evolution are summarized to exemplify the expansion of knowledge pedigree and the paradigm shift to the omics-based approaches, which update our awareness about plant genome evolution and enable the molecular breeding of medicinal plants and the sustainable utilization of plant pharmaceutical resources.

Repetitive sequence analysis and karyotyping reveals centromere-associated DNA sequences in radish (Raphanus sativus L.)

BACKGROUND

Radish (Raphanus sativus L., 2n = 2x = 18) is a major root vegetable crop especially in eastern Asia. Radish root contains various nutritions which play an important role in strengthening immunity. Repetitive elements are primary components of the genomic sequence and the most important factors in genome size variations in higher eukaryotes. To date, studies about repetitive elements of radish are still limited. To better understand genome structure of radish, we undertook a study to evaluate the proportion of repetitive elements and their distribution in radish.

RESULTS

We conducted genome-wide characterization of repetitive elements in radish with low coverage genome sequencing followed by similarity-based cluster analysis. Results showed that about 31% of the genome was composed of repetitive sequences. Satellite repeats were the most dominating elements of the genome. The distribution pattern of three satellite repeat sequences (CL1, CL25, and CL43) on radish chromosomes was characterized using fluorescence in situ hybridization (FISH). CL1 was predominantly located at the centromeric region of all chromosomes, CL25 located at the subtelomeric region, and CL43 was a telomeric satellite. FISH signals of two satellite repeats, CL1 and CL25, together with 5S rDNA and 45S rDNA, provide useful cytogenetic markers to identify each individual somatic metaphase chromosome. The centromere-specific histone H3 (CENH3) has been used as a marker to identify centromere DNA sequences. One putative CENH3 (RsCENH3) was characterized and cloned from radish. Its deduced amino acid sequence shares high similarities to those of the CENH3s in Brassica species. An antibody against B. rapa CENH3, specifically stained radish centromeres. Immunostaining and chromatin immunoprecipitation (ChIP) tests with anti-BrCENH3 antibody demonstrated that both the centromere-specific retrotransposon (CR-Radish) and satellite repeat (CL1) are directly associated with RsCENH3 in radish.

CONCLUSIONS

Proportions of repetitive elements in radish were estimated and satellite repeats were the most dominating elements. Fine karyotyping analysis was established which allow us to easily identify each individual somatic metaphase chromosome. Immunofluorescence- and ChIP-based assays demonstrated the functional significance of satellite and centromere-specific retrotransposon at centromeres. Our study provides a valuable basis for future genomic studies in radish.