Comparison of dinucleotide frequency and codon usage in Toxoplasma and Plasmodium: evolutionary implications

Deciphering the co-adaptation of codon usage between respiratory coronaviruses and their human host uncovers candidate therapeutics for COVID-19

Coronavirus disease 2019 (COVID-19) has caused thousands of deaths worldwide and has become an urgent public health concern. The extraordinary interhuman transmission of this disease has urged scientists to examine the various facets of its pathogenic agent, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Herein, based on publicly available genomic data, we analyzed the codon usage co-adaptation profiles of SARS-CoV-2 and other respiratory coronaviruses (CoVs) with their human host, identified CoV-responsive human genes and their functional roles on the basis of both the relative synonymous codon usage (RSCU)-based correlation of viral genes with human genes and differential gene expression analysis, and predicted potential drugs for COVID-19 treatment based on these genes. The relatively high codon adaptation index (CAI) values (>0.70) signposted the gene expressivity efficiency of CoVs in human. The ENc-GC3 plot indicated that SARS-CoV-2 genome was under strict selection pressure while SARS-CoV and MERS-CoV were under selection and mutational pressures. The RSCU-based correlation analysis indicated that the viral genomes shared similar codons with a panoply of human genes. The merging of RSCU-based correlation data and SARS-CoV-2-responsive differentially expressed genes allowed the identification of human genes potentially affected by SARS-CoV-2 infection. Functional enrichment analysis indicated that these genes were enriched in biological processes and pathways related to host response to viral infection and immune response. Using the drug-gene interaction database, we screened a list of drugs that could target these genes as potential COVID-19 therapeutics. Our findings not only will contribute in vaccine development but also provide a useful set of drugs that could guide practitioners in strategical monitoring of COVID-19. We recommend practitioners to scrupulously screen this list of predicted drugs in order to authenticate those qualified for treating COVID-19 symptoms.

Comprehensive analysis of prokaryotic mechanosensation genes: their characteristics in codon usage

In the present study, we examined GC nucleotide composition, relative synonymous codon usage (RSCU), effective number of codons (ENC), codon adaptation index (CAI) and gene length for 308 prokaryotic mechanosensitive ion channel (MSC) genes from six evolutionary groups: Euryarchaeota, Actinobacteria, Alphaproteobacteria, Betaproteobacteria, Firmicutes, and Gammaproteobacteria. Results showed that: (1) a wide variation of overrepresentation of nucleotides exists in the MSC genes; (2) codon usage bias varies considerably among the MSC genes; (3) both nucleotide constraint and gene length play an important role in shaping codon usage of the bacterial MSC genes; and (4) synonymous codon usage of prokaryotic MSC genes is phylogenetically conserved. Knowledge of codon usage in prokaryotic MSC genes may benefit from the study of the MSC genes in eukaryotes in which few MSC genes have been identified and functionally analysed.

Toxoplasma gondii expresses two distinct lactate dehydrogenase homologous genes during its life cycle in intermediate hosts

Two Toxoplasma gondii genes were characterized that are differentially expressed during the parasite's life cycle. The genes named LDH1 and LDH2, respectively, encode polypeptides similar to the enzyme lactate dehydrogenase (LDH; L-lactate:NAD+ oxidoreductase, EC 1.1.1.27) from a variety of organisms. They show 64.0% nucleotide identity in the coding region and both have an intron at the same relative position. The deduced amino acid sequences of LDH1 and LDH2 share 71.1% identity. LDH1 and LDH2 are most similar to an LDH of Plasmodium falciparum (46.5% and 48.5% amino acid identities, respectively). The mRNA of LDH2 was only detected in the bradyzoite stage, while the mRNA of LDH1 was detected in both the bradyzoite and tachyzoite stages. However, by isoelectric focusing and immunoblot analysis, only one LDH isoform was found to be expressed in each stage. Furthermore, the expression of a reporter gene carrying chloramphenicol acetyltransferase (CAT) coding sequence and the putative LDH2 promoter sequence was significantly up-regulated by growing parasites in tissue culture in media with alkaline pH (pH 8.2, a condition known to induce the expression of bradyzoite-specific antigens), while the expression of a CAT reporter construct carrying the putative LDH1 promoter sequence was down-regulated by similar treatment. These results indicate that LDH expression is developmentally regulated in T. gondii and suggest a possible correlation between stage conversion and alteration in carbohydrate or energy metabolism in this parasite.

Possible implications of CpG avoidance in the flatworm Schistosoma mansoni

We report the analysis of the biases of CpG, TpG, and CpA of all the DNA sequences data from the Trematode Schistosoma mansoni. Our results show CpG avoidance whereas TpG and CpA frequencies are over the expected values. These characteristics are similar to the biases displayed by methylated genomes, but in platyhelminths 5mC has not been detected by biochemical methods. The possible implications of this CpG shortage are discussed.

Analysis of dinucleotide frequency and codon usage in the phylum Apicomplexa

Dinucleotide frequency (DiF) and codon usage (cu) were analysed in gene sequences from four parasitic protozoa, Babesia bovis, Theileria parva, Toxoplasma gondii and Eimeria tenella, of the phylum Apicomplexa. In keeping with the 'genome hypothesis', cu was found to be non-random and species specific in these organisms, although cu among members of the same subclass was found to be very similar. Several low-usage (lu) codons were identified, and the usage of lu codons appears to be related to the taxonomic position of the organisms under study. A comparison of the observed/expected DiF ratios obtained from gene coding regions revealed a low frequency of the TA and CG dinucleotides in all organisms studied. A comparison of these DiF ratios with those found in rRNA-encoding genes and in introns, showed that in the parasites, B. bovis and Th. parva (representing the piroplasms), the low frequency of dinucleotides appeared to be the result of coding pressure alone. In T. gondii and E. tenella (representing the coccidia), however, coding pressure could not completely explain differences in DiF.

Cloning, expression and nucleotide sequence of the gene fragment encoding an antigenic portion of the nucleoside triphosphate hydrolase of Toxoplasma gondii

Toxoplasma gondii expresses high levels of an active nucleoside triphosphate hydrolase (NTPase; EC 3.6.1.3) with several unique properties. It has been detected as a circulating antigen in mice, making it an ideal candidate for diagnostic tests for toxoplasmosis. A cDNA library constructed from T. gondii poly(A)+RNA was made in lambda gt11. One hundred thousand members of this library were immunoscreened with a rabbit polyclonal antibody to the purified NTPase. Six positive clones were subcloned into the plasmid, pGEX-IN, and the inserts were restriction-mapped. All clones had identical partial restriction enzyme maps. One insert was subcloned into M13mp18 and sequencing by the deoxy/dideoxy method showed an NTPase-encoding gene (ntp) fragment of 571 bp. The insert was also purified, radiolabelled, and used to hybridize to Northern blots of tachyzoite RNA and quantitative Southern blots of tachyzoite DNA. Northern blotting revealed that the NTPase mRNA was in great abundance and had a length of about 2800 nucleotides. Southern blotting showed a gene copy number of between one and five, and the possibility that ntp is tandemly repeated over a large length of DNA. The NTPase was expressed as a glutathione S-transferase (EC 2.5.1.18) fusion protein of about 50 kDa, which reacted with polyclonal rabbit antibody on Western blotting.