The chloroplast and mitochondrial C-to-U RNA editing in Arabidopsis thaliana shows signals of adaptation

Convergent evolution of allele-specific gene expression that leads to non-small cell lung cancer in different human populations

Phenotypical innovations during evolution are caused by novel mutations, which are usually heterozygous at the beginning. The gene expressions on two alleles of these mutation sites are not necessarily identical, leading to flexible allele-specific regulation in cell systems. We retrieve the transcriptome data of normal and non-small cell lung cancer (NSCLC) tissues from 47 African Americans (AA) and 50 European Americans (EA). We analyze the differentially expressed genes (DEGs) in NSCLC as well as the tumor-specific mutations. Expression and mutation profiles show convergent evolution in AA and EA populations. The tumor-specific mutations are poorly overlapped, but many of them are located in the same genes, mainly oncogenes and tumor suppressor genes. The DEGs in tumors are majorly caused by the mutated alleles rather than normal alleles. The relative expressions of mutated alleles are highly correlated between AA and EA. The differential expression in NSCLC is predominantly mediated by the mutated alleles on heterozygous sites. This molecular mechanism underlying NSCLC oncogenesis is conserved across different human populations, exhibiting convergent evolution. We present this novel angle that differential expression analysis should be performed separately for different alleles. Our ideas should greatly benefit the cancer community.

Assembly and comparative analysis of the first complete mitochondrial genome of the invasive water hyacinth, Eichhornia crassipes

Eichhornia crassipes is an aquatic plant in tropical and subtropical regions, renowned for its notorious invasive tendencies. In this study, we assembled the complete mitogenome of E. crassipes into a single circle molecule of 397,361 bp. The mitogenome has 58 unique genes, including 37 protein-coding genes (PCGs), 18 tRNA genes, three rRNA genes, and 47 % GC content. Sixteen (6.93 %) homologous fragments, ranging from 31 bp to 8548 bp, were identified, indicating the transfer of genetic material from chloroplasts to mitochondria. In addition, we detected positive selection in six PCGs (ccmB, ccmC, ccmFC, nad3, nad4 and sdh4), along with the identification of 782 RNA editing sites across 37 mt-PCGs. These findings suggest a potential contribution to the robust adaptation of this invasive plant to the stressful environment. Lastly, we inferred that phylogenetic conflicts of E. crassipes between the plastome and mitogenome may be attributed to the difference in nucleotide substitution rates between the two organelle genomes. In conclusion, our study provided vital genomic resources for further understanding the invasive mechanism of this species and exploring the dynamic evolution of mitogenomes within the monocot clade.

RNA editing in host lncRNAs as potential modulator in SARS-CoV-2 variants-host immune response dynamics

Both host and viral RNA editing plays a crucial role in host's response to infection, yet our understanding of host RNA editing remains limited. In this study of in-house generated RNA sequencing (RNA-seq) data of 211 hospitalized COVID-19 patients with PreVOC, Delta, and Omicron variants, we observed a significant differential editing frequency and patterns in long non-coding RNAs (lncRNAs), with Delta group displaying lower RNA editing compared to PreVOC/Omicron patients. Notably, multiple transcripts of UGDH-AS1 and NEAT1 exhibited high editing frequencies. Expression of ADAR1/APOBEC3A/APOBEC3G and differential abundance of repeats were possible modulators of differential editing across patient groups. We observed a shift in crucial infection-related pathways wherein the pathways were downregulated in Delta compared to PreVOC and Omicron. Our genomics-based evidence suggests that lncRNA editing influences stability, miRNA binding, and expression of both lncRNA and target genes. Overall, the study highlights the role of lncRNAs and how editing within host lncRNAs modulates the disease severity.

The first complete chloroplast genome of Thalictrum fargesii: insights into phylogeny and species identification

Introduction

Thalictrum fargesii is a medicinal plant belonging to the genus Thalictrum of the Ranunculaceae family and has been used in herbal medicine in the Himalayan regions of China and India. This species is taxonomically challenging because of its morphological similarities to other species within the genus. Thus, herbal drugs from this species are frequently adulterated, substituted, or mixed with other species, thereby endangering consumer safety.

Methods

The present study aimed to sequence and assemble the entire chloroplast (cp) genome of T. fargesii using the Illumina HiSeq 2500 platform to better understand the genomic architecture, gene composition, and phylogenetic relationships within the Thalictrum.

Results and discussion

The cp genome was 155,929 bp long and contained large single-copy (85,395 bp) and small single-copy (17,576 bp) regions that were segregated by a pair of inverted repeat regions (26,479 bp) to form a quadripartite structure. The cp genome contains 133 genes, including 88 protein-coding genes (PCGs), 37 tRNA genes, and 8 rRNA genes. Additionally, this genome contains 64 codons that encode 20 amino acids, the most preferred of which are alanine and leucine. We identified 68 SSRs, 27 long repeats, and 242 high-confidence C-to-U RNA-editing sites in the cp genome. Moreover, we discovered seven divergent hotspot regions in the cp genome of T. fargesii, among which ndhD-psaC and rpl16-rps3 may be useful for developing molecular markers for identifying ethnodrug species and their contaminants. A comparative study with eight other species in the genus revealed that pafI and rps19 had highly variable sites in the cp genome of T. fargesii. Additionally, two special features, (i) the shortest length of the ycf1 gene at the IRA-SSC boundary and (ii) the distance between the rps19 fragment and trnH at the IRA-LSC junction, distinguish the cp genome of T. fargesii from those of other species within the genus. Furthermore, phylogenetic analysis revealed that T. fargesii was closely related to T. tenue and T. petaloidium.

Conclusion

Considering all these lines of evidence, our findings offer crucial molecular and evolutionary information that could play a significant role in further species identification, evolution, and phylogenetic studies on T. fargesii.

The first A-to-I RNA editome of hemipteran species Coridius chinensis reveals overrepresented recoding and prevalent intron editing in early-diverging insects

BACKGROUND

Metazoan adenosine-to-inosine (A-to-I) RNA editing resembles A-to-G mutation and increases proteomic diversity in a temporal-spatial manner, allowing organisms adapting to changeable environment. The RNA editomes in many major animal clades remain unexplored, hampering the understanding on the evolution and adaptation of this essential post-transcriptional modification.

METHODS

We assembled the chromosome-level genome of Coridius chinensis belonging to Hemiptera, the fifth largest insect order where RNA editing has not been studied yet. We generated ten head RNA-Seq libraries with DNA-Seq from the matched individuals.

RESULTS

We identified thousands of high-confidence RNA editing sites in C. chinensis. Overrepresentation of nonsynonymous editing was observed, but conserved recoding across different orders was very rare. Under cold stress, the global editing efficiency was down-regulated and the general transcriptional processes were shut down. Nevertheless, we found an interesting site with "conserved editing but non-conserved recoding" in potassium channel Shab which was significantly up-regulated in cold, serving as a candidate functional site in response to temperature stress.

CONCLUSIONS

RNA editing in C. chinensis largely recodes the proteome. The first RNA editome in Hemiptera indicates independent origin of beneficial recoding during insect evolution, which advances our understanding on the evolution, conservation, and adaptation of RNA editing.

Liver cancer-specific mutations in functional domains of ADAR2 lead to the elevation of coding and non-coding RNA editing in multiple tumor-related genes

Mutation is the major cause of phenotypic innovations. Apart from DNA mutations, the alteration on RNA such as the ADAR-mediated A-to-I RNA editing could also shape the phenotype. These two layers of variations have not been systematically combined to study their collective roles in cancers. We collected the high-quality transcriptomes of ten hepatocellular carcinoma (HCC) and the matched control samples. We systematically identified HCC-specific mutations in the exonic regions and profiled the A-to-I RNA editome in each sample. All ten HCC samples had mutations in the CDS of ADAR2 gene (dsRNA-binding domain or catalytic domain). The consequence of these mutations converged to the elevation of ADAR2 efficiency as reflected by the global increase of RNA editing levels in HCC. The up-regulated editing sites (UES) were enriched in the CDS and UTR of oncogenes and tumor suppressor genes (TSG), indicating the possible roles of these target genes in HCC oncogenesis. We present the mutation-ADAR2-UES-oncogene/TSG-HCC axis that explains how mutations at different layers would finally lead to abnormal phenotype. In the light of central dogma, our work provides novel insights into how to fully take advantage of the transcriptome data to decipher the consequence of mutations.

SARS-CoV-2 continuously optimizes its codon usage to adapt to human lung environment

Viruses need to utilize the resources from host cells to reproduce themselves. RNA translation rate, which is largely determined by codon usage, is the rate-limiting step across the life cycle of viruses. Adapting to the codon usage of hosts would help virus better proliferate. We retrieved the time-course mutation profile of millions of world-wide SARS-CoV-2 sequences. For synonymous mutations, we defined whether a mutation elevate or reduce the relative synonymous codon usage (RSCU). We found that if a synonymous mutation in SARS-CoV-2 increases the RSCU (calculated from human lungs), denoted as delta RSCU > 0, then this mutation is positively selected because the allele frequency (AF) of this mutation increases with time, and vice versa. The results suggest that in SARS-CoV-2, the synonymous mutations that increase codon optimality are beneficial to the virus and are favored by natural selection. For the first time, we used the dynamics of allele frequency to demonstrate that SARS-CoV-2 is continuously optimizing its codon usage to adapt to human lungs. Nevertheless, adaptation to other human tissues cannot be excluded. These results warn us that under this global pandemic, synonymous mutations in SARS-CoV-2 should not be automatically ignored since they indeed change the fitness of the virus.

A-to-I RNA editing shows dramatic up-regulation in osteosarcoma and broadly regulates tumor-related genes by altering microRNA target regions

A-to-I RNA editing is a prevalent type of RNA modification in animals. The dysregulation of RNA editing has led to multiple human cancers. However, the role of RNA editing has never been studied in osteosarcoma, a complex bone cancer with unknown molecular basis. We retrieved the RNA-sequencing data from 24 primary osteosarcoma patients and 3 healthy controls. We systematically profiled the RNA editomes in these samples and quantitatively identified reliable differential editing sites (DES) between osteosarcoma and normal samples. RNA editing efficiency is dramatically increased in osteosarcoma, presumably due to the significant up-regulation of editing enzymes ADAR1 and ADAR2. Up-regulated DES in osteosarcoma are enriched in 3'UTRs. Strikingly, such 3'UTR sites are further enriched in microRNA binding regions of gene EMP2 and other oncogenes, abolishing the microRNA suppression on target genes. Accordingly, the expression of these tumor-promoting genes is elevated in osteosarcoma. There might be an RNA editing-dependent pathway leading to osteosarcoma. We expanded our knowledge on the potential roles of RNA editing in oncogenesis. Based on these molecular features, our work is valuable for future prognosis and diagnosis of osteosarcoma.

Analysis of the RNA Editing Sites and Orthologous Gene Function of Transcriptome and Chloroplast Genomes in the Evolution of Five Deutzia Species

In this study, the chloroplast genomes and transcriptomes of five Deutzia genus species were sequenced, characterized, combined, and analyzed. A phylogenetic tree was constructed, including 32 other chloroplast genome sequences of Hydrangeoideae species. The results showed that the five Deutzia chloroplast genomes were typical circular genomes 156,860-157,025 bp in length, with 37.58-37.6% GC content. Repeat analysis showed that the Deutzia species had 41-45 scattered repeats and 199-201 simple sequence repeats. Comparative genomic and pi analyses indicated that the genomes are conservative and that the gene structures are stable. According to the phylogenetic tree, Deutzia species appear to be closely related to Kirengeshoma palmata and Philadelphus. By combining chloroplast genomic and transcriptomic analyses, 29-31 RNA editing events and 163-194 orthologous genes were identified. The ndh, rpo, rps, and atp genes had the most editing sites, and all RNA editing events were of the C-to-U type. Most of the orthologous genes were annotated to the chloroplast, mitochondria, and nucleus, with functions including energy production and conversion, translation, and protein transport. Genes related to the biosynthesis of monoterpenoids and flavonoids were also identified from the transcriptome of Deutzia spp. Our results will contribute to further studies of the genomic information and potential uses of the Deutzia spp.

Differential adaptive RNA editing signals between insects and plants revealed by a new measurement termed haplotype diversity

BACKGROUND

C-to-U RNA editing in plants is believed to confer its evolutionary adaptiveness by reversing unfavorable DNA mutations. This "restorative hypothesis" has not yet been tested genome-wide. In contrast, A-to-I RNA editing in insects like Drosophila and honeybee is already known to benefit the host by increasing proteomic diversity in a spatial-temporal manner (namely "diversifying hypothesis").

METHODS

We profiled the RNA editomes of multiple tissues of Arabidopsis thaliana, Drosophila melanogaster, and Apis melifera. We unprecedentedly defined the haplotype diversity (HD) of RNA molecules based on nonsynonymous editing events (recoding sites).

RESULTS

Signals of adaptation is confirmed in Arabidopsis by observing higher frequencies and levels at nonsynonymous editing sites over synonymous sites. Compared to A-to-I recoding sites in Drosophila, the C-to-U recoding sites in Arabidopsis show significantly lower HD, presumably due to the stronger linkage between C-to-U events.

CONCLUSIONS

C-to-U RNA editing in Arabidopsis is adaptive but it is not designed for diversifying the proteome like A-to-I editing in Drosophila. Instead, C-to-U recoding sites resemble DNA mutations. Our observation supports the restorative hypothesis of plant C-to-U editing which claims that editing is used for fixing unfavorable genomic sequences.

Comparative analysis of the chloroplast genomes of Rosa species and RNA editing analysis

BACKGROUND

The genus Rosa (Rosaceae) contains approximately 200 species, most of which have high ecological and economic values. Chloroplast genome sequences are important for studying species differentiation, phylogeny, and RNA editing.

RESULTS

In this study, the chloroplast genomes of three Rosa species, Rosa hybrida, Rosa acicularis, and Rosa rubiginosa, were assembled and compared with other reported Rosa chloroplast genomes. To investigate the RNA editing sites in R. hybrida (commercial rose cultivar), we mapped RNA-sequencing data to the chloroplast genome and analyzed their post-transcriptional features. Rosa chloroplast genomes presented a quadripartite structure and had highly conserved gene order and gene content. We identified four mutation hotspots (ycf3-trnS, trnT-trnL, psbE-petL, and ycf1) as candidate molecular markers for differentiation in the Rosa species. Additionally, 22 chloroplast genomic fragments with a total length of 6,192 bp and > 90% sequence similarity with their counterparts were identified in the mitochondrial genome, representing 3.96% of the chloroplast genome. Phylogenetic analysis including all sections and all subgenera revealed that the earliest divergence in the chloroplast phylogeny roughly distinguished species of sections Pimpinellifoliae and Rosa and subgenera Hulthemia. Moreover, DNA- and RNA-sequencing data revealed 19 RNA editing sites, including three synonymous and 16 nonsynonymous, in the chloroplast genome of R. hybrida that were distributed among 13 genes.

CONCLUSIONS

The genome structure and gene content of Rosa chloroplast genomes are similar across various species. Phylogenetic analysis based on the Rosa chloroplast genomes has high resolution. Additionally, a total of 19 RNA editing sites were validated by RNA-Seq mapping in R. hybrida. The results provide valuable information for RNA editing and evolutionary studies of Rosa and a basis for further studies on genomic breeding of Rosa species.

The efficacy and safety of SARS-CoV-2 vaccines mRNA1273 and BNT162b2 might be complicated by rampant C-to-U RNA editing

The SARS-CoV-2 RNA vaccines are smartly designed to increase the synonymous codon usage by introducing multiple U-to-C mutations. This design would elevate the translation efficiency of vaccine RNAs. However, we found evidence to reason that the designed cytidines might be converted to uridines again by C-to-U RNA deamination in host cells. This C-to-U mechanism might be a main factor that affects the efficacy and safety of RNA vaccines.

Evidence Supporting That C-to-U RNA Editing Is the Major Force That Drives SARS-CoV-2 Evolution

Mutations of DNA organisms are introduced by replication errors. However, SARS-CoV-2, as an RNA virus, is additionally subjected to rampant RNA editing by hosts. Both resources contributed to SARS-CoV-2 mutation and evolution, but the relative prevalence of the two origins is unknown. We performed comparative genomic analyses at intra-species (world-wide SARS-CoV-2 strains) and inter-species (SARS-CoV-2 and RaTG13 divergence) levels. We made prior predictions of the proportion of each mutation type (nucleotide substitution) under different scenarios and compared the observed versus the expected. C-to-T alteration, representing C-to-U editing, is far more abundant that all other mutation types. Derived allele frequency (DAF) as well as novel mutation rate of C-to-T are the highest in SARS-CoV-2 population, and C-T substitution dominates the divergence sites between SARS-CoV-2 and RaTG13. This is compelling evidence suggesting that C-to-U RNA editing is the major source of SARS-CoV-2 mutation. While replication errors serve as a baseline of novel mutation rate, the C-to-U editing has elevated the mutation rate for orders of magnitudes and accelerates the evolution of the virus.

Genome-wide investigation and functional analysis of RNA editing sites in wheat

Wheat is an important cereal and half of the world population consumed it. Wheat faces environmental stresses and different techniques (CRISPR, gene silencing, GWAS, etc.) were used to enhance its production but RNA editing (RESs) is not fully explored in wheat. RNA editing has a special role in controlling environmental stresses. The genome-wide identification and functional characterization of RESs in different types of wheat genotypes was done. We employed six wheat genotypes by RNA-seq analyses to achieve RESs. The findings revealed that RNA editing events occurred on all chromosomes equally. RNA editing sites were distributed randomly and 10–12 types of RESs were detected in wheat genotypes. Higher number of RESs were detected in drought-tolerant genotypes. A-to-I RNA editing (2952, 2977, 1916, 2576, 3422, and 3459) sites were also identified in six wheat genotypes. Most of the genes were found to be engaged in molecular processes after a Gene Ontology analysis. PPR (pentatricopeptide repeat), OZ1 (organelle zinc-finger), and MORF/RIP gene expression levels in wheat were also examined. Normal growth conditions diverge gene expression of these three different gene families, implying that normal growth conditions for various genotypes can modify RNA editing events and have an impact on gene expression levels. While the expression of PPR genes was not change. We used Variant Effect Predictor (VEP) to annotate RNA editing sites, and Local White had the highest RESs in the CDS region of the protein. These findings will be useful for prediction of RESs in other crops and will be helpful in drought tolerance development in wheat.

ORTHOSKIM: in silico sequence capture from genomic and transcriptomic libraries for phylogenomic and barcoding applications

Low-coverage whole genome shotgun sequencing (or genome skimming) has emerged as a cost-effective method for acquiring genomic data in non-model organisms. This method provides sequence information on chloroplast genome (cpDNA), mitochondrial genome (mtDNA) and nuclear ribosomal regions (rDNA), which are over-represented within cells. However, numerous bioinformatic challenges remain to accurately and rapidly obtain such data in organisms with complex genomic structures and rearrangements, in particular for mtDNA in plants or for cpDNA in some plant families. Here we introduce the pipeline ORTHOSKIM, which performs in silico capture of targeted sequences from genomic and transcriptomic libraries without assembling whole organelle genomes. ORTHOSKIM proceeds in three steps: 1) global sequence assembly, 2) mapping against reference sequences, and 3) target sequence extraction; importantly it also includes a range of quality control tests. Different modes are implemented to capture both coding and non-coding regions of cpDNA, mtDNA and rDNA sequences, along with predefined nuclear sequences (e.g. ultra-conserved elements) or collections of single-copy ortholog genes. Moreover, aligned DNA matrices are produced for phylogenetic reconstructions, by performing multiple alignments of the captured sequences. While ORTHOSKIM is suitable for any eukaryote, a case study is presented here, using 114 genome-skimming libraries and 4 RNAseq libraries obtained for two plant families, Primulaceae and Ericaceae, the latter being a well-known problematic family for cpDNA assemblies. ORTHOSKIM recovered with high success rates cpDNA, mtDNA and rDNA sequences, well suited to accurately infer evolutionary relationships within these families. ORTHOSKIM is released under a GPL-3 license and is available at: https://github.com/cpouchon/ORTHOSKIM.

Comparative genomic analysis of a naturally born serpentized pig reveals putative mutations related to limb and bone development

BACKGROUND

It is believed that natural selection acts on the phenotypical changes caused by mutations. Phenotypically, from fishes to amphibians to reptiles, the emergence of limbs greatly facilitates the landing of ancient vertebrates, but the causal mutations and evolutionary trajectory of this process remain unclear.

RESULTS

We serendipitously obtained a pig of limbless phenotype. Mutations specific to this handicapped pig were identified using genome re-sequencing and comparative genomic analysis. We narrowed down the causal mutations to particular chromosomes and even several candidate genes and sites, such like a mutation-containing codon in gene BMP7 (bone morphogenetic protein) which was conserved in mammals but variable in lower vertebrates.

CONCLUSIONS

We parsed the limbless-related mutations in the light of evolution. The limbless pig shows phenocopy of the clades before legs were evolved. Our findings might help deduce the emergence of limbs during vertebrate evolution and should be appealing to the broad community of human genetics and evolutionary biology.

Genome-wide expression changes mediated by A-to-I RNA editing correlate with hepatic oncogenesis

Background

Adenosine-to-inosine (A-to-I) RNA editing is one of the most prevalent RNA modifications in the animal kingdom. Since inosine is recognized as guanosines, the A-to-I process mimics A-to-G DNA mutations but can be controlled in a more flexible manner compared to DNA alterations.

Methods

We parsed the transcriptomes and translatomes of liver cancer and normal tissues from ten patients. We profiled the landscape of the A-to-I RNA editome in these samples and interrogated whether the A-to-I processes participated in the gene expression regulation in oncogenesis.

Results

Globally, editing activity was enhanced in all tumor samples compared to that in normal samples. Accordingly, expression of the gene encoding the RNA editing enzyme ADAR (adenosine deaminase acting on RNA) was elevated. Two intronic self-editing sites in ADAR mRNAs controlled its splicing pattern and may regulate its translation efficiency (TE). Moreover, the expression of oncogenes was generally upregulated in tumors, whereas tumor suppressor genes (TSG) were downregulated, possibly due to alterations to microRNA binding sites or RNA splicing defects caused by A-to-I editing.

Conclusions

A-to-I RNA editing plays a crucial role in the oncogenesis of liver cancer. ADAR regulates its own expression via self-editing, and it also affects global transcriptomes and translatomes involving cancer-related genes by RNA editing and changing their expression patterns.

Retrieving the deleterious mutations before extinction: genome-wide comparison of shared derived mutations in liver cancer and normal population

STUDY PURPOSE

Deleterious mutations would be rapidly purged from natural populations along with the extinction of their carriers. The currently observed mutations in existing species are mostly neutral. The inaccessibility of deleterious mutations impedes the functional studies on how these mutations affect the fitness at individual level.

STUDY DESIGN

The connection between the deleterious genotype and the non-adaptive phenotype could be bridged by sequencing the genome before extinction. Although this approach is no longer feasible for evolutionary biologists, it is feasible for cancer biologists by profiling the mutations in tumour samples which are so deleterious that the carriers hardly live.

RESULTS

By comparing the derived mutation profile between normal populations and patients with liver cancer, we found that the shared mutations, which are highly deleterious, are suppressed to low allele frequencies in normal populations and tissues, but show remarkably high frequency in tumours. The density of shared mutations is negatively correlated with gene conservation and expression levels.

CONCLUSIONS

Deleterious mutations are suppressed in functionally important genes as well as in normal populations. This work deepened our understanding on how natural selection act on deleterious mutations by analogising the cancer evolution to species evolution, which are essentially the same molecular process but at different time scales.

Synonymous mutations that regulate translation speed might play a non-negligible role in liver cancer development

Background

Synonymous mutations do not change the protein sequences. Automatically, they have been regarded as neutral events and are ignored in the mutation-based cancer studies. However, synonymous mutations will change the codon optimality, resulting in altered translational velocity.

Methods

We fully utilized the transcriptome and translatome of liver cancer and normal tissue from ten patients. We profiled the mutation spectrum and examined the effect of synonymous mutations on translational velocity.

Results

Synonymous mutations that increase the codon optimality significantly enhanced the translational velocity, and were enriched in oncogenes. Meanwhile, synonymous mutations decreasing codon optimality slowed down translation, and were enriched in tumor suppressor genes. These synonymous mutations significantly contributed to the translational changes in tumor samples compared to normal samples.

Conclusions

Synonymous mutations might play a role in liver cancer development by altering codon optimality and translational velocity. Synonymous mutations should no longer be ignored in the genome-wide studies.

Detection and Analysis of C-to-U RNA Editing in Rice Mitochondria-Encoded ORFs

Cytidine to uridine (C-to-U) RNA editing is an important type of substitutional RNA modification and is almost omnipresent in plant chloroplasts and mitochondria. In rice mitochondria, 491 C-to-U editing sites have been identified previously, and case studies have elucidated the function of several C-to-U editing sites in rice, but the functional consequence of most C-to-U alterations needs to be investigated further. Here, by means of Sanger sequencing and publicly available RNA-seq data, we identified a total of 569 C-to-U editing sites in rice mitochondria-encoded open reading frames (ORFs), 85.41% of these editing sites were observed on the first or the second base of a codon, resulting in the alteration of encoded amino acid. Moreover, we found some novel editing sites and several inaccurately annotated sites which may be functionally important, based on the highly conserved amino acids encoded by these edited codons. Finally, we annotated all 569 C-to-U RNA editing sites in their biological context. More precise information about C-to-U editing sites in rice mitochondria-encoded ORFs will facilitate our investigation on the function of C-to-U editing events in rice and also provide a valid benchmark from rice for the analysis of mitochondria C-to-U editing in other plant species.

The Genome and Transcriptome Analysis of the Vigna mungo Chloroplast

Vigna mungo is cultivated in approximately 5 million hectares worldwide. The chloroplast genome of this species has not been previously reported. In this study, we sequenced the genome and transcriptome of the V. mungo chloroplast. We identified many positively selected genes in the photosynthetic pathway (e.g., rbcL, ndhF, and atpF) and RNA polymerase genes (e.g., rpoC2) from the comparison of the chloroplast genome of V. mungo, temperate legume species, and tropical legume species. Our transcriptome data from PacBio isoform sequencing showed that the 51-kb DNA inversion could affect the transcriptional regulation of accD polycistronic. Using Illumina deep RNA sequencing, we found RNA editing of clpP in the leaf, shoot, flower, fruit, and root tissues of V. mungo. We also found three G-to-A RNA editing events that change guanine to adenine in the transcripts transcribed from the adenine-rich regions of the ycf4 gene. The edited guanine bases were found particularly in the chloroplast genome of the Vigna species. These G-to-A RNA editing events were likely to provide a mechanism for correcting DNA base mutations. The V. mungo chloroplast genome sequence and the analysis results obtained in this study can apply to phylogenetic studies and chloroplast genome engineering.

Systematic analysis reveals cis and trans determinants affecting C-to-U RNA editing in Arabidopsis thaliana

Background

C-to-U RNA editing is prevalent in the mitochondrial and chloroplast genes in plants. The biological functions of a fraction of C-to-U editing sites are continuously discovered by case studies. However, at genome-wide level, the cis and trans determinants affecting the occurrence or editing levels of these C-to-U events are relatively less studied. What is known is that the PPR (pentatricopeptide repeat) proteins are the main trans-regulatory elements responsible for the C-to-U conversion, but other determinants especially the cis-regulatory elements remain largely uninvestigated.

Results

By analyzing the transcriptome and translatome data in Arabidopsis thaliana roots and shoots, combined with RNA-seq data from hybrids of Arabidopsis thaliana and Arabidopsis lyrata, we perform genome-wide investigation on the cis elements and trans-regulatory elements that potentially affect C-to-U editing events. An upstream guanosine or double-stranded RNA (dsRNA) regions are unfavorable for editing events. Meanwhile, many genes including the transcription factors may indirectly play regulatory roles in trans.

Conclusions

The 5-prime thymidine facilitates editing and dsRNA structures prevent editing in cis. Many transcription factors affect editing in trans. Although the detailed molecular mechanisms underlying the cis and trans regulation remain to be experimentally verified, our findings provide novel aspects in studying the botanical C-to-U RNA editing events.

Mutations in m6A consensus motifs are suppressed in the m6A modified genes in human cancer cells

N⁶-methyladenosine (m⁶A) is the most prevalent type of RNA modification. METTL3 in the methyltransferase complex is the core enzyme responsible for methylation. METTL3 selectively catalyzes the adenosines centered in the RRAC motif. Functional studies established that m⁶A could enhance the translation efficiency (TE) of modified genes by recruiting reader protein YTHDF1 and other initiation factors. We downloaded the m⁶A peaks in HeLa cells from a previous study and defined the m⁶A modified genes and sites. Ancestral mutations in the genic region fixed in the HeLa cell samples were defined using their mRNA-Seq data and the alignment between human and mouse genomes. Furthermore, in the small interfering (si)-METTL3 sample, the calculated TE foldchange of all genes was compared to that in the negative control. The TE of m⁶A genes was globally down-regulated in si-METTL3 versus control compared to the non-m⁶A genes. In m⁶A modified genes, RRAC motif mutations were suppressed compared to mutations in non-motif regions or non-m⁶A genes. Among the m⁶A genes, a fraction RRAC motif mutations negatively correlated with the TE foldchange (si-METTL3 versus control). The TE of m⁶A modified genes was enhanced in HeLa cells. RRAC motif mutations could potentially prevent methylation of adenosines and consequently abolish the enhanced translation. Such mutations in the RRAC motif might be deleterious. Accordingly, we observed lower fractions of mutations in RRAC motifs than in other regions. This prevention of mutations in the RRAC motif could be a strategy adopted by cancer cells to maintain the elevated translation of particular genes.

Disruption of a Rice Chloroplast-Targeted Gene OsHMBPP Causes a Seedling-Lethal Albino Phenotype

BACKGROUND

Chloroplast development is coordinately regulated by plastid- and nuclear-encoding genes. Although many regulators have been reported to be involved in chloroplast development, new factors remain to be identified, given the complexity of this process.

RESULTS

In this study, we characterized a rice mutant lethal albinic seedling 1(las1)form of a 4-hydroxy-3-methylbut-2-enyl diphosphate reductase (OsHMBPP) that was targeted to the chloroplasts. The LAS1 mutation caused the albino lethal phenotype in seedlings. Transmission electron microscopy indicated that las1 were defective in early chloroplast development. LAS1 is preferentially expressed in leaves, implying its role in controlling chloroplast development. The expression levels of many chloroplast-encoded genes were altered significantly in las1. The expression levels of nuclear-encoded gene involved in Chl biosynthesis were also decreased in las1. We further investigated plastidic RNA editing in las1 and found that the edit efficiency of four chloroplast genes were markly altered. Compared with WT, las1 exhibited defective in biogenesis of chloroplast ribosomes.

CONCLUSIONS

Our results show that LAS1/OsHMBPP plays an essential role in the early chloroplast development in rice.

Salinity effects on nad3 gene RNA editing of wild barley mitochondria

Nad complex plays a very important role during cellular respiration. nad3 (nad dehydrogenase subunit 3) is one of the biggest subunits in this complex. Four cDNAs of nad3 gene were characterized in Hordeum vulgare subsp. spontaneum at exposed to four periods of 500 mM salinity, 0 h or control (accession no. MN066165), after 2 h (accession no. MN066166), after 12 h (accession no. MN066167) and after 24 h (accession no. MN066168) using RNA-seq raw data. Seventeen RNA editing sites were found in positions (or nucleotide nos. C5, C39, C44, C61, C62, C79, C80, C147, C185, C190, C191, C208, C209, C275, C317, C344, C349) within the nad3 coding region. These alterations represent differential editing at four exposure times. The maximum editing rate was revealed 2 and 12 h after salinity exposure. However, these edits were disrupted after 24 h probably due to the initiation of program cell death machinery. We found that RNA editing not only improved protein function but also may improve codon bias by altering the nucleotide without any change in amino acid. Characterization of pentatricopeptide repeat-containing protein At4g13650 (PPRSp1) in wild barley helped us to understand the behavior of editing sites C190 and C191 under salinity. Position - 6 in cis-element upstream editing sites of C155, C190 and C191 may be vital to the editing process in these sites by PPRSp1 protein. The differential editing of this gene under salinity led to a relationship between RNA editing and cellular respiration regulation.

Genome-wide analysis on the maize genome reveals weak selection on synonymous mutations

Background

Synonymous mutations are able to change the tAI (tRNA adaptation index) of a codon and consequently affect the local translation rate. Intuitively, one may hypothesize that those synonymous mutations which increase the tAI values are favored by natural selection.

Results

We use the maize (Zea mays) genome to test our assumption. The first supporting evidence is that the tAI-increasing synonymous mutations have higher fixed-to-polymorphic ratios than the tAI-decreasing ones. Next, the DAF (derived allele frequency) or MAF (minor allele frequency) of the former is significantly higher than the latter. Moreover, similar results are obtained when we investigate CAI (codon adaptation index) instead of tAI.

Conclusion

The synonymous mutations in the maize genome are not strictly neutral. The tAI-increasing mutations are positively selected while those tAI-decreasing ones undergo purifying selection. This selection force might be weak but should not be automatically ignored.

Reduced C-to-U RNA editing rates might play a regulatory role in stress response of Arabidopsis

C-to-U RNA editing is prevalent in the mitochondrial and chloroplast genes in plants. The C-to-U editing rates are constantly very high. During genome evolution, those edited cytidines are likely to be replaced with thymidines at the DNA level. C-to-U editing events are suggested to be designed for reversing the unfavorable T-to-C DNA mutations. Despite the existing theory showing the importance of editing mechanisms, few studies have investigated the genome-wide adaptive signals of the C-to-U editome or the potential function of C-to-U editing events in the stress response. By analyzing the transcriptome and translatome data of normal and heat-shocked Arabidopsis thaliana and the RNA-seq from cold-stressed plants, combined with genome-wide comparison of mitochondrial/chloroplast genes and nuclear genes from multiple aspects, we present the conservational and translational features of each gene and depict the dynamic mitochondrial/chloroplast C-to-U RNA editome. We found that the tAI (tRNA adaptation index) and basic translation levels are lower for mitochondrial/chloroplast genes than for nuclear genes. Interestingly, although we found adaptive signals for the global C-to-U RNA editome in mitochondrial/chloroplast genes, the C-to-U (T) alteration would usually cause a reduction in the codon tAI value. Moreover, the C-to-U editing rates are significantly reduced under heat or cold stress when compared to the normal condition. This reduction is irrelevant to the temperature-sensitive RNA structures. Several cases have illustrated that under heat stress, the reduced C-to-U editing rates alleviate ribosome stalling and consequently facilitate the local translation. Our study reveals that in Arabidopsis thaliana the mitochondrial/chloroplast C-to-U RNA editing rates are reduced under heat or cold stress. This reduction is associated with the alleviation of decreased tAI/translation rate of edited codons. The regulation of C-to-U editing rates could be the tradeoff between quantity and quality. We profile the dynamic change of C-to-U RNA editome under heat stress and propose a potential role of editing sites in the heat response. Our work should be appealing to the plant physiologists as well as the RNA editing community.