Journal of Molecular Evolution最新文献

The genetic architecture of multiple sclerosis is complicated. Additionally, the disease incidence varies per population or per geographical region. A recent study gives convincing explanations about the north-south incidence gradient of multiple sclerosis in Europe, by analyzing ancient and modern human genomes. Interestingly, the evidence shows that multiple sclerosis associated immunogenetic variants underwent positive selection in Asian and European populations. Lifestyle and pathogen infections probably shaped the overall multiple sclerosis risk. These results complete the findings of previous studies that showed that a high percentage of the autoimmunity associated genetic variants are under selection pressure.

Adenosine-to-inosine (A-to-I) RNA editing recodes the genetic information. Apart from diversifying the proteome, another tempting advantage of RNA recoding is to correct deleterious DNA mutation and restore ancestral allele. Solid evidences for beneficial restorative editing are very rare in animals. By searching for "convergent recoding" under a phylogenetic context, we proposed this term for judging the potential restorative functions of particular editing site. For the well-known mammalian Gln>Arg (Q>R) recoding site, its ancestral state in vertebrate genomes was the pre-editing Gln, and all 470 available mammalian genomes strictly avoid other three equivalent ways to achieve Arg in protein. The absence of convergent recoding from His>Arg, or synonymous mutations on Gln codons, could be attributed to the strong maintenance on editing motif and structure, but the absence of direct A-to-G mutation is extremely unexpected. With similar ideas, we found cases of convergent recoding in Drosophila genus, reducing the possibility of their restorative function. In summary, we defined an interesting scenario of convergent recoding, the occurrence of which could be used as preliminary judgements for whether a recoding site has a sole restorative role. Our work provides novel insights to the natural selection and evolution of RNA editing.

Kinesins are eukaryotic microtubule motor proteins subdivided into conserved families with distinct functional roles. While many kinesin families are widespread in eukaryotes, each organismal lineage maintains a unique kinesin repertoire composed of many families with distinct numbers of genes. Previous genomic surveys indicated that land plant kinesin repertoires differ markedly from other eukaryotes. To determine when repertoires diverged during plant evolution, we performed robust phylogenomic analyses of kinesins in 24 representative plants, two algae, two animals, and one yeast. These analyses show that kinesin repertoires expand and contract coincident with major shifts in the biology of algae and land plants. One kinesin family and five subfamilies, each defined by unique domain architectures, emerged in the green algae. Four of those kinesin groups expanded in ancestors of modern land plants, while six other kinesin groups were lost in the ancestors of pollen-bearing plants. Expansions of different kinesin families and subfamilies occurred in moss and angiosperm lineages. Other kinesin families remained stable and did not expand throughout plant evolution. Collectively these data support a radiation of kinesin domain architectures in algae followed by differential positive and negative selection on kinesins families and subfamilies in different lineages of land plants.

Genome size variation in eukaryotes has myriad effects on organismal biology from the genomic to whole-organism level. Large genome size may be associated with lower selection efficiency because lower effective population sizes allow fixation of deleterious mutations via genetic drift, increasing genome size and decreasing selection efficiency. Because of a hypothesized negative relationship between genome size and recombination rate per base pair, increased genome size could also increase the effect of linked selection in the genome, decreasing the efficiency with which natural selection can fix or remove mutations. We used a transcriptomic dataset of 15 and a subset of six Neotropical salamander species ranging in genome size from 12 to 87 pg to study the relationship between genome size and efficiency of selection. We estimated dN/dS of salamanders with small and large genomes and tested for relaxation of selection in the larger genomes. Contrary to our expectations, we did not find a significant relationship between genome size and selection efficiency or strong evidence for higher dN/dS values in species with larger genomes for either species set. We also found little evidence for relaxation of selection in species with larger genomes. A positive correlation between genome size and range size (a proxy of population size) in this group disagrees with predictions of stronger drift in species with larger genomes. Our results highlight the complex interactions between the many forces shaping genomic variation in organisms with genomic gigantism.

Gene duplication followed by nucleotide differentiation is one of the simplest mechanisms to develop new functions for genes. However, the evolutionary processes underlying the divergence of multigene families remain controversial. We used multigene families found within the diversity of toxic proteins in centipede venom to test two hypotheses related to venom evolution: the two-speed mode of venom evolution and the rapid accumulation of variation in exposed residues (RAVER) model. The two-speed mode of venom evolution proposes that different types of selection impact ancient and younger venomous lineages with negative selection being the predominant form in ancient lineages and positive selection being the dominant form in younger lineages. The RAVER hypothesis proposes that, instead of different types of selection acting on different ages of venomous lineages, the different types of selection will selectively contribute to amino acid variation based on whether the residue is exposed to the solvent where it can potentially interact directly with toxin targets. This hypothesis parallels the longstanding understanding of protein evolution that suggests that residues found within the structural or active regions of the protein will be under negative or purifying selection, and residues that do not form part of these areas will be more prone to positive selection. To test these two hypotheses, we compared the venom of 26 centipedes from the order Scolopendromorpha from six currently recognized species from across North America using both transcriptomics and proteomics. We first estimated their phylogenetic relationships and uncovered paraphyly among the genus Scolopendra and evidence for cryptic diversity among currently recognized species. Using our phylogeny, we then characterized the diverse venom components from across the identified clades using a combination of transcriptomics and proteomics. We conducted selection-based analyses in the context of predicted three-dimensional properties of the venom proteins and found support for both hypotheses. Consistent with the two-speed hypothesis, we found a prevalence of negative selection across all proteins. Consistent with the RAVER hypothesis, we found evidence of positive selection on solvent-exposed residues, with structural and less-exposed residues showing stronger signal for negative selection. Through the use of phylogenetics, transcriptomics, proteomics, and selection-based analyses, we were able to describe the evolution of venom from an ancient venomous lineage and support principles of protein evolution that directly relate to multigene family evolution.

Quantum mechanical calculations are used to explore the thermodynamics of possible prebiotic synthesis of the building blocks of nucleic acids. Different combinations of D-ribofuranose (Ribf) and N-(2-aminoethyl)-glycine (AEG) (trifunctional connectors (TCs)); the nature of the Ribf, its anomeric form, and its ring puckering (conformation); and the nature of the nucleobases (recognition units (RUs)) are considered. The combinatorial explosion of possible nucleosides has been drastically reduced on physicochemical grounds followed by a detailed thermodynamic evaluation of alternative synthetic pathways. The synthesis of nucleosides containing N-(2-aminoethyl)-glycine (AEG) is predicted to be thermodynamically favored suggesting a possible role of AEG as a component of an ancestral proto-RNA that may have preceded today's nucleic acids. A new pathway for the building of free nucleotides (exemplified by 5'-uridine monophosphate (UMP)) and of AEG dipeptides is proposed. This new pathway leads to a spontaneous formation of free UMP assisted by an AEG nucleoside in an aqueous environment. This appears to be a workaround to the "water problem" that prohibits the synthesis of nucleotides in water.

By looking for a lack of homologs in a reference database of 27 well-annotated proteomes of primates and 52 well-annotated proteomes of other mammals, 170 putative human-specific proteins were identified. While most of them are deemed uncertain, 2 are known at the protein level and 23 at the transcript level, according to UniProt. Interestingly, 23 of these 25 proteins are found to be encoded or to have close homologs in an open reading frame of a long noncoding human RNA. However, half of them are predicted to be at least 80% globular, with a single structural domain, according to IUPred, and with at least 80% of ordered residues, according to flDPnn. Strikingly, there is a near-complete lack of structural knowledge about these proteins, with no tertiary structure presently available in the Protein Data Bank and a fair prediction for one of them in the AlphaFold Protein Structure Database. Moreover, knowledge about the function of these possibly key proteins remains scarce.

In the present work, we carried out a comparative genomic analysis to trace the evolutionary trajectory of the bacterial species that make up the Liquorilactobacillus genus, from the identification of genes and speciation/adaptation mechanisms in their unique characteristics to the identification of the pattern grouping these species. We present phylogenetic relationships between Liquorilactobacillus and related taxa such as Bacillus, basal lactobacilli and Ligilactobacillus, highlighting evolutionary divergences and lifestyle transitions across different taxa. The species of this genus share a core genome of 1023 genes, distributed in all COGs, which made it possible to characterize it as Liquorilactobacillus sensu lato: few amino acid auxotrophy, low genes number for resistance to antibiotics and general and specific cellular reprogramming mechanisms for environmental responses. These species were divided into four clades, with diversity being enhanced mainly by the diversity of genes involved in sugar metabolism. Clade 1 presented lower (< 70%) average amino acid identity with the other clades, with exclusive or absent genes, and greater distance in the genome compared to clades 2, 3 and 4. The data pointed to an ancestor of clades 2, 3 and 4 as being the origin of the genus Ligilactobacillus, while the species of clade 1 being closer to the ancestral Bacillus. All these traits indicated that the species of clade 1 could be soon separated in a distinct genus.

Cetaceans and pinnipeds are lineages of mammals that have independently returned to the aquatic environment, acquiring varying degrees of dependence on it while sharing adaptations for underwater living. Here, we focused on one critical adaptation from both groups, their ability to withstand the ischemia and reperfusion experienced during apnea diving, which can lead to the production of reactive oxygen species (ROS) and subsequent oxidative damage. Previous studies have shown that cetaceans and pinnipeds possess efficient antioxidant enzymes that protect against ROS. In this study, we investigated the molecular evolution of key antioxidant enzyme genes (CAT, GPX3, GSR, PRDX1, PRDX3, and SOD1) and the ROS-producing gene XDH, in cetaceans and pinnipeds lineages. We used the ratio of non-synonymous (dN) to synonymous (dS) substitutions as a measure to identify signatures of adaptive molecular evolution in these genes within and between the two lineages. Additionally, we performed protein modeling and variant impact analyzes to assess the functional consequences of observed mutations. Our findings revealed distinct selective regimes between aquatic and terrestrial mammals in five of the examined genes, including divergences within cetacean and pinniped lineages, between ancestral and recent lineages and between crowns groups. We identified specific sites under positive selection unique to Cetacea and Pinnipedia, with one site showing evidence of convergent evolution in species known for their long and deep-diving capacities. Notably, many sites under adaptive selection exhibited radical changes in amino acid properties, with some being damaging mutations in human variations, but with no apparent detrimental impacts on aquatic mammals. In conclusion, our study provides insights into the adaptive changes that have occurred in the antioxidant systems of aquatic mammals throughout their evolutionary history. We observed both distinctive features within each group of Cetacea and Pinnipedia and instances of convergence. These findings highlight the dynamic nature of the antioxidant system in response to challenges of the aquatic environment and provide a foundation for further investigations into the molecular mechanisms underlying these adaptations.

Brucellosis is a notifiable disease induced by a facultative intracellular Brucella pathogen. In this study, eight Brucella abortus and eighteen Brucella melitensis strains from Egypt were annotated and compared with RB51 and REV1 vaccines respectively. RAST toolkit in the BV-BRC server was used for annotation, revealing genome length of 3,250,377 bp and 3,285,803 bp, 3289 and 3323 CDS, 48 and 49 tRNA genes, the same number of rRNA (3) genes, 583 and 586 hypothetical proteins, 2697 and 2726 functional proteins for B. abortus and B. melitensis respectively. B. abortus strains exhibit a similar number of candidate genes, while B. melitensis strains showed some differences, especially in the SRR19520422 Faiyum strain. Also, B. melitensis clarified differences in antimicrobial resistance genes (KatG, FabL, MtrA, MtrB, OxyR, and VanO-type) in SRR19520319 Faiyum and (Erm C and Tet K) in SRR19520422 Faiyum strain. Additionally, the whole genome phylogeny analysis proved that all B. abortus strains were related to vaccinated animals and all B. melitensis strains of Menoufia clustered together and closely related to Gharbia, Dameitta, and Kafr Elshiek. The Bowtie2 tool identified 338 (eight B. abortus) and 4271 (eighteen B. melitensis) single nucleotide polymorphisms (SNPs) along the genomes. These variants had been annotated according to type and impact. Moreover, thirty candidate genes were predicted and submitted at GenBank (24 in B. abortus) and (6 in B. melitensis). This study contributes significant insights into genetic variation, virulence factors, and vaccine-related associations of Brucella pathogens, enhancing our knowledge of brucellosis epidemiology and evolution in Egypt.