Highly preserved roles of Brassica MIR172 in polyploid Brassicas: ectopic expression of variants of Brassica MIR172 accelerates floral transition

Single-cell transcriptome analysis dissects lncRNA-associated gene networks in Arabidopsis

The plant genome produces an extremely large collection of long noncoding RNAs (lncRNAs) that are generally expressed in a context-specific manner and have pivotal roles in regulation of diverse biological processes. Here, we mapped the transcriptional heterogeneity of lncRNAs and their associated gene regulatory networks at single-cell resolution. We generated a comprehensive cell atlas at the whole-organism level by integrative analysis of 28 published single-cell RNA sequencing (scRNA-seq) datasets from juvenile Arabidopsis seedlings. We then provided an in-depth analysis of cell-type-related lncRNA signatures that show expression patterns consistent with canonical protein-coding gene markers. We further demonstrated that the cell-type-specific expression of lncRNAs largely explains their tissue specificity. In addition, we predicted gene regulatory networks on the basis of motif enrichment and co-expression analysis of lncRNAs and mRNAs, and we identified putative transcription factors orchestrating cell-type-specific expression of lncRNAs. The analysis results are available at the single-cell-based plant lncRNA atlas database (scPLAD; https://biobigdata.nju.edu.cn/scPLAD/). Overall, this work demonstrates the power of integrative single-cell data analysis applied to plant lncRNA biology and provides fundamental insights into lncRNA expression specificity and associated gene regulation.

Transgressive and parental dominant gene expression and cytosine methylation during seed development in Brassica napus hybrids

KEY MESSAGE

Transcriptomic and epigenomic profiling of gene expression and small RNAs during seed and seedling development reveals expression and methylation dominance levels with implications on early stage heterosis in oilseed rape. The enhanced performance of hybrids through heterosis remains a key aspect in plant breeding; however, the underlying mechanisms are still not fully elucidated. To investigate the potential role of transcriptomic and epigenomic patterns in early expression of hybrid vigor, we investigated gene expression, small RNA abundance and genome-wide methylation in hybrids from two distant Brassica napus ecotypes during seed and seedling developmental stages using next-generation sequencing. A total of 31117, 344, 36229 and 7399 differentially expressed genes, microRNAs, small interfering RNAs and differentially methylated regions were identified, respectively. Approximately 70% of the differentially expressed or methylated features displayed parental dominance levels where the hybrid followed the same patterns as the parents. Via gene ontology enrichment and microRNA-target association analyses during seed development, we found copies of reproductive, developmental and meiotic genes with transgressive and paternal dominance patterns. Interestingly, maternal dominance was more prominent in hypermethylated and downregulated features during seed formation, contrasting to the general maternal gamete demethylation reported during gametogenesis in angiosperms. Associations between methylation and gene expression allowed identification of putative epialleles with diverse pivotal biological functions during seed formation. Furthermore, most differentially methylated regions, differentially expressed siRNAs and transposable elements were in regions that flanked genes without differential expression. This suggests that differential expression and methylation of epigenomic features may help maintain expression of pivotal genes in a hybrid context. Differential expression and methylation patterns during seed formation in an F₁ hybrid provide novel insights into genes and mechanisms with potential roles in early heterosis.

Identification of microRNA and analysis of target genes in Panax ginseng

Objective

Ginsenosides, polysaccharides and phenols, the main active ingredients in Panax ginseng, are not different significantly in content between 3 and 5 years old of ginsengs called Yuan ginseng and more than ten years old ones called Shizhu ginseng. The responsible chemical compounds cannot fully explain difference in efficacy between them. According to reports in Lonicerae Japonicae Flos (Jinyinhua in Chinese) and Glycyrrhizae Radix et Rhizoma (Gancao in Chinese), microRNA may play a role in efficacy, so we identified microRNAs in P. ginseng at the different growth years and analyzed their target genes.

Methods

Using high-throughput sequencing, the RNA-Seq, small RNA-Seq and degradome databases of P. ginseng were constructed. The differentially expressed microRNAs was identified by qRT-PCR.

Results

A total of 63,875 unigenes and 24,154,579 small RNA clean reads were obtained from the roots of P. ginseng. From these small RNAs, 71 miRNA families were identified by bioinformatics target prediction software, including 34 conserved miRNAs, 37 non-conserved miRNA families, as well as 179 target genes of 17 known miRNAs. Through degradome sequencing and computation, we finally verified 13 targets of eight miRNAs involved in transcription, energy metabolism, biological stress and disease resistance, suggesting the significance of miRNAs in the development of P. ginseng. Consistently, major miRNA targets exhibited tissue specificity and complexity in expression patterns.

Conclusion

Differential expression microRNAs were found in different growth years of ginsengs (Shizhu ginseng and Yuan ginseng), and the regulatory roles and functional annotations of miRNA targets in P. ginseng need further investigation.

A prescient evolutionary model for genesis, duplication and differentiation of MIR160 homologs in Brassicaceae

MicroRNA160 is a class of nitrogen-starvation responsive genes which governs establishment of root system architecture by down-regulating AUXIN RESPONSE FACTOR genes (ARF10, ARF16 and ARF17) in plants. The high copy number of MIR160 variants discovered by us from land plants, especially polyploid crop Brassicas, posed questions regarding genesis, duplication, evolution and function. Absence of studies on impact of whole genome and segmental duplication on retention and evolution of MIR160 homologs in descendent plant lineages prompted us to undertake the current study. Herein, we describe ancestry and fate of MIR160 homologs in Brassicaceae in context of polyploidy driven genome re-organization, copy number and differentiation. Paralogy amongst Brassicaceae MIR160a, MIR160b and MIR160c was inferred using phylogenetic analysis of 468 MIR160 homologs from land plants. The evolutionarily distinct MIR160a was found to represent ancestral form and progenitor of MIR160b and MIR160c. Chronology of evolutionary events resulting in origin and diversification of genomic loci containing MIR160 homologs was delineated using derivatives of comparative synteny. A prescient model for causality of segmental duplications in establishment of paralogy in Brassicaceae MIR160, with whole genome duplication accentuating the copy number increase, is being posited in which post-segmental duplication events viz. differential gene fractionation, gene duplications and inversions are shown to drive divergence of chromosome segments. While mutations caused the diversification of MIR160a, MIR160b and MIR160c, duplicated segments containing these diversified genes suffered gene rearrangements via gene loss, duplications and inversions. Yet the topology of phylogenetic and phenetic trees were found congruent suggesting similar evolutionary trajectory. Over 80% of Brassicaceae genomes and subgenomes showed a preferential retention of single copy each of MIR160a, MIR160b and MIR160c suggesting functional relevance. Thus, our study provides a blue-print for reconstructing ancestry and phylogeny of MIRNA gene families at genomics level and analyzing the impact of polyploidy on organismal complexity. Such studies are critical for understanding the molecular basis of agronomic traits and deploying appropriate candidates for crop improvement.

Interplay between miRNAs and lncRNAs: Mode of action and biological roles in plant development and stress adaptation

Plants employ sophisticated mechanisms to control developmental processes and to cope with environmental changes at transcriptional and post-transcriptional levels. MicroRNAs (miRNAs) and long noncoding RNAs (lncRNAs), two classes of endogenous noncoding RNAs, are key regulators of gene expression in plants. Recent studies have identified the interplay between miRNAs and lncRNAs as a novel regulatory layer of gene expression in plants. On one hand, miRNAs target lncRNAs for the production of phased small interfering RNAs (phasiRNAs). On the other hand, lncRNAs serve as origin of miRNAs or regulate the accumulation or activity of miRNAs at transcription and post-transcriptional levels. Theses lncRNA-miRNA interplays are crucial for plant development, physiology and responses to biotic and abiotic stresses. In this review, we summarize recent advances in the biological roles, interaction mechanisms and computational predication methods of the interplay between miRNAs and lncRNAs in plants.

Homeologs of Brassica SOC1, a central regulator of flowering time, are differentially regulated due to partitioning of evolutionarily conserved transcription factor binding sites in promoters

Evolution of Brassica genome post-polyploidization reveals asymmetrical genome fractionation and copy number variation. Herein, we describe the impact of promoter divergence among SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1) homeologs on expression and function in Brassica spp. SOC1, a regulated floral pathway integrator, is conserved as 3 redundant homeologs in diploid Brassicas. Even with high sequence identity within coding regions (92.8-100%), the spatio-temporal expression patterns of 9 SOC1 homologs in B. juncea and B. nigra indicates regulatory divergence. While LF and MF2 SOC1 homeologs are upregulated during floral transition, MF1 is barely expressed. Also, MF2 homeolog levels do not decline post-flowering, unlike LF. To investigate the underlying source of divergence, we analyzed the sequence and phylogeny of all reported (22) and isolated (21) upstream regions of Brassica SOC1. Full length upstream regions (4712-19189 bp) reveal 5 ubiquitously conserved ancestral Blocks, harboring binding sites of 18 TFs (TFBSs) characterized in Arabidopsis thaliana. The orthologs of these TFBSs are differentially conserved among Brassica SOC1 homeologs, imparting expression divergence. No crucial TFBSs are exclusively lost from LF_SOC1 promoter, while MF1_SOC1 has lost NF-Y binding site crucial for SOC1 activation by CONSTANS. MF2_SOC1 homeologs have lost important TFBSs (SEP3, AP1 and SMZ), responsible for SOC1 repression post-flowering. BjuAALF_SOC1 promoter (proximal 2 kb) shows ubiquitous reporter expression in B. juncea cv. Varuna transgenics, while BjuAAMF1_SOC1 promoter shows absence of reporter expression, validating the impact of TFBS divergence. Conservation of the original primary protein sequence is discovered in B. rapa homeologs (46) of 18 TFs. Co-regulation pattern of these TFs appeared similar for B. rapa LF and MF2 SOC1 homeologs; MF1 shows significant variation. Strong regulatory association is recorded for AP1, AP2, SEP3, FLC and CONSTANS/NF-Y, highlighting their importance in homeolog-specific SOC1 regulation. Correlation of B. juncea AP1, AP2 and FLC expression with SOC1 homeologs also complies with the TFBS differences. We thus conclude that redundant SOC1 loci contribute differentially to cumulative expression of SOC1 due to divergent selection of ancestral TFBSs.

Besides and Beyond Flowering: Other roles of EuAP2 Genes in Plant Development

EuAP2 genes are well-known for their role in flower development, a legacy of the founding member of this subfamily of transcription factors, whose mutants lacked petals in Arabidopsis. However, studies of euAP2 genes in several species have accumulated evidence highlighting the diverse roles of euAP2 genes in other aspects of plant development. Here, we emphasize other developmental roles of euAP2 genes in various species and suggest a shift from regarding euAP2 genes as just flowering genes to consider the global role they may be playing in plant development. We hypothesize that their almost universal expression profile and pleiotropic effects of their mutation suggest their involvement in fundamental plant development processes.

miR172 Regulates both Vegetative and Reproductive Development in the Perennial Woody Plant Jatropha curcas

Jatropha curcas is a promising feedstock for biofuel production because its oil is highly suitable for processing bio-jet fuels and biodiesel. However, Jatropha exhibits a long juvenile stage in subtropical areas. miR172, a conserved small non-protein-coding RNA molecule with 21 nucleotides, regulates a wide range of developmental processes. To date, however, no studies have examined the function of miR172 in Jatropha. There are five miR172 precursors encoding two mature miR172s in Jatropha, which are expressed in all tissues, with the highest expression level in leaves, and the levels are up-regulated with age. Overexpression of JcmiR172a resulted in early flowering, abnormal flowers, and altered leaf morphology in transgenic Arabidopsis and Jatropha. The expression levels of miR172 target genes were down-regulated, and the flower identity genes were up-regulated in the JcmiR172a-overexpressing transgenic plants. Interestingly, we showed that JcmiR172 might be involved in regulation of stem vascular development through manipulating the expression of cellulose and lignin biosynthesis genes. Overexpression of JcmiR172a enhanced xylem development and reduced phloem and pith development. This study helped elucidate the functions of miR172 in perennial plants, a known age-related miRNA involved in the regulation of perennial plant phase change and organ development.