Allopolyploid origin of Mediterranean species in Helictotrichon (Poaceae) and its consequences for karyotype repatterning and homogenisation of rDNA repeat units

Contrasting evolutionary origins of two mountain endemics: Saxifraga wahlenbergii (Western Carpathians) and S. styriaca (Eastern Alps)

BACKGROUND

The Carpathians and the Alps are the largest mountain ranges of the European Alpine System and important centres of endemism. Among the distinctive endemic species of this area is Saxifraga wahlenbergii, a Western Carpathians member of the speciose genus Saxifraga. It was frequently considered a taxonomically isolated Tertiary palaeopolyploid and palaeoendemic, for which the closest relatives could not yet be traced. A recently described narrow endemic of the Eastern Alps, S. styriaca, was hypothesized to be closely related to S. wahlenbergii based on shared presence of peculiar glandular hairs. To elucidate the origin and phylogenetic relationships of both species we studied nuclear and plastid DNA markers based on multiple accessions and analysed the data in a wide taxonomic context. We applied Sanger sequencing, followed by targeted next-generation sequencing (NGS) for a refined analysis of nrITS variants to detect signatures of ancient hybridization. The ITS data were used to estimate divergence times of different lineages using a relaxed molecular clock.

RESULTS

We demonstrate divergent evolutionary histories for the two mountain endemics. For S. wahlenbergii we revealed a complicated hybrid origin. Its maternal parent belongs to a Western Eurasian lineage of high mountain taxa grouped in subsect. Androsaceae and is most likely the widespread S. androsacea. The putative second parent was most likely S. adscendens, which belongs to the distantly related subsect. Tridactylites. While Sanger sequencing of nrITS only showed S. adscendens-related variants in S. wahlenbergii, our NGS screening revealed presence of sequences from both lineages with clear predominance of the paternal over the maternal lineage.

CONCLUSIONS

Saxifraga styriaca was unambiguously assigned to subsect. Androsaceae and is not the sister taxon of S. wahlenbergii. Accordingly, the similarity of the glandular hairs observed in both taxa rests on parallelism and both species do not constitute an example of a close evolutionary link between the floras of the Western Carpathians and Eastern Alps. With the origin of its paternal, S. adscendens-like ITS DNA estimated to ca. 4.7 Ma, S. wahlenbergii is not a relict of the mid-Tertiary climate optimum. Its hybrid origin is much younger and most likely took place in the Pleistocene.

Glacial refugia and speciation in a group of wind-pollinated and -dispersed, endemic Alpine species of Helictotrichon (Poaceae)

In the Alps phylogeographic studies indicate for small insect-pollinated herbs that climatic fluctuations caused significant population migrations and fragmentations into glacial refugia at the periphery of the Alps. Here we investigate whether this holds also for wind-pollinated and -dispersed species. We therefore analysed the phylogeographic pattern (nuclear and chloroplast dataset) of a clade of the four species of the Helictotrichon parlatorei species group (Poaceae) endemic to the Alps. In contrast to earlier findings for small insect-pollinated herbs no clear barriers to gene flow could be detected in this species group. Instead a few haplotypes are widespread across the entire Alpine region. While the complete absence of a phylogeographic structure in the plastid dataset hints towards very efficient long distance seed dispersal, the moderate phylogeographic structure in the nuclear dataset indicates at least some spatial restriction to pollen dispersal. Rare haplotypes cluster solely in the Western and Southern central Alps and thereby suggest this to be the area of origin for the H. parlatorei species group from where expansion occurred following the presence of calcareous bedrock into the Eastern Alps. We thus conclude that the inclusion of taxa with complementary life-history traits is vital in understanding the glacial history of the Alpine flora.

Dynamics of tandemly repeated DNA sequences during evolution of diploid and tetraploid botiid loaches (Teleostei: Cobitoidea: Botiidae)

Polyploidization has played an important role in the evolution of vertebrates, particularly at the base of Teleostei-an enormously successful ray-finned fish group with additional genome doublings on lower taxonomic levels. The investigation of post-polyploid genome dynamics might provide important clues about the evolution and ecology of respective species and can help to decipher the role of polyploidy per se on speciation. Few studies have attempted to investigate the dynamics of repetitive DNA sequences in the post-polyploid genome using molecular cytogenetic tools in fishes, though recent efforts demonstrated their usefulness. The demonstrably monophyletic freshwater loach family Botiidae, branching to evolutionary diploid and tetraploid lineages separated >25 Mya, offers a suited model group for comparing the long-term repetitive DNA evolution. For this, we integrated phylogenetic analyses with cytogenetical survey involving Giemsa- and Chromomycin A3 (CMA3)/DAPI stainings and fluorescence in situ hybridization with 5S/45S rDNA, U2 snDNA and telomeric probes in representative sample of 12 botiid species. The karyotypes of all diploids were composed of 2n = 50 chromosomes, while majority of tetraploids had 2n = 4x = 100, with only subtle interspecific karyotype differences. The exceptional karyotype of Botia dario (2n = 4x = 96) suggested centric fusions behind the 2n reduction. Variable patterns of FISH signals revealed cases of intraspecific polymorphisms, rDNA amplification, variable degree of correspondence with CMA3+ sites and almost no phylogenetic signal. In tetraploids, either additivity or loci gain/loss was recorded. Despite absence of classical interstitial telomeric sites, large blocks of interspersed rDNA/telomeric regions were found in diploids only. We uncovered different molecular drives of studied repetitive DNA classes within botiid genomes as well as the advanced stage of the re-diploidization process in tetraploids. Our results may contribute to link genomic approach with molecular cytogenetic analyses in addressing the origin and mechanism of this polyploidization event.

Karyotype evolution in Phalaris (Poaceae): The role of reductional dysploidy, polyploidy and chromosome alteration in a wide-spread and diverse genus

Karyotype characteristics can provide valuable information on genome evolution and speciation, in particular in taxa with varying basic chromosome numbers and ploidy levels. Due to its worldwide distribution, remarkable variability in morphological traits and the fact that ploidy change plays a key role in its evolution, the canary grass genus Phalaris (Poaceae) is an excellent study system to investigate the role of chromosomal changes in species diversification and expansion. Phalaris comprises diploid species with two basic chromosome numbers of x = 6 and 7 as well as polyploids based on x = 7. To identify distinct karyotype structures and to trace chromosome evolution within the genus, we apply fluorescence in situ hybridisation (FISH) of 5S and 45S rDNA probes in four diploid and four tetraploid Phalaris species of both basic numbers. The data agree with a dysploid reduction from x = 7 to x = 6 as the result of reciprocal translocations between three chromosomes of an ancestor with a diploid chromosome complement of 2n = 14. We recognize three different genomes in the genus: (1) the exclusively Mediterranean genome A based on x = 6, (2) the cosmopolitan genome B based on x = 7 and (3) a genome C based on x = 7 and with a distribution in the Mediterranean and the Middle East. Both auto- and allopolyploidy of genomes B and C are suggested for the formation of tetraploids. The chromosomal divergence observed in Phalaris can be explained by the occurrence of dysploidy, the emergence of three different genomes, and the chromosome rearrangements accompanied by karyotype change and polyploidization. Mapping the recognized karyotypes on the existing phylogenetic tree suggests that genomes A and C are restricted to sections Phalaris and Bulbophalaris, respectively, while genome B occurs across all taxa with x = 7.

Peterson, Soreng RJ, Paszko B. Molecular phylogenetics of cool-season grasses in the subtribes Agrostidinae, Anthoxanthinae, Aveninae, Brizinae, Calothecinae, Koeleriinae and Phalaridinae (Poaceae, Pooideae, Poeae, Poeae chloroplast group 1)

Circumscriptions of and relationships among many genera and suprageneric taxa of the diverse grass tribe Poeae remain controversial. In an attempt to clarify these, we conducted phylogenetic analyses of >2400 new DNA sequences from two nuclear ribosomal regions (ITS, including internal transcribed spacers 1 and 2 and the 5.8S gene, and the 3'-end of the external transcribed spacer (ETS)) and five plastid regions (matK, trnL-trnF, atpF-atpH, psbK-psbI, psbA-rps19-trnH), and of more than 1000 new and previously published ITS sequences, focused particularly on Poeae chloroplast group 1 and including broad and increased species sampling compared to previous studies. Deep branches in the combined plastid and combined ITS+ETS trees are generally well resolved, the trees are congruent in most aspects, branch support across the trees is stronger than in trees based on only ITS and fewer plastid regions, and there is evidence of conflict between data partitions in some taxa. In plastid trees, a strongly supported clade corresponds to Poeae chloroplast group 1 and includes Agrostidinae p.p., Anthoxanthinae, Aveninae s.str., Brizinae, Koeleriinae (sometimes included in Aveninae s.l.), Phalaridinae and Torreyochloinae. In the ITS+ETS tree, a supported clade includes these same tribes as well as Sesleriinae and Scolochloinae. Aveninae s.str. and Sesleriinae are sister taxa and form a clade with Koeleriinae in the ITS+ETS tree whereas Aveninae s.str. and Koeleriinae form a clade and Sesleriinae is part of Poeae chloroplast group 2 in the plastid tree. All species of Trisetum are part of Koeleriinae, but the genus is polyphyletic. Koeleriinae is divided into two major subclades: one comprises Avellinia, Gaudinia, Koeleria, Rostraria, Trisetaria and Trisetum subg. Trisetum, and the other Calamagrostis/Deyeuxia p.p. (multiple species from Mexico to South America), Peyritschia, Leptophyllochloa, Sphenopholis, Trisetopsis and Trisetum subg. Deschampsioidea. Graphephorum, Trisetum cernuum, T. irazuense and T. macbridei fall in different clades of Koeleriinae in plastid vs. nuclear ribosomal trees, and are likely of hybrid origin. ITS and matK trees identify a third lineage of Koeleriinae corresponding to Trisetum subsect. Sibirica, and affinities of Lagurus ovatus with respect to Aveninae s.str. and Koeleriinae are incongruent in nuclear ribosomal and plastid trees, supporting recognition of Lagurus in its own subtribe. A large clade comprises taxa of Agrostidinae, Brizinae and Calothecinae, but neither Agrostidinae nor Calothecinae are monophyletic as currently circumscribed and affinities of Brizinae differ in plastid and nuclear ribosomal trees. Within this clade, one newly identified lineage comprises Calamagrostis coarctata, Dichelachne, Echinopogon (Agrostidinae p.p.) and Relchela (Calothecinae p.p.), and another comprises Chascolytrum (Calothecinae p.p.) and Deyeuxia effusa (Agrostidinae p.p.). Within Agrostidinae p.p., the type species of Deyeuxia and Calamagrostis s.str. are closely related, supporting classification of Deyeuxia as a synonym of Calamagrostis s.str. Furthermore, the two species of Ammophila are not sister taxa and are nested among different groups of Calamagrostis s.str., supporting their classification in Calamagrostis. Agrostis, Lachnagrostis and Polypogon form a clade and species of each are variously intermixed in plastid and nuclear ribosomal trees. Additionally, all but one species from South America classified in Deyeuxia sect. Stylagrostis resolve in Holcinae p.p. (Deschampsia). The current phylogenetic results support recognition of the latter species in Deschampsia, and we also demonstrate Scribneria is part of this clade. Moreover, Holcinae is not monophyletic in its current circumscription because Deschampsia does not form a clade with Holcus and Vahlodea, which are sister taxa. The results support recognition of Deschampsia in its own subtribe Aristaveninae. Substantial further changes to the classification of these grasses will be needed to produce generic circumscriptions consistent with phylogenetic evidence. The following 15 new combinations are made: Calamagrostis × calammophila, C. breviligulata, C. breviligulata subsp. champlainensis, C. × don-hensonii, Deschampsia aurea, D. bolanderi, D. chrysantha, D. chrysantha var. phalaroides, D. eminens, D. eminens var. fulva, D. eminens var. inclusa, D. hackelii, D. ovata, and D. ovata var. nivalis. D. podophora; the new name Deschampsia parodiana is proposed; the new subtribe Lagurinae is described; and a second-step lectotype is designated for the name Deyeuxia phalaroides.

Genome evolution in alpine oat-like grasses through homoploid hybridization and polyploidy

Hybridization and polyploidization can radically impact genome organization from sequence level to chromosome structure. As a result, often in response to environmental change and species isolation, the development of novel traits can arise and will tend to result in the formation of homoploid or polyploid hybrid species. In this study we focus on evidence of hybridization and polyploidization by ascertaining the species parentage of the endemic alpine Helictotrichon parlatorei group. This group comprises five taxa; the diploids H. parlatorei, Helictotrichon setaceum subsp. setaceum and subsp. petzense, their putative hybrid Helictotrichon ×krischae and the hexaploid Helictotrichon sempervirens. For molecular analyses, cloned nuclear Topoisomerase VI genes of H. sempervirens and H. ×krischae were sequenced and compared with sequences of the diploids to estimate the evolutionary history in this group. In addition, detailed chromosome studies were carried out including fluorescence in situ hybridization (FISH) with 5S and 45S ribosomal and satellite DNA probes, and fluorochrome staining with chromomycin and DAPI. Two distinct types of Topoisomerase VI sequences were identified. One of them (SET) occurs in both subspecies of H. setaceum, the other (PAR) in H. parlatorei. Both types were found in H. ×krischae and H. sempervirens Karyotypes of H. parlatorei and H. setaceum could be distinguished by chromosomes with a clearly differentiated banding pattern of ribosomal DNAs. Both patterns occurred in the hybrid H. ×krischae Hexaploid H. sempervirens shares karyotype features with diploid H. parlatorei, but lacks the expected chromosome characteristics of H. setaceum, possibly an example of beginning diploidization after polyploidization. The geographic origin of the putative parental species and their hybrids and the possible biogeographical spread through the Alps are discussed.

Mitosis and Interphase of the Highly Polyploid Palm Voanioalagerardii (2n = 606 ± 3)

The endemic, highly polyploid, monotypic Madagascan palm genus Voanioala (2n ≈ 606) was studied with regard to mitotic stages and interphase. Features of the cell cycle, morphology and sizes of metaphase chromosomes, fluorochrome banding patterns, and silver staining of NORs of such an extremely high polyploid organism are reported for the first time. On a whole, karyokinesis appears to be stable and efficient. A comparison with closely related palm taxa reveals that V. gerardii is 38-ploid, and comparison with the closely related genera Butia, Cocos (coconut) and Jubaea shows that Voanioala has lost ∼ 35% of its DNA amount subsequent to polyploidization and has suppressed between 74 and 88% of the original nucleolar organizers. About 10 active NORs are present in the nuclei. An auto- or allopolyploid origin of Voanioala is discussed with respect to currently available nuclear gene data. The biogeographic relations to Jubaeopsis, a closely related, monotypic, apparently likewise relict palm genus from eastern mainland South Africa are discussed. From a cytogenetic point of view, a common polyploid ancestor of both genera is most likely, but the available molecular phylogenetic data are not univocal.

Plastid phylogenomics of the cool-season grass subfamily: clarification of relationships among early-diverging tribes

Whole plastid genomes are being sequenced rapidly from across the green plant tree of life, and phylogenetic analyses of these are increasing resolution and support for relationships that have varied among or been unresolved in earlier single- and multi-gene studies. Pooideae, the cool-season grass lineage, is the largest of the 12 grass subfamilies and includes important temperate cereals, turf grasses and forage species. Although numerous studies of the phylogeny of the subfamily have been undertaken, relationships among some 'early-diverging' tribes conflict among studies, and some relationships among subtribes of Poeae have not yet been resolved. To address these issues, we newly sequenced 25 whole plastomes, which showed rearrangements typical of Poaceae. These plastomes represent 9 tribes and 11 subtribes of Pooideae, and were analysed with 20 existing plastomes for the subfamily. Maximum likelihood (ML), maximum parsimony (MP) and Bayesian inference (BI) robustly resolve most deep relationships in the subfamily. Complete plastome data provide increased nodal support compared with protein-coding data alone at nodes that are not maximally supported. Following the divergence of Brachyelytrum, Phaenospermateae, Brylkinieae-Meliceae and Ampelodesmeae-Stipeae are the successive sister groups of the rest of the subfamily. Ampelodesmeae are nested within Stipeae in the plastome trees, consistent with its hybrid origin between a phaenospermatoid and a stipoid grass (the maternal parent). The core Pooideae are strongly supported and include Brachypodieae, a Bromeae-Triticeae clade and Poeae. Within Poeae, a novel sister group relationship between Phalaridinae and Torreyochloinae is found, and the relative branching order of this clade and Aveninae, with respect to an Agrostidinae-Brizinae clade, are discordant between MP and ML/BI trees. Maximum likelihood and Bayesian analyses strongly support Airinae and Holcinae as the successive sister groups of a Dactylidinae-Loliinae clade.

Polyploidy and hybridization as main factors of speciation: complex reticulate evolution within the grass genus Helictochloa

To study the origin and evolution of naturally occurring polyploids, we performed phylogenetic analyses of nuclear ribosomal DNA spacers combined with molecular cytogenetics in 55 accessions of 27 taxa of the oat genus Helictochloa. A complex pattern of reticulate evolution was revealed with many diploid species and extensive polyploidy up to 20x. Altogether 11 groups of internal transcribed spacer (ITS) sequences can be distinguished. Sequences from 1-3 different ITS lineages were detected in polyploids. Cytogenetic data allow reconstruction of 8 basic monoploid chromosome sets. Six of these genomes occur in different combinations in the polyploid species. Two genomes are only found in diploids. Our sequence and karyological data highlight the occurrence of autopolyploidy and allopolyploidy, provide new information about the evolutionary history of taxa, and allow a more accurate systematic treatment of the concerned species. The geographical distribution of the 11 ITS lineages distinguished is highly structured and points to an origin of the genus in western Asia, presumably in grasslands like steppes or mountain steppes and meadows. The evolutionary basal lineages are of Asian, Minor Asian and east Mediterranean distribution and are present also in North America. The western and central parts of the Mediterranean and northern Europe harbor the modern lineages.

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.

Similar patterns of rDNA evolution in synthetic and recently formed natural populations of Tragopogon (Asteraceae) allotetraploids

BACKGROUND

Tragopogon mirus and T. miscellus are allotetraploids (2n = 24) that formed repeatedly during the past 80 years in eastern Washington and adjacent Idaho (USA) following the introduction of the diploids T. dubius, T. porrifolius, and T. pratensis (2n = 12) from Europe. In most natural populations of T. mirus and T. miscellus, there are far fewer 35S rRNA genes (rDNA) of T. dubius than there are of the other diploid parent (T. porrifolius or T. pratensis). We studied the inheritance of parental rDNA loci in allotetraploids resynthesized from diploid accessions. We investigate the dynamics and directionality of these rDNA losses, as well as the contribution of gene copy number variation in the parental diploids to rDNA variation in the derived tetraploids.

RESULTS

Using Southern blot hybridization and fluorescent in situ hybridization (FISH), we analyzed copy numbers and distribution of these highly reiterated genes in seven lines of synthetic T. mirus (110 individuals) and four lines of synthetic T. miscellus (71 individuals). Variation among diploid parents accounted for most of the observed gene imbalances detected in F1 hybrids but cannot explain frequent deviations from repeat additivity seen in the allotetraploid lines. Polyploid lineages involving the same diploid parents differed in rDNA genotype, indicating that conditions immediately following genome doubling are crucial for rDNA changes. About 19% of the resynthesized allotetraploid individuals had equal rDNA contributions from the diploid parents, 74% were skewed towards either T. porrifolius or T. pratensis-type units, and only 7% had more rDNA copies of T. dubius-origin compared to the other two parents. Similar genotype frequencies were observed among natural populations. Despite directional reduction of units, the additivity of 35S rDNA locus number is maintained in 82% of the synthetic lines and in all natural allotetraploids.

CONCLUSIONS

Uniparental reductions of homeologous rRNA gene copies occurred in both synthetic and natural populations of Tragopogon allopolyploids. The extent of these rDNA changes was generally higher in natural populations than in the synthetic lines. We hypothesize that locus-specific and chromosomal changes in early generations of allopolyploids may influence patterns of rDNA evolution in later generations.