Transcriptional analysis of nucleolar dominance in polyploid plants: biased expression/silencing of progenitor rRNA genes is developmentally regulated in Brassica

Transgenerational epigenetic inheritance during plant evolution and breeding

Plants can program and reprogram their genomes to create genetic variation and epigenetic modifications, leading to phenotypic plasticity. Although consequences of genetic changes are comprehensible, the basis for transgenerational inheritance of epigenetic variation is elusive. This review addresses contributions of external (environmental) and internal (genomic) factors to the establishment and maintenance of epigenetic memory during plant evolution, crop domestication, and modern breeding. Dynamic and pervasive changes in DNA methylation and chromatin modifications provide a diverse repertoire of epigenetic variation potentially for transgenerational inheritance. Elucidating and harnessing epigenetic inheritance will help us develop innovative breeding strategies and biotechnological tools to improve crop yield and resilience in the face of environmental challenges. Beyond plants, epigenetic principles are shared across sexually reproducing organisms including humans with relevance to medicine and public health.

Transcriptional Silencing of 35S rDNA in Tragopogon porrifolius Correlates with Cytosine Methylation in Sequence-Specific Manner

Despite the widely accepted involvement of DNA methylation in the regulation of rDNA transcription, the relative participation of different cytosine methylation pathways is currently described only for a few model plants. Using PacBio, Bisulfite, and RNA sequencing; PCR; Southern hybridizations; and FISH, the epigenetic consequences of rDNA copy number variation were estimated in two T. porrifolius lineages, por1 and por2, the latter with more than twice the rDNA copy numbers distributed approximately equally between NORs on chromosomes A and D. The lower rDNA content in por1 correlated with significantly reduced (>90%) sizes of both D-NORs. Moreover, two (L and S) prominent rDNA variants, differing in the repetitive organization of intergenic spacers, were detected in por2, while only the S-rDNA variant was detected in por1. Transcriptional activity of S-rDNA in por1 was associated with secondary constriction of both A-NORs. In contrast, silencing of S-rDNA in por2 was accompanied by condensation of A-NORs, secondary constriction on D-NORs, and L-rDNA transcriptional activity, suggesting (i) bidirectional nucleolar dominance and (ii) association of S-rDNAs with A-NORs and L-rDNAs with D-NORs in T. porrifolius. Each S- and L-rDNA array was formed of several sub-variants differentiating both genetically (specific SNPs) and epigenetically (transcriptional efficiency and cytosine methylation). The most significant correlations between rDNA silencing and methylation were detected for symmetric CWG motifs followed by CG motifs. No correlations were detected for external cytosine in CCGs or asymmetric CHHs, where methylation was rather position-dependent, particularly for AT-rich variants. We conclude that variations in rDNA copy numbers in plant diploids can be accompanied by prompt epigenetic responses to maintain an appropriate number of active rDNAs. The methylation dynamics of CWGs are likely to be the most responsible for regulating silent and active rDNA states.

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

Arabidopsis ASYMMETRIC LEAVES2 and Nucleolar Factors Are Coordinately Involved in the Perinucleolar Patterning of AS2 Bodies and Leaf Development

Arabidopsis ASYMMETRIC LEAVES2 (AS2) plays a key role in the formation of flat symmetric leaves. AS2 represses the expression of the abaxial gene ETTIN/AUXIN RESPONSE FACTOR3 (ETT/ARF3). AS2 interacts in vitro with the CGCCGC sequence in ETT/ARF3 exon 1. In cells of leaf primordia, AS2 localizes at peripheral regions of the nucleolus as two AS2 bodies, which are partially overlapped with chromocenters that contain condensed 45S ribosomal DNA repeats. AS2 contains the AS2/LOB domain, which consists of three sequences conserved in the AS2/LOB family: the zinc finger (ZF) motif, the ICG sequence including the conserved glycine residue, and the LZL motif. AS2 and the genes NUCLEOLIN1 (NUC1), RNA HELICASE10 (RH10), and ROOT INITIATION DEFECTIVE2 (RID2) that encode nucleolar proteins coordinately act as repressors against the expression of ETT/ARF3. Here, we examined the formation and patterning of AS2 bodies made from as2 mutants with amino acid substitutions in the ZF motif and the ICG sequence in cells of cotyledons and leaf primordia. Our results showed that the amino acid residues next to the cysteine residues in the ZF motif were essential for both the formation of AS2 bodies and the interaction with ETT/ARF3 DNA. The conserved glycine residue in the ICG sequence was required for the formation of AS2 bodies, but not for the DNA interaction. We also examined the effects of nuc1, rh10, and rid2 mutations, which alter the metabolism of rRNA intermediates and the morphology of the nucleolus, and showed that more than two AS2 bodies were observed in the nucleolus and at its periphery. These results suggested that the patterning of AS2 bodies is tightly linked to the morphology and functions of the nucleolus and the development of flat symmetric leaves in plants.

New insights on the evolution of nucleolar dominance in newly resynthesized hexaploid wheat Triticum zhukovskyi

Nucleolar dominance (ND) is a widespread epigenetic phenomenon in hybridizations where nucleolus transcription fails at the nucleolus organizer region (NOR). However, the dynamics of NORs during the formation of Triticum zhukovskyi (GGAu Au Am Am ), another evolutionary branch of allohexaploid wheat, remains poorly understood. Here, we elucidated genetic and epigenetic changes occurring at the NOR loci within the Am , G, and D subgenomes during allopolyploidization by synthesizing hexaploid wheat GGAu Au Am Am and GGAu Au DD. In T. zhukovskyi, Au genome NORs from T. timopheevii (GGAu Au ) were lost, while the second incoming NORs from T. monococcum (Am Am ) were retained. Analysis of the synthesized T. zhukovskyi revealed that rRNA genes from the Am genome were silenced in F1 hybrids (GAu Am ) and remained inactive after genome doubling and subsequent self-pollinations. We observed increased DNA methylation accompanying the inactivation of NORs in the Am genome and found that silencing of NORs in the S1 generation could be reversed by a cytidine methylase inhibitor. Our findings provide insights into the ND process during the evolutionary period of T. zhukovskyi and highlight that inactive rDNA units may serve as a 'first reserve' in the form of R-loops, contributing to the successful evolution of T. zhukovskyi.

Reaching for the off switch in nucleolar dominance

Nucleolus organizer regions (NORs) are eukaryotic chromosomal loci where ribosomal RNA (rRNA) genes are clustered, typically in hundreds to thousands of copies. Transcription of these rRNA genes by RNA polymerase I and processing of their transcripts results in the formation of the nucleolus, the sub-nuclear domain in which ribosomes are assembled. Approximately 90 years ago, cytogenetic observations revealed that NORs inherited from the different parents of an interspecific hybrid sometimes differ in morphology at metaphase. Fifty years ago, those chromosomal differences were found to correlate with differences in rRNA gene transcription and the phenomenon became known as nucleolar dominance. Studies of the past 30 years have revealed that nucleolar dominance results from selective rRNA gene silencing, involving repressive chromatin modifications, and occurs in pure species as well as hybrids. Recent evidence also indicates that silencing depends on the NOR in which an rRNA gene is located, and not on the gene's sequence. In this perspective, we discuss how our thinking about nucleolar dominance has shifted over time from the kilobase scale of individual genes to the megabase scale of NORs and chromosomes and questions that remain unanswered in the search for a genetic and biochemical understanding of the off switch.

Switch them off or not: selective rRNA gene repression in grasses

Nucleolar dominance (ND) is selective epigenetic silencing of 35-48S rDNA loci. In allopolyploids, it is frequently manifested at the cytogenetic level by the inactivation of nucleolar organiser region(s) (NORs) inherited from one or several evolutionary ancestors. Grasses are ecologically and economically one of the most important land plant groups, which have frequently evolved through hybridisation and polyploidisation events. Here we review common and unique features of ND phenomena in this monocot family from cytogenetic, molecular, and genomic perspectives. We highlight recent advances achieved by using an allotetraploid model grass, Brachypodium hybridum, where ND commonly occurs at a population level, and we cover modern genomic approaches that decipher structural features of core arrays of NORs.

Structure and Methylation of 35S rDNA in Allopolyploids Anemone multifida (2n = 4x = 32, BBDD) and Anemone baldensis (2n = 6x = 48, AABBDD) and Their Parental Species Show Evidence of Nucleolar Dominance

Transcriptional silencing of 35S rDNA loci inherited from one parental species is occurring relatively frequently in allopolyploids. However, molecular mechanisms by which it is selected for transcriptional silencing remain unclear. We applied NGS, silver staining and bisulfite sequencing to study the structure, expression and methylation landscape of 35S rDNA in two allopolyploids of common origin, allotetraploid Anemone multifida (2n = 4x = 32, genome composition BBDD) and allohexaploid A. baldensis (2n = 6x = 48, AABBDD), and their genome donors, A. sylvestris (2n = 16, AA), A. cylindrica (2n = 16, BB) and A. parviflora (2n = 16, DD). The size of the recovered 35S rDNA units varied from 10,489 bp in A. cylindrica to 12,084 bp in A. sylvestris. Anemone showed an organization typical of most ribosomal 35S rDNA composed of NTS, ETS, rRNA genes, TTS and TIS with structural features of plant IGS sequences and all functional elements needed for rRNA gene activity. The NTS was more variable than the ETS and consisted of SRs which are highly variable among Anemone. Five to six CpG-rich islands were found within the ETS. CpG island located adjacent to the transcription initiation site (TIS) was highly variable regarding the sequence size and methylation level and exhibited in most of the species lower levels of methylation than CpG islands located adjacent to the 18S rRNA gene. Our results uncover hypomethylation of A. sylvestris- and A. parviflora-derived 35S rDNA units in allopolyploids A. multifida and A. baldensis. Hypomethylation of A. parviflora-derived 35S rDNA was more prominent in A. baldensis than in A. multifida. We showed that A. baldensis underwent coupled A. sylvestris-derived 35S rDNA array expansion and A. parviflora-derived 35S rDNA copy number decrease that was accompanied by lower methylation level of A. sylvestris-derived 35S rDNA units in comparison to A. parviflora-derived 35S rDNA units. These observations suggest that in A. baldensis nucleolar dominance is directed toward A. sylvestris-derived chromosomes. This work broadens our current knowledge of the 35S rDNA organization in Anemone and provides evidence of the progenitor-specific 35S rDNA methylation in nucleolar dominance.

Rapid and Synchronous Breeding of Cytoplasmic Male Sterile and Maintainer Line Through Mitochondrial DNA Rearrangement Using Doubled Haploid Inducer in Brassica napus

When homozygously fertile plants were induced using doubled haploid (DH) induction lines Y3380 and Y3560, the morphology of the induced F1 generation was basically consistent with the female parent, but the fertility was separated, showing characteristics similar to cytoplasmic male sterile (CMS) and maintainer lines. In this study, the morphology, fertility, ploidy, and cytoplasm genotype of the induced progeny were identified, and the results showed that the sterile progeny was polima cytoplasm sterile (pol CMS) and the fertile progeny was nap cytoplasm. The molecular marker and test-cross experimental results showed that the fertile progeny did not carry the restorer gene of pol CMS and the genetic distance between the female parent and the offspring was 0.002. This suggested that those inductions which produced sterile and fertile progeny were coordinated to CMS and maintainer lines. Through the co-linearity analysis of the mitochondrial DNA (mtDNA), it was found that the rearrangement of mtDNA by DH induction was the key factor that caused the transformation of fertility (nap) into sterility (pol). Also, when heterozygous females were induced with DH induction lines, the induction F2 generation also showed the segregation of fertile and sterile lines, and the genetic distance between sterile and fertile lines was approximately 0.075. Therefore, the induction line can induce different types of female parents, and the breeding of the sterile line and the maintainer line can be achieved through the rapid synchronization of sister crosses and self-crosses. The induction of DH inducer in B. napus can provide a new model for the innovation of germplasm resources and open up a new way for its application.

Transcriptome shock in interspecific F1 allotriploid hybrids between Brassica species

Interspecific hybridization drives the evolution of angiosperms and can be used to introduce novel alleles for important traits or to activate heterosis in crop breeding. Hybridization brings together gene expression networks from two different species, potentially causing global alterations of gene expression in the F1 plants which is called 'transcriptome shock'. Here, we explored such a transcriptome shock in allotriploid Brassica hybrids. We generated interspecific F1 allotriploid hybrids between the allotetraploid species Brassica napus and three accessions of the diploid species Brassica rapa. RNA-seq of the F1 hybrids and the parental plants revealed that 26.34-30.89% of genes were differentially expressed between the parents. We also analyzed expression level dominance and homoeolog expression bias between the parents and the F1 hybrids. The expression-level dominance biases of the Ar, An, and Cn subgenomes was genotype and stage dependent, whereas significant homoeolog expression bias was observed among three subgenomes from different parents. Furthermore, more genes were involved in trans regulation than in cis regulation in allotriploid F1 hybrids. Our findings provide new insights into the transcriptomic responses of cross-species hybrids and hybrids showing heterosis, as well as a new method for promoting the breeding of desirable traits in polyploid Brassica species.

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.

Nonadditive gene expression is correlated with nonadditive phenotypic expression in interspecific triploid hybrids of willow (Salix spp.)

Many studies have highlighted the complex and diverse basis for heterosis in inbred crops. Despite the lack of a consensus model, it is vital that we turn our attention to understanding heterosis in undomesticated, heterozygous, and polyploid species, such as willow (Salix spp.). Shrub willow is a dedicated energy crop bred to be fast-growing and high yielding on marginal land without competing with food crops. A trend in willow breeding is the consistent pattern of heterosis in triploids produced from crosses between diploid and tetraploid species. Here, we test whether differentially expressed genes are associated with heterosis in triploid families derived from diploid Salix purpurea, diploid Salix viminalis, and tetraploid Salix miyabeana parents. Three biological replicates of shoot tips from all family progeny and parents were collected after 12 weeks in the greenhouse and RNA extracted for RNA-Seq analysis. This study provides evidence that nonadditive patterns of gene expression are correlated with nonadditive phenotypic expression in interspecific triploid hybrids of willow. Expression-level dominance was most correlated with heterosis for biomass yield traits and was highly enriched for processes involved in starch and sucrose metabolism. In addition, there was a global dosage effect of parent alleles in triploid hybrids, with expression proportional to copy number variation. Importantly, differentially expressed genes between family parents were most predictive of heterosis for both field and greenhouse collected traits. Altogether, these data will be used to progress models of heterosis to complement the growing genomic resources available for the improvement of heterozygous perennial bioenergy crops.

Fine structure and transcription dynamics of bread wheat ribosomal DNA loci deciphered by a multi-omics approach

Three out of four RNA components of ribosomes are encoded by 45S ribosomal DNA (rDNA) loci, which are organized as long head-to-tail tandem arrays of nearly identical units, spanning several megabases of sequence. Due to this structure, the rDNA loci are the major sources of gaps in genome assemblies, and gene copy number, sequence composition, and expression status of particular arrays remain elusive, especially in complex genomes harboring multiple loci. Here we conducted a multi-omics study to decipher the 45S rDNA loci in hexaploid bread wheat. Coupling chromosomal genomics with optical mapping, we reconstructed individual rDNA arrays, enabling locus-specific analyses of transcription activity and methylation status from RNA- and bisulfite-sequencing data. We estimated a total of 6,650 rDNA units in the bread wheat genome, with approximately 2,321, 3,910, 253, and 50 gene copies located in short arms of chromosomes 1B, 6B, 5D, and 1A, respectively. Only 1B and 6B loci contributed substantially to rRNA transcription at a roughly 2:1 ratio. The ratio varied among five tissues analyzed (embryo, coleoptile, root tip, primary leaf, mature leaf), being the highest (2.64:1) in mature leaf and lowest (1.72:1) in coleoptile. Cytosine methylation was considerably higher in CHG context in the silenced 5D locus as compared with the active 1B and 6B loci. In conclusion, a fine genomic organization and tissue-specific expression of rDNA loci were deciphered, for the first time, in a complex polyploid species. The results are discussed in the context of wheat evolution and transcription regulation.

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.

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

Feasible strategies for studying the involvement of DNA methylation and histone acetylation in the stress-induced formation of quality-related metabolites in tea (Camellia sinensis)

Tea plants are subjected to multiple stresses during growth, development, and postharvest processing, which affects levels of secondary metabolites in leaves and influences tea functional properties and quality. Most studies on secondary metabolism in tea have focused on gene, protein, and metabolite levels, whereas upstream regulatory mechanisms remain unclear. In this review, we exemplify DNA methylation and histone acetylation, summarize the important regulatory effects that epigenetic modifications have on plant secondary metabolism, and discuss feasible research strategies to elucidate the underlying specific epigenetic mechanisms of secondary metabolism regulation in tea. This information will help researchers investigate the epigenetic regulation of secondary metabolism in tea, providing key epigenetic data that can be used for future tea genetic breeding.

Trajectories of Homoeolog-Specific Expression in Allotetraploid Tragopogon castellanus Populations of Independent Origins

Polyploidization can have a significant ecological and evolutionary impact by providing substantially more genetic material that may result in novel phenotypes upon which selection may act. While the effects of polyploidization are broadly reviewed across the plant tree of life, the reproducibility of these effects within naturally occurring, independently formed polyploids is poorly characterized. The flowering plant genus Tragopogon (Asteraceae) offers a rare glimpse into the intricacies of repeated allopolyploid formation with both nascent (< 90 years old) and more ancient (mesopolyploids) formations. Neo- and mesopolyploids in Tragopogon have formed repeatedly and have extant diploid progenitors that facilitate the comparison of genome evolution after polyploidization across a broad span of evolutionary time. Here, we examine four independently formed lineages of the mesopolyploid Tragopogon castellanus for homoeolog expression changes and fractionation after polyploidization. We show that expression changes are remarkably similar among these independently formed polyploid populations with large convergence among expressed loci, moderate convergence among loci lost, and stochastic silencing. We further compare and contrast these results for T. castellanus with two nascent Tragopogon allopolyploids. While homoeolog expression bias was balanced in both nascent polyploids and T. castellanus, the degree of additive expression was significantly different, with the mesopolyploid populations demonstrating more non-additive expression. We suggest that gene dosage and expression noise minimization may play a prominent role in regulating gene expression patterns immediately after allopolyploidization as well as deeper into time, and these patterns are conserved across independent polyploid lineages.

Histone H3K27 dimethylation landscapes contribute to genome stability and genetic recombination during wheat polyploidization

Bread wheat (Triticum aestivum) is an allohexaploid that was formed via two allopolyploidization events. Growing evidence suggests histone modifications are involved in the response to 'genomic shock' and environmental adaptation during polyploid formation and evolution. However, the role of histone modifications, especially histone H3 lysine-27 dimethylation (H3K27me2), in genome evolution remains elusive. Here we analyzed H3K27me2 and H3K27me3 profiles in hexaploid wheat and its tetraploid and diploid relatives. Although H3K27me3 levels were relatively stable among wheat species with different ploidy levels, H3K27me2 intensities increased concurrent with increased ploidy levels, and H3K27me2 peaks were colocalized with massively amplified DTC transposons (CACTA family) in euchromatin, which may silence euchromatic transposons to maintain genome stability during polyploid wheat evolution. Consistently, the distribution of H3K27me2 is mutually exclusive with another repressive histone mark, H3K9me2, that mainly silences transposons in heterochromatic regions. Remarkably, the regions with low H3K27me2 levels (named H3K27me2 valleys) were associated with the formation of DNA double-strand breaks in genomes of wheat, maize (Zea mays) and Arabidopsis. Our results provide a comprehensive view of H3K27me2 and H3K27me3 distributions during wheat evolution, which support roles for H3K27me2 in silencing euchromatic transposons to maintain genome stability and in modifying genetic recombination landscapes. These genomic insights may empower breeding improvement of crops.

Analysis of Transcriptional Changes in Different Brassica napus Synthetic Allopolyploids

Allopolyploidy is an evolutionary and mechanistically intriguing process involving the reconciliation of two or more sets of diverged genomes and regulatory interactions, resulting in new phenotypes. In this study, we explored the gene expression patterns of eight F2 synthetic Brassica napus using RNA sequencing. We found that B. napus allopolyploid formation was accompanied by extensive changes in gene expression. A comparison between F2 and the parent shows a certain proportion of differentially expressed genes (DEG) and activation\silent gene, and the two genomes (female parent (AA)\male parent (CC) genomes) showed significant differences in response to whole-genome duplication (WGD); non-additively expressed genes represented a small portion, while Gene Ontology (GO) enrichment analysis showed that it played an important role in responding to WGD. Besides, genome-wide expression level dominance (ELD) was biased toward the AA genome, and the parental expression pattern of most genes showed a high degree of conservation. Moreover, gene expression showed differences among eight individuals and was consistent with the results of a cluster analysis of traits. Furthermore, the differential expression of waxy synthetic pathways and flowering pathway genes could explain the performance of traits. Collectively, gene expression of the newly formed allopolyploid changed dramatically, and this was different among the selfing offspring, which could be a prominent cause of the trait separation. Our data provide novel insights into the relationship between the expression of differentially expressed genes and trait segregation and provide clues into the evolution of allopolyploids.

To Be or Not to Be Expressed: The First Evidence of a Nucleolar Dominance Tissue-Specificity in Brachypodium hybridum

Nucleolar dominance (ND) is an epigenetic, developmentally regulated phenomenon that describes the selective inactivation of 35S rDNA loci derived from one progenitor of a hybrid or allopolyploid. The presence of ND was documented in an allotetraploid grass, Brachypodium hybridum (genome composition DDSS), which is a polyphyletic species that arose from crosses between two putative ancestors that resembled the modern B. distachyon (DD) and B. stacei (SS). In this work, we investigated the developmental stability of ND in B. hybridum genotype 3-7-2 and compared it with the reference genotype ABR113. We addressed the question of whether the ND is established in generative tissues such as pollen mother cells (PMC). We examined condensation of rDNA chromatin by fluorescence in situ hybridization employing state-of-art confocal microscopy. The transcription of rDNA homeologs was determined by reverse-transcription cleaved amplified polymorphic sequence analysis. In ABR113, the ND was stable in all tissues analyzed (primary and adventitious root, leaf, and spikes). In contrast, the 3-7-2 individuals showed a strong upregulation of the S-genome units in adventitious roots but not in other tissues. Microscopic analysis of the 3-7-2 PMCs revealed extensive decondensation of the D-genome loci and their association with the nucleolus in meiosis. As opposed, the S-genome loci were always highly condensed and localized outside the nucleolus. These results indicate that genotype-specific loss of ND in B. hybridum occurs probably after fertilization during developmental processes. This finding supports our view that B. hybridum is an attractive model to study ND in grasses.

The fate of 35S rRNA genes in the allotetraploid grass Brachypodium hybridum

Nucleolar dominance (ND) consists of the reversible silencing of 35S/45S rDNA loci inherited from one of the ancestors of an allopolyploid. The molecular mechanisms by which one ancestral rDNA set is selected for silencing remain unclear. We applied a combination of molecular (Southern blot hybridization and reverse-transcription cleaved amplified polymorphic sequence analysis), genomic (analysis of variants) and cytogenetic (fluorescence in situ hybridization) approaches to study the structure, expression and epigenetic landscape of 35S rDNA in an allotetraploid grass that exhibits ND, Brachypodium hybridum (genome composition DDSS), and its putative progenitors, Brachypodium distachyon (DD) and Brachypodium stacei (SS). In progenitor genomes, B. stacei showed a higher intragenomic heterogeneity of rDNA compared with B. distachyon. In all studied accessions of B. hybridum, there was a reduction in the copy number of S homoeologues, which was accompanied by their inactive transcriptional status. The involvement of DNA methylation in CG and CHG contexts in the silencing of the S-genome rDNA loci was revealed. In the B. hybridum allotetraploid, ND is stabilized towards the D-genome units, irrespective of the polyphyletic origin of the species, and does not seem to be influenced by homoeologous 35S rDNA ratios and developmental stage.

Homoeolog expression bias and expression level dominance (ELD) in four tissues of natural allotetraploid Brassica napus

Background

Allopolyploidy is widespread in angiosperms, and they can coordinate two or more different genomes through genetic and epigenetic modifications to exhibit stronger vigor and adaptability. To explore the changes in homologous gene expression patterns in the natural allotetraploid Brassica napus (A_nA_nC_nC_n) relative to its two diploid progenitors, B. rapa (A_rA_r) and B. oleracea (C_oC_o), after approximately 7500 years of domestication, the global gene pair expression patterns in four major tissues (stems, leaves, flowers and siliques) of these three species were analyzed using an RNA sequencing approach.

Results

The results showed that the ‘transcriptomic shock’ phenomenon was alleviated in natural B. napus after approximately 7500 years of natural domestication, and most differentially expressed genes (DEGs) in B. napus were downregulated relative to those in its two diploid progenitors. The KEGG analysis indicated that three pathways related to photosynthesis were enriched in both comparison groups (A_nA_nC_nC_n vs A_rA_r and A_nA_nC_nC_n vs C_oC_o), and these pathways were all downregulated in four tissues of B. napus. In addition, homoeolog expression bias and expression level dominance (ELD) in B. napus were thoroughly studied through analysis of expression levels of 27,609 B. rapa-B. oleracea orthologous gene pairs. The overwhelming majority of gene pairs (an average of 86.7%) in B. napus maintained their expression pattern in two diploid progenitors, and approximately 78.1% of the gene pairs showed expression bias with a preference toward the A subgenome. Overall, an average of 48, 29.7 and 22.3% homologous gene pairs exhibited additive expression, ELD and transgressive expression in B. napus, respectively. The ELD bias varies from tissue to tissue; specifically, more gene pairs in stems and siliques showed ELD-A, whereas the opposite was observed in leaves and flowers. More transgressive upregulation, rather than downregulation, was observed in gene pairs of B. napus.

Conclusions

In general, these results may provide a comprehensive understanding of the changes in homologous gene expression patterns in natural B. napus after approximately 7500 years of evolution and domestication and may enhance our understanding of allopolyploidy.

Regulation of Nucleolar Dominance in Drosophila melanogaster

In eukaryotic genomes, ribosomal RNA (rRNA) genes exist as tandemly repeated clusters, forming ribosomal DNA (rDNA) loci. Each rDNA locus typically contains hundreds of rRNA genes to meet the high demand of ribosome biogenesis. Nucleolar dominance is a phenomenon whereby individual rDNA loci are entirely silenced or transcribed, and is believed to be a mechanism to control rRNA dosage. Nucleolar dominance was originally noted to occur in interspecies hybrids, and has been shown to occur within a species (i.e., nonhybrid context). However, studying nucleolar dominance within a species has been challenging due to the highly homogenous sequence across rDNA loci. By utilizing single nucleotide polymorphisms between X rDNA and Y rDNA loci in males, as well as sequence variations between two X rDNA loci in females, we conducted a thorough characterization of nucleolar dominance throughout development of Drosophila melanogaster. We demonstrate that nucleolar dominance is a developmentally regulated program that occurs in nonhybrid, wild-type D. melanogaster, where Y rDNA dominance is established during male embryogenesis, whereas females normally do not exhibit dominance between two X rDNA loci. By utilizing various chromosomal complements (e.g., X/Y, X/X, X/X/Y) and a chromosome rearrangement, we show that the short arm of the Y chromosome including the Y rDNA likely contains information that instructs the state of nucleolar dominance. Our study begins to reveal the mechanisms underlying the selection of rDNA loci for activation/silencing in nucleolar dominance in the context of nonhybrid D. melanogaster.

Analysis on the Meiosis-Related Gene (Dmc1, Ph1) Expression in Autotriploid Carassius auratus

Triploid is usually considered to be unable to perform normal meiosis due to the abnormal behavior of the three sets of chromosomes. But autotriploid Carassius auratus in the Dongting water system (3n = 150, abbreviated as 3nCC) can perform normal meiosis. In artificial autotriploid Carassius auratus (3n = 150, abbreviated as 3nRR), female individuals undergo normal meiosis and produce mature gametes, while male individuals cannot. To better understand the effects of triploidization on meiosis in fish, we study the structure, methylation level, and expression level of meiosis-related genes (Dmc1, Ph1) in diploid Carassius auratus (2n = 100, abbreviated as 2nCC), Carassius auratus red var.(2n = 100, abbreviated as RCC), 3nCC and 3nRR. The results show that, compared with their diploid ancestors (2nCC and RCC), Dmc1 and Ph1 genes are hypomethylated in all 3nCC and female 3nRR, while are hypermethylated in male 3nRR. Correspondingly, Dmc1 and Ph1 genes are highly expressed in all 3nCC and female 3nRR, while are lowly expressed in male 3nRR. These results indicate that high expression of meiosis-related genes can contribute to restoration of bivalent pairing during meiosis in autotriploid Carassius auratus. This study provides new insights into the effect of DNA methylation on the fertility in triploid fish.

The Hordeum bulbosum 25S-18S rDNA region: comparison with Hordeum vulgare and other Triticeae

Hordeum vulgare and Hordeum bulbosum are two closely related barley species, which share a common H genome. H. vulgare has two nucleolar organizer regions (NORs), while the NOR of H. bulbosum is only one. We sequenced the 2.5 kb 25S-18S region in the rDNA of H. bulbosum and compared it to the same region in H. vulgare as well as to the other Triticeae. The region includes an intergenic spacer (IGS) with a number of subrepeats, a promoter, and an external transcribed spacer (5'ETS). The IGS of H. bulbosum downstream of 25S rRNA contains two 143-bp repeats and six 128-bp repeats. In contrast, the IGS in H. vulgare contains an array of seven 79-bp repeats and a varying number of 135-bp repeats. The 135-bp repeats in H. vulgare and the 128-bp repeats in H. bulbosum show similarity. Compared to H. vulgare, the 5'ETS of H. bulbosum is shorter. Additionally, the 5'ETS regions in H. bulbosum and H. vulgare diverged faster than in other Triticeae genera. Alignment of the Triticeae promoter sequences suggests that in Hordeum, as in diploid Triticum, transcription starts with guanine and not with adenine as it is in many other plants.

Intragenomic heterogeneity of intergenic ribosomal DNA spacers in Cucurbita moschata is determined by DNA minisatellites with variable potential to form non-canonical DNA conformations

The intergenic spacer (IGS) of rDNA is frequently built of long blocks of tandem repeats. To estimate the intragenomic variability of such knotty regions, we employed PacBio sequencing of the Cucurbita moschata genome, in which thousands of rDNA copies are distributed across a number of loci. The rRNA coding regions are highly conserved, indicating intensive interlocus homogenization and/or high selection pressure. However, the IGS exhibits high intragenomic structural diversity. Two repeated blocks, R1 (300-1250 bp) and R2 (290-643 bp), account for most of the IGS variation. They exhibit minisatellite-like features built of multiple periodically spaced short GC-rich sequence motifs with the potential to adopt non-canonical DNA conformations, G-quadruplex-folded and left-handed Z-DNA. The mutual arrangement of these motifs can be used to classify IGS variants into five structural families. Subtle polymorphisms exist within each family due to a variable number of repeats, suggesting the coexistence of an enormous number of IGS variants. The substantial length and structural heterogeneity of IGS minisatellites suggests that the tempo of their divergence exceeds the tempo of the homogenization of rDNA arrays. As frequently occurring among plants, we hypothesize that their instability may influence transcription regulation and/or destabilize rDNA units, possibly spreading them across the genome.

Allele-specific expression variation at different ploidy levels in Squalius alburnoides

Allopolyploid plants are long known to be subject to a homoeolog expression bias of varying degree. The same phenomenon was only much later suspected to occur also in animals based on studies of single selected genes in an allopolyploid vertebrate, the Iberian fish Squalius alburnoides. Consequently, this species became a good model for understanding the evolution of gene expression regulation in polyploid vertebrates. Here, we analyzed for the first time genome-wide allele-specific expression data from diploid and triploid hybrids of S. alburnoides and compared homoeolog expression profiles of adult livers and of juveniles. Co-expression of alleles from both parental genomic types was observed for the majority of genes, but with marked homoeolog expression bias, suggesting homoeolog specific reshaping of expression level patterns in hybrids. Complete silencing of one allele was also observed irrespective of ploidy level, but not transcriptome wide as previously speculated. Instead, it was found only in a restricted number of genes, particularly ones with functions related to mitochondria and ribosomes. This leads us to hypothesize that allelic silencing may be a way to overcome intergenomic gene expression interaction conflicts, and that homoeolog expression bias may be an important mechanism in the achievement of sustainable genomic interactions, mandatory to the success of allopolyploid systems, as in S. alburnoides.

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.

Differential expression of HPG-axis genes in autotetraploids derived from red crucian carp Carassius auratus red var., ♀ × blunt snout bream Megalobrama amblycephala, ♂

Autotetraploid fish (4n = 200, abbreviated as 4nRR), which reach sexual maturity at 1 year of age, were derived from the whole genome duplication of red crucian carp Carassius auratus red var. (RCC; 2n = 100) and possess four sets of chromosomes from RCC. The histological features of the gonads showed that the RCC and 4nRR both possessed normal gonadal structure and could arrive at maturation. To understand the expression characteristics of genes related to reproductive development in the autotetraploid fish, we analysed the nucleotide sequence and expression characteristics of the gnrh2, gthb and gthr genes, which are the pivotal genes of the hypothalamic-pituitary-gonadal (HPG) axis. We found that the gnrh2, gthb and gthr genes in 4nRR share remarkable homology with RCC, but there were obvious differences in expression levels between 4nRR and RCC. These results demonstrate that autotetraploidization can lead to gene expression changes. This study provides insights into the molecular mechanism underlying the reproductive development of autotetraploid fish and is expected to be of great significance for subsequent research on polyploidization.

Coregulation of ribosomal RNA with hundreds of genes contributes to phenotypic variation

Ribosomal repeats occupy 5% of a plant genome, yet there has been little study of their diversity in the modern age of genomics. Ribosomal copy number and expression variation present an opportunity to tap a novel source of diversity. In the present study, we estimated the ribosomal DNA (rDNA) copy number and ribosomal RNA (rRNA) expression for a population of maize inbred lines and investigated the potential role of rDNA and rRNA dosage in regulating global gene expression. Extensive variation was found in both ribosomal DNA copy number and ribosomal RNA expression among maize inbred lines. However, rRNA abundance was not consistent with the copy number of the rDNA. We have not found that the rDNA gene dosage has a regulatory role in gene expression; however, thousands of genes are identified to be coregulated with rRNA expression, including genes participating in ribosome biogenesis and other functionally relevant pathways. We further investigated the potential roles of copy number and the expression level of rDNA on agronomic traits and found that both correlated with flowering time but through different regulatory mechanisms. This comprehensive analysis suggested that rRNA expression variation is a valuable source of functional diversity that affects gene expression variation and field-based phenotypic changes.

Homoeolog expression bias and expression level dominance in resynthesized allopolyploid Brassica napus

Background

Allopolyploids require rapid genetic and epigenetic modifications to reconcile two or more sets of divergent genomes. To better understand the fate of duplicate genes following genomic mergers and doubling during allopolyploid formation, in this study, we explored the global gene expression patterns in resynthesized allotetraploid Brassica napus (AACC) and its diploid parents B. rapa (AA) and B. oleracea (CC) using RNA sequencing of leaf transcriptomes.

Results

We found that allopolyploid B. napus formation was accompanied by extensive changes (approximately one-third of the expressed genes) in the parental gene expression patterns (‘transcriptome shock’). Interestingly, the majority (85%) of differentially expressed genes (DEGs) were downregulated in the allotetraploid. Moreover, the homoeolog expression bias (relative contribution of homoeologs to the transcriptome) and expression level dominance (total expression level of both homoeologs) were thoroughly investigated by monitoring the expression of 23,766 B. oleracea-B. rapa orthologous gene pairs. Approximately 36.5% of the expressed gene pairs displayed expression bias with a slight preference toward the A-genome. In addition, 39.6, 4.9 and 9.0% of the expressed gene pairs exhibited expression level dominance (ELD), additivity expression and transgressive expression, respectively. The genome-wide ELD was also biased toward the A-genome in the resynthesized B. napus. To explain the ELD phenomenon, we compared the individual homoeolog expression levels relative to those of the diploid parents and found that ELD in the direction of the higher-expression parent can be explained by the downregulation of homoeologs from the dominant parent or upregulation of homoeologs from the nondominant parent; however, ELD in the direction of the lower-expression parent can be explained only by the downregulation of the nondominant parent or both homoeologs. Furthermore, Gene Ontology (GO) enrichment analysis suggested that the alteration in the gene expression patterns could be a prominent cause of the phenotypic variation between the newly formed B. napus and its parental species.

Conclusions

Collectively, our data provide insight into the rapid repatterning of gene expression at the beginning of Brassica allopolyploidization and enhance our knowledge of allopolyploidization processes.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4966-5) contains supplementary material, which is available to authorized users.

Development of Brassica oleracea-nigra monosomic alien addition lines: genotypic, cytological and morphological analyses

We report the development and characterization of Brassica oleracea - nigra monosomic alien addition lines (MAALs) to dissect the Brassica B genome. Brassica nigra (2n = 16, BB) represents the diploid Brassica B genome which carries many useful genes and traits for breeding but received limited studies. To dissect the B genome from B. nigra, the triploid F₁ hybrid (2n = 26, CCB) obtained previously from the cross B. oleracea var. alboglabra (2n = 18, CC) × B. nigra was used as the maternal parent and backcrossed successively to parental B. oleracea. The progenies in BC₁ to BC₃ generations were analyzed by the methods of FISH and SSR markers to screen the monosomic alien addition lines (MAALs) with each of eight different B-genome chromosomes added to C genome (2n = 19, CC + 1B_1-8), and seven different MAALs were established, except for the one with chromosome B2 which existed in one triple addition. Most of these MAALs were distinguishable morphologically from each other, as they expressed the characters from B. nigra differently and at variable extents. The alien chromosome remained unpaired as a univalent in 86.24% pollen mother cells at diakinesis or metaphase I, and formed a trivalent with two C-genome chromosomes in 13.76% cells. Transmission frequency of all the added chromosomes was far higher through the ovules (averagely 14.40%) than the pollen (2.64%). The B1, B4 and B5 chromosomes were transmitted by female at much higher rates (22.38-30.00%) than the other four (B3, B6, B7, B8) (5.04-8.42%). The MAALs should be valuable for exploiting the genome structure and evolution of B. nigra.

Dominance and Sexual Dimorphism Pervade the Salix purpurea L. Transcriptome

The heritability of gene expression is critical in understanding heterosis and is dependent on allele-specific regulation by local and remote factors in the genome. We used RNA-Seq to test whether variation in gene expression among F1 and F2 intraspecific Salix purpurea progeny is attributable to cis- and trans-regulatory divergence. We assessed the mode of inheritance based on gene expression levels and allele-specific expression for F1 and F2 intraspecific progeny in two distinct tissue types: shoot tip and stem internode. In addition, we explored sexually dimorphic patterns of inheritance and regulatory divergence among F1 progeny individuals. We show that in S. purpurea intraspecific crosses, gene expression inheritance largely exhibits a maternal dominant pattern, regardless of tissue type or pedigree. A significantly greater number of cis- and trans-regulated genes coincided with upregulation of the maternal parent allele in the progeny, irrespective of the magnitude, whereas the paternal allele was higher expressed for genes showing cis × trans or compensatory regulation. Importantly, consistent with previous genetic mapping results for sex in shrub willow, we have delimited sex-biased gene expression to a 2 Mb pericentromeric region on S. purpurea chr15 and further refined the sex determination region. Altogether, our results offer insight into the inheritance of gene expression in S. purpurea as well as evidence of sexually dimorphic expression which may have contributed to the evolution of dioecy in Salix.

Epigenetics of the preferential silencing of Brachypodium stacei-originated 35S rDNA loci in the allotetraploid grass Brachypodium hybridum

Nucleolar dominance (ND), initially described as ‘differential amphiplasty’, is a phenomenon observed in some plant and animal allopolyploids and hybrids in which the selective suppression of the activity of 35S rRNA gene loci that have been inherited from one of the two or more ancestral genomes occurs. Although more than 80 years have passed since the discovery of ND, there is still a significant lack in our understanding of the mechanisms that determine this phenomenon. Here, we aimed to investigate the epigenetic status of 35S rRNA gene loci in the monocotyledonous Brachypodium hybridum, which is an allotetraploid that has resulted from a cross between B. distachyon and B. stacei. We revealed that the repressed B. stacei-inherited rDNA loci are characterised by a high level of DNA methylation. The global hypomethylation of B. hybridum nuclear DNA induced by 5-azacytidine, however, seems to be insufficient for the transcriptional reactivation of these loci, which indicates that factors other than DNA methylation are behind the suppression of B. stacei-originated loci. We also showed that the transcriptionally active and silenced fractions of rRNA genes that had been inherited from B. distachyon occupy different domains within the chromocentres adjacent to the nucleolus, depending on their epigenetic status.

MOLECULAR CHARACTERIZATION OF PgiC IN A TETRAPLOID PLANT AND ITS DIPLOID RELATIVES

Allotetraploid plant species combine the genomes of related diploid species, but little is known about whether homologous genes from the diploid genomes are expressed, how they interact, or whether they evolve differently when in a common tetraploid nucleus. Polyploidy may lead to gene silencing, but few molecular characterizations of silenced genes encoding enzymes in polyploids and related diploids have been reported. Here we describe the PgiC genes in the tetraploid Clarkia gracilis and related diploid species, which are native from California to southern Washington. PgiC encodes the cytosolic isozyme of phosphoglucose isomerase (PGIC; EC 5.3.1.9). The gene was duplicated in the basal stock of Clarkia and now both genes, PgiC1 and PgiC2, are active in about half of the diploid species, whereas only PgiC1 is active in the others. Clarkia gracilis was found to have three PgiC genes: two PgiC1s and a PgiC2. Reverse-transcriptase-polymerase chain reaction (RT-PCR) experiments, starting with mRNAs prepared from seedling leaves of C. gracilis, showed that the three genes are expressed. Analysis of their sequences showed they are evolving at similar rates to their homologues and that they have the same intron-exon structure. The presence of an expressed PgiC2 in C. gracilis was unexpected because all related diploid species, including one identified as a parent, have only active PgiC1s. The donor of the PgiC2 is now presumed extinct, but parsimony analysis identified its phylogenetic position. None of the PgiC genes that were active when C. gracilis arose were silenced. A possible example of gene conversion involving a 300-nuclectide region of one PgiC1 and PgiC2 was identified, but it probably occurred in the diploid parental species rather than in C. gracilis. PgiC2 is the first known example of an active locus in a tetraploid plant species that is no longer expressed in its diploid relatives.

Epistatic and allelic interactions control expression of ribosomal RNA gene clusters in Arabidopsis thaliana

BACKGROUND

Ribosomal RNA (rRNA) accounts for the majority of the RNA in eukaryotic cells, and is encoded by hundreds to thousands of nearly identical gene copies, only a subset of which are active at any given time. In Arabidopsis thaliana, 45S rRNA genes are found in two large ribosomal DNA (rDNA) clusters and little is known about the contribution of each to the overall transcription pattern in the species.

RESULTS

By taking advantage of genome sequencing data from the 1001 Genomes Consortium, we characterize rRNA gene sequence variation within and among accessions. Notably, variation is not restricted to the pre-rRNA sequences removed during processing, but it is also present within the highly conserved ribosomal subunits. Through linkage mapping we assign these variants to a particular rDNA cluster unambiguously and use them as reporters of rDNA cluster-specific expression. We demonstrate that rDNA cluster-usage varies greatly among accessions and that rDNA cluster-specific expression and silencing is controlled via genetic interactions between entire rDNA cluster haplotypes (alleles).

CONCLUSIONS

We show that rRNA gene cluster expression is controlled via complex epistatic and allelic interactions between rDNA haplotypes that apparently regulate the entire rRNA gene cluster. Furthermore, the sequence polymorphism we discovered implies that the pool of rRNA in a cell may be heterogeneous, which could have functional consequences.

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.

Nucleolar Dominance and Repression of 45S Ribosomal RNA Genes in Hybrids between Xenopus borealis and X. muelleri (2n = 36)

Nucleolar dominance is a dramatic disruption in the formation of nucleoli and the expression of ribosomal RNA (rRNA) genes, characteristic of some plant and animal hybrids. Here, we report that F1 hybrids produced from reciprocal crosses between 2 sister species of Xenopus clawed frogs, X. muelleri and X. borealis, undergo nucleolar dominance somewhat distinct from a pattern previously reported in hybrids between phylogenetically more distant Xenopus species. Patterns of nucleolar development, 45S rRNA expression, and gene copy inheritance were investigated using a combination of immunostaining, pyrosequencing, droplet digital PCR, flow cytometry, and epigenetic inhibition. In X. muelleri × X. borealis hybrids, typically only 1 nucleolus is formed, and 45S rRNA genes are predominantly expressed from 1 progenitor's alleles, X. muelleri, regardless of the cross-direction. These changes are accompanied by an extensive (∼80%) loss of rRNA gene copies in the hybrids relative to their parents, with the transcriptionally underdominant variant (X. borealis) being preferentially lost. Chemical treatment of hybrid larvae with a histone deacetylase inhibitor resulted in a partial derepression of the underdominant variant. Together, these observations shed light on the genetic and epigenetic basis of nucleolar dominance as an underappreciated manifestation of genetic conflicts within a hybrid genome.

Rapid genomic changes in allopolyploids of Carassius auratus red var. (♀) × Megalobrama amblycephala (♂)

To better understand genomic changes in the early generations after polyploidisation, we examined the chromosomal consequences of genomic merger in allotetraploid hybrids (4 nF₁) (AABB, 4n = 148) of Carassius auratus red var. (RCC) (AA, 2n = 100) (♀) × Megalobrama amblycephala (BSB) (BB, 2n = 48) (♂). Complete loss of the paternal 5S rDNA sequence and the expected number of maternal chromosomal loci were found in 4 nF₁, suggesting directional genomic changes occurred in the first generations after polyploidisation. Recent studies have reported instability of newly established allotetraploid genomes. To assess this in the newly formed 4 nF₁ genome, we performed fluorescence in situ hybridisation on an allotetraploid gynogenetic hybrid (4 nG) (AABB, 4n = 148) and an allopentaploid hybrid (5 nH) (AABBB, 5n = 172) from 4 nF₁ (♀) × BSB (♂) with 5S rDNA gene and centromere probes from RCC, the original diploid parent. The expected numbers of maternal chromosomal loci were found in 4 nG, while chromosomal locus deletions and chromosome recombinations were detected in 5 nH. These observations suggest that abnormal meiosis did not lead to obvious genomic changes in the newly established allotetraploid genomes, but hybridisation with the original diploid parent resulted in obvious genomic changes in the newly established allotetraploid genomes, as was found for the maternal genome.

Genome-Wide Gene/Genome Dosage Imbalance Regulates Gene Expressions in Synthetic Brassica napus and Derivatives (AC, AAC, CCA, CCAA)

Gene/genome dosage balance is an essential evolutionary mechanism for organisms to ensure a normal function, but the underlying causes of dosage-imbalance regulation remain poorly understood. Herein, the serial Brassica hybrids/polyploids (AC, AAC, CCA, CCAA) with different copies of A and C subgenomes from the same two parents of Brassica rapa and Brassica oleracea were synthesized to investigate the effects of genome dosages on gene expressions and interactions by using RNA-Seq. The expression changes of A- and C-subgenome genes were consistent with dosage alterations. Dosage-dependent and -independent genes were grouped according to the correlations between dosage variations and gene expressions. Expression levels of dosage-dependent genes were strongly correlated with dosage changes and mainly contributed to dosage effects, while those of dosage-independent genes gave weak correlations with dosage variations and mostly facilitated dosage compensation. More protein-protein interactions were detected for dosage-independent genes than dosage-dependent ones, as predicted by the dosage balance hypothesis. Dosage-dependent genes more likely impacted the expressions by trans effects, whereas dosage-independent genes preferred to play by cis effects. Furthermore, dosage-dependent genes were mainly associated with the basic biological processes to maintain the stability of the growth and development, while dosage-independent genes were more enriched in the stress response related processes to accelerate adaptation. The present comprehensive analysis of gene expression dependent/independent on dosage alterations in Brassica polyploids provided new insights into gene/genome dosage-imbalance regulation of gene expressions.

Analysis of Different Ploidy and Parent-Offspring Genomic DNA Methylation in the Loach Misgurnus anguillicaudatus

In this study, we selected natural polyploidy loach (diploid, triploid and tetraploid) and hybrid F₁ generation obverse cross (4 × 2) and inverse cross (2 × 4) by diploids and tetraploids as the research model. The MSAP (methylation-sensitive amplified polymorphism) reaction system was established by our laboratory to explore methylation levels and pattern diversification features at the whole genome level of the polyploidy loach. The results showed that the total methylation and full methylation rates decreased on increased ploidy individuals; moreover, the hemimethylation rate showed no consistent pattern. Compared with diploid loach, the methylation patterns of tetraploid sites changed 68.17%, and the methylation patterns of triploid sites changed 73.05%. The proportion of hypermethylation genes is significantly higher than the proportion of demethylation genes. The methylation level of reciprocal cross F₁ generation is lower than the male diploid and higher than the female tetraploid. The hemimethylation and total methylation rate of the cross hybrid F₁ generation is significantly higher than the orthogonal F₁ generation (p < 0.01). After readjusting, the methylation pattern of genome DNA of reciprocal hybrids changed 69.59% and 72.83%, respectively.

HANDS2: accurate assignment of homoeallelic base-identity in allopolyploids despite missing data

Characterization of homoeallelic base-identity in allopolyploids is difficult since homeologous subgenomes are closely related and becomes further challenging if diploid-progenitor data is missing. We present HANDS2, a next-generation sequencing-based tool that enables highly accurate (>90%) genome-wide discovery of homeolog-specific base-identity in allopolyploids even in the absence of a diploid-progenitor. We applied HANDS2 to the transcriptomes of various cruciferous plants belonging to genus Brassica. Our results suggest that the three C genomes in Brassica are more similar to each other than the three A genomes, and provide important insights into the relationships between various Brassica tetraploids and their diploid-progenitors at a single-base resolution.

Genetic and Epigenetic Changes in Oilseed Rape (Brassica napus L.) Extracted from Intergeneric Allopolyploid and Additions with Orychophragmus

Allopolyploidization with the merger of the genomes from different species has been shown to be associated with genetic and epigenetic changes. But the maintenance of such alterations related to one parental species after the genome is extracted from the allopolyploid remains to be detected. In this study, the genome of Brassica napus L. (2n = 38, genomes AACC) was extracted from its intergeneric allohexaploid (2n = 62, genomes AACCOO) with another crucifer Orychophragmus violaceus (2n = 24, genome OO), by backcrossing and development of alien addition lines. B. napus-type plants identified in the self-pollinated progenies of nine monosomic additions were analyzed by the methods of amplified fragment length polymorphism, sequence-specific amplified polymorphism, and methylation-sensitive amplified polymorphism. They showed modifications to certain extents in genomic components (loss and gain of DNA segments and transposons, introgression of alien DNA segments) and DNA methylation, compared with B. napus donor. The significant differences in the changes between the B. napus types extracted from these additions likely resulted from the different effects of individual alien chromosomes. Particularly, the additions which harbored the O. violaceus chromosome carrying dominant rRNA genes over those of B. napus tended to result in the development of plants which showed fewer changes, suggesting a role of the expression levels of alien rRNA genes in genomic stability. These results provided new cues for the genetic alterations in one parental genome that are maintained even after the genome becomes independent.

Polyploidy and the proteome

Although major advances have been made during the past 20 years in our understanding of the genetic and genomic consequences of polyploidy, our knowledge of polyploidy and the proteome is in its infancy. One of our goals is to stimulate additional study, particularly broad-scale proteomic analyses of polyploids and their progenitors. Although it may be too early to generalize regarding the extent to which transcriptomic data are predictive of the proteome of polyploids, it is clear that the proteome does not always reflect the transcriptome. Despite limited data, important observations on the proteomes of polyploids are emerging. In some cases, proteomic profiles show qualitatively and/or quantitatively non-additive patterns, and proteomic novelty has been observed. Allopolyploids generally combine the parental contributions, but there is evidence of parental dominance of one contributing genome in some allopolyploids. Autopolyploids are typically qualitatively identical to but quantitatively different from their parents. There is also evidence of parental legacy at the proteomic level. Proteomes clearly provide insights into the consequences of genomic merger and doubling beyond what is obtained from genomic and/or transcriptomic data. Translating proteomic changes in polyploids to differences in morphology and physiology remains the holy grail of polyploidy--this daunting task of linking genotype to proteome to phenotype should emerge as a focus of polyploidy research in the next decade. This article is part of a Special Issue entitled: Plant Proteomics--a bridge between fundamental processes and crop production, edited by Dr. Hans-Peter Mock.

Interpopulation hybridization generates meiotically stable rDNA epigenetic variants in allotetraploid Tragopogon mirus

Uniparental silencing of 35S rRNA genes (rDNA), known as nucleolar dominance (ND), is common in interspecific hybrids. Allotetraploid Tragopogon mirus composed of Tragopogon dubius (d) and Tragopogon porrifolius (p) genomes shows highly variable ND. To examine the molecular basis of such variation, we studied the genetic and epigenetic features of rDNA homeologs in several lines derived from recently and independently formed natural populations. Inbred lines derived from T. mirus with a dominant d-rDNA homeolog transmitted this expression pattern over generations, which may explain why it is prevalent among natural populations. In contrast, lines derived from the p-rDNA dominant progenitor were meiotically unstable, frequently switching to co-dominance. Interpopulation crosses between progenitors displaying reciprocal ND resulted in d-rDNA dominance, indicating immediate suppression of p-homeologs in F1 hybrids. Original p-rDNA dominance was not restored in later generations, even in those segregants that inherited the corresponding parental rDNA genotype, thus indicating the generation of additional p-rDNA and d-rDNA epigenetic variants. Despite preserved intergenic spacer (IGS) structure, they showed altered cytosine methylation and chromatin condensation patterns, and a correlation between expression, hypomethylation of RNA Pol I promoters and chromatin decondensation was apparent. Reversion of such epigenetic variants occurred rarely, resulting in co-dominance maintained in individuals with distinct genotypes. Generally, interpopulation crosses may generate epialleles that are not present in natural populations, underlying epigenetic dynamics in young allopolyploids. We hypothesize that highly expressed variants with distinct IGS features may induce heritable epigenetic reprogramming of the partner rDNA arrays, harmonizing the expression of thousands of genes in allopolyploids.

Evidence of the occurrence of structural chromosome changes at the initial diploid diversification of the autopolyploid Turnera sidoides L. (Passifloraceae) complex

Turnera sidoides is an autopolyploid complex of obligate outcrossing perennial herbs. It includes five subspecies and five morphotypes in which diploid to octoploid cytotypes were found. Based on phenetic analyses of the complex and karyotype data of polyploid cytotypes, it has been hypothesized that morphological and chromosome differentiation of T. sidoides occurred at the diploid level. To test this hypothesis, we present the first detailed chromosome analysis of diploid populations of three subspecies and four morphotypes. CMA(+)/DAPI(-) bands were restricted to secondary constrictions (except in the andino morphotype) and varied in number and position among taxa. By contrast, DAPI staining was uniform in all the materials investigated. The number and position of 45S rDNA loci were coincident with the CMA(+)/DAPI(-) bands associated with secondary constrictions. Only one pair of 5S rDNA loci was detected in all the taxa (except in subsp. holosericea), but its position was variable. The identified chromosome markers varied among the three subspecies analyzed, but they were more conserved among the morphotypes of subsp. pinnatifida. Cluster analysis of these chromosome markers supports the current taxonomic arrangement of diploids and demonstrates that structural chromosome changes would have led or accompanied the initial differentiation of T. sidoides at the diploid level.

Cytomolecular Analysis of Ribosomal DNA Evolution in a Natural Allotetraploid Brachypodium hybridum and Its Putative Ancestors-Dissecting Complex Repetitive Structure of Intergenic Spacers

Nucleolar dominance is an epigenetic phenomenon associated with nuclear 35S rRNA genes and consists in selective suppression of gene loci inherited from one of the progenitors in the allopolyploid. Our understanding of the exact mechanisms that determine this process is still fragmentary, especially in case of the grass species. This study aimed to shed some light on the molecular basis of this genome-specific inactivation of 35S rDNA loci in an allotetraploid Brachypodium hybridum (2n = 30), which arose from the interspecific hybridization between two diploid ancestors that were very similar to modern B. distachyon (2n = 10) and B. stacei (2n = 20). Using fluorescence in situ hybridization with 25S rDNA and chromosome-specific BAC clones as probes we revealed that the nucleolar dominance is present not only in meristematic root-tip cells but also in differentiated cell fraction of B. hybridum. Additionally, the intergenic spacers (IGSs) from both of the putative ancestors and the allotetraploid were sequenced and analyzed. The presumptive transcription initiation sites, spacer promoters and repeated elements were identified within the IGSs. Two different length variants, 2.3 and 3.5 kb, of IGSs were identified in B. distachyon and B. stacei, respectively, however only the IGS that had originated from B. distachyon-like ancestor was present in the allotetraploid. The amplification pattern of B. hybridum IGSs suggests that some genetic changes occurred in inactive B. stacei-like rDNA loci during the evolution of the allotetraploid. We hypothesize that their preferential silencing is an effect of structural changes in the sequence rather than just the result of the sole inactivation at the epigenetic level.

Cytonuclear Coordination Is Not Immediate upon Allopolyploid Formation in Tragopogon miscellus (Asteraceae) Allopolyploids

Allopolyploids, formed by hybridization and chromosome doubling, face the immediate challenge of having duplicated nuclear genomes that interact with the haploid and maternally inherited cytoplasmic (plastid and mitochondrial) genomes. Most of our knowledge of the genomic consequences of allopolyploidy has focused on the fate of the duplicated nuclear genes without regard to their potential interactions with cytoplasmic genomes. As a step toward understanding the fates of nuclear-encoded subunits that are plastid-targeted, here we examine the retention and expression of the gene encoding the small subunit of Ribulose-1, 5-bisphosphate carboxylase/oxygenase (Rubisco; rbcS) in multiple populations of allotetraploid Tragopogon miscellus (Asteraceae). These polyploids formed recently (~80 years ago) and repeatedly from T. dubius and T. pratensis in the northwestern United States. Examination of 79 T. miscellus individuals from 10 natural populations, as well as 25 synthetic allotetraploids, including reciprocally formed plants, revealed a low percentage of naturally occurring individuals that show a bias in either gene (homeolog) loss (12%) or expression (16%), usually toward maintaining the maternal nuclear copy of rbcS. For individuals showing loss, seven retained the maternally derived rbcS homeolog only, while three had the paternally derived copy. All of the synthetic polyploid individuals examined (S0 and S1 generations) retained and expressed both parental homeologs. These results demonstrate that cytonuclear coordination does not happen immediately upon polyploid formation in Tragopogon miscellus.