Autotetraploid rice methylome analysis reveals methylation variation of transposable elements and their effects on gene expression.

Polyploidy, or whole-genome duplication (WGD), serves as a key innovation in plant evolution and is an important genomic feature for all eukaryotes. Neopolyploids have to overcome difficulties in meiosis, genomic alterations, changes of gene expression, and epigenomic reorganization. However, the underlying mechanisms for these processes are poorly understood. One of the most interesting aspects is that genome doubling events increase the dosage of all genes. Unlike allopolyploids entangled by both hybridization and polyploidization, autopolyploids, especially artificial lines, in relatively uniform genetic background offer a model system to understand mechanisms of genome-dosage effects. To investigate DNA methylation effects in response to WGD rather than hybridization, we produced autotetraploid rice with its diploid donor, Oryza sativa ssp. indica cv. Aijiaonante, both of which were independently self-pollinated over 48 generations, and generated and compared their comprehensive transcriptomes, base pair-resolution methylomes, and siRNAomes. DNA methylation variation of transposable elements (TEs) was observed as widespread in autotetraploid rice, in which hypermethylation of class II DNA transposons was predominantly noted in CHG and CHH contexts. This was accompanied by changes of 24-nt siRNA abundance, indicating the role of the RNA-directed DNA methylation pathway. Our results showed that the increased methylation state of class II TEs may suppress the expression of neighboring genes in autotetraploid rice that has obtained double alleles, leading to no significant differences in transcriptome alterations for most genes from its diploid donor. Collectively, our findings suggest that chromosome doubling induces methylation variation in TEs that affect gene expression and may become a "genome shock" response factor to help neoautopolyploids adapt to genome-dosage effects.

Rapid diversification of five Oryza AA genomes associated with rice adaptation.

Comparative genomic analyses among closely related species can greatly enhance our understanding of plant gene and genome evolution. We report de novo-assembled AA-genome sequences for Oryza nivara, Oryza glaberrima, Oryza barthii, Oryza glumaepatula, and Oryza meridionalis. Our analyses reveal massive levels of genomic structural variation, including segmental duplication and rapid gene family turnover, with particularly high instability in defense-related genes. We show, on a genomic scale, how lineage-specific expansion or contraction of gene families has led to their morphological and reproductive diversification, thus enlightening the evolutionary process of speciation and adaptation. Despite strong purifying selective pressures on most Oryza genes, we documented a large number of positively selected genes, especially those genes involved in flower development, reproduction, and resistance-related processes. These diversifying genes are expected to have played key roles in adaptations to their ecological niches in Asia, South America, Africa and Australia. Extensive variation in noncoding RNA gene numbers, function enrichment, and rates of sequence divergence might also help account for the different genetic adaptations of these rice species. Collectively, these resources provide new opportunities for evolutionary genomics, numerous insights into recent speciation, a valuable database of functional variation for crop improvement, and tools for efficient conservation of wild rice germplasm.

De novo transcriptome sequencing of Camellia sasanqua and the analysis of major candidate genes related to floral traits.

Camellia sasanqua is one of the most famous horticultural plants in Camellia (Theaceae) due to its aesthetic appeal as landscape plant. Knowledge regarding the genetic basis of flowering time, floral aroma and color in C. sasanqua is limited, but is essential to breed new varieties with desired floral traits. Here, we described the de novo transcriptome of young leaves, flower buds and flowers of C. sasanqua. A total of 60,127 unigenes were functionally annotated based on the sequence similarity. After analysis, we found that two floral integrator genes, SOC1 and AP1, in flowering time pathway showed evidence of gene family expansion. Compared with 1-deoxy-D-xylulose-5-phosphate pathway, some genes in the mevalonate pathway were most highly expressed, suggesting that this might represent the major pathway for terpenoid biosynthesis related to floral aroma in C. sasanqua. In flavonoid biosynthesis pathway, PAL, CHI, DFR and ANS showing significantly higher expression levels in flowers and flower buds might have important role in regulation of floral color. The top five most transcription factors (TFs) families in C. sasanqua transcriptome were MYB, MIKC, C3H, FAR1 and HD-ZIP, many of which have a direct relationship with floral traits. In addition, we also identified 33,540 simple sequence repeats (SSRs) in the C. sasanqua transcriptome. Collectively, the C. sasanqua transcriptome dataset generated from this study along with the SSR markers provide a new resource for the identification of novel regulatory transcripts and will accelerate the genetic improvement of C. sasanqua breeding programs.

Transcriptome Analysis Identifies Candidate Genes Related to Triacylglycerol and Pigment Biosynthesis and Photoperiodic Flowering in the Ornamental and Oil-Producing Plant, Camellia reticulata (Theaceae).

Camellia reticulata, which is native to Southwest China, is famous for its ornamental flowers and high-quality seed oil. However, the lack of genomic information for this species has largely hampered our understanding of its key pathways related to oil production, photoperiodic flowering process and pigment biosynthesis. Here, we first sequenced and characterized the transcriptome of a diploid C. reticulata in an attempt to identify genes potentially involved in triacylglycerol biosynthesis (TAGBS), photoperiodic flowering, flavonoid biosynthesis (FlaBS), carotenoid biosynthesis (CrtBS) pathways. De novo assembly of the transcriptome provided a catalog of 141,460 unigenes with a total length of ~96.1 million nucleotides (Mnt) and an N50 of 1080 nt. Of them, 22,229 unigenes were defined as differentially expressed genes (DEGs) across five sequenced tissues. A large number of annotated genes in C. reticulata were found to have been duplicated, and differential expression patterns of these duplicated genes were commonly observed across tissues, such as the differential expression of SOC1_a, SOC1_b, and SOC1_c in the photoperiodic flowering pathway. Up-regulation of SAD_a and FATA genes and down-regulation of FAD2_a gene in the TAGBS pathway in seeds may be relevant to the ratio of monounsaturated fatty acid (MUFAs) to polyunsaturated fatty acid (PUFAs) in seed oil. MYBF1, a transcription regulator gene of the FlaBS pathway, was found with great sequence variation and alteration of expression patterns, probably resulting in functionally evolutionary differentiation in C. reticulata. MYBA1_a and some anthocyanin-specific biosynthetic genes in the FlaBS pathway were highly expressed in both flower buds and flowers, suggesting important roles of anthocyanin biosynthesis in flower development. Besides, a total of 40,823 expressed sequence tag simple sequence repeats (EST-SSRs) were identified in the C. reticulata transcriptome, providing valuable marker resources for further basic and applied researches on this economically important Camellia plant.

Genome-wide identification and comparative expression analysis reveal a rapid expansion and functional divergence of duplicated genes in the WRKY gene family of cabbage, Brassica oleracea var. capitata.

WRKY transcription factors (TFs), one of the ten largest TF families in higher plants, play important roles in regulating plant development and resistance. To date, little is known about the WRKY TF family in Brassica oleracea. Recently, the completed genome sequence of cabbage (B. oleracea var. capitata) allows us to systematically analyze WRKY genes in this species. A total of 148 WRKY genes were characterized and classified into seven subgroups that belong to three major groups. Phylogenetic and synteny analyses revealed that the repertoire of cabbage WRKY genes was derived from a common ancestor shared with Arabidopsis thaliana. The B. oleracea WRKY genes were found to be preferentially retained after the whole-genome triplication (WGT) event in its recent ancestor, suggesting that the WGT event had largely contributed to a rapid expansion of the WRKY gene family in B. oleracea. The analysis of RNA-Seq data from various tissues (i.e., roots, stems, leaves, buds, flowers and siliques) revealed that most of the identified WRKY genes were positively expressed in cabbage, and a large portion of them exhibited patterns of differential and tissue-specific expression, demonstrating that these gene members might play essential roles in plant developmental processes. Comparative analysis of the expression level among duplicated genes showed that gene expression divergence was evidently presented among cabbage WRKY paralogs, indicating functional divergence of these duplicated WRKY genes.

CandiSSR: An Efficient Pipeline used for Identifying Candidate Polymorphic SSRs Based on Multiple Assembled Sequences.

Simple sequence repeats (SSRs), also known as microsatellites, are ubiquitous short tandem duplications commonly found in genomes and/or transcriptomes of diverse organisms. They represent one of the most powerful molecular markers for genetic analysis and breeding programs because of their high mutation rate and neutral evolution. However, traditionally experimental screening of the SSR polymorphic status and their subsequent applicability to genetic studies are extremely labor-intensive and time-consuming. Thankfully, the recently decreased costs of next generation sequencing and increasing availability of large genome and/or transcriptome sequences have provided an excellent opportunity and sources for large-scale mining this type of molecular markers. However, current tools are limited. Thus we here developed a new pipeline, CandiSSR, to identify candidate polymorphic SSRs (PolySSRs) based on the multiple assembled sequences. The pipeline allows users to identify putative PolySSRs not only from the transcriptome datasets but also from multiple assembled genome sequences. In addition, two confidence metrics including standard deviation and missing rate of the SSR repetitions are provided to systematically assess the feasibility of the detected PolySSRs for subsequent application to genetic characterization. Meanwhile, primer pairs for each identified PolySSR are also automatically designed and further evaluated by the global sequence similarities of the primer-binding region, ensuring the successful rate of the marker development. Screening rice genomes with CandiSSR and subsequent experimental validation showed an accuracy rate of over 90%. Besides, the application of CandiSSR has successfully identified a large number of PolySSRs in the Arabidopsis genomes and Camellia transcriptomes. CandiSSR and the PolySSR marker sources are publicly available at: http://www.plantkingdomgdb.com/CandiSSR/index.html.