Beginner's Guide on the Use of PAML to Detect Positive Selection

The fast-evolving FIKK kinase family of Plasmodium falciparum can be inhibited by a single compound

Of 250 Plasmodium species, 6 infect humans, with P. falciparum causing over 95% of 600,000 annual malaria-related deaths. Its pathology arises from host cell remodelling driven by over 400 exported parasite proteins, including the FIKK kinase family. About one million years ago, a bird-infecting Plasmodium species crossed into great apes and a single non-exported FIKK kinase gained an export element. This led to a rapid expansion into 15-21 atypical, exported Ser/Thr effector kinases. Here, using genomic and proteomic analyses, we demonstrate FIKK differentiation via changes in subcellular localization, expression timing and substrate motifs, which supports an individual important role in host-pathogen interactions. Structural data and AlphaFold2 predictions reveal fast-evolving loops in the kinase domain that probably enabled rapid functional diversification for substrate preferences. One FIKK evolved exclusive tyrosine phosphorylation, previously thought absent in Plasmodium. Despite divergence of substrate preferences, the atypical ATP binding pocket is conserved and we identified a single compound that inhibits all FIKKs. A pan-specific inhibitor could reduce resistance development and improve malaria control strategies.

Positive selection and relaxed purifying selection contribute to rapid evolution of sex-biased genes in green seaweed Ulva

BACKGROUND

The evolution of differences in gamete size and number between sexes is a cornerstone of sexual selection theories. The green macroalga Ulva, with incipient anisogamy and parthenogenetic gametes, provides a unique system to investigate theoretical predictions regarding the evolutionary pressures that drive the transition from isogamy to anisogamy, particularly in relation to gamete size differentiation and sexual selection. Its minimal gamete dimorphism and facultative parthenogenesis enable a rare window into early evolutionary steps toward anisogamy.

RESULTS

By analyzing the expression profiles of sex-biased genes (SBGs) during gametogenesis, we found that SBGs evolve faster than unbiased genes, driven by higher rates of non-synonymous substitution (dN), indicating that SBGs are under stronger selective pressures. Mating type minus-biased genes (mt-BGs) exhibit higher dN/dS values than mating type plus-biased genes (mt+BGs), suggesting stronger selective pressures on mt-BGs, although this difference was not statistically significant (P = 0.08). Using branch-site and RELAX models, we found positive selection and relaxed purifying selection acting on a significant proportion of SBGs, particularly those associated with flagella function.

CONCLUSIONS

This study highlights the selective pressures shaping anisogamy and provides insights into the molecular mechanisms underlying its evolution. The faster evolution of SBGs, particularly mt-BGs, and the positive selection on genes associated with motility, such as those related to flagella function, suggest the importance of enhanced gamete motility in the transition to anisogamy. These findings contribute to our understanding of sexual selection and the evolutionary forces that drive the differentiation of gamete size and number between sexes.

Characterization of five Neisseria homoserine dehydrogenases with diverse coenzyme specificities reveals adaptive evolution of the hom6 genes

Several dehydrogenases using NAD+ or NADP+ as a coenzyme have been characterized, but the molecular evolutionary mechanisms underlying differential coenzyme preferences of these dehydrogenases are largely unknown. In this study, we performed biochemical, kinetic, and phylogenetic analyses of five monofunctional homoserine dehydrogenases (HSDs) encoded by hom6 genes from different Neisseria species (Neisseria elongata, Neisseria animalis, Neisseria dumasiana, Neisseria iguanae, and Neisseria shayeganii) with key amino acid substitutions related to coenzyme specificities. These HSDs were over-expressed in Escherichia coli and purified to apparent homogeneity. Kinetic analysis demonstrated that, in contrast to the NAD+-dependent Neisseria gonorrhoeae HSD (NgHSD), N. elongata HSD (NeHSD) was a NADP+-dependent enzyme, with an approximately 189-fold preference for NADP+ over NAD+. The other four HSDs exhibited NAD+/NADP+ dual coenzyme specificities. Furthermore, other biochemical properties of NeHSD were comparable to NgHSD, the NAD+-dependent homolog. Computational and site-directed mutagenesis studies suggested that Arg45 of NeHSD was a key residue for NADP+ binding. Phylogenetic analysis of Neisseria hom6 genes and positive selection analysis using the branch-site model resulted in the identification of at least four positively selected sites with Bayes empirical Bayes posterior probabilities >0.95. Among these, Leu45 (amino acid numbering according to NgHSD) was implicated in coenzyme specificity. Therefore, we concluded that the coenzyme specificity changes in the HSDs from different Neisseria species were driven by adaptive evolution. These findings significantly advance our understanding regarding the molecular evolution of hom6 genes from Neisseria and provide a foundation for investigating evolutionary mechanisms of coenzyme utilization in other dehydrogenase-coding genes.

Evolutionary trajectories and subfunctionalization of 2 key methyltransferase regulator subfamilies in plants

DNA methylation, a conserved epigenetic mark in both plants and animals, plays a critical role in growth, development, and adaptability. This study explores the origin, evolution, and functional diversification of 2 key methyltransferase regulators, DNAJ-domain-containing protein 1/2/3 (DNAJ1/2/3) and SU(VAR)3-9 HOMOLOG 1/3 (SUVH1/3), in plants. By analyzing genomic data from 21 algae and 86 land plants, we discovered that DNAJ1/2/3 originated within Magnoliopsida, while SUVH1/3 emerged in ferns and evolved through retrotransposition. Both protein families have undergone multiple duplication events and positive selection throughout plant evolution, resulting in their expansion and functional divergence. In dicotyledons, DNAJ1/2/3 diverged into 3 subclades, whereas SUVH1/3 underwent a common duplication event in its ancestral lineage, resulting in 2 subgroups. Structural domain analysis revealed that the evolution of PHD fingers in DNAJ1/2/3 and AT domains in SUVH1/3, under selective pressure, enhanced their interaction capabilities and contributed to the formation of complexes involved in DNA methylation and demethylation regulation. Expression profile analysis across various plant taxa demonstrated tissue-specific expression patterns, with higher expression levels observed in meristematic tissues and active cell regions. These findings elucidate the evolutionary patterns of DNAJ1/2/3 and SUVH1/3 and offer insights into their functional diversification in plants.

The gene characteristics and adaptive evolution of the tumor necrosis factor superfamily (TNFSF) in miiuy croaker, Miichthysmiiuy

The tumor necrosis factor superfamily (TNFSF) is crucial in regulating immune responses, with its members mediating various biological functions through key signaling pathways. However, the gene characteristics of this family and their comparative and evolutionary analysis across species remain limited. In this study, 12 TNFSF genes were identified in the genome-wide of miiuy croaker. Analyses were conducted on evolutionary relationships, conserved motifs, gene duplication, and selection pressure. Conserved motif analyses revealed that the C-terminal motifs of vertebrate TNFSF proteins were more conserved than the N-terminus. Sequence alignment and conservation analysis identified an unrecognized helix structure within the TNF homology domain, which exhibited structural conservation among vertebrates. Synteny and selection pressure analyses indicated that the TNFSF in miiuy croaker exhibited tandem and segmental duplication events. Evolutionary selection pressures may contributed to the functional differentiation of this family. These findings could enhance the understanding of TNFSF gene characteristics and evolutionary relationships, and provide new insights for studying immune-related TNFSF genes.

Selection Pressure Regulates the Evolution-Structure-Function Paradigm of Monocyte Chemoattractant Protein Family

Monocyte chemoattractant proteins (MCPs) are involved in monocyte trafficking during severe inflammation by modulating the chemokine-glycosaminoglycan-receptor signaling axis. MCPs comprise a family of four chemokines (CCL2, CCL7, CCL8, and CCL13/12) that exhibit differential expression patterns in mammals, functional diversity, and receptor/glycosaminoglycan (GAG) binding promiscuity. In this context, the evolution-structure-function paradigm of MCP chemokines in mammals was established by assessing phylogeny, functional divergence, selection pressure, and coevolution in correlation with structural and surface characteristics. Comprehensive analyses were performed using an array of evolutionary and structural bioinformatic methods including molecular dynamics simulations. Our findings demonstrate that substitutions in receptor/GAG-interacting residues mediate episodic diversification and functional diversity in MCP chemokines. Additionally, a balanced interplay of selection pressures has driven the functional changes observed among MCP paralogs, with positive selection at various receptor/GAG-binding sites contributing to their promiscuous receptor/GAG interactions. Meanwhile, processes like purifying selection and coevolution maintain the classical chemokine structure and preserve the ancestral functions of MCP chemokines. Overall, this study suggests that selection pressure on sites within the N-terminal region [N-loop and 310-helix] and 40S loop of MCP chemokines alters surface properties to fine-tune the molecular interactions and functional characteristics without altering the overall chemokine structure.

Characterization of the Pyropia katadae plastid genome and comparative analysis within Bangiales species

Bangiales species live in intertidal regions and suffer from stresses from a variable environment, making them suitable research objects for studying how plants adapt to nature. In this study, the plastid genome of Pyropia katadae was sequenced and compared with those of 11 previously reported Bangiales species. The Py. katadae plastid genome was 193,531 bp long and contained a single-copy region (LSC) of 149,821 bp, a small single-copy region (SSC) of 34,732 bp, and two direct repeats (DRs) 4489 bp long in between. Furthermore, we compared the plastid genomes of 12 Bangiales species. Among the 12 Bangiales plastid genomes, Bangia fuscopurpurea harbored the largest plastid genome size (196,913 bp), and Py. perforata harbored the smallest (189,789 bp). Phylogenetic analyses of shared genes indicated that Py. katadae clustered with Py. yezoensis into a single clade with a high bootstrap value. An overall high degree of similarity in gene content and arrangement among the Bangiales plastid genomes was observed. The size of the plastid genomes and that of the repeats of Pyropia were positively correlated, demonstrating that the repeats were essential for changes in the plastid genome size over a short evolutionary time. The presence of DR or approximate DR regions in most Bangiales plastid genomes indicates the existence of DR regions in their last common ancestor. The different shortened lengths of identical DR regions showed that each species experienced species-specific evolutionary events, which might cause variations in the sequences and the loss of genes. The two steps of fragment reversal could generate the DRs of Bangiales species from an ancestor in common with Florideophyceae. We identified positive selection sites in eight genes that appeared to be essential for Bangiales species to adapt to diverse environments. Our results provide essential genetic data for an in-depth understanding of the evolution and phylogeny of Bangiales species.

Evolution of the Complex Growth Hormone Gene Cluster in Macaques

In higher primates, unlike other mammals, the GH gene locus is complex, comprising several GH-like genes, resulting from gene duplication and divergent evolution, expressed in pituitary and placenta. There are 5 genes in this GH gene cluster in human and 5 to 7 in apes and most Old-World monkeys, but in macaques the cluster has expanded further. Here the nature and evolution of the GH locus in this important primate genus is explored. Analysis of genomic data for Macaca fascicularis (crab-eating macaque) revealed that the GH gene cluster in this species is variable, with at least 5 different haplotypes, comprising 11 to 14 GH-like genes. Gene-number heterozygosity was also detected in Macaca mulatta (rhesus macaque) with 9 to 13 genes. Analysis of genomic data for other macaque species revealed GH gene clusters containing 8 to 14 GH-like genes, but gene-number heterozygosity was not detected. Expression of GH-like genes in pituitary and placenta was examined for Macaca fascicularis. This analysis has established that the complexity of the GH gene cluster increased during the evolution of macaques, by gene duplication and divergent evolution, and that these processes continue within at least 2 extant species. Analysis of rate of sequence change, and distribution of substitutions within the 3D structure, shows that for at least 1 GH-like gene (GH2), the changes reflect positive selection, implying adaptive biological change. Whether this involves changes in physiological (endocrine) function or response to viral or other pathogenic challenge is not yet clear.

Diversification of social complexity following a major evolutionary transition in bees

How social complexity evolved remains a long-standing enigma. In most animal groups, social complexity is typically classified into a few discrete classes. This approach is oversimplified and constrains our inference of social evolution to a narrow trajectory consisting of transitions between classes. Such categorical classifications also limit quantitative studies on the molecular and environmental drivers of social complexity. The recent accumulation of relevant quantitative data has set the stage to overcome these limitations. Here, we propose a data-driven, high-dimensional approach for studying the full diversity of social phenotypes. We curated and analyzed a comprehensive dataset encompassing 17 social traits across 80 species and studied the evolution of social complexity in bees. We found that honey bees, stingless bees, and bumble bees underwent a major evolutionary transition ∼80 mya, inconsistent with the stepwise progression of the social ladder conceptual framework. This major evolutionary transition was followed by a phase of substantial phenotypic diversification of social complexity. Other bee lineages display a continuum of social complexity, ranging from solitary to simple societies, but do not reach the levels of social complexity seen in honey bees, stingless bees, and bumble bees. Bee evolution, therefore, provides a remarkable demonstration of a macroevolutionary process in which a major transition removed biological constraints and opened novel evolutionary opportunities, driving the exploration of the landscape of social phenotypes. Our approach can be extended to incorporate additional data types and readily applied to illuminate the evolution of social complexity in other animal groups.

Evolution of canonical circadian clock genes underlies unique sleep strategies of marine mammals for secondary aquatic adaptation

To satisfy the needs of sleeping underwater, marine mammals, including cetaceans, sirenians, and pinnipeds, have evolved an unusual form of sleep, known as unihemispheric slow-wave sleep (USWS), in which one brain hemisphere is asleep while the other is awake. All aquatic cetaceans have only evolved USWS without rapid eye movement (REM) sleep, whereas aquatic sirenians and amphibious pinnipeds display both bihemispheric slow-wave sleep (BSWS) and USWS, as well as REM sleep. However, the molecular genetic changes underlying USWS remain unknown. The present study investigated the evolution of eight canonical circadian genes and found that positive selection occurred mainly within cetacean lineages. Furthermore, convergent evolution was observed in lineages with USWS at three circadian clock genes. Remarkably, in vitro assays showed that cetacean-specific mutations increased the nuclear localization of zebrafish clocka, and enhanced the transcriptional activation activity of Clocka and Bmal1a. In vivo, transcriptome analysis showed that the overexpression of the cetacean-specific mutant clocka (clocka-mut) caused the upregulation of the wakefulness-promoting glutamatergic genes and the differential expression of multiple genes associated with sleep regulation. In contrast, the GABAergic and cholinergic pathways, which play important roles in promoting sleep, were downregulated in the bmal1a-mut-overexpressing zebrafish. Concordantly, sleep time of zebrafish overexpressing clocka-mut and bmal1a-mut were significantly less than the zebrafish overexpressing the wild-type genes, respectively. These findings support our hypothesis that canonical circadian clock genes may have evolved adaptively to enhance circadian regulation ability relating to sleep in cetaceans and, in turn, contribute to the formation of USWS.

Comparative genomics uncovers evolutionary drivers of locust migratory adaptation

BACKGROUND

Locust migration is one of the main causes of locust plagues. While existing research has highlighted the adaptive migratory capabilities of locusts, the evolutionary patterns of their migration remain elusive. This study aims to explore these evolutionary patterns of locust migratory behavior at the genomic level. To achieve this, we conducted comparative genomics analysis using genomic data from 10 locust species with diverse migratory tendencies.

RESULTS

We identified 1064 genes showing signatures of positive selection in five migratory locust species using a dN/dS model. The BUSTED-PH model revealed 116 genes associated with migratory phenotypes. Gene ontology enrichment analysis indicated that these genes were predominantly related to metabolism and mitochondria-related pathways through both methods. Additionally, the evolutionary rate (RER) analysis between migratory and non-migratory locusts revealed significant divergence in energy metabolism pathways. Notably, of the genes analyzed, the SETX gene consistently showed evidence of positive selection across all five migratory species.

CONCLUSIONS

The findings suggest that the evolution of migratory behavior is associated with increased selective pressure on metabolism and mitochondria-related pathways. Hundreds of genes undergo selective changes during repetitive transitions to migratory behavior. These findings enhance our understanding of the genetic and phenotypic relationships underlying different locust migratory behaviors, providing important data for understanding the biological mechanisms behind locust outbreaks.

Comparative analysis of the mitochondrial genomes of the soft-shelled turtles Palea steindachneri and Pelodiscus axenaria and phylogenetic implications for Trionychia

Soft-shelled turtles, or Trionychia, are an enigmatic and fascinating group due to their specific morphological features and ecological adaptations. Based on mitochondrial DNA (mtDNA) and/or nuclear markers, previous studies showed the incongruent phylogenetic topologies within Trionychia (e.g., the Palea and its closely related species). In order to resolve the equivocal relationships and obtain some "genome-level" common evolutionary characters of soft-shelled turtles, in this study, we assembled and annotated the complete mitochondrial genomes of Palea steindachneri and Pelodiscus axenaria, both naturally distributed in Asia. The sizes of the two mitochondrial genomes were 16,811 bp and 17,143 bp, respectively. Typical vertebrate animal mtDNA features were observed, such as the usual gene components and arrangements (37 genes with a non-coding control region) and the A + T biased nucleotide compositions on the light strand (61.5% and 62.7%, respectively). All conserved blocks common to the vertebrates control region except for the extended terminal associated sequences (ETAS2) were found in the two soft-shelled turtles. The ω ratio averaged over all sites of each protein-coding gene (PCG) was below 1, which indicated purifying selection at the gene-wide level. However, a positive selection site at the 350-codon position in the cytb gene was detected, as estimated by Bayes empirical Bayes (BEB) analysis. Compared with the gene subsets, the mitogenomes provided the most robust phylogenetic resolution. The monophyly of the clades Amydona, Gigantaesuarochelys, and Apalonia was well supported. Topology discrepancies were observed among different datasets (e.g., the positions of Lissemys and Palea), reflecting the heterogeneous phylogenetic signals in the soft-shelled turtle mitogenomes. Precise date estimation based on Bayesian relaxed clock analyses indicated that the crown group age of extant Trionychia was approximately 115.84 Ma (95% HPD: 91.33-142.18 Ma). Paleoclimate changes, especially the Eocene - Oligocene transition, could be responsible for the speciation in these groups. Our results reiterated the necessity and effectiveness of incorporating entire mitochondrial genomes to delineate phylogenetic relationships in chelonian phylogeny studies.

Multi-Omics Profiling of Lipid Variation and Regulatory Mechanisms in Poultry Breast Muscles

This study aimed to elucidate the genetic basis of lipid composition in the breast muscles of poultry, including AA broilers, dwarf guinea fowl, quails, and pigeons, and the impact of artificial selection on lipid traits. By employing lipidomics and transcriptomic sequencing, the research analyzed the chest muscle tissues of these four poultry. A total of 1542 lipid molecules were identified, with 711 showing significant differences among species. These lipids primarily belonged to subclasses such as TG, PC, Phosphatidylethanolamine (PE), Ceramides (Cer), and Diglyceride (DG), with each species demonstrating distinct profiles in these subclasses. Additionally, 5790 orthologous genes were identified, with 763, 767, 24, and 8 genes in AA broilers, dwarf guinea fowl, quails, and pigeons, respectively, exhibiting positive selection (Ka/Ks > 1). Notably, 114 genes related to lipid metabolism displayed significant differential expression, particularly between AA broilers and dwarf guinea fowl. The findings revealed that the metabolic pathways of PC and LPC lipid molecules in the glycerophospholipid pathway, as well as TG lipid molecules in the glycerolipid pathway, exhibited marked interspecies differences, potentially contributing to variations in breast muscle lipid composition. These results provide a solid foundation for understanding the lipid composition and molecular regulatory mechanisms in diverse poultry, offering valuable insights for further research in poultry lipid metabolism and artificial breeding programs.

The Silene latifolia genome and its giant Y chromosome

In many species with sex chromosomes, the Y is a tiny chromosome. However, the dioecious plant Silene latifolia has a giant ~550-megabase Y chromosome, which has remained unsequenced so far. We used a long- and short-read hybrid approach to obtain a high-quality male genome. Comparative analysis of the sex chromosomes with their homologs in outgroups showed that the Y is highly rearranged and degenerated. Recombination suppression between X and Y extended in several steps and triggered a massive accumulation of repeats on the Y as well as in the nonrecombining pericentromeric region of the X, leading to giant sex chromosomes. Using sex phenotype mutants, we identified candidate sex-determining genes on the Y in locations consistent with their favoring recombination suppression events 11 and 5 million years ago.

Writers, readers, and erasers of N6-Methyladenosine (m6A) methylomes in oilseed rape: identification, molecular evolution, and expression profiling

BACKGROUND

m6A RNA modifications are the most prevalent internal modifications in eukaryotic mRNAs and are crucial for plant growth and development, as well as for responses to biotic or abiotic stresses. The modification is catalyzed by writers, removed by erasers, and decoded by various m6A-binding proteins, which are readers. Brassica napus is a major oilseed crop. The dynamic regulation of m6A modifications by writers, erasers, and readers offers potential targets for improving the quality of this crop.

RESULTS

In this study, we identified 92 m6A-regulatory genes in B. napus, including 13 writers, 29 erasers, and 50 readers. A phylogenetic analysis revealed that they could be further divided into four, three, and two clades, respectively. The distribution of protein motifs and gene structures among members of the same clade exhibited notable similarity. During the course of evolution, whole genome duplication (WGD) and segmental duplication were the primary drivers of the expansion of m6A-related gene families. The genes were subjected to rigorous purification selection. Additionally, several sites under positive selection were identified in the proteins. RNA-seq and quantitative real-time PCR (qRT-PCR) expression analyses revealed that the identified Bnam6As exhibit tissue-specific expression patterns, as well as their expression patterns in response to various abiotic and biotic stresses. The 2000 bp sequence upstream of Bnam6As contained a number of cis-acting elements that regulate plant growth and environmental response. Furthermore, the protein interaction network revealed their interactions with a number of proteins of significant functional importance.

CONCLUSION

The identification of m6A modifiers in oilseed rape and their molecular evolution and expression profiling have revealed potential functions and molecular mechanisms of m6A, thus establishing a foundation for further functional validation and molecular breeding.

High-quality genome of black wolfberry (Lycium ruthenicum Murr.) provides insights into the genetics of anthocyanin biosynthesis regulation

Black wolfberry (Lycium ruthenicum Murr.) is an important plant for ecological preservation. In addition, its fruits are rich in anthocyanins and have important edible and medicinal value. However, a high-quality chromosome-level genome for this species is not yet available, and the regulatory mechanisms involved in the biosynthesis of anthocyanins are unclear. In this study, haploid material was used to assemble a high-quality chromosome-level reference genome of Lycium ruthenicum, resulting in a genome size of 2272 Mb with contig N50 of 92.64 Mb, and 38 993 annotated gene models. In addition, the evolution of this genome and large-scale variations compared with the Ningxia wolfberry Lycium barbarum were determined. Importantly, homology annotation identified 86 genes involved in the regulatory pathway of anthocyanin biosynthesis, five of which [LrCHS1 (evm.TU.Chr05.295), LrCHS2 (evm.TU.Chr09.488), LrAOMT (evm.TU.Chr09.809), LrF3'5'H (evm.TU.Chr06.177), and LrAN2.1 (evm.TU.Chr05.2618)] were screened by differential expression analysis and correlation analysis using a combination of transcriptome and metabolome testing. Overexpression of these genes could significantly up- or downregulate anthocyanin-related metabolites. These results will help accelerate the functional genomic research of L. ruthenicum, and the elucidation of the genes involved in anthocyanin synthesis will be beneficial for breeding new varieties and further exploring its ecological conservation potential.

A single gene orchestrates androgen variation underlying male mating morphs in ruffs

Androgens are pleiotropic and play pivotal roles in the formation and variation of sexual phenotypes. We show that differences in circulating androgens between the three male mating morphs in ruff sandpipers are linked to 17-beta hydroxysteroid dehydrogenase 2 (HSD17B2), encoded by a gene within the supergene that determines the morphs. Low-testosterone males had higher HSD17B2 expression in blood than high-testosterone males, as well as in brain areas related to social behaviors and testosterone production. Derived HSD17B2 isozymes, which are absent in high-testosterone males but preferentially expressed in low-testosterone males, converted testosterone to androstenedione faster than the ancestral isozyme. Thus, a combination of evolutionary changes in regulation, sequence, and structure of a single gene introduces endocrine variation underlying reproductive phenotypes.

Genetic characterization of immune adaptor molecule MyD88 in Culex pipiens complex (Diptera: Culicidae) mosquitoes from China

Mosquitoes of the Culex (Cx.) pipiens complex are vectors of severe diseases including West Nile fever by West Nile virus, Japanese encephalitis by Japanese encephalitis virus, and Lymphatic filariasis by filarial nematode Wuchereria bancrofti. As a major portion of mosquito immune system, the Toll pathway implicates in response against infections of mosquito-borne pathogens and biocontrol agents. The genetic diversity of immune-related molecules is expected to be a feasible and effective introduction to expand our knowledge of the mosquito-microbe interplay. However, a comprehensive description is currently lacking regarding the genetic characteristic of the Toll pathway molecules in Cx. pipiens complex mosquitoes. In the present study, genetic changes in Cx. pipiens complex MyD88 (Myeloid differentiation primary response protein 88) were analyzed as a precedent for the Toll pathway molecules in this taxon. MyD88 is a critical adaptor of the pathway transducing signals from TIR-containing receptors to downstream death domain-containing molecules. Our results revealed that adaptive selection has influenced the genetic changes of the molecule, giving rise to acceleration of diversity at a number of amino acid sites. The adaptively selected sites lie in the death domain, intermediate domain, and C-terminal extension. The characteristics of the genetic changes shed insights into the prominent molecular-level structural basis and the involvement strategy of the adaptor in the arms race against exogenous challenges. This finding would be beneficial for further exploration and deeper understanding of the mosquitoes' vectorial capacity and facilitating the effectiveness and sustainability of the biocontrol agents.

Adaptive evolution of stress response genes in parasites aligns with host niche diversity

BACKGROUND

Stress responses are key the survival of parasites and, consequently, also the evolutionary success of these organisms. Despite this importance, our understanding of the evolution of molecular pathways dealing with environmental stressors in parasitic animals remains limited. Here, we tested the link between adaptive evolution of parasite stress response genes and their ecological diversity and species richness. We comparatively investigated antioxidant, heat shock, osmoregulatory, and behaviour-related genes (foraging) in two model parasitic flatworm lineages with contrasting ecological diversity, Cichlidogyrus and Kapentagyrus (Platyhelminthes: Monopisthocotyla), through whole-genome sequencing of 11 species followed by in silico exon bait capture as well as phylogenetic and codon analyses.

RESULTS

We assembled the sequences of 48 stress-related genes and report the first foraging (For) gene orthologs in flatworms. We found duplications of heat shock (Hsp) and oxidative stress genes in Cichlidogyrus compared to Kapentagyrus. We also observed positive selection patterns in genes related to mitochondrial protein import (Hsp) and behaviour (For) in species of Cichlidogyrus infecting East African cichlids-a host lineage under adaptive radiation. These patterns are consistent with a potential adaptation linked to a co-radiation of these parasites and their hosts. Additionally, the absence of cytochrome P450 and kappa and sigma-class glutathione S-transferases in monogenean flatworms is reported, genes considered essential for metazoan life.

CONCLUSIONS

This study potentially identifies the first molecular function linked to a flatworm radiation. Furthermore, the observed gene duplications and positive selection indicate the potentially important role of stress responses for the ecological adaptation of parasite species.

Evolutionary analysis of genes associated with the sense of balance in semi-aquatic mammals

BACKGROUND

Semi-aquatic mammals represent a transitional phase in the evolutionary spectrum between terrestrial and aquatic mammals. The sense of balance is crucial for mammalian locomotion, and in semi-aquatic mammals, the structural foundation of this sense (the vestibular system) shows distinct morphological adaptations to both aquatic and terrestrial environments compared to their terrestrial counterparts. Despite this, the precise molecular mechanisms driving these adaptations remain elusive. Our study endeavors to unravel the genetic components associated with the sense of balance in semi-aquatic mammals and to examine the evolutionary trajectories of these genes, shed light on the molecular mechanisms underlying the adaptive evolution of balance perception in semi-aquatic mammals.

RESULTS

We selected 42 mammal species across 20 orders, 38 families, and 42 genera for analysis. We analyzed a comprehensive set of 116 genes related to the vestibular system's development or function. Our findings indicate that 27 of these genes likely experienced adaptive evolution in semi-aquatic mammals. Particularly, genes such as SLC26A2, SOX10, MYCN, and OTX1 are implicated in collectively orchestrating morphological adaptations in the semicircular canals to suit semi-aquatic environments. Additionally, genes associated with otolith development, including SLC26A2, OC90, and OTOP1, likely regulate otolith sensitivity across various locomotor modes. Moreover, genes linked to vestibular disorders, such as GJB2, GJB6, and USH1C, may provide a molecular foundation for averting vertigo amidst intricate locomotor scenarios in semi-aquatic mammals.

CONCLUSIONS

Our research offers insights into the molecular mechanisms underlying the evolution of the sense of balance in semi-aquatic mammals, while also providing a new research direction for the adaptive evolution of mammals undergoing a secondary transition to an aquatic lifestyle.

Mitochondrial Genome Characteristics Reveal Evolution of Danxiaorchis yangii and Phylogenetic Relationships

Danxiaorchis yangii is a fully mycoheterotrophic orchid that lacks both leaves and roots, belonging to the genus Danxiaorchis in the subtribe Calypsoinae. In this study, we assembled and annotated its mitochondrial genome (397,867 bp, GC content: 42.70%), identifying 55 genes, including 37 protein-coding genes (PCGs), 16 tRNAs, and 2 rRNAs, and conducted analyses of relative synonymous codon usage (RSCU), repeat sequences, horizontal gene transfers (HGTs), and gene selective pressure (dN/dS). Additionally, we sequenced and assembled its plastome, which has a reduced size of 110,364 bp (GC content: 36.60%), comprising 48 PCGs, 26 tRNAs, and 4 rRNAs. We identified 64 potential chloroplast DNA fragments transferred to the mitogenome. Phylogenomic analysis focusing on 33 mitogenomes, with Vitis vinifera as the outgroup, indicated that D. yangii is grouped as follows: D. yangii + ((Dendrobium wilsonii + Dendrobium wilsonii henanense) + Phalaenopsis aphrodite). Phylogenetic analysis based on 83 plastid PCGs from these species showed that D. yangii is grouped as follows: (D. yangii + Pha. aphrodite) + (Den. wilsonii + Den. henanense). Gene selective pressure analysis revealed that most mitochondrial and plastid genes in D. yangii are under purifying selection, ensuring functional stability, and certain genes may have undergone positive selection or adaptive evolution, reflecting the species' adaptation to specific ecological environments. Our study provides valuable data on the plastomes and mitogenomes of D. yangii and lays the groundwork for future research on genetic variation, evolutionary relationships, and the breeding of orchids.

Uncovering selection pressures on the IRF gene family in bats' immune system

Unlike other mammals, bats serve as natural reservoirs for several highly pathogenic viruses without exhibiting symptoms of infection. Recent research has explored the complex mechanisms underlying the balance between bats' antiviral defenses and their pathological responses. However, the evolution of the molecular drivers behind bats' antiviral strategies remains largely unknown. Interferon regulatory factors (IRFs) are essential transcription factors that bind to DNA and regulate the expression of numerous genes involved in antiviral defense, inflammation, immune cell differentiation, apoptosis, and oncogenesis. Our research focused on members of the IRF family, using 17 bat species and four terrestrial mammals available in GenBank. We employed CodeML to detect signs of positive selection through three different models. Statistically significant results were obtained for the IRF-1, IRF-4, IRF-5, IRF-6, and IRF-9 genes, which are known to play pivotal roles in various regulation mechanisms. Specifically, IRF-4 and IRF-5 are key in modulating the inflammatory response, while IRF-1 is essential for antiviral defense in bats, and IRF-9 regulates genes activated by type I interferon. Although the role of IRF-6 in these mechanisms requires further investigation in bats, all these genes show signs of positive selection, suggesting an optimization of the processes they regulate. These findings highlight the adaptive role of IRF elements in enhancing, among other things, the bat immune system, potentially improving their resilience and efficacy. Our study not only provides new genetic insights into bats but also underscores the remarkable molecular evolution within this unique group of mammals.

Positive Selection on Mammalian Immune Genes-Effects of Gene Function and Selective Constraint

Genome-wide analyses of various taxa have repeatedly shown that immune genes are important targets of positive selection. However, little is known about what factors determine which immune genes are under positive selection. To address this question, we here focus on the mammalian immune system and investigate the importance of gene function and other factors such as gene expression, protein-protein interactions, and overall selective constraint as determinants of positive selection. We compiled a list of >1,100 immune genes that were divided into six functional categories and analyzed using data from rodents. Genes encoding proteins that are in direct interactions with pathogens, such as pattern recognition receptors (PRRs), are often expected to be key targets of positive selection. We found that categories containing cytokines, cytokine receptors, and other cell surface proteins involved in, for example, cell-cell interactions were at least as important targets as PRRs, with three times higher rate of positive selection than nonimmune genes. The higher rate of positive selection on cytokines and cell surface proteins was partly an effect of these categories having lower selective constraint. Nonetheless, cytokines had a higher rate of positive selection than nonimmune genes even at a given level of selective constraint, indicating that gene function per se can also be a determinant of positive selection. These results have broad implications for understanding the causes of positive selection on immune genes, specifically the relative importance of host-pathogen coevolution versus other processes.

Diversification and recurrent adaptation of the synaptonemal complex in Drosophila

The synaptonemal complex (SC) is a protein-rich structure essential for meiotic recombination and faithful chromosome segregation. Acting like a zipper to paired homologous chromosomes during early prophase I, the complex is a symmetrical structure where central elements are connected on two sides by the transverse filaments to the chromatin-anchoring lateral elements. Despite being found in most major eukaryotic taxa implying a deeply conserved evolutionary origin, several components of the complex exhibit unusually high rates of sequence turnover. This is puzzlingly exemplified by the SC of Drosophila, where the central elements and transverse filaments display no identifiable homologs outside of the genus. Here, we exhaustively examine the evolutionary history of the SC in Drosophila taking a comparative phylogenomic approach with high species density to circumvent obscured homology due to rapid sequence evolution. Contrasting starkly against other genes involved in meiotic chromosome pairing, SC genes show significantly elevated rates of coding evolution due to a combination of relaxed constraint and recurrent, widespread positive selection. In particular, the central element cona and transverse filament c(3)G have diversified through tandem and retro-duplications, repeatedly generating paralogs with novel germline activity. In a striking case of molecular convergence, c(3)G paralogs that independently arose in distant lineages evolved under positive selection to have convergent truncations to the protein termini and elevated testes expression. Surprisingly, the expression of SC genes in the germline is prone to change suggesting recurrent regulatory evolution which, in many species, resulted in high testes expression even though Drosophila males are achiasmic. Overall, our study recapitulates the poor conservation of SC components, and further uncovers that the lack of conservation extends to other modalities including copy number, genomic locale, and germline regulation. Considering the elevated testes expression in many Drosophila species and the common ancestor, we suggest that the activity of SC genes in the male germline, while still poorly understood, may be a prime target of constant evolutionary pressures driving repeated adaptations and innovations.

Lineage classification and selective site identification of Orthoebolavirus zairense

As the high pathogenic species of Filoviridae virus family, Orthoebolavirus zairense (EBOV) shows frequent outbreaks in human in recently years since its first emerging in 1976 in Democratic Republic of the Congo (COD), bringing ongoing risks and burden on public health safety. Here, the phylogenetic relationship among major outbreaks was analyzed. The results showed that EBOV isolates could be divided into four lineages according to spatial and temporal epidemics. Then, the positive selection sites (PSSs) were detected on all proteins of the EBOV, exhibiting lineage characteristic. Particularly, sites in GP and VP24 were identified to be significantly under positive selection, and partial of which were maintained in the latest isolates in 2021. GP and L were found to have high variability between lineages. Substitutions including F443L and F443S in GP, as well as F1610L and I1951V in L could be characteristic of the two large outbreaks in COD (2018) and West Africa (2014), respectively. Further, substitutions of significant PSSs in VP24 and L proteins were visualized for analysis of structural changes, which may affect EBOV pathogenesis. In summary, our results gains insights in genetic characteristic and adaptive evolution of EBOV, which could facilitate gene functional research against EBOV.

Into the Blue: Exploring genetic mechanisms behind the evolution of baleen whales

Marine ecosystems are ideal for studying evolutionary adaptations involved in lineage diversification due to few physical barriers and reduced opportunities for strict allopatry compared to terrestrial ecosystems. Cetaceans (whales, dolphins, and porpoises) are a diverse group of mammals that successfully adapted to various habitats within the aquatic environment around 50 million years ago. While the overall adaptive transition from terrestrial to fully aquatic species is relatively well understood, the radiation of modern whales is still unclear. Here high-quality genomes derived from previously published data were used to identify genomic regions that potentially underpinned the diversification of baleen whales (Balaenopteridae). A robust molecular phylogeny was reconstructed based on 10,159 single copy and complete genes for eight mysticetes, seven odontocetes and two cetacean outgroups. Analysis of positive selection across 3,150 genes revealed that balaenopterids have undergone numerous idiosyncratic and convergent genomic variations that may explain their diversification. Genes associated with aging, survival and homeostasis were enriched in all species. Additionally, positive selection on genes involved in the immune system were disclosed for the two largest species, blue and fin whales. Such genes can potentially be ascribed to their morphological evolution, allowing them to attain greater length and increased cell number. Further evidence is presented about gene regions that might have contributed to the extensive anatomical changes shown by cetaceans, including adaptation to distinct environments and diets. This study contributes to our understanding of the genomic basis of diversification in baleen whales and the molecular changes linked to their adaptive radiation, thereby enhancing our understanding of cetacean evolution.

Animal models in neuroscience with alternative approaches: Evolutionary, biomedical, and ethical perspectives

Animal models have been a crucial tool in neuroscience research for decades, providing insights into the biomedical and evolutionary mechanisms of the nervous system, disease, and behavior. However, their use has raised concerns on several ethical, clinical, and scientific considerations. The welfare of animals and the 3R principles (replacement, reduction, refinement) are the focus of the ethical concerns, targeting the importance of reducing the stress and suffering of these models. Several laws and guidelines are applied and developed to protect animal rights during experimenting. Concurrently, in the clinic and biomedical fields, discussions on the relevance of animal model findings on human organisms have increased. Latest data suggest that in a considerable amount of time the animal model results are not translatable in humans, costing time and money. Alternative methods, such as in vitro (cell culture, microscopy, organoids, and micro physiological systems) techniques and in silico (computational) modeling, have emerged as potential replacements for animal models, providing more accurate data in a minimized cost. By adopting alternative methods and promoting ethical considerations in research practices, we can achieve the 3R goals while upholding our responsibility to both humans and other animals. Our goal is to present a thorough review of animal models used in neuroscience from the biomedical, evolutionary, and ethical perspectives. The novelty of this research lies in integrating diverse points of views to provide an understanding of the advantages and disadvantages of animal models in neuroscience and in discussing potential alternative methods.

Molecular evolution and genetic diversity of defective chorion 1 in Anastrepha fraterculus and Anastrepha obliqua (Diptera, Tephritidae)

The family Tephritidae comprises numerous fruit fly species, some of which are economically significant, such as several in the genus Anastrepha. Most pest species in this genus belong to the fraterculus group, characterized by closely related species that are difficult to differentiate due to recent divergence and gene flow. Identifying genetic markers for their study is paramount for understanding the group's evolution and eventual phytosanitary control. Because there is variation in eggshell morphology among species in the genus, the study of the rapidly evolving defective chorion 1 (dec-1) gene, which is crucial for chorion formation and reproduction, could provide relevant information for Anastrepha differentiation. We compared transcriptome sequences of dec-1 from two of the most important pest species in the genus, Anastrepha fraterculus and Anastrepha obliqua to dec-1 sequences from Anastrepha ludens, which was used for structure prediction. Furthermore, we amplified a conserved exon across populations of these species. These data revealed three alternative transcripts in A. fraterculus and A. obliqua, consistent with patterns found in other Tephritidae; we obtained orthologous sequences for these other tephritids from NCBI to investigate patterns of selection affecting this gene at different hierarchical levels using different methods. These analyses show a general pattern of purifying selection across the whole gene and throughout its history at different hierarchical levels, from populations to more distantly related species. That notwithstanding, we still found evidence of positive and episodic diversifying selection at different levels. Different parts of the gene have shown distinct evolutionary rates, which were associated with the diverse proproteins produced by posttranslational changes of DEC-1, with proproteins that are incorporated in the chorion earlier in egg formation being in general more conserved than others that are incorporated later. This correlation appears more evident in certain lineages, including the branch that separates Anastrepha, as well as other internal branches that differentiate species within the genus. Our data showed that this gene shows remarkable variation across its different exons, which has proven to be informative at different evolutionary levels. These changes hold promise not only for studying differentiation in Anastrepha but also for the eventual management of selected pest species.

Reticulate evolution: Detection and utility in the phylogenomics era

Phylogenomics has enriched our understanding that the Tree of Life can have network-like or reticulate structures among some taxa and genes. Two non-vertical modes of evolution - hybridization/introgression and horizontal gene transfer - deviate from a strictly bifurcating tree model, causing non-treelike patterns. However, these reticulate processes can produce similar patterns to incomplete lineage sorting or recombination, potentially leading to ambiguity. Here, we present a brief overview of a phylogenomic workflow for inferring organismal histories and compare methods for distinguishing modes of reticulate evolution. We discuss how the timing of coalescent events can help disentangle introgression from incomplete lineage sorting and how horizontal gene transfer events can help determine the relative timing of speciation events. In doing so, we identify pitfalls of certain methods and discuss how to extend their utility across the Tree of Life. Workflows, methods, and future directions discussed herein underscore the need to embrace reticulate evolutionary patterns for understanding the timing and rates of evolutionary events, providing a clearer view of life's history.

Phylogenomics and adaptive evolution of hydrophytic umbellifers (tribe Oenantheae, Apioideae) revealed from chloroplast genomes

BACKGROUND

Tribe Oenantheae consists mainly of aquatic species within the Apioideae. The unique morphology and habitat distinguish this group from other Apioideae groups. However, the genomic information of these group species has not been widely developed, and the molecular mechanisms of adaptive evolution remain unclear.

RESULTS

We provide comparative analyses on 30 chloroplast genomes of this tribe representing five genera to explore the molecular variation response to plant adaptations. The Oenantheae chloroplast genomes presented typical quadripartite structures, with sizes ranging from 153,024 bp to 155,006 bp. Gene content and order were highly conserved with no significant expansion or contraction observed. Seven regions (rps16 intron-trnK, rpoB-trnC, trnE-trnT-psbD, petA-psbJ, ndhF-rpl32-trnL, ycf1a-rps15, and ycf1a gene) were identified as remarkable candidate DNA markers for future studies on species identification, biogeography, and phylogeny of tribe Oenantheae. Our study elucidated the relationships among the genera of tribe Oenantheae and subdivided the genera of Sium and Oenanthe. However, relationships among the Oenanthe I clade remain to be further clarified. Eight positively selected genes (accD, rbcL, rps8, ycf1a, ycf1b, ycf2, ndhF, and ndhK) were persuasively detected under site models tests, and these genes might have played roles in Oenantheae species adaptation to the aquatic environments.

CONCLUSIONS

Our results provide sufficient molecular markers for the subsequent molecular studies of the tribe Oenantheae, and promote the understanding of the adaptation of the Oenantheae species to aquatic environments.

Phylogenetics of Lepidonotopodini (Macellicephalinae, Polynoidae, Annelida) and Comparative Mitogenomics of Shallow-Water vs. Deep-Sea Scaleworms (Aphroditiformia)

Within Polynoidae, a diverse aphroditiform family, the subfamily Macellicephalinae comprises anchialine cave-dwelling and deep-sea scaleworms. In this study, Lepidonotopodinae is synonymized with Macellicephalinae, and the tribe Lepidonotopodini is applied to a well-supported clade inhabiting deep-sea chemosynthetic-based ecosystems. Newly sequenced "genome skimming" data for 30 deep-sea polynoids and the comparatively shallow living Eulagisca gigantea is used to bioinformatically assemble their mitogenomes. When analyzed with existing scaleworm mitogenomes, deep-sea scaleworms exhibit increased gene order rearrangement events compared to shallow-water relatives. Additionally, comparative analyses of shallow-water vs. deep-sea polynoid substitution rates in mitochondrial protein-coding genes show an overall relaxed purifying selection and a positive selection of several amino acid sites in deep-sea species, indicating that polynoid mitogenomes have undergone selective pressure to evolve metabolic adaptations suited to deep-sea environments. Furthermore, the inclusion of skimming data for already known Lepidonotopodini species allowed for an increased coverage of DNA data and a representation of the taxa necessary to create a more robust phylogeny using 18 genes, as opposed to the six genes previously used. The phylogenetic results support the erection of Cladopolynoe gen. nov., Mamiwata gen. nov., Photinopolynoe gen. nov., Stratigos gen. nov., and Themis gen. nov., and emended diagnoses for Branchinotogluma, Branchipolynoe, Lepidonotopodium, and Levensteiniella.

Deciphering the complex organelle genomes of two Rhododendron species and insights into adaptive evolution patterns in high-altitude

BACKGROUND

The genomes within organelles are crucial for physiological functions such as respiration and photosynthesis and may also contribute to environmental adaptation. However, the limited availability of genetic resources, particularly mitochondrial genomes, poses significant challenges for in-depth investigations.

RESULTS

Here, we explored various assembly methodologies and successfully reconstructed the complex organelle genomes of two Rhododendron species: Rhododendron nivale subsp. boreale and Rhododendron vialii. The mitogenomes of these species exhibit various conformations, as evidenced by long-reads mapping. Notably, only the mitogenome of R. vialii can be depicted as a singular circular molecule. The plastomes of both species conform to the typical quadripartite structure but exhibit elongated inverted repeat (IR) regions. Compared to the high similarity between plastomes, the mitogenomes display more obvious differences in structure, repeat sequences, and codon usage. Based on the analysis of 58 organelle genomes from angiosperms inhabiting various altitudes, we inferred the genetic adaptations associated with high-altitude environments. Phylogenetic analysis revealed partial inconsistencies between plastome- and mitogenome-derived phylogenies. Additionally, evolutionary lineage was determined to exert a greater influence on codon usage than altitude. Importantly, genes such as atp4, atp9, mttB, and clpP exhibited signs of positive selection in several high-altitude species, suggesting a potential link to alpine adaptation.

CONCLUSIONS

We tested the effectiveness of different organelle assembly methods for dealing with complex genomes, while also providing and validating high-quality organelle genomes of two Rhododendron species. Additionally, we hypothesized potential strategies for high-altitude adaptation of organelles. These findings offer a reference for the assembly of complex organelle genomes, while also providing new insights and valuable resources for understanding their adaptive evolution patterns.

Selection on synonymous codon usage in soybean (Glycine max) WRKY genes

The WRKY transcription factor gene family in soybean [Glycine max (L.) Merr.] (GmWRKY) is critical for the plant's development and stress responses. This study examines the evolutionary dynamics of the GmWRKY gene family, focusing on its synonymous codon usage bias (CUB) in a comprehensive set of 179 coding sequences. CUB was analyzed using various indices, revealing a preference for A/T-ending codons and relatively low codon bias. Codon adaptation index (CAI) analysis suggested that these genes are optimized for efficient translation despite relatively low bias, reflecting a balance between codon diversity and translation efficiency. Neutrality and NC plots indicated that selective forces dominate over mutational forces in shaping codon usage, while selection signature analysis showed purifying selection being prevalent across the gene family. However, episodic positive selection was also detected in certain clades, highlighting potential adaptive diversification in response to environmental stress. Additionally, promoter binding site analysis uncovered correlations between codon usage and transcriptional regulation, indicating a context-dependent relationship between CUB and gene expression. Phylogenetic analysis identified 11 well-supported clades in the modern GmWRKY gene family and ancestral sequence reconstruction revealed more relaxed codon preferences and reduced selection constraints in modern GmWRKY genes, potentially linked to neofunctionalization and adaptation to environmental changes. These findings provide a framework for optimizing gene expression in transgenic soybean crops with resilience. Further functional validation of positively selected genes is recommended to elucidate their role in stress responses.

Increased Evolutionary Rate in the Z chromosome of Sympatric and Allopatric Species of Morpho Butterflies

Divergent evolution of genomes among closely related species is shaped by both neutral processes and ecological forces, such as local adaptation and reinforcement. These factors can drive accelerated evolution of sex chromosomes relative to autosomes. Comparative genomic analyses between allopatric and sympatric species with overlapping or divergent ecological niches offer insights into reinforcement and ecological specialization on genome evolution. In the butterfly genus Morpho, several species coexist in sympatry, with specialization across forest strata and temporal niches. We analyzed the genomes of eight Morpho species, along with previously published genomes of three others, to compare chromosomal rearrangements and signs of positive selection in the Z chromosome vs. autosomes. We found extensive chromosomal rearrangements in Z chromosome, particularly in sympatric species with similar ecological niches, suggesting a role for inversions in restricting gene flow at a postzygotic level. Z-linked genes also exhibited significantly higher dN/dS ratios than autosomal genes across the genus, with pronounced differences in closely related species living in sympatry. Additionally, we examined the evolution of eight circadian clock genes, detecting positive selection in Period, located on the Z chromosome. Our findings suggest that the Z chromosome evolves more rapidly than autosomes, particularly among closely related species, raising questions about its role in prezygotic and postzygotic isolation mechanisms.

Lost in translation? Evidence for a muted proteomic response to thermal stress in a stenothermal Antarctic fish and possible evolutionary mechanisms

Stenothermal Antarctic notothenioid fishes are noteworthy for their history of isolation in extreme cold and their corresponding lack of the canonical heat shock response. Despite extensive transcriptomic studies, the mechanistic basis for stenothermy has not been fully elucidated. Given that the proteome better represents an organism's physiology, the possibility exists that some aspects of stenothermy arise posttranscriptionally. Here, Antarctic emerald rockcod (Trematomus bernacchii) were sampled after exposure to chronic and/or acute high temperatures, followed by a thorough assessment of proteomic responses in the brain, gill, and kidney. Few cellular stress response proteins were induced, and overall responses were modest in terms of the numbers of differentially expressed proteins and their fold changes. Inconsistencies in protein induction across treatments and tissues are suggestive of dysregulation, rather than an adaptive response. Changes in regulation of the translational machinery in Antarctic notothenioids could explain these patterns. Some components of translational regulatory pathways are highly conserved [e.g., Ser-52, eukaryotic translation initiation factor 2α (eIF2α)], but other proteins comprising the cellular "integrated stress response," specifically, the eIF2α kinases general control nonderepressible 2 (GCN2) and PKR-like endoplasmic reticulum kinase (PERK), may have evolved along different trajectories in Antarctic fishes. Taken together, these observations suggest a novel hypothesis for stenothermy and the absence of a coordinated cellular stress response in Antarctic fishes.NEW & NOTEWORTHY Antarctic fishes have some of the lowest known heat tolerances among vertebrates, but the molecular mechanisms underlying this pattern are not fully understood. By combining detailed analyses of protein expression patterns in several tissues under various heat treatments with a broader evolutionary perspective, this study offers a novel hypothesis to explain the narrow range of temperature tolerance in this extraordinary group of fishes.

An integrative genomic toolkit for studying the genetic, evolutionary, and molecular underpinnings of eusociality in insects

While genomic resources for social insects have vastly increased over the past two decades, we are still far from understanding the genetic and molecular basis of eusociality. Here, we briefly review three scientific advancements that, when integrated, can be highly synergistic for advancing our knowledge of the genetics and evolution of eusocial traits. Population genomics provides a natural way to quantify the strength of natural selection on coding and regulatory sequences, highlighting genes that have undergone adaptive evolution during the evolution or maintenance of eusociality. Genome-wide association studies (GWAS) can be used to characterize the complex genetic architecture underlying eusocial traits and identify candidate causal variants. Concurrently, CRISPR/Cas9 enables the precise manipulation of gene function to both validate genotype-phenotype associations and study the molecular biology underlying interesting traits. While each approach has its own advantages and disadvantages, which we discuss herein, we argue that their combination will ultimately help us better understand the genetics and evolution of eusocial behavior. Specifically, by triangulating across these three different approaches, researchers can directly identify and study loci that have a causal association with key phenotypes and have evidence of positive selection over the relevant timescales associated with the evolution and maintenance of eusociality in insects.

Evolution of Sensory Receptors

Sensory receptors are at the interface between an organism and its environment and thus represent key sites for biological innovation. Here, we survey major sensory receptor families to uncover emerging evolutionary patterns. Receptors for touch, temperature, and light constitute part of the ancestral sensory toolkit of animals, often predating the evolution of multicellularity and the nervous system. In contrast, chemoreceptors exhibit a dynamic history of lineage-specific expansions and contractions correlated with the disparate complexity of chemical environments. A recurring theme includes independent transitions from neurotransmitter receptors to sensory receptors of diverse stimuli from the outside world. We then provide an overview of the evolutionary mechanisms underlying sensory receptor diversification and highlight examples where signatures of natural selection are used to identify novel sensory adaptations. Finally, we discuss sensory receptors as evolutionary hotspots driving reproductive isolation and speciation, thereby contributing to the stunning diversity of animals.

Genomic analysis and mechanisms exploration of a stress tolerance and high-yield pullulan producing strain

Pullulan is a kind of natural polymer, which is widely used in medicine and food because of its solubility, plasticity, edible, non-toxicity and good biocompatibility. It is of great significance to improve the yield of pullulan by genetic modification of microorganisms. It was previously reported that Aureobasidium melanogenum TN3-1 isolated from honey-comb could produce high-yield of pullulan, but the molecular mechanisms of its production of pullulan had not been completely solved. In this study, the reported strains of Aureobasidium spp. were further compared and analyzed at genome level. It was found that genome duplication and genome genetic variations might be the crucial factors for the high yield of pullulan and stress resistance. This particular phenotype may be the result of adaptive evolution, which can adapt to its environment through genetic variation and adaptive selection. In addition, the TN3-1 strain has a large genome, and the special regulatory sequences of its specific genes and promoters may ensure a unique characteristics. This study is a supplement of the previous studies, and provides basic data for the research of microbial genome modification in food and healthcare applications.

Evolutionary characterization of the establishment of H6 influenza viruses in domestic geese in China: implications for the position of the host in the ecosystem

Geese, both wild and domestic, are generally considered part of the natural reservoir for influenza A viruses. The highly pathogenic H5 Goose/Guangdong avian influenza virus lineage that is still causing outbreaks worldwide was first detected in domestic geese in 1996. However, while wild geese might have a somewhat restricted role in the influenza ecosystem, the role of domestic geese is little studied. Here, 109 H6 viruses isolated from domestic geese during 2001-2018 in southern China had their phylogeny, evolutionary dynamics, and molecular signatures characterized to examine the role of domestic geese. Our findings demonstrated that all geese H6 viruses were derived from H6 viruses established in ducks and that they subsequently formed three distinct hemagglutinin lineages. Rapid evolution of the hemagglutinin genes was not detected after the duck-to-goose transmissions of H6 viruses that then circulated in geese. Despite long-term persistence in geese, H6 viruses were rarely observed to transmit back to ducks or terrestrial poultry and never exchanged genes with viruses from other subtypes. Most geese H6 viruses maintained the primary molecular signatures of their duck precursors. This study raises the possibility that, rather than being part of the natural reservoir, domestic geese might be more like an aberrant host species for influenza A viruses, and perhaps a "dead-end" host.

Independent repeated mutations within the alphaviruses Ross River virus and Barmah Forest virus indicates convergent evolution and past positive selection in ancestral populations despite ongoing purifying selection

Ross River virus (RRV) and Barmah Forest virus (BFV) are arthritogenic arthropod-borne viruses (arboviruses) that exhibit generalist host associations and share distributions in Australia and Papua New Guinea (PNG). Using stochastic mapping and discrete-trait phylogenetic analyses, we profiled the independent evolution of RRV and BFV signature mutations. Analysis of 186 RRV and 88 BFV genomes demonstrated their viral evolution trajectories have involved repeated selection of mutations, particularly in the nonstructural protein 1 (nsP1) and envelope 3 (E3) genes suggesting convergent evolution. Convergent mutations in the nsP1 genes of RRV (residues 248 and 441) and BFV (residues 297 and 447) may be involved with catalytic enzyme mechanisms and host membrane interactions during viral RNA replication and capping. Convergent E3 mutations (RRV site 59 and BFV site 57) may be associated with enzymatic furin activity and cleavage of E3 from protein precursors assisting viral maturation and infectivity. Given their requirement to replicate in disparate insect and vertebrate hosts, convergent evolution in RRV and BFV may represent a dynamic link between their requirement to selectively 'fine-tune' intracellular host interactions and viral replicative enzymatic processes. Despite evidence of evolutionary convergence, selection pressure analyses did not reveal any RRV or BFV amino acid sites under strong positive selection and only weak positive selection for nonstructural protein sites. These findings may indicate that their alphavirus ancestors were subject to positive selection events which predisposed ongoing pervasive convergent evolution, and this largely supports continued purifying selection in RRV and BFV populations during their replication in mosquito and vertebrate hosts.

Reversions mask the contribution of adaptive evolution in microbiomes

When examining bacterial genomes for evidence of past selection, the results depend heavily on the mutational distance between chosen genomes. Even within a bacterial species, genomes separated by larger mutational distances exhibit stronger evidence of purifying selection as assessed by dN/dS, the normalized ratio of nonsynonymous to synonymous mutations. Here, we show that the classical interpretation of this scale dependence, weak purifying selection, leads to problematic mutation accumulation when applied to available gut microbiome data. We propose an alternative, adaptive reversion model with opposite implications for dynamical intuition and applications of dN/dS. Reversions that occur and sweep within-host populations are nearly guaranteed in microbiomes due to large population sizes, short generation times, and variable environments. Using analytical and simulation approaches, we show that adaptive reversion can explain the dN/dS decay given only dozens of locally fluctuating selective pressures, which is realistic in the context of Bacteroides genomes. The success of the adaptive reversion model argues for interpreting low values of dN/dS obtained from long timescales with caution as they may emerge even when adaptive sweeps are frequent. Our work thus inverts the interpretation of an old observation in bacterial evolution, illustrates the potential of mutational reversions to shape genomic landscapes over time, and highlights the importance of studying bacterial genomic evolution on short timescales.

Rapid Loss of Nutritional Symbionts in an Endemic Hawaiian Herbivore Radiation Is Associated with Plant Galling Habit

Insect herbivores frequently cospeciate with symbionts that enable them to survive on nutritionally unbalanced diets. While ancient symbiont gain and loss events have been pivotal for insect diversification and feeding niche specialization, evidence of recent events is scarce. We examine the recent loss of nutritional symbionts (in as little as 1 MY) in sap-feeding Pariaconus, an endemic Hawaiian insect genus that has undergone adaptive radiation, evolving various galling and free-living ecologies on a single host-plant species, Metrosideros polymorpha within the last ∼5 MY. Using 16S rRNA sequencing, we investigated the bacterial microbiomes of 19 Pariaconus species and identified distinct symbiont profiles associated with specific host-plant ecologies. Phylogenetic analyses and metagenomic reconstructions revealed significant differences in microbial diversity and functions among psyllids with different host-plant ecologies. Within a few millions of years, Pariaconus species convergently evolved the closed-gall habit twice. This shift to enclosed galls coincided with the loss of the Morganella-like symbiont that provides the essential amino acid arginine to free-living and open-gall sister species. After the Pariaconus lineage left Kauai and colonized younger islands, both open- and closed-gall species lost the Dickeya-like symbiont. This symbiont is crucial for synthesizing essential amino acids (phenylalanine, tyrosine, and lysine) as well as B vitamins in free-living species. The recurrent loss of these symbionts in galling species reinforces evidence that galls are nutrient sinks and, combined with the rapidity of the evolutionary timeline, highlights the dynamic role of insect-symbiont relationships during the diversification of feeding ecologies. We propose new Candidatus names for the novel Morganella-like and Dickeya-like symbionts.

Diverse signatures of convergent evolution in cactus-associated yeasts

Many distantly related organisms have convergently evolved traits and lifestyles that enable them to live in similar ecological environments. However, the extent of phenotypic convergence evolving through the same or distinct genetic trajectories remains an open question. Here, we leverage a comprehensive dataset of genomic and phenotypic data from 1,049 yeast species in the subphylum Saccharomycotina (Kingdom Fungi, Phylum Ascomycota) to explore signatures of convergent evolution in cactophilic yeasts, ecological specialists associated with cacti. We inferred that the ecological association of yeasts with cacti arose independently approximately 17 times. Using a machine learning-based approach, we further found that cactophily can be predicted with 76% accuracy from both functional genomic and phenotypic data. The most informative feature for predicting cactophily was thermotolerance, which we found to be likely associated with altered evolutionary rates of genes impacting the cell envelope in several cactophilic lineages. We also identified horizontal gene transfer and duplication events of plant cell wall-degrading enzymes in distantly related cactophilic clades, suggesting that putatively adaptive traits evolved independently through disparate molecular mechanisms. Notably, we found that multiple cactophilic species and their close relatives have been reported as emerging human opportunistic pathogens, suggesting that the cactophilic lifestyle-and perhaps more generally lifestyles favoring thermotolerance-might preadapt yeasts to cause human disease. This work underscores the potential of a multifaceted approach involving high-throughput genomic and phenotypic data to shed light onto ecological adaptation and highlights how convergent evolution to wild environments could facilitate the transition to human pathogenicity.

Evolution of rhodopsin in flatfishes (Pleuronectiformes) is associated with depth and migratory behavior

Visual signals are involved in many fitness-related tasks and are therefore essential for survival in many species. Aquatic organisms are ideal systems to study visual evolution, as the high diversity of spectral properties in aquatic environments generates great potential for adaptation to different light conditions. Flatfishes are an economically important group, with over 800 described species distributed globally, including halibut, flounder, sole, and turbot. The diversity of flatfish species and wide array of environments they occupy provides an excellent opportunity to understand how this variation translates to molecular adaptation of vision genes. Using models of molecular evolution, we investigated how the light environments inhabited by different flatfish lineages have shaped evolution in the rhodopsin gene, which is responsible for mediating dim-light visual transduction. We found strong evidence for positive selection in rhodopsin, and this was correlated with both migratory behavior and several fundamental aspects of habitat, including depth and freshwater/marine evolutionary transitions. We also identified several mutations that likely affect the wavelength of peak absorbance of rhodopsin, and outline how these shifts in absorbance correlate with the response to the light spectrum present in different habitats. This is the first study of rhodopsin evolution in flatfishes that considers their extensive diversity, and our results highlight how ecologically-driven molecular adaptation has occurred across this group in response to transitions to novel light environments.

Genomic Gigantism is not Associated with Reduced Selection Efficiency in Neotropical Salamanders

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.

Distinct prophage gene profiles of Staphylococcus aureus strains from atopic dermatitis patients and healthy individuals

Staphylococcus aureus strains exhibit varying associations with atopic dermatitis (AD), but the genetic determinants underpinning the pathogenicity are yet to be fully characterized. To reveal the genetic differences between S. aureus strains from AD patients and healthy individuals (HE), we developed and employed a random forest classifier to identify potential marker genes responsible for their phenotypic variations. The classifier was able to effectively distinguish strains from AD and HE. We also uncovered strong links between certain marker genes and phage functionalities, with phage holin emerging as the most pivotal differentiating factor. Further examination of S. aureus gene content highlighted the genetic diversity and functional implications of prophages in driving differentiation between strains from AD and HE. The HE group exhibited greater gene content diversity, largely influenced by their prophages. While strains from both AD and HE universally housed prophages, those in the HE group were distinctively higher at the strain level. Moreover, although prophages in the HE group exhibited variously higher enrichment of differential functions, the AD group displayed a notable enrichment of virulence factors within their prophages, underscoring the important contribution of prophages to the pathogenesis of AD-associated strains. Overall, prophages significantly shape the genetic and functional profiles of S. aureus strains, shedding light on their pathogenic potential and elucidating the mechanisms behind the phenotypic variations in AD and HE environments.

IMPORTANCE

Through a nuanced exploration of Staphylococcus aureus strains obtained from atopic dermatitis (AD) patients and healthy controls (HE), our research unveils pivotal genetic determinants influencing their pathogenic associations. Utilizing a random forest classifier, we illuminate distinct marker genes, with phage holin emerging as a critical differential factor, revealing the profound impact of prophages on genetic and pathogenic profiles. HE strains exhibited a diverse gene content, notably shaped by unique, heightened prophages. Conversely, AD strains emphasized a pronounced enrichment of virulence factors within prophages, signifying their key role in AD pathogenesis. This work crucially highlights prophages as central architects of the genetic and functional attributes of S. aureus strains, providing vital insights into pathogenic mechanisms and phenotypic variations, thereby paving the way for targeted AD therapeutic approaches and management strategies by demystifying specific genetic and pathogenic mechanisms.

Comparative Genomic Analysis of the Pattern of Evolution of Male and Female Reproductive Proteins in Seed Beetles

Male seminal fluid proteins often show signs of positive selection and divergent evolution, believed to reflect male-female coevolution. Yet, our understanding of the predicted concerted evolution of seminal fluid proteins and female reproductive proteins is limited. We sequenced, assembled, and annotated the genome of two species of seed beetles allowing a comparative analysis of four closely related species of these herbivorous insects. We compare the general pattern of evolution in genes encoding seminal fluid proteins and female reproductive proteins with those in digestive protein genes and well-conserved reference genes. We found that female reproductive proteins showed an overall ratio of nonsynonymous to synonymous substitutions (ω) similar to that of conserved genes, while seminal fluid proteins and digestive proteins exhibited higher overall ω values. Further, seminal fluid proteins and digestive proteins showed a higher proportion of sites putatively under positive selection, and explicit tests showed no difference in relaxed selection between protein types. Evolutionary rate covariation analyses showed that evolutionary rates among seminal fluid proteins were on average more closely correlated with those in female reproductive proteins than with either digestive or conserved genes. Gene expression showed the expected negative covariation with ω values, except for male-biased genes where this negative relationship was reversed. In conclusion, seminal fluid proteins showed relatively rapid evolution and signs of positive selection. In contrast, female reproductive proteins evolved at a lower rate under selective constraints, on par with genes known to be well conserved. Although our findings provide support for concerted evolution of seminal fluid proteins and female reproductive proteins, they also suggest that these two classes of proteins evolve under partly distinct selective regimes.

Molecular evolution of the mammalian kinetochore complex

Mammalian centromeres are satellite-rich chromatin domains that serve as sites for kinetochore complex assembly. Centromeres are highly variable in sequence and satellite organization across species, but the processes that govern the co-evolutionary dynamics between rapidly evolving centromeres and their associated kinetochore proteins remain poorly understood. Here, we pursue a course of phylogenetic analyses to investigate the molecular evolution of the complete kinetochore complex across primate and rodent species with divergent centromere repeat sequences and features. We show that many protein components of the core centromere associated network (CCAN) harbor signals of adaptive evolution, consistent with their intimate association with centromere satellite DNA and roles in the stability and recruitment of additional kinetochore proteins. Surprisingly, CCAN and outer kinetochore proteins exhibit comparable rates of adaptive divergence, suggesting that changes in centromere DNA can ripple across the kinetochore to drive adaptive protein evolution within distant domains of the complex. Our work further identifies kinetochore proteins subject to lineage-specific adaptive evolution, including rapidly evolving proteins in species with centromere satellites characterized by higher-order repeat structure and lacking CENP-B boxes. Thus, features of centromeric chromatin beyond the linear DNA sequence may drive selection on kinetochore proteins. Overall, our work spotlights adaptively evolving proteins with diverse centromere-associated functions, including centromere chromatin structure, kinetochore protein assembly, kinetochore-microtubule association, cohesion maintenance, and DNA damage response pathways. These adaptively evolving kinetochore protein candidates present compelling opportunities for future functional investigations exploring how their concerted changes with centromere DNA ensure the maintenance of genome stability.

Comparative genomics of the closely related fungal genera Cryptococcus and Kwoniella reveals karyotype dynamics and suggests evolutionary mechanisms of pathogenesis

In exploring the evolutionary trajectories of both pathogenesis and karyotype dynamics in fungi, we conducted a large-scale comparative genomic analysis spanning the Cryptococcus genus, encompassing both global human fungal pathogens and nonpathogenic species, and related species from the sister genus Kwoniella. Chromosome-level genome assemblies were generated for multiple species, covering virtually all known diversity within these genera. Although Cryptococcus and Kwoniella have comparable genome sizes (about 19.2 and 22.9 Mb) and similar gene content, hinting at preadaptive pathogenic potential, our analysis found evidence of gene gain (via horizontal gene transfer) and gene loss in pathogenic Cryptococcus species, which might represent evolutionary signatures of pathogenic development. Genome analysis also revealed a significant variation in chromosome number and structure between the 2 genera. By combining synteny analysis and experimental centromere validation, we found that most Cryptococcus species have 14 chromosomes, whereas most Kwoniella species have fewer (11, 8, 5, or even as few as 3). Reduced chromosome number in Kwoniella is associated with formation of giant chromosomes (up to 18 Mb) through repeated chromosome fusion events, each marked by a pericentric inversion and centromere loss. While similar chromosome inversion-fusion patterns were observed in all Kwoniella species with fewer than 14 chromosomes, no such pattern was detected in Cryptococcus. Instead, Cryptococcus species with less than 14 chromosomes showed reductions primarily through rearrangements associated with the loss of repeat-rich centromeres. Additionally, Cryptococcus genomes exhibited frequent interchromosomal translocations, including intercentromeric recombination facilitated by transposons shared between centromeres. Overall, our findings advance our understanding of genetic changes possibly associated with pathogenicity in Cryptococcus and provide a foundation to elucidate mechanisms of centromere loss and chromosome fusion driving distinct karyotypes in closely related fungal species, including prominent global human pathogens.

Study on the detection rate, genetic polymorphism, viral load, persistent infection capacity, and pathogenicity of human papillomavirus type 81

Human papillomavirus (HPV) type 81 has recently become one of the most common low-risk HPV types; however, literature focusing on it is limited. This study aimed to analyze the reasons for the increased detection rate of HPV81 and investigate its evolving pathogenicity. We analyzed the detection rates and trends of HPV81 in 229 061 exfoliated cervical cell samples collected from 2014 to 2023; collected samples of HPV81 single infections from two different time periods; and analyzed the allele frequencies, positive selection, viral load, persistent infection capacity, and pathogenicity of E6 and E7 genotypes. We found that the detection rate of HPV81 ranked first among the low-risk types in exfoliated cervical cells and exhibited a significantly increasing trend (p < 0.001). The frequency of the E6 prototype allele of HPV81 (n = 317) was significantly increased (p = 0.018) and demonstrated the strongest adaptive capacity. The viral load and persistent infection capacity of the E6 prototype were significantly higher than those of the mutants, thus serving as key drivers for increasing the detection rate of HPV81 and enhancing its pathogenicity. The viral load was positively correlated with persistent infection capacity and pathogenicity. Persistent infection was a crucial factor in the pathogenicity of HPV81. Successful adaptive evolution of HPV81 is accompanied by enhanced pathogenicity.