A genome-wide transcriptional activity survey of rice transposable element-related genes

Identification and characterization of RuvBL DNA helicase genes for tolerance against abiotic stresses in bread wheat (Triticum aestivum L.) and related species

Recombination UVB (sensitivity) like (RuvBL) helicase genes represent a conserved family of genes, which are known to be involved in providing tolerance against abiotic stresses like heat and drought. We identified nine wheat RuvBL genes, one each on nine different chromosomes, belonging to homoeologous groups 2, 3, and 4. The lengths of genes ranged from 1647 to 2197 bp and exhibited synteny with corresponding genes in related species including Ae. tauschii, Z. mays, O. sativa, H. vulgare, and B. distachyon. The gene sequences were associated with regulatory cis-elements and transposable elements. Two genes, namely TaRuvBL1a-4A and TaRuvBL1a-4B, also carried targets for a widely known miRNA, tae-miR164. Gene ontology revealed that these genes were closely associated with ATP-dependent formation of histone acetyltransferase complex. Analysis of the structure and function of RuvBL proteins revealed that the proteins were localized mainly in the cytoplasm. A representative gene, namely TaRuvBL1a-4A, was also shown to be involved in protein-protein interactions with ten other proteins. On the basis of phylogeny, RuvBL proteins were placed in two sub-divisions, namely RuvBL1 and RuvBL2, which were further classified into clusters and sub-clusters. In silico studies suggested that these genes were differentially expressed under heat/drought. The qRT-PCR analysis confirmed that expression of TaRuvBL genes differed among wheat cultivars, which differed in the level of thermotolerance. The present study advances our understanding of the biological role of wheat RuvBL genes and should help in planning future studies on RuvBL genes in wheat including use of RuvBL genes in breeding thermotolerant wheat cultivars.

LTR Retroelements and Bird Adaptation to Arid Environments

TEs are known to be among the main drivers in genome evolution, leading to the generation of evolutionary advantages that favor the success of organisms. The aim of this work was to investigate the TE landscape in bird genomes to look for a possible relationship between the amount of specific TE types and environmental changes that characterized the Oligocene era in Australia. Therefore, the mobilome of 29 bird species, belonging to a total of 11 orders, was analyzed. Our results confirmed that LINE retroelements are not predominant in all species of this evolutionary lineage and highlighted an LTR retroelement dominance in species with an Australian-related evolutionary history. The bird LTR retroelement expansion might have happened in response to the Earth’s dramatic climate changes that occurred about 30 Mya, followed by a progressive aridification across most of Australian landmasses. Therefore, in birds, LTR retroelement burst might have represented an evolutionary advantage in the adaptation to arid/drought environments.

Genome-wide identification and characterization of exapted transposable elements in the large genome of sunflower (Helianthus annuus L.)

Transposable elements (TEs) are an important source of genome variability, playing many roles in the evolution of eukaryotic species. Besides well-known phenomena, TEs may undergo the exaptation process and generate the so-called exapted transposable element genes (ETEs). Here we present a genome-wide survey of ETEs in the large genome of sunflower (Helianthus annuus L.), in which the massive amount of TEs, provides a significant source for exaptation. A library of sunflower TEs was used to build TE-specific Hidden Markov Model profiles, to search for all available sunflower gene products. In doing so, 20 016 putative ETEs were identified and further investigated for the characteristics that distinguish TEs from genes, leading to the validation of 3530 ETEs. The analysis of ETEs transcription patterns under different stress conditions showed a differential regulation triggered by treatments mimicking biotic and abiotic stress; furthermore, the distribution of functional domains of differentially regulated ETEs revealed a relevant presence of domains involved in many aspects of cellular functions. A comparative genomic investigation was performed including species representative of Asterids and appropriate outgroups: the bulk of ETEs that resulted were specific to the sunflower, while few ETEs presented orthologues in the genome of all analyzed species, making the hypothesis of a conserved function. This study highlights the crucial role played by exaptation, actively contributing to species evolution.

MKKK20 works as an upstream triple-kinase of MKK3-MPK6-MYC2 module in Arabidopsis seedling development

The mitogen-activated protein kinase (MAPK) cascade is involved in several signal transduction processes in eukaryotes. Here, we report a mechanistic function of MAP kinase kinase kinase 20 (MKKK20) in light signal transduction pathways. We show that MKKK20 acts as a negative regulator of photomorphogenic growth at various wavelengths of light. MKKK20 not only regulates the expression of light signaling pathway regulatory genes but also gets regulated by the same pathway genes. The atmyc2 mkkk20 double mutant analysis shows that MYC2 works downstream to MKKK20 in the regulation of photomorphogenic growth. MYC2 directly binds to the promoter of MKKK20 to modulate its expression. The protein-protein interaction study indicates that MKKK20 physically interacts with MYC2, and this interaction likely suppresses the MYC2-mediated promotion of MKKK20 expression. Further, the protein phosphorylation studies demonstrate that MKKK20 works as the upstream kinase of MKK3-MPK6-MYC2 module in photomorphogenesis.

Environmental stress and transposons in plants

Transposons were once thought to be junk repetitive DNA in the genome. However, their importance gradually became apparent as it became clear that they regulate gene expression, which is essential for organisms to survive, and that they are important factors in the driving force of evolution. Since there are multiple transposons in the genomes of all organisms, transposons have likely been activated and increased in copy number throughout their long history. This review focuses on environmental stress as a factor in transposon activation, paying particular attention to transposons in plants that are activated by environmental stresses. It is now known that plants respond to environmental stress in various ways, and correspondingly, many transposons respond to stress. The relationship between environmental stress and transposons is reviewed, including the mechanisms of their activation and the effects of transposon activation on host plants.

Differential nuclease sensitivity profiling uncovers a drought responsive change in maize leaf chromatin structure for two large retrotransposon derivatives, Uloh and Vegu

Plant chromatin dynamics are generally recognized as playing a role in the genomic response to environmental stress. Although stress-induced transcriptional activities of LTR-retrotransposons have been reported, whether the stress response can be detected at the level of chromatin structure for LTR-retrotransposons is largely unknown. Using differential nuclease sensitivity profiling, we identified that two out of 29 maize LTR-retrotransposon families change their chromatin structure in response to drought stress in leaf tissue. The two LTR-retrotransposon families, uloh and vegu, are classified as nonautonomous LTR-retrotransposons. Differently from other families, the chromatin structure of these two families shifted from more open in normal conditions to more closed following drought stress. Although uloh and vegu lack sequence similarity, most of them shared an intriguing feature of having a new and uncharacterized insertion of a DNA sequence near one side of an LTR. In the uloh family, nine members with a strong drought response also exhibited a drought-induced reduction of published H3K4me3 histone modification in the inserted DNA region, implicating this modification in the chromatin structural changes. Our results provide new insight into how LTR-retrotransposons can alter their chromatin structure following stress response in plants.

Genome-wide discovery of natural variation in pre-mRNA splicing and prioritising causal alternative splicing to salt stress response in rice

Pre-mRNA splicing is an essential step for the regulation of gene expression. In order to specifically capture splicing variants in plants for genome-wide association studies (GWAS), we developed a software tool to quantify and visualise Variations of Splicing in Population (VaSP). VaSP can quantify splicing variants from short-read RNA-seq datasets and discover genotype-specific splicing (GSS) events, which can be used to prioritise causal pre-mRNA splicing events in GWAS. We applied our method to an RNA-seq dataset with 328 samples from 82 genotypes from a rice diversity panel exposed to optimal and saline growing conditions. In total, 764 significant GSS events were identified in salt stress conditions. GSS events were used as markers for a GWAS with the shoot Na⁺ accumulation, which identified six GSS events in five genes significantly associated with the shoot Na⁺ content. Two of these genes, OsNUC1 and OsRAD23 emerged as top candidate genes with splice variants that exhibited significant divergence between the variants for shoot growth under salt stress conditions. VaSP is a versatile tool for alternative splicing analysis in plants and a powerful tool for prioritising candidate causal pre-mRNA splicing and corresponding genomic variations in GWAS.

Transcriptional regulatory network of the light signaling pathways

The developmental program by which plants respond is tightly controlled by a complex cascade in which photoreceptors perceive and transduce the light signals that drive signaling processes and direct the transcriptional reprogramming, yielding specific cellular responses. The molecular mechanisms involved in the transcriptional regulation include light-regulated nuclear localization (the phytochromes and UVR8) and nuclear accumulation (the cryptochrome, cry2) of photoreceptors. This regulatory cascade also includes master regulatory transcription factors (TFs) that bridge photoreceptor activation with chromatin remodeling and regulate the expression of numerous light-responsive genes. Light signaling-related TFs often function as signal convergence points in concert with TFs in other signaling pathways to integrate complex endogenous and environmental cues that help the plant adapt to the surrounding environment. Increasing evidence suggests that chromatin modifications play a critical role in regulating light-responsive gene expression and provide an additional layer of light signaling regulation. Here, we provide an overview of our current knowledge of the transcriptional regulatory network involved in the light response, particularly the roles of TFs and chromatin in regulating light-responsive gene expression.

Transcriptional activity of PIF and Pong-like Class II transposable elements in Triticeae

Background

Transposable elements are major contributors to genome size and variability, accounting for approximately 70–80% of the maize, barley, and wheat genomes. PIF and Pong-like elements belong to two closely-related element families within the PIF/Harbinger superfamily of Class II (DNA) transposons. Both elements contain two open reading frames; one encodes a transposase (ORF2) that catalyzes transposition of the functional elements and their related non-autonomous elements, while the function of the second is still debated. In this work, we surveyed for PIF- and Pong-related transcriptional activity in 13 diploid Triticeae species, all of which have been previously shown to harbor extensive within-genome diversity of both groups of elements.

Results

The results revealed that PIF elements have considerable transcriptional activity in Triticeae, suggesting that they can escape the initial levels of plant cell control and are regulated at the post-transcriptional level. Phylogenetic analysis of 156 PIF cDNA transposase fragments along with 240 genomic partial transposase sequences showed that most, if not all, PIF clades are transcriptionally competent, and that multiple transposases coexisting within a single genome have the potential to act simultaneously. In contrast, we did not detect any transcriptional activity of Pong elements in any sample.

Conclusions

The lack of Pong element transcription shows that even closely related transposon families can exhibit wide variation in their transposase transcriptional activity within the same genome.

Electronic supplementary material

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

Transgenerational epimutations induced by multi-generation drought imposition mediate rice plant's adaptation to drought condition

Epigenetic mechanisms are crucial mediators of appropriate plant reactions to adverse environments, but their involvement in long-term adaptation is less clear. Here, we established two rice epimutation accumulation lines by applying drought conditions to 11 successive generations of two rice varieties. We took advantage of recent technical advances to examine the role of DNA methylation variations on rice adaptation to drought stress. We found that multi-generational drought improved the drought adaptability of offspring in upland fields. At single-base resolution, we discovered non-random appearance of drought-induced epimutations. Moreover, we found that a high proportion of drought-induced epimutations maintained their altered DNA methylation status in advanced generations. In addition, genes related to transgenerational epimutations directly participated in stress-responsive pathways. Analysis based on a cluster of drought-responsive genes revealed that their DNA methylation patterns were affected by multi-generational drought. These results suggested that epigenetic mechanisms play important roles in rice adaptations to upland growth conditions. Epigenetic variations have morphological, physiological and ecological consequences and are heritable across generations, suggesting that epigenetics can be considered an important regulatory mechanism in plant long-term adaptation and evolution under adverse environments.

Moving through the Stressed Genome: Emerging Regulatory Roles for Transposons in Plant Stress Response

The recognition of a positive correlation between organism genome size with its transposable element (TE) content, represents a key discovery of the field of genome biology. Considerable evidence accumulated since then suggests the involvement of TEs in genome structure, evolution and function. The global genome reorganization brought about by transposon activity might play an adaptive/regulatory role in the host response to environmental challenges, reminiscent of McClintock's original 'Controlling Element' hypothesis. This regulatory aspect of TEs is also garnering support in light of the recent evidences, which project TEs as "distributed genomic control modules." According to this view, TEs are capable of actively reprogramming host genes circuits and ultimately fine-tuning the host response to specific environmental stimuli. Moreover, the stress-induced changes in epigenetic status of TE activity may allow TEs to propagate their stress responsive elements to host genes; the resulting genome fluidity can permit phenotypic plasticity and adaptation to stress. Given their predominating presence in the plant genomes, nested organization in the genic regions and potential regulatory role in stress response, TEs hold unexplored potential for crop improvement programs. This review intends to present the current information about the roles played by TEs in plant genome organization, evolution, and function and highlight the regulatory mechanisms in plant stress responses. We will also briefly discuss the connection between TE activity, host epigenetic response and phenotypic plasticity as a critical link for traversing the translational bridge from a purely basic study of TEs, to the applied field of stress adaptation and crop improvement.

DNA methylation changes facilitated evolution of genes derived from Mutator-like transposable elements

Background

Mutator-like transposable elements, a class of DNA transposons, exist pervasively in both prokaryotic and eukaryotic genomes, with more than 10,000 copies identified in the rice genome. These elements can capture ectopic genomic sequences that lead to the formation of new gene structures. Here, based on whole-genome comparative analyses, we comprehensively investigated processes and mechanisms of the evolution of putative genes derived from Mutator-like transposable elements in ten Oryza species and the outgroup Leersia perieri, bridging ~20 million years of evolutionary history.

Results

Our analysis identified thousands of putative genes in each of the Oryza species, a large proportion of which have evidence of expression and contain chimeric structures. Consistent with previous reports, we observe that the putative Mutator-like transposable element-derived genes are generally GC-rich and mainly derive from GC-rich parental sequences. Furthermore, we determine that Mutator-like transposable elements capture parental sequences preferentially from genomic regions with low methylation levels and high recombination rates. We explicitly show that methylation levels in the internal and terminated inverted repeat regions of these elements, which might be directed by the 24-nucleotide small RNA-mediated pathway, are different and change dynamically over evolutionary time. Lastly, we demonstrate that putative genes derived from Mutator-like transposable elements tend to be expressed in mature pollen, which have undergone de-methylation programming, thereby providing a permissive expression environment for newly formed/transposable element-derived genes.

Conclusions

Our results suggest that DNA methylation may be a primary mechanism to facilitate the origination, survival, and regulation of genes derived from Mutator-like transposable elements, thus contributing to the evolution of gene innovation and novelty in plant genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-0954-8) contains supplementary material, which is available to authorized users.

Mutator and MULE Transposons

The Mutator system of transposable elements (TEs) is a highly mutagenic family of transposons in maize. Because they transpose at high rates and target genic regions, these transposons can rapidly generate large numbers of new mutants, which has made the Mutator system a favored tool for both forward and reverse mutagenesis in maize. Low copy number versions of this system have also proved to be excellent models for understanding the regulation and behavior of Class II transposons in plants. Notably, the availability of a naturally occurring locus that can heritably silence autonomous Mutator elements has provided insights into the means by which otherwise active transposons are recognized and silenced. This chapter will provide a review of the biology, regulation, evolution and uses of this remarkable transposon system, with an emphasis on recent developments in our understanding of the ways in which this TE system is recognized and epigenetically silenced as well as recent evidence that Mu-like elements (MULEs) have had a significant impact on the evolution of plant genomes.

Structural Diversity of a Novel LTR Retrotransposon, RTPOSON, in the Genus Oryza

Retrotransposons with long terminal repeats (LTRs) are the most abundant transposable elements in plant genomes. A novel LTR retrotransposon named RTPOSON primarily occurs in the genus Oryza and in several species of the Poaceae family. RTPOSON has been identified in the Ty1-copia group of retrotransposons because two of its open reading frames encode an uncharacterized protein and UBN2_2 and zinc knuckle, respectively. More than 700 RTPOSONs were identified in Oryza genomes; 127 RTPOSONs with LTRs and gag-pol elements were classified into three subgroups. The subgroup RTPOSON_sub3 had the smallest DNA size and 97% (32/33) of RTPOSON elements from Oryza punctata are classified in this group. The insertion time of these RTPOSONs varied and their proliferation occurred within the last 8 Mya, with two bursting periods within the last 1.5–5.0 Mya. A total of 37 different orthologous insertions of RTPOSONs, with different nested transposable elements and gene fragments, were identified by comparing the genomes of ssp. japonica cv. Nipponbare and ssp. indica cv. 93–11. A part of intact RTPOSON elements was evolved independently after the divergence of indica and japonica. In addition, intact RTPOSONs and homologous fragments were preferentially retained or integrated in genic regions. This novel LTR retrotransposon, RTPOSON, might have an impact on genome evolution, genic innovation, and genetic variation.

A highly specific microRNA-mediated mechanism silences LTR retrotransposons of strawberry

Small RNAs are involved in a plethora of functions in plant genomes. In general, transcriptional gene silencing is mediated by 24-nucleotide siRNAs and is required for maintaining transposable elements in a silenced state. However, microRNAs are not commonly associated with transposon silencing. In this study, we performed small RNA transcriptome and degradome analyses of the Rosaceae model plant Fragaria vesca (the woodland strawberry) at the genome-wide level, and identified miRNA families and their targets. We report a highly specific mechanism of LTR retrotransposon silencing mediated by an abundant, ubiquitously expressed miRNA (fve-miR1511) generated from a single locus. This miRNA specifically targets LTR retroelements, silencing them post-transcriptionally by perfectly pairing to the highly conserved primer binding site for methionyl initiator tRNA that is essential for reverse transcription. We investigated the possible origins of this miRNA, and present evidence that the pre-miR1511 hairpin structure probably derived from a locus coding for tRNA(iM) (et) through a single microinversion event. Our study shows that this miRNA targets retrotransposons specifically and constitutively, and contributes to features such as genome stability, size and architecture in a far more direct way than previously thought.

What makes up plant genomes: The vanishing line between transposable elements and genes

The ultimate source of evolution is mutation. As the largest component in plant genomes, transposable elements (TEs) create numerous types of mutations that cannot be mimicked by other genetic mechanisms. When TEs insert into genomic sequences, they influence the expression of nearby genes as well as genes unlinked to the insertion. TEs can duplicate, mobilize, and recombine normal genes or gene fragments, with the potential to generate new genes or modify the structure of existing genes. TEs also donate their transposase coding regions for cellular functions in a process called TE domestication. Despite the host defense against TE activity, a subset of TEs survived and thrived through discreet selection of transposition activity, target site, element size, and the internal sequence. Finally, TEs have established strategies to reduce the efficacy of host defense system by increasing the cost of silencing TEs. This review discusses the recent progress in the area of plant TEs with a focus on the interaction between TEs and genes.

Discovery of novel genes derived from transposable elements using integrative genomic analysis

Complex eukaryotes contain millions of transposable elements (TEs), comprising large fractions of their nuclear genomes. TEs consist of structural, regulatory, and coding sequences that are ordinarily associated with transposition, but that occasionally confer on the organism a selective advantage and may thereby become exapted. Exapted transposable element genes (ETEs) are known to play critical roles in diverse systems, from vertebrate adaptive immunity to plant development. Yet despite their evident importance, most ETEs have been identified fortuitously and few systematic searches have been conducted, suggesting that additional ETEs may await discovery. To explore this possibility, we develop a comprehensive systematic approach to searching for ETEs. We use TE-specific conserved domains to identify with high precision genes derived from TEs and screen them for signatures of exaptation based on their similarities to reference sets of known ETEs, conventional (non-TE) genes, and TE genes across diverse genetic attributes including repetitiveness, conservation of genomic location and sequence, and levels of expression and repressive small RNAs. Applying this approach in the model plant Arabidopsis thaliana, we discover a surprisingly large number of novel high confidence ETEs. Intriguingly, unlike known plant ETEs, several of the novel ETE families form tandemly arrayed gene clusters, whereas others are relatively young. Our results not only identify novel TE-derived genes that may have practical applications but also challenge the notion that TE exaptation is merely a relic of ancient life, instead suggesting that it may continue to fundamentally drive evolution.

Transcriptomic dissection of the rice-Burkholderia glumae interaction

BACKGROUND

Bacterial panicle blight caused by the bacterium Burkholderia glumae is an emerging disease of rice in the United States. Not much is known about this disease, the disease cycle or any source of disease resistance. To understand the interaction between rice and Burkholderia glumae, we used transcriptomics via next-generation sequencing (RNA-Seq) and bioinformatics to identify differentially expressed transcripts between resistant and susceptible interactions and formulate a model for rice resistance to the disease.

RESULTS

Using inoculated young seedlings as sample tissues, we identified unique transcripts involved with resistance to bacterial panicle blight, including a PIF-like ORF1 and verified differential expression of some selected genes using qRT-PCR. These transcripts, which include resistance genes of the NBS-LRR type, kinases, transcription factors, transporters and expressed proteins with functions that are not known, have not been reported in other pathosystems including rice blast or bacterial blight. Further, functional annotation analysis reveals enrichment of defense response and programmed cell death (biological processes); ATP and protein binding (molecular functions); and mitochondrion-related (cell component) transcripts in the resistant interaction.

CONCLUSION

Taken together, we formulated a model for rice resistance to bacterial panicle blight that involves an activation of previously unknown resistance genes and their activation partners upon challenge with B. glumae. Other interesting findings are that 1) though these resistance transcripts were up-regulated upon inoculation in the resistant interaction, some of them were already expressed in the water-inoculated control from the resistant genotype, but not in the water- and bacterium-inoculated samples from the susceptible genotype; 2) rice may have co-opted an ORF that was previously a part of a transposable element to aid in the resistance mechanism; and 3) resistance may have existed immediately prior to rice domestication.

Low temperature-responsive changes in the anther transcriptome's repeat sequences are indicative of stress sensitivity and pollen sterility in rice strains

Genome-wide transcriptome analyses using microarray probes containing genes and repeat sequences have been performed to examine responses to low temperatures in rice (Oryza sativa). We focused particularly on the rice anther at the booting stage, because a low temperature at this stage can result in pollen abortion. The five rice strains examined in this study showed different pollen fertilities due to a low-temperature treatment during the booting stage. The microarray analyses demonstrated that the low-temperature stress caused genome-wide changes in the transcriptional activities not only of genes but also of repeat sequences in the rice anther. The degree of the temperature-responsive changes varied among the five rice strains. Interestingly, the low-temperature-sensitive strains revealed more changes in the transcriptome when compared with the tolerant strains. The expression patterns of the repeat sequences, including miniature inverted-repeat transposable elements, transposons, and retrotransposons, were correlated with the pollen fertilities of the five strains, with the highest correlation coefficient being 0.979. Even in the low-temperature-sensitive strains, the transcriptomes displayed distinct expression patterns. The elements responding to the low temperatures were evenly distributed throughout the genome, and the major cis-motifs involved in temperature-responsive changes were undetectable from the upstream sequences in the corresponding repeats. The genome-wide responses of transcription to the temperature shift may be associated with chromatin dynamics, which facilitates environmental plasticity. A genome-wide analysis using repeat sequences suggested that stress tolerance could be conferred by insensitivity to the stimuli.

Genome-wide identification and analysis of Japonica and Indica cultivar-preferred transcripts in rice using 983 Affymetrix array data

BACKGROUND

Accumulation of genome-wide transcriptome data provides new insight on a genomic scale which cannot be gained by analyses of individual data. The majority of rice (O. sativa) species are japonica and indica cultivars. Genome-wide identification of genes differentially expressed between japonica and indica cultivars will be very useful in understanding the domestication and evolution of rice species.

RESULTS

In this study, we analyzed 983 of the 1866 entries in the Affymetrix array data in the public database: 595 generated from indica and 388 from japonica rice cultivars. To discover differentially expressed genes in each cultivar, we performed significance analysis of microarrays for normalized data, and identified 490 genes preferentially expressed in japonica and 104 genes in indica. Gene Ontology analyses revealed that defense response-related genes are significantly enriched in both cultivars, indicating that japonica and indica might be under strong selection pressure for these traits during domestication. In addition, 36 (34.6%) of 104 genes preferentially expressed in indica and 256 (52.2%) of 490 genes preferentially expressed in japonica were annotated as genes of unknown function. Biotic stress overview in the MapMan toolkit revealed key elements of the signaling pathway for defense response in japonica or indica eQTLs.

CONCLUSIONS

The percentage of screened genes preferentially expressed in indica was 4-fold higher (34.6%) and that in japonica was 5-fold (52.2%) higher than expected (11.1%), suggesting that genes of unknown function are responsible for the novel traits that distinguish japonica and indica cultivars. The identification of 10 functionally characterized genes expressed preferentially in either japonica or indica highlights the significance of our candidate genes during the domestication of rice species. Functional analysis of the roles of individual components of stress-mediated signaling pathways will shed light on potential molecular mechanisms to improve disease resistance in rice.

Gene body methylation shows distinct patterns associated with different gene origins and duplication modes and has a heterogeneous relationship with gene expression in Oryza sativa (rice)

Whole-genome duplication (WGD) has been recurring and single-gene duplication is also widespread in angiosperms. Recent whole-genome DNA methylation maps indicate that gene body methylation (i.e. of coding regions) has a functional role. However, whether gene body methylation is related to gene origins and duplication modes has yet to be reported. In rice (Oryza sativa), we computed a body methylation level (proportion of methylated CpG within coding regions) for each gene in five tissues. Body methylation levels follow a bimodal distribution, but show distinct patterns associated with transposable element-related genes; WGD, tandem, proximal and transposed duplicates; and singleton genes. For pairs of duplicated genes, divergence in body methylation levels increases with physical distance and synonymous (Ks) substitution rates, and WGDs show lower divergence than single-gene duplications of similar Ks levels. Intermediate body methylation tends to be associated with high levels of gene expression, whereas heavy body methylation is associated with lower levels of gene expression. The biological trends revealed here are consistent across five rice tissues, indicating that genes of different origins and duplication modes have distinct body methylation patterns, and body methylation has a heterogeneous relationship with gene expression and may be related to survivorship of duplicated genes.

The SLEEPER genes: a transposase-derived angiosperm-specific gene family

BACKGROUND

DAYSLEEPER encodes a domesticated transposase from the hAT-superfamily, which is essential for development in Arabidopsis thaliana. Little is known about the presence of DAYSLEEPER orthologs in other species, or how and when it was domesticated. We studied the presence of DAYSLEEPER orthologs in plants and propose a model for the domestication of the ancestral DAYSLEEPER gene in angiosperms.

RESULTS

Using specific BLAST searches in genomic and EST libraries, we found that DAYSLEEPER-like genes (hereafter called SLEEPER genes) are unique to angiosperms. Basal angiosperms as well as grasses (Poaceae) and dicotyledonous plants possess such putative orthologous genes, but SLEEPER-family genes were not found in gymnosperms, mosses and algae. Most species contain more than one SLEEPER gene. All SLEEPERs contain a C2H2 type BED-zinc finger domain and a hATC dimerization domain. We designated 3 motifs, partly overlapping the BED-zinc finger and dimerization domain, which are hallmark features in the SLEEPER family. Although SLEEPER genes are structurally conserved between species, constructs with SLEEPER genes from grapevine and rice did not complement the daysleeper phenotype in Arabidopsis, when expressed under control of the DAYSLEEPER promoter. However these constructs did cause a dominant phenotype when expressed in Arabidopsis. Rice plant lines with an insertion in the RICESLEEPER1 or 2 locus displayed phenotypic abnormalities, indicating that these genes are functional and important for normal development in rice. We suggest a model in which we hypothesize that an ancestral hAT transposase was retrocopied and stably integrated in the genome during early angiosperm evolution. Evidence is also presented for more recent retroposition events of SLEEPER genes, such as an event in the rice genome, which gave rise to the RICESLEEPER1 and 2 genes.

CONCLUSIONS

We propose the ancestral SLEEPER gene was formed after a process of retro-transposition during the evolution of the first angiosperms. It may have acquired an important function early on, as mutation of two SLEEPER genes in rice, like the daysleeper mutant in A. thaliana gave a developmental phenotype indicative of their importance for normal plant development.

A gene family derived from transposable elements during early angiosperm evolution has reproductive fitness benefits in Arabidopsis thaliana

The benefits of ever-growing numbers of sequenced eukaryotic genomes will not be fully realized until we learn to decipher vast stretches of noncoding DNA, largely composed of transposable elements. Transposable elements persist through self-replication, but some genes once encoded by transposable elements have, through a process called molecular domestication, evolved new functions that increase fitness. Although they have conferred numerous adaptations, the number of such domesticated transposable element genes remains unknown, so their evolutionary and functional impact cannot be fully assessed. Systematic searches that exploit genomic signatures of natural selection have been employed to identify potential domesticated genes, but their predictions have yet to be experimentally verified. To this end, we investigated a family of domesticated genes called MUSTANG (MUG), identified in a previous bioinformatic search of plant genomes. We show that MUG genes are functional. Mutants of Arabidopsis thaliana MUG genes yield phenotypes with severely reduced plant fitness through decreased plant size, delayed flowering, abnormal development of floral organs, and markedly reduced fertility. MUG genes are present in all flowering plants, but not in any non-flowering plant lineages, such as gymnosperms, suggesting that the molecular domestication of MUG may have been an integral part of early angiosperm evolution. This study shows that systematic searches can be successful at identifying functional genetic elements in noncoding regions and demonstrates how to combine systematic searches with reverse genetics in a fruitful way to decipher eukaryotic genomes.

The genomes of complex organisms are mostly made up not of ordinary genes but of transposable elements. Transposable elements have been called “selfish DNA” because they normally persist by copying themselves, not by helping the organism to survive or reproduce. Yet transposable elements can help organisms to evolve; for instance, transposable element genes sometimes acquire new functions that do benefit the organism. Because they are difficult to distinguish from transposable elements, little is known about these “domesticated genes.” Although studies have attempted to identify them computationally, the predictions have not been verified experimentally. Here, we examine some of the first domesticated genes to be predicted computationally, the MUSTANG family of plant genes. We show that the predictions were correct: MUSTANGs are, like ordinary genes, functional. MUSTANG mutations result in serious defects in how plants grow, flower, and reproduce. Since they are present only in flowering plants, MUSTANG probably originated when flowers first evolved, perhaps taking on a key role. This study is important both because it shows that MUSTANG is critical to plant fitness and because, in the future, a similar approach can be used to find additional domesticated genes and to better understand how transposable elements contribute to evolution.

Genotype-dependent Burst of Transposable Element Expression in Crowns of Hexaploid Wheat (Triticum aestivum L.) during Cold Acclimation

The expression of 1,613 transposable elements (TEs) represented in the Affymetrix Wheat Genome Chip was examined during cold treatment in crowns of four hexaploid wheat genotypes that vary in tolerance to cold and in flowering time. The TE expression profiles showed a constant level of expression throughout the experiment in three of the genotypes. In winter Norstar, the most cold-hardy of the four genotypes, a subset of the TEs showed a burst of expression after vernalization saturation was achieved. About 47% of the TEs were expressed, and both Class I (retrotransposons) and Class II (DNA transposons) types were well represented. Gypsy and Copia were the most represented among the retrotransposons while CACTA and Mariner were the most represented DNA transposons. The data suggests that the Vrn-A1 region plays a role in the stage-specific induction of TE expression in this genotype.

Diversification of hAT transposase paralogues in the sugarcane genome

Transposons are abundant components of eukaryotic genomes, and play important role in genome evolution. The knowledge about these elements should contribute to the understanding of their impact on the host genomes. The hAT transposon superfamily is one of the best characterized superfamilies in diverse organisms, nevertheless, a detailed study of these elements was never carried in sugarcane. To address this question we analyzed 32 cDNAs similar to that of hAT superfamily of transposons previously identified in the sugarcane transcriptome. Our results revealed that these hAT-like transposases cluster in one highly homogeneous and other more heterogeneous lineage. We present evidences that support the hypothesis that the highly homogeneous group is a domesticated transposase while the remainder of the lineages are composed of transposon units. The first is common to grasses, clusters significantly with domesticated transposases from Arabidopsis, rice and sorghum and is expressed in different tissues of two sugarcane cultivars analyzed. In contrast, the more heterogeneous group represents at least two transposon lineages. We recovered five genomic versions of one lineage, characterizing a novel transposon family with conserved DDE motif, named SChAT. These results indicate the presence of at least three distinct lineages of hAT-like transposase paralogues in sugarcane genome, including a novel transposon family described in Saccharum and a domesticated transposase. Taken together, these findings permit to follow the diversification of some hAT transposase paralogues in sugarcane, aggregating knowledge about the co-evolution of transposons and their host genomes.

Transposable elements in Phyllostachys pubescens (Poaceae) genome survey sequences and the full-length cDNA sequences, and their association with simple-sequence repeats

Phyllostachys pubescens is a woody bamboo with the highest ecological, economic, and cultural values of all bamboos in Asia. There is more genomic data available for P. pubescens than for any other bamboo species, including 2.12-Mb genome survey sequences (GSS) and 11.4-Mb full-length cDNA sequences (FL-cDNAs) currently deposited in GenBank. Analysis of these sequences revealed that transposable elements (TEs) are abundant, diverse and polyphyletic in the P. pubescens genome, of which Ty3-gypsy and Ty1-copia are the two most abundant families. Phylogenic analysis showed that both elements probably arose before the Bambusoideae separated from the other Poaceae subfamilies. We found evidence that the distribution of some intragenic TEs correlated with transcript profiles, of which Mutator elements preferred to insert in the transcripts of transcription factors. Additionally, we found that the abundance of SSRs in TEs (4.56%) was significantly higher than in GSS (0.098%) and in FL-cDNAs (2.60%) in P. pubescens genome, and TA/AT and CT/AG repeats were found to be intimately associated with En/Spm and Mutator elements, respectively. Our data provide a glimpse of the structure and evolution of P. pubescens genome, although large-scale sequencing of the genome would be required to fully understand the architecture of the P. pubescens genome.

Retrogenes in rice (Oryza sativa L. ssp. japonica) exhibit correlated expression with their source genes

Gene duplication occurs by either DNA- or RNA-based processes; the latter duplicates single genes via retroposition of messenger RNA. The expression of a retroposed gene copy (retrocopy) is expected to be uncorrelated with its source gene because upstream promoter regions are usually not part of the retroposition process. In contrast, DNA-based duplication often encompasses both the coding and the intergenic (promoter) regions; hence, expression is often correlated, at least initially, between DNA-based duplicates. In this study, we identified 150 retrocopies in rice (Oryza sativa L. ssp japonica), most of which represent ancient retroposition events. We measured their expression from high-throughput RNA sequencing (RNAseq) data generated from seven tissues. At least 66% of the retrocopies were expressed but at lower levels than their source genes. However, the tissue specificity of retrogenes was similar to their source genes, and expression between retrocopies and source genes was correlated across tissues. The level of correlation was similar between RNA- and DNA-based duplicates, and they decreased over time at statistically indistinguishable rates. We extended these observations to previously identified retrocopies in Arabidopsis thaliana, suggesting they may be general features of the process of retention of plant retrogenes.

Cell-type specific analysis of translating RNAs in developing flowers reveals new levels of control

Combining translating ribosome affinity purification with RNA-seq for cell-specific profiling of translating RNAs in developing flowers.

Cell type comparisons of cell type-specific hormone responses, promoter motifs, coexpressed cognate binding factor candidates, and splicing isoforms.

Widespread post-transcriptional regulation at both the intron splicing and translational stages.

A new class of noncoding RNAs associated with polysomes.

What constitutes a differentiated cell type? How much do cell types differ in their transcription of genes? The development and functions of tissues rely on constant interactions among distinct and nonequivalent cell types. Answering these questions will require quantitative information on transcriptomes, proteomes, protein–protein interactions, protein–nucleic acid interactions, and metabolomes at cellular resolution. The systems approaches emerging in biology promise to explain properties of biological systems based on genome-wide measurements of expression, interaction, regulation, and metabolism. To facilitate a systems approach, it is essential first to capture such components in a global manner, ideally at cellular resolution.

Recently, microarray analysis of transcriptomes has been extended to a cellular level of resolution by using laser microdissection or fluorescence-activated sorting (for review, see Nelson et al, 2008). These methods have been limited by stresses associated with cellular separation and isolation procedures, and biases associated with mandatory RNA amplification steps. A newly developed method, translating ribosome affinity purification (TRAP; Zanetti et al, 2005; Heiman et al, 2008; Mustroph et al, 2009), circumvents these problems by epitopetagging a ribosomal protein in specific cellular domains to selectively purify polysomes. We combined TRAP with deep sequencing, which we term TRAP-seq, to provide cell-level spatiotemporal maps for Arabidopsis early floral development at single-base resolution.

Flower development in Arabidopsis has been studied extensively and is one of the best understood aspects of plant development (for review, see Krizek and Fletcher, 2005). Genetic analysis of homeotic mutants established the ABC model, in which three classes of regulatory genes, A, B and C, work in a combinatorial manner to confer organ identities of four whorls (Coen and Meyerowitz, 1991). Each class of regulatory gene is expressed in a specific and evolutionarily conserved domain, and the action of the class A, B and C genes is necessary for specification of organ identity (Figure 1A).

Using TRAP-seq, we purified cell-specific translating mRNA populations, which we and others call the translatome, from the A, B and C domains of early developing flowers, in which floral patterning and the specification of floral organs is established. To achieve temporal specificity, we used a floral induction system to facilitate collection of early stage flowers (Wellmer et al, 2006). The combination of TRAP-seq with domain-specific promoters and this floral induction system enabled fine spatiotemporal isolation of translating mRNA in specific cellular domains, and at specific developmental stages.

Multiple lines of evidence confirmed the specificity of this approach, including detecting the expression in expected domains but not in other domains for well-studied flower marker genes and known physiological functions (Figures 1B–D and 2A–C). Furthermore, we provide numerous examples from flower development in which a spatiotemporal map of rigorously comparable cell-specific translatomes makes possible new views of the properties of cell domains not evident in data obtained from whole organs or tissues, including patterns of transcription and cis-regulation, new physiological differences among cell domains and between flower stages, putative hormone-active centers, and splicing events specific for flower domains (Figure 2A–D). Such findings may provide new targets for reverse genetics studies and may aid in the formulation and validation of interaction and pathway networks.

Beside cellular heterogeneity, the transcriptome is regulated at several steps through the life of mRNA molecules, which are not directly available through traditional transcriptome profiling of total mRNA abundance. By comparing the translatome and transcriptome, we integratively profiled two key posttranscriptional control points, intron splicing and translation state. From our translatome-wide profiling, we (i) confirmed that both posttranscriptional regulation control points were used by a large portion of the transcriptome; (ii) identified a number of cis-acting features within the coding or noncoding sequences that correlate with splicing or translation state; and (iii) revealed correlation between each regulation mechanism and gene function. Our transcriptome-wide surveys have highlighted target genes transcripts of which are probably under extensive posttranscriptional regulation during flower development.

Finally, we reported the finding of a large number of polysome-associated ncRNAs. About one-third of all annotated ncRNA in the Arabidopsis genome were observed co-purified with polysomes. Coding capacity analysis confirmed that most of them are real ncRNA without conserved ORFs. The group of polysome-associated ncRNA reported in this study is a potential new addition to the expanding riboregulator catalog; they could have roles in translational regulation during early flower development.

Determining both the expression levels of mRNA and the regulation of its translation is important in understanding specialized cell functions. In this study, we describe both the expression profiles of cells within spatiotemporal domains of the Arabidopsis thaliana flower and the post-transcriptional regulation of these mRNAs, at nucleotide resolution. We express a tagged ribosomal protein under the promoters of three master regulators of flower development. By precipitating tagged polysomes, we isolated cell type-specific mRNAs that are probably translating, and quantified those mRNAs through deep sequencing. Cell type comparisons identified known cell-specific transcripts and uncovered many new ones, from which we inferred cell type-specific hormone responses, promoter motifs and coexpressed cognate binding factor candidates, and splicing isoforms. By comparing translating mRNAs with steady-state overall transcripts, we found evidence for widespread post-transcriptional regulation at both the intron splicing and translational stages. Sequence analyses identified structural features associated with each step. Finally, we identified a new class of noncoding RNAs associated with polysomes. Findings from our profiling lead to new hypotheses in the understanding of flower development.

Relationship between gene compactness and base composition in rice and human genome

In human, highly expressed genes contain shorter and fewer introns and these have been attributed to selection for economy in transcription and translation. On the other hand, in plants, it has been shown that highly expressed genes tend to be longer than lowly expressed genes. Here, in this study, we analyzed compositional influence on genome organization in both rice and human. We demonstrated that, in GC rich rice genes, highly expressed genes are less compact than lowly expressed genes. In GC-poor class, there is no difference in gene compactness between highly and lowly expressed genes. However, the scenario is different for human as there is no influence of GC composition on gene compactness due to their expression levels. We also reported that, highly expressed rice GC-rich pre-mRNA tend to form less stable secondary structure than that of lowly expressed genes. However, on removing intronic sequences, highly expressed mRNA form a stable secondary structure as compared to lowly expressed GC-rich genes. We suggest that in GC-rich rice genes long introns are under selection for enhancing transcriptional efficiency by modulating pre-mRNA secondary structural stability. Thus evolutionary mechanisms behind genome organization are different between these two genomes (human and rice).

High-resolution mapping of epigenetic modifications of the rice genome uncovers interplay between DNA methylation, histone methylation, and gene expression

We present high-resolution maps of DNA methylation and H3K4 di- and trimethylation of two entire chromosomes and two fully sequenced centromeres in rice (Oryza sativa) shoots and cultured cells. This analysis reveals combinatorial interactions between these epigenetic modifications and chromatin structure and gene expression. Cytologically densely stained heterochromatin had less H3K4me2 and H3K4me3 and more methylated DNA than the less densely stained euchromatin, whereas centromeres had a unique epigenetic composition. Most transposable elements had highly methylated DNA but no H3K4 methylation, whereas more than half of protein-coding genes had both methylated DNA and di- and/or trimethylated H3K4. Methylation of DNA but not H3K4 was correlated with suppressed transcription. By contrast, when both DNA and H3K4 were methylated, transcription was only slightly reduced. Transcriptional activity was positively correlated with the ratio of H3K4me3/H3K4me2: genes with predominantly H3K4me3 were actively transcribed, whereas genes with predominantly H3K4me2 were transcribed at moderate levels. More protein-coding genes contained all three modifications, and more transposons contained DNA methylation in shoots than cultured cells. Differential epigenetic modifications correlated to tissue-specific expression between shoots and cultured cells. Collectively, this study provides insights into the rice epigenomes and their effect on gene expression and plant development.

Retrotransposon insertion polymorphisms in six rice genes and their evolutionary history

Retrotransposons are abundant in higher plant genomes. Although retrotransposons associated with plant genes have been identified, little is known about their evolutionary conservation at the level of species and subspecies. In the present study, we investigated the phylogenetic distribution of long terminal repeat (LTR) retrotransposon, long interspersed nuclear element (LINE) and short interspersed nuclear element (SINE) insertions in six genes in 95 cultivated and wild rice genotypes. These six genes are likely to be functional based on nonsynonymous (Ka) to synonymous (Ks) substitution ratios which were found to be significantly <1. Different conservation patterns of these retrotransposons in genes were observed in cultivated and wild rice species. Four out of seven retrotransposon insertions appear to predate the ancestral Oryza AA genome. Two of these insertions in genes 4 and 5 occurred early in the evolutionary history of Oryza. Two retrotransposon insertions in gene 1 arose after the divergence of Asian cultivated rice from its wild ancestor. Furthermore, the retrotransposon insertion in gene 3 appears to have occurred in the ancestral lineage leading to temperate japonicas. Conservation of retrotransposon insertions in genes in specific groups, species, and lineages might be related to their specific function.

Analysis of genes associated with retrotransposons in the rice genome

Retrotransposons comprise a significant fraction of the rice genome. Despite their prevalence, the effects of retrotransposon insertions are not well understood, especially with regard to how they affect the expression of genes. In this study, we identified one-sixth of rice genes as being associated with retrotransposons, with insertions either in the gene itself or within its putative promoter region. Among genes with insertions in the promoter region, the likelihood of the gene being expressed was shown to be directly proportional to the distance of the retrotransposon from the translation start site. In addition, retrotransposon insertions in the transcribed region of the gene were found to be positively correlated with the presence of alternative splicing forms. Furthermore, preferential association of retrotransposon insertions with genes in several functional classes was identified. Some of the retrotransposons that are part of full-length cDNA (fl-cDNA) contribute splice sites and give rise to novel exons. Several interesting trends concerning the effects of retrotransposon insertions on gene expression were identified. Taken together, our data suggests that retrotransposon association with genes have a role in gene regulation. The data presented in this study provides a foundation for experimental studies to determine the role of retrotransposons in gene regulation.