Molecular biology and evolution最新文献

Although synonymous mutations are commonly assumed neutral or nearly so, recent years have seen reports of fitness effects of synonymous mutations detected under laboratory conditions. In a previous study, we used genome editing to construct thousands of yeast mutants each carrying a synonymous or nonsynonymous mutation in one of 21 genes, and discovered that most synonymous and most nonsynonymous mutations are deleterious. A concern was raised that this observation could be caused by the fitness effects of potential CRISPR/Cas9 off-target edits and/or secondary mutations, and an experiment that would be refractory to such effects was proposed. Using genome sequencing, we here show that no CRISPR/Cas9 off-target editing occurred, although some mutants did carry secondary mutations. Analysis of mutants with negligible effects from secondary mutations and new data collected from the proposed experiment confirms the original conclusion. These findings, along with other reports of fitness effects of synonymous mutations from both case and systematic studies, necessitate a paradigm shift from assuming (near) neutrality of synonymous mutations.

Functional innovation at the protein level is a key source of evolutionary novelties. The constraints on functional innovations are likely to be highly specific in different proteins, which are shaped by their unique histories and the extent of global epistasis that arises from their structures and biochemistries. These contextual nuances in the sequence-function relationship have implications both for a basic understanding of the evolutionary process and for engineering proteins with desirable properties. Here, we have investigated the molecular basis of novel function in a model member of an ancient, conserved, and biotechnologically relevant protein family. These Major Facilitator Superfamily sugar porters are a functionally diverse group of proteins that are thought to be highly plastic and evolvable. By dissecting a recent evolutionary innovation in an α-glucoside transporter from the yeast Saccharomyces eubayanus, we show that the ability to transport a novel substrate requires high-order interactions between many protein regions and numerous specific residues proximal to the transport channel. To reconcile the functional diversity of this family with the constrained evolution of this model protein, we generated new, state-of-the-art genome annotations for 332 Saccharomycotina yeast species spanning ∼400 My of evolution. By integrating phylogenetic and phenotypic analyses across these species, we show that the model yeast α-glucoside transporters likely evolved from a multifunctional ancestor and became subfunctionalized. The accumulation of additive and epistatic substitutions likely entrenched this subfunction, which made the simultaneous acquisition of multiple interacting substitutions the only reasonably accessible path to novelty.

Comparative genomics approaches seek to associate molecular evolution with the evolution of phenotypes across a phylogeny. Many of these methods lack the ability to analyze non-ordinal categorical traits with more than two categories. To address this limitation, we introduce an expansion to RERconverge that associates shifts in evolutionary rates with the convergent evolution of categorical traits. The categorical RERconverge expansion includes methods for performing categorical ancestral state reconstruction, statistical tests for associating relative evolutionary rates with categorical variables, and a new method for performing phylogeny-aware permutations, "permulations", on categorical traits. We demonstrate our new method on a three-category diet phenotype, and we compare its performance to binary RERconverge analyses and two existing methods for comparative genomic analyses of categorical traits: phylogenetic simulations and a phylogenetic signal based method. We present an analysis of how the categorical permulations scale with the number of species and the number of categories included in the analysis. Our results show that our new categorical method outperforms phylogenetic simulations at identifying genes and enriched pathways significantly associated with the diet phenotypes and that the categorical ancestral state reconstruction drives an improvement in our ability to capture diet-related enriched pathways compared to binary RERconverge when implemented without user input on phenotype evolution. The categorical expansion to RERconverge will provide a strong foundation for applying the comparative method to categorical traits on larger data sets with more species and more complex trait evolution than have previously been analyzed.

Comparative genomic studies of social insects suggest that changes in gene regulation are associated with evolutionary transitions in social behavior, but the activity of predicted regulatory regions has not been tested empirically. We used self-transcribing active regulatory region sequencing, a high-throughput enhancer discovery tool, to identify and measure the activity of enhancers in the socially variable sweat bee, Lasioglossum albipes. We identified over 36,000 enhancers in the L. albipes genome from 3 social and 3 solitary populations. Many enhancers were identified in only a subset of L. albipes populations, revealing rapid divergence in regulatory regions within this species. Population-specific enhancers were often proximal to the same genes across populations, suggesting compensatory gains and losses of regulatory regions may preserve gene activity. We also identified 1,182 enhancers with significant differences in activity between social and solitary populations, some of which are conserved regulatory regions across species of bees. These results indicate that social trait variation in L. albipes is associated with the fine-tuning of ancient enhancers as well as lineage-specific regulatory changes. Combining enhancer activity with population genetic data revealed variants associated with differences in enhancer activity and identified a subset of differential enhancers with signatures of selection associated with social behavior. Together, these results provide the first empirical map of enhancers in a socially flexible bee and highlight links between cis-regulatory variation and the evolution of social behavior.

The evolution of sex chromosomes can involve recombination suppression sometimes involving structural changes, such as inversions, allowing subsequent rearrangements, including inversions and gene transpositions. In the two major genus Salix clades, Salix and Vetrix, almost all species are dioecious, and sex-linked regions have evolved on chromosome 7 and 15, with either male or female heterogamety. We used chromosome conformation capture (Hi-C) and PacBio HiFi (high-fidelity) reads to assemble chromosome-level, gap-free X and Y chromosomes from both clades, S. triandra (15XY system), a basal species in the Vetrix clade, and the Salix clade species S. mesnyi (7XY system). Combining these with other available genome assemblies, we found inversions within the sex-linked regions, which are likely to be pericentromeric and probably recombined rarely in the ancestral species, before sex-linkage evolved. The Y-linked regions in all 15XY and 7XY species include partial duplicates containing exon 1 of an ARR17-like gene similar to male-determining factors in other Salicaceae species. We also found duplicates of a Y-specific gene, which we named MSF. The derived Salix clade 7XY chromosome systems appear to have evolved when these two genes transposed from the 15Y to the 7Y. Additionally, the 7Y chromosomes in S. dunnii and S. chaenomeloides probably evolved from the ancestral 7X of the Salix clade, involving a similar transposition, and loss of the ancestral 7Y. We suggest that pericentromeric regions that recombine infrequently may facilitate the evolution of sex linkage.

Adaptation to high-altitude hypoxia is characterized by systemic and organ-specific metabolic changes. This study investigates whether intestinal metabolic rewiring is a contributing factor to hypoxia adaptation. We conducted a longitudinal analysis over 108 days, with seven time points, examining fecal metabolomic data from a cohort of 46 healthy male adults traveling from Chongqing (a.s.l. 243 m) to Lhasa (a.s.l. 3,658 m) and back. Our findings reveal that short-term hypoxia exposure significantly alters intestinal metabolic pathways, particularly those involving purines, pyrimidines, and amino acids. A notable observation was the significantly reduced level of intestinal uric acid, the end product of purine metabolism, during acclimatization (also called acclimation) and additional two long-term exposed cohorts (Han Chinese and Tibetans) residing in Shigatse, Xizang (a.s.l. 4,700 m), suggesting that low intestinal uric acid levels facilitate adaptation to high-altitude hypoxia. Integrative analyses with gut metagenomic data showed consistent trends in intestinal uric acid levels and the abundance of key uric acid-degrading bacteria, predominantly from the Lachnospiraceae family. The sustained high abundance of these bacteria in the long-term resident cohorts underscores their essential role in maintaining low intestinal uric acid levels. Collectively, these findings suggest that the rewiring of intestinal uric acid metabolism, potentially orchestrated by gut bacteria, is crucial for enhancing human resilience and adaptability in extreme environments.

For a set of binary unrooted subtrees generating all binary unrooted trees compatible with them, i.e. generating their stand, is one of the classical problems in phylogenetics. Here, we introduce Gentrius-an efficient algorithm to tackle this task. The algorithm has a direct application in practice. Namely, Gentrius generates phylogenetic terraces-topologically distinct, equally scoring trees due to missing data. Despite stand generation being computationally intractable, we showed on simulated and biological datasets that Gentrius generates stands with millions of trees in feasible time. We exemplify that depending on the distribution of missing data across species and loci and the inferred phylogeny, the number of equally optimal terrace trees varies tremendously. The strict consensus tree computed from them displays all the branches unaffected by the pattern of missing data. Thus, by solving the problem of stand generation, in practice Gentrius provides an important systematic assessment of phylogenetic trees inferred from incomplete data. Furthermore, Gentrius can aid theoretical research by fostering understanding of tree space structure imposed by missing data.

The transcriptional regulatory network (TRN) in bacteria is thought to rapidly evolve in response to selection pressures, modulating transcription factor (TF) activities and interactions. In order to probe the limits and mechanisms surrounding the short-term adaptability of the TRN, we generated, evolved, and characterized knockout (KO) strains in Escherichia coli for 11 regulators selected based on measured growth impact on glucose minimal media. All but one knockout strain (Δlrp) were able to recover growth and did so requiring few convergent mutations. We found that the TF knockout adaptations could be divided into four categories: (i) Strains (ΔargR, ΔbasR, Δlon, ΔzntR, and Δzur) that recovered growth without any regulator-specific adaptations, likely due to minimal activity of the regulator on the growth condition, (ii) Strains (ΔcytR, ΔmlrA, and ΔybaO) that recovered growth without TF-specific mutations but with differential expression of regulators with overlapping regulons to the KO'ed TF, (iii) Strains (Δcrp and Δfur) that recovered growth using convergent mutations within their regulatory networks, including regulated promoters and connected regulators, and (iv) Strains (Δlrp) that were unable to fully recover growth, seemingly due to the broad connectivity of the TF within the TRN. Analyzing growth capabilities in evolved and unevolved strains indicated that growth adaptation can restore fitness to diverse substrates often despite a lack of TF-specific mutations. This work reveals the breadth of TRN adaptive mechanisms and suggests these mechanisms can be anticipated based on the network and functional context of the perturbed TFs.

Transposable elements (TEs) compose a substantial proportion of the largest eukaryotic genomes. TE diversity has been hypothesized to be negatively correlated with genome size, yet empirical demonstrations of such a relationship in a phylogenetic context are largely lacking. Furthermore, the most abundant type of TEs in genomes varies across groups, and it is not clear if there are patterns of TE activity consistent with genome size among different taxa with large genome sizes. We use low-coverage sequencing of 16 species of Neotropical salamanders, which vary ∼7-fold in genome size, to estimate TE relative abundance and diversity for each species. We also compare the divergence of copies of each TE superfamily to estimate patterns of TE activity in each species. We find a negative relationship between TE diversity and genome size, which is consistent with the hypothesis that either competition among TEs or reduced selection against ectopic recombination may result in lower diversity in the largest genomes. We also find divergent activity patterns in the largest versus the smallest genomes, suggesting that the history of TE activity may explain differences in genome size. Our results suggest that both TE diversity and relative abundance may be predictable, at least within taxonomic groups.

Mollusk-hunting (molluscivorous) cone snails belong to a monophyletic group in Conus, a genus of venomous marine snails. The molluscivorous lineage evolved from ancestral worm-hunting (vermivorous) snails ∼18 Ma. To enable the shift to a molluscivorous lifestyle, molluscivorous cone snails must solve biological problems encountered when hunting other gastropods, namely: (i) preventing prey escape and (ii) overcoming the formidable defense of the prey in the form of the molluscan shell, a problem unique to molluscivorous Conus. Here, we show that χ-conotoxins, peptides exclusively expressed in the venoms of molluscivorous Conus, provide solutions to the above problems. Injecting χ-conotoxins into the gastropod mollusk Aplysia californica results in impaired locomotion and uncoordinated hyperactivity. Impaired locomotion impedes escape, and a hyperactive snail will likely emerge from its shell, negating the protection the shell provides. Thus, χ-conotoxins are an evolutionary innovation that accompanied the emergence of molluscivory in Conus and provide solutions to problems posed by hunting other snails.