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

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.

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.

Hepatitis B virus (HBV) codon adapts well to the gene expression profile of liver cancer: an evolutionary explanation for HBV's oncogenic role

Due to the evolutionary arms race between hosts and viruses, viruses must adapt to host translation systems to rapidly synthesize viral proteins. Highly expressed genes in hosts have a codon bias related to tRNA abundance, the primary RNA translation rate determinant. We calculated the relative synonymous codon usage (RSCU) of three hepatitis viruses (HAV, HBV, and HCV), SARS-CoV-2, 30 human tissues, and hepatocellular carcinoma (HCC). After comparing RSCU between viruses and human tissues, we calculated the codon adaptation index (CAI) of viral and human genes. HBV and HCV showed the highest correlations with HCC and the normal liver, while SARS-CoV-2 had the strongest association with lungs. In addition, based on HCC RSCU, the CAI of HBV and HCV genes was the highest. HBV and HCV preferentially adapt to the tRNA pool in HCC, facilitating viral RNA translation. After an initial trigger, rapid HBV/HCV translation and replication may change normal liver cells into HCC cells. Our findings reveal a novel perspective on virus-mediated oncogenesis.

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.

SARS-CoV-2 competes with host mRNAs for efficient translation by maintaining the mutations favorable for translation initiation

During SARS-CoV-2 proliferation, the translation of viral RNAs is usually the rate-limiting step. Understanding the molecular details of this step is beneficial for uncovering the origin and evolution of SARS-CoV-2 and even for controlling the pandemic. To date, it is unclear how SARS-CoV-2 competes with host mRNAs for ribosome binding and efficient translation. We retrieved the coding sequences of all human genes and SARS-CoV-2 genes. We systematically profiled the GC content and folding energy of each CDS. Considering that some fixed or polymorphic mutations exist in SARS-CoV-2 and human genomes, all algorithms and analyses were applied to both pre-mutate and post-mutate versions. In SARS-CoV-2 but not human, the 5-prime end of CDS had lower GC content and less RNA structure than the 3-prime part, which was favorable for ribosome binding and efficient translation initiation. Globally, the fixed and polymorphic mutations in SARS-CoV-2 had created an even lower GC content at the 5-prime end of CDS. In contrast, no similar patterns were observed for the fixed and polymorphic mutations in human genome. Compared with human RNAs, the SARS-CoV-2 RNAs have less RNA structure in the 5-prime end and thus are more favorable of fast translation initiation. The fixed and polymorphic mutations in SARS-CoV-2 are further amplifying this advantage. This might serve as a strategy for SARS-CoV-2 to adapt to the human host.

Compelling Evidence Suggesting the Codon Usage of SARS-CoV-2 Adapts to Human After the Split From RaTG13

SARS-CoV-2 needs to efficiently make use of the resources from hosts in order to survive and propagate. Among the multiple layers of regulatory network, mRNA translation is the rate-limiting step in gene expression. Synonymous codon usage usually conforms with tRNA concentration to allow fast decoding during translation. It is acknowledged that SARS-CoV-2 has adapted to the codon usage of human lungs so that the virus could rapidly proliferate in the lung environment. While this notion seems to nicely explain the adaptation of SARS-CoV-2 to lungs, it is unable to tell why other viruses do not have this advantage. In this study, we retrieve the GTEx RNA-seq data for 30 tissues (belonging to over 17 000 individuals). We calculate the RSCU (relative synonymous codon usage) weighted by gene expression in each human sample, and investigate the correlation of RSCU between the human tissues and SARS-CoV-2 or RaTG13 (the closest coronavirus to SARS-CoV-2). Lung has the highest correlation of RSCU to SARS-CoV-2 among all tissues, suggesting that the lung environment is generally suitable for SARS-CoV-2. Interestingly, for most tissues, SARS-CoV-2 has higher correlations with the human samples compared with the RaTG13-human correlation. This difference is most significant for lungs. In conclusion, the codon usage of SARS-CoV-2 has adapted to human lungs to allow fast decoding and translation. This adaptation probably took place after SARS-CoV-2 split from RaTG13 because RaTG13 is less perfectly correlated with human. This finding depicts the trajectory of adaptive evolution from ancestral sequence to SARS-CoV-2, and also well explains why SARS-CoV-2 rather than other viruses could perfectly adapt to human lung environment.

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.

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.

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.