Synonymous Substitution Rates Research Articles

Pervasive relaxed selection in termite genomes.

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

Death and resurrection of an inhibitory Na+ binding site in a human epithelial Na+ channel subunit

Epithelial Na+ channels (ENaC) transport and sense Na+ in a variety of tissues. These channels selectively transport Na+ and are also inhibited by it, an effect called Na+ self-inhibition. ENaCs are heterotrimers comprising α, β, and γ subunits. A fourth ENaC subunit δ, related by gene duplication to α, can replace the α subunit and confer unique functional properties. Among these are Na+ self-inhibition, present in αβγ, but absent in δβγ. We previously identified a putative Na+ binding site centered around mouse αD365. Functional data provided evidence that a homologous site was present in the Xenopus α and δ subunits, and structural data provided evidence for a cation bound at the homologous site in the human α subunit (D338). We hypothesized that 1) the α/δ ancestral subunit had an inhibitory Na+ binding site, 2) that the site was lost on the δ subunit lineage, and 3) that the human δ subunit site could be resurrected. We determined a phylogenetic tree from >800 α and δ subunit sequences, and examined substitution at the position of the key aspartate. Among the α subunits, the key aspartate was omnipresent, absent only in the sloth. Among the δ subunits, the key aspartate was absent in >60 species, including primates and other mammalian orders, and some perching birds. By contrasting nested models of site evolution, we found that site retention was 67-fold more favorable for α than for δ (p=0.02). When we measured non-synonymous to synonymous substitution rate ratios (dN/dS), we found that the α subunit was under purifying selection (p < 0.001) while the δ subunit was not (p = 0.4), providing evidence for a loss of selection pressure after the divergence of the α and δ subunits. Results for other residues at the putative inhibitory Na+ binding site were similar. The treeshrew ( Tupaia belangerie) was the closest primate relative with a putative Na+ self-inhibition site in its δ subunit. We confirmed the site’s functionality by complementary expression with human β and γ subunits in Xenopus oocytes and performing two-electrode voltage clamp experiments. Mutating the homologous aspartate in treeshrew δ significantly reduced Na+ self-inhibition (p < 0.0001). To confer Na+ self-inhibition to human δ, we mutated residues at the putative site to their treeshrew δ equivalents. Restoring the key aspartate in human δ (dP314D) resulted in discernable Na+ self-inhibition (p=0.003), and treeshrew-ifying 2 additional sites (G136M and L322E) conferred Na+ self-inhibition in human δ like that measured in human α. We then examined ENaC subunit expression in tissues from several species. In contrast to Xenopus laevis, we did not detect δ subunit transcripts in kidneys from the treeshrew, chicken, or opossum. Our study provides evidence that the site centered around human aD338 is a bona fide allosteric Na+ binding site, and that selection pressure maintained this site in the α subunit but not the δ subunit, possibly due to differential renal expression. DK125439. This is the full abstract presented at the American Physiology Summit 2024 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.

Dynamic changes in the plastid and mitochondrial genomes of the angiosperm Corydalis pauciovulata (Papaveraceae)

BackgroundCorydalis DC., the largest genus in the family Papaveraceae, comprises > 465 species. Complete plastid genomes (plastomes) of Corydalis show evolutionary changes, including syntenic arrangements, gene losses and duplications, and IR boundary shifts. However, little is known about the evolution of the mitochondrial genome (mitogenome) in Corydalis. Both the organelle genomes and transcriptomes are needed to better understand the relationships between the patterns of evolution in mitochondrial and plastid genomes.ResultsWe obtained complete plastid and mitochondrial genomes from Corydalis pauciovulata using a hybrid assembly of Illumina and Oxford Nanopore Technologies reads to assess the evolutionary parallels between the organelle genomes. The mitogenome and plastome of C. pauciovulata had sizes of 675,483 bp and 185,814 bp, respectively. Three ancestral gene clusters were missing from the mitogenome, and expanded IR (46,060 bp) and miniaturized SSC (202 bp) regions were identified in the plastome. The mitogenome and plastome of C. pauciovulata contained 41 and 67 protein-coding genes, respectively; the loss of genes was a plastid-specific event. We also generated a draft genome and transcriptome for C. pauciovulata. A combination of genomic and transcriptomic data supported the functional replacement of acetyl-CoA carboxylase subunit β (accD) by intracellular transfer to the nucleus in C. pauciovulata. In contrast, our analyses suggested a concurrent loss of the NADH-plastoquinone oxidoreductase (ndh) complex in both the nuclear and plastid genomes. Finally, we performed genomic and transcriptomic analyses to characterize DNA replication, recombination, and repair (DNA-RRR) genes in C. pauciovulata as well as the transcriptomes of Liriodendron tulipifera and Nelumbo nuicifera. We obtained 25 DNA-RRR genes and identified their structure in C. pauciovulata. Pairwise comparisons of nonsynonymous (dN) and synonymous (dS) substitution rates revealed that several DNA-RRR genes in C. pauciovulata have higher dN and dS values than those in N. nuicifera.ConclusionsThe C. pauciovulata genomic data generated here provide a valuable resource for understanding the evolution of Corydalis organelle genomes. The first mitogenome of Papaveraceae provides an example that can be explored by other researchers sequencing the mitogenomes of related plants. Our results also provide fundamental information about DNA-RRR genes in Corydalis and their related rate variation, which elucidates the relationships between DNA-RRR genes and organelle genome stability.

Auto-phylo v2 and auto-phylo-pipeliner: building advanced, flexible, and reusable pipelines for phylogenetic inferences, estimation of variability levels and identification of positively selected amino acid sites.

The vast amount of genome sequence data that is available, and that is predicted to drastically increase in the near future, can only be efficiently dealt with by building automated pipelines. Indeed, the Earth Biogenome Project will produce high-quality reference genome sequences for all 1.8 million named living eukaryote species, providing unprecedented insight into the evolution of genes and gene families, and thus on biological issues. Here, new modules for gene annotation, further BLAST search algorithms, further multiple sequence alignment methods, the adding of reference sequences, further tree rooting methods, the estimation of rates of synonymous and nonsynonymous substitutions, and the identification of positively selected amino acid sites, have been added to auto-phylo (version 2), a recently developed software to address biological problems using phylogenetic inferences. Additionally, we present auto-phylo-pipeliner, a graphical user interface application that further facilitates the creation and running of auto-phylo pipelines. Inferences on S-RNase specificity, are critical for both cross-based breeding and for the establishment of pollination requirements. Therefore, as a test case, we develop an auto-phylo pipeline to identify amino acid sites under positive selection, that are, in principle, those determining S-RNase specificity, starting from both non-annotated Prunus genomes and sequences available in public databases.

Complete Mitochondrial Genome and Phylogenetic Position of Nurudea zhengii Ren (Insecta, Hemiptera, Aphididae).

Nurudea zhengii Ren was identified by aphid morphological characteristics as well as the gall shape and host plant species, and placed in the tribe Fordini (Hemiptera, Aphididae, Eriosomatinae). Here, its whole genome was firstly sequenced by a genome-skimming method and its mitochondrial genome (mitogenome) was assembled to examine its genetic variation and phylogenetic position. The complete mitogenome of Nurudea zhengii is 15,392bp in length, and consists of 13 protein-coding genes, 22 tRNAs, two rRNAs and one D-loop region. The gene order follows the mitogenomes of the other Rhus gall aphids, and similarly has an AT bias with the content of 83.9%. The majority strand is A-skewed and C-skewed, and shows opposite skewness for G-skewed in the minority strands. The ratios of nonsynonymous to synonymous substitution rates of protein-coding genes are lower than one except for ATP8, which indicated that ATP8 was undergoing positive selection. Phylogenetic analysis among the Rhus gall aphids based on 13 protein-coding genes and two rRNA genes showed that N. zhengii was sister to N. shiraii, and then clustered with N. yanoniella as a group with high support value. The two species, N. shiraii and N. yanoniella, share the same host plant Rhus chinensis, while the host of N. zhengii is R. hypoleuca. However, the phylogenetic relationship indicated that the taxa sharing the same host plant were not absolutely clustered as the closest taxa at least at species level.

Characterization of the complete mitochondrial genome of Paecilomyces variotii and comparative evolutionary mitochondriomics of 36 fungi.

Paecilomyces variotii has important economic value in stimulating crop growth, biodegradation, and other aspects. Up to now, there are no research reports on its mitochondrial genome. The mitochondrial genome of Paecilomyces variotii was determined with the next-generation sequencing method (Illumina, NovaSeq), and its characteristics were analyzed using various bioinformatics approaches. The length of complete mitochondrial genome sequence of P. variotii is 40,965bp and consists of 14 protein-coding genes, 2 ribosomal RNA genes, 1 ribosomal protein S3 gene, 26 transport RNA genes. The results of phylogenetics analysis using Bayesian inference and Maximum likelihood methods showed that P. variotii belongs to the Eurotiales order in the Thermoascaceae family, and 9 genera within the Eurotiomycetes class were effectively distinguished with high support rates (bootstrap value > 92% and posterior probabilities > 99%). The analysis of synonymous substitution rates and nonsynonymous substitution rates indicated that the Ka/Ks values of the 14 PCGs in the mitochondrial genomes of the two orders in the Eurotiomycetes class ranged from 0 to 0.4333. This study revealed the structural and sequence information characteristics of the mitochondrial genome of P. variotii, and the phylogenetic results strongly support its classification within the family Thermoascaceae, consistent with traditional morphological taxonomy studies. The 14 PCGs in the mitochondrial genomes of the two orders in the Eurotiomycetes class are subject to strong purifying (negative) selection. The results of this research provides an important molecular basis for the development of genomics, evolutionary genetics and molecular markers of P. variotii in the future.

Combining extracellular matrix proteome and phosphoproteome of chickpea and meta-analysis reveal novel proteoforms and evolutionary significance of clade-specific wall-associated events in plant.

Extracellular matrix (ECM) plays central roles in cell architecture, innate defense and cell wall integrity (CWI) signaling. During transition to multicellularity, modular domain structures of ECM proteins and proteoforms have evolved due to continuous adaptation across taxonomic clades under different ecological niche. Although this incredible diversity has to some extent been investigated at protein level, extracellular phosphorylation events and molecular evolution of ECM proteoform families remains unexplored. We developed matrisome proteoform atlas in a grain legume, chickpea and performed meta-analyses of 74 plant matrisomes. MS/MS analysis identified 1,424 proteins and 315 phosphoproteins involved in diverse functions. Cross-species ECM protein network identified proteoforms associated with CWI maintenance system. Phylogenetic characterization of eighteen matrix protein families highlighted the role of taxon-specific paralogs and orthologs. Novel information was acquired on gene expansion and loss, co-divergence, sub functionalization and neofunctionalization during evolution. Modular networks of matrix protein families and hub proteins showed higher diversity across taxonomic clades than among organs. Furthermore, protein families differ in nonsynonymous to synonymous substitution rates. Our study pointed towards the matrix proteoform functionality, sequence divergence variation, interactions between wall remodelers and molecular evolution using a phylogenetic framework. This is the first report on comprehensive matrisome proteoform network illustrating presence of CWI signaling proteins in land plants.

Genome copy number predicts extreme evolutionary rate variation in plant mitochondrial DNA

Nuclear and organellar genomes can evolve at vastly different rates despite occupying the same cell. In most bilaterian animals, mitochondrial DNA (mtDNA) evolves faster than nuclear DNA, whereas this trend is generally reversed in plants. However, in some exceptional angiosperm clades, mtDNA substitution rates have increased up to 5,000-fold compared with closely related lineages. The mechanisms responsible for this acceleration are generally unknown. Because plants rely on homologous recombination to repair mtDNA damage, we hypothesized that mtDNA copy numbers may predict evolutionary rates, as lower copy numbers may provide fewer templates for such repair mechanisms. In support of this hypothesis, we found that copy number explains 47% of the variation in synonymous substitution rates of mtDNA across 60 diverse seed plant species representing ~300 million years of evolution. Copy number was also negatively correlated with mitogenome size, which may be a cause or consequence of mutation rate variation. Both relationships were unique to mtDNA and not observed in plastid DNA. These results suggest that homologous recombinational repair plays a role in driving mtDNA substitution rates in plants and may explain variation in mtDNA evolution more broadly across eukaryotes. Our findings also contribute to broader questions about the relationships between mutation rates, genome size, selection efficiency, and the drift-barrier hypothesis.

Genome assembly provides insights into the genome evolution of Baccaurea ramiflora Lour.

Baccaurea ramiflora Lour., an evergreen tree of the Baccaurea genus of the Phyllanthaceae family, is primarily distributed in South Asia, Southeast Asia, and southern China, including southern Yunnan Province. It is a wild or semi-cultivated tree species with ornamental, edible, and medicinal value, exhibiting significant development potential. In this study, we present the whole-genome sequencing of B. ramiflora, employing a combination of PacBio SMRT and Illumina HiSeq 2500 sequencing techniques. The assembled genome size was 975.8 Mb, with a contig N50 of 509.33 kb and the longest contig measuring 7.74 Mb. The genome comprises approximately 73.47% highly repetitive sequences, of which 52.1% are long terminal repeat–retrotransposon sequences. A total of 29,172 protein-coding genes were predicted, of which 25,980 (89.06%) have been annotated, Additionally, 3452 non-coding RNAs were identified. Comparative genomic analysis revealed a close relationship between B. ramiflora and the Euphorbiaceae family, with both being sister groups that diverged approximately 59.9 million years ago. During the evolutionary process, B. ramiflora exhibited positive selection in 278 candidate genes. Synonymous substitution rate and collinearity analysis demonstrated that B. ramiflora underwent a single ancient genome-wide triploidization event, without recent genome-wide duplication events. This high-quality B. ramiflora genome provides a valuable resource for basic research and tree improvement programs focusing on the Phyllanthaceae family.

Characterization of Complete Mitochondrial Genome of Badri Breed of Bos indicus (Bovidae: Bovinae): Selection Pressure and Comparative Analysis.

High-altitude mammals are often subject to specific environmental obstacles, which exert selective pressure on their physiological and morphological traits, hence driving their evolutionary processes. It is anticipated that these circumstances will lead to the adaptive evolution of protein-coding genes (PCGs) in the mitochondrial genome, which play a crucial role in the oxidative phosphorylation system. In this study, we have generated the complete mitochondrial genome of the Badri breed of Bos indicus inhabiting a high-altitude environment to test the signatures of adaptive evolution on PCGs and their phylogenetic relationships. The complete mitogenome of the Badri breed is 16,339bp and most tRNAs showed typical clover-leaf secondary structure with a few exceptions, like trnS1 and trnS2 without DHU arm and trnK without DHU loop. Comparative analysis of PCGs indicated that cox1 is the most conserved, while atp6 is the most variable gene. Moreover, the ratios of non-synonymous to synonymous substitution rates indicated the purifying selection (Ka/Ks < 1) in the protein-coding genes that shape the diversity in mitogenome of Bos indicus. Furthermore, Branch-site model (BSM) suggested thatcox1, cox2, nad3, nad4L, and nad6 underwent stronger purifying selection (ω < 1) than other PCGs in 15 breeds of 4 species, including Badri. BSM also detected 10 positive sites in PCGs and one in 13 PCGs concatenated dataset. Additional analyses in Datamonkey indicated 11 positive sites and 23 purifying sites in the concatenated dataset, a relaxation of selection strength innad3, and no evidence of episodic diversifying selection in any PCGs. Phylogeny revealed the sister relationship of the Badri with other breeds of Bos indicus as well as Bos frontalis (Gayal-2). The mitogenome of the Badri breed is an important genomic resource for conservation genetics of this species and also contributes to the understanding of the adaptive evolution of mitochondrial protein coding genes.

Genome-wide analysis of the peanut CaM/CML gene family reveals that the AhCML69 gene is associated with resistance to Ralstonia solanacearum

BackgroundCalmodulins (CaMs)/CaM-like proteins (CMLs) are crucial Ca2+-binding sensors that can decode and transduce Ca2+ signals during plant development and in response to various stimuli. The CaM/CML gene family has been characterized in many plant species, but this family has not yet been characterized and analyzed in peanut, especially for its functions in response to Ralstonia solanacearum. In this study, we performed a genome-wide analysis to analyze the CaM/CML genes and their functions in resistance to R. solanacearum.ResultsHere, 67, 72, and 214 CaM/CML genes were identified from Arachis duranensis, Arachis ipaensis, and Arachis hypogaea, respectively. The genes were divided into nine subgroups (Groups I-IX) with relatively conserved exon‒intron structures and motif compositions. Gene duplication, which included whole-genome duplication, tandem repeats, scattered repeats, and unconnected repeats, produced approximately 81 pairs of homologous genes in the AhCaM/CML gene family. Allopolyploidization was the main reason for the greater number of AhCaM/CML members. The nonsynonymous (Ka) versus synonymous (Ks) substitution rates (less than 1.0) suggested that all homologous pairs underwent intensive purifying selection pressure during evolution. AhCML69 was constitutively expressed in different tissues of peanut plants and was involved in the response to R. solanacearum infection. The AhCML69 protein was localized in the cytoplasm and nucleus. Transient overexpression of AhCML69 in tobacco leaves increased resistance to R. solanacearum infection and induced the expression of defense-related genes, suggesting that AhCML69 is a positive regulator of disease resistance.ConclusionsThis study provides the first comprehensive analysis of the AhCaM/CML gene family and potential genetic resources for the molecular design and breeding of peanut bacterial wilt resistance.

Comparative Mitogenome Analyses of Fifteen Ramshorn Snails and Insights into the Phylogeny of Planorbidae (Gastropoda: Hygrophila).

Ramshorn snails from the family Planorbidae are important freshwater snails due to their low trophic level, and some of them act as intermediate hosts for zoonotic trematodes. There are about 250 species from 40 genera of Planorbidae, but only 14 species from 5 genera (Anisus, Biomphalaria, Bulinus, Gyraulus, and Planorbella) have sequenced complete mitochondrial genomes (mitogenomes). In this study, we sequenced and assembled a high-quality mitogenome of a ramshorn snail, Polypylis sp. TS-2018, which represented the first mitogenome of the genus. The mitogenome of Polypylis sp. TS-2018 is 13,749 bp in length, which is shorter than that of most gastropods. It contains 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and 2 ribosomal RNA (rRNA). We compared mitogenome characteristics, selection pressure, and gene rearrangement among all of the available mitogenomes of ramshorn snails. We found that the nonsynonymous and synonymous substitution rates (Ka/Ks) of most PCGs indicated purifying and negative selection, except for atp8 of Anisus, Biomphalaria, and Gyraulus, which indicated positive selection. We observed that transpositions and reverse transpositions occurred on 10 tRNAs and rrnS, which resulted in six gene arrangement types. We reconstructed the phylogenetic trees using the sequences of PCGs and rRNAs and strongly supported the monophyly of each genus, as well as three tribes in Planorbidae. Both the gene rearrangement and phylogenetic results suggested that Polypylis had a close relationship with Anisus and Gyraulus, while Bulinus was the sister group to all of the other genera. Our results provide useful data for further investigation of species identification, population genetics, and phylogenetics among ramshorn snails.

Remarkable mitochondrial genome heterogeneity in Meniocus linifolius (Brassicaceae).

Detailed analyses of 16 genomes identified a remarkable acceleration of mutation rate, hence mitochondrial sequence and structural heterogeneity, in Meniocus linifolius (Brassicaceae). The powerhouse, mitochondria, in plantsfeature high levels of structural variation, while the encoded genes are normally conserved. However, the substitution rates and spectra of mitochondria DNA within the Brassicaceae, a family with substantial scientific and economic importance, have not been adequately deciphered. Here, by analyzing three newly assembled and 13 known mitochondrial genomes (mitogenomes), we report the highly variable genome structure and mutation rates in Brassicaceae. The genome sizes and GC contents are 196,604bp and 46.83%, 288,122bp and 44.79%, and 287,054bp and 44.93%, for Meniocus linifolius (Mli), Crucihimalaya lasiocarpa (Cla), and Lepidium sativum (Lsa), respectively. In total, 29, 33, and 34 protein-coding genes (PCGs) and 14, 18, and 18 tRNAs are annotated for Mli, Cla, and Lsa, respectively, while all mitogenomes contain one complete circular molecule with three rRNAs and abundant RNA editingsites.The Mli mitogenome features four conformations likely mediatedby the two pairs of long repeats, while at the same time seems to have an unusual evolutionary history due to higher GC content, loss of more genes and sequences, but having more repeats and plastid DNA insertions. Corroborating with these, an ambiguous phylogenetic position with long branch length and elevated synonymous substitution rate in nearly all PCGs are observed for Mli. Taken together, our results reveal a high level of mitogenome heterogeneity at the family level and provide valuable resources for further understanding the evolutionary pattern of organelle genomes in Brassicaceae.

Characterization of Four Complete Mitogenomes of Monolepta Species and Their Related Phylogenetic Implications.

Monolepta is one of the diverse genera in the subfamily Galerucinae, including 708 species and 6 sub-species worldwide. To explore the information on the mitogenome characteristics and phylogeny of the section "Monoleptites", especially the genus Monolepta, we obtained the newly completed mitochondrial genomes (mitogenomes) of four Monolepta species using high-throughput sequencing technology. The lengths of these four new mitochondrial genomes are 16,672 bp, 16,965 bp, 16,012 bp, and 15,866 bp in size, respectively. All four mitochondrial genomes include 22 transfer RNA genes (tRNAs), 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), and one control region, which is consistent with other Coleoptera. The results of the nonsynonymous with synonymous substitution rates showed that ND6 had the highest evolution rate, while COI displayed the lowest evolution rate. The substitution saturation of three datasets (13 PCGs_codon1, 13 PCGs_codon2, 13 PCGs_codon3) showed that there was no saturation across all datasets. Phylogenetic analyses based on three datasets (ND1, 15 genes of mitogenomes, and 13 PCGs_AA) were carried out using maximum likelihood (ML) and Bayesian inference (BI) methods. The results showed that mitogenomes had a greater capacity to resolve the main clades than the ND1 gene at the suprageneric and species levels. The section "Monoleptites" was proven to be a monophyletic group, while Monolepta was a non-monophyletic group. Based on ND1 data, the newly sequenced species whose antennal segment 2 was shorter than 3 were split into several clades, while, based on the mitogenomic dataset, the four newly sequenced species had close relationships with Paleosepharia. The species whose antennal segment 2 was as long as 3 were split into two clades, which indicated that the characteristic of "antennal segment 2 as long as 3" of the true "Monolepta" evolved multiple times in several subgroups. Therefore, to explore the relationships among the true Monolepta, the most important thing is to perform a thorough revision of Monolepta and related genera in the future.

Large-scale Pan Genomic Analysis of Mycobacterium tuberculosis Reveals Key Insights Into Molecular Evolutionary Rate of Specific Processes and Functions.

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis (TB), an infectious disease that is a major killer worldwide. Due to selection pressure caused by the use of antibacterial drugs, Mtb is characterised by mutational events that have given rise to multi drug resistant (MDR) and extensively drug resistant (XDR) phenotypes. The rate at which mutations occur is an important factor in the study of molecular evolution, and it helps understand gene evolution. Within the same species, different protein-coding genes evolve at different rates. To estimate the rates of molecular evolution of protein-coding genes, a commonly used parameter is the ratio dN/dS, where dN is the rate of non-synonymous substitutions and dS is the rate of synonymous substitutions. Here, we determined the estimated rates of molecular evolution of select biological processes and molecular functions across 264 strains of Mtb. We also investigated the molecular evolutionary rates of core genes of Mtb by computing the dN/dS values, and estimated the pan genome of the 264 strains of Mtb. Our results show that the cellular amino acid metabolic process and the kinase activity function evolve at a significantly higher rate, while the carbohydrate metabolic process evolves at a significantly lower rate for M. tuberculosis. These high rates of evolution correlate well with Mtb physiology and pathogenicity. We further propose that the core genome of M. tuberculosis likely experiences varying rates of molecular evolution which may drive an interplay between core genome and accessory genome during M. tuberculosis evolution.

The most exposed regions of SARS-CoV-2 structural proteins are subject to strong positive selection and gene overlap may locally modify this behavior.

The SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) pandemic that emerged in 2019 has been an unprecedented event in international science, as it has been possible to sequence millions of genomes, tracking their evolution very closely. This has enabled various types of secondary analyses of these genomes, including the measurement of their sequence selection pressure. In this work, we have been able to measure the selective pressure of all the described SARS-CoV-2 genes, even analyzed by sequence regions, and we show how this type of analysis allows us to separate the genes between those subject to positive selection (usually those that code for surface proteins or those exposed to the host immune system) and those subject to negative selection because they require greater conservation of their structure and function. We have also seen that when another gene with an overlapping reading frame appears within a gene sequence, the overlapping sequence between the two genes evolves under a stronger purifying selection than the average of the non-overlapping regions of the main gene. We propose this type of analysis as a useful tool for locating and analyzing all the genes of a viral genome when an adequate number of sequences are available.IMPORTANCEWe have analyzed the selection pressure of all severe acute respiratory syndrome coronavirus 2 genes by means of the nonsynonymous (Ka) to synonymous (Ks) substitution rate. We found that protein-coding genes are exposed to strong positive selection, especially in the regions of interaction with other molecules (host receptor and genome of the virus itself). However, overlapping coding regions are more protected and show negative selection. This suggests that this measure could be used to study viral gene function as well as overlapping genes.

Concern regarding H3-subtype avian influenza virus.

The H3-subtype of avian influenza virus (AIV) is one of the most frequently detected low pathogenic avian influenza virus (LPAIV) subtypes in birds and fowls, causing substantial economic loss to the poultry industry. Most importantly, besides poultry, mammals could also be infected with it, such as swines, canines, equines, felines, and humans, posing a serious public health threat. This allows the virus to persist widely in poultry and wild birds for a long time, where it may mix with other subtypes, providing conditions for viral recombination or reassortment. Currently, the monitoring of H3-subtype AIV is inadequate, and there is a lack of effective prevention and control measures for H3-subtype AIV. Here, the epidemiology, phylogeny, and genetic variation of H3-subtype AIV were analyzed, and nonsynonymous and synonymous substitution rates (dN/dS) were calculated. Through these steps, we aimed to clarify the current epidemiological feature and evolutionary characteristics of H3-subtype AIV, and provide an operative reference for future scientific control of H3-subtype AIV.

The pace of mitochondrial molecular evolution varies with seasonal migration distance.

Animals that engage in long-distance seasonal migration experience strong selective pressures on their metabolic performance and life history, with potential consequences for molecular evolution. Species with slow life histories typically show lower rates of synonymous substitution (dS) than "fast" species. Previous research suggests long-distance seasonal migrants have a slower life history strategy than short-distance migrants, raising the possibility that rates of molecular evolution may covary with migration distance. Additionally, long-distance migrants may face strong selection on metabolically-important mitochondrial genes due to their long-distance flights. Using over 1,000 mitochondrial genomes, we assessed the relationship between migration distance and mitochondrial molecular evolution in 39 boreal-breeding migratory bird species. We show that migration distance correlates negatively with dS, suggesting that the slow life history associated with long-distance migration is reflected in rates of molecular evolution. Mitochondrial genes in every study species exhibited evidence of purifying selection, but the strength of selection was greater in short-distance migrants, contrary to our predictions. This result may indicate effects of selection for cold tolerance on mitochondrial evolution among species overwintering at high latitudes. Our study demonstrates that the pervasive correlation between life history and molecular evolutionary rates exists in the context of differential adaptations to seasonality.

Genomic estimates of mutation and substitution rates contradict the evolutionary speed hypothesis of the latitudinal diversity gradient.

The latitudinal diversity gradient (LDG) refers to a decrease in biodiversity from the equator to the poles. The evolutionary speed hypothesis, backed by the metabolic theory of ecology, asserts that nucleotide mutation and substitution rates per site per year are higher and thereby speciation rates are higher at higher temperatures, generating the LDG. However, prior empirical investigations of the relationship between the temperature and mutation or substitution rate were based on a few genes and the results were mixed. We here revisit this relationship using genomic data. No significant correlation between the temperature and mutation rate is found in 13 prokaryotes or in 107 eukaryotes. An analysis of 234 diverse trios of bacterial taxa indicates that the synonymous substitution rate is not significantly associated with the growth temperature. The same data, however, reveal a significant negative association between the nonsynonymous substitution rate and temperature, which is explainable by a larger fraction of detrimental nonsynonymous mutations at higher temperatures due to a stronger demand for protein stability. We conclude that the evolutionary speed hypothesis of the LDG is unsupported by genomic data and advise that future mechanistic studies of the LDG should focus on other hypotheses.

Characterization of complete mitochondrial genome of three Horse flies of the genus Tabanus (Diptera: Tabanidae): comparative analysis.

Tabanidae (Horse-Flies or Deer-Flies) are one of the most economically important as well as medically significant haematophagous insect family within the order Diptera. Members of this group are also responsible for the mortality of substantial number of live-stock every year. Due to their pathogen transmission potential and vector competencies makes them an important insect group to study. Till now, mitochondrial genome of 18 species of tabanids were available. The complete mitogenome of three species T. diversifrons (15,809bp), T. rubidus (15,878bp) and T. tenens (15,872bp) were generated by Next generation sequencing method. They consist 37 genes, with a positive AT skew and a negative GC skew. The gene order of these three species is similar to the typical gene arrangement of infra-order Tabanomorpha. Most of the tRNAs showed typical clover-leaf secondary structure except trnS1, which lacks the DHU arm. The sliding window analysis showed that thenad4L is the most conserved while atp8, and nad6 are the most variable genes. Moreover, the ratios of non-synonymous to synonymous substitution rates indicated that all PCGs under the purifying selection. Phylogeny revealed Chrysops and Haematopota are monophyletic while species of Hybomitra are nested within the polyphyletic clade of Tabanus. T. diversifrons exhibits sister relationship with Atylotus miser. Two morphologically divergent species T. rubidus and T. tenens are found to be genetically similar and indistinguishable by mitochondrial genome. The hypervariable genes like atp8 and nad6 can be used as molecular markers for the identification of recently diverged lineages of family Tabanidae. Further, to address uncertainties arising from the two morphological divergent species, it is imperative to obtain data from nuclear gene markers.