Elimination of 5S DNA unit classes in newly formed allopolyploids of the genera Aegilops and Triticum

Evolution and origin of bread wheat

Bread wheat (Triticum aestivum, genome BBAADD) is a young hexaploid species formed only 8,500-9,000 years ago through hybridization between a domesticated free-threshing tetraploid progenitor, genome BBAA, and Aegilops tauschii, the diploid donor of the D subgenome. Very soon after its formation, it spread globally from its cradle in the fertile crescent into new habitats and climates, to become a staple food of humanity. This extraordinary global expansion was probably enabled by allopolyploidy that accelerated genetic novelty through the acquisition of new traits, new intergenomic interactions, and buffering of mutations, and by the attractiveness of bread wheat's large, tasty, and nutritious grain with high baking quality. New genome sequences suggest that the elusive donor of the B subgenome is a distinct (unknown or extinct) species rather than a mosaic genome. We discuss the origin of the diploid and tetraploid progenitors of bread wheat and the conflicting genetic and archaeological evidence on where it was formed and which species was its free-threshing tetraploid progenitor. Wheat experienced many environmental changes throughout its evolution, therefore, while it might adapt to current climatic changes, efforts are needed to better use and conserve the vast gene pool of wheat biodiversity on which our food security depends.

Prospects and challenges of epigenomics in crop improvement

BACKGROUND

The advent of high-throughput epigenome mapping techniques has ushered in a new era of multiomics with powerful tools now available to map and record genomic output at different levels. Integrating the different components of the epigenome from these multiomics measures allows investigations of cis-regulatory elements on a genome-scale. Mapping of chromatin state, chromatin accessibility dynamics, and higher-order chromatin structure enables a new level of understanding of cell fate determination, identity and function in normal growth and development, disease resistance, and yield.

OBJECTIVE

In this paper, the recent advances in epigenomics research of rice, maize, and wheat are reviewed, and the development trends of epigenomics of major crops in the coming years are projected.

METHODS

We highlight the role of epigenomics in regulating growth and development and identifying potential distal cis-regulatory elements in three major crops, and discuss the prospects and challenges for new epigenetics-mediated breeding technologies in crop improvement.

CONCLUSION

In this review, we summarize and analyze recent epigenomic advances in three major crops epigenomics and discuss possibilities and challenges for future research in the field.

Mosaic Arrangement of the 5S rDNA in the Aquatic Plant Landoltia punctata (Lemnaceae)

Duckweeds are a group of monocotyledonous aquatic plants in the Araceae superfamily, represented by 37 species divided into five genera. Duckweeds are the fastest growing flowering plants and are distributed around the globe; moreover, these plants have multiple applications, including biomass production, wastewater remediation, and making pharmaceutical proteins. Dotted duckweed (Landoltia punctata), the sole species in genus Landoltia, is one of the most resilient duckweed species. The ribosomal DNA (rDNA) encodes the RNA components of ribosomes and represents a significant part of plant genomes but has not been comprehensively studied in duckweeds. Here, we characterized the 5S rDNA genes in L. punctata by cloning and sequencing 25 PCR fragments containing the 5S rDNA repeats. No length variation was detected in the 5S rDNA gene sequence, whereas the nontranscribed spacer (NTS) varied from 151 to 524 bp. The NTS variants were grouped into two major classes, which differed both in nucleotide sequence and the type and arrangement of the spacer subrepeats. The dominant class I NTS, with a characteristic 12-bp TC-rich sequence present in 3-18 copies, was classified into four subclasses, whereas the minor class II NTS, with shorter, 9-bp nucleotide repeats, was represented by two identical sequences. In addition to these diverse subrepeats, class I and class II NTSs differed in their representation of cis-elements and the patterns of predicted G-quadruplex structures, which may influence the transcription of the 5S rDNA. Similar to related duckweed species in the genus Spirodela, L. punctata has a relatively low rDNA copy number, but in contrast to Spirodela and the majority of other plants, the arrangement of the 5S rDNA units demonstrated an unusual, heterogeneous pattern in L. punctata, as revealed by analyzing clones containing double 5S rDNA neighboring units. Our findings may further stimulate the research on the evolution of the plant rDNA and discussion of the molecular forces driving homogenization of rDNA repeats in concerted evolution.

Chromosome diversity in species of the genus Arachis, revealed by FISH and CMA/DAPI banding, and inferences about their karyotype differentiation

The species of the genus Arachis (Leguminosae) are ordered into nine sections. The assignment of genome types in this genus has been based on cross-compatibility analysis and molecular cytogenetic studies. The latter has also allowed karyotypically establishing well-defined genomes and reassigning the genome of several species. However, most of these studies have been focused mainly on the sections Arachis and Rhizomatosae. To increase the knowledge about the chromosome diversity of the whole genus, here we performed a detailed karyotype characterization of representative species of most of the sections and genomes of Arachis. This characterization included chromosome morphology, CMA/DAPI chromosome banding, and chromosome marker localization (rDNAloci and one satDNA sequence) by fluorescent in situ hybridization (FISH). Based on the data obtained and other previously published data, we established the karyotype similarities by cluster analysis and defined eleven karyotype groups. The grouping was partly coincident with the traditional genome assignment, except for some groups and some individual species. Karyotype similarities among some genomes were also found. The main characteristics of each karyotype group of Arachis were summarized. Together, our results provide information that may be beneficial for future cytogenetic and evolutionary studies, and also contribute to the identification of interspecific hybrids.

Fine organization of genomic regions tagged to the 5S rDNA locus of the bread wheat 5B chromosome

BACKGROUND

The multigene family encoding the 5S rRNA, one of the most important structurally-functional part of the large ribosomal subunit, is an obligate component of all eukaryotic genomes. 5S rDNA has long been a favored target for cytological and phylogenetic studies due to the inherent peculiarities of its structural organization, such as the tandem arrays of repetitive units and their high interspecific divergence. The complex polyploid nature of the genome of bread wheat, Triticum aestivum, and the technically difficult task of sequencing clusters of tandem repeats mean that the detailed organization of extended genomic regions containing 5S rRNA genes remains unclear. This is despite the recent progress made in wheat genomic sequencing. Using pyrosequencing of BAC clones, in this work we studied the organization of two distinct 5S rDNA-tagged regions of the 5BS chromosome of bread wheat.

RESULTS

Three BAC-clones containing 5S rDNA were identified in the 5BS chromosome-specific BAC-library of Triticum aestivum. Using the results of pyrosequencing and assembling, we obtained six 5S rDNA- containing contigs with a total length of 140,417 bp, and two sets (pools) of individual 5S rDNA sequences belonging to separate, but closely located genomic regions on the 5BS chromosome. Both regions are characterized by the presence of approximately 70-80 copies of 5S rDNA, however, they are completely different in their structural organization. The first region contained highly diverged short-type 5S rDNA units that were disrupted by multiple insertions of transposable elements. The second region contained the more conserved long-type 5S rDNA, organized as a single tandem array. FISH using probes specific to both 5S rDNA unit types showed differences in the distribution and intensity of signals on the chromosomes of polyploid wheat species and their diploid progenitors.

CONCLUSION

A detailed structural organization of two closely located 5S rDNA-tagged genomic regions on the 5BS chromosome of bread wheat has been established. These two regions differ in the organization of both 5S rDNA and the neighboring sequences comprised of transposable elements, implying different modes of evolution for these regions.

Asymmetric epigenetic modification and elimination of rDNA sequences by polyploidization in wheat

rRNA genes consist of long tandem repeats clustered on chromosomes, and their products are important functional components of the ribosome. In common wheat (Triticum aestivum), rDNA loci from the A and D genomes were largely lost during the evolutionary process. This biased DNA elimination may be related to asymmetric transcription and epigenetic modifications caused by the polyploid formation. Here, we observed both sets of parental nucleolus organizing regions (NORs) were expressed after hybridization, but asymmetric silencing of one parental NOR was immediately induced by chromosome doubling, and reversing the ploidy status could not reactivate silenced NORs. Furthermore, increased CHG and CHH DNA methylation on promoters was accompanied by asymmetric silencing of NORs. Enrichment of H3K27me3 and H3K9me2 modifications was also observed to be a direct response to increased DNA methylation and transcriptional inactivation of NOR loci. Both A and D genome NOR loci with these modifications started to disappear in the S4 generation and were completely eliminated by the S7 generation in synthetic tetraploid wheat. Our results indicated that asymmetric epigenetic modification and elimination of rDNA sequences between different donor genomes may lead to stable allopolyploid wheat with increased differentiation and diversity.

Transposon variation by order during allopolyploidisation between Brassica oleracea and Brassica rapa

Although many studies have shown that transposable element (TE) activation is induced by hybridisation and polyploidisation in plants, much less is known on how different types of TE respond to hybridisation, and the impact of TE-associated sequences on gene function. We investigated the frequency and regularity of putative transposon activation for different types of TE, and determined the impact of TE-associated sequence variation on the genome during allopolyploidisation. We designed different types of TE primers and adopted the Inter-Retrotransposon Amplified Polymorphism (IRAP) method to detect variation in TE-associated sequences during the process of allopolyploidisation between Brassica rapa (AA) and Brassica oleracea (CC), and in successive generations of self-pollinated progeny. In addition, fragments with TE insertions were used to perform Blast2GO analysis to characterise the putative functions of the fragments with TE insertions. Ninety-two primers amplifying 548 loci were used to detect variation in sequences associated with four different orders of TE sequences. TEs could be classed in ascending frequency into LTR-REs, TIRs, LINEs, SINEs and unknown TEs. The frequency of novel variation (putative activation) detected for the four orders of TEs was highest from the F1 to F2 generations, and lowest from the F2 to F3 generations. Functional annotation of sequences with TE insertions showed that genes with TE insertions were mainly involved in metabolic processes and binding, and preferentially functioned in organelles. TE variation in our study severely disturbed the genetic compositions of the different generations, resulting in inconsistencies in genetic clustering. Different types of TE showed different patterns of variation during the process of allopolyploidisation.

High penetrance of a pan-canina type rDNA family in intersection Rosa hybrids suggests strong selection of bivalent chromosomes in the section Caninae

All dogroses (Rosa sect. Caninae) are characterized by the peculiar canina meiosis in which genetic material is unevenly distributed between female and male gametes. The pan-canina rDNA family (termed beta) appears to be conserved in all dogroses analyzed so far. Here, we have studied rDNAs in experimental hybrids obtained from open pollination of F1 plants derived from 2 independent intersectional crosses between the pentaploid dogrose species (2n = 5x = 35) Rosa rubiginosa as female parent (producing 4x egg cells due to the unique asymmetrical canina meiosis) and the tetraploid (2n = 4x = 28) garden rose R. hybrida 'André Brichet' as male parent (producing 2x pollen after normal meiosis). We analyzed the structure of rDNA units by molecular methods [CAPS and extensive sequencing of internal transcribed spacers (ITS)] and determined the number of loci on chromosomes by FISH. FISH showed that R. rubiginosa and 'André Brichet' harbored 5 and 4 highly heteromorphic rDNA loci, respectively. In the second generation of hybrid lines, we observed a reduced number of loci (4 and 5 instead of the expected 6). In R. rubiginosa and 'André Brichet', 2-3 major ITS types were found which is consistent with a weak homogenization pressure maintaining high diversity of ITS types in this genus. In contrast to expectation (the null hypothesis of Mendelian inheritance of ITS families), we observed reduced ITS diversity in some individuals of the second generation which might derive from self-fertilization or from a backcross to R. rubiginosa. In these individuals, the pan-canina beta family appeared to be markedly enriched, while the paternal families were lost or diminished in copies. Although the mechanism of biased meiotic transmission of certain rDNA types is currently unknown, we speculate that the bivalent-forming chromosomes carrying the beta rDNA family exhibit extraordinary pairing efficiency and/or are subjected to strong selection in Caninae polyploids.

What does the nr5S DNA multigene family tell us about the genomic relationship between Dasypyrum breviaristatum and D. villosum (Triticeae: Poaceae)?

The genus Dasypyrum contains two species: the annual and widespread D. villosum (2x = 2n = 14) and the perennial and generally rare D. breviaristatum (2x = 2n = 14 and 4x = 2n = 28). The origin of the latter and its genome constitution have been subject of several studies. There is agreement that the genome of the diploid D. villosum (VV) is different from the diploid cytotype of D. breviaristatum (VbVb), but there is no agreement of the constitution of the tetraploid cytotype, specifically whether is it an autotetraploid or an allotetraploid. This is a long-standing disagreement that this study aims to resolve using the 5S nrDNA as a genomic marker. Our studies suggest that the 4x D. breviaristatum is an allotetraploid (VVVbVb).

Genomic asymmetry in allopolyploid plants: wheat as a model

The evolvement of duplicated gene loci in allopolyploid plants has become the subject of intensive studies. Most duplicated genes remain active in neoallopolyploids contributing either to a favourable effect of an extra gene dosage or to the build-up of positive inter-genomic interactions when genes or regulation factors on homoeologous chromosomes are divergent. However, in a small number of loci (about 10%), genes of only one genome are active, while the homoeoalleles on the other genome(s) are either eliminated or partially or completely suppressed by genetic or epigenetic means. For several traits, the retention of controlling genes is not random, favouring one genome over the other(s). Such genomic asymmetry is manifested in allopolyploid wheat by the control of various morphological and agronomical traits, in the production of rRNA and storage proteins, and in interaction with pathogens. It is suggested that the process of cytological diploidization leading to exclusive intra-genomic meiotic pairing and, consequently, to complete avoidance of inter-genomic recombination, has two contrasting effects. Firstly, it provides a means for the fixation of positive heterotic inter-genomic interactions and also maintains genomic asymmetry resulting from loss or silencing of genes. The possible mechanisms and evolutionary advantages of genomic asymmetry are discussed.

Chromosomal characterization in native populations of Elymus scabrifolius from Argentina through classical and molecular cytogenetics (FISH–GISH)

The karyotype of Elymus scabrifolius (Döll) J.H. Hunz. (2n = 4x = 28) was investigated by DAPI staining and in situ hybridization. All the accessions studied presented a symmetric and uniform karyotype constituted by 9m+2m–sm+3sm. DAPI stain showed 1–7 conspicuous bands in all the chromosomes and polymorphisms between accessions. FISH experiments carried out with 45S rDNA as probe (pTa71) showed strong hybridization signals on the metacentric SAT-chromosome pair 8; the submetacentric SAT-chromosome pair 13 presented weaker hybridization. FISH using pSc119.2 clone as probe identified five chromosome pairs. Then, the combination of chromosome morphology, DAPI-staining, and FISH enabled the accurate identification of each chromosome pair in E. scabrifolius. Genomic in situ hybridization (GISH) experiments using Hordeum DNA as probe on mitotic metaphases confirmed unequivocally the presence of the H genome in E. scabrifolius, allowing us to observe six uniformly labeled chromosome pairs and two chromosome pairs with only one arm labeled. The remaining six chromosome pairs were weakly labeled. The rehybridization of FISH slides with Hordeum DNA as probe allow us to assign the genomic provenance of most of the chromosomes in the studied accessions. Moreover, intergenomic rearrangement was detected between genome H and the still unknown progenitor genome.

Identification of a major QTL controlling the content of B-type starch granules in Aegilops

Starch within the endosperm of most species of the Triticeae has a unique bimodal granule morphology comprising large lenticular A-type granules and smaller near-spherical B-type granules. However, a few wild wheat species (Aegilops) are known to lack B-granules. Ae. peregrina and a synthetic tetraploid Aegilops with the same genome composition (SU) were found to differ in B-granule number. The synthetic tetraploid had normal A- and B-type starch granules whilst Ae. peregrina had only A-granules because the B-granules failed to initiate. A population segregating for B-granule number was generated by crossing these two accessions and was used to study the genetic basis of B-granule initiation. A combination of Bulked Segregant Analysis and QTL mapping identified a major QTL located on the short arm of chromosome 4S that accounted for 44.4% of the phenotypic variation. The lack of B-granules in polyploid Aegilops with diverse genomes suggests that the B-granule locus has been lost several times independently during the evolution of the Triticeae. It is proposed that the B-granule locus is susceptible to silencing during polyploidization and a model is presented to explain the observed data based on the assumption that the initiation of B-granules is controlled by a single major locus per haploid genome.