High and uneven levels of 45S rDNA site-number variation across wild populations of a diploid plant genus (Anacyclus, Asteraceae)

Cytomolecular diversity among Vigna Savi (Leguminosae) subgenera

The genus Vigna (Leguminosae) comprises about 150 species grouped into five subgenera. The present study aimed to improve the understanding of karyotype diversity and evolution in Vigna, using new and previously published data through different cytogenetic and DNA content approaches. In the Vigna subgenera, we observed a random distribution of rDNA patterns. The 35S rDNA varied in position, from terminal to proximal, and in number, ranging from one (V. aconitifolia, V. subg. Ceratotropis) to seven pairs (V. unguiculata subsp. unguiculata, V. subg. Vigna). On the other hand, the number of 5S rDNA was conserved (one or two pairs), except for V. radiata (V. subg. Ceratotropis), which had three pairs. Genome size was relatively conserved within the genus, ranging from 1C = 0.43 to 0.70 pg in V. oblongifolia and V. unguiculata subsp. unguiculata, respectively, both belonging to V. subg. Vigna. However, we observed a positive correlation between DNA content and the number of 35S rDNA sites. In addition, data from chromosome-specific BAC-FISH suggest that the ancestral 35S rDNA locus is conserved on chromosome 6 within Vigna. Considering the rapid diversification in the number and position of rDNA sites, such conservation is surprising and suggests that additional sites may have spread out from this ancestral locus.

Unravelling plant diversification: Intraspecific genetic differentiation in hybridizing Anacyclus species in the western Mediterranean Basin

PREMISE

The interfertile species Anacyclus clavatus, A. homogamos, and A. valentinus represent a plant complex coexisting in large anthropic areas of the western Mediterranean Basin with phenotypically mixed populations exhibiting a great floral variation. The goal of this study was to estimate the genetic identity of each species, to infer the role of hybridization in the observed phenotypic diversity, and to explore the effect of climate on the geographic distribution of species and genetic clusters.

METHODS

We used eight nuclear microsatellites to genotype 585 individuals from 31 populations of three Anacyclus species for population genetic analyses by using clustering algorithms based on Bayesian models and ordination methods. In addition, we used ecological niche models and niche overlap analyses for both the species and genetic clusters. We used an expanded data set, including 721 individuals from 129 populations for ecological niche models of the genetic clusters.

RESULTS

We found a clear correspondence between species and genetic clusters, except for A. clavatus that included up to three genetic clusters. We detected individuals with admixed genetic ancestry in A. clavatus and in mixed populations. Ecological niche models predicted similar distributions for species and genetic clusters. For the two specific genetic clusters of A. clavatus, ecological niche models predicted remarkably different areas.

CONCLUSIONS

Gene flow between Anacyclus species likely explains phenotypic diversity in contact areas. In addition, we suggest that introgression could be involved in the origin of one of the two A. clavatus genetic clusters, which also showed ecological differentiation.

Comparative cytogenetics of three economically important Piper L. species from the Brazilian Amazon

The species Piper hispidinervum, Piper aduncum, and Piper affinis hispidinervum have essential oils with high levels of safrole, dillapiole, and sarisan, respectively. Safrole is important for pharmaceutical and chemical industries, while dillapiole and sarisan are promising compounds to control insects and fungi. These species are very similar morphologically and their taxonomy is controversial. Divergent hypotheses consider P. aduncum and P. hispidinervum either as a single species or as distinct taxa, while P. affinis hispidinervum is inferred to be a natural hybrid or a chemotype of P. hispidinervum. Delimiting the taxonomic boundaries would be helpful for germplasm conservation and breeding programs. This study aimed to undertake a detailed analysis of P. aduncum, P. hispidinervum, and P. affinis hispidinervum karyotype and rDNA sites. Genomic in situ hybridization (GISH) was used to establish genomic homology among species and to test the natural hybridization hypothesis for origin of P. affinis hispidinervum. Karyotype traits were similar for all three species: 2n = 26 small chromosomes, predominantly metacentric. All three species exhibited CMA+ bands on the secondary constriction of chromosome pair 4. A size-heteromorphic 35S rDNA site was co-localized with the CMA+ band. A 5S rDNA site was located in the proximal region of chromosome pair 7. The patterns of genomic hybridization revealed that the repetitive DNA fraction of the species is highly similar in terms of proportion of genome, sequence type, and distribution. Our findings did not allow us to differentiate the three species and point to the importance of deeper genomic studies to elucidate the taxonomic controversy.

Variation in the Number and Position of rDNA Loci Contributes to the Diversification and Speciation in Nigella (Ranunculaceae)

Nigella is a small genus belonging to the Ranunculaceae family which is presumably originated and distributed in Aegean and the adjacent Western-Irano-Turanian region. Comparative repeat analysis of N. sativa, N. damascena and N. bucharica was performed using low-pass Illumina genomic reads followed by karyotyping and FISH mapping of seven Nigella species using the in silico identified repeats and ribosomal DNA (rDNA) probes. High- and moderate-copy repeat sequences occupy 57.52, 59.01, and 64.73% of N. sativa, N. damascena and N. bucharica genomes, respectively. Roughly, half of the genomes are retrotransposons (class I transposons), while DNA transposons (class II transposons) contributed to only about 2% of the genomes. The analyzed Nigella species possess large genomes of about 7.4 to 12.4 Gbp/1C. Only two satellite repeats in N. sativa, one in N. damascena and four in N. bucharica were identified, which were mostly (peri)centromeric and represented about 1% of each genome. A high variation in number and position of 45S rDNA loci were found among Nigella species. Interestingly, in N. hispanica, each chromosome revealed at least one 45S rDNA site and one of them occurs in hemizygous condition. Based on the chromosome numbers, genome size and (peri)centromeric satellites, three karyotype groups were observed: Two with 2n = 2x = 12 and a karyotype formula of 10m + 2t (including N. sativa, N. arvensis, N. hispanica as the first group and N. damascena and N. orientalis as the second group) and a more distant group with 2n = 2x = 14 and a karyotype formula of 8m + 2st + 4t (including N. integrifolia and N. bucharica). These karyotype groups agreed with the phylogenetic analysis using ITS and rbcL sequences. We conclude that variation in (peri)centromeric sequences, number and localization of rDNA sites as well as chromosome number (dysploidy) are involved in the diversification of the genus Nigella.

Uneven Levels of 5S and 45S rDNA Site Number and Loci Variations across Wild Chrysanthemum Accessions

Ribosomal DNA (rDNA) is an excellent cytogenetic marker owing to its tandem arrangement and high copy numbers. However, comparative studies have focused more on the number of rDNA site variations within the Chrysanthemum genus, and studies on the types of rDNA sites with the same experimental procedures at the species levels are lacking. To further explore the number and types of rDNA site variations, we combined related data to draw ideograms of the rDNA sites of Chrysanthemum accessions using oligonucleotide fluorescence in situ hybridization (Oligo-FISH). Latent variations (such as polymorphisms of 45S rDNA sites and co-localized 5S-45S rDNA) also occurred among the investigated accessions. Meanwhile, a significant correlation was observed between the number of 5S rDNA sites and chromosome number. Additionally, the clumped and concentrated geographical distribution of different ploidy Chrysanthemum accessions may significantly promote the karyotype evolution. Based on the results above, we identified the formation mechanism of rDNA variations. Furthermore, these findings may provide a reliable method to examine the sites and number of rDNA variations among Chrysanthemum and its related accessions and allow researchers to further understand the evolutionary and phylogenetic relationships of the Chrysanthemum genus.

The Nuclear 35S rDNA World in Plant Systematics and Evolution: A Primer of Cautions and Common Misconceptions in Cytogenetic Studies

The ubiquitous presence of rRNA genes in nuclear, plastid, and mitochondrial genomes has provided an opportunity to use genomic markers to infer patterns of molecular and organismic evolution as well as to assess systematic issues throughout the tree of life. The number, size, location, and activity of the 35S rDNA cistrons in plant karyotypes have been used as conventional cytogenetic landmarks. Their scrutiny has been useful to infer patterns of chromosomal evolution and the data have been used as a proxy for assessing species discrimination, population differentiation and evolutionary relationships. The correct interpretation of rDNA markers in plant taxonomy and evolution is not free of drawbacks given the complexities derived from the lability of the genetic architecture, the diverse patterns of molecular change, and the fate and evolutionary dynamics of the rDNA units in hybrids and polyploid species. In addition, the terminology used by independent authors is somewhat vague, which often complicates comparisons. To date, no efforts have been reported addressing the potential problems and limitations involved in generating, utilizing, and interpreting the data from the 35S rDNA in cytogenetics. This review discusses the main technical and conceptual limitations of these rDNA markers obtained by cytological and karyological experimental work, in order to clarify biological and evolutionary inferences postulated in a systematic and phylogenetic context. Also, we provide clarification for some ambiguity and misconceptions in terminology usually found in published work that may help to improve the usage of the 35S ribosomal world in plant evolution.

Variation in Ribosomal DNA in the Genus Trifolium (Fabaceae)

The genus Trifolium L. is characterized by basic chromosome numbers 8, 7, 6, and 5. We conducted a genus-wide study of ribosomal DNA (rDNA) structure variability in diploids and polyploids to gain insight into evolutionary history. We used fluorescent in situ hybridization to newly investigate rDNA variation by number and position in 30 Trifolium species. Evolutionary history among species was examined using 85 available sequences of internal transcribed spacer 1 (ITS1) of 35S rDNA. In diploid species with ancestral basic chromosome number (x = 8), one pair of 5S and 26S rDNA in separate or adjacent positions on a pair of chromosomes was prevalent. Genomes of species with reduced basic chromosome numbers were characterized by increased number of signals determined on one pair of chromosomes or all chromosomes. Increased number of signals was observed also in diploids Trifolium alpestre and Trifolium microcephalum and in polyploids. Sequence alignment revealed ITS1 sequences with mostly single nucleotide polymorphisms, and ITS1 diversity was greater in diploids with reduced basic chromosome numbers compared to diploids with ancestral basic chromosome number (x = 8) and polyploids. Our results suggest the presence of one 5S rDNA site and one 26S rDNA site as an ancestral state.

Using target capture to address conservation challenges: Population-level tracking of a globally-traded herbal medicine

The promotion of responsible and sustainable trade in biological resources is widely proposed as one solution to mitigate current high levels of global biodiversity loss. Various molecular identification methods have been proposed as appropriate tools for monitoring global supply chains of commercialized animals and plants. Here, we demonstrate the efficacy of target capture genomic barcoding in identifying and establishing the geographic origin of samples traded as Anacyclus pyrethrum, a medicinal plant assessed as globally vulnerable in the IUCN Red List of Threatened Species. Samples collected from national and international supply chains were identified through target capture sequencing of 443 low-copy nuclear makers and compared to results derived from genome skimming of plastome and DNA barcoding of standard plastid regions and ITS. Both target capture and genome skimming provided approximately 3.4 million reads per sample, but target capture largely outperformed standard plant barcodes and entire plastid genome sequences. We were able to discern the geographical origin of Anacyclus samples collected in Moroccan, Indian and Sri Lankan markets, differentiating between plant materials originally harvested from diverse populations in Algeria and Morocco. Dropping costs of analysing samples enables the potential of target capture to routinely identify commercialized plant species and determine their geographic origin. It promises to play an important role in monitoring and regulation of plant species in trade, supporting biodiversity conservation efforts, and in ensuring that plant products are unadulterated, contributing to consumer protection.

Molecular Cytogenetics of Eurasian Species of the Genus Hedysarum L. (Fabaceae)

The systematic knowledge on the genus Hedysarum L. (Fabaceae: Hedysareae) is still incomplete. The species from the section Hedysarum are valuable forage and medicinal resources. For eight Hedysarum species, we constructed the integrated schematic map of their distribution within Eurasia based on currently available scattered data. For the first time, we performed cytogenomic characterization of twenty accessions covering eight species for evaluating genomic diversity and relationships within the section Hedysarum. Based on the intra- and interspecific variability of chromosomes bearing 45S and 5S rDNA clusters, four main karyotype groups were detected in the studied accessions: (1) H.arcticum, H. austrosibiricum, H. flavescens, H. hedysaroides, and H. theinum (one chromosome pair with 45S rDNA and one pair bearing 5S rDNA); (2) H. alpinum and one accession of H. hedysaroides (one chromosome pair with 45S rDNA and two pairs bearing 5S rDNA); (3) H. caucasicum (one chromosome pair with 45S rDNA and one chromosome pair bearing 5S rDNA and 45S rDNA); (4) H. neglectum (two pairs with 45S rDNA and one pair bearing 5S rDNA). The species-specific chromosomal markers detected in karyotypes of H. alpinum, H. caucasicum, and H. neglectum can be useful in taxonomic studies of this section.

Hybridization and cryptic speciation in the Iberian endemic plant genus Phalacrocarpum (Asteraceae-Anthemideae)

Understanding the role and impact of reticulation in phylogenetic inquiry has improved with extended use of high throughput sequencing data. Yet, due to the dynamism of genomes over evolutionary time, disentangling old hybridization events remains a serious challenge. Phalacrocarpum (DC.) Willk. is one of the 27 Iberian endemic plant genera, currently considered monotypic but including three subspecies. Its uncertain phylogenetic relationships within tribe Anthemideae (Asteraceae) point to an Early Miocene divergence from its sister group, and its persistent taxonomic instability has been proposed to be due to hybridization. We aim at understanding the evolutionary history of this genus using SNPs called from a genotyping-by-sequencing (GBS) analysis, Sanger sequences-from three plastid DNA regions (psbJ-petA, petB-petD, trnH-psbA) and the nuclear ribosomal ITS regions (cloned)-as well as leaf morphometric multivariate analysis. SNP data and Sanger sequences strongly support the unforeseen existence of a cryptic species in the eastern populations of P. oppositifolium subsp. anomalum. Broad molecular and morphometric patterns of variation found in conflictive populations from the Sanabria Valley region convincingly identify a recent previously undocumented hybrid zone. By contrast, evidence is less conclusive on relationships between subspecies hoffmannseggii, oppositifolium and a second conflictive group distributed along the Galician-Portuguese border (Orense massifs). Although genetic clustering analysis of SNP data suggests that the former subspecies was the maternal progenitor in hybridization events that gave rise to the other two groups, we found considerable uniqueness of ITS ribotypes and plastid haplotypes in them. This result, in the context of Pleistocene climatically-driven range shifts in NW Iberian Peninsula, can be due to periods of isolation, genetic bottlenecks and drift superimposed on old hybridization events. Our study confirms the idea that unravelling old hybridization events may be compromised by the suite of evolutionary processes accumulated subsequently, particularly in areas with a history of climatic instability.

Genome size variation at constant chromosome number is not correlated with repetitive DNA dynamism in Anacyclus (Asteraceae)

BACKGROUND AND AIMS

Changes in the amount of repetitive DNA (dispersed and tandem repeats) are considered the main contributors to genome size variation across plant species in the absence of polyploidy. However, the study of repeatome dynamism in groups showing contrasting genomic features and complex evolutionary histories is needed to determine whether other processes underlying genome size variation may have been overlooked. The main aim here was to elucidate which mechanism best explains genome size evolution in Anacyclus (Asteraceae).

METHODS

Using data from Illumina sequencing, we analysed the repetitive DNA in all species of Anacyclus, a genus with a reticulate evolutionary history, which displays significant genome size and karyotype diversity albeit presenting a stable chromosome number.

KEY RESULTS

By reconstructing ancestral genome size values, we inferred independent episodes of genome size expansions and contractions during the evolution of the genus. However, analysis of the repeatome revealed a similar DNA repeat composition across species, both qualitative and quantitative. Using comparative methods to study repeatome dynamics in the genus, we found no evidence for repeat activity causing genome size variation among species.

CONCLUSIONS

Our results, combined with previous cytogenetic data, suggest that genome size differences in Anacyclus are probably related to chromosome rearrangements involving losses or gains of chromosome fragments, possibly associated with homoploid hybridization. These could represent balanced rearrangements that do not disrupt gene dosage in merged genomes, for example via chromosome segment exchanges.

The Mendelian inheritance of gynomonoecy: insights from Anacyclus hybridizing species

PREMISE

Gynomonoecy is an infrequent sexual system in angiosperms, although widely represented within the Asteraceae family. Currently, the hypothesis of two nuclear loci controling gynomonoecy is the most accepted. However, the genic interactions are still uncertain. Anacyclus clavatus, A. homogamos, and A. valentinus differ in their sexual system and floral traits. Here, we investigate the inheritance of gynomonoecy in this model system to understand its prevalence in the family.

METHODS

We selected six natural populations (two per species) for intra- and interspecific experimental crosses, and generated a total of 1123 individuals from the F₁ generation, F₂ , and backcrosses for sexual system characterization. The frequency of gynomonoecy observed for each cross was tested to fit different possible hypotheses of genic interaction. Additionally, the breeding system and the degree of reproductive isolation between these species were assessed.

RESULTS

Complementary epistasis, in which two dominant alleles are required for trait expression, explained the frequencies of gynomonoecy observed across all generations. The heterozygosity inferred in Anacyclus valentinus, as well as its lower and variable seed set, is congruent with its hybrid origin.

CONCLUSIONS

In our model system gynomonoecy is controlled by complementary epistasis of two genes. A common origin of this sexual system in Asteraceae, in which genic duplications, mutations, and hybridization between lineages played a key role, is hypothesized whereas independent evolutionary pathways and possibly diverse underlying genetic factors are suggested for gynomonoecy expression in other angiosperm families.

Dissimilar molecular and morphological patterns in an introgressed peripheral population of a sand dune species (Armeria pungens, Plumbaginaceae)

Introgression is a poorly understood evolutionary outcome of hybridisation because it may remain largely undetected whenever it involves the transfer of small parts of the genome from one species to another. Aiming to understand the early stages of this process, a putative case from the southernmost border of the Armeria pungens range from its congener A. macrophylla is revisited following the discovery of a subpopulation that does not show phenotypic signs of introgression and resembles typical A. pungens. We analysed morphometrics, nuclear ribosomal DNA ITS and plastid DNA (trnL-trnF) sequences, genome size, 45S and 5S rDNA loci-FISH data and nrDNA IGS sequences. Within the study site, most individuals match morphologies of either of the two hybridising species, particularly the new subpopulation, with intermediate phenotypes being scarce. This pattern does not fully fit molecular evidence revealing two ITS ribotypes co-occurring intragenomically in most plants from the study site and one single plastid haplotype. Genome size and structural features of the IGS sequences both indicate that A. pungens from the study site is genetically more similar to its sympatric congener than to the remainder of its conspecifics. Introgression of A. macrophylla into A. pungens and plastid capture explain all the evidence analysed. However, important features to understand the origin and fate of the introgressed population, such as the degree and direction of introgression, which are important for understanding early stages of hybridisation in plants with low reproductive barriers, should be addressed with new data.

Comparative Dissection of Three Giant Genomes: Allium cepa, Allium sativum, and Allium ursinum

Knowledge of the fascinating world of DNA repeats is continuously being enriched by newly identified elements and their hypothetical or well-established biological relevance. Genomic approaches can be used for comparative studies of major repeats in any group of genomes, regardless of their size and complexity. Such studies are particularly fruitful in large genomes, and useful mainly in crop plants where they provide a rich source of molecular markers or information on indispensable genomic components (e.g., telomeres, centromeres, or ribosomal RNA genes). Surprisingly, in Allium species, a comprehensive comparative study of repeats is lacking. Here we provide such a study of two economically important species, Allium cepa (onion), and A. sativum (garlic), and their distantly related A. ursinum (wild garlic). We present an overview and classification of major repeats in these species and have paid specific attention to sequence conservation and copy numbers of major representatives in each type of repeat, including retrotransposons, rDNA, or newly identified satellite sequences. Prevailing repeats in all three studied species belonged to Ty3/gypsy elements, however they significantly diverged and we did not detect them in common clusters in comparative analysis. Actually, only a low number of clusters was shared by all three species. Such conserved repeats were for example 5S and 45S rDNA genes and surprisingly a specific and quite rare Ty1/copia lineage. Species-specific long satellites were found mainly in A. cepa and A. sativum. We also show in situ localization of selected repeats that could potentially be applicable as chromosomal markers, e.g., in interspecific breeding.

Healthy Roots and Leaves: Comparative Genome Structure of Horseradish and Watercress

Horseradish (Armoracia rusticana) and watercress (Nasturtium officinale) are economically important cruciferous vegetable species with limited genomic resources. We used comparative chromosome painting to identify the extent of chromosomal collinearity between horseradish and watercress, and to reconstruct the origin and evolution of the two tetraploid genomes (2n = 4x = 32). Our results show that horseradish and watercress genomes originated from a common ancestral (n = 8) genome, structurally resembling the Ancestral Crucifer Karyotype (n = 8), which, however, contained two unique translocation chromosomes (AK6/8 and AK8/6). Except for a 2.4-Mb unequal chromosome translocation in watercress, both genomes are structurally identical. The structural similarity of the two parental subgenomes might suggest an autotetraploid origin of horseradish and watercress genomes. The subgenome stasis, apart from the single-chromosome translocation, indicates that homeologous recombination played a limited role in postpolyploid evolution in both tetraploid genomes. The octoploid genome of one-rowed watercress (N. microphyllum, 2n = 8x = 64), structurally mirroring the tetraploid horseradish and watercress genomes, originated via autopolyploidization from the immediate tetraploid predecessor of watercress or hybridization between this and another now-extinct tetraploid Nasturtium species. These comparative cytogenomic maps in horseradish and watercress represent a first stepping stone for future whole-genome sequencing efforts and genetic improvement of both crop species.

Karyotype heterogeneity in Philodendron s.l. (Araceae) revealed by chromosome mapping of rDNA loci

Philodendron s.l. (Araceae) has been recently focus of taxonomic and phylogenetic studies, but karyotypic data are limited to chromosome numbers and a few published genome sizes. In this work, karyotypes of 34 species of Philodendron s.l. (29 species of Philodendron and five of Thaumatophyllum), ranging from 2n = 28 to 36 chromosomes, were analyzed by fluorescence in situ hybridization (FISH) with rDNA and telomeric probes, aiming to understand the evolution of the karyotype diversity of the group. Philodendron presented a high number variation of 35S rDNA, ranging from two to 16 sites, which were mostly in the terminal region of the short arms, with nine species presenting heteromorphisms. In the case of Thaumatophyllum species, we observed a considerably lower variation, which ranged from two to four terminal sites. The distribution of the 5S rDNA clusters was more conserved, with two sites for most species, being preferably located interstitially in the long chromosome arms. For the telomeric probe, while exclusively terminal sites were observed for P. giganteum (2n = 30) chromosomes, P. callosum (2n = 28) presented an interstitial distribution associated with satellite DNA. rDNA sites of the analyzed species of Philodendron s.l. species were randomly distributed considering the phylogenetic context, probably due to rapid evolution and great diversity of these genomes. The observed heteromorphisms suggest the accumulation of repetitive DNA in the genomes of some species and the occurrence of chromosomal rearrangements along the karyotype evolution of the group.

Inter- and intraspecific hypervariability in interstitial telomeric-like repeats (TTTAGGG)n in Anacyclus (Asteraceae)

Background and Aims

Interstitial telomeric repeat (ITR) sites, consisting of tandem repeats of telomeric motifs localized at intrachromosomal sites, have been reported in a few unrelated organisms including plants. However, the causes for the occurrence of ITRs outside of the chromosomal termini are not fully understood. One possible explanation are the chromosomal rearrangements involving telomeric sites, which could also affect the location of other structural genome elements, such as the 45S rDNA. Taking advantage of the high dynamism in 45S rDNA loci previously found in Anacyclus (Asteraceae, Anthemideae), the occurrence and patterns of variation of ITRs were explored in this genus with the aim of finding common underlying causes.

Methods

In total, 132 individuals from 44 populations of nine species were analysed by fluorescence in situ hybridization using an Arabidopsis-type telomeric sequence as a probe.

Key results

Variable presence of ITR sites was detected in six out of nine species of Anacyclus, ranging from two to 45 sites and showing contrasting chromosomal locations and a differential presence of the ITR site on homologous chromosome pairs. At the intraspecific level, the ranges were as large as 0-12 ITR sites. Although only 26 % of the total observed ITR sites were localized in chromosomes bearing 45S rDNA loci, all cases of interstitial 45S rDNA reported in a previous work co-occurred with ITRs in close proximity in the same chromosome arms.

Conclusions

High levels of ITR polymorphism within a single species have not been previously reported in plants and suggest that this pattern might have been overlooked due to insufficient sampling. Although ancient Robertsonian translocations or the amplification of terminal 45S rDNA sites cannot, on their own, explain all of the levels of variability in ITRs reported here, there are suggestions that they may have been involved in the evolutionary history of this genus or its ancestors in Anthemideae.