Microsynteny and phylogenetic analysis of tandemly organised miRNA families across five members of Brassicaceae reveals complex retention and loss history

The non-canonically organized members of MIR395 gene family in Brassica juncea are associated with developmentally regulated, sulfate-stress responsive bidirectional promoters that exhibit orientation-dependent differential transcriptional activity

Several MICRORNA genes belonging to same family or different families are often found in homologous or non-homologous clusters. Among the various classes, head-to-head arranged genes form one of the largest categories of non-canonically organized genes. Such head-to-head arranged, non-canonically organized genes possibly share cis-regulatory region with the intergenic sequence having the potential to function as bi-directional promoter (BDP). The transcriptional regulation of head-to-head arranged genes, especially with bidirectional promoters, remains an enigma. In the past, bidirectional promoters have been characterized for a small set of protein-coding gene pairs in plants; however, to the best of our knowledge, no such study has been carried so far for MICRORNA genes. The present study thus functionally characterizes bidirectional promoters associated with members of MIR395 family, which is evolutionary conserved and is most frequently occurring cluster across plant kingdom. In Arabidopsis thaliana, the MIR395 gene family contains six members with two head-to-head arranged gene pairs- MIR395A-B and MIR395E-F. This organization was found to be conserved at seven loci for MIR395A-B, and eleven loci for MIR395E-F in five Brassica sps. Sequence analysis of the putative bidirectional promoters revealed variation in length, GC content and distribution of strict TATA-box. Comparatively higher level of conservation at both the ends of the bidirectional promoters, corresponding to ca. 250 bp upstream of 5'end of the respective MIRNA precursor, was observed. These conserved regions harbour several abiotic stress (nutrient, salt, drought) and hormone (ABA, ethylene) responsive cis-motifs. Functional characterization of putative bidirectional promoters associated with MIR395A-B and MIR395E-F from Arabidopsis and their respective orthologs from Brassica juncea (Bj_A08 MIR395A-B, Bj_B03 MIR395A-B, Bj_A07.1 MIR395E-F and Bj_A07.2 MIR395E-F) was carried out using a dual-reporter vector with β-glucuronidase (GUS) and Green Fluorescent Protein (GFP). Analysis of transcriptional regulation of the two reporter genes - GUS and GFP during developmental stages confirmed their bidirectional nature. Orientation-dependent differential reporter activity indicated asymmetric nature of the promoters. Comparison of the reporter activity amongst orthologs, paralogs and homeologs revealed regulatory diversification, an outcome expected in polyploid genomes. Interestingly, reporter gene activities driven by selected bidirectional promoters were also observed in anther and siliques apart vegetative tissues indicating role of miR395 in anther and fruit development. Finally, we evaluated the activity of reporter genes driven under transcriptional regulation of bidirectional promoters under normal and sulfate-deprived conditions which revealed asymmetric inducibility under sulfate-starvation, in agreement with the known role of miR395 in sulfate homeostasis.

Complex origin, evolution, and diversification of non-canonically organized OVATE-OFP and OVATE-Like OFP gene pair across Embryophyta

Ovate Family Proteins (OFP) is a plant-specific gene family of negative transcriptional regulators. Till-date, a handful of in-silico studies have provided glimpses into family size, expansion patterns, and genic features across all major plant lineages. A major lacuna exists in understanding origin of organisation complexity of members such as those arranged in a head-to-head manner which may lead to transcriptional co-regulation via a common bi-directional promoter. To address this gap, we investigated the origin, organization and evolution of two head-to-head arranged gene pairs of homologs of AtOFP2-AtOFP17, and, AtOFP4-AtOFP20 across Archaeplastida. The ancestral forms of AtOFP2, AtOFP4, AtOFP17, and AtOFP20 are likely to have evolved in last common ancestors of Embryophyta (land plants) given their complete absence in Rhodophyta and Chlorophyta. The OFP gene family originated and expanded in Bryophyta, including protein variants with complete (OVATE-OFP) or partial (OVATE-Like OFP) OVATE domain; with head-to-head organization present only in Spermatophyta (gymnosperms and angiosperms). Ancestral State Reconstruction revealed the origin of head-to-head organized gene pair in gymnosperms, with both genes being OVATE-OFP (homologs of AtOFP2/4). Phylogenetic reconstruction and copy number analysis suggests the presence of a single copy of the head-to-head arranged pair of OFP2/4 (OVATE)-OFP17/20 (OVATE-Like) in all angiosperms except Brassicaceae, and a duplication event in last common ancestor of core Brassicaceae approximately 32-54 MYA leading to origin of AtOFP2-AtOFP17 and AtOFP4-AtOFP20 as paralogs. Synteny analysis of genomic regions harbouring homologs of AtOFP2-AtOFP17, AtOFP4-AtOFP20 and AtOFP2/4-AtOFP17/20 across angiosperms suggested ancestral nature of AtOFP2-AtOFP17 gene pair. The present study thus establishes the orthology and evolutionary history of two non-canonically organised gene pairs with variation in their OVATE domain. The non-canonical organisation, atleast in Brassicaceae, has the potential of generating complex transcriptional regulation mediated via a common bi-directional promoter. The study thus lays down a framework to understand evolution of gene and protein structure, transcriptional regulation and function across a phylogenetic lineage through comparative analyses.

Epigenetic control of transposons during plant reproduction: From meiosis to hybrid seeds

The regulation of transposable elements (TEs) requires overlapping epigenetic modifications that must be reinforced every cell division and generation. In plants, this is achieved by multiple pathways including small RNAs, DNA methylation, and repressive histone marks that act together to control TE expression and activity throughout the entire life cycle. However, transient TE activation is observed during reproductive transitions as a result of epigenome reprogramming, thus providing windows of opportunity for TE proliferation and epigenetic novelty. Ultimately, these events may originate complex TE-driven transcriptional networks or cell-to-cell communication strategies via mobile small RNAs. In this review, we discuss recent findings and current understanding of TE regulation during sexual plant reproduction, and its implications for fertility, early seed development, and epigenetic inheritance.

Analysis of Homologous Regions of Small RNAs MIR397 and MIR408 Reveals the Conservation of Microsynteny among Rice Crop-Wild Relatives

MIRNAs are small non-coding RNAs that play important roles in a wide range of biological processes in plant growth and development. MIR397 (involved in drought, low temperature, and nitrogen and copper (Cu) starvation) and MIR408 (differentially expressed in response to environmental stresses such as copper, light, mechanical stress, dehydration, cold, reactive oxygen species, and drought) belong to conserved MIRNA families that either negatively or positively regulate their target genes. In the present study, we identified the homologs of MIR397 and MIR408 in Oryza sativa and its six wild progenitors, three non-Oryza species, and one dicot species. We analyzed the 100 kb segments harboring MIRNA homologs from 11 genomes to obtain a comprehensive view of their community evolution around these loci in the farthest (distant) relatives of rice. Our study showed that mature MIR397 and MIR408 were highly conserved among all Oryza species. Comparative genomics analyses also revealed that the microsynteny of the 100 kb region surrounding MIRNAs was only conserved in Oryza spp.; disrupted in Sorghum, maize, and wheat; and completely lost in Arabidopsis. There were deletions, rearrangements, and translocations within the 100 kb segments in Oryza spp., but the overall microsynteny of the region was maintained. The phylogenetic analyses of the precursor regions of all MIRNAs under study revealed a bimodal clade of common origin. This comparative analysis of miRNA involved in abiotic stress tolerance in plants provides a powerful tool for future Oryza research. Crop wild relatives (CWRs) offer multiple traits with potential to decrease the amount of yield loss owing to biotic and abiotic stresses. Using a comparative genomics approach, the exploration of CWRs as a source of tolerance to these stresses by understanding their evolution can be further used to leverage their yield potential.

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.

Characterization of novel pollen-expressed transcripts reveals their potential roles in pollen heat stress response in Arabidopsis thaliana

Arabidopsis pollen transcriptome analysis revealed new intergenic transcripts of unknown function, many of which are long non-coding RNAs, that may function in pollen-specific processes, including the heat stress response. The male gametophyte is the most heat sensitive of all plant tissues. In recent years, long noncoding RNAs (lncRNAs) have emerged as important components of cellular regulatory networks involved in most biological processes, including response to stress. While examining RNAseq datasets of developing and germinating Arabidopsis thaliana pollen exposed to heat stress (HS), we identified 66 novel and 246 recently annotated intergenic expressed loci (XLOCs) of unknown function, with the majority encoding lncRNAs. Comparison with HS in cauline leaves and other RNAseq experiments indicated that 74% of the 312 XLOCs are pollen-specific, and at least 42% are HS-responsive. Phylogenetic analysis revealed that 96% of the genes evolved recently in Brassicaceae. We found that 50 genes are putative targets of microRNAs and that 30% of the XLOCs contain small open reading frames (ORFs) with homology to protein sequences. Finally, RNAseq of ribosome-protected RNA fragments together with predictions of periodic footprint of the ribosome P-sites indicated that 23 of these ORFs are likely to be translated. Our findings indicate that many of the 312 unknown genes might be functional and play a significant role in pollen biology, including the HS response.

Independent recurrent evolution of MICRORNA genes converging onto similar non-canonical organisation across green plant lineages is driven by local and segmental duplication events in species, family and lineages

The relationship between evolutionary history, organisation and transcriptional regulation of genes are intrinsically linked. These have been well studied in canonically organised protein-coding genes but not of MIRNAs. In the present study, we investigated the non-canonical arrangement of MIRNAs across taxonomic boundaries from algae to angiosperms employing a combination of genome organization, phylogeny and synteny. We retrieved the complete dataset of MIRNA from twenty-five species to identify and classify based on organisational patterns. The median size of cluster was between 2-5 kb and between 1-20 % of all MIRNAs are organized in head-to-head (with bidirectional promoter), head-to-tail (tandem), and overlapping manner. Although majority of the clusters are composed of MIRNA homologs, 25% of all clusters comprises of non-homologous genes with a potential of generating functional and regulatory complexity. A comparison of phylogeny and organizational patterns revealed that multiple independent events, some of which are species-specific, and some ancient, in different lineages, are responsible for non-canonical organization. Detailed investigation of MIR395 family across the plants revealed a complex origin of non-canonical arrangement through ancient and recent, segmental and local duplications; analysis of MIR399 family revealed major expansion occurred prior to monocot-dicot split, with few lineage-specific events. Evolution of "convergent" organization pattern of non-canonical arrangement originating from independent loci through recurrent event highlights our poor understanding of evolutionary process of MIRNA genes. The present investigation thus paves way for comparative functional genomics to understand the role of non-canonical organization on transcriptional regulation and regulatory diversity in MIRNA gene families.

Identification and Analysis of OVATE Family Members from Genome of the Early Land Plants Provide Insights into Evolutionary History of OFP Family and Function

Mosses, liverworts, hornworts and lycophytes represent transition stages between the aquatic to terrestrial/land plants. Several morphological and adaptive novelties driven by genomic components including emergence and expansion of new or existing gene families have played a critical role during and after the transition, and contributed towards successful colonization of terrestrial ecosystems. It is crucial to decipher the evolutionary transitions and natural selection on the gene structure and function to understand the emergence of phenotypic and adaptive diversity. Plants at the "transition zone", between aquatic and terrestrial ecosystem, are also the most vulnerable because of climate change and may contain clues for successful mitigation of the challenges of climate change. Identification and comparative analyses of such genetic elements and gene families are few in mosses, liverworts, hornworts and lycophytes. Ovate family proteins (OFPs) are plant-specific transcriptional repressors and are acknowledged for their roles in important growth and developmental processes in land plants, and information about the functional aspects of OFPs in early land plants is fragmentary. As a first step towards addressing this gap, a comprehensive in silico analysis was carried out utilizing publicly available genome sequences of Marchantia polymorpha (Mp), Physcomitrella patens (Pp), Selaginella moellendorffii (Sm) and Sphagnum fallax (Sf). Our analysis led to the identification of 4 MpOFPs, 19 PpOFPs, 6 SmOFPs and 3 SfOFPs. Cross-genera analysis revealed a drastic change in the structure and physiochemical properties in OFPs suggesting functional diversification and genomic plasticity during the evolutionary course. Knowledge gained from this comparative analysis will form the framework towards deciphering and dissection of their developmental and adaptive role/s in early land plants and could provide insights into evolutionary strategies adapted by land plants.

A segmental duplication in the common ancestor of Brassicaceae is responsible for the origin of the paralogs KCS6-KCS5, which are not shared with other angiosperms

Novel morphological structures allowed adaptation to dry conditions in early land plants. The cuticle, one such novelty, plays diverse roles in tolerance to abiotic and biotic stresses and plant development. Cuticular waxes represent a major constituent of the cuticle and are comprised of an assortment of chemicals that include, among others, very long chain fatty acids (VLCFAs). Members of the β-ketoacyl coenzyme A synthases (KCS) gene family code for enzymes that are essential for fatty acid biosynthesis. The gene KCS6 (CUT1) is known to be a key player in the production of VLCFA precursors essential for the synthesis of cuticular waxes in the model plant Arabidopsis thaliana (Brassicaceae). Despite its functional importance, relatively little is known about the evolutionary history of KCS6 or its paralog KCS5 in Brassicaceae or beyond. This lacuna becomes important when we extrapolate understanding of mechanisms gained from the model plant to its containing clades Brassicaceae, flowering plants, or beyond. The Brassicaceae, with several sequenced genomes and a known history of paleoploidy, mesopolyploidy and neopolyploidy, offer a system in which to study the evolution and diversification of the KCS6-KCS5 paralogy. Our phylogenetic analyses across green plants, combined with comparative genomic, microsynteny and evolutionary rates analyses across nine genomes of Brassicaceae, reveal that (1) the KCS6-KCS5 paralogy arose as the result of a large segmental duplication in the ancestral Brassicaceae, (2) the KCS6-KCS5 lineage is represented by a single copy in other flowering plant lineages, (3) the duplicated segments undergo different degrees of retention and loss, and (4) most of the genes in the KCS6 and KCS5 gene blocks (including KCS6 and KCS5 themselves) are under purifying selection. The last also true for most members of the KCS gene family in Brassicaceae, except for KCS8, KCS9 and KCS17, which are under positive selection and may be undergoing functional evolution, meriting further investigation. Overall, our results clearly establish that the ancestral KCS6/5 gene duplicated in the Brassicaceae lineage. It is possible that any specialized functions of KCS5 found in Brassicaceae are either part of a set of KCS6/5 gene functions in the rest of the flowering plants, or unique to Brassicaceae.

Transposon-derived small RNAs triggered by miR845 mediate genome dosage response in Arabidopsis

Chromosome dosage has substantial effects on reproductive isolation and speciation in both plants and animals, but the underlying mechanisms are largely obscure 1 . Transposable elements in animals can regulate hybridity through maternal small RNA 2 , whereas small RNAs in plants have been postulated to regulate dosage response via neighboring imprinted genes3,4. Here we show that a highly conserved microRNA in plants, miR845, targets the tRNAMet primer-binding site (PBS) of long terminal repeat (LTR) retrotransposons in Arabidopsis pollen, and triggers the accumulation of 21-22-nucleotide (nt) small RNAs in a dose-dependent fashion via RNA polymerase IV. We show that these epigenetically activated small interfering RNAs (easiRNAs) mediate hybridization barriers between diploid seed parents and tetraploid pollen parents (the 'triploid block'), and that natural variation for miR845 may account for 'endosperm balance' allowing the formation of triploid seeds. Targeting of the PBS with small RNA is a common mechanism for transposon control in mammals and plants, and provides a uniquely sensitive means to monitor chromosome dosage and imprinting in the developing seed.

High-Throughput Sequencing of Small RNA Transcriptomes in Maize Kernel Identifies miRNAs Involved in Embryo and Endosperm Development

Maize kernel development is a complex biological process that involves the temporal and spatial expression of many genes and fine gene regulation at a transcriptional and post-transcriptional level, and microRNAs (miRNAs) play vital roles during this process. To gain insight into miRNA-mediated regulation of maize kernel development, a deep-sequencing technique was used to investigate the dynamic expression of miRNAs in the embryo and endosperm at three developmental stages in B73. By miRNA transcriptomic analysis, we characterized 132 known miRNAs and six novel miRNAs in developing maize kernel, among which, 15 and 14 miRNAs were commonly differentially expressed between the embryo and endosperm at 9 days after pollination (DAP), 15 DAP and 20 DAP respectively. Conserved miRNA families such as miR159, miR160, miR166, miR390, miR319, miR528 and miR529 were highly expressed in developing embryos; miR164, miR171, miR393 and miR2118 were highly expressed in developing endosperm. Genes targeted by those highly expressed miRNAs were found to be largely related to a regulation category, including the transcription, macromolecule biosynthetic and metabolic process in the embryo as well as the vitamin biosynthetic and metabolic process in the endosperm. Quantitative reverse transcription-PCR (qRT-PCR) analysis showed that these miRNAs displayed a negative correlation with the levels of their corresponding target genes. Importantly, our findings revealed that members of the miR169 family were highly and dynamically expressed in the developing kernel, which will help to exploit new players functioning in maize kernel development.