Transcription of nuclear organellar DNA in a model plant system

Plastid DNA is a major source of nuclear genome complexity and of RNA genes in the orphan crop moringa

BACKGROUND

Unlike Transposable Elements (TEs) and gene/genome duplication, the role of the so-called nuclear plastid DNA sequences (NUPTs) in shaping the evolution of genome architecture and function remains poorly studied. We investigate here the functional and evolutionary fate of NUPTs in the orphan crop Moringa oleifera (moringa), featured by the highest fraction of plastid DNA found so far in any plant genome, focusing on (i) any potential biases in their distribution in relation to specific nuclear genomic features, (ii) their contribution to the emergence of new genes and gene regions, and (iii) their impact on the expression of target nuclear genes.

RESULTS

In agreement with their potential mutagenic effect, NUPTs are underrepresented among structural genes, although their overall transcription levels and broadness were only lower when involved exonic regions; the occurrence of plastid DNA generally did not result in a broader expression, except among those affected in introns by older NUPTs. In contrast, we found a strong enrichment of NUPTs among specific superfamilies of retrotransposons and several classes of RNA genes, including those participating in the protein biosynthetic machinery (i.e., rRNA and tRNA genes) and a specific class of regulatory RNAs. A significant fraction of NUPT RNA genes was found to be functionally expressed, thus potentially contributing to the nuclear pool.

CONCLUSIONS

Our results complete our view of the molecular factors driving the evolution of nuclear genome architecture and function, and support plastid DNA in moringa as a major source of (i) genome complexity and (ii) the nuclear pool of RNA genes.

Cryptic promoter activation occurs by at least two different mechanisms in the Arabidopsis genome

In gene-trap screening of plant genomes, promoterless reporter constructs are often expressed without trapping of annotated gene promoters. The molecular basis of this phenomenon, which has been interpreted as the trapping of cryptic promoters, is poorly understood. Here, we found that cryptic promoter activation occurs by at least two different mechanisms using Arabidopsis gene-trap lines in which a firefly luciferase (LUC) open reading frame (ORF) without an apparent promoter sequence was expressed from intergenic regions: one mechanism is 'cryptic promoter capturing', in which the LUC ORF captured pre-existing promoter-like chromatin marked by H3K4me3 and H2A.Z, and the other is 'promoter de novo origination', in which the promoter chromatin was newly formed near the 5' end of the inserted LUC ORF. The latter finding raises a question as to how the inserted LUC ORF sequence is involved in this phenomenon. To examine this, we performed a model experiment with chimeric LUC genes in transgenic plants. Using Arabidopsis psaH1 promoter-LUC constructs, we found that the functional core promoter region, where transcription start sites (TSSs) occur, cannot simply be determined by the upstream nor core promoter sequences; rather, its positioning proximal to the inserted LUC ORF sequence was more critical. This result suggests that the insertion of the coding sequence alters the local distribution of TSSs in the plant genome. The possible impact of the two types of cryptic promoter activation mechanisms on plant genome evolution and endosymbiotic gene transfer is discussed.

A Repertory of Rearrangements and the Loss of an Inverted Repeat Region in Passiflora Chloroplast Genomes

Chloroplast genomes (cpDNA) in angiosperms are usually highly conserved. Although rearrangements have been observed in some lineages, such as Passiflora, the mechanisms that lead to rearrangements are still poorly elucidated. In the present study, we obtained 20 new chloroplast genomes (18 species from the genus Passiflora, and Dilkea retusa and Mitostemma brevifilis from the family Passifloraceae) in order to investigate cpDNA evolutionary history in this group. Passiflora cpDNAs vary in size considerably, with ∼50 kb between shortest and longest. Large inverted repeat (IR) expansions were identified, and at the extreme opposite, the loss of an IR was detected for the first time in Passiflora, a rare event in angiosperms. The loss of an IR region was detected in Passiflora capsularis and Passiflora costaricensis, a species in which occasional biparental chloroplast inheritance has previously been reported. A repertory of rearrangements such as inversions and gene losses were detected, making Passiflora one of the few groups with complex chloroplast genome evolution. We also performed a phylogenomic study based on all the available cp genomes and our analysis implies that there is a need to reconsider the taxonomic classifications of some species in the group.

Kozak Sequence Acts as a Negative Regulator for De Novo Transcription Initiation of Newborn Coding Sequences in the Plant Genome

The manner in which newborn coding sequences and their transcriptional competency emerge during the process of gene evolution remains unclear. Here, we experimentally simulated eukaryotic gene origination processes by mimicking horizontal gene transfer events in the plant genome. We mapped the precise position of the transcription start sites (TSSs) of hundreds of newly introduced promoterless firefly luciferase (LUC) coding sequences in the genome of Arabidopsis thaliana cultured cells. The systematic characterization of the LUC-TSSs revealed that 80% of them occurred under the influence of endogenous promoters, while the remainder underwent de novo activation in the intergenic regions, starting from pyrimidine-purine dinucleotides. These de novo TSSs obeyed unexpected rules; they predominantly occurred ∼100 bp upstream of the LUC inserts and did not overlap with Kozak-containing putative open reading frames (ORFs). These features were the output of the immediate responses to the sequence insertions, rather than a bias in the screening of the LUC gene function. Regarding the wild-type genic TSSs, they appeared to have evolved to lack any ORFs in their vicinities. Therefore, the repulsion by the de novo TSSs of Kozak-containing ORFs described above might be the first selection gate for the occurrence and evolution of TSSs in the plant genome. Based on these results, we characterized the de novo type of TSS identified in the plant genome and discuss its significance in genome evolution.

Deficiency of mitoribosomal S10 protein affects translation and splicing in Arabidopsis mitochondria

The ribosome is not only a protein-making machine, but also a regulatory element in protein synthesis. This view is supported by our earlier data showing that Arabidopsis mitoribosomes altered due to the silencing of the nuclear RPS10 gene encoding mitochondrial ribosomal protein S10 differentially translate mitochondrial transcripts compared with the wild-type. Here, we used ribosome profiling to determine the contribution of transcriptional and translational control in the regulation of protein synthesis in rps10 mitochondria compared with the wild-type ones. Oxidative phosphorylation system proteins are preferentially synthesized in wild-type mitochondria but this feature is lost in the mutant. The rps10 mitoribosomes show slightly reduced translation efficiency of most respiration-related proteins and at the same time markedly more efficiently synthesize ribosomal proteins and MatR and TatC proteins. The mitoribosomes deficient in S10 protein protect shorter transcript fragments which exhibit a weaker 3-nt periodicity compared with the wild-type. The decrease in the triplet periodicity is particularly drastic for genes containing introns. Notably, splicing is considerably less effective in the mutant, indicating an unexpected link between the deficiency of S10 and mitochondrial splicing. Thus, a shortage of the mitoribosomal S10 protein has wide-ranging consequences on mitochondrial gene expression.

Pod-shattering characteristics differences between two groups of soybeans are associated with specific changes in gene expression

Soybean is an economically important leguminous crop, and pod dehiscence of soybean could cause huge yield loss. In this study, we measured fruit-cracking forces and percentages of dehisced pods for ten soybean accessions, then separated them into two groups as shattering-sensitive (SS) and shattering-resistant (SR) soybeans. Pod transcriptomes from these two groups were analyzed, and 225 differentially expressed genes (DEGs) were identified between SS and SR soybeans. Some of these DEGs have been previously reported to be associated with pod dehiscence in soybean. The expression patterns of selected DEGs were validated by real-time quantitative reverse transcription PCR, which confirmed the expression changes found in RNA-seq analysis. We also de novo identified 246 soybean pod-long intergenic ncRNAs (lincRNAs), 401 intronic lncRNAs, and 23 antisense lncRNAs from these transcriptomes. Furthermore, genes and lincRNAs co-expression network analysis showed that there are distinct expression patterns between SS and SR soybeans in some co-expression modules. In conclusion, we systematically investigated potential genes and molecular pathways as candidates for differences in soybean pod dehiscence and will provide a useful resource for molecular breeding of soybeans.

Intergenomic gene transfer in diploid and allopolyploid Gossypium

BACKGROUND

Intergenomic gene transfer (IGT) between nuclear and organellar genomes is a common phenomenon during plant evolution. Gossypium is a useful model to evaluate the genomic consequences of IGT for both diploid and polyploid species. Here, we explore IGT among nuclear, mitochondrial, and plastid genomes of four cotton species, including two allopolyploids and their model diploid progenitors (genome donors, G. arboreum: A₂ and G. raimondii: D₅).

RESULTS

Extensive IGT events exist for both diploid and allotetraploid cotton (Gossypium) species, with the nuclear genome being the predominant recipient of transferred DNA followed by the mitochondrial genome. The nuclear genome has integrated 100 times more foreign sequences than the mitochondrial genome has in total length. In the nucleus, the integrated length of chloroplast DNA (cpDNA) was between 1.87 times (in diploids) to nearly four times (in allopolyploids) greater than that of mitochondrial DNA (mtDNA). In the mitochondrion, the length of nuclear DNA (nuDNA) was typically three times than that of cpDNA. Gossypium mitochondrial genomes integrated three nuclear retrotransposons and eight chloroplast tRNA genes, and incorporated chloroplast DNA prior to divergence between the diploids and allopolyploid formation. For mitochondrial chloroplast-tRNA genes, there were 2-6 bp conserved microhomologies flanking their insertion sites across distantly related genera, which increased to 10 bp microhomologies for the four cotton species studied. For organellar DNA sequences, there are source hotspots, e.g., the atp6-trnW intergenic region in the mitochondrion and the inverted repeat region in the chloroplast. Organellar DNAs in the nucleus were rarely expressed, and at low levels. Surprisingly, there was asymmetry in the survivorship of ancestral insertions following allopolyploidy, with most numts (nuclear mitochondrial insertions) decaying or being lost whereas most nupts (nuclear plastidial insertions) were retained.

CONCLUSIONS

This study characterized and compared intracellular transfer among nuclear and organellar genomes within two cultivated allopolyploids and their ancestral diploid cotton species. A striking asymmetry in the fate of IGTs in allopolyploid cotton was discovered, with numts being preferentially lost relative to nupts. Our results connect intergenomic gene transfer with allotetraploidy and provide new insight into intracellular genome evolution.

Starch content differences between two sweet potato accessions are associated with specific changes in gene expression

Sweet potato (Ipomoea batatas (L.) Lam.) is one of the most important root crops in the world. Initial formation and development of storage roots (SRs) are key factors affecting its yields. In order to study the molecular mechanism and regulatory networks of the SRs development process, we have analyzed root transcriptomes between the high and low starch content sweet potato accessions at three different developmental stages. In this study, we assembled 46,840 unigenes using Illumina paired-end sequencing reads and identified differentially expressed genes (DEGs) between two accessions. The numbers of DEGs were increased with the development of SRs, indicating that the difference between two accessions is enlarging with the maturation. DEGs were mainly enriched in starch biosynthesis, plant hormones regulatory, and genetic information processing pathways. Then, expression patterns of DEGs that are most significant and starch biosynthesis related were validated using qRT-PCR. Our results provide valuable resources to future study on molecular mechanisms of SRs development and candidate genes for starch content improvement in sweet potato.

The Roles of Mitochondrion in Intergenomic Gene Transfer in Plants: A Source and a Pool

Intergenomic gene transfer (IGT) is continuous in the evolutionary history of plants. In this field, most studies concentrate on a few related species. Here, we look at IGT from a broader evolutionary perspective, using 24 plants. We discover many IGT events by assessing the data from nuclear, mitochondrial and chloroplast genomes. Thus, we summarize the two roles of the mitochondrion: a source and a pool. That is, the mitochondrion gives massive sequences and integrates nuclear transposons and chloroplast tRNA genes. Though the directions are opposite, lots of likenesses emerge. First, mitochondrial gene transfer is pervasive in all 24 plants. Second, gene transfer is a single event of certain shared ancestors during evolutionary divergence. Third, sequence features of homologies vary for different purposes in the donor and recipient genomes. Finally, small repeats (or micro-homologies) contribute to gene transfer by mediating recombination in the recipient genome.

Seed weight differences between wild and domesticated soybeans are associated with specific changes in gene expression

KEY MESSAGE

Our study systematically explored potential genes and molecular pathways as candidates for differences in seed weight resulting from soybean domestication. In addition, potential contributions of lncRNAs to seed weight were also investigated. Soybeans have a long history of domestication in China, and there are several significant phenotypic differences between cultivated and wild soybeans, for example, seeds of cultivars are generally larger and heavier than those from wild accessions. We analyzed seed transcriptomes from thirteen soybean samples, including six landraces and seven wild accessions using strand-specific RNA sequencing. Differentially expressed genes related to seed weight were identified, and some of their homologs were associated with seed development in Arabidopsis. We also identified 1251 long intergenic noncoding RNAs (lincRNAs), 243 intronic RNAs and 81 antisense lncRNAs de novo from these soybean transcriptomes. We then profiled the expression patterns of lncRNAs in cultivated and wild soybean seeds, and found that transcript levels of a number of lncRNAs were sample-specific. Moreover, gene transcript and lincRNA co-expression network analysis showed that some soybean lincRNAs might have functional roles as they were hubs of co-expression modules. In conclusion, this study systematically explored potential genes and molecular pathways as candidates for differences in seed weight resulting from soybean domestication, and will provide a useful future resource for molecular breeding of soybeans.

Structure of Plastid Genomes of Photosynthetic Eukaryotes

This review presents current views on the plastid genomes of higher plants and summarizes data on the size, structural organization, gene content, and other features of plastid DNAs. Special emphasis is placed on the properties of organization of land plant plastid genomes (nucleoids) that distinguish them from bacterial genomes. The prospects of genetic engineering of chloroplast genomes are discussed.

Transposable elements (TEs) contribute to stress-related long intergenic noncoding RNAs in plants

Noncoding RNAs have been extensively described in plant and animal transcriptomes by using high-throughput sequencing technology. Of these noncoding RNAs, a growing number of long intergenic noncoding RNAs (lincRNAs) have been described in multicellular organisms, however the origins and functions of many lincRNAs remain to be explored. In many eukaryotic genomes, transposable elements (TEs) are widely distributed and often account for large fractions of plant and animal genomes yet the contribution of TEs to lincRNAs is largely unknown. By using strand-specific RNA-sequencing, we profiled the expression patterns of lincRNAs in Arabidopsis, rice and maize, and identified 47 611 and 398 TE-associated lincRNAs (TE-lincRNAs), respectively. TE-lincRNAs were more often derived from retrotransposons than DNA transposons and as retrotransposon copy number in both rice and maize genomes so did TE-lincRNAs. We validated the expression of these TE-lincRNAs by strand-specific RT-PCR and also demonstrated tissue-specific transcription and stress-induced TE-lincRNAs either after salt, abscisic acid (ABA) or cold treatments. For Arabidopsis TE-lincRNA11195, mutants had reduced sensitivity to ABA as demonstrated by longer roots and higher shoot biomass when compared to wild-type. Finally, by altering the chromatin state in the Arabidopsis chromatin remodelling mutant ddm1, unique lincRNAs including TE-lincRNAs were generated from the preceding untranscribed regions and interestingly inherited in a wild-type background in subsequent generations. Our findings not only demonstrate that TE-associated lincRNAs play important roles in plant abiotic stress responses but lincRNAs and TE-lincRNAs might act as an adaptive reservoir in eukaryotes.

Whole-Genome Sequencing Reveals a New Genospecies of Methylobacterium sp. GXS13, Isolated from Vitis vinifera L. Xylem Sap

The whole-genome sequence of a new genospecies of Methylobacterium sp., named GXS13 and isolated from grapevine xylem sap, is reported and demonstrates potential for methylotrophy, cytokinin synthesis, and cell wall modification. In addition, biosynthetic gene clusters were identified for cupriachelin, carotenoid, and acyl-homoserine lactone using the antiSMASH server.

The 'dark matter' in the plant genomes: non-coding and unannotated DNA sequences associated with open chromatin

Sequencing of complete plant genomes has become increasingly more routine since the advent of the next-generation sequencing technology. Identification and annotation of large amounts of noncoding but functional DNA sequences, including cis-regulatory DNA elements (CREs), have become a new frontier in plant genome research. Genomic regions containing active CREs bound to regulatory proteins are hypersensitive to DNase I digestion and are called DNase I hypersensitive sites (DHSs). Several recent DHS studies in plants illustrate that DHS datasets produced by DNase I digestion followed by next-generation sequencing (DNase-seq) are highly valuable for the identification and characterization of CREs associated with plant development and responses to environmental cues. DHS-based genomic profiling has opened a door to identify and annotate the 'dark matter' in sequenced plant genomes.