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.

Gene Rearrangements in the Mitochondrial Genomes of Whiteflies (Hemiptera: Aleyrodinae): Plesiomorphies, Synapomorphies and Autapomorphies

Mitochondrial genome rearrangements have been used for defining historical relationships, but there have been incidences of convergences at different taxonomic levels. Here, we sequenced complete mitogenome of Aleurodicus rugioperculatus (Aleyrodidae: Aleurodicinae) to examine gene rearrangements and phylogenetic relationships within the family Aleyrodidae. We identified five gene blocks (I-V) in the whitefly ancestor that are shared plesiomorphies retained in different whitefly lineages. Gene block I is conserved in all whiteflies except three species (Tetraleurodes acaciae and two Bemisia species). Conversely, we detected 83 derived gene boundaries within the family. Mapping these gene boundaries onto a phylogenetic tree revealed that 16 were symplesiomorphies for two subfamilies; 9 were synapomorphies at different taxonomic levels, and 28 autapomorphies for individual species. Bayesian Inference (BI) and Maximum Likelihood (ML) phylogenetic analyses yielded similar topologies supporting the monophyly of Aleyrodinae and Aleurodicinae. Exclusion of PCG third codon positions from phylogenetic analyses improved both node support and consistency with prior analyses. To understand the significance of gene order convergence on phylogeny of the whiteflies, more species-level data is required.

Variation in synonymous nucleotide composition among genomes of sarbecoviruses and consequences for the origin of COVID-19

The subgenus Sarbecovirus includes two human viruses, SARS-CoV and SARS-CoV-2, respectively responsible for the SARS epidemic and COVID-19 pandemic, as well as many bat viruses and two pangolin viruses.Here, the synonymous nucleotide composition (SNC) of Sarbecovirus genomes was analysed by examining third codon-positions, dinucleotides, and degenerate codons. The results show evidence for the eigth following groups: (i) SARS-CoV related coronaviruses (SCoVrC including many bat viruses from China), (ii) SARS-CoV-2 related coronaviruses (SCoV2rC; including five bat viruses from Cambodia, Thailand and Yunnan), (iii) pangolin viruses, (iv) three bat viruses showing evidence of recombination between SCoVrC and SCoV2rC genomes, (v) two highly divergent bat viruses from Yunnan, (vi) the bat virus from Japan, (vii) the bat virus from Bulgaria, and (viii) the bat virus from Kenya. All these groups can be diagnosed by specific nucleotide compositional features except the one concerned by recombination between SCoVrC and SCoV2rC. In particular, SCoV2rC genomes are characterised by the lowest percentages of cyosine and highest percentages of uracil at third codon-positions, whereas the genomes of pangolin viruses exhibit the highest percentages of adenine at third codon-positions. I suggest that latitudinal and taxonomic differences in the imbalanced nucleotide pools available in host cells during viral replication can explain the seven groups of SNC here detected among Sarbecovirus genomes. A related effect due to hibernating bats is also considered. I conclude that the two independent host switches from Rhinolophus bats to pangolins resulted in convergent mutational constraints and that SARS-CoV-2 emerged directly from a horseshoe bat virus.

Comparative Mitogenomics and Phylogenetic Analyses of Pentatomoidea (Hemiptera: Heteroptera)

Pentatomoidea is the largest superfamily of Pentatomomorpha; however, the phylogenetic relationships among pentatomoid families have been debated for a long time. In the present study, we gathered the mitogenomes of 55 species from eight common families (Acanthosomatidae, Cydnidae, Dinidoridae, Scutelleridae, Tessaratomidae, Plataspidae, Urostylididae and Pentatomidae), including 20 newly sequenced mitogenomes, and conducted comparative mitogenomic studies with an emphasis on the structures of non-coding regions. Heterogeneity in the base composition, and contrasting evolutionary rates were encountered among the mitogenomes in Pentatomoidea, especially in Urostylididae, which may lead to unstable phylogenetic topologies. When the family Urostylididae is excluded in taxa sampling or the third codon positions of protein coding genes are removed, phylogenetic analyses under site-homogenous models could provide more stable tree topologies. However, the relationships between families remained the same in all PhyloBayes analyses under the site-heterogeneous mixture model CAT + GTR with different datasets and were recovered as (Cydnidae + (((Tessaratomidae + Dinidoridae) + (Plataspidae + Scutelleridae)) + ((Acanthosomatidae + Urostylididae) + Pentatomidae)))). Our study showed that data optimizing strategies after heterogeneity assessments based on denser sampling and the use of site-heterogeneous mixture models are essential for further analysis of the phylogenetic relationships of Pentatomoidea.

Molecular Evolution of Alphabaculovirus genomes: Evidence of Mutational bias and Natural selection

Codon usage is a reflection of evolutionary adaptation to environmental pressure. The pattern of usage may be unique to species of viruses, genomes of the same species or genes within the same genome. Here we have analysed the overall nucleotide composition and the nucleotides at different codon positions in the genomes of 6 Alphabaculoviruses. Principle Component Analysis (PCA) based on Relative Synonymous Codon Usage (RSCU) of all Open Reading Frames (ORFs) was employed to investigate the pattern of the codon usage. The results suggest the Alphabaculovirus genomes, except that of Agrotis Ipsilon mNPV (AgipNPV), are predominantly under an influence of a neutral mutation that bias toward A/T. The majority of the ORFs, except those of the AgipNPV, cluster at the same location in the 2-dimensional PCA map with one prominent outlier that has been identified as a P6.9 gene. The six Alpha-baculovirus P6.9 genes have a high G/C content, dissimilar to the majority of the ORFs. The G/C content is found to be significantly high at the 2 nd codon position, suggesting the influence of natural selection and perhaps reflecting its functional conservation in DNA packaging as well as its evolutionary relation to Protamine.

Nanopore whole genome sequencing and partitioned phylogenetic analysis supports a new salmonid alphavirus genotype (SAV7)

Salmon pancreas disease virus, more commonly known as salmonid alphavirus (SAV), is a single-stranded positive sense RNA virus and the causative agent of pancreas disease and sleeping disease in salmonids. In this study, a unique strain of SAV previously isolated from ballan wrasse was subjected to whole genome sequencing using nanopore sequencing. In order to accurately examine the evolutionary history of this strain in comparison to other SAV strains, a partitioned phylogenetic analysis was performed to account for variation in the rate of evolution for both individual genes and codon positions. Partitioning the genome alignments almost doubled the observed branch lengths in the phylogenetic tree when compared to the more common approach of applying one model of substitution across the genome and significantly increased the statistical fit of the best-fitting models of nucleotide substitution. Based on the genomic data, a valid case can be made for the viral strain examined in this study to be considered a new SAV genotype. In addition, this study adds to a growing number of studies in which SAV has been found to infect non-salmonid fish, and as such we have suggested that the viral species name be amended to the more inclusive ‘piscine alphavirus’.

HCV Genetic Diversity Can Be Used to Infer Infection Recency and Time since Infection

HIV-1 genetic diversity can be used to infer time since infection (TSI) and infection recency. We adapted this approach for HCV and identified genomic regions with informative diversity. We included 72 HCV/HIV-1 coinfected participants of the Swiss HIV Cohort Study, for whom reliable estimates of infection date and viral sequences were available. Average pairwise diversity (APD) was calculated over each codon position for the entire open reading frame of HCV. Utilizing cross validation, we evaluated the correlation of APD with TSI, and its ability to infer TSI via a linear model. We additionally studied the ability of diversity to classify infections as recent (infected for <1 year) or chronic, using receiver-operator-characteristic area under the curve (ROC-AUC) in 50 patients whose infection could be unambiguously classified as either recent or chronic. Measuring HCV diversity over third or all codon positions gave similar performances, and notable improvement over first or second codon positions. APD calculated over the entire genome enabled classification of infection recency (ROC-AUC = 0.76). Additionally, APD correlated with TSI (R2 = 0.33) and could predict TSI (mean absolute error = 1.67 years). Restricting the region over which APD was calculated to E2-NS2 further improved accuracy (ROC-AUC = 0.85, R2 = 0.54, mean absolute error = 1.38 years). Genetic diversity in HCV correlates with TSI and is a proxy for infection recency and TSI, even several years post-infection.

The role of selection in the evolution of marine turtles mitogenomes

Sea turtles are the only extant chelonian representatives that inhabit the marine environment. One key to successful colonization of this habitat is the adaptation to different energetic demands. Such energetic requirement is intrinsically related to the mitochondrial ability to generate energy through oxidative phosphorylation (OXPHOS) process. Here, we estimated Testudines phylogenetic relationships from 90 complete chelonian mitochondrial genomes and tested the adaptive evolution of 13 mitochondrial protein-coding genes of sea turtles to determine how natural selection shaped mitochondrial genes of the Chelonioidea clade. Complete mitogenomes showed strong support and resolution, differing at the position of the Chelonioidea clade in comparison to the turtle phylogeny based on nuclear genomic data. Codon models retrieved a relatively increased dN/dS (ω) on three OXPHOS genes for sea turtle lineages. Also, we found evidence of positive selection on at least three codon positions, encoded by NADH dehydrogenase genes (ND4 and ND5). The accelerated evolutionary rates found for sea turtles on COX2, ND1 and CYTB and the molecular footprints of positive selection found on ND4 and ND5 genes may be related to mitochondrial molecular adaptation to stress likely resulted from a more active lifestyle in sea turtles. Our study provides insight into the adaptive evolution of the mtDNA genome in sea turtles and its implications for the molecular mechanism of oxidative phosphorylation.

Demographic history and adaptive synonymous and nonsynonymous variants of nuclear genes in Rhododendron oldhamii (Ericaceae)

Demographic events are important in shaping the population genetic structure and exon variation can play roles in adaptive divergence. Twelve nuclear genes were used to investigate the species-level phylogeography of Rhododendron oldhamii, test the difference in the average GC content of coding sites and of third codon positions with that of surrounding non-coding regions, and test exon variants associated with environmental variables. Spatial expansion was suggested by R2 index of the aligned intron sequences of all genes of the regional samples and sum of squared deviations statistic of the aligned intron sequences of all genes individually and of all genes of the regional and pooled samples. The level of genetic differentiation was significantly different between regional samples. Significantly lower and higher average GC contents across 94 sequences of the 12 genes at third codon positions of coding sequences than that of surrounding non-coding regions were found. We found seven exon variants associated strongly with environmental variables. Our results demonstrated spatial expansion of R. oldhamii in the late Pleistocene and the optimal third codon position could end in A or T rather than G or C as frequent alleles and could have been important for adaptive divergence in R. oldhamii.

Interrogating Genomic-Scale Data to Resolve Recalcitrant Nodes in the Spider Tree of Life

Genome-scale data sets are converging on robust, stable phylogenetic hypotheses for many lineages; however, some nodes have shown disagreement across classes of data. We use spiders (Araneae) as a system to identify the causes of incongruence in phylogenetic signal between three classes of data: exons (as in phylotranscriptomics), noncoding regions (included in ultraconserved elements [UCE] analyses), and a combination of both (as in UCE analyses). Gene orthologs, coded as amino acids and nucleotides (with and without third codon positions), were generated by querying published transcriptomes for UCEs, recovering 1,931 UCE loci (codingUCEs). We expected that congeners represented in the codingUCE and UCEs data would form clades in the presence of phylogenetic signal. Noncoding regions derived from UCE sequences were recovered to test the stability of relationships. Phylogenetic relationships resulting from all analyses were largely congruent. All nucleotide data sets from transcriptomes, UCEs, or a combination of both recovered similar topologies in contrast with results from transcriptomes analyzed as amino acids. Most relationships inferred from low-occupancy data sets, containing several hundreds of loci, were congruent across Araneae, as opposed to high occupancy data matrices with fewer loci, which showed more variation. Furthermore, we found that low-occupancy data sets analyzed as nucleotides (as is typical of UCE data sets) can result in more congruent relationships than high occupancy data sets analyzed as amino acids (as in phylotranscriptomics). Thus, omitting data, through amino acid translation or via retention of only high occupancy loci, may have a deleterious effect in phylogenetic reconstruction.