Translational adaptation of human viruses to the tissues they infect

mRNA Technology and Mucosal Immunization

Current mRNA vaccines are mainly administered via intramuscular injection, which induces good systemic immunity but limited mucosal immunity. Achieving mucosal immunity through mRNA vaccination could diminish pathogen replication at the entry site and reduce interhuman transmission. However, delivering mRNA vaccines to mucosae faces challenges like mRNA degradation, poor entry into cells, and reactogenicity. Encapsulating mRNA in extracellular vesicles may protect the mRNA and reduce reactogenicity, making mucosal mRNA vaccines possible. Plant-derived extracellular vesicles from edible fruits have been investigated as mRNA carriers. Studies in animals show that mRNA vehiculated in orange-derived extracellular vesicles can elicit both systemic and mucosal immune responses when administered by the oral, nasal, or intramuscular routes. Once lyophilized, these products show remarkable stability. The optimization of mRNA to improve translation efficiency, immunogenicity, reactogenicity, and stability can be obtained through adjustments of the 5'cap region, poly-A tail, codons selection, and the use of nucleoside analogues. Recent studies have also proposed self-amplifying RNA vaccines containing an RNA polymerase as well as circular mRNA constructs. Data from parenterally primed animals demonstrate the efficacy of nasal immunization with non-adjuvanted protein, and studies in humans indicate that the combination of a parenteral vaccine with the natural exposure of mucosae to the same antigen provides protection and reduces transmission. Hence, mucosal mRNA vaccination would be beneficial at least in organisms pre-treated with parenteral vaccines. This practice could have wide applications for the treatment of infectious diseases.

Proteolytic cleavage and inactivation of the TRMT1 tRNA modification enzyme by SARS-CoV-2 main protease

Nonstructural protein 5 (Nsp5) is the main protease of SARS-CoV-2 that cleaves viral polyproteins into individual polypeptides necessary for viral replication. Here, we show that Nsp5 binds and cleaves human tRNA methyltransferase 1 (TRMT1), a host enzyme required for a prevalent post-transcriptional modification in tRNAs. Human cells infected with SARS-CoV-2 exhibit a decrease in TRMT1 protein levels and TRMT1-catalyzed tRNA modifications, consistent with TRMT1 cleavage and inactivation by Nsp5. Nsp5 cleaves TRMT1 at a specific position that matches the consensus sequence of SARS-CoV-2 polyprotein cleavage sites, and a single mutation within the sequence inhibits Nsp5-dependent proteolysis of TRMT1. The TRMT1 cleavage fragments exhibit altered RNA binding activity and are unable to rescue tRNA modification in TRMT1-deficient human cells. Compared to wild-type human cells, TRMT1-deficient human cells infected with SARS-CoV-2 exhibit reduced levels of intracellular viral RNA. These findings provide evidence that Nsp5-dependent cleavage of TRMT1 and perturbation of tRNA modification patterns contribute to the cellular pathogenesis of SARS-CoV-2 infection.

Ten computational challenges in human virome studies

In recent years, substantial advancements have been achieved in understanding the diversity of the human virome and its intricate roles in human health and diseases. Despite this progress, the field of human virome research remains nascent, primarily hindered by the absence of effective methods, particularly in the domain of computational tools. This perspective systematically outlines ten computational challenges spanning various types of virome studies. These challenges arise due to the vast diversity of viromes, the absence of a universal marker gene in viral genomes, the low abundance of virus populations, the remote or minimal homology of viral proteins to known proteins, and the highly dynamic and heterogeneous nature of viromes. For each computational challenge, we discuss the underlying reasons, current research progress, and potential solutions. The resolution of these challenges necessitates ongoing collaboration among computational scientists, virologists, and multidisciplinary experts. In essence, this perspective serves as a comprehensive guide for directing computational efforts in human virome studies.

Establishment of immortalized Egyptian Rousettus bat cell lines

The Egyptian Rousettus bat (Rousettus aegyptiacus) is a common fruit bat species that is distributed mainly in Africa and the Middle East. Bats serve as reservoir hosts for numerous pathogens. Human activities, such as hunting bats for food, managing vermin, and causing habitat loss, elevate the likelihood of transmission of bat pathogens to humans and other animals. Consequently, bat cell lines play a crucial role as research materials for investigating viral pathogens. However, the inherent limitation of finite cell division in primary cells necessitates the use of immortalized cells derived from various bat tissues. Herein, we successfully established six fibroblast cell lines derived from an infant bat heart and lungs and an elderly bat heart. Three of the six cell lines, called K4DT cells, were transduced by a combination of cell cycle regulators, mutant cyclin-dependent kinase 4, cyclin D1, and human telomerase reverse transcriptase. The other three cell lines, named SV40 cells, were transfected with simian virus 40 large T antigen. Transgene protein expression was detected in the transduced cells. All three K4DT cell lines and one lung-derived SV40 cell line were virtually immortalized and nearly maintained the normal diploid karyotypes. However, the two other heart-derived SV40 cell lines had aberrant karyotypes and the young bat-derived cell line stopped proliferating at approximately 40 population doublings. These bat cell lines are valuable for studying pathogen genomics and biology.

The human papillomavirus late life cycle and links to keratinocyte differentiation

Regulation of human papillomavirus (HPV) gene expression is tightly linked to differentiation of the keratinocytes the virus infects. HPV late gene expression is confined to the cells in the upper layers of the epithelium where the virus capsid proteins are synthesized. As these proteins are highly immunogenic, and the upper epithelium is an immune-privileged site, this spatial restriction aids immune evasion. Many decades of work have contributed to the current understanding of how this restriction occurs at a molecular level. This review will examine what is known about late gene expression in HPV-infected lesions and will dissect the intricacies of late gene regulation. Future directions for novel antiviral approaches will be highlighted.

The Typical tRNA Co-Expresses Multiple 5' tRNA Halves Whose Sequences and Abundances Depend on Isodecoder and Isoacceptor and Change with Tissue Type, Cell Type, and Disease

Transfer RNA-derived fragments (tRFs) are noncoding RNAs that arise from either mature transfer RNAs (tRNAs) or their precursors. One important category of tRFs comprises the tRNA halves, which are generated through cleavage at the anticodon. A given tRNA typically gives rise to several co-expressed 5'-tRNA halves (5'-tRHs) that differ in the location of their 3' ends. These 5'-tRHs, even though distinct, have traditionally been treated as indistinguishable from one another due to their near-identical sequences and lengths. We focused on co-expressed 5'-tRHs that arise from the same tRNA and systematically examined their exact sequences and abundances across 10 different human tissues. To this end, we manually curated and analyzed several hundred human RNA-seq datasets from NCBI's Sequence Run Archive (SRA). We grouped datasets from the same tissue into their own collection and examined each group separately. We found that a given tRNA produces different groups of co-expressed 5'-tRHs in different tissues, different cell lines, and different diseases. Importantly, the co-expressed 5'-tRHs differ in their sequences, absolute abundances, and relative abundances, even among tRNAs with near-identical sequences from the same isodecoder or isoacceptor group. The findings suggest that co-expressed 5'-tRHs that are produced from the same tRNA or closely related tRNAs have distinct, context-dependent roles. Moreover, our analyses show that cell lines modeling the same tissue type and disease may not be interchangeable when it comes to experimenting with tRFs.

Comprehensive quantitative modeling of translation efficiency in a genome-reduced bacterium

Translation efficiency has been mainly studied by ribosome profiling, which only provides an incomplete picture of translation kinetics. Here, we integrated the absolute quantifications of tRNAs, mRNAs, RNA half-lives, proteins, and protein half-lives with ribosome densities and derived the initiation and elongation rates for 475 genes (67% of all genes), 73 with high precision, in the bacterium Mycoplasma pneumoniae (Mpn). We found that, although the initiation rate varied over 160-fold among genes, most of the known factors had little impact on translation efficiency. Local codon elongation rates could not be fully explained by the adaptation to tRNA abundances, which varied over 100-fold among tRNA isoacceptors. We provide a comprehensive quantitative view of translation efficiency, which suggests the existence of unidentified mechanisms of translational regulation in Mpn.

"Transfer" of power: The intersection of DNA virus infection and tRNA biology

Transfer RNAs (tRNAs) are at the heart of the molecular biology central dogma, functioning to decode messenger RNAs into proteins. As obligate intracellular parasites, viruses depend on the host translation machinery, including host tRNAs. Thus, the ability of a virus to fine-tune tRNA expression elicits the power to impact the outcome of infection. DNA viruses commonly upregulate the output of RNA polymerase III (Pol III)-dependent transcripts, including tRNAs. Decades after these initial discoveries we know very little about how mature tRNA pools change during viral infection, as tRNA sequencing methodology has only recently reached proficiency. Here, we review perturbation of tRNA biogenesis by DNA virus infection, including an emerging player called tRNA-derived fragments (tRFs). We discuss how tRNA dysregulation shifts the power landscape between the host and virus, highlighting the potential for tRNA-based antivirals as a future therapeutic.

Human genes with codon usage bias similar to that of the nonstructural protein 1 gene of influenza A viruses are conjointly involved in the infectious pathogenesis of influenza A viruses

Molecular mechanisms of the non-structural protein 1 (NS1) in influenza A-induced pathological changes remain ambiguous. This study explored the pathogenesis of human infection by influenza A viruses (IAVs) through identifying human genes with codon usage bias (CUB) similar to NS1 gene of these viruses based on the relative synonymous codon usage (RSCU). CUB of the IAV subtypes H1N1, H3N2, H3N8, H5N1, H5N2, H5N8, H7N9 and H9N2 was analyzed and the correlation of RSCU values of NS1 sequences with those of the human genes was calculated. The CUB of NS1 was uneven and codons ending with A/U were preferred. The ENC-GC3 and neutrality plots suggested natural selection as the main determinant for CUB. The RCDI, CAI and SiD values showed that the viruses had a high degree of adaptability to human. A total of 2155 human genes showed significant RSCU-based correlation (p < 0.05 and r > 0.5) with NS1 coding sequences and was considered as human genes with CUB similar to NS1 gene of IAV subtypes. Differences and similarities in the subtype-specific human protein–protein interaction (PPI) networks and their functions were recorded among IAVs subtypes, indicating that NS1 of each IAV subtype has a specific pathogenic mechanism. Processes and pathways involved in influenza, transcription, immune response and cell cycle were enriched in human gene sets retrieved based on the CUB of NS1 gene of IAV subtypes. The present work may advance our understanding on the mechanism of NS1 in human infections of IAV subtypes and shed light on the therapeutic options.

Distinct codon usage bias evolutionary patterns between weakly and strongly virulent respiratory viruses

Human respiratory viruses are of vastly different virulence, giving rise to symptoms ranging from common cold to severe pneumonia or even death. Although this most likely impacts molecular evolution of the corresponding viruses, the specific differences in their evolutionary patterns remain largely unknown. By comparing structural and nonstructural genes within respiratory viruses, greater similarities in codon usage bias (CUB) between nonstructural genes and humans were observed in weakly virulent viruses, whereas in strongly virulent viruses, it was structural genes whose CUBs were more similar to that of humans. Further comparisons between genomes of weakly and strongly virulent coronaviruses revealed greater similarities in CUBs between strongly virulent viruses and humans. Finally, using phylogenetic independent contrasts, dissimilation of viral CUB from that of humans was observed in SARS-CoV-2. Our work revealed distinct CUB evolutionary patterns between weakly and strongly virulent viruses, a previously unrecognized interaction between CUB and virulence in respiratory viruses.

Causes and Consequences of Purifying Selection on SARS-CoV-2

Owing to a lag between a deleterious mutation's appearance and its selective removal, gold-standard methods for mutation rate estimation assume no meaningful loss of mutations between parents and offspring. Indeed, from analysis of closely related lineages, in SARS-CoV-2, the Ka/Ks ratio was previously estimated as 1.008, suggesting no within-host selection. By contrast, we find a higher number of observed SNPs at 4-fold degenerate sites than elsewhere and, allowing for the virus's complex mutational and compositional biases, estimate that the mutation rate is at least 49-67% higher than would be estimated based on the rate of appearance of variants in sampled genomes. Given the high Ka/Ks one might assume that the majority of such intrahost selection is the purging of nonsense mutations. However, we estimate that selection against nonsense mutations accounts for only ∼10% of all the "missing" mutations. Instead, classical protein-level selective filters (against chemically disparate amino acids and those predicted to disrupt protein functionality) account for many missing mutations. It is less obvious why for an intracellular parasite, amino acid cost parameters, notably amino acid decay rate, is also significant. Perhaps most surprisingly, we also find evidence for real-time selection against synonymous mutations that move codon usage away from that of humans. We conclude that there is common intrahost selection on SARS-CoV-2 that acts on nonsense, missense, and possibly synonymous mutations. This has implications for methods of mutation rate estimation, for determining times to common ancestry and the potential for intrahost evolution including vaccine escape.

Age and gender differences in ACE2 and TMPRSS2 expressions in oral epithelial cells

Background

SARS-CoV-2, which has brought a huge negative impact on the world since the end of 2019, is reported to invade cells using the spike (S) protein to bind to angiotensin-converting enzyme II (ACE2) receptors on human cells while the transmembrane protease serine 2 (TMPRSS2) is the key protease that activates the S protein, which greatly facilitates the entry of SARS-CoV-2 into target cells. In our previous study, it was observed that the positive rate of SARS-CoV-2 nucleic acids in saliva was higher in male and the elderly COVID-19 patients, suggesting that the susceptibility of oral tissues to SARS-CoV-2 may be related to gender and age. This research aimed to further investigate the SARS-CoV-2 susceptibility in oral tissues and influencing factors from the perspective of ACE2 and TMPRSS2, which were two proteins closely associated with SARS-CoV-2 infection.

Methods

Immunofluorescence was used to find the localization of ACE2 and TMPRSS2 in oral mucosal tissues. Transcriptomic sequencing data of several datasets were then collected to analysis the relationship between the expressions of ACE2 and TMPRSS2 with the age and gender of patients. Furthermore, oral tissues from patients with different ages and genders were collected. Immunohistochemistry staining, qRT-PCR and western blot were performed to explore the relationship between expression levels of ACE2 and TMPRSS2 and patient age as well as gender.

Results

The results showed that the two proteins were able to be co-expressed in the epithelial cells of oral tissues, and their expression levels were higher in the relatively elderly group than those in relatively younger group. Male oral epithelial cells exhibited higher level of TMPRSS2.

Conclusions

Our findings comprehensively confirmed the existence of ACE2 and TMPRSS2 in oral tissues and clarify the relationship between the expression levels with human age and gender for the first time, providing evidence for possible entry routes of SARS-CoV-2 and the influencing factors of SARS-CoV-2 colonization in oral cavity. Thus, the oral mucosa might be at potential risk of infection by SARS-CoV-2, especially in male or elderly patients. Using saliva to detect the nucleic acids of SARS-CoV-2 may be more accurate for elder male COVID-19 patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-03037-4.

The Functional Consequences of the Novel Ribosomal Pausing Site in SARS-CoV-2 Spike Glycoprotein RNA

The SARS-CoV-2 Spike glycoprotein (S protein) acquired a unique new 4 amino acid -PRRA- insertion sequence at amino acid residues (aa) 681-684 that forms a new furin cleavage site in S protein as well as several new glycosylation sites. We studied various statistical properties of the -PRRA- insertion at the RNA level (CCUCGGCGGGCA). The nucleotide composition and codon usage of this sequence are different from the rest of the SARS-CoV-2 genome. One of such features is two tandem CGG codons, although the CGG codon is the rarest codon in the SARS-CoV-2 genome. This suggests that the insertion sequence could cause ribosome pausing as the result of these rare codons. Due to population variants, the Nextstrain divergence measure of the CCU codon is extremely large. We cannot exclude that this divergence might affect host immune responses/effectiveness of SARS-CoV-2 vaccines, possibilities awaiting further investigation. Our experimental studies show that the expression level of original RNA sequence "wildtype" spike protein is much lower than for codon-optimized spike protein in all studied cell lines. Interestingly, the original spike sequence produces a higher titer of pseudoviral particles and a higher level of infection. Further mutagenesis experiments suggest that this dual-effect insert, comprised of a combination of overlapping translation pausing and furin sites, has allowed SARS-CoV-2 to infect its new host (human) more readily. This underlines the importance of ribosome pausing to allow efficient regulation of protein expression and also of cotranslational subdomain folding.

Transcription, mRNA export and immune evasion shape the codon usage of viruses

The nucleotide composition, dinucleotide composition, and codon usage of many viruses differs from their hosts. These differences arise because viruses are subject to unique mutation and selection pressures that do not apply to host genomes; however, the molecular mechanisms that underlie these evolutionary forces are unclear. Here, we analysed the patterns of codon usage in 1,520 vertebrate-infecting viruses, focusing on parameters known to be under selection and associated with gene regulation. We find that GC content, dinucleotide content, and splicing and m6A modification-related sequence motifs are associated with the type of genetic material (DNA or RNA), strandedness, and replication compartment of viruses. In an experimental follow-up, we find that the effects of GC content on gene expression depend on whether the genetic material is delivered to the cell as DNA or mRNA, whether it is transcribed by endogenous or exogenous RNA polymerase, and whether transcription takes place in the nucleus or cytoplasm. Our results suggest that viral codon usage cannot be explained by a simple adaptation to the codon usage of the host - instead, it reflects the combination of multiple selective and mutational pressures, including the need for efficient transcription, export, and immune evasion.