Phylogenetic relationships among the polyploid and diploid Aegilops species inferred from the nuclear 5S rDNA sequences (Poaceae: Triticeae)

5S Ribosomal DNA of Genus Solanum: Molecular Organization, Evolution, and Taxonomy

The Solanum genus, being one of the largest among high plants, is distributed worldwide and comprises about 1,200 species. The genus includes numerous agronomically important species such as Solanum tuberosum (potato), Solanum lycopersicum (tomato), and Solanum melongena (eggplant) as well as medical and ornamental plants. The huge Solanum genus is a convenient model for research in the field of molecular evolution and structural and functional genomics. Clear knowledge of evolutionary relationships in the Solanum genus is required to increase the effectiveness of breeding programs, but the phylogeny of the genus is still not fully understood. The rapidly evolving intergenic spacer region (IGS) of 5S rDNA has been successfully used for inferring interspecific relationships in several groups of angiosperms. Here, combining cloning and sequencing with bioinformatic analysis of genomic data available in the SRA database, we evaluate the molecular organization and diversity of IGS for 184 accessions, representing 137 species of the Solanum genus. It was found that the main mechanisms of IGS molecular evolution was step-wise accumulation of single base substitution or short indels, and that long indels and multiple base substitutions, which arose repeatedly during evolution, were mostly not conserved and eliminated. The reason for this negative selection seems to be association between indels/multiple base substitutions and pseudogenization of 5S rDNA. Comparison of IGS sequences allowed us to reconstruct the phylogeny of the Solanum genus. The obtained dendrograms are mainly congruent with published data: same major and minor clades were found. However, relationships between these clades and position of some species (S. cochoae, S. clivorum, S. macrocarpon, and S. spirale) were different from those of previous results and require further clarification. Our results show that 5S IGS represents a convenient molecular marker for phylogenetic studies on the Solanum genus. In particular, the simultaneous presence of several structural variants of rDNA in the genome enables the detection of reticular evolution, especially in the largest and economically most important sect. Petota. The origin of several polyploid species should be reconsidered.

Identification of 5S and 45S rDNA sites in Chrysanthemum species by using oligonucleotide fluorescence in situ hybridization (Oligo-FISH)

Fluorescence in situ hybridization (FISH) is a conventional method used to visualize the distribution of DNA elements within a genome. To examine the relationships within the Chrysanthemum genus, ribosomal DNA (rDNA), a popular cytogenetic marker, was utilized as a probe for FISH within this genus. Based on the genome data of Chrysanthemum nankingense, C. seticuspe and its allied genera in the Compositae(Asteraceae), we explored rDNA sequences to design oligonucleotide probes and perform oligonucleotide fluorescence in situ hybridization (Oligo-FISH) in eight Chrysanthemum accessions. The results showed that the majority of 5S rDNA signals were located in subterminal chromosome regions and that the number of 5S rDNA sites might be tightly associated with ploidy. For 45S rDNA sites, the number and intensity of signals differed from those of previously investigated Chrysanthemum resources. These findings may provide an optimally reliable method of examining the chromosome composition and structural variation of Chrysanthemum and its related species and allow researchers to understand the evolutionary history and phylogenetic relationships of Chrysanthemum.

The Utility of Graph Clustering of 5S Ribosomal DNA Homoeologs in Plant Allopolyploids, Homoploid Hybrids, and Cryptic Introgressants

INTRODUCTION

Ribosomal DNA (rDNA) loci have been widely used for identification of allopolyploids and hybrids, although few of these studies employed high-throughput sequencing data. Here we use graph clustering implemented in the RepeatExplorer (RE) pipeline to analyze homoeologous 5S rDNA arrays at the genomic level searching for hybridogenic origin of species. Data were obtained from more than 80 plant species, including several well-defined allopolyploids and homoploid hybrids of different evolutionary ages and from widely dispersed taxonomic groups.

RESULTS

(i) Diploids show simple circular-shaped graphs of their 5S rDNA clusters. In contrast, most allopolyploids and other interspecific hybrids exhibit more complex graphs composed of two or more interconnected loops representing intergenic spacers (IGS). (ii) There was a relationship between graph complexity and locus numbers. (iii) The sequences and lengths of the 5S rDNA units reconstituted in silico from k-mers were congruent with those experimentally determined. (iv) Three-genomic comparative cluster analysis of reads from allopolyploids and progenitor diploids allowed identification of homoeologous 5S rRNA gene families even in relatively ancient (c. 1 Myr) Gossypium and Brachypodium allopolyploids which already exhibit uniparental partial loss of rDNA repeats. (v) Finally, species harboring introgressed genomes exhibit exceptionally complex graph structures.

CONCLUSION

We found that the cluster graph shapes and graph parameters (k-mer coverage scores and connected component index) well-reflect the organization and intragenomic homogeneity of 5S rDNA repeats. We propose that the analysis of 5S rDNA cluster graphs computed by the RE pipeline together with the cytogenetic analysis might be a reliable approach for the determination of the hybrid or allopolyploid plant species parentage and may also be useful for detecting historical introgression events.

Identification of New Leaf Rust Resistance Loci in Wheat and Wild Relatives by Array-Based SNP Genotyping and Association Genetics

Leaf rust caused by Puccinia triticina is the most widespread rust disease of wheat. As pathogen populations are constantly evolving, identification of novel sources of resistance is necessary to maintain disease resistance and stay ahead of this plant-pathogen evolutionary arms race. The wild genepool of wheat is a rich source of genetic diversity, accounting for 44% of the Lr genes identified. Here we performed a genome-wide association study (GWAS) on a diverse germplasm of 385 accessions, including 27 different Triticum and Aegilops species. Genetic characterization using the wheat 90 K array and subsequent filtering identified a set of 20,501 single nucleotide polymorphic (SNP) markers. Of those, 9,570 were validated using exome capture and mapped onto the Chinese Spring reference sequence v1.0. Phylogenetic analyses illustrated four major clades, clearly separating the wild species from the T. aestivum and T. turgidum species. GWAS was conducted using eight statistical models for infection types against six leaf rust isolates and leaf rust severity rated in field trials for 3-4 years at 2-3 locations in Canada. Functional annotation of genes containing significant quantitative trait nucleotides (QTNs) identified 96 disease-related loci associated with leaf rust resistance. A total of 21 QTNs were in haplotype blocks or within flanking markers of at least 16 known Lr genes. The remaining significant QTNs were considered loci that putatively harbor new Lr resistance genes. Isolation of these candidate genes will contribute to the elucidation of their role in leaf rust resistance and promote their usefulness in marker-assisted selection and introgression.

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.

Molecular Cytogenetic Mapping of Satellite DNA Sequences in Aegilops geniculata and Wheat

Fluorescence in situ hybridization (FISH) provides an efficient system for cytogenetic analysis of wild relatives of wheat for individual chromosome identification, elucidation of homoeologous relationships, and for monitoring alien gene transfers into wheat. This study is aimed at developing cytogenetic markers for chromosome identification of wheat and Aegilops geniculata (2n = 4x = 28, UgUgMgMg) using satellite DNAs obtained from flow-sorted chromosome 5Mg. FISH was performed to localize the satellite DNAs on chromosomes of wheat and selected Aegilops species. The FISH signals for satellite DNAs on chromosome 5Mg were generally associated with constitutive heterochromatin regions corresponding to C-band-positive chromatin including telomeric, pericentromeric, centromeric, and interstitial regions of all the 14 chromosome pairs of Ae. geniculata. Most satellite DNAs also generated FISH signals on wheat chromosomes and provided diagnostic chromosome arm-specific cytogenetic markers that significantly improved chromosome identification in wheat. The newly identified satellite DNA CL36 produced localized Mg genome chromosome-specific FISH signals in Ae. geniculata and in the M genome of the putative diploid donor species Ae. comosa subsp. subventricosa but not in Ae. comosa subsp. comosa, suggesting that the Mg genome of Ae. geniculata was probably derived from subsp. subventricosa.

Genome reorganization in F1 hybrids uncovers the role of retrotransposons in reproductive isolation

Interspecific hybridization leads to new interactions among divergent genomes, revealing the nature of genetic incompatibilities having accumulated during and after the origin of species. Conflicts associated with misregulation of transposable elements (TEs) in hybrids expectedly result in their activation and genome-wide changes that may be key to species boundaries. Repetitive genomes of wild wheats have diverged under differential dynamics of specific long terminal repeat retrotransposons (LTR-RTs), offering unparalleled opportunities to address the underpinnings of plant genome reorganization by selfish sequences. Using reciprocal F1 hybrids between three Aegilops species, restructuring and epigenetic repatterning was assessed at random and LTR-RT sequences with amplified fragment length polymorphism and sequence-specific amplified polymorphisms as well as their methylation-sensitive counterparts, respectively. Asymmetrical reorganization of LTR-RT families predicted to cause conflicting interactions matched differential survival of F1 hybrids. Consistent with the genome shock model, increasing divergence of merged LTR-RTs yielded higher levels of changes in corresponding genome fractions and lead to repeated reorganization of LTR-RT sequences in F1 hybrids. Such non-random reorganization of hybrid genomes is coherent with the necessary repression of incompatible TE loci in support of hybrid viability and indicates that TE-driven genomic conflicts may represent an overlooked factor supporting reproductive isolation.

Diversity within the genus Elymus (Poaceae: Triticeae) II: analyses of variation within 5S nrDNA restrict membership in the genus to species with StH genomes

The genus Elymus is a repository for a large number of species that have been difficult to classify by traditional techniques due to their remarkable levels of polymorphism. Following the genome analyses of Yen and Yang (Genus Elymus 5:58-362, 2013), we used sequences of the nr5SDNA to investigate diversity within those 24 species having St and H haplomes (Baum et al. Mol Genet Genomics 290:329-42, 2015) and for which the genome status was known. The present work extends this analysis to include eight species for which there was no information on genomic status. Our results show that these eight have nr5SDNA sequences that can be assigned to unit classes of orthologous sequences found in St and H haplomes, suggesting that the presence of St and H haplomes is characteristic of the genus. We then carried out a set of canonical discriminant analyses based on 247 DNA new sequences from these 8 species plus the 1054 sequences previously identified from 24 Elymus species. Sequences were analyzed to answer the following questions: Do the species integrate or are they different? Are the tetraploids different from the higher-ploid species? Are the species united within sections, or the same within regions? How do the species fare when divided according to sections? The main results of the canonical discriminant analyses are that the species are united within the tetraploids and within the hexaploids, within each region and within each section. In addition, a series of classificatory discriminant analyses showed that the identification tests are different, although not sufficiently useful for the discrimination of all the species. We also demonstrate the power of our approach by showing that the voucher for Elymus mobilis is not Elymus at all, but Leymus.

Diversity within the genus Elymus (Poaceae: Triticeae) as investigated by the analysis of the nr5S rDNA variation in species with St and H haplomes

The genus Elymus ("Ryegrass") is a repository for a range of species with a variety of haplome contents; hence the pejorative name "dustbin" genus. We have analyzed 1,059 sequences from 128 accessions representing 24 species to investigate the relationships among the StH haplomes-containing species described by Yen and Yang (Genus Elymus Beijing 5:58-362, 2013). Sequences were assigned to "unit classes" of orthologous sequences and subjected to a suite of analyses including BLAST (Basic Local Alignment Search Tool) searches, phylogenetic analysis and population genetic analysis to estimate species diversity. Our results support the genome analyses in Yen and Yang (Genus Elymus Beijing 5:58-362, 2013), i.e., genomic constitution StStHH including variants restricted to Elymus. Population genetic analysis of the 5S nrDNA sequence data revealed that the within-species variance component is roughly ±89 %; thus, we were unable to identify molecular markers capable to separate the 24 species analyzed. Separate phylogenetic analyses of the two unit classes and of all the data exhibit a trend only of the species to cluster on the phylograms. Finally, the analysis provides evidence for the multiple origins of American and Eurasian species.

Contrasting evolutionary trajectories of multiple retrotransposons following independent allopolyploidy in wild wheats

Transposable elements (TEs) are expectedly central to genome evolution. To assess the impact of TEs in driving genome turnover, we used allopolyploid genomes, showing considerable deviation from the predicted additivity of their diploid progenitors and thus having undergone major restructuring. Genome survey sequencing was used to select 17 putatively active families of long terminal repeat retrotransposons. Genome-wide TE insertions were genotyped with sequence-specific amplified polymorphism (SSAP) in diploid progenitors and their derived polyploids, and compared with changes in random sequences to assess restructuring of four independent Aegilops allotetraploid genomes. Generally, TEs with different evolutionary trajectories from those of random sequences were identified. Thus, TEs presented family-specific and species-specific dynamics following polyploidy, as illustrated by Sabine showing proliferation in particular polyploids, but massive elimination in others. Contrasting with that, only a few families (BARE1 and Romani) showed proliferation in all polyploids. Overall, TE divergence between progenitors was strongly correlated with the degree of restructuring in polyploid TE fractions. TE families present evolutionary trajectories that are decoupled from genome-wide changes after allopolyploidy and have a pervasive impact on their restructuring.

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).

What does the 5S rRNA multigene family tell us about the origin of the annual Triticeae (Poaceae)?

We have investigated the complex relationships among the annual genera within the tribe Triticeae through phylogenetic analyses of the 5S rRNA multigene family. Cloned sequences were assigned to groups of orthologous sequences, called unit classes, that were subjected to several analyses including BLAST (Basic Local Alignment Search Tool) searches to assess possible ancestral relationships with perennial genera; phylogenetic analyses using parsimony (Pars), maximum likelihood (ML), and Bayesian methods; and minimum reticulation networks from the Pars, ML, and Bayesian trees. In this study, we included genera with both annual and perennial species, such as Dasypyrum, Hordeum, and Secale. BLAST pointed to Pseudoroegneria (carrier of the St genome) and possibly Thinopyrum (carrier of the J genome) as the potential next of kin. However, Thinopyrum and Pseudoroegneria have never fallen together on the individual trees with the former generally associated with Crithopsis, Aegilops, Triticum, and Dasypyrum, while the latter is usually associated with the rest of the genera within Triticeae. The "long" unit classes placed Dasypyrum breviaristatum together with Dasypyrum villosum, whereas the "short" unit classes put them far apart on the trees. None of the gene trees alone was able to summarize the complex relationships among the genera, in line with previous results in the Triticeae. However, the application of tools designed to display phylogenetic networks was able to depict the complex links among the genera based on the short and the long gene trees, including the close link between Thinopyrum and Pseudoroegneria suggested by the phylogenetic analyses. In addition, our analyses provide support for the hypothesis that at least some annual Triticeae taxa are derived from their perennial relatives.

Evolutionary dynamics of retrotransposons assessed by high-throughput sequencing in wild relatives of wheat

Transposable elements (TEs) represent a major fraction of plant genomes and drive their evolution. An improved understanding of genome evolution requires the dynamics of a large number of TE families to be considered. We put forward an approach bypassing the required step of a complete reference genome to assess the evolutionary trajectories of high copy number TE families from genome snapshot with high-throughput sequencing. Low coverage sequencing of the complex genomes of Aegilops cylindrica and Ae. geniculata using 454 identified more than 70% of the sequences as known TEs, mainly long terminal repeat (LTR) retrotransposons. Comparing the abundance of reads as well as patterns of sequence diversity and divergence within and among genomes assessed the dynamics of 44 major LTR retrotransposon families of the 165 identified. In particular, molecular population genetics on individual TE copies distinguished recently active from quiescent families and highlighted different evolutionary trajectories of retrotransposons among related species. This work presents a suite of tools suitable for current sequencing data, allowing to address the genome-wide evolutionary dynamics of TEs at the family level and advancing our understanding of the evolution of nonmodel genomes.