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

PlantC2U: deep learning of cross-species sequence landscapes predicts plastid C-to-U RNA editing in plants

In plants, C-to-U RNA editing mainly occurs in plastid and mitochondrial transcripts, which contributes to a complex transcriptional regulatory network. More evidence reveals that RNA editing plays critical roles in plant growth and development. However, accurate detection of RNA editing sites using transcriptome sequencing data alone is still challenging. In the present study, we develop PlantC2U, which is a convolutional neural network, to predict plastid C-to-U RNA editing based on the genomic sequence. PlantC2U achieves >95% sensitivity and 99% specificity, which outperforms the PREPACT tool, random forests, and support vector machines. PlantC2U not only further checks RNA editing sites from transcriptome data to reduce possible false positives, but also assesses the effect of different mutations on C-to-U RNA editing based on the flanking sequences. Moreover, we found the patterns of tissue-specific RNA editing in the mangrove plant Kandelia obovata, and observed reduced C-to-U RNA editing rates in the cold stress response of K. obovata, suggesting their potential regulatory roles in plant stress adaptation. In addition, we present RNAeditDB, available online at https://jasonxu.shinyapps.io/RNAeditDB/. Together, PlantC2U and RNAeditDB will help researchers explore the RNA editing events in plants and thus will be of broad utility for the plant research community.

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.

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.

C-to-U RNA deamination is the driving force accelerating SARS-CoV-2 evolution

Understanding the molecular mechanism underlying the rampant mutation of SARS-CoV-2 would help us control the COVID-19 pandemic. The APOBEC-mediated C-to-U deamination is a major mutation type in the SARS-CoV-2 genome. However, it is unclear whether the novel mutation rate u is higher for C-to-U than for other mutation types, and what the detailed driving force is. By analyzing the time course SARS-CoV-2 global population data, we found that C-to-U has the highest novel mutation rate u among all mutation types and that this u is still increasing with time (du/dt > 0). Novel C-to-U events, rather than other mutation types, have a preference over particular genomic regions. A less local RNA structure is correlated with a high novel C-to-U mutation rate. A cascade model nicely explains the du/dt > 0 for C-to-U deamination. In SARS-CoV-2, the RNA structure serves as the molecular basis of the extremely high and continuously accelerating C-to-U deamination rate. This mechanism is the driving force of the mutation, adaptation, and evolution of SARS-CoV-2. Our findings help us understand the dynamic evolution of the virus mutation rate.

Exploring the RNA Editing Events and Their Potential Regulatory Roles in Tea Plant (Camellia sinensis L.)

RNA editing is a post-transcriptional modification process that alters the RNA sequence relative to the genomic blueprint. In plant organelles (namely, mitochondria and chloroplasts), the most common type is C-to-U, and the absence of C-to-U RNA editing results in abnormal plant development, such as etiolation and albino leaves, aborted embryonic development and retarded seedling growth. Here, through PREP, RES-Scanner, PCR and RT-PCR analyses, 38 and 139 RNA editing sites were identified from the chloroplast and mitochondrial genomes of Camellia sinensis, respectively. Analysis of the base preference around the RNA editing sites showed that in the -1 position of the edited C had more frequent occurrences of T whereas rare occurrences of G. Three conserved motifs were identified at 25 bases upstream of the RNA editing site. Structural analyses indicated that the RNA secondary structure of 32 genes, protein secondary structure of 37 genes and the three-dimensional structure of 5 proteins were altered due to RNA editing. The editing level analysis of matK and ndhD in six tea cultivars indicated that matK-701 might be involved in the color change of tea leaves. Furthermore, 218 PLS-CsPPR proteins were predicted to interact with the identified RNA editing sites. In conclusion, this study provides comprehensive insight into RNA editing events, which will facilitate further study of the RNA editing phenomenon of the tea plant.

COVID-19 and its genomic variants: Molecular pathogenesis and therapeutic interventions

Coronavirus disease-19 (COVID-19), caused by a β-coronavirus and its genomic variants, is associated with substantial morbidities and mortalities globally. The COVID-19 virus and its genomic variants enter host cells upon binding to the angiotensin converting enzyme 2 receptors that are expressed in a variety of tissues, but predominantly in the lungs, heart, and blood vessels. Patients afflicted with COVID-19 may be asymptomatic or present with critical symptoms possibly due to diverse lifestyles, immune responses, aging, and underlying medical conditions. Geriatric populations, especially men in comparison to women, with immunocompromised conditions, are most vulnerable to severe COVID-19 associated infections, complications, and mortalities. Notably, whereas immunomodulation, involving nutritional consumption, is essential to protecting an individual from COVID-19, immunosuppression is detrimental to a person with this aggressive disease. As such, immune health is inversely correlated to COVID-19 severity and resulting consequences. Advances in genomic and proteomic technologies have helped us to understand the molecular events underlying symptomatology, transmission and, pathogenesis of COVID-19 and its genomic variants. Accordingly, there has been development of a variety of therapeutic interventions, ranging from mask wearing to vaccination to medication. This review summarizes the current understanding of molecular pathogenesis of COVID-19, effects of comorbidities on COVID-19, and prospective therapeutic strategies for the prevention and treatment of this contagious disease.

Gene expression variation explains maize seed germination heterosis

BACKGROUND

Heterosis has been extensively utilized in plant breeding, however, the underlying molecular mechanism remains largely elusive. Maize (Zea mays), which exhibits strong heterosis, is an ideal material for studying heterosis.

RESULTS

In this study, there is faster imbibition and development in reciprocal crossing Zhengdan958 hybrids than in their parent lines during seed germination. To investigate the mechanism of heterosis of maize germination, comparative transcriptomic analyses were conducted. The gene expression patterns showed that 1324 (47.27%) and 1592 (66.44%) of the differential expression genes between hybrids and either parental line display parental dominance up or higher levels in the reciprocal cross of Zhengdan958, respectively. Notably, these genes were mainly enriched in metabolic pathways, including carbon metabolism, glycolysis/gluconeogenesis, protein processing in endoplasmic reticulum, etc. CONCLUSION: Our results provide evidence for the higher expression level genes in hybrid involved in metabolic pathways acting as main contributors to maize seed germinating heterosis. These findings provide new insights into the gene expression variation of maize embryos and improve the understanding of maize seed germination heterosis.

Large-scale omics data reveal the cooperation of mutation-circRNA-miRNA-target gene network in liver cancer oncogenesis

Aims: Clarifying the initial trigger of the differentially expressed genes in cancers helps researchers understand the cellular system as a whole network. Materials & methods: We retrieve the transcriptome and translatome of tumor and normal tissues from ten liver cancer patients and define differentially expressed genes and tumor-specific mutations. We associate the oncogenesis with the mutations by target prediction and experimental verification. Results: Upregulated genes have tumor-specific mutations in 3'UTRs that abolish the binding of miRNAs. For downregulated genes, their corresponding miRNAs are mutually targeted by two circRNAs, with mutations in base-pairing regions. Transfection experiments support the oncogenic role of these mutations. Conclusions: The tumor-specific mutations serve as the initial trigger of liver cancer. The mutation-circRNA-miRNA-target gene chain is completed.

Direct in vivo observation of the effect of codon usage bias on gene expression in Arabidopsis hybrids

Hybrids are the perfect materials to study cis regulatory elements because the two parental alleles are subjected to identical trans environments. There has been a debate on whether synonymous codon usage could affect gene expression. In vitro experiments found that luciferase genes with enhanced codon optimality showed elevated mRNA expression. However, the underlying mechanism is still unclear, and no direct evidence is observed to support this notion. By mapping the RNA-seq data of hybrids of Arabidopsis thaliana and Arabidopsis lyrata, we quantified the allele-specific reads and estimated the relative expression of orthologous genes. We focused on orthologous genes with dN = 0 and dS > 0, which means that they only differ in synonymous codon usage. We found that orthologous genes with higher codon optimality in A. thaliana tend to have higher expression levels of the A. thaliana allele. Codon usage bias could influence gene expression. This phenomenon is not only found in in vitro experiments but also supported by in vivo observations. Therefore, synonymous mutations could have a broad impact from multiple aspects and should not be automatically ignored in genomic studies. KEY MESSAGE: In Arabidopsis hybrids, alleles with higher codon optimality tend to have higher expression levels.

Mutation profiling of a limbless pig reveals genome-wide regulation of RNA processing related to bone development

Mutation is the basis of phenotypic changes and serves as the source of natural selection. The development of limbs has been the milestone in vertebrate evolution. Several limb and bone-related genes were verified experimentally, but other indirect and regulatory factors of limb development remained untested, especially very few cases were observed in natural environment. We report a naturally born serpentized pig without hindlimbs. Whole genome sequencing followed by comparative genomic analysis revealed multiple interesting patterns on the handicapped pig-specific mutations. Although the bone-related genes are not directly subjected to mutations, other regulatory factors such as the RNA deaminase genes Adar are damaged in the handicapped pig, leading to the abolished A-to-I deamination in many functional, conserved genes as well as the bone-related genes. This is a precious case that the limbless phenotype is observed in naturally born non-model organisms. Our study broadened the generality of the limbless phenotype across mammals and extended the regulation of hindlimb development to other non-bone-related genes. Our knowledge of limb and bone-related mutations and regulation would also contribute to human genetics.

Context-dependent and -independent selection on synonymous mutations revealed by 1,135 genomes of Arabidopsis thaliana

BACKGROUND

Synonymous mutations do not alter the amino acids and therefore are regarded as neutral for a long time. However, they do change the tRNA adaptation index (tAI) of a particular codon (independent of its context), affecting the tRNA availability during translation. They could also change the isoaccepting relationship with its neighboring synonymous codons in particular context, which again affects the local translation process. Evidence of selection pressure on synonymous mutations has emerged.

RESULTS

The proposed selection patterns on synonymous mutations are never formally and systematically tested in plant species. We fully take advantage of the SNP data from 1,135 A. thaliana lines, and found that the synonymous mutations that increase tAI or the isoaccepting mutations in isoaccepting codon context tend to have higher derived allele frequencies (DAF) compared to other synonymous mutations of the opposite effects.

CONCLUSIONS

Synonymous mutations are not strictly neutral. The synonymous mutations that increase tAI or the isoaccepting mutations in isoaccepting codon context are likely to be positively selected. We propose the concept of context-dependent and -independent selection on synonymous mutations. These concepts broaden our knowledge of the functional consequences of synonymous mutations, and should be appealing to phytologists and evolutionary biologists.

Maize kernel development

Maize (Zea mays) is a leading cereal crop in the world. The maize kernel is the storage organ and the harvest portion of this crop and is closely related to its yield and quality. The development of maize kernel is initiated by the double fertilization event, leading to the formation of a diploid embryo and a triploid endosperm. The embryo and endosperm are then undergone independent developmental programs, resulting in a mature maize kernel which is comprised of a persistent endosperm, a large embryo, and a maternal pericarp. Due to the well-characterized morphogenesis and powerful genetics, maize kernel has long been an excellent model for the study of cereal kernel development. In recent years, with the release of the maize reference genome and the development of new genomic technologies, there has been an explosive expansion of new knowledge for maize kernel development. In this review, we overviewed recent progress in the study of maize kernel development, with an emphasis on genetic mapping of kernel traits, transcriptome analysis during kernel development, functional gene cloning of kernel mutants, and genetic engineering of kernel traits.