Concerted evolution of 18-5.8-26S rDNA repeats in Nicotiana allotetraploids

Evolution of the Araliaceae family involved rapid diversification of the Asian Palmate group and Hydrocotyle specific mutational pressure

The Araliaceae contain many valuable species in medicinal and industrial aspects. We performed intensive phylogenomics using the plastid genome (plastome) and 45S nuclear ribosomal DNA sequences. A total of 66 plastome sequences were used, 13 of which were newly assembled in this study, 12 from new sequences, and one from existing data. While Araliaceae plastomes showed conserved genome structure, phylogenetic reconstructions based on four different plastome datasets revealed phylogenetic discordance within the Asian Palmate group. The divergence time estimation revealed that splits in two Araliaceae subfamilies and the clades exhibiting phylogenetic discordances in the Asian Palmate group occurred at two climatic optima, suggesting that global warming events triggered species divergence, particularly the rapid diversification of the Asian Palmate group during the Middle Miocene. Nucleotide substitution analyses indicated that the Hydrocotyloideae plastomes have undergone accelerated AT-biased mutations (C-to-T transitions) compared with the Aralioideae plastomes, and the acceleration may occur in their mitochondrial and nuclear genomes as well. This implies that members of the genus Hydrocotyle, the only aquatic plants in the Araliaceae, have experienced a distinct evolutionary history from the other species. We also discussed the intercontinental disjunction in the genus Panax and proposed a hypothesis to complement the previously proposed hypothesis. Our results provide the evolutionary trajectory of Araliaceae and advance our current understanding of the evolution of Araliaceae species.

Analysis of the complete plastomes and nuclear ribosomal DNAs from Euonymus hamiltonianus and its relatives sheds light on their diversity and evolution

Euonymus hamiltonianus and its relatives (Celastraceae family) are used for ornamental and medicinal purposes. However, species identification in Euonymus is difficult due to their morphological diversity. Using plastid genome (plastome) data, we attempt to reveal phylogenetic relationship among Euonymus species and develop useful markers for molecular identification. We assembled the plastome and nuclear ribosomal DNA (nrDNA) sequences from five Euonymus lines collected from South Korea: three Euonymus hamiltonianus accessions, E. europaeus, and E. japonicus. We conducted an in-depth comparative analysis using ten plastomes, including other publicly available plastome data for this genus. The genome structures, gene contents, and gene orders were similar in all Euonymus plastomes in this study. Analysis of nucleotide diversity revealed six divergence hotspots in their plastomes. We identified 339 single nucleotide polymorphisms and 293 insertion or deletions among the four E. hamiltonianus plastomes, pointing to abundant diversity even within the same species. Among 77 commonly shared genes, 9 and 33 were identified as conserved genes in the genus Euonymus and E. hamiltonianus, respectively. Phylogenetic analysis based on plastome and nrDNA sequences revealed the overall consensus and relationships between plastomes and nrDNAs. Finally, we developed six barcoding markers and successfully applied them to 31 E. hamiltonianus lines collected from South Korea. Our findings provide the molecular basis for the classification and molecular taxonomic criteria for the genus Euonymus (at least in Korea), which should aid in more objective classification within this genus. Moreover, the newly developed markers will be useful for understanding the species delimitation of E. hamiltonianus and closely related species.

Intact ribosomal DNA arrays of Potentilla origin detected in Erythronium nucleus suggest recent eudicot-to-monocot horizontal transfer

During our initial phylogenetic study of the monocot genus Erythronium (Liliaceae), we observed peculiar eudicot-type internal transcribed spacer (ITS) sequences in a dataset derived from genomic DNA of Erythronium dens-canis. This raised the possibility of horizontal transfer of a eudicot alien ribosomal DNA (rDNA) into the Erythronium genome. In this work we aimed to support this hypothesis by carrying out genomic, molecular, and cytogenetic analyses. Genome skimming coupled by PacBio HiFi sequencing of a bacterial artificial chromosome clone derived from flow-sorted nuclei was used to characterise the alien 45S rDNA. Integration of alien rDNA in the recipient genome was further proved by Southern blotting and fluorescence in situ hybridization using specific probes. Alien rDNA, nested among Potentilla species in phylogenetic analysis, likely entered the Erythronium lineage in the common ancestor of E. dens-canis and E. caucasicum. Transferred eudicot-type rDNA preserved its tandemly arrayed feature on a single chromosome and was found to be transcribed in the monocot host, albeit much less efficiently than the native counterpart. This study adds a new example to the rarely documented nuclear-to-nuclear jumps of DNA between eudicots and monocots while holding the scientific community continually in suspense about the mode of DNA transfer.

Global Patterns of subgenome evolution in organelle-targeted genes of six allotetraploid angiosperms

Whole-genome duplications (WGDs) are a prominent process of diversification in eukaryotes. The genetic and evolutionary forces that WGD imposes on cytoplasmic genomes are not well understood, despite the central role that cytonuclear interactions play in eukaryotic function and fitness. Cellular respiration and photosynthesis depend on successful interaction between the 3,000+ nuclear-encoded proteins destined for the mitochondria or plastids and the gene products of cytoplasmic genomes in multi-subunit complexes such as OXPHOS, organellar ribosomes, Photosystems I and II, and Rubisco. Allopolyploids are thus faced with the critical task of coordinating interactions between the nuclear and cytoplasmic genes that were inherited from different species. Because the cytoplasmic genomes share a more recent history of common descent with the maternal nuclear subgenome than the paternal subgenome, evolutionary “mismatches” between the paternal subgenome and the cytoplasmic genomes in allopolyploids might lead to the accelerated rates of evolution in the paternal homoeologs of allopolyploids, either through relaxed purifying selection or strong directional selection to rectify these mismatches. We report evidence from six independently formed allotetraploids that the subgenomes exhibit unequal rates of protein-sequence evolution, but we found no evidence that cytonuclear incompatibilities result in altered evolutionary trajectories of the paternal homoeologs of organelle-targeted genes. The analyses of gene content revealed mixed evidence for whether the organelle-targeted genes are lost more rapidly than the non-organelle-targeted genes. Together, these global analyses provide insights into the complex evolutionary dynamics of allopolyploids, showing that the allopolyploid subgenomes have separate evolutionary trajectories despite sharing the same nucleus, generation time, and ecological context.

Tracing the Evolution of the Angiosperm Genome from the Cytogenetic Point of View

Cytogenetics constitutes a branch of genetics that is focused on the cellular components, especially chromosomes, in relation to heredity and genome structure, function and evolution. The use of modern cytogenetic approaches and the latest microscopes with image acquisition and processing systems enables the simultaneous two- or three-dimensional, multicolour visualisation of both single-copy and highly-repetitive sequences in the plant genome. The data that is gathered using the cytogenetic methods in the phylogenetic background enable tracing the evolution of the plant genome that involve changes in: (i) genome sizes; (ii) chromosome numbers and morphology; (iii) the content of repetitive sequences and (iv) ploidy level. Modern cytogenetic approaches such as FISH using chromosome- and genome-specific probes have been widely used in studies of the evolution of diploids and the consequences of polyploidy. Nowadays, modern cytogenetics complements analyses in other fields of cell biology and constitutes the linkage between genetics, molecular biology and genomics.

Chromosome Evolution in the Family Solanaceae

This review summarizes and discusses the knowledge of cytogenetics in Solanaceae, the tomato family, its current applications, and prospects for making progress in fundamental systematic botany and plant evolution. We compile information on basic chromosome features (number, size, morphology) and molecular cytogenetics (chromosome banding and rDNA patterns). These data were mapped onto the Solanaceae family tree to better visualize the changes in chromosome features and evaluate them in a phylogenetic context. We conclude that chromosomal features are important in understanding the evolution of the family, especially in delimiting clades, and therefore it is necessary to continue producing this type of data. The potential for future applications in plant biology is outlined. Finally, we provide insights into understanding the mechanisms underlying Solanaceae's diversification that could substantially contribute to developing new approaches for future research.

Personal Perspectives on Plant Ribosomal RNA Genes Research: From Precursor-rRNA to Molecular Evolution

The history of rDNA research started almost 90 years ago when the geneticist, Barbara McClintock observed that in interphase nuclei of maize the nucleolus was formed in association with a specific region normally located near the end of a chromosome, which she called the nucleolar organizer region (NOR). Cytologists in the twentieth century recognized the nucleolus as a common structure in all eukaryotic cells, using both light and electron microscopy and biochemical and genetic studies identified ribosomes as the subcellular sites of protein synthesis. In the mid- to late 1960s, the synthesis of nuclear-encoded rRNA was the only system in multicellular organisms where transcripts of known function could be isolated, and their synthesis and processing could be studied. Cytogenetic observations of NOR regions with altered structure in plant interspecific hybrids and detailed knowledge of structure and function of rDNA were prerequisites for studies of nucleolar dominance, epistatic interactions of rDNA loci, and epigenetic silencing. In this article, we focus on the early rDNA research in plants, performed mainly at the dawn of molecular biology in the 60 to 80-ties of the last century which presented a prequel to the modern genomic era. We discuss - from a personal view - the topics such as synthesis of rRNA precursor (35S pre-rRNA in plants), processing, and the organization of 35S and 5S rDNA. Cloning and sequencing led to the observation that the transcribed and processed regions of the rRNA genes vary enormously, even between populations and species, in comparison with the more conserved regions coding for the mature rRNAs. Epigenetic phenomena and the impact of hybridization and allopolyploidy on rDNA expression and homogenization are discussed. This historical view of scientific progress and achievements sets the scene for the other articles highlighting the immense progress in rDNA research published in this special issue of Frontiers in Plant Science on "Molecular organization, evolution, and function of ribosomal DNA."

Power and Weakness of Repetition - Evaluating the Phylogenetic Signal From Repeatomes in the Family Rosaceae With Two Case Studies From Genera Prone to Polyploidy and Hybridization (Rosa and Fragaria)

Plant genomes consist, to a considerable extent, of non-coding repetitive DNA. Several studies showed that phylogenetic signals can be extracted from such repeatome data by using among-species dissimilarities from the RepeatExplorer2 pipeline as distance measures. Here, we advanced this approach by adjusting the read input for comparative clustering indirectly proportional to genome size and by summarizing all clusters into a main distance matrix subjected to Neighbor Joining algorithms and Principal Coordinate Analyses. Thus, our multivariate statistical method works as a "repeatomic fingerprint," and we proved its power and limitations by exemplarily applying it to the family Rosaceae at intrafamilial and, in the genera Fragaria and Rosa, at the intrageneric level. Since both taxa are prone to hybridization events, we wanted to show whether repeatome data are suitable to unravel the origin of natural and synthetic hybrids. In addition, we compared the results based on complete repeatomes with those from ribosomal DNA clusters only, because they represent one of the most widely used barcoding markers. Our results demonstrated that repeatome data contained a clear phylogenetic signal supporting the current subfamilial classification within Rosaceae. Accordingly, the well-accepted major evolutionary lineages within Fragaria were distinguished, and hybrids showed intermediate positions between parental species in data sets retrieved from both complete repeatomes and rDNA clusters. Within the taxonomically more complicated and particularly frequently hybridizing genus Rosa, we detected rather weak phylogenetic signals but surprisingly found a geographic pattern at a population scale. In sum, our method revealed promising results at larger taxonomic scales as well as within taxa with manageable levels of reticulation, but success remained rather taxon specific. Since repeatomes can be technically easy and comparably inexpensively retrieved even from samples of rather poor DNA quality, our phylogenomic method serves as a valuable alternative when high-quality genomes are unavailable, for example, in the case of old museum specimens.

Analysis of Codon Usage Patterns of Six Sequenced Brachypodium distachyon Lines Reveals a Declining CG Skew of the CDSs from the 5'-ends to the 3'-ends

Brachypodium distachyon, a new monocotyledonous model plant, has received wide attention in biological research due to its small genome and numerous genetic resources. Codon usage bias is an important feature of genes and genomes, and it can be used in transgenic and evolutionary studies. In this study, the nucleotide compositions and patterns of codon usage bias were calculated using Codon W. Additionally, an ENC plot, Parity rule 2 and correspondence analyses were used to explore the major factors influencing codon usage bias patterns. The numbers of hydrogen bonds and skews were used to analyze the GC trend in the 5'-ends of the coding sequences. The results showed that minor differences in the codon usage bias patterns were revealed by the ENC plot, Parity rule 2 and correspondence analyses. The analyses of the CG-skew and the number of hydrogen bonds showed a declining trend in the number of cytosines at the 5'-ends of the CDSs (from the 5'-ends to the 3'-ends), indicating that GC may play a major role in codon usage bias. In addition, our results laid a foundation for the study of codon usage bias patterns in Brachypodium genus and suggested that the GC plays a major role in determining these patterns.

Chromosome inheritance and meiotic stability in allopolyploid Brassica napus

Homoeologous recombination, aneuploidy, and other genetic changes are common in resynthesized allopolyploid Brassica napus. In contrast, the chromosomes of cultivars have long been considered to be meiotically stable. To gain a better understanding of the underlying mechanisms leading to stabilization in the allopolyploid, the behavior of chromosomes during meiosis can be compared by unambiguous chromosome identification between resynthesized and natural B. napus. Compared with natural B. napus, resynthesized lines show high rates of nonhomologous centromere association, homoeologous recombination leading to translocation, homoeologous chromosome replacement, and association and breakage of 45S rDNA loci. In both natural and resynthesized B. napus, we observed low rates of univalents, A–C bivalents, and early sister chromatid separations. Reciprocal homoeologous chromosome exchanges and double reductions were photographed for the first time in meiotic telophase I. Meiotic errors were non-uniformly distributed across the genome in resynthesized B. napus, and in particular homoeologs sharing synteny along their entire length exhibited multivalents at diakinesis and polysomic inheritance at telophase I. Natural B. napus appeared to resolve meiotic errors mainly by suppressing homoeologous pairing, resolving nonhomologous centromere associations and 45S rDNA associations before diakinesis, and reducing homoeologous cross-overs.

Plastid genomics of Nicotiana (Solanaceae): insights into molecular evolution, positive selection and the origin of the maternal genome of Aztec tobacco (Nicotiana rustica)

Species of the genus Nicotiana (Solanaceae), commonly referred to as tobacco plants, are often cultivated as non-food crops and garden ornamentals. In addition to the worldwide production of tobacco leaves, they are also used as evolutionary model systems due to their complex development history tangled by polyploidy and hybridization. Here, we assembled the plastid genomes of five tobacco species: N. knightiana, N. rustica, N. paniculata, N. obtusifolia and N. glauca. De novo assembled tobacco plastid genomes had the typical quadripartite structure, consisting of a pair of inverted repeat (IR) regions (25,323-25,369 bp each) separated by a large single-copy (LSC) region (86,510-86,716 bp) and a small single-copy (SSC) region (18,441-18,555 bp). Comparative analyses of Nicotiana plastid genomes with currently available Solanaceae genome sequences showed similar GC and gene content, codon usage, simple sequence and oligonucleotide repeats, RNA editing sites, and substitutions. We identified 20 highly polymorphic regions, mostly belonging to intergenic spacer regions (IGS), which could be suitable for the development of robust and cost-effective markers for inferring the phylogeny of the genus Nicotiana and family Solanaceae. Our comparative plastid genome analysis revealed that the maternal parent of the tetraploid N. rustica was the common ancestor of N. paniculata and N. knightiana, and the later species is more closely related to N. rustica. Relaxed molecular clock analyses estimated the speciation event between N. rustica and N. knightiana appeared 0.56 Ma (HPD 0.65-0.46). Biogeographical analysis supported a south-to-north range expansion and diversification for N. rustica and related species, where N. undulata and N. paniculata evolved in North/Central Peru, while N. rustica developed in Southern Peru and separated from N. knightiana, which adapted to the Southern coastal climatic regimes. We further inspected selective pressure on protein-coding genes among tobacco species to determine if this adaptation process affected the evolution of plastid genes. These analyses indicate that four genes involved in different plastid functions, including DNA replication (rpoA) and photosynthesis (atpB, ndhD and ndhF), came under positive selective pressure as a result of specific environmental conditions. Genetic mutations in these genes might have contributed to better survival and superior adaptations during the evolutionary history of tobacco species.

Chromosome instabilities in resynthesized Brassica napus revealed by FISH

Brassica napus is an allopolyploid plant, derived from spontaneous hybridization between Brassica rapa and Brassica oleracea. Intensive breeding has led to a significant reduction in genetic and phenotypic diversity within this species. Newly resynthesized hybrids from progenitor species may restore some diversity in B. napus, but they often are chromosomally and phenotypically unstable. Using fluorescence in situ hybridization, we tested chromosome constitutions in a range of new allopolyploids resynthesized from various parental species. A majority of these allopolyploids were euploid, with the expected chromosome numbers and constitutions, but deviations were also identified. We detected a low level of intergenomic rearrangements in analyzed hybrids and a high level of changes in rDNA loci. Our study revealed a significant effect of maternal cross combination on loss of 35S rDNA loci, especially when B. rapa was the maternal parent. The studied lines were characterized by diversified of pollen viability. In the analyzed hybrids, the erucic acid level in the seed oil ranged from 0 to 43.4% and total glucosinolate content in seeds ranged from 24.3 to 119.2 μmol g⁻¹. Our study shows that cytogenetic analysis of B. napus resynthesized hybrids would be useful in breeding for the selection of lines with important agricultural characters and genetically stable stock seed production.

Molecular cytogenetic characterisation of Elytrigia ×mucronata, a natural hybrid of E. intermedia and E. repens (Triticeae, Poaceae)

BACKGROUND

Interspecific hybridisation resulting in polyploidy is one of the major driving forces in plant evolution. Here, we present data from the molecular cytogenetic analysis of three cytotypes of Elytrigia ×mucronata using sequential fluorescence (5S rDNA, 18S rDNA and pSc119.2 probes) and genomic in situ hybridisation (four genomic probes of diploid taxa, i.e., Aegilops, Dasypyrum, Hordeum and Pseudoroegneria).

RESULTS

The concurrent presence of Hordeum (descended from E. repens) and Dasypyrum + Aegilops (descended from E. intermedia) chromosome sets in all cytotypes of E. ×mucronata confirmed the assumed hybrid origin of the analysed plants. The following different genomic constitutions were observed for E. ×mucronata. Hexaploid plants exhibited three chromosome sets from Pseudoroegneria and one chromosome set each from Aegilops, Hordeum and Dasypyrum. Heptaploid plants harboured the six chromosome sets of the hexaploid plants and an additional Pseudoroegneria chromosome set. Nonaploid cytotypes differed in their genomic constitutions, reflecting different origins through the fusion of reduced and unreduced gametes. The hybridisation patterns of repetitive sequences (5S rDNA, 18S rDNA, and pSc119.2) in E. ×mucronata varied between and within cytotypes. Chromosome alterations that were not identified in the parental species were found in both heptaploid and some nonaploid plants.

CONCLUSIONS

The results confirmed that both homoploid hybridisation and heteroploid hybridisation that lead to the coexistence of four different haplomes within single hybrid genomes occur in Elytrigia allopolyploids. The chromosomal alterations observed in both heptaploid and some nonaploid plants indicated that genome restructuring occurs during and/or after the hybrids arose. Moreover, a specific chromosomal translocation detected in one of the nonaploids indicated that it was not a primary hybrid. Therefore, at least some of the hybrids are fertile. Hybridisation in Triticeae allopolyploids clearly and significantly contributes to genomic diversity. Different combinations of parental haplomes coupled with chromosomal alterations may result in the establishment of unique lineages, thus providing raw material for selection.

Ribosomal DNA loci derived from Brachypodium stacei are switched off for major parts of the life cycle of Brachypodium hybridum

Nucleolar dominance is an epigenetic phenomenon that occurs in some plant and animal allopolyploids and hybrids, whereby only one ancestral set of 35S rRNA genes retains the ability to form the nucleolus while the rDNA loci derived from the other progenitor are transcriptionally silenced. There is substantial evidence that nucleolar dominance is regulated developmentally. This study focuses upon the establishment and/or maintenance of nucleolar dominance during different stages of development in the model grass allotetraploid Brachypodium hybridum. Fluorescence in situ hybridization with a 25S rDNA probe to cells in three-dimensional cytogenetic preparations showed that nucleolar dominance is present not only in root meristematic and differentiated cells of this species, but also in male meiocytes at prophase I, tetrads of microspores, and different embryonic tissues. The inactive state of Brachypodium stacei-originated rDNA loci was confirmed by silver staining. Only B. distachyon-derived 35S rDNA loci formed nucleoli in the aforementioned tissues, whereas B. stacei-like loci remained highly condensed and thus transcriptionally suppressed. The establishment of nucleolar dominance during earlier stages of B. hybridum embryo development cannot be ruled out. However, we propose that gradual pseudogenization of B. stacei-like loci in the evolution of the allotetraploid seems to be more likely.

Parental transposable element loads influence their dynamics in young Nicotiana hybrids and allotetraploids

The genomic shock hypothesis suggests that allopolyploidy is associated with genome changes driven by transposable elements, as a response to imbalances between parental insertion loads. To explore this hypothesis, we compared three allotetraploids, Nicotiana arentsii, N. rustica and N. tabacum, which arose over comparable time frames from hybridisation between increasingly divergent diploid species. We used sequence-specific amplification polymorphism (SSAP) to compare the dynamics of six transposable elements in these allopolyploids, their diploid progenitors and in corresponding synthetic hybrids. We show that element-specific dynamics in young Nicotiana allopolyploids reflect their dynamics in diploid progenitors. Transposable element mobilisation is not concomitant with immediate genome merger, but occurs within the first generations of allopolyploid formation. In natural allopolyploids, such mobilisations correlate with imbalances in the repeat profile of the parental species, which increases with their genetic divergence. Other restructuring leading to locus loss is immediate, nonrandom and targeted at specific subgenomes, independently of cross orientation. The correlation between transposable element mobilisation in allopolyploids and quantitative imbalances in parental transposable element loads supports the genome shock hypothesis proposed by McClintock.

The fate of ribosomal RNA genes in spontaneous polyploid dogrose hybrids [Rosa L. sect. Caninae (DC.) Ser.] exhibiting non-symmetrical meiosis

Dogroses represent an exceptional system for studying the effects of genome doubling and hybridization: their asymmetrical meiosis enables recombination in bi-parentally inherited chromosomes but prevents it in maternally inherited ones. We employed fluorescent in situ hybridization, genome skimming, amplicon sequencing of genomic and cDNA as well as conventional cloning of nuclear ribosomal DNA in two phylogenetically distinct pentaploid (2n = 5x = 35) species, Rosa canina and Rosa inodora, and their naturally occurring reciprocal hybrids, Rosa dumalis (5x) and Rosa agrestis (5x, 6x). Both progenitor species differed in composition, meiotic behaviour and expression of rDNA loci: R. canina (five 18S and 5-8 5S loci) was dominated by the Canina ribotypes, but R. inodora (four 18S loci and 7-8 5S loci) by the Rubiginosa ribotype. The co-localized 5S/18S loci occurred on either bivalent-forming (R. canina) or univalent-forming (R. inodora) chromosomes. Ribosomal DNA loci were additively inherited; however, the Canina ribotypes were dominantly expressed, even in genotypes with relatively low copy number of these genes. Moreover, we observed rDNA homogenization towards the paternally transmitted Canina ribotype in 6x R. agrestis. The here-observed variation in arrangement and composition of rDNA types between R. canina and R. inodora suggests the involvement of different genomes in bivalent formation. This results supports the hypothesis that the asymmetrical meiosis arose at least twice by independent ancient hybridization events.

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.

Ecological Roles and Biological Activities of Specialized Metabolites from the Genus Nicotiana

Species of Nicotiana grow naturally in different parts of the world and have long been used both medicinally and recreationally by human societies. More recently in our history, Nicotiana tabacum has attracted interest as one of the most economically important industrial crops. Nicotiana species are frequently investigated for their bioactive natural products, and the ecological role of their specialized metabolites in responses to abiotic stress or biotic stress factors like pathogens and herbivores. The interest of tobacco companies in genetic information as well as the success of a few wild tobacco species as experimental model organisms have resulted in growing knowledge about the molecular biology and ecology of these plants and functional studies of the plant's natural products. Although a large number of reviews and books on biologically active natural products already exists, mostly from N. tabacum, we focus our attention on the ecological roles and biological activity of natural products, versus products from cured and processed material, in this Review. The studied compounds include alkaloids, aromatic compounds, flavonoids, volatiles, sesquiterpenoids, diterpenes alcohols, and sugar esters from trichomes of the plants, and recently characterized acyclic hydroxygeranyllinalool diterpene glycosides (HGL-DTGs). In this Review (1800s-2017), we describe the above-mentioned classes of natural products, emphasizing their biological activities and functions as they have been determined either in bioassay-guided purification approaches or in bioassays with plants in which the expression of specific biosynthetic genes has been genetically manipulated. Additionally, a review on the history, taxonomy, ecology, and medicinal application of different Nicotiana species growing around the globe presented in this Review may be of interest for pharmacognosists, natural products, and ecological chemists.

Comparative genetic analysis of the 45S rDNA intergenic spacers from three Saccharum species

The 45S ribosomal DNA (rDNA) units are separated by an intergenic spacer (IGS) containing the signals for transcription and processing of rRNAs. For the first time, we sequenced and analyzed the entire IGS region from three original species within the genus Saccharum, including S. spontaneum, S. robustum, and S. officinarum in this study. We have compared the IGS organization within three original species of the genus Saccharum. The IGS of these three original species showed similar overall organizations comprised of putative functional elements needed for rRNA gene activity as well as a non-transcribed spacer (NTS), a promoter region, and an external transcribed spacer (ETS). The variability in length of the IGS sequences was assessed at the individual, intraspecies, and interspecies levels of the genus Saccharum, including S. spontaneum, S. robustum, and S. officinarum. The ETS had greater similarity than the NTS across species, but nevertheless exhibited variation in length. Within the IGS of the Saccharum species, base substitutions and copy number variation of sub-repeat were causes of the divergence in IGS sequences. We also identified a significant number of methylation sites. Furthermore, fluorescent in situ hybridization (FISH) co-localization of IGS and pTa71 probes was detected on all representative species of the genus Saccharum tested. Taken together, the results of this study provide a better insight into the structure and organization of the IGS in the genus Saccharum.

Homoeologous chromosome pairing across the eukaryote phylogeny

During the past quarter century, molecular phylogenetic inferences have significantly resolved evolutionary relationships spanning the eukaryotic tree of life. With improved phylogenies in hand, the focus of systematics will continue to expand from estimating species relationships toward examining the evolution of specific, fundamental traits across the eukaryotic tree. Undoubtedly, this will expose knowledge gaps in the evolution of key traits, particularly with respect to non-model lineages. Here, we examine one such trait across eukaryotes-the regulation of homologous chromosome pairing during meiosis-as an illustrative example. Specifically, we present an overview of the breakdown of homologous chromosome pairing in model eukaryotes and provide a discussion of various meiotic aberrations that result in the failure of homolog recognition, with a particular focus on lineages with a history of hybridization and polyploidization, across major eukaryotic clades. We then explore what is known about these processes in natural and non-model eukaryotic taxa, thereby exposing disparities in our understanding of this key trait among non-model groups.

Higher-order organisation of extremely amplified, potentially functional and massively methylated 5S rDNA in European pikes (Esox sp.)

Background

Pikes represent an important genus (Esox) harbouring a pre-duplication karyotype (2n = 2x = 50) of economically important salmonid pseudopolyploids. Here, we have characterized the 5S ribosomal RNA genes (rDNA) in Esox lucius and its closely related E. cisalpinus using cytogenetic, molecular and genomic approaches. Intragenomic homogeneity and copy number estimation was carried out using Illumina reads. The higher-order structure of rDNA arrays was investigated by the analysis of long PacBio reads. Position of loci on chromosomes was determined by FISH. DNA methylation was analysed by methylation-sensitive restriction enzymes.

Results

The 5S rDNA loci occupy exclusively (peri)centromeric regions on 30–38 acrocentric chromosomes in both E. lucius and E. cisalpinus. The large number of loci is accompanied by extreme amplification of genes (>20,000 copies), which is to the best of our knowledge one of the highest copy number of rRNA genes in animals ever reported. Conserved secondary structures of predicted 5S rRNAs indicate that most of the amplified genes are potentially functional. Only few SNPs were found in genic regions indicating their high homogeneity while intergenic spacers were more heterogeneous and several families were identified. Analysis of 10–30 kb-long molecules sequenced by the PacBio technology (containing about 40% of total 5S rDNA) revealed that the vast majority (96%) of genes are organised in large several kilobase-long blocks. Dispersed genes or short tandems were less common (4%). The adjacent 5S blocks were directly linked, separated by intervening DNA and even inverted. The 5S units differing in the intergenic spacers formed both homogeneous and heterogeneous (mixed) blocks indicating variable degree of homogenisation between the loci. Both E. lucius and E. cisalpinus 5S rDNA was heavily methylated at CG dinucleotides.

Conclusions

Extreme amplification of 5S rRNA genes in the Esox genome occurred in the absence of significant pseudogenisation suggesting its recent origin and/or intensive homogenisation processes. The dense methylation of units indicates that powerful epigenetic mechanisms have evolved in this group of fish to silence amplified genes. We discuss how the higher-order repeat structures impact on homogenisation of 5S rDNA in the genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3774-7) contains supplementary material, which is available to authorized users.

Evolutional dynamics of 45S and 5S ribosomal DNA in ancient allohexaploid Atropa belladonna

BACKGROUND

Polyploid hybrids represent a rich natural resource to study molecular evolution of plant genes and genomes. Here, we applied a combination of karyological and molecular methods to investigate chromosomal structure, molecular organization and evolution of ribosomal DNA (rDNA) in nightshade, Atropa belladonna (fam. Solanaceae), one of the oldest known allohexaploids among flowering plants. Because of their abundance and specific molecular organization (evolutionarily conserved coding regions linked to variable intergenic spacers, IGS), 45S and 5S rDNA are widely used in plant taxonomic and evolutionary studies.

RESULTS

Molecular cloning and nucleotide sequencing of A. belladonna 45S rDNA repeats revealed a general structure characteristic of other Solanaceae species, and a very high sequence similarity of two length variants, with the only difference in number of short IGS subrepeats. These results combined with the detection of three pairs of 45S rDNA loci on separate chromosomes, presumably inherited from both tetraploid and diploid ancestor species, example intensive sequence homogenization that led to substitution/elimination of rDNA repeats of one parent. Chromosome silver-staining revealed that only four out of six 45S rDNA sites are frequently transcriptionally active, demonstrating nucleolar dominance. For 5S rDNA, three size variants of repeats were detected, with the major class represented by repeats containing all functional IGS elements required for transcription, the intermediate size repeats containing partially deleted IGS sequences, and the short 5S repeats containing severe defects both in the IGS and coding sequences. While shorter variants demonstrate increased rate of based substitution, probably in their transition into pseudogenes, the functional 5S rDNA variants are nearly identical at the sequence level, pointing to their origin from a single parental species. Localization of the 5S rDNA genes on two chromosome pairs further supports uniparental inheritance from the tetraploid progenitor.

CONCLUSIONS

The obtained molecular, cytogenetic and phylogenetic data demonstrate complex evolutionary dynamics of rDNA loci in allohexaploid species of Atropa belladonna. The high level of sequence unification revealed in 45S and 5S rDNA loci of this ancient hybrid species have been seemingly achieved by different molecular mechanisms.

Gene conversion events and variable degree of homogenization of rDNA loci in cultivars of Brassica napus

BACKGROUND AND AIMS

Brassica napus (AACC, 2n = 38, oilseed rape) is a relatively recent allotetraploid species derived from the putative progenitor diploid species Brassica rapa (AA, 2n = 20) and Brassica oleracea (CC, 2n = 18). To determine the influence of intensive breeding conditions on the evolution of its genome, we analysed structure and copy number of rDNA in 21 cultivars of B. napus, representative of genetic diversity.

METHODS

We used next-generation sequencing genomic approaches, Southern blot hybridization, expression analysis and fluorescence in situ hybridization (FISH). Subgenome-specific sequences derived from rDNA intergenic spacers (IGS) were used as probes for identification of loci composition on chromosomes.

KEY RESULTS

Most B. napus cultivars (18/21, 86 %) had more A-genome than C-genome rDNA copies. Three cultivars analysed by FISH ('Darmor', 'Yudal' and 'Asparagus kale') harboured the same number (12 per diploid set) of loci. In B. napus 'Darmor', the A-genome-specific rDNA probe hybridized to all 12 rDNA loci (eight on the A-genome and four on the C-genome) while the C-genome-specific probe showed weak signals on the C-genome loci only. Deep sequencing revealed high homogeneity of arrays suggesting that the C-genome genes were largely overwritten by the A-genome variants in B. napus 'Darmor'. In contrast, B. napus 'Yudal' showed a lack of gene conversion evidenced by additive inheritance of progenitor rDNA variants and highly localized hybridization signals of subgenome-specific probes on chromosomes. Brassica napus 'Asparagus kale' showed an intermediate pattern to 'Darmor' and 'Yudal'. At the expression level, most cultivars (95 %) exhibited stable A-genome nucleolar dominance while one cultivar ('Norin 9') showed co-dominance.

CONCLUSIONS

The B. napus cultivars differ in the degree and direction of rDNA homogenization. The prevalent direction of gene conversion (towards the A-genome) correlates with the direction of expression dominance indicating that gene activity may be needed for interlocus gene conversion.

Correlation analysis of the mRNA and miRNA expression profiles in the nascent synthetic allotetraploid Raphanobrassica

Raphanobrassica is an allopolyploid species derived from inter-generic hybridization that combines the R genome from R. sativus and the C genome from B. oleracea var. alboglabra. In the present study, we used a high-throughput sequencing method to identify the mRNA and miRNA profiles in Raphanobrassica and its parents. A total of 33,561 mRNAs and 283 miRNAs were detected, 9,209 mRNAs and 134 miRNAs were differentially expressed respectively, 7,633 mRNAs and 39 miRNAs showed ELD expression, 5,219 mRNAs and 57 miRNAs were non-additively expressed in Raphanobrassica. Remarkably, differentially expressed genes (DEGs) were up-regulated and maternal bias was detected in Raphanobrassica. In addition, a miRNA-mRNA interaction network was constructed based on reverse regulated miRNA-mRNAs, which included 75 miRNAs and 178 mRNAs, 31 miRNAs were non-additively expressed target by 13 miRNAs. The related target genes were significantly enriched in the GO term 'metabolic processes'. Non-additive related target genes regulation is involved in a range of biological pathways, like providing a driving force for variation and adaption in this allopolyploid. The integrative analysis of mRNA and miRNA profiling provides more information to elucidate gene expression mechanism and may supply a comprehensive and corresponding method to study genetic and transcription variation of allopolyploid.

Critical factors in the establishment of allopolyploids

PREMISE OF THE STUDY

The growth and spread of new polyploid populations have been explained in terms of fitness advantages over their diploid progenitors. However, a fitness advantage is not sufficient to insure the establishment of a polyploid; it must also overcome the obstacles of demographic stochasticity and minority disadvantage. Several studies have addressed the population dynamics of autopolyploids, but the present study is the first to consider allopolyploids, which are affected by more factors than autopolyploids.

METHODS

We constructed a population dynamic model of four types of plants (two parent species, hybrids, allopolyploids) that also included an explicit breeding system.

KEY RESULTS

The numbers of plants of each type were the most important factors determining whether the new allopolyploid would become established. More polyploid plants greatly increased the likelihood of polyploid persistence. More plants of the parent species and more hybrids resulted in more polyploids being produced. The model parameters with the most effect on polyploid establishment were potential population size (K), individual plant fecundity, and niche separation (α). The most important breeding system parameters were selfing rates, which mitigated minority disadvantage imposed by pollen limitation.

CONCLUSIONS

The importance of population sizes, and the parameters that controlled them, in overcoming demographic stochasticity parallels the well-recognized role of propagule pressure in determining the success of invasive species. We modeled the establishment of a new allopolyploid; analogous considerations would affect the establishment of a new autopolyploid. The critical role of population sizes in polyploid establishment should be more widely recognized.

Molecular and Cytogenetic Characterization of Wild Musa Species

The production of bananas is threatened by rapid spreading of various diseases and adverse environmental conditions. The preservation and characterization of banana diversity is essential for the purposes of crop improvement. The world's largest banana germplasm collection maintained at the Bioversity International Transit Centre (ITC) in Belgium is continuously expanded by new accessions of edible cultivars and wild species. Detailed morphological and molecular characterization of the accessions is necessary for efficient management of the collection and utilization of banana diversity. In this work, nuclear DNA content and genomic distribution of 45S and 5S rDNA were examined in 21 diploid accessions recently added to ITC collection, representing both sections of the genus Musa. 2C DNA content in the section Musa ranged from 1.217 to 1.315 pg. Species belonging to section Callimusa had 2C DNA contents ranging from 1.390 to 1.772 pg. While the number of 45S rDNA loci was conserved in the section Musa, it was highly variable in Callimusa species. 5S rRNA gene clusters were found on two to eight chromosomes per diploid cell. The accessions were genotyped using a set of 19 microsatellite markers to establish their relationships with the remaining accessions held at ITC. Genetic diversity done by SSR genotyping platform was extended by phylogenetic analysis of ITS region. ITS sequence data supported the clustering obtained by SSR analysis for most of the accessions. High level of nucleotide diversity and presence of more than two types of ITS sequences in eight wild diploids pointed to their origin by hybridization of different genotypes. This study significantly expands the number of wild Musa species where nuclear genome size and genomic distribution of rDNA loci is known. SSR genotyping identified Musa species that are closely related to the previously characterized accessions and provided data to aid in their classification. Sequence analysis of ITS region provided further information about evolutionary relationships between individual accessions and suggested that some of analyzed accessions were interspecific hybrids and/or backcross progeny.

Nonadditive gene expression in polyploids

Allopolyploidy involves hybridization and duplication of divergent parental genomes and provides new avenues for gene expression. The expression levels of duplicated genes in polyploids can show deviation from parental additivity (the arithmetic average of the parental expression levels). Nonadditive expression has been widely observed in diverse polyploids and comprises at least three possible scenarios: (a) The total gene expression level in a polyploid is similar to that of one of its parents (expression-level dominance); (b) total gene expression is lower or higher than in both parents (transgressive expression); and (c) the relative contribution of the parental copies (homeologs) to the total gene expression is unequal (homeolog expression bias). Several factors may result in expression nonadditivity in polyploids, including maternal-paternal influence, gene dosage balance, cis- and/or trans-regulatory networks, and epigenetic regulation. As our understanding of nonadditive gene expression in polyploids remains limited, a new generation of investigators should explore additional phenomena (i.e., alternative splicing) and use other high-throughput "omics" technologies to measure the impact of nonadditive expression on phenotype, proteome, and metabolome.

Evolutionary origins and dynamics of octoploid strawberry subgenomes revealed by dense targeted capture linkage maps

Whole-genome duplications are radical evolutionary events that have driven speciation and adaptation in many taxa. Higher-order polyploids have complex histories often including interspecific hybridization and dynamic genomic changes. This chromosomal reshuffling is poorly understood for most polyploid species, despite their evolutionary and agricultural importance, due to the challenge of distinguishing homologous sequences from each other. Here, we use dense linkage maps generated with targeted sequence capture to improve the diploid strawberry (Fragaria vesca) reference genome and to disentangle the subgenomes of the wild octoploid progenitors of cultivated strawberry, Fragaria virginiana and Fragaria chiloensis. Our novel approach, POLiMAPS (Phylogenetics Of Linkage-Map-Anchored Polyploid Subgenomes), leverages sequence reads to associate informative interhomeolog phylogenetic markers with linkage groups and reference genome positions. In contrast to a widely accepted model, we find that one of the four subgenomes originates with the diploid cytoplasm donor F. vesca, one with the diploid Fragaria iinumae, and two with an unknown ancestor close to F. iinumae. Extensive unidirectional introgression has converted F. iinumae-like subgenomes to be more F. vesca-like, but never the reverse, due either to homoploid hybridization in the F. iinumae-like diploid ancestors or else strong selection spreading F. vesca-like sequence among subgenomes through homeologous exchange. In addition, divergence between homeologous chromosomes has been substantially augmented by interchromosomal rearrangements. Our phylogenetic approach reveals novel aspects of the complicated web of genetic exchanges that occur during polyploid evolution and suggests a path forward for unraveling other agriculturally and ecologically important polyploid genomes.

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.

Doubling down on genomes: polyploidy and crop plants

Polyploidy, or whole genome multiplication, is ubiquitous among angiosperms. Many crop species are relatively recent allopolyploids, resulting from interspecific hybridization and polyploidy. Thus, an appreciation of the evolutionary consequences of (allo)polyploidy is central to our understanding of crop plant domestication, agricultural improvement, and the evolution of angiosperms in general. Indeed, many recent insights into plant biology have been gleaned from polyploid crops, including, but not limited to wheat, tobacco, sugarcane, apple, and cotton. A multitude of evolutionary processes affect polyploid genomes, including rapid and substantial genome reorganization, transgressive gene expression alterations, gene fractionation, gene conversion, genome downsizing, and sub- and neofunctionalization of duplicate genes. Often these genomic changes are accompanied by heterosis, robustness, and the improvement of crop yield, relative to closely related diploids. Historically, however, the genome-wide analysis of polyploid crops has lagged behind those of diploid crops and other model organisms. This lag is partly due to the difficulties in genome assembly, resulting from the genomic complexities induced by combining two or more evolutionarily diverged genomes into a single nucleus and by the significant size of polyploid genomes. In this review, we explore the role of polyploidy in angiosperm evolution, the domestication process and crop improvement. We focus on the potential of modern technologies, particularly next-generation sequencing, to inform us on the patterns and processes governing polyploid crop improvement and phenotypic change subsequent to domestication.

When fathers are instant losers: homogenization of rDNA loci in recently formed Cardamine × schulzii trigenomic allopolyploid

Recently formed allopolyploids represent an excellent system to study the impacts of hybridization and genomic duplication on genome structure and evolution. Here we explored the 35SrRNA genes (rDNA) in the Cardamine × schulzii allohexaploid that was formed by two subsequent hybridization events within the past c. 150 yr. The rDNA loci were analyzed by cloning, next generation sequencing (NGS), RT-PCR and FISH methods. The primary C. × insueta triploid hybrid derived from C. rivularis (♀) and C. amara (♂) had gene ratios highly skewed towards maternal sequences. Similarly, C. × schulzii, originating from the secondary hybridization event involving C. × insueta (♀) and C. pratensis (♂), showed a reduction in paternal rDNA homeologs despite an excess of chromosomes inherited from C. pratensis. We also identified novel rDNA loci in C. × schulzii, suggesting that lost loci might be slowly reinstalled by translocation (but not recombination) of genes from partner genomes. Prevalent clonal propagation of allopolyploids, C. × insueta and C. × schulzii, indicates that concerted evolution of rDNA may occur in the absence of extensive meiotic cycles. Adoption of NGS in rDNA variant analysis is highly informative for deciphering the evolutionary histories of allopolyploid species with ongoing homogenization processes.

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.

The widespread crucifer species Cardamine flexuosa is an allotetraploid with a conserved subgenomic structure

The origin of Cardamine flexuosa (Wavy Bittercress) has been a conundrum for more than six decades. Here we identify its parental species, analyse its genome structure in comparison to parental genomes and describe intergenomic structural variations in C. flexuosa. Genomic in situ hybridization (GISH) and comparative chromosome painting (CCP) uncovered the parental genomes and the chromosome composition of C. flexuosa and its presumed diploid progenitors. Cardamine flexuosa is an allotetraploid (2n = 4x = 32), originating from two diploid species, Cardamine amara and Cardamine hirsuta (2n = 2x = 16). The two parental species display almost perfectly conserved chromosomal collinearity for seven out of the eight chromosomes. A 13 Mb pericentric inversion distinguishes chromosome CA1 from CH1. A comparative cytomolecular map was established for C. flexuosa by CCP/GISH. Whereas conserved chromosome collinearity between the C. amara and C. hirsuta subgenomes might have promoted intergenomic rearrangements through homeologous recombination, only one reciprocal translocation between two homeologues has occurred since the origin of C. flexuosa. The genome of C. flexuosa demonstrates that allopolyploids can maintain remarkably stable subgenomes over 10(4) -10(5)  yr throughout a wide distribution range. By contrast, the rRNA genes underwent genome-specific elimination towards a diploid-like number of loci.

Host adaptation and speciation through hybridization and polyploidy in Phytophthora

It is becoming increasingly evident that interspecific hybridization is a common event in phytophthora evolution. Yet, the fundamental processes underlying interspecific hybridization and the consequences for its ecological fitness and distribution are not well understood. We studied hybridization events in phytophthora clade 8b. This is a cold-tolerant group of plant pathogenic oomycetes in which six host-specific species have been described that mostly attack winter-grown vegetables. Hybrid characterization was done by sequencing and cloning of two nuclear (ITS and Ypt1) and two mitochondrial loci (Cox1 and Nadh1) combined with DNA content estimation using flow cytometry. Three different mtDNA haplotypes were recovered among the presumed hybrid isolates, dividing the hybrids into three types, with different parental species involved. In the nuclear genes, additivity, i.e. the presence of two alleles coming from different parents, was detected. Hybrid isolates showed large variations in DNA content, which was positively correlated with the additivity in nuclear loci, indicating allopolyploid hybridization followed by a process of diploidization. Moreover, indications of homeologous recombination were found in the hybrids by cloning ITS products. The hybrid isolates have been isolated from a range of hosts that have not been reported previously for clade 8b species, indicating that they have novel pathogenic potential. Next to this, DNA content measurements of the non-hybrid clade 8b species suggest that polyploidy is a common feature of this clade. We hypothesize that interspecific hybridization and polyploidy are two linked phenomena in phytophthora, and that these processes might play an important and ongoing role in the evolution of this genus.

Genomic characterization of interspecific hybrids between the scallops Argopecten purpuratus and A. irradians irradians

The Peruvian scallop (Argopecten purpuratus) has been introduced to China and has successfully been hybridized with the bay scallop (A. irradians irradians). The F1 hybrids of these two scallops exhibited a large increase in production traits and some other interesting new characteristics. To understand the genetic basis of this heterosis, nuclear gene and partial mtDNA sequences, and genomic in situ hybridization (GISH) were employed to analyze the genomic organization of the hybrids. Amplification of the ribosomal DNA internal transcribed spacer (ITS) showed that the parental ITS sequences were present in all the hybrid individuals, illustrating that the hybrid offspring inherited nuclear DNA from both parents. Sequence analyses of the ITS region further confirmed that the hybrids harbored alleles from their parents; some recombinant variants were also detected, which revealed some alterations in the nuclear genetic material of the hybrids. The analysis of mitochondrial 16S rDNA showed that the hybrids possessed sequences that were identical to the 16S rDNA of the female parents, proving a matrilineal inheritance of mitochondrial genes in scallops. In addition, GISH clearly discriminated between the parental chromosomes and indicated a combination of haploid genomes of duplex parents in the hybrids. The genetic analyses in our study illustrated that the F1 hybrids inherited nuclear material from both parents and cytoplasmic genetic material maternally, and some variations occurred in the genome, which might contribute to a further understanding of crossbreeding and heterosis in scallop species.

Next generation sequencing analysis reveals a relationship between rDNA unit diversity and locus number in Nicotiana diploids

BACKGROUND

Tandemly arranged nuclear ribosomal DNA (rDNA), encoding 18S, 5.8S and 26S ribosomal RNA (rRNA), exhibit concerted evolution, a pattern thought to result from the homogenisation of rDNA arrays. However rDNA homogeneity at the single nucleotide polymorphism (SNP) level has not been detailed in organisms with more than a few hundred copies of the rDNA unit. Here we study rDNA complexity in species with arrays consisting of thousands of units.

METHODS

We examined homogeneity of genic (18S) and non-coding internally transcribed spacer (ITS1) regions of rDNA using Roche 454 and/or Illumina platforms in four angiosperm species, Nicotiana sylvestris, N. tomentosiformis, N. otophora and N. kawakamii. We compared the data with Southern blot hybridisation revealing the structure of intergenic spacer (IGS) sequences and with the number and distribution of rDNA loci.

RESULTS AND CONCLUSIONS

In all four species the intragenomic homogeneity of the 18S gene was high; a single ribotype makes up over 90% of the genes. However greater variation was observed in the ITS1 region, particularly in species with two or more rDNA loci, where >55% of rDNA units were a single ribotype, with the second most abundant variant accounted for >18% of units. IGS heterogeneity was high in all species. The increased number of ribotypes in ITS1 compared with 18S sequences may reflect rounds of incomplete homogenisation with strong selection for functional genic regions and relaxed selection on ITS1 variants. The relationship between the number of ITS1 ribotypes and the number of rDNA loci leads us to propose that rDNA evolution and complexity is influenced by locus number and/or amplification of orphaned rDNA units at new chromosomal locations.

Analysis of two abundant, highly related satellites in the allotetraploid Nicotiana arentsii using double-strand conformation polymorphism analysis and sequencing

• Allopolyploidy, a driving force in plant evolution, can induce rapid structural changes in parental subgenomes. Here, we examined the fate of homologous subtelomeric satellites in intrasection allotetraploid Nicotiana arentsii formed from N. undulata and N. wigandioides progenitors < 200,000 yr ago. • We cloned and sequenced a number of monomers from progenitors and the allotetraploid. Structural features of both cloned and genomic monomers were studied using double-strand conformation polymorphism analysis. • Two homologous satellites were isolated from N. undulata (called NUNSSP) and N. wigandioides (NWISSP). While the NUNSSP monomers were highly homogeneous in nucleotide sequences, the NWISSP monomers formed two separate clades. Likewise, the genomic NUNSSP monomers showed less DNA conformation heterogeneity than NWISSP monomers, with distinct conformations. While both satellites predominantly occupy subtelomeric positions, a fraction of the NWISSP repeats was found in an intercalary location, supporting the hypothesis that dispersion prevents the repeats becoming homogeneous. Sequence, structural and chromosomal features of the parental satellites were faithfully inherited by N. arentsii. • Our study revealed that intergenomic homogenization of subtelomeric satellite repeats does not occur in N. arentsii allotetraploid. We propose that the sequence and structural divergence of subtelomeric satellites may render allopolyploid chromosomes less vulnerable to intergenomic exchanges.

Karyotype analysis and genetic variation of a mutant in Siraitia grosvenorii

This study analyzed the differences in karyotype and genetic variation between a mutant and wild-type Siraitia grosvenorii. Genetic variation included changes in genome and gene expression by SRAP molecular markers. Results showed that wild-type S. grosvenorii was diploid, with a chromosome number of 2n = 2x = 28, whereas the mutant was tetraploid with a chromosome number of 2n = 4x = 56. 4573 DNA bands were obtained using 189 different primer combinations, 577 of which were polymorphic, averaging 3.1 bands for each primer pair, while 1998 pairs were identical. There were no apparent differences on bands amplified by most primer pairs. After comparing the diploid and tetraploid strains, the data generally indicated that the polymorphism would be quite low. 2917 cDNA bands were generated using 133 primer combinations, and stable and clearly differential fragments were sorted out, cloned and sequenced. Ninety-two differentially expressed fragments were successfully sequenced. Sequence analysis showed that most fragments had significant homologous nucleotide sequences with resistant to stress and photosynthesis genes, including ribulose-1,5-bisphosphate carboxylase/oxygenase, phosphoenolpyruvate carboxykinase, pyruvate kinase, peroxisomal membrane transporter, NBS-LRR type resistance protein, protein phosphatase and others. The results revealed that the tetraploid strain has more resistant and photosynthesis ability than its diploid relatives, which providing reference information and resources for molecular breeding and seedless Luohanguo.

The ITS1-5.8S-ITS2 sequence region in the Musaceae: structure, diversity and use in molecular phylogeny

Genes coding for 45S ribosomal RNA are organized in tandem arrays of up to several thousand copies and contain 18S, 5.8S and 26S rRNA units separated by internal transcribed spacers ITS1 and ITS2. While the rRNA units are evolutionary conserved, ITS show high level of interspecific divergence and have been used frequently in genetic diversity and phylogenetic studies. In this work we report on the structure and diversity of the ITS region in 87 representatives of the family Musaceae. We provide the first detailed information on ITS sequence diversity in the genus Musa and describe the presence of more than one type of ITS sequence within individual species. Both Sanger sequencing of amplified ITS regions and whole genome 454 sequencing lead to similar phylogenetic inferences. We show that it is necessary to identify putative pseudogenic ITS sequences, which may have negative effect on phylogenetic reconstruction at lower taxonomic levels. Phylogenetic reconstruction based on ITS sequence showed that the genus Musa is divided into two distinct clades – Callimusa and Australimusa and Eumusa and Rhodochlamys. Most of the intraspecific banana hybrids analyzed contain conserved parental ITS sequences, indicating incomplete concerted evolution of rDNA loci. Independent evolution of parental rDNA in hybrids enables determination of genomic constitution of hybrids using ITS. The observation of only one type of ITS sequence in some of the presumed interspecific hybrid clones warrants further study to confirm their hybrid origin and to unravel processes leading to evolution of their genomes.

Contrasting patterns of the 5S and 45S rDNA evolutions in the Byblis liniflora complex (Byblidaceae)

To clarify the evolutionary dynamics of ribosomal RNA genes (rDNAs) in the Byblis liniflora complex (Byblidaceae), we investigated the 5S and 45S rDNA genes through (1) chromosomal physical mapping by fluorescence in situ hybridization (FISH) and (2) phylogenetic analyses using the nontranscribed spacer of 5S rDNA (5S-NTS) and the internal transcribed spacer of 45S rDNA (ITS). In addition, we performed phylogenetic analyses based on rbcL and trnK intron. The complex was divided into 2 clades: B. aquatica-B. filifolia and B. guehoi-B. liniflora-B. rorida. Although members of the complex had conservative symmetric karyotypes, they were clearly differentiated on chromosomal rDNA distribution patterns. The sequence data indicated that ITS was almost homogeneous in all taxa in which two or four 45S rDNA arrays were frequently found at distal regions of chromosomes in the somatic karyotype. ITS homogenization could have been prompted by relatively distal 45S rDNA positions. In contrast, 2-12 5S rDNA arrays were mapped onto proximal/interstitial regions of chromosomes, and some paralogous 5S-NTS were found in the genomes harboring 4 or more arrays. 5S-NTS sequence type-specific FISH analysis showed sequence heterogeneity within and between some 5S rDNA arrays. Interlocus homogenization may have been hampered by their proximal location on chromosomes. Chromosomal location may have affected the contrasting evolutionary dynamics of rDNAs in the B. liniflora complex.