DNA methylation and Mycoplasma genomes

Extreme Genomic CpG Deficiency in SARS-CoV-2 and Evasion of Host Antiviral Defense

Wild mammalian species, including bats, constitute the natural reservoir of betacoronavirus (including SARS, MERS, and the deadly SARS-CoV-2). Different hosts or host tissues provide different cellular environments, especially different antiviral and RNA modification activities that can alter RNA modification signatures observed in the viral RNA genome. The zinc finger antiviral protein (ZAP) binds specifically to CpG dinucleotides and recruits other proteins to degrade a variety of viral RNA genomes. Many mammalian RNA viruses have evolved CpG deficiency. Increasing CpG dinucleotides in these low-CpG viral genomes in the presence of ZAP consistently leads to decreased viral replication and virulence. Because ZAP exhibits tissue-specific expression, viruses infecting different tissues are expected to have different CpG signatures, suggesting a means to identify viral tissue-switching events. The author shows that SARS-CoV-2 has the most extreme CpG deficiency in all known betacoronavirus genomes. This suggests that SARS-CoV-2 may have evolved in a new host (or new host tissue) with high ZAP expression. A survey of CpG deficiency in viral genomes identified a virulent canine coronavirus (alphacoronavirus) as possessing the most extreme CpG deficiency, comparable with that observed in SARS-CoV-2. This suggests that the canine tissue infected by the canine coronavirus may provide a cellular environment strongly selecting against CpG. Thus, viral surveys focused on decreasing CpG in viral RNA genomes may provide important clues about the selective environments and viral defenses in the original hosts.

Improving Phylogenetic Signals of Mitochondrial Genes Using a New Method of Codon Degeneration

Recovering deep phylogeny is challenging with animal mitochondrial genes because of their rapid evolution. Codon degeneration decreases the phylogenetic noise and bias by aiming to achieve two objectives: (1) alleviate the bias associated with nucleotide composition, which may lead to homoplasy and long-branch attraction, and (2) reduce differences in the phylogenetic results between nucleotide-based and amino acid (AA)-based analyses. The discrepancy between nucleotide-based analysis and AA-based analysis is partially caused by some synonymous codons that differ more from each other at the nucleotide level than from some nonsynonymous codons, e.g., Leu codon TTR in the standard genetic code is more similar to Phe codon TTY than to synonymous CTN codons. Thus, nucleotide similarity conflicts with AA similarity. There are many such examples involving other codon families in various mitochondrial genetic codes. Proper codon degeneration will make synonymous codons more similar to each other at the nucleotide level than they are to nonsynonymous codons. Here, I illustrate a "principled" codon degeneration method that achieves these objectives. The method was applied to resolving the mammalian basal lineage and phylogenetic position of rheas among ratites. The codon degeneration method was implemented in the user-friendly and freely available DAMBE software for all known genetic codes (genetic codes 1 to 33).

Bioinformatics and Translation Elongation

Codon usage depends on mutation bias, tRNA-mediated selection, and the need for high efficiency and accuracy in translation. One codon in a synonymous codon family is often strongly over-used, especially in highly expressed genes, which often leads to a high dN/dS ratio because dS is very small. Many different codon usage indices have been proposed to measure codon usage and codon adaptation. Sense codon could be misread by release factors and stop codons misread by tRNAs, which also contribute to codon usage in rare cases. This chapter outlines the conceptual framework on codon evolution, illustrates codon-specific and gene-specific codon usage indices, and presents their applications. A new index for codon adaptation that accounts for background mutation bias (Index of Translation Elongation) is presented and contrasted with codon adaptation index (CAI) which does not consider background mutation bias. They are used to re-analyze data from a recent paper claiming that translation elongation efficiency matters little in protein production. The reanalysis disproves the claim.

Bioinformatics and Drug Discovery

Bioinformatic analysis can not only accelerate drug target identification and drug candidate screening and refinement, but also facilitate characterization of side effects and predict drug resistance. High-throughput data such as genomic, epigenetic, genome architecture, cistromic, transcriptomic, proteomic, and ribosome profiling data have all made significant contribution to mechanismbased drug discovery and drug repurposing. Accumulation of protein and RNA structures, as well as development of homology modeling and protein structure simulation, coupled with large structure databases of small molecules and metabolites, paved the way for more realistic protein-ligand docking experiments and more informative virtual screening. I present the conceptual framework that drives the collection of these high-throughput data, summarize the utility and potential of mining these data in drug discovery, outline a few inherent limitations in data and software mining these data, point out news ways to refine analysis of these diverse types of data, and highlight commonly used software and databases relevant to drug discovery.

The evolution of genomic GC content undergoes a rapid reversal within the genus Plasmodium

The genome of the malarial parasite Plasmodium falciparum is extremely AT rich. This bias toward a low GC content is a characteristic of several, but not all, species within the genus Plasmodium. We compared 4283 orthologous pairs of protein-coding sequences between Plasmodium falciparum and the less AT-biased Plasmodium vivax. Our results indicate that the common ancestor of these two species was also extremely AT rich. This means that, although there was a strong bias toward A+T during the early evolution of the ancestral Plasmodium lineage, there was a subsequent reversal of this trend during the more recent evolution of some species, such as P. vivax. Moreover, we show that not only is the P. vivax genome losing its AT richness, it is actually gaining a very significant degree of GC richness. This example illustrates the potential volatility of nucleotide content during the course of molecular evolution. Such reversible fluxes in nucleotide content within lineages could have important implications for phylogenetic reconstruction based on molecular sequence data.

The effect of mutation and selection on codon adaptation in Escherichia coli bacteriophage

Studying phage codon adaptation is important not only for understanding the process of translation elongation, but also for reengineering phages for medical and industrial purposes. To evaluate the effect of mutation and selection on phage codon usage, we developed an index to measure selection imposed by host translation machinery, based on the difference in codon usage between all host genes and highly expressed host genes. We developed linear and nonlinear models to estimate the C→T mutation bias in different phage lineages and to evaluate the relative effect of mutation and host selection on phage codon usage. C→T-biased mutations occur more frequently in single-stranded DNA (ssDNA) phages than in double-stranded DNA (dsDNA) phages and affect not only synonymous codon usage, but also nonsynonymous substitutions at second codon positions, especially in ssDNA phages. The host translation machinery affects codon adaptation in both dsDNA and ssDNA phages, with a stronger effect on dsDNA phages than on ssDNA phages. Strand asymmetry with the associated local variation in mutation bias can significantly interfere with codon adaptation in both dsDNA and ssDNA phages.

DAMBE5: a comprehensive software package for data analysis in molecular biology and evolution

Since its first release in 2001 as mainly a software package for phylogenetic analysis, data analysis for molecular biology and evolution (DAMBE) has gained many new functions that may be classified into six categories: 1) sequence retrieval, editing, manipulation, and conversion among more than 20 standard sequence formats including MEGA, NEXUS, PHYLIP, GenBank, and the new NeXML format for interoperability, 2) motif characterization and discovery functions such as position weight matrix and Gibbs sampler, 3) descriptive genomic analysis tools with improved versions of codon adaptation index, effective number of codons, protein isoelectric point profiling, RNA and protein secondary structure prediction and calculation of minimum folding energy, and genomic skew plots with optimized window size, 4) molecular phylogenetics including sequence alignment, testing substitution saturation, distance-based, maximum parsimony, and maximum-likelihood methods for tree reconstructions, testing the molecular clock hypothesis with either a phylogeny or with relative-rate tests, dating gene duplication and speciation events, choosing the best-fit substitution models, and estimating rate heterogeneity over sites, 5) phylogeny-based comparative methods for continuous and discrete variables, and 6) graphic functions including secondary structure display, optimized skew plot, hydrophobicity plot, and many other plots of amino acid properties along a protein sequence, tree display and drawing by dragging nodes to each other, and visual searching of the maximum parsimony tree. DAMBE features a graphic, user-friendly, and intuitive interface and is freely available from http://dambe.bio.uottawa.ca (last accessed April 16, 2013).

In vivo function of airway epithelial TLR2 in host defense against bacterial infection

Decreased Toll-like receptor 2 (TLR2) expression has been reported in patients with chronic obstructive pulmonary disease and in a murine asthma model, which may predispose the hosts to bacterial infections, leading to disease exacerbations. Since airway epithelial cells serve as the first line of respiratory mucosal defense, the present study aimed to reveal the role of airway epithelial TLR2 signaling to lung bacterial [i.e., Mycoplasma pneumoniae (Mp)] clearance. In vivo TLR2 gene transfer via intranasal inoculation of adenoviral vector was performed to reconstitute TLR2 expression in airway epithelium of TLR2(-/-) BALB/c mice, with or without ensuing Mp infection. TLR2 and lactotransferrin (LTF) expression in airway epithelial cells and lung Mp load were assessed. Adenovirus-mediated TLR2 gene transfer to airway epithelial cells of TLR2(-/-) mice reconstituted 30-40% TLR2 expression compared with TLR2(+/+) cells. Such airway epithelial TLR2 reconstitution in TLR2(-/-) mice significantly reduced lung Mp load (an appropriate 45% reduction), coupled with elevated LTF expression. LTF expression in mice was shown to be mainly dependent on TLR2 signaling in response to Mp infection. Exogenous human LTF protein dose-dependently decreased lung bacterial load in Mp-infected TLR2(-/-) mice. In addition, human LTF protein directly dose-dependently decreased Mp levels in vitro. These data indicate that reconstitution of airway epithelial TLR2 signaling in TLR2(-/-) mice significantly restores lung defense against bacteria (e.g., Mp) via increased lung antimicrobial protein LTF production. Our findings may offer a deliverable approach to attenuate bacterial infections in airways of asthma or chronic obstructive pulmonary disease patients with impaired TLR2 function.

Information-theoretic indices and an approximate significance test for testing the molecular clock hypothesis with genetic distances

Distance-based phylogenetic methods are widely used in biomedical research. However, distance-based dating of speciation events and the test of the molecular clock hypothesis are relatively underdeveloped. Here I develop an approximate test of the molecular clock hypothesis for distance-based trees, as well as information-theoretic indices that have been used frequently in model selection, for use with distance matrices. The results are in good agreement with the conventional sequence-based likelihood ratio test. Among the information-theoretic indices, AICu is the most consistent with the sequence-based likelihood ratio test. The confidence in model selection by the indices can be evaluated by bootstrapping. I illustrate the usage of the indices and the approximate significance test with both empirical and simulated sequences. The tests show that distance matrices from protein gel electrophoresis and from genome rearrangement events do not violate the molecular clock hypothesis, and that the evolution of the third codon position conforms to the molecular clock hypothesis better than the second codon position in vertebrate mitochondrial genes. I outlined evolutionary distances that are appropriate for phylogenetic reconstruction and dating.

Toll-like receptor 2 down-regulation in established mouse allergic lungs contributes to decreased mycoplasma clearance

RATIONALE

Respiratory Mycoplasma pneumoniae (Mp) infection is involved in asthma pathobiology, but whether the established allergic airway inflammation compromises lung innate immunity and subsequently predisposes patients with asthma to Mp infection remains unknown.

OBJECTIVES

To test whether the established allergic airway inflammation compromises host innate immunity (e.g., Toll-like receptor 2 [TLR2]) to hinder the elimination of Mp from the lungs.

METHODS

We used mouse models of ovalbumin (OVA)-induced allergic airway inflammation with an ensuing Mp infection, and cultures of mouse primary lung dendritic cells (DCs) and bone marrow-derived DCs.

MEASUREMENTS AND MAIN RESULTS

Lung Mp clearance in allergic mice and TLR2 and IL-6 levels in lung cells, including DCs as well as cultured primary lung DCs and bone marrow-derived DCs, were assessed. The established OVA-induced allergic airway inflammation, or the prominent Th2 cytokines IL-4 and IL-13, inhibited TLR2 expression and IL-6 production in lung cells, including lung DCs, and eventually led to impaired host defense against Mp. Studies in IL-6 knockout mice indicated that IL-6 directly promoted Mp clearance from the lungs. IL-4- and IL-13-induced suppression of TLR2 was mediated by inhibiting nuclear factor-kappaB activation through signal transducer and activator of transcription 6 (STAT6) signaling pathway.

CONCLUSIONS

The established OVA-induced allergic airway inflammation impairs TLR2 expression and host defense cytokine (e.g., IL-6) production, and subsequently delays lung bacterial clearance. This could offer novel therapeutic strategies to reinstate TLR2 activation by using TLR2 ligands and/or blocking IL-4 and IL-13 to ameliorate persisting respiratory bacterial infections in allergic lungs.