Molecular biology and evolution最新文献

Even though gene duplication is a key source of new genes and evolutionary innovation, it is unclear how duplicates survive the period immediately following gene duplication, in which both copies are functionally redundant. In the absence of epigenetic silencing, the abundance of the gene product would double after gene duplication, which would often have deleterious effects. However, recent duplicates exhibit low expression levels, which could be at least partially explained by high levels of promoter methylation. What evolutionary paths lead to duplicate hypermethylation, and does it affect both duplicates, or only one? Here, we compare levels of promoter methylation in 10 human and 16 mouse tissues, between singletons and duplicates and among human-mouse orthologs of different kinds (one-to-one, one-to-many, many-to-one, and many-to-many). Our results indicate that: (1) on average, duplicates are more methylated than singletons in mouse, but less methylated than singletons in human, (2) recently duplicated genes tend to exhibit high levels of promoter methylation, (3) genes that undergo duplication tend to be highly methylated before duplication, (4) after gene duplication, one of the copies (the daughter copy, i.e. the one that relocates to a new genomic context) tends to undergo an additional increase in promoter methylation, whereas the other (the parental copy, which remains in the original genomic location) tends to retain pre-duplication methylation levels, and (5) daughter copies tend to be lowly expressed. These observations support a model in which daughter copies are repressed via promoter hypermethylation and can thus survive the filter of purifying selection until both copies diverge functionally.

Understanding metabolic plasticity of animal evolution is a fundamental challenge in evolutionary biology. Owing to the diversification of insect wing morphology and dynamic energy requirements, the molecular adaptation mechanisms underlying the metabolic pathways in wing evolution remain largely unknown. This study reveals the pivotal role of the duplicated Apolipoprotein D (ApoD) gene in lipid and energy homeostasis in the lepidopteran wing. ApoD underwent significant expansion in insects, with gene duplication and consistent retention observed in Lepidoptera. Notably, duplicated ApoD2 was highly expressed in lepidopteran wings and encoded a unique C-terminal tail, conferring distinct ligand-binding properties. Using Bombyx mori as a model organism, we integrated evolutionary analysis, multiomics, and in vivo functional experiments to elucidate the way duplicated ApoD2 mediates lipid trafficking and homeostasis via the AMP-activated protein kinase pathway in wings. Moreover, we revealed the specific expression and functional divergence of duplicated ApoD as a key mechanism regulating lipid homeostasis in the lepidopteran wing. These findings highlight an evolutionary scenario in which neofunctionalization conferred a novel role of ApoD in shaping adaptive lipid metabolic regulatory networks during wing phenotypic evolution. Overall, we provide in vivo evidence for the functional differentiation of duplicate genes in shaping adaptive metabolic regulatory networks during phenotypic evolution.

Caspase-1 and related proteases are key players in inflammation and innate immunity. Here, we characterize the evolutionary history of caspase-1 and its close relatives across 19 primates and 21 rodents, focusing on differences that may cause discrepancies between humans and animal studies. While caspase-1 has been retained in all these taxa, other members of the caspase-1 subfamily (caspase-4, caspase-5, caspase-11, and caspase-12 and CARD16, 17, and 18) each have unique evolutionary trajectories. Caspase-4 is found across simian primates, whereas we identified multiple pseudogenization and gene loss events in caspase-5, caspase-11, and the CARDs. Because caspase-4 and caspase-11 are both key players in the noncanonical inflammasome pathway, we expected that these proteins would be likely to evolve rapidly. Instead, we found that these two proteins are largely conserved, whereas caspase-4's close paralog, caspase-5, showed significant indications of positive selection, as did primate caspase-1. Caspase-12 is a nonfunctional pseudogene in humans. We find this extends across most primates, although many rodents and some primates retain an intact, and likely functional, caspase-12. In mouse laboratory lines, we found that 50% of common strains carry nonsynonymous variants that may impact the functions of caspase-11 and caspase-12 and therefore recommend specific strains to be used (and avoided). Finally, unlike rodents, primate caspases have undergone repeated rounds of gene conversion, duplication, and loss leading to a highly dynamic proinflammatory caspase repertoire. Thus, we uncovered many differences in the evolution of primate and rodent proinflammatory caspases and discuss the potential implications of this history for caspase gene functions.

While whole-genome sequencing has been applied extensively to investigate the genetic diversity of global populations, ethnic minority groups in Pakistan are generally underrepresented. In particular, little is known about the genetic origin and highland adaptation of the Pamirian Wakhi people. According to Chinese historical records, the geographical location and language usage of Wakhi may be closely related to Xinjiang Tajiks. In this study, based on high-coverage (∼30×) whole-genome sequencing of eight Wakhi and 25 Xinjiang Tajik individuals, we performed data analyses together with worldwide populations to gain insights into their genetic composition, demography, and adaptive evolution to the highland environment. The Wakhi derived more than 85% of their ancestry from West Eurasian populations (European ∼44.5%, South Asian ∼42.2%) and 10% from East Eurasian populations (Siberian ∼6.0%, East Asian ∼4.3%). Modeling the admixture history of the Wakhi indicated that the early West-East admixture occurred ∼3,875 to 2,250 years ago and that the recent admixture occurred ∼750 to 375 years ago. We identified selection signatures across EGLN3, in particular, a distinctive evolutionary signature was observed, and a certain underlying selected haplotype showed higher frequency (87.5%) in the Wakhi than in nearby Xinjiang Tajiks and other highlanders. Interestingly, we found high-frequency archaic sequences in the Wakhi genome, which overlapped with several genes related to cellular signaling transduction, including MAGI2, previously associated with high-altitude adaptation. Our analysis indicates that the Wakhi are distinct from the Xinjiang Tajiks and Tajikistan Tajiks and sheds light on the Wakhi's ancestral origin and genetic basis of high-altitude adaptation.

Bats possess a range of distinctive characteristics, including flight, echolocation, impressive longevity, and the ability to harbor various zoonotic pathogens. Additionally, they account for the second-highest species diversity among mammalian orders, yet their phylogenetic relationships and demographic history remain underexplored. Here, we generated de novo assembled genomes for 17 bat species and 2 of their mammalian relatives (the Amur hedgehog and Chinese mole shrew), with 12 genomes reaching chromosome-level assembly. Comparative genomics and ChIP-seq assays identified newly gained genomic regions in bats potentially linked to the regulation of gene activity and expression. Notably, some antiviral infection-related gene under positive selection exhibited the activity of suppressing cancer, evidencing the linkage between virus tolerance and cancer resistance in bats. By integrating published bat genome assemblies, phylogenetic reconstruction established the proximity of noctilionoid bats to vesper bats. Interestingly, we found 2 distinct patterns of ancient population dynamics in bats and population changes since the last glacial maximum does not reflect species phylogenetic relationships. These findings enriched our understanding of adaptive mechanisms and demographic history of bats.

Pangenomes can facilitate a deeper understanding of genome complexity. Using de novo phased long-read assemblies of eight representative goat breeds, we constructed a graph-based pangenome of goats (Capra hircus) and discovered 113-Mb autosomal novel sequences. Combining this multi-assembly pangenome with low-coverage PacBio HiFi sequences, we constructed a long-read structural variations (SVs) database containing 59,325 SV deletions, 84,910 SV insertions, and 24,954 other complex SV alleles. This resource allowed reliable graph-based genotyping from short reads of 79 wild and 1,148 worldwide domestic goats. Selection signal analysis of SV captured a novel immune-related domestication locus containing the galectin-9 gene and extra copies of the ruminant-specific galectin-9-like genes (LGALS9L), which have high tissue specificity. A segmental duplication in domestic goats generates three additional LGALS9L copies. Ancient goat genome sequences show a gradual increase in frequency of this duplication from the Neolithic to the present. Two other newly detected SVs also have higher selection signals than adjacent SNPs, a truncated-LINE1 deletion in EDAR2 associated with cashmere production and a VNTR-related insertion in PAPSS2 linked to high-altitude adaptation. In summary, the multi-assembly goat pangenome and long-read SV database facilitates detecting complex variations that are important in evolution and selection.

Rates of molecular evolution vary greatly among even closely related species. Although theory predicts that antagonistic interactions between species increase rates of molecular evolution, predictions for how mutualism affects evolutionary rates are mixed. We compared rates of molecular evolution between (i) mutualistic and non-mutualistic legumes, (ii) an independent set of symbiotic rhizobia and their non-symbiotic close relatives, and (iii) symbiotic and non-symbiotic clades within Ensifer, a diverse genus of bacteria with various lifestyles. We assembled transcriptomes de novo for 12 legume species and calculated dN/dS ratios at orthologous genes in all species to determine if genes in mutualistic plants evolve faster or slower than in their non-mutualistic relatives. We also calculated dN/dS ratios in genes known to be important for symbiosis. We found that mutualists have higher rates of molecular evolution genome-wide compared to non-mutualistic legumes, but this pattern did not hold in symbiosis genes. We next calculated dN/dS ratios in 14 bacteria species across the proteobacteria phylogeny that differ in whether they associate mutualistically with plants, using published data. In most pairs, symbiotic rhizobia show higher dN/dS values compared to their non-symbiotic relatives. Within a bacterial genus with many well-characterized mutualist species (Ensifer), we calculated dN/dS ratios in symbiotic and non-symbiotic clades and found that symbiotic lineages have higher rates of molecular evolution genome-wide, but not at genes on the symbiotic plasmid pSymB. Our results suggest that although mutualism between legumes and rhizobia is associated with elevated rates of molecular evolution genome-wide, symbiosis genes may be evolutionarily stagnant.

Phylogenomics reveals reticulate evolution to be widespread across taxa, but whether reticulation is due to low statistical power or it is a true evolutionary pattern remains a field of study. Here, we investigate the phylogeny and quantify reticulation in the Drosophila saltans species group, a Neotropical clade of the subgenus Sophophora comprising 23 species whose relationships have long been problematic. Phylogenetic analyses revealed conflicting topologies between the X chromosome, autosomes and the mitochondria. We extended the ABBA-BABA test of asymmetry in phylogenetic discordance to cases where no "true" species tree could be inferred, and applied our new test (called 2A2B) to whole genome data and to individual loci. We used four strategies, two based on our new assemblies using either conserved genes or ≥50 kb-long syntenic blocks with conserved collinearity across Neotropical Sophophora, and two consisted of windows from pseudo-reference genomes aligned to either an ingroup or outgroup species. Evidence for reticulation varied among the strategies, being lowest in the synteny-based approach, where it did not exceed ∼7% of the blocks in the most conflicting species quartets. High incidences of reticulation were restricted to three nodes on the tree that coincided with major paleogeographical events in South America. Our results identify possible technical biases in quantifying reticulate evolution and indicate that episodic rapid radiations have played a major role in the evolution of a largely understudied Neotropical clade.

Malate dehydrogenases (MalDHs) (EC.1.1.1.37), which are involved in the conversion of oxaloacetate to pyruvate in the tricarboxylic acid cycle, are a relevant model for the study of enzyme evolution and adaptation. Likewise, a recent study showed that Methanococcales, a major lineage of Archaea, is a good model to study the molecular processes of proteome thermoadaptation in prokaryotes. Here, we use ancestral sequence reconstruction and paleoenzymology to characterize both ancient and extant MalDHs. We observe a good correlation between inferred optimal growth temperatures and experimental optimal temperatures for activity (A-Topt). In particular, we show that the MalDH present in the ancestor of Methanococcales was hyperthermostable and had an A-Topt of 80 °C, consistent with a hyperthermophilic lifestyle. This ancestor gave rise to two lineages with different thermal constraints: one remained hyperthermophilic, while the other underwent several independent adaptations to colder environments. Surprisingly, the enzymes of the first lineage have retained a thermoresistant behavior (i.e. strong thermostability and high A-Topt), whereas the ancestor of the second lineage shows a strong thermostability, but a reduced A-Topt. Using mutants, we mimic the adaptation trajectory toward mesophily and show that it is possible to significantly reduce the A-Topt without altering the thermostability of the enzyme by introducing a few mutations. Finally, we reveal an unexpected link between thermostability and the ability to resist γ-irradiation-induced unfolding.

In recent years, advances in image processing and machine learning have fueled a paradigm shift in detecting genomic regions under natural selection. Early machine learning techniques employed population-genetic summary statistics as features, which focus on specific genomic patterns expected by adaptive and neutral processes. Though such engineered features are important when training data are limited, the ease at which simulated data can now be generated has led to the recent development of approaches that take in image representations of haplotype alignments and automatically extract important features using convolutional neural networks. Digital image processing methods termed α-molecules are a class of techniques for multiscale representation of objects that can extract a diverse set of features from images. One such α-molecule method, termed wavelet decomposition, lends greater control over high-frequency components of images. Another α-molecule method, termed curvelet decomposition, is an extension of the wavelet concept that considers events occurring along curves within images. We show that application of these α-molecule techniques to extract features from image representations of haplotype alignments yield high true positive rate and accuracy to detect hard and soft selective sweep signatures from genomic data with both linear and nonlinear machine learning classifiers. Moreover, we find that such models are easy to visualize and interpret, with performance rivaling those of contemporary deep learning approaches for detecting sweeps.