Genomic analysis and phylogenetic characterization of Himalayan snow trout, Schizothorax esocinus based on mitochondrial protein-coding genes.

BACKGROUND: Mitochondrial DNA (mtDNA) has become a significant tool for exploring genetic diversity and delineating evolutionary links across diverse taxa. Within the group of cold-water fish species that are native to the Indian Himalayan region, Schizothorax esocinus holds particular importance due to its ecological significance and is potentially vulnerable to environmental changes. This research aims to clarify the phylogenetic relationships within the Schizothorax genus by utilizing mitochondrial protein-coding genes. METHODS: Standard protocols were followed for the isolation of DNA from S. esocinus. For the amplification of mtDNA, overlapping primers were used, and then subsequent sequencing was performed. The genetic features were investigated by the application of bioinformatic approaches. These approaches covered the evaluation of nucleotide composition, codon usage, selective pressure using nonsynonymous substitution /synonymous substitution (Ka/Ks) ratios, and phylogenetic analysis. RESULTS: The study specifically examined the 13 protein-coding genes of Schizothorax species which belongs to the Schizothoracinae subfamily. Nucleotide composition analysis showed a bias towards A + T content, consistent with other cyprinid fish species, suggesting evolutionary conservation. Relative Synonymous Codon Usage highlighted leucine as the most frequent (5.18%) and cysteine as the least frequent (0.78%) codon. The positive AT-skew and the predominantly negative GC-skew indicated the abundance of A and C. Comparative analysis revealed significant conservation of amino acids in multiple genes. The majority of amino acids were hydrophobic rather than polar. The purifying selection was revealed by the genetic distance and Ka/Ks ratios. Phylogenetic study revealed a significant genetic divergence between S. esocinus and other Schizothorax species with interspecific K2P distances ranging from 0.00 to 8.87%, with an average of 5.76%. CONCLUSION: The present study provides significant contributions to the understanding of mitochondrial genome diversity and genetic evolution mechanisms in Schizothoracinae, hence offering vital insights for the development of conservation initiatives aimed at protecting freshwater fish species.

Comparative Mitogenomic Analyses of Darkling Beetles (Coleoptera: Tenebrionidae) Provide Evolutionary Insights into tRNA-like Sequences.

Tenebrionidae is widely recognized owing to its species diversity and economic importance. Here, we determined the mitochondrial genomes (mitogenomes) of three Tenebrionidae species (Melanesthes exilidentata, Anatolica potanini, and Myladina unguiculina) and performed a comparative mitogenomic analysis to characterize the evolutionary characteristics of the family. The tenebrionid mitogenomes were highly conserved with respect to genome size, gene arrangement, base composition, and codon usage. All protein-coding genes evolved under purifying selection. The largest non-coding region (i.e., control region) showed several unusual features, including several conserved repetitive fragments (e.g., A+T-rich regions, G+C-rich regions, Poly-T tracts, TATA repeat units, and longer repetitive fragments) and tRNA-like structures. These tRNA-like structures can bind to the appropriate anticodon to form a cloverleaf structure, although base-pairing is not complete. We summarized the quantity, types, and conservation of tRNA-like sequences and performed functional and evolutionary analyses of tRNA-like sequences with various anticodons. Phylogenetic analyses based on three mitogenomic datasets and two tree inference methods largely supported the monophyly of each of the three subfamilies (Stenochiinae, Pimeliinae, and Lagriinae), whereas both Tenebrioninae and Diaperinae were consistently recovered as polyphyletic. We obtained a tenebrionid mitogenomic phylogeny: (Lagriinae, (Pimeliinae, ((Tenebrioninae + Diaperinae), Stenochiinae))). Our results provide insights into the evolution and function of tRNA-like sequences in tenebrionid mitogenomes and contribute to our general understanding of the evolution of Tenebrionidae.

Analysis of the Compositional Features and Codon Usage Pattern of Genes Involved in Human Autophagy.

Autophagy plays an intricate role in paradigmatic human pathologies such as cancer, and neurodegenerative, cardiovascular, and autoimmune disorders. Autophagy regulation is performed by a set of autophagy-related (ATG) genes, first recognized in yeast genome and subsequently identified in other species, including humans. Several other genes have been identified to be involved in the process of autophagy either directly or indirectly. Studying the codon usage bias (CUB) of genes is crucial for understanding their genome biology and molecular evolution. Here, we examined the usage pattern of nucleotide and synonymous codons and the influence of evolutionary forces in genes involved in human autophagy. The coding sequences (CDS) of the protein coding human autophagy genes were retrieved from the NCBI nucleotide database and analyzed using various web tools and software to understand their nucleotide composition and codon usage pattern. The effective number of codons (ENC) in all genes involved in human autophagy ranges between 33.26 and 54.6 with a mean value of 45.05, indicating an overall low CUB. The nucleotide composition analysis of the autophagy genes revealed that the genes were marginally rich in GC content that significantly influenced the codon usage pattern. The relative synonymous codon usage (RSCU) revealed 3 over-represented and 10 under-represented codons. Both natural selection and mutational pressure were the key forces influencing the codon usage pattern of the genes involved in human autophagy.

Fowl adenovirus strains 1/A and 11/D isolated from birds with reovirus infection.

Inclusion body hepatitis (IBH) is, in some cases, a fatal disease affecting fowl by adenovirus strains which are subdivided into 5 species (A-E). In the current study, we investigated sequences from the Loop L1 region of the hexon gene of sequences of adenovirus field stains 1/A and 11/D isolated from a poultry flock co-infected with IBH and avian reoviruses ARVs. In early 2021, an epidemiologic survey highlighted the coinfection adenoviruses with other viruses (orthoreovirus infection) as being particularly deleterious within the poultry industry. Here, we investigated the Loop L1 HVR1-4 region of the hexon gene with relative synonymous codon usage (RSCU) designation and RSCU inclusive of all the mutations. These are the first results that have been presented on fowl adenovirus species A and D with simultaneous reovirus infection in 38-days old broiler chickens in Poland.

Relative synonymous codon usage of ORF1ab in SARS-CoV-2 and SARS-CoV.

BACKGROUND: COVID-19, as a novel coronavirus disease caused by new coronavirus SARS-CoV-2, spreads all over the world, and brings harm to human in many countries. Humans suffered a lot from both SARS-CoV-2 now and by SARS-CoV in the year 2003. It is important to understand the differences and the relationships between these two types of viruses. OBJECTIVE: To compare relative synonymous codon usage of ORF1ab gene in SARS-CoV-2 and SARS-CoV, relative synonymous codon usage of their genomes are studied in this paper from the bioinformatics perspective. METHODS: The ORF1ab gene, which is an important non-structural polyprotein coding gene and now used for nucleic acid detection markers in many measurement method, in both SARS-CoV-2 (30 strains) and SARS-CoV (20 strains) are considered to be the research object in the present paper. The relative synonymous codon usage values of the ORF1ab gene are calculated to characterize the differences and the evolutionary characteristics among 50 strains. RESULTS: There is a significant difference between SARS-CoV and SARS-CoV-2 when the relative synonymous codon usage value of ORF1ab genes is concerned. The results suggest that codon usage pattern of SARS-CoV is more similar to human than that of the SARS-CoV-2, and that the inner difference in SARS-CoV-2 strains is larger than that of SARS-CoV, which denote the larger diversity exits in the SARS-CoV-2 virus. CONCLUSION: These results show that the relative synonymous codon usage values in the coronavirus could be used for further research on their evolutionary phenomenon.

Translational adaptation of human viruses to the tissues they infect.

Viruses need to hijack the translational machinery of the host cell for a productive infection to happen. However, given the dynamic landscape of tRNA pools among tissues, it is unclear whether different viruses infecting different tissues have adapted their codon usage toward their tropism. Here, we collect the coding sequences of 502 human-infecting viruses and determine that tropism explains changes in codon usage. Using the tRNA abundances across 23 human tissues from The Cancer Genome Atlas (TCGA), we build an in silico model of translational efficiency that validates the correspondence of the viral codon usage with the translational machinery of their tropism. For instance, we detect that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is specifically adapted to the upper respiratory tract and alveoli. Furthermore, this correspondence is specifically defined in early viral proteins. The observed tissue-specific translational efficiency could be useful for the development of antiviral therapies and vaccines.

Leucine encoding codon TTG shows an inverse relationship with GC content in genes involved in neurodegeneration with iron accumulation.

We determined various forces involved in shaping codon usage of the genes linked to brain iron accumulation and infantile neuroaxonal dystrophy. The analysis paved the way for determining the forces responsible for composition, expression level, physical properties and codon bias of a gene. An interesting observation related to composition was that, on all the three codon positions, any two of the four nucleotides had similar compositions. CpG, TpA, and GpT dinucleotides were underrepresented with the overrepresentation of TpG dinucleotide. CpG and TpA containing codons ATA, CTA, TCG, and GCG were underrepresented, while TpG dinucleotide containing codon CTG was overrepresented, indicative of compositional constraints importance. GC ending codons were favored when the genome is GC rich, except leucine encoding codon TTG, which exhibits an inverse relationship with GC content. Nucleotide disproportions are found associated with the physical properties of proteins. The values of CAI and ENc are suggestive of low codon bias in genes. Considering the results of neutrality analysis, parity analysis, underrepresentation of TpA and CpG codons, and over-representation of TpG codons, the correlation between the compositional constraints and skew relationships with protein properties suggested the role of all the three selectional, mutational and compositional forces in shaping codon usage with the dominance of selectional pressure.

Codon bias analyses on thyroid carcinoma genes.

BACKGROUND: Thyroid carcinoma is one of the most common cancers in the world. Although the genetics of thyroid carcinoma was intensively studied, new mechanisms could be involved in its development as the codon bias. In this paper, we studied the codon bias of thyroid-cancer genes, considering not only the sequences but also the synonymous mutations. METHODS: Different measures and statistical analyses were employed to characterize the thyroid-cancer genes. We considered classical measures as RSCU and ENC, the compositional and protein characteristics, but also the codon bias landscape via the %MinMax algorithm. RESULTS: The compositional analyses highlighted two groups of thyroid cancer genes according to the GC% and GC3% content. The ENC did not show a clear codon bias in the genes. Differently, the RSCU analyses showed interesting codons that could play an important role in the development of thyroid cancer as the codon Ser-tcG. Furthermore, interesting synonymous mutations were detected that could affect the codon bias. The codon bias landscape detected genes enriched in rare codons as AKAP9 and KTN1. A cluster analysis based on %MinMax classified the thyroid cancer genes in four different groups according to the distribution of rare/frequent codons in the sequence. CONCLUSIONS: This is the first study that analyzed the codon bias in thyroid cancer genes based also on synonymous mutations. This study provided different hints that should be further investigated by wet-lab validation and that it could open new scenarios in the understanding the molecular mechanisms involved in thyroid cancer development based on codon bias.

Natural Selection Shapes Codon Usage in the Human Genome.

Synonymous codon usage has been identified as a determinant of translational efficiency and mRNA stability in model organisms and human cell lines. However, whether natural selection shapes human codon content to optimize translation efficiency is unclear. Furthermore, aside from those that affect splicing, synonymous mutations are typically ignored as potential contributors to disease. Using genetic sequencing data from nearly 200,000 individuals, we uncover clear evidence that natural selection optimizes codon content in the human genome. In deriving intolerance metrics to quantify gene-level constraint on synonymous variation, we discover that dosage-sensitive genes, DNA-damage-response genes, and cell-cycle-regulated genes are particularly intolerant to synonymous variation. Notably, we illustrate that reductions in codon optimality in BRCA1 can attenuate its function. Our results reveal that synonymous mutations most likely play an underappreciated role in human variation.

Comprehensive analysis of the codon usage patterns of polyprotein of Zika virus.

Zika virus (ZIKV) is a mosquito-borne virus in the family Flaviviridae, and the massive outbreak of ZIKV has endangered public health. Codon usage patterns of viruses reflect a series of evolutionary changes that enable viruses to shape their survival rates and fitness toward the external environment and, most importantly, their hosts. In this study, 90 ZIKV isolates were used for a comprehensive analysis on the codon usage patterns. The overall codon usage among ZIKV strains is similar and slightly biased. The value of effective number of codons (ENC) showed that the overall extent of codon usage bias in ZIKV is relatively low. Nucleotide analysis showed that the overall codon usage is biased toward A- and G-ending codons. The phylogenetic analysis indicated that their independent evolutionary origins from a common ancestor. The RSCU analysis showed that the codon usage pattern of ZIKV is more similar to that of Homo sapiens. Correlation analysis, Correspondence analysis, ENC-GC3S plot, and PR2 plot indicated that the codon usage patterns of the viruses are not only influenced by mutational pressure but also by natural selection, but neutrality plot analysis showed that the latter plays a major role. These results built the base for further research on the molecular evolution of ZIKV.

Codon usage of human hepatitis C virus clearance genes in relation to its expression.

Hepatitis C virus (HCV) infection is among the leading causes of hepatocellular carcinoma and liver cirrhosis globally, with a high economic burden. The disease progression is well established, but less is known about the spontaneous HCV infection clearance. This study tries to establish the relationship between codon biasness and expression of HCV clearance candidate genes in normal and HCV infected liver tissues. A total of 112 coding sequences comprising 151 679 codons were subjected to the computation of codon indices, namely relative synonymous codon usage, an effective number of codon (Nc), frequency of optimal codon, codon adaptation index, codon bias index, and base compositions. Codon indices report of GC3s, GC12, hydropathicity, and aromaticity implicates both mutational and translational selection in the candidate gene set. This was further correlated with the differentially expressed genes among the selected genes using BioGPS. A significant correlation is observed between the gene expression of normal liver and cancerous liver tissues with codon bias (Nc). Gene expression is also correlated with relative codon bias values, indicating that CCL5, APOA2, CD28, IFITM1, and TNFSF4 genes have higher expression. These results are quite encouraging in selecting the high responsive genes in HCV clearance. However, there could be additional genes which could also orchestrate the clearance role with the above mentioned first line of defensive genes.

Is there codon usage bias for poly-Q stretches in the human proteome?

We have analyzed codon usage for poly-Q stretches of different lengths for the human proteome. First, we have obtained that all long poly-Q stretches in Protein Data Bank (PDB) belong to the disordered regions. Second, we have found the bias for codon usage for glutamine homo-repeats in the human proteome. In the cases when the same codon is used for poly-Q stretches only CAG triplets are found. Similar results are obtained for human proteins with glutamine homo-repeats associated with diseases. Moreover, for proteins associated with diseases (from the HraDis database), the fraction of proteins for which the same codon is used for glutamine homo-repeats is less (22%) than for proteins from the human proteome (26%). We have demonstrated for poly-Q stretches in the human proteome that in some cases (28) the splicing sites correspond to the homo-repeats and in 11 cases, these sites appear at the C -terminal part of the homo-repeats with statistical significance 10 -8 .