The divergence between SARS-CoV-2 and RaTG13 might be overestimated due to the extensive RNA modification

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.

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.

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.

Retrospect of the Two-Year Debate: What Fuels the Evolution of SARS-CoV-2: RNA Editing or Replication Error?

Mutation is one of the mechanisms of the evolutionary divergence of an organism. Under this global COVID-19 pandemic, the fast evolution of SARS-CoV-2 became one of the most worrying issues. Some researchers believed that the hosts' RNA deamination systems (APOBECs and ADARs) are the major source of mutations and have driven the evolution of SARS-CoV-2. However, apart from RNA editing, the RDRP (RNA-dependent RNA polymerase)-mediated replication errors may also contribute to the mutation of SARS-CoV-2 (just like the single-nucleotide polymorphisms/variations in eukaryotes caused by DNA replication errors). Unfortunately, it is technically unable to distinguish RNA editing and replication errors (SNPs) in this RNA virus. Here comes a fundamental question: we indeed observed the fast evolution of SARS-CoV-2, but what exactly fuels its evolution: RNA editing or replication errors? This debate lasts for 2 years. In this piece, we will retrospect the 2-year debate on RNA editing versus SNPs.

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.

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.

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.

The Sponge Interaction Between Circular RNA and microRNA Serves as a Fast-Evolving Mechanism That Suppresses Non-small Cell Lung Cancer (NSCLC) in Humans

Non-small cell lung cancer (NSCLC) is one of the most lethal cancer types in the world. Currently, the molecular mechanisms and pathways underlying NSCLC oncogenesis are poorly understood. Using multiple Omics data, we systematically explored the differentially expressed circular RNAs (circRNAs) in NSCLC. We also investigated potential microRNA sponges (that absorb circRNAs) in NSCLC and downstream target genes with experimental verifications. hsa_circ_0003497 was down-regulated in NSCLC and played an inhibitory role in tumorigenesis. In contrast, miR-197-3p was up-regulated in NSCLC. hsa_circ_0003497 directly interacts with miR-197-3p and releases a target gene of miR-197-3p termed CTNND1 (a known tumor suppressor gene). Evolutionary analysis reveals fast evolution of this hsa_circ_0003497-miR-197-3p-CTNND1-NSCLC axis in mammals. This work clarified the biological functions and molecular mechanisms of how hsa_circ_0003497 suppresses NSCLC through miR-197-3p and CTNND1. We discovered molecular markers for the prognosis of NSCLC and provided potential intervention targets for its treatment.

uORF-Mediated Translational Regulation of ATF4 Serves as an Evolutionarily Conserved Mechanism Contributing to Non-Small-Cell Lung Cancer (NSCLC) and Stress Response

Diseases and environmental stresses are two distinct challenges for virtually all living organisms. In light of evolution, cellular responses to diseases and stresses might share similar molecular mechanisms, but the detailed regulation pathway is not reported yet.

We obtained the transcriptomes and translatomes from several NSCLC (non-small-cell lung cancer) patients as well as from different species under normal or stress conditions. We found that the translation level of gene ATF4 is remarkably enhanced in NSCLC due to the reduced number of ribosomes binding to its upstream open reading frames (uORFs). We also showed the evolutionary conservation of this uORF-ATF4 regulation in the stress response of other species. Molecular experiments showed that knockdown of ATF4 reduced the cell growth rate while overexpression of ATF4 enhanced cell growth, especially for the ATF4 allele with mutated uORFs. Population genetics analyses in multiple species verified that the mutations that abolish uATGs (start codon of uORFs) are highly deleterious, suggesting the functional importance of uORFs.

Our study proposes an evolutionarily conserved pattern that enhances the ATF4 translation by uORFs upon stress or disease. We generalized the concept of cellular response to diseases and stresses. These two biological processes may share similar molecular mechanisms.

Nothing in SARS-CoV-2 makes sense except in the light of RNA modification?

The expression pattern of RNA deaminases determines the mutation and evolution of SARS-CoV-2.

Evidence for selection on SARS-CoV-2 RNA translation revealed by the evolutionary dynamics of mutations in UTRs and CDSs

RNA translation is the rate-limiting step when cells synthesize proteins. Elevating translation efficiency (TE) is intuitively beneficial. Particularly, when viruses invade host cells, how to compete with endogenous RNAs for efficient translation is a major issue to be resolved. We collected millions of worldwide SARS-CoV-2 sequences during the past year and traced the dynamics of allele frequency of every mutation. We defined adaptive and deleterious mutations according to the rise and fall of their frequencies along time. For 5ʹUTR and synonymous mutations in SARS-CoV-2, the selection on TE is evident near start codons. Adaptive mutations generally decrease GC content while deleterious mutations increase GC content. This trend fades away with increasing distance to start codons. Mutations decreasing GC content near start codons would unravel the complex RNA structure and facilitate translation initiation, which are beneficial to SARS-CoV-2, and vice versa. During this evolutionary arms race between human and virus, SARS-CoV-2 tries to improve its cis elements to compete with host RNAs for rapid translation.

Inverted repeats in coronavirus SARS-CoV-2 genome manifest the evolution events

The world faces a great unforeseen challenge through the COVID-19 pandemic caused by coronavirus SARS-CoV-2. The virus genome structure and evolution are positioned front and center for further understanding insights on vaccine development, monitoring of transmission trajectories, and prevention of zoonotic infections of new coronaviruses. Of particular interest are genomic elements Inverse Repeats (IRs), which maintain genome stability, regulate gene expressions, and are the targets of mutations. However, little research attention is given to the IR content analysis in the SARS-CoV-2 genome. In this study, we propose a geometric analysis method and using the method to investigate the distributions of IRs in SARS-CoV-2 and its related coronavirus genomes. The method represents each genomic IR sequence pair as a single point and constructs the geometric shape of the genome using the IRs. Thus, the IR shape can be considered as the signature of the genome. The genomes of different coronaviruses are then compared using the constructed IR shapes. The results demonstrate that SARS-CoV-2 genome, specifically, has an abundance of IRs, and the IRs in coronavirus genomes show an increase during evolution events.

Fast evolution of SARS-CoV-2 driven by deamination systems in hosts

As an RNA virus, the fast evolution of SARS-CoV-2 is driven by the extensive RNA deamination by the host cells.

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.

Natural selection on synonymous mutations in SARS-CoV-2 and the impact on estimating divergence time

To adapt to human host environment, synonymous mutations in SARS-CoV-2 are shaped by tRNA selection, energy cost and RNA structure.

The percentages of SARS-CoV-2 protein similarity and identity with SARS-CoV and BatCoV RaTG13 proteins can be used as indicators of virus origin

There are three types of proteins in coronaviruses: nonstructural, structural, and accessory proteins. Coronavirus proteins are essential for viral replication and for the binding and invasion of hosts and the regulation of host cell metabolism and immunity. This study investigated the amino acid sequence similarity and identity percentages of 10 proteins in SARS-CoV-2, SARS-CoV and the Rhinolophus affinis bat coronavirus (BatCoV RaTG13). The investigated proteins were the 1ab polyprotein, spike protein, orf3a, the envelope protein, the membrane protein, orf6, orf7a, orf7b, orf8, and the nucleocapsid protein. The online sequence alignment service of The European Molecular Biology Open Software Suite (EMBOSS) was used to determine the percentages of protein similarity and identity in the three viruses. The results showed that the similarity and identity percentages of the SARS-CoV-2 and BatCoV RaTG13 proteins were both greater than 95%, while the identity and similarity percentages of SARS-CoV-2 and SARS-CoV were both greater than 38%. The proteins of SARS-CoV-2 and BatCoV RaTG13 have high identity and similarity compared to those of SARS-CoV-2 and SARS-CoV.

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.

Clade GR and clade GH isolates of SARS-CoV-2 in Asia show highest amount of SNPs

Clades are monophyletic groups composed of a common ancestor and all its lineal descendants. As the propensity of virulence of a disease depends upon the type of clade the virus belongs to and it causes different fatality rates of disease in different countries, so the clade-wise analysis of SARS-CoV-2 isolates collected from different countries can illuminate the actual evolutionary relationships between them. In this study, 1566 SARS-CoV-2 genome sequences across ten Asian countries are collected, clustered, and characterized based on the clade they belong to. The isolates are compared to the Wuhan reference sequence” hCoV-19/Wuhan/WIV04/19″ to identify the mutations that occurred at different protein regions. Structural changes in amino acids due to mutations lead to functional instability of the proteins. Detailed clade-wise functional assessments are carried out to quantify the stability and vulnerability of the mutations occurring in SARS-CoV-2 genomes which can shade light on personalized prevention and treatment of the disease and encourage towards the invention of clade-specific vaccines.

A Comprehensive Analysis of cis-Acting RNA Elements in the SARS-CoV-2 Genome by a Bioinformatics Approach

The emergence of a new coronavirus (CoV), severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for severe respiratory disease in humans termed coronavirus disease of 2019 (COVID-19), became a new global threat for health and the economy. The SARS-CoV-2 genome is about a 29,800-nucleotide-long plus-strand RNA that can form functionally important secondary and higher-order structures called cis-acting RNA elements. These elements can interact with viral proteins, host proteins, or other RNAs and be involved in regulating translation and replication processes of the viral genome and encapsidation of the virus. However, the cis-acting RNA elements and their biological roles in SARS-CoV-2 as well as their comparative analysis in the closely related viral genome have not been well explored, which is very important to understand the molecular mechanism of viral infection and pathogenies. In this study, we used a bioinformatics approach to identify the cis-acting RNA elements in the SARS-CoV-2 genome. Initially, we aligned the full genomic sequence of six different CoVs, and a phylogenetic analysis was performed to understand their evolutionary relationship. Next, we predicted the cis-acting RNA elements in the SARS-CoV-2 genome using the structRNAfinder tool. Then, we annotated the location of these cis-acting RNA elements in different genomic regions of SARS-CoV-2. After that, we analyzed the sequence conservation patterns of each cis-acting RNA element among the six CoVs. Finally, the presence of cis-acting RNA elements across different CoV genomes and their comparative analysis was performed. Our study identified 12 important cis-acting RNA elements in the SARS-CoV-2 genome; among them, Corona_FSE, Corona_pk3, and s2m are highly conserved across most of the studied CoVs, and Thr_leader, MAT2A_D, and MS2 are uniquely present in SARS-CoV-2. These RNA structure elements can be involved in viral translation, replication, and encapsidation and, therefore, can be potential targets for better treatment of COVID-19. It is imperative to further characterize these cis-acting RNA elements experimentally for a better mechanistic understanding of SARS-CoV-2 infection and therapeutic intervention.

Mutational insights into the envelope protein of SARS-CoV-2

The ongoing mutations in the structural proteins of SARS-CoV-2 are the major impediment for prevention and control of the COVID-19 disease. Presently we focused on evolution of the envelope (E) protein, one of the most enigmatic and less studied protein among the four structural proteins (S, E, M and N) associated with multitude of immunopathological functions of SARS-CoV-2. In the present study, we comprehensively analyzed 81,818 high quality E protein sequences of SARS-CoV-2 globally available in the GISAID database as of 20 August 2020. Compared to Wuhan reference strain, our mutational analysis explored only 1.2 % (982/81818) mutant strains undergoing a total of 115 unique amino acid (aa) substitutions in the E protein, highlighting the fact that most (98.8 %) of the E protein of SARS-CoV-2 strains are highly conserved. Moreover, we found 58.77 % (134 of 228) nucleotides (nt) positions of SARS-CoV-2 E gene encountering a total of 176 unique nt-level mutations globally, which may affect the efficacy of real time RT-PCR-based molecular detection of COVID-19. Importantly, higher aa variations observed in the C-terminal domain (CTD) of the E protein, particularly at Ser⁵⁵-Phe⁵⁶, Arg⁶⁹ and the C-terminal end (DLLV: 72–75) may alter the binding of SARS-CoV-2 Envelope protein to tight junction-associated PALS1 and thus could play a key role in COVID-19 pathogenesis. Furthermore, this study revealed the V25A mutation in the transmembrane domain which is a key factor for the homopentameric conformation of E protein. Our analysis also observed a triple cysteine motif harboring mutation (L39M, A41S, A41V, C43F, C43R, C43S, C44Y, N45R) which may hinder the binding of E protein with spike glycoprotein. These results therefore suggest the continuous monitoring of the structural proteins including the envelope protein of SARS-CoV-2 since the number of genome sequences from across the world are continuously increasing.

Positive selection within the genomes of SARS-CoV-2 and other Coronaviruses independent of impact on protein function

BACKGROUND

The emergence of a novel coronavirus (SARS-CoV-2) associated with severe acute respiratory disease (COVID-19) has prompted efforts to understand the genetic basis for its unique characteristics and its jump from non-primate hosts to humans. Tests for positive selection can identify apparently nonrandom patterns of mutation accumulation within genomes, highlighting regions where molecular function may have changed during the origin of a species. Several recent studies of the SARS-CoV-2 genome have identified signals of conservation and positive selection within the gene encoding Spike protein based on the ratio of synonymous to nonsynonymous substitution. Such tests cannot, however, detect changes in the function of RNA molecules.

METHODS

Here we apply a test for branch-specific oversubstitution of mutations within narrow windows of the genome without reference to the genetic code.

RESULTS

We recapitulate the finding that the gene encoding Spike protein has been a target of both purifying and positive selection. In addition, we find other likely targets of positive selection within the genome of SARS-CoV-2, specifically within the genes encoding Nsp4 and Nsp16. Homology-directed modeling indicates no change in either Nsp4 or Nsp16 protein structure relative to the most recent common ancestor. These SARS-CoV-2-specific mutations may affect molecular processes mediated by the positive or negative RNA molecules, including transcription, translation, RNA stability, and evasion of the host innate immune system. Our results highlight the importance of considering mutations in viral genomes not only from the perspective of their impact on protein structure, but also how they may impact other molecular processes critical to the viral life cycle.

GC usage of SARS-CoV-2 genes might adapt to the environment of human lung expressed genes

Understanding how SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus 2) efficiently reproduces itself by taking resources from the human host could facilitate the development of drugs against the virus. SARS-CoV-2 translates its own proteins by using the host tRNAs, so that its GC or codon usage should fit that of the host cells. It is necessary to study both the virus and human genomes in the light of evolution and adaptation. The SARS-CoV-2 virus has significantly lower GC content and GC3 as compared to human. However, when we selected a set of human genes that have similar GC properties to SARS-CoV-2, we found that these genes were enriched in particular pathways. Moreover, these human genes have the codon composition perfectly correlated with the SARS-CoV-2, and were extraordinarily highly expressed in human lung tissues, demonstrating that the SARS-CoV-2 genes have similar GC usage as compared to the lung expressed human genes. RSCU (relative synonymous codon usage) and CAI (codon adaptation index) profiles further support the matching between SARS-CoV-2 and lungs. Our study indicates that SARS-CoV-2 might have adapted to the human lung environment by observing the high correlation between GC usage of SARS-CoV-2 and human lung genes, which suggests the GC content of SARS-CoV-2 is optimized to take advantage of human lung tissues.