Search

The role of balancing selection is a long-standing evolutionary puzzle. Balancing selection is a crucial evolutionary process that maintains genetic variation (polymorphism) over extended periods of time; however, detecting it poses a significant challenge. Building upon the Polymorphism-aware phylogenetic Models (PoMos) framework rooted in the Moran model, we introduce a PoMoBalance model. This novel approach is designed to disentangle the interplay of mutation, genetic drift, and directional selection (GC-biased gene conversion), along with the previously unexplored balancing selection pressures on ultra-long timescales comparable with species divergence times by analyzing multi-individual genomic and phylogenetic divergence data. Implemented in the open-source RevBayes Bayesian framework, PoMoBalance offers a versatile tool for inferring phylogenetic trees as well as quantifying various selective pressures. The novel aspect of our approach in studying balancing selection lies in polymorphism-aware phylogenetic models' ability to account for ancestral polymorphisms and incorporate parameters that measure frequency-dependent selection, allowing us to determine the strength of the effect and exact frequencies under selection. We implemented validation tests and assessed the model on the data simulated with SLiM and a custom Moran model simulator. Real sequence analysis of Drosophila populations reveals insights into the evolutionary dynamics of regions subject to frequency-dependent balancing selection, particularly in the context of sex-limited color dimorphism in Drosophila erecta. [ABSTRACT FROM AUTHOR]

Different frequencies amongst codons that encode the same amino acid (i.e. synonymous codons) have been observed in multiple species. Studies focused on uncovering the forces that drive such codon usage showed that a combined effect of mutational biases and translational selection works to produce different frequencies of synonymous codons. However, only few have been able to measure and distinguish between these forces that may leave similar traces on the coding regions. Here, we have developed a codon model that allows the disentangling of mutation, selection on amino acids and synonymous codons, and GC-biased gene conversion (gBGC) which we employed on an extensive dataset of 415 chordates and 191 arthropods. We found that chordates need 15 more synonymous codon categories than arthropods to explain the empirical codon frequencies, which suggests that the extent of codon usage can vary greatly between animal phyla. Moreover, methylation at CpG sites seems to partially explain these patterns of codon usage in chordates but not in arthropods. Despite the differences between the two phyla, our findings demonstrate that in both, GC-rich codons are disfavored when mutations are GC-biased, and the opposite is true when mutations are AT-biased. This indicates that selection on the genomic coding regions might act primarily to stabilize its GC/AT content on a genome-wide level. Our study shows that the degree of synonymous codon usage varies considerably among animals, but is likely governed by a common underlying dynamic. [ABSTRACT FROM AUTHOR]

The European Union has committed to achieving carbon neutrality by 2050 and green hydrogen has been chosen as a priority vector for reaching that goal. Accordingly, Portugal has drafted a National Hydrogen Strategy laying out the various steps for the development of a green hydrogen economy. One element of this strategy is the development of a gigawatt-scale hydrogen production facility powered by dedicated renewable electricity sources. This work presents an analysis of the technical and economic feasibility of a facility consisting of a gigawatt-scale polymer electrolyte membrane electrolyser powered by solar photovoltaic and wind electricity, using the energy analysis model EnergyPLAN. Different capacities and modes of operation of the electrolyser are considered, including the complementary use of grid electricity as well as different combinations of renewable power, resulting in a total of 72 different configurations. An economic analysis is conducted addressing the related annualised capital expenditures, maintenance, and variable costs, to allow for the determination of the levelised cost of hydrogen for the different configurations. This analysis shows the conditions required for maximising annual hydrogen production at the lowest levelised cost of hydrogen. The best options consist of an electrolyser powered by a combination of solar photovoltaic and wind, with limited exchanges with the electricity grid, and a levelised cost of hydrogen in the range 3.13–3.48 EUR/kg. [ABSTRACT FROM AUTHOR]

For over a decade, experimental evolution has been combined with high‐throughput sequencing techniques. In so‐called Evolve‐and‐Resequence (E&R) experiments, populations are kept in the laboratory under controlled experimental conditions where their genomes are sampled and allele frequencies monitored. However, identifying signatures of adaptation in E&R datasets is far from trivial, and it is still necessary to develop more efficient and statistically sound methods for detecting selection in genome‐wide data. Here, we present Bait‐ER – a fully Bayesian approach based on the Moran model of allele evolution to estimate selection coefficients from E&R experiments. The model has overlapping generations, a feature that describes several experimental designs found in the literature. We tested our method under several different demographic and experimental conditions to assess its accuracy and precision, and it performs well in most scenarios. Nevertheless, some care must be taken when analysing trajectories where drift largely dominates and starting frequencies are low. We compare our method with other available software and report that ours has generally high accuracy even for trajectories whose complexity goes beyond a classical sweep model. Furthermore, our approach avoids the computational burden of simulating an empirical null distribution, outperforming available software in terms of computational time and facilitating its use on genome‐wide data. We implemented and released our method in a new open‐source software package that can be accessed at https://doi.org/10.5281/zenodo.7351736. [ABSTRACT FROM AUTHOR]

The availability of population genomic data through new sequencing technologies gives unprecedented opportunities for estimating important evolutionary forces such as genetic drift, selection and mutation biases across organisms. Yet, analytical methods that can handle polymorphisms jointly with sequence divergence across species are rare and not easily accessible to empiricists.We implemented polymorphism‐aware phylogenetic models (PoMos), an alternative approach for species tree estimation, in the Bayesian phylogenetic software RevBayes. PoMos naturally account for incomplete lineage sorting, which is known to cause difficulties for phylogenetic inference in species radiations, and scale well with genome‐wide data. Simultaneously, PoMos can estimate mutation and selection biases.We have applied our methods to resolve the complex phylogenetic relationships of a young radiation of Chorthippus grasshoppers, based on coding sequences. In addition to establishing a well‐supported species tree, we found a mutation bias favouring AT alleles and selection bias promoting the fixation of GC alleles, the latter consistent with GC‐biased gene conversion. The selection bias is two orders of magnitude lower than genetic drift, validating the critical role of nearly neutral evolutionary processes in species radiation.PoMos offer a wide range of models to reconstruct phylogenies and can be easily combined with existing models in RevBayes—for example, relaxed clock and divergence time estimation—offering new insights into the evolutionary processes underlying molecular evolution and, ultimately, species diversification. [ABSTRACT FROM AUTHOR]

Despite the importance of natural selection in species' evolutionary history, phylogenetic methods that take into account population-level processes typically ignore selection. The assumption of neutrality is often based on the idea that selection occurs at a minority of loci in the genome and is unlikely to compromise phylogenetic inferences significantly. However, genome-wide processes like GC-bias and some variation segregating at the coding regions are known to evolve in the nearly neutral range. As we are now using genome-wide data to estimate species trees, it is natural to ask whether weak but pervasive selection is likely to blur species tree inferences. We developed a polymorphism-aware phylogenetic model tailored for measuring signatures of nucleotide usage biases to test the impact of selection in the species tree. Our analyses indicate that although the inferred relationships among species are not significantly compromised, the genetic distances are systematically underestimated in a node-height-dependent manner: that is, the deeper nodes tend to be more underestimated than the shallow ones. Such biases have implications for molecular dating. We dated the evolutionary history of 30 worldwide fruit fly populations, and we found signatures of GC-bias considerably affecting the estimated divergence times (up to 23%) in the neutral model. Our findings call for the need to account for selection when quantifying divergence or dating species evolution. [ABSTRACT FROM AUTHOR]

The SARS-CoV-2 pandemic has led to unprecedented, nearly real-time genetic tracing due to the rapid community sequencing response. Researchers immediately leveraged these data to infer the evolutionary relationships among viral samples and to study key biological questions, including whether host viral genome editing and recombination are features of SARS-CoV-2 evolution. This global sequencing effort is inherently decentralized and must rely on data collected by many labs using a wide variety of molecular and bioinformatic techniques. There is thus a strong possibility that systematic errors associated with lab—or protocol—specific practices affect some sequences in the repositories. We find that some recurrent mutations in reported SARS-CoV-2 genome sequences have been observed predominantly or exclusively by single labs, co-localize with commonly used primer binding sites and are more likely to affect the protein-coding sequences than other similarly recurrent mutations. We show that their inclusion can affect phylogenetic inference on scales relevant to local lineage tracing, and make it appear as though there has been an excess of recurrent mutation or recombination among viral lineages. We suggest how samples can be screened and problematic variants removed, and we plan to regularly inform the scientific community with our updated results as more SARS-CoV-2 genome sequences are shared (https://virological.org/t/issues-with-sars-cov-2-sequencing-data/473 and https://virological.org/t/masking-strategies-for-sars-cov-2-alignments/480). We also develop tools for comparing and visualizing differences among very large phylogenies and we show that consistent clade- and tree-based comparisons can be made between phylogenies produced by different groups. These will facilitate evolutionary inferences and comparisons among phylogenies produced for a wide array of purposes. Building on the SARS-CoV-2 Genome Browser at UCSC, we present a toolkit to compare, analyze and combine SARS-CoV-2 phylogenies, find and remove potential sequencing errors and establish a widely shared, stable clade structure for a more accurate scientific inference and discourse. Author summary: SARS-CoV-2 genome sequences have been produced by hundreds of labs across the world. Idiosyncratic data generation or processing has the potential to inject non-random errors into genome sequences provided by individual lab groups. Here we show that these sites can be detected and removed by identifying variants that appear to reoccur many times across a phylogeny and are associated with specific lab groups. We show that the identified variants are consistent with expectations for recurrent error. These sites may produce spurious signatures of natural selection and viral recombination. We also provide methods for comparisons and visualization of extremely large phylogenies. [ABSTRACT FROM AUTHOR]

Typical avian eyes are phenotypically engineered for photopic vision (daylight). In contrast, the highly derived eyes of the barn owl (Tyto alba) are adapted for scotopic vision (dim light). The dramatic modifications distinguishing barn owl eyes from other birds include: 1) shifts in frontal orientation to improve binocularity, 2) rod-dominated retina, and 3) enlarged corneas and lenses. Some of thesefeaturesparallelmammalianeyepatterns,whicharehypothesizedtohaveinitiallyevolvedinnocturnalenvironments.Here,we used an integrative approach combining phylogenomics and functional phenotypes of 211 eye-development genes across 48 avian genomesrepresentingmostavianorders,includingthestemlineageofthescotopic-adaptedbarnowl.Overall,weidentified25eyedevelopment genes that coevolved under intensified or relaxed selection in the retina, lens, cornea, and optic nerves of the barn owl. The agtpbp1 gene, which is associated with the survival of photoreceptor populations, was pseudogenized in the barn owl genome. Our results further revealed that barn owl retinal genes responsible for the maintenance, proliferation, and differentiation of photoreceptors experienced an evolutionary relaxation. Signatures of relaxed selection were also observed in the lens and cornea morphology-associated genes, suggesting that adaptive evolution in these structures was essentially structural. Four eyedevelopment genes (ephb1, phactr4, prph2, and rs1) evolved in positive association with the orbit convergence in birds and under relaxedselectioninthebarnowllineage,likelycontributingtoanincreasedrelianceonbinocularvisioninthebarnowl.Moreover,we found evidence of coevolutionary interactions among genes that are expressed in the retina, lens, and optic nerve, suggesting synergetic adaptive events. Our study disentangles the genomic changes governing the binocularity and low-light perception adaptations of barn owls to nocturnal environments while revealing the molecular mechanisms contributing to the shift from the typical avian photopic vision to the more-novel scotopic-adapted eye. [ABSTRACT FROM AUTHOR]

As multi-individual population-scale data become available, more complex modeling strategies are needed to quantify genome-wide patterns of nucleotide usage and associated mechanisms of evolution. Recently, the multivariate neutral Moran model was proposed. However, it was shown insufficient to explain the distribution of alleles in great apes. Here, we propose a new model that includes allelic selection. Our theoretical results constitute the basis of a new Bayesian framework to estimate mutation rates and selection coefficients from population data. We apply the new framework to a great ape dataset, where we found patterns of allelic selection that match those of genome-wide GC-biased gene conversion (gBGC). In particular, we show that great apes have patterns of allelic selection that vary in intensity--a feature that we correlated with great apes' distinct demographies. We also demonstrate that the AT/GC toggling effect decreases the probability of a substitution, promoting more polymorphisms in the base composition of great ape genomes. We further assess the impact of GC-bias in molecular analysis, and find that mutation rates and genetic distances are estimated under bias when gBGC is not properly accounted for. Our results contribute to the discussion on the tempo and mode of gBGC evolution, while stressing the need for gBGC-aware models in population genetics and phylogenetics. [ABSTRACT FROM AUTHOR]

Motivation Determining whether a trait and phylogeny share some degree of phylogenetic signal is a flagship goal in evolutionary biology. Signatures of phylogenetic signal can assist the resolution of a broad range of evolutionary questions regarding the tempo and mode of phenotypic evolution. However, despite the considerable number of strategies to measure it, few and limited approaches exist for categorical traits. Here, we used the concept of Shannon entropy and propose the δ statistic for evaluating the degree of phylogenetic signal between a phylogeny and categorical traits. Results We validated δ as a measure of phylogenetic signal: the higher the δ -value the higher the degree of phylogenetic signal between a given tree and a trait. Based on simulated data we proposed a threshold-based classification test to pinpoint cases of phylogenetic signal. The assessment of the test's specificity and sensitivity suggested that the δ approach should only be applied to 20 or more species. We have further tested the performance of δ in scenarios of branch length and topology uncertainty, unbiased and biased trait evolution and trait saturation. Our results showed that δ may be applied in a wide range of phylogenetic contexts. Finally, we investigated our method in 14 360 mammalian gene trees and found that olfactory receptor genes are significantly associated with the mammalian activity patterns, a result that is congruent with expectations and experiments from the literature. Our application shows that δ can successfully detect molecular signatures of phenotypic evolution. We conclude that δ represents a useful measure of phylogenetic signal since many phenotypes can only be measured in categories. Availability and implementation https://github.com/mrborges23/delta%5fstatistic. Supplementary information Supplementary data are available at Bioinformatics online. [ABSTRACT FROM AUTHOR]

The Asian clam, Corbicula fluminea, is an invasive alien species (IAS) originally from Asia that has spread worldwide causing major ecological and economic impacts in aquatic ecosystems. Here, we evaluated C. fluminea genetic (using COI mtDNA, CYTb mtDNA and 18S rDNA gene markers), morphometric and sperm morphology variation in Portuguese freshwater ecosystems. The COI marker revealed a single haplotype, which belongs to the Asian FW5 invasive lineage, suggesting a common origin for all the 13 Portuguese C. fluminea populations analysed. Morphometric analyses showed differences between the populations colonizing the North (with the exception of the Lima River) and the Centre/South ecosystems. The sperm morphology examination revealed the presence of biflagellate sperm, a distinctive character of the invasive androgenetic lineages. The low genetic variability of the Portuguese C. fluminea populations and the pattern of sperm morphology have been illuminating for understanding the demographic history of this invasive species. We hypothesize that these populations were derived from a unique introductory event of a Corbicula fluminea FW5 invasive androgenic lineage in the Tejo River, which subsequently dispersed to other Portuguese freshwater ecosystems. The C. fluminea asexual reproductive mode may have assisted these populations to become highly invasive despite the low genetic diversity. [ABSTRACT FROM AUTHOR]

Background: The wide range of complex photic systems observed in birds exemplifies one of their key evolutionary adaptions, a well-developed visual system. However, genomic approaches have yet to be used to disentangle the evolutionary mechanisms that govern evolution of avian visual systems. Results: We performed comparative genomic analyses across 48 avian genomes that span extant bird phylogenetic diversity to assess evolutionary changes in the 17 representatives of the opsin gene family and five plumage coloration genes. Our analyses suggest modern birds have maintained a repertoire of up to 15 opsins. Synteny analyses indicate that PARA and PARIE pineal opsins were lost, probably in conjunction with the degeneration of the parietal organ. Eleven of the 15 avian opsins evolved in a non-neutral pattern, confirming the adaptive importance of vision in birds. Visual conopsins sw1, sw2 and lw evolved under negative selection, while the dim-light RH1 photopigment diversified. The evolutionary patterns of sw1 and of violet/ultraviolet sensitivity in birds suggest that avian ancestors had violet-sensitive vision. Additionally, we demonstrate an adaptive association between the RH2 opsin and the MC1R plumage color gene, suggesting that plumage coloration has been photic mediated. At the intra-avian level we observed some unique adaptive patterns. For example, barn owl showed early signs of pseudogenization in RH2, perhaps in response to nocturnal behavior, and penguins had amino acid deletions in RH2 sites responsible for the red shift and retinal binding. These patterns in the barn owl and penguins were convergent with adaptive strategies in nocturnal and aquatic mammals, respectively. Conclusions: We conclude that birds have evolved diverse opsin adaptations through gene loss, adaptive selection and coevolution with plumage coloration, and that differentiated selective patterns at the species level suggest novel photic pressures to influence evolutionary patterns of more-recent lineages. [ABSTRACT FROM AUTHOR]

Polymorphism-aware phylogenetic models (PoMo) constitute an alternative approach for species tree estimation from genome-wide data. PoMo builds on the standard substitution models of DNA evolution but expands the classic alphabet of the four nucleotide bases to include polymorphic states. By doing so, PoMo accounts for ancestral and current intra-population variation, while also accommodating population-level processes ruling the substitution process (e.g. genetic drift, mutations, allelic selection). PoMo has shown to be a valuable tool in several phylogenetic applications but a proof of statistical consistency (and identifiability, a necessary condition for consistency) is lacking. Here, we prove that PoMo is identifiable and, using this result, we further show that the maximum a posteriori (MAP) tree estimator of PoMo is a consistent estimator of the species tree. We complement our theoretical results with a simulated data set mimicking the diversity observed in natural populations exhibiting incomplete lineage sorting. We implemented PoMo in a Bayesian framework and show that the MAP tree easily recovers the true tree for typical numbers of sites that are sampled in genome-wide analyses. [ABSTRACT FROM AUTHOR]

Background: Based on evolutionary patterns of the vertebrate eye, Walls (1942) hypothesized that early placental mammals evolved primarily in nocturnal habitats. However, not only Eutheria, but all mammals show photic characteristics (i.e. dichromatic vision, rod-dominated retina) suggestive of a scotopic eye design. Results: Here, we used integrative comparative genomic and phylogenetic methodologies employing the photoreceptive opsin gene family in 154 mammals to test the likelihood of a nocturnal period in the emergence of all mammals. We showed that mammals possess genomic patterns concordant with a nocturnal ancestry. The loss of the RH2, VA, PARA, PARIE and OPN4x opsins in all mammals led us to advance a probable and most-parsimonious hypothesis of a global nocturnal bottleneck that explains the loss of these genes in the emerging lineage (> > 215.5 million years ago). In addition, ancestral character reconstruction analyses provided strong evidence that ancestral mammals possessed a nocturnal lifestyle, ultra-violet-sensitive vision, low visual acuity and low orbit convergence (i.e. panoramic vision). Conclusions: Overall, this study provides insight into the evolutionary history of the mammalian eye while discussing important ecological aspects of the photic paleo-environments ancestral mammals have occupied. [ABSTRACT FROM AUTHOR]

Summary: The phylogenetic signal, frequently used to identify signatures of adaptive evolution or important associations between genes and phenotypes, measures the tendency for recently diverged species to resemble each other more than distantly related species. An example of such a measure is the δ statistic, which uses Shannon entropy to measure the degree of phylogenetic signal between a categorical trait and a phylogeny. In this study, we refined this statistic to account for tree uncertainty, resulting in more accurate assessments of phylogenetic associations. In addition, we provided a more accessible and computationally efficient implementation of the δ statistic that will facilitate its use by the evolutionary community., Availability and Implementation: github.com/diogo-s-ribeiro/delta-statistic., (© The Author(s) 2023. Published by Oxford University Press.)