Codon usage bias and phylogenetic analysis of mitochondrial ND1 gene in pisces, aves, and mammals

Five new mitogenomes sequences of Calidridine sandpipers (Aves: Charadriiformes) and comparative mitogenomics of genus Calidris

Background

The genus Calidris (Charadriiformes, Scolopacidae) includes shorebirds known as dunlin, knots, and sanderlings. The relationships between species nested within Calidris, including Eurynorynchus, Limicola and Aphriza, are not well-resolved.

Methods

Samples were collected from Xiaoyangkou, Rudong County, Jiangsu Province, China. Mitogenomes were sequenced using the Illumina Novaseq 6000 platform for PE 2 × 150 bp sequencing, and then checked for PCR products. Protein-coding genes were determined using an Open Reading Frame Finder. tRNAscan-SE, MITOS, and ARWEN were used to confirm tRNA and rRNA annotations. Bioinformatic analyses were conducted using DnaSP 5.1 and MEGA X. Phylogenic trees were constructed using maximum likelihood and Bayesian analyses.

Results

We sequenced and annotated the mitogenome of five species and obtained four complete mitogenomes and one nearly complete mitogenome. Circular mitogenomes displayed moderate size variation, with a mean length of 16,747 bp, ranging from 16,642 to 16,791 bp. The mitogenome encoded a control region and a typical set of 37 genes containing two rRNA genes, 13 protein-coding genes, and 22 tRNA genes. There were four start codons, four stop codons, and one incomplete stop codon (T-). The nucleotide composition was consistently AT-biased. The average uncorrected pairwise distances revealed heterogeneity in the evolutionary rate for each gene; the COIII had a slow evolutionary rate, whereas the ATP8 gene had a fast rate. dN/dS analysis indicated that the protein-coding genes were under purifying selection. The genetic distances between species showed that the greatest genetic distance was between Eurynorhynchus pygmeus and Limicola falcinellus (22.5%), and the shortest was between E. pygmeus and Calidris ruficollis (12.8%). Phylogenetic trees revealed that Calidris is not a monophyletic genus, as species from the genera Eurynorynchus and Limicola were nested within Calidris. The molecular data obtained in this study are valuable for research on the taxonomy, population genetics, and evolution of birds in the genus Calidris.

Characterisation of the Complete Mitochondrial Genome of Critically Endangered Mustela lutreola (Carnivora: Mustelidae) and Its Phylogenetic and Conservation Implications

In this paper, a complete mitochondrial genome of the critically endangered European mink Mustela lutreola L., 1761 is reported. The mitogenome was 16,504 bp in length and encoded the typical 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes, and harboured a putative control region. The A+T content of the entire genome was 60.06% (A > T > C > G), and the AT-skew and GC-skew were 0.093 and −0.308, respectively. The encoding-strand identity of genes and their order were consistent with a collinear gene order characteristic for vertebrate mitogenomes. The start codons of all protein-coding genes were the typical ATN. In eight cases, they were ended by complete stop codons, while five had incomplete termination codons (TA or T). All tRNAs had a typical cloverleaf secondary structure, except tRNA^Ser(AGC) and tRNA^Lys, which lacked the DHU stem and had reduced DHU loop, respectively. Both rRNAs were capable of folding into complex secondary structures, containing unmatched base pairs. Eighty-one single nucleotide variants (substitutions and indels) were identified. Comparative interspecies analyses confirmed the close phylogenetic relationship of the European mink to the so-called ferret group, clustering the European polecat, the steppe polecat and the black-footed ferret. The obtained results are expected to provide useful molecular data, informing and supporting effective conservation measures to save M. lutreola.

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.

Analysis of codon usage bias in mitochondrial CO gene among platyhelminthes

The phenomenon of non-uniform usage of the synonymous codons, where some codons are given more preference to others, is known as codon usage bias (CUB). CUB is known to be determined by two major evolutionary forces i.e. mutation pressure and selection. We used various approaches to understand the codon usage pattern in mitochondrial CO (MT-CO) genes involved in complex IV of the respiratory chain (RC) as no work was reported yet. Our present study revealed that CUB was relatively high and the coding sequences were rich in A and T. Correspondence analysis further indicated that A/T compositional properties under mutational pressure might be affecting the codon usage pattern and was different in different classes for MT-CO gene. A highly significant correlation between A% and A3%, T% and T3%, G% and G3%, C% and C3%, GC% and GC3% in all the classes indicated that compositional constraints under mutational pressure and natural selection might affect the CUB. Neutrality plot indicated that both natural selection and mutational bias affected the CUB, where, natural selection played the major role as compared to mutational pressure.

Understanding the codon usage patterns of mitochondrial CO genes among Amphibians

A universal phenomenon of using synonymous codons unequally in coding sequences known as codon usage bias (CUB) is observed in all forms of life. Mutation and natural selection drive CUB in many species but the relative role of evolutionary forces varies across species, genes and genomes. We studied the CUB in mitochondrial (mt) CO genes from three orders of Amphibia using bioinformatics approach as no work was reported yet. We observed that CUB of mt CO genes of Amphibians was weak across different orders. Order Caudata had higher CUB followed by Gymnophiona and Anura for all genes and CUB also varied across genes. Nucleotide composition analysis showed that CO genes were AT-rich. The AT content in Caudata was higher than that in Gymnophiona while Anura showed the least content. Multiple investigations namely nucleotide composition, correspondence analysis, parity plot analysis showed that the interplay of mutation pressure and natural selection caused CUB in these genes. Neutrality plot suggested the involvement of natural selection was more than the mutation pressure. The contribution of natural selection was higher in Anura than Gymnophiona and the lowest in Caudata. The codons CGA, TGA, AAA were found to be highly favoured by nature across all genes and orders.

Codon usage pattern and its influencing factors for mitochondrial CO genes among different classes of Arthropoda

Analysis of codon usage bias (CUB) is very much important in perceiving the knowledge of molecular biology, the discovery of a new gene, designing of transgenes and evolution of gene. In this study, we analyzed compositional features and codon usage of MT-CO (COI, COII and COIII) genes among the classes of Arthropoda to explore the pattern of CUB as no research work was reported yet. Nucleotide composition analysis in CO genes suggested that the genes were AT-rich in all the four classes of Arthropoda. CUB was low in all the classes of Arthropoda for MT-CO genes as revealed from a high effective number of codons (ENC). We also found that the evolutionary forces namely mutation pressure and natural selection were the key influencing factors in CUB among MT-CO genes as revealed by correlation analysis between overall nucleotide composition and nucleotide composition at the 3rd codon position. Correspondence analysis suggested that the pattern of CUB was different among the classes of Arthropoda. Further, it was revealed from the neutrality plot that natural selection had a dominant role while mutation pressure exhibited a minor role in structuring the pattern of codon usage in all the classes of Arthropoda across COI, COII and COIII genes.

Comparative analysis of codon usage between Gossypium hirsutum and G. barbadense mitochondrial genomes

Gossypium hirsutum and G. barbadense mitochondrial genomes were analyzed to understand the factors shaping codon usage. While most analyses of codon usage suggest minimal to no bias, nucleotide composition, specifically GC content, was significantly correlated with codon usage. In general, both mitochondrial genomes favor codons that end in A or U, with a secondary preference for pyrimidine rich codons. These observations are similar to previous reports of codon usage in cotton nuclear genomes, possibly suggestive of a general bias spanning genomic compartment. Although evidence for codon usage bias is weak for most genes, we identified six genes (i.e. atp8, atp9, sdh3, sdh4, mttB and rpl2) with significant nonrandom codon usage. In general, we find multiple factors that influence cotton mitochondrial genome codon usage, which may include selection in a subset of genes.

Compositional properties and codon usage pattern of mitochondrial ATP gene in different classes of Arthropoda

Codon usage bias (CUB) is defined as the usage of synonymous codons unequally for an amino acid in a gene transcript. It is influenced by both mutation pressure and natural selection and is a species-specific property. In our current study, we used bioinformatic methods to investigate the coding sequences of mitochondrial adenosine triphosphate gene (MT-ATP) in different classes of arthropoda to know the codon usage pattern of the gene as no work was described earlier. The analysis of compositional properties suggested that the gene is AT rich. The effective number of codons revealed the CUB of both ATP6 and ATP8 gene was moderate. Heat map showed that the codons ending with AT were negatively associated with GC3 while the codons ending with GC were positively associated with GC3 in all the classes of arthropoda. Correspondence study revealed that the pattern of codon usage of ATP6 and ATP8 genes differed across classes. Neutrality plot suggested the codon usage bias of these two genes in phylum arthropoda was influenced by both mutation pressure and natural selection.

A cross-talk on compositional dynamics and codon usage patterns of mitochondrial CYB gene in Echinodermata

Codon usage bias (CUB) refers to a phenomenon in which some synonymous codons are used in mature mRNA at a higher frequency than other members codifying the same amino acid. CUB is mainly determined by mutation pressure and natural selection. We used bioinformatic tools to analyze the protein coding sequences of mitochondrial CYB gene in different classes of Echinodermata to understand the patterns of codon usage. The ENC values of CYB gene in five different classes of Echinodermata were 41.64, 30.33, 43.63, 41.11, and 41.33, which suggested that the CUB of this gene was low. The relative synonymous codon usage (RSCU) values showed that the patterns of over-represented and under-represented codons were different among different classes. Correspondence analysis indicated that the plots of CYB gene were different across classes, suggesting that the pattern of codon usage was also different among five classes under study. Highly significant correlation (p < .01) between overall nucleotide composition and its 3rd codon position indicated that both mutational pressure and natural selection had an influence on the codon usage bias of CYB gene. Furthermore, PR-2 bias plot analysis showed that both mutation pressure and natural selection might have affected the pattern of codon usage in CYB gene of Echinodermata.

The mitochondrial genome of red-necked phalarope Phalaropus lobatus (Charadriiformes: Scolopacidae) and phylogeny analysis among Scolopacidae

The red-necked phalarope is a wonderful species with specific morphological characters and lifestyles. Mitochondrial genomes, encoding necessary proteins involved in the system of energy metabolism, are important for the evolution and adaption of species. In this study, we determined the complete mitogenome sequence of Phalaropus lobatus (Charadriiformes: Scolopacidae). The circular genome is 16714 bp in size, containing 13 PCGs, two ribosomal RNAs and 22 tRNAs and a high AT-rich control region. The AT skew and GC skew of major strand is positive and negative respectively. Most of PCGs are biased towards A-rich except ND1. A codon usage analysis shows that 3 start codons (ATG, GTG and ATA), 4 stop codons (TAA, TAG, AGG, AGA) and two incomplete terminate codons (T-). Twenty two transfer RNAs have the typical cloverleaf structure, and a total of ten base pairs are mismatched throughout the nine tRNA genes. The phylogenetic tree based on 13 PCGs and 2 rRNA genes indicates that monophyly of the family and genus Phalaropus is close to genus Xenus plus Tringa. The analysis of selective pressure shows 13 protein-coding genes are evolving under the purifying selection and P. lobatus is different from other Scolopacidae species on the selective pressure of gene ND4. This study helps us know the inherent mechanism of mitochondrial structure and natural selection.