Nucleotide Usage Research Articles

Mycobacteria that cause tuberculosis have retained ancestrally acquired genes for the biosynthesis of chemically diverse terpene nucleosides.

Mycobacterium tuberculosis (Mtb) releases the unusual terpene nucleoside 1-tuberculosinyladenosine (1-TbAd) to block lysosomal function and promote survival in human macrophages. Using conventional approaches, we found that genes Rv3377c and Rv3378c, but not Rv3376, were necessary for 1-TbAd biosynthesis. Here, we introduce linear models for mass spectrometry (limms) software as a next-generation lipidomics tool to study the essential functions of lipid biosynthetic enzymes on a whole-cell basis. Using limms, whole-cell lipid profiles deepened the phenotypic landscape of comparative mass spectrometry experiments and identified a large family of approximately 100 terpene nucleoside metabolites downstream of Rv3378c. We validated the identity of previously unknown adenine-, adenosine-, and lipid-modified tuberculosinol-containing molecules using synthetic chemistry and collisional mass spectrometry, including comprehensive profiling of bacterial lipids that fragment to adenine. We tracked terpene nucleoside genotypes and lipid phenotypes among Mycobacterium tuberculosis complex (MTC) species that did or did not evolve to productively infect either human or nonhuman mammals. Although 1-TbAd biosynthesis genes were thought to be restricted to the MTC, we identified the locus in unexpected species outside the MTC. Sequence analysis of the locus showed nucleotide usage characteristic of plasmids from plant-associated bacteria, clarifying the origin and timing of horizontal gene transfer to a pre-MTC progenitor. The data demonstrated correlation between high level terpene nucleoside biosynthesis and mycobacterial competence for human infection, and 2 mechanisms of 1-TbAd biosynthesis loss. Overall, the selective gain and evolutionary retention of tuberculosinyl metabolites in modern species that cause human TB suggest a role in human TB disease, and the newly discovered molecules represent candidate disease-specific biomarkers.

Characterization of the mitochondrial genomes of the Mexican endemic bats Corynorhinus mexicanus and Corynorhinus leonpaniaguae (Chiroptera: Vespertilionidae).

The genus Corynorhinus is composed of four recognized species: C. rafinesquii, C. townsendii, C. mexicanus, and C. leonpaniaguae, the latter two being endemic to Mexico. According to the IUCN, C. mexicanus is considered "Near Threatened", as its populations are dwindling and habitats are affected by anthropogenic disturbance. Corynorhinus leonpaniaguae has not been assigned to an IUCN Red List risk category due to its recent description. In this study, the mitochondrial genomes of C. mexicanus and C. leonpaniaguae were assembled and characterized in detail. The mitochondrial genomes (mtDNA) of C. mexicanus and C. leonpaniaguae have lengths of 16,470 and 16,581bp respectively, with a predominant nucleotide usage of adenine (31.670% and 31.729%, respectively) and thymine (26.15% and 26.18%, respectively). The mtDNA of C. mexicanus and C. leonpaniaguae is composed of 37 coding and non-coding elements: 22 transfer RNAs (tRNA), 13 protein-coding genes (PCGs), two ribosomal RNAs and a non-coding region, the control region, which has a length of 933bp and 1,149bp, respectively. All tRNAs exhibited a cloverleaf secondary structure, with the exception of trn-Ser1 which showed a deletion of the dihydrouridine arm in the two species. All PCGs are subjected to purifying selection, with atp8 being the gene showing the highest Ka/Ks value. These are the first whole mitogenomic resources developed for C. mexicanus and C. leonpaniaguae and enhance our knowledge of the ecology of these species and aid in their conservation.

The systematic codon usage bias has an important effect on genetic adaption in native species

Random mutations increase genetic variety and natural selection enhances adaption over generations. Codon usage biases (CUB) provide clues about the genome adaptation mechanisms of native species and extremophile species. Significant numbers of gene (CDS) of nine classes of endangered, native species, including extremophiles and mesophiles were utilised to compute CUB. Codon usage patterns differ among the lineages of endangered and extremophiles with native species. Polymorphic usage of nucleotides with codon burial suggests parallelism of native species within relatively confined taxonomic groups. Utilizing the deviation pattern of CUB of endangered and native species, I present a calculation parameter to estimate the extinction risk of endangered species. Species diversity and extinction risk are both positively associated with the propensity of random mutation in CDS (Coding DNA sequence). Codon bias tenet profoundly selected and it governs to adaptive evolution of native species.

Synthetic Introns Identify the Novel RNA Splicing Factor GPATCH8 As Required for Mis-Splicing Induced By SF3B1 Mutations

Disclosures Mutations in the core RNA splicing factor SF3B1 are common across MDS, CLL, and clonal hematopoiesis. Prior studies have elucidated that mutations in SF3B1 result in neomorphic widespread changes in splicing due to usage of aberrant intronic branchpoint nucleotides. However, the molecular basis by which mutant SF3B1 induces mis-splicing is not established and therapeutic means to correct mis-splicing due to mutant SF3B1 do not exist. Recently, the aberrant gain-of-function splicing activity of mutant SF3B1 was harnessed to regulate expression of proteins in a mutant selective manner using synthetic version of endogenous mRNA sequences uniquely recognized by mutant SF3B1. Here, we utilized this technology to engineer a synthetic intron derived from MAP3K7 that interrupts the coding sequence of the fluorescent protein mEmerald such that mutant SF3B1 produces lower mEmerald expression than SF3B1 wild-type (WT) cells (Fig.A). This fluorescent splicing reporter was then used to perform positive enrichment whole genome CRISPR screens to identify genes whose deletion corrects SF3B1 mutant aberrant splicing. GPATCH8 was the single most robust hit. GPATCH8 knockout strongly corrected mis-splicing of both synthetic and endogenous MAP3K7 introns (Fig.A). A previously unexplored protein, GPATCH8 has domains characteristic of RNA splicing factors including a G-patch motif, thought to be important in activating RNA helicases. Interestingly, recent work suggests that mutations in SF3B1 result in mis-splicing by disrupting physical interaction of SF3B1 to the RNA helicase DHX15 via the G-patch domain containing protein SUGP1. We therefore compared SUGP1 to GPATCH8 activity, evaluating first their transcriptome-wide binding sites using TRIBE-seq. We found that both proteins predominantly target intronic sequences over other genomic regions suggesting roles in splicing. RNA-seq of SF3B1 mutant and WT cells with or without GPATCH8/SUGP1 deletion revealed that ~33% of SF3B1 mutant splicing alterations are corrected by GPATCH8 deletion, while ~60% of SF3B1 mutant splicing alterations are recapitulated upon SUGP1 deletion. GPATCH8 regulated splicing events have stronger branchpoints indicating a potential role for GPATCH8 in 3' splice site recognition. To further understand the connection between GPATCH8 and mutant SF3B1, we performed mass spectrometry studies of immunoprecipitated endogenous GPATCH8 in SF3B1 WT and mutant knockin cells. This revealed strong interaction between GPATCH8 and DHX15. GPATCH8 interacts with DHX15 at the same sites previously shown to be occupied by SUGP1 in crystal structures of the SUGP1/DHX15 interaction. These data and further biochemical studies elucidated that GPATCH8 competes with SUGP1 for interaction to DHX15. As such, deletion of GPATCH8 corrects SF3B1 mutant mis-splicing by enhancing the interaction of SUGP1 and DHX15 to the mutant SF3b complex. Given that GPATCH8 is required for a large proportion of the splicing alterations induced by mutant SF3B1, we next evaluated the phenotypic effect of correcting SF3B1 mutant mis-splicing. We first tested the impact of anti-Gpatch8 shRNAs in bone marrow from WT mice as well as animals with conditional knockin of Sf3b1 K700E, K666N, or R625H mutations. While each Sf3b1 mutation results in impaired colony formation of hematopoietic precursors in methylcellulose assays, silencing of Gpatch8 rescued colony formation from hematopoietic precursors of each of the different Sf3b1 mutant mice and was tolerated by normal hematopoietic precursors (Fig.B). Moreover, while CRISPR base-edited knockin of the SF3B1 K700E mutation in adult CD34 + cells impaired erythroid development, this erythroid differentiation defect was rescued by GPATCH8 deletion in the same cells (Fig.B). These data identify GPATCH8 as a novel RNA splicing factor involved in quality control of RNA branchpoint selection whose expression is required for mis-splicing by the different mutant forms of SF3B1. GPATCH8 antagonizes the activity of SUGP1 by competing for interaction with DHX15. These findings suggest that disrupting GPATCH8/DHX15 interaction could have important therapeutic benefit for the multitude of SF3B1 mutant hematopoietic diseases. This study also demonstrates the power of synthetic intronic splicing assays for discovery of trans factors and druggable proteins required by leukemia-associated mutant RNA splicing factors. Abdel-Wahab: AbbVie, Inc.: Consultancy, Membership on an entity's Board of Directors or advisory committees, Speakers Bureau; Loxo/Lilly: Consultancy; Nurix Therapeutics: Research Funding; Minovia Therapeutics: Research Funding; Amphista Therapeutics: Consultancy; AstraZeneca: Consultancy; Harmonic Discovery: Current holder of stock options in a privately-held company. Figure 1

Comprehensive genome analysis of Pseudomonas sp. SWRIQ11, a new plant growth-promoting bacterium that alleviates salinity stress in olive.

The study presents the genome analysis of a new Pseudomonas sp. (SWRIQ11), which can alleviate salinity stress effects on growth of olive seedlings in greenhouse study. The strain SWRIQ11 can tolerate salinity up to 6%, produce siderophores, indole acetic acid (IAA), aminocyclopropane-1-carboxylate (ACC) deaminase, and has the phosphate-solubilizing capability. The SWRIQ11 genome contained an assembly size of 6,196,390 bp with a GC content of 60.1%. According to derived indices based on whole-genome sequences for species delineation, including tetra nucleotide usage patterns (TETRA), genome-to-genome distance (GGDC), and average nucleotide identity (ANI), Pseudomonas sp. SWRIQ11 can be considered a novel species candidate. The phylogenetic analysis revealed SWRIQ11 clusters with Pseudomonas tehranensis SWRI196 in the same clade. The SWRIQ11 genome was rich in genes related to stress sensing, signaling, and response, chaperones, motility, attachments, colonization, and enzymes for degrading plant-derived carbohydrates. Furthermore, the genes for production of exopolysaccharides, osmoprotectants, phytohormones, and ACC deaminase, ion homeostasis, nutrient acquisition, and antioxidant defenses were identified in the SWRIQ11 genome. The results of genome analysis (identification of more than 825 CDSs related to plant growth-promoting and stress-alleviating traits in the SWRIQ11 genome which is more than 15% of its total CDSs) are in accordance with laboratory and greenhouse experiments assigning the Pseudomonas sp. SWRIQ11 as a halotolerant plant growth-promoting bacterium (PGPB). This research highlights the potential safe application of this new PGPB species in agriculture as a potent biofertilizer.

Consequences of asymmetric mutational pressure for the dynamic of linear B-cell epitopes repertoire of influenza a virus neuraminidase rearrangement

Full-length nucleotide sequences of avian influenza A virus neuraminidase coding region (20,631 sequences) were analyzed and compared with those isolated from viruses infecting human and swine (63,750 sequences). If in fourfold degenerate sites there is asymmetric A-bias that may be more or less asymmetric depending on the type of neuraminidase and the host, than in twofold degenerate sites from third codon positions there is a strong asymmetric U-bias in coding regions of N4, N5, and N8 isolated from viruses infecting birds, as well as in those of N1 and N2 isolated from viruses infecting human, swine, and birds, while in coding regions of N9 isolated from birds, there is surprisingly strong C-bias, and in sequences of N3, N6, and N7 the usage of C is quite close to the level of U. Revealed stabilization of both U and C in twofold degenerate sites is the evidence of frequent changes in mutational pressure direction. Asymmetric mutational pressure was one of the sources of amino acid replacements that resulted in an equal percentage of sites with appeared and disappeared linear B-cell epitopes in N1, N2, N4, and N5 (33.62–35.33% vs. 32.41–36.45%, respectively), and controlled by the immune pressure it resulted in a stronger tendency to disappear for B-cell epitopes of N3, N6, N7, N8, and N9 of avian viruses (8.74–28.77% vs. 28.96–38.89%). The lack of correlation between nucleotide usages in fourfold and twofold degenerate sites for three nucleotides, except U, is a strong evidence of mutational pressure theory.

Analysis for codon usage bias in membrane anchor of nonstructural protein 5A from BVDV.

The nonstructural protein 5A (NS5A) of the bovine viral diarrhea virus (BVDV) is a monotopic membrane protein. This protein can anchor to the cell membrane by an in-plane amphipathic ⍺-helix, which participates in the viral replication complex. In this study, the effects of synonymous codon usage pattern of NS5A and the overall transfer RNA (tRNA) abundance in cells on the formation of the in-plane membrane anchor of NS5A were analyzed, based on NS5A coding sequences of different BVDV genotypes. BVDV NS5A coding sequences represent the most potential for BVDV genotyping. Moreover, the nucleotide usage of BVDV NS5A dominates the genotype-specific pattern of synonymous codon usage. There is an obvious relationship between synonymous codon usage bias and the spatial conformation of the in-plane membrane anchor. Furthermore, the overall tRNA abundance profiling displays that codon positions with a high level of tRNA abundance are more than ones with a low level of tRNA abundance in the in-plane membrane anchor, implying that high translation speed probably acts on the spatial conformation of in-plane membrane anchor of BVDV NS5A. These results give a new opinion on the effect of codon usage bias in the formation of the in-plane membrane anchor of BVDV NS5A.

Distal nucleotides affect the rate of stop codon read‐through

BackgroundA key step in gene expression is the recognition of the stop codon to terminate translation at the correct position. However, it has been observed that ribosomes can misinterpret the stop codon and continue the translation in the 3′UTR region. This phenomenon is called stop codon read‐through (SCR). It has been suggested that these events would occur on a programmed basis, but the underlying mechanisms are still not well understood.MethodsHere, we present a strategy for the comprehensive identification of SCR events in the Drosophila melanogaster transcriptome by evaluating the ribosomal density profiles. The associated ribosomal leak rate was estimated for every event identified. A statistical characterization of the frequency of nucleotide use in the proximal region to the stop codon in the sequences associated to SCR events was performed.ResultsThe results show that the nucleotide usage pattern in transcripts with the UGA codon is different from the pattern for those transcripts ending in the UAA codon, suggesting the existence of at least two mechanisms that could alter the translational termination process. Furthermore, a linear regression models for each of the three stop codons was developed, and we show that the models using the nucleotides at informative positions outperforms those models that consider the entire sequence context to the stop codon.ConclusionsWe report that distal nucleotides can affect the SCR rate in a stop‐codon dependent manner.

Analysis of 3.5 million SARS-CoV-2 sequences reveals unique mutational trends with consistent nucleotide and codon frequencies

BackgroundSince the onset of the SARS-CoV-2 pandemic, bioinformatic analyses have been performed to understand the nucleotide and synonymous codon usage features and mutational patterns of the virus. However, comparatively few have attempted to perform such analyses on a considerably large cohort of viral genomes while organizing the plethora of available sequence data for a month-by-month analysis to observe changes over time. Here, we aimed to perform sequence composition and mutation analysis of SARS-CoV-2, separating sequences by gene, clade, and timepoints, and contrast the mutational profile of SARS-CoV-2 to other comparable RNA viruses.MethodsUsing a cleaned, filtered, and pre-aligned dataset of over 3.5 million sequences downloaded from the GISAID database, we computed nucleotide and codon usage statistics, including calculation of relative synonymous codon usage values. We then calculated codon adaptation index (CAI) changes and a nonsynonymous/synonymous mutation ratio (dN/dS) over time for our dataset. Finally, we compiled information on the types of mutations occurring for SARS-CoV-2 and other comparable RNA viruses, and generated heatmaps showing codon and nucleotide composition at high entropy positions along the Spike sequence.ResultsWe show that nucleotide and codon usage metrics remain relatively consistent over the 32-month span, though there are significant differences between clades within each gene at various timepoints. CAI and dN/dS values vary substantially between different timepoints and different genes, with Spike gene on average showing both the highest CAI and dN/dS values. Mutational analysis showed that SARS-CoV-2 Spike has a higher proportion of nonsynonymous mutations than analogous genes in other RNA viruses, with nonsynonymous mutations outnumbering synonymous ones by up to 20:1. However, at several specific positions, synonymous mutations were overwhelmingly predominant.ConclusionsOur multifaceted analysis covering both the composition and mutation signature of SARS-CoV-2 gives valuable insight into the nucleotide frequency and codon usage heterogeneity of SARS-CoV-2 over time, and its unique mutational profile compared to other RNA viruses.

Nucleotide and codon usage biases involved in the evolution of African swine fever virus: A comparative genomics analysis.

Since African swine fever virus (ASFV) replication is closely related to its host's machinery, codon usage of viral genome can be subject to selection pressures. A better understanding of codon usage can give new insights into viral evolution. We implemented information entropy and revealed that the nucleotide usage pattern of ASFV is significantly associated with viral isolation factors (region and time), especially the usages of thymine and cytosine. Despite the domination of adenine and thymine in the viral genome, we found that mutation pressure alters the overall codon usage pattern of ASFV, followed by selective forces from natural selection. Moreover, the nucleotide skew index at the gene level indicates that nucleotide usages influencing synonymous codon bias of ASFV are significantly correlated with viral protein hydropathy. Finally, evolutionary plasticity is proved to contribute to the weakness in synonymous codons with A- or T-end serving as optimal codons of ASFV, suggesting that fine-tuning translation selection plays a role in synonymous codon usages of ASFV for adapting host. Taken together, ASFV is subject to evolutionary dynamics on nucleotide selections and synonymous codon usage, and our detailed analysis offers deeper insights into the genetic characteristics of this newly emerging virus around the world.

Germline mutations directions are different between introns of the same gene: case study of the gene coding for amyloid-beta precursor protein.

Amyloid-beta precursor protein (APP) is highly conserved in mammals. This feature allowed us to compare nucleotide usage biases in fourfold degenerated sites along the length of its coding region for 146 species of mammals and birds in search of fragments with significant deviations. Even though cytosine usage has the highest value in fourfold degenerated sites in APP coding region from all tested placental mammals, in contrast to marsupial mammals with the bias toward thymine usage, the most frequent germline and somatic mutations in human APP coding region are C to T and G to A transitions. The same mutational AT-pressure is characteristic for germline mutations in introns of human APP gene. However, surprisingly, there are several exceptional introns with deviations in germline mutations rates. The most of those introns surround exons with exceptional biases in nucleotide usage in fourfold degenerated sites. Existence of such fragments in exons 4 and 5, as well as in exon 14, can be connected with the presence of lncRNA genes in complementary strand of DNA. Exceptional nucleotide usage bias in exons 16 and 17 that contain a sequence encoding amyloid-beta peptides can be explained either by the presence of yet unmapped lncRNA(s), or by the autonomous expression of a short mRNA that encodes just C-terminal part of the APP providing an alternative source of amyloid-beta peptides. This hypothesis is supported by the increased rate of T to C transitions in introns 16-17 and 17-18 of Human APP gene relatively to other introns.

Genetic features of bovine viral diarrhea virus subgenotype 1c in newborn calves at nucleotide and synonymous codon usages.

Bovine viral diarrhea virus (BVDV), serving as an important pathogen for newborn calves, poses threat to reproductive and economic losses in the cattle industry. To survey the infection rate and genetic diversity of BVDV in newborn calves in northern China, a total of 676 sera samples of newborn calves were collected from four provinces between 2021 and 2022. All sera samples were individually detected for BVDV infection by RT-PCR and ELISA. Our results showed that the overall serological rate was 9.76% (66/676) and the average positive rate of BVDV RNA was 8.14% (55/676) in the newborn calves. Eight BVDV strains were successfully isolated from RT-PCR positive sera samples, and four isolates displayed the cytopathic effect (CPE). Based on phylogenetic tree at the genome level, the eight strains were classified into subgenotype 1c. Moreover, the BVDV isolates had a close genetic relationship with the GSTZ strain at either nucleotide or codon usage level. Interestingly, in comparison of synonymous codon usage patterns between the BVDV isolates with CPE and ones without CPE, there were four synonymous codons (UCG, CCC, GCA, and AAC) which displayed the significant differences (p < 0.05) at codon usage pattern, suggesting that synonymous codon usage bias might play a role in BVDV-1c biotypes. In addition, the usage of synonymous codons containing CpG dinucleotides was suppressed by the BVDV-1c isolates, reflecting one of strategies of immune evasion of BVDV to its host. Taken together, our study provided data for monitoring and vaccination strategies of BVDV for newborn calves in northern China.

Evolutionary dynamics of codon usages for peste des petits ruminants virus.

Peste des petits ruminants virus (PPRV) is an important agent of contagious, acute and febrile viral diseases in small ruminants, while its evolutionary dynamics related to codon usage are still lacking. Herein, we adopted information entropy, the relative synonymous codon usage values and similarity indexes and codon adaptation index to analyze the viral genetic features for 45 available whole genomes of PPRV. Some universal, lineage-specific, and gene-specific genetic features presented by synonymous codon usages of the six genes of PPRV that encode N, P, M, F, H and L proteins reflected evolutionary plasticity and independence. The high adaptation of PPRV to hosts at codon usages reflected high viral gene expression, but some synonymous codons that are rare in the hosts were selected in high frequencies in the viral genes. Another obvious genetic feature was that the synonymous codons containing CpG dinucleotides had weak tendencies to be selected in viral genes. The synonymous codon usage patterns of PPRV isolated during 2007–2008 and 2013–2014 in China displayed independent evolutionary pathway, although the overall codon usage patterns of these PPRV strains matched the universal codon usage patterns of lineage IV. According to the interplay between nucleotide and synonymous codon usages of the six genes of PPRV, the evolutionary dynamics including mutation pressure and natural selection determined the viral survival and fitness to its host.

In Silico Analyses of the Role of Codon Usage at the Hemagglutinin Cleavage Site in Highly Pathogenic Avian Influenza Genesis.

A vast diversity of 16 influenza hemagglutinin (HA) subtypes are found in birds. Interestingly, viruses from only two subtypes, H5 and H7, have so far evolved into highly pathogenic avian influenza viruses (HPAIVs) following insertions or substitutions at the HA cleavage site by the viral polymerase. The mechanisms underlying this striking subtype specificity are still unknown. Here, we compiled a comprehensive dataset of 20,488 avian influenza virus HA sequences to investigate differences in nucleotide and amino acid usage at the HA cleavage site between subtypes and how these might impact the genesis of HPAIVs by polymerase stuttering and realignment. We found that sequences of the H5 and H7 subtypes stand out by their high purine content at the HA cleavage site. In addition, fewer substitutions were necessary in H5 and H7 HAs than in HAs from other subtypes to acquire an insertion-prone HA cleavage site sequence, as defined based on in vitro and in vivo data from the literature. Codon usage was more favorable for HPAIV genesis in sequences of viruses isolated from species or geographical regions in which HPAIV genesis is more frequently observed in nature. The results of the present analyses suggest that the subtype restriction of HPAIV genesis to H5 and H7 influenza viruses might be due to the particular codon usage at the HA cleavage site in these subtypes.

A Study on microRNAs Targeting the Genes Overexpressed in Lung Cancer and their Codon Usage Patterns

As reported by WHO in 2018, there were 2.09 million victims of lung cancer and 1.76 million fatalities worldwide. Tobacco remains the biggest hazard in causing this lethal disease. To execute the computational analysis, the overexpressed lung cancer genes were retrieved from literature and subsequently their complete coding sequences (CDS) were downloaded. The mature microRNA sequences of human were extracted from miRBASE. The 7mer-m8 perfect seed match between miRNAs and mRNAs was found. Following filtration, 7 genes were selected that possessed binding sites for maximum miRNAs viz., MUC5B (miR-4479, miR-1227-5p, miR-3940-5p, miR-604, miR-4455, miR-4267, miR-6750-3p, miR-4530, miR-5587-5p, miR-4508, miR-4534, miR-4443, miR-4253, miR-1321, miR-4655-5p, miR-4297, miR-4296, miR-1268a, miR-3178, miR-4750-3p, miR-1306-3p, miR-1268b, miR-3656, miR-1233-3p, miR-6804-5p), MUC16 (miR-4456, miR-1205, miR-665, miR-6808-3p, miR-1279, miR-4257, miR-1227-5p, miR-888-3p, miR-4455, miR-4267, miR-4294, miR-1275, miR-4288, miR-1178-5p, miR-4314, miR-6829-3p, miR-548av-5p, miR-1294, miR-5587-5p, miR-3622b-5p, miR-1273f, miR-4770, miR-4327, miR-4318, miR-4531, miR-4534, miR-4443, miR-7106-5p, miR-3125, miR-3650, miR-4325, miR-4266, miR-7976, miR-1290, miR-4500, miR-7160-5p, miR-4291, miR-1306-3p, miR-6130, miR-4430, miR-4725-5p, miR-4441, miR-6077, miR-1304-5p, miR-7515, miR-3182, miR-6134), COL1A1 (miR-3665, miR-1227-5p, miR-6132, miR-2861, miR-4530, miR-3155b, miR-3155a, miR-1292-3p, miR-4497), COL5A1 (miR-7162-5p, miR-3665, miR-6809-3p, miR-4313, miR-4531, miR-4532, miR-3155b, miR-4323, miR-1207-3p, miR-4260, miR-6071, miR-4710, miR-7162-5p), CELSR2 (miR-7150, miR-4308, miR-6132, miR-4770, miR-4534, miR-4492, miR-3960, miR-3178, miR-4291, miR-563), COL7A1 (miR-665, miR-6730-3p, miR-1227-5p, miR-4265, miR-6829-3p, miR-4297, miR-4532, miR-3181, miR-4310, miR-4441, miR-4497, miR-1237-3p), and FAT2 (miR-4267, miR-1275, miR-4770, miR-1825, miR-6895-5p, miR-4535, miR-4493, miR-940, miR-6861-3p, miR-4310, miR-4710, miR-4447, miR-4472). The miRNA-target site and their flank regions were compared with respect to site accessibility, translational rate, and relationship between RSCU and tRNAs. Higher accessibilities to miRNA-binding regions and lower translational rates indicated that miRNAs' binding to their respective targets might be efficient. The presence of rare codons might further augment miRNA targeting. The codon usage bias study of the genes related to lung cancer revealed non-uniform usage of nucleotides and comparatively higher GC content. Lower biasness prevailed in the genes and selective constraint mostly governed them. Lastly, the functionalities of target genes were also revealed. The silencing characteristic of miRNAs might be exploited to design miRNA-mediated therapy that might potentially repress the overexpressed genes in carcinoma.

Characteristics of the nucleotide composition and de novo structure of 28S rRNA gene of intestinal flukes (Echinostomatidae) and genetic distance among species (suborder Echinostomata)

The family Echinostomatidae includes four genera, Echinostoma, Hypoderaeum, Echinoparyphium, and Artyfechinostomum, which are of public epidemiological and pathogenic importance. The morphological similarity of species requires differentiating methods, and the use of molecular genetic markers is mainly encountered. The 28S rRNA gene sequence in the family Echinostomatidae is an important marker for the analysis of the interspecies and phylogenetic relationships. However, information about the structure and composition of the complete 28S rRNA gene in Echinostomatidae is limited. In this study, the complete 28S rRNA sequences (3,881–3,883 bp) of four species, including Echinostoma revolutum (designated Erev-MSD15-TH), E. malayanum (Emal-EMI3-TH), E. miyagawai (Emiy-RED11-TH), and Hypoderaeum conoideum (Hcon-RED42-TH) were obtained and analyzed for their genetic characteristics. All had complex secondary structures, containing numerous branches, “hairpins”, and “loops” that provide stability of the rRNA gene. The de novo structural pattern and the "skew" value in the nucleotide use of E. revolutum, E. malayanum, E. miyagawai, and H. conoideum were analyzed. The nucleotide usage of A, T, G, C and A+T and G+C of 28S rRNA in all species of the family Echinostomatidae was similar, in which A and T usage were almost equal, and G and C were slightly different in the range of 7–8%, so the skewness for AT was very low negative value (–0.029 to –0.045), and for GC was not very high (0.163 to 0.172). The genetic distance estimated among strains within a species (intra-specific level) was very low in both the Echinostomatidae (0.05%) and the Fasciolidae families (0–0.03%), and that among different species (inter-specific level) was higher. However, among species in different families (inter-familial level) the genetic distance was highest, at 3.80–4.86%, as seen between Echinostoma and Fasciola species. The 28S rRNA gene dataset obtained from four species of the intestinal flukes (family Echinostomatidae) in this study contribute to the understandings of genetic markers in the suborder Echinostomata for use in diagnostic, phylogenetic, taxonomic, epidemiological, and population genetics studies.

Bias at the third nucleotide of codon pairs in virus and host genomes

Genomes of different sizes and complexity can be compared using common features. Most genomes contain open reading frames, and most genomes use the same genetic code. Redundancy in the genetic code means that different biases in the third nucleotide position of a codon exist in different genomes. However, the nucleotide composition of viruses can be quite different from host nucleotide composition making it difficult to assess the relevance of these biases. Here we show that grouping codons of a codon-pair according to the GC content of the first two nucleotide positions of each codon reveals patterns in nucleotide usage at the third position of the 1st codon. Differences between the observed and expected biases occur predominantly when the first two nucleotides of the 2nd codon are both S (strong, G or C) or both W (weak, A or T), not a mixture of strong and weak. The data indicates that some codon pairs are preferred because of the strength of the interactions between the codon and anticodon, the adjacent tRNAs and the ribosome. Using base-pairing strength and third position bias facilitates the comparison of genomes of different size and nucleotide composition and reveals patterns not previously described.

Nucleotide Usage Biases Distort Inferences of the Species Tree.

Despite the importance of natural selection in species’ evolutionary history, phylogenetic methods that take into account population-level processes typically ignore selection. The assumption of neutrality is often based on the idea that selection occurs at a minority of loci in the genome and is unlikely to compromise phylogenetic inferences significantly. However, genome-wide processes like GC-bias and some variation segregating at the coding regions are known to evolve in the nearly neutral range. As we are now using genome-wide data to estimate species trees, it is natural to ask whether weak but pervasive selection is likely to blur species tree inferences. We developed a polymorphism-aware phylogenetic model tailored for measuring signatures of nucleotide usage biases to test the impact of selection in the species tree. Our analyses indicate that although the inferred relationships among species are not significantly compromised, the genetic distances are systematically underestimated in a node-height-dependent manner: that is, the deeper nodes tend to be more underestimated than the shallow ones. Such biases have implications for molecular dating. We dated the evolutionary history of 30 worldwide fruit fly populations, and we found signatures of GC-bias considerably affecting the estimated divergence times (up to 23%) in the neutral model. Our findings call for the need to account for selection when quantifying divergence or dating species evolution.

Compositional features and pattern of codon usage for mitochondrial CO genes among reptiles

The phenomenon of non-random occurrence of synonymous nucleotide triplets (codons) in the coding sequences of genes is the codon usage bias (CUB). In this study, we used bioinformatic tool kit to analyze the compositional pattern and CUB of mitogenes namely COI, COII and COIII across different orders of reptiles. Estimation of overall base composition in the protein-coding sequences of COI, COII and COIII genes of the reptilian orders revealed an uneven usage of nucleotides. The overall count of A nucleotide was found to be the highest while the overall count of G nucleotide was the least. The CO genes across the three reptilian orders were prominently AT biased. Comparison of the GC proportion at each codon position displayed that GC1 percentage ranked the highest in all the three CO genes of the reptilian orders. SCUO values indicated weaker CUB, while considerable variation of SCUO values existed in the three CO genes across the studied reptiles. Relative synonymous codon usage (RSCU) values indicated that mostly the A ending codons were preferred. Based on the parameters namely neutrality plot, mutational responsive index and translational selection, we could conclude that natural selection was the major evolutionary force in COI, COII and COIII genes in the studied reptilian orders. However, correspondence analysis, parity plot and correlation studies indicated the existence of mutation pressure as well on the CO genes.

Codon Usage Bias in Autophagy-Related Gene 13 in Eukaryotes: Uncovering the Genetic Divergence by the Interplay Between Nucleotides and Codon Usages.

Synonymous codon usage bias is a universal characteristic of genomes across various organisms. Autophagy-related gene 13 (atg13) is one essential gene for autophagy initiation, yet the evolutionary trends of the atg13 gene at the usages of nucleotide and synonymous codon remains unexplored. According to phylogenetic analyses for the atg13 gene of 226 eukaryotic organisms at the nucleotide and amino acid levels, it is clear that their nucleotide usages exhibit more genetic information than their amino acid usages. Specifically, the overall nucleotide usage bias quantified by information entropy reflected that the usage biases at the first and second codon positions were stronger than those at the third position of the atg13 genes. Furthermore, the bias level of nucleotide ‘G’ usage is highest, while that of nucleotide ‘C’ usage is lowest in the atg13 genes. On top of that, genetic features represented by synonymous codon usage exhibits a species-specific pattern on the evolution of the atg13 genes to some extent. Interestingly, the codon usages of atg13 genes in the ancestor animals (Latimeria chalumnae, Petromyzon marinus, and Rhinatrema bivittatum) are strongly influenced by mutation pressure from nucleotide composition constraint. However, the distributions of nucleotide composition at different codon positions in the atg13 gene display that natural selection still dominates atg13 codon usages during organisms’ evolution.