Polymorphism-Aware Species Trees with Advanced Mutation Models, Bootstrap, and Rate Heterogeneity

Pollinator shifts despite hybridisation in the Cape's hyperdiverse heathers (Erica, Ericaceae)

Interrogating the ecological and geographic factors that influence population divergence dynamics can reveal why some groups of organisms diversify more prolifically than others. One such group is the heathers (Erica, Ericaceae), the largest plant genus in the Cape Floristic Region. We study Erica abietina, a highly variable species complex with four subspecies differing in geographic range, habitat and pollination syndrome. We test for population differentiation, hybridisation, introgression and pollinator-driven divergence using genotyping-by-sequencing on samples across the entire distribution. We find five variably distinct genetic groups, with one subspecies comprising two independent lineages that are geographically isolated and occur on different soil types. Phylogenetic analysis suggests two independent shifts between bird and insect pollination, with accompanying genetic divergence. However, for one pair of populations with different pollinators, we uncover several individuals of hybrid origin at a site of sympatry. These results suggest that floral differentiation driven by divergent selection acts in concert with geographic isolation to maintain reproductive isolation and promote speciation. Our investigations reveal a highly dynamic system whose diversity has been shaped by a variety of interacting forces. We suggest that such a system could be a model for much of the diversification of the Cape flora.

Polymorphism-Aware Models in RevBayes: Species Trees, Disentangling Balancing Selection, and GC-Biased Gene Conversion

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.

Genetic characterization of the first Deltacoronavirus from wild birds around Qinghai Lake

Deltacoronavirus, widely distributed among pigs and wild birds, pose a significant risk of cross-species transmission, including potential human epidemics. Metagenomic analysis of bird samples from Qinghai Lake, China in 2021 reported the presence of Deltacoronavirus. A specific gene fragment of Deltacoronavirus was detected in fecal samples from wild birds at a positive rate of 5.94% (6/101). Next-generation sequencing (NGS) identified a novel Deltacoronavirus strain, which was closely related to isolates from the United Arab Emirates (2018), China (2022), and Poland (2023). Subsequently the strain was named A/black-headed gull/Qinghai/2021(BHG-QH-2021) upon confirmation of the Cytochrome b gene of black-headed gull in the sample. All available genome sequences of avian Deltacoronavirus, including the newly identified BHG-QH-2021 and 5 representative strains of porcine Deltacoronavirus (PDCoV), were classified according to ICTV criteria. In contrast to Coronavirus HKU15, which infects both mammals and birds and shows the possibility of cross-species transmission from bird to mammal host, our analysis revealed that BHG-QH-2021 is classified as Putative species 4. Putative species 4 has been reported to infect 5 species of birds but not mammals, suggesting that cross-species transmission of Putative species 4 is more prevalent among birds. Recombination analysis traced BHG-QH-2021 origin to dut148cor1 and MW01_1o strains, with MW01_1o contributing the S gene. Surprisingly, SwissModle prediction showed that the optimal template for receptor-binding domain (RBD) of BHG-QH-2021 is derived from the human coronavirus 229E, a member of the Alphacoronavirus, rather than the anticipated RBD structure of PDCoV of Deltacoronavirus. Further molecular docking analysis revealed that substituting the loop 1-2 segments of HCoV-229E significantly enhanced the binding capability of BHG-QH-2021 with human Aminopeptidase N (hAPN), surpassing its native receptor-binding domain (RBD). Most importantly, this finding was further confirmed by co-immunoprecipitation experiment that loop 1-2 segments of HCoV-229E enable BHG-QH-2021 RBD binding to hAPN, indicating that the loop 1-2 segment of the RBD in Putative species 4 is a probable key determinant for the virus ability to spill over into humans. Our results summarize the phylogenetic relationships among known Deltacoronavirus, reveal an independent putative avian Deltacoronavirus species with inter-continental and inter-species transmission potential, and underscore the importance of continuous surveillance of wildlife Deltacoronavirus.

ncOrtho: efficient and reliable identification of miRNA orthologs

MicroRNAs (miRNAs) are post-transcriptional regulators that finetune gene expression via translational repression or degradation of their target mRNAs. Despite their functional relevance, frameworks for the scalable and accurate detection of miRNA orthologs are missing. Consequently, there is still no comprehensive picture of how miRNAs and their associated regulatory networks have evolved. Here we present ncOrtho, a synteny informed pipeline for the targeted search of miRNA orthologs in unannotated genome sequences. ncOrtho matches miRNA annotations from multi-tissue transcriptomes in precision, while scaling to the analysis of hundreds of custom-selected species. The presence-absence pattern of orthologs to 266 human miRNA families across 402 vertebrate species reveals four bursts of miRNA acquisition, of which the most recent event occurred in the last common ancestor of higher primates. miRNA families are rarely modified or lost, but notable exceptions for both events exist. miRNA co-ortholog numbers faithfully indicate lineage-specific whole genome duplications, and miRNAs are powerful markers for phylogenomic analyses. Their exceptionally low genetic diversity makes them suitable to resolve clades where the phylogenetic signal is blurred by incomplete lineage sorting of ancestral alleles. In summary, ncOrtho allows to routinely consider miRNAs in evolutionary analyses that were thus far reserved to protein-coding genes.

Phylogenomic comparative methods: Accurate evolutionary inferences in the presence of gene tree discordance

Phylogenetic comparative methods have long been a mainstay of evolutionary biology, allowing for the study of trait evolution across species while accounting for their common ancestry. These analyses typically assume a single, bifurcating phylogenetic tree describing the shared history among species. However, modern phylogenomic analyses have shown that genomes are often composed of mosaic histories that can disagree both with the species tree and with each other-so-called discordant gene trees. These gene trees describe shared histories that are not captured by the species tree, and therefore that are unaccounted for in classic comparative approaches. The application of standard comparative methods to species histories containing discordance leads to incorrect inferences about the timing, direction, and rate of evolution. Here, we develop two approaches for incorporating gene tree histories into comparative methods: one that constructs an updated phylogenetic variance-covariance matrix from gene trees, and another that applies Felsenstein's pruning algorithm over a set of gene trees to calculate trait histories and likelihoods. Using simulation, we demonstrate that our approaches generate much more accurate estimates of tree-wide rates of trait evolution than standard methods. We apply our methods to two clades of the wild tomato genus Solanum with varying rates of discordance, demonstrating the contribution of gene tree discordance to variation in a set of floral traits. Our approaches have the potential to be applied to a broad range of classic inference problems in phylogenetics, including ancestral state reconstruction and the inference of lineage-specific rate shifts.

Range-wide population genomics of common seadragons shows secondary contact over a former barrier and insights on illegal capture

BACKGROUND

Common seadragons (Phyllopteryx taeniolatus, Syngnathidae) are an emblem of the diverse endemic fauna of Australia's southern rocky reefs, the newly recognized "Great Southern Reef." A lack of assessments spanning this global biodiversity hotspot in its entirety is currently hampering an understanding of the factors that have contributed to its diversity. The common seadragon has a wide range across Australia's entire temperate south and includes a geogenetic break over a former land bridge, which has called its status as a single species into question. As a popular aquarium display that sells for high prices, common seadragons are also vulnerable to illegal capture.

RESULTS

Here, we provide range-wide nuclear sequences (986 variable Ultraconserved Elements) for 198 individuals and mitochondrial genomes for 140 individuals to assess species status, identify genetic units and their diversity, and trace the source of two poached individuals. Using published data of the other two seadragon species, we found that lineages of common seadragons have diverged relatively recently (< 0.63 Ma). Within common seadragons, we found pronounced genetic structure, falling into three major groups in the western, central, and eastern parts of the range. While populations across the Bassian Isthmus were divergent, there is also evidence for secondary contact since the passage opened. We found a strong cline of genetic diversity from the range center tapering symmetrically towards the range peripheries. Based on their genetic similarities, the poached individuals were inferred to have originated from around Albany in southwestern Australia.

CONCLUSIONS

We conclude that common seadragons constitute a single species with strong geographic structure but coherence through gene flow. The low genetic diversity on the east and west coasts is concerning given that these areas are projected to face fast climate change. Our results suggest that in addition to their life history, geological events and demographic expansions have all played a role in shaping populations in the temperate south. These insights are an important step towards understanding the historical determinants of the diversity of species endemic to the Great Southern Reef.

The North Asian Genus Kolhymamnicola Starobogatov and Budnikova 1976 (Gastropoda: Amnicolidae), Its Extended Diagnosis, Distribution, and Taxonomic Relationships

The taxonomic position and phylogenetic affinities of the endemic North Asian genus Kolhymamnicola Starobogatov and Budnikova, 1976 (Gastropoda: Amnicolidae) remain unknown. To resolve this, we studied key morpho-anatomical characteristics of Kolhymamnicola snails and performed a molecular phylogenetic analysis based on sequences of COI mtDNA, 16S rRNA, and 18S rRNA genes. In terms of protoconch microsculpture, operculum, radular teeth, and gill complex morphology, Kolhymamnicola snails do not differ significantly from the North American genera Amnicola Gould and Haldeman, 1840 and Taylorconcha Hershler et al., 1994, and the European genus Marstoniopsis van Regteren Altena 1936. The bifid penis found in Kolhymamnicola is similar to that in the genus Marstoniopsis. The female reproductive anatomy has some features shared by Kolhymamnicola and Taylorconcha (absence of bursa copulatrix, single seminal receptacle in rs2′ position, and ventral channel). The molecular analysis has revealed Taylorconcha as the closest relative to Kolhymamnicola; the COI-based genetic distance between them amounted to 0.113. We discuss the possible time of divergence of these two genera, as well as of European Marstoniopsis and the Baikal Lake endemic family Baicaliidae. The last common ancestor of these groups was widely distributed in Miocene–Pliocene in the Holarctic waterbodies. Recent Kolhymamnicola snails are distributed in Northern Asia, including lakes of the Baikal rift zone. We rank the Baicaliidae as a family rather than a subfamily of Amnicolidae based on their distinct, unique morpho-anatomical characteristics and highly supported separate position on the molecular tree. The tribe Erhaiini Davis and Kuo, 1985 is elevated to the rank of the family, with 3–4 recent genera included. The family Palaeobaicaliidae Sitnikova et Vinarski fam. nov. is established to embrace the Cretaceous North Asian gastropods conchologically similar to the recent Baicaliidae and Pyrgulidae.

Bactrocera dorsalis in the Indian Ocean: A tale of two invasions

An increasing number of invasive fruit fly pests are colonizing new grounds. With this study, we aimed to uncover the invasion pathways of the oriental fruit fly, Bactrocera dorsalis into the islands of the Indian Ocean. By using genome-wide SNP data and a multipronged approach consisting of PCA, ancestry analysis, phylogenetic inference, and kinship networks, we were able to resolve two independent invasion pathways. A western invasion pathway involved the stepping-stone migration of B. dorsalis from the east African coast into the Comoros, along Mayotte and into Madagascar with a decreasing genetic diversity. The Mascarene islands (Reunion and Mauritius), on the contrary, were colonized directly from Asia and formed a distinct cluster. The low nucleotide diversity suggests that only a few genotypes invaded the Mascarenes. The presence of many long runs of homozygosity (ROH) in the introduced populations is indicative of population bottlenecks, with evidence of a more severe bottleneck for populations along the western migration pathway than on the Mascarene islands. More strict phytosanitary regulations are recommended in order to prevent the further spread of B. dorsalis.

Origin and diversification of the saguaro cactus (Carnegiea gigantea): a within-species phylogenomic analysis

Reconstructing accurate historical relationships within a species poses numerous challenges, not least in many plant groups in which gene flow is high enough to extend well beyond species boundaries. Nonetheless, the extent of tree-like history within a species is an empirical question on which it is now possible to bring large amounts of genome sequence to bear. We assess phylogenetic structure across the geographic range of the saguaro cactus, an emblematic member of Cactaceae, a clade known for extensive hybridization and porous species boundaries. Using 200 Gb of whole genome resequencing data from 20 individuals sampled from 10 localities, we assembled two data sets comprising 150,000 biallelic single nucleotide polymorphisms (SNPs) from protein coding sequences. From these we inferred within-species trees and evaluated their significance and robustness using five qualitatively different inference methods. Despite the low sequence diversity, large census population sizes, and presence of wide-ranging pollen and seed dispersal agents, phylogenetic trees were well resolved and highly consistent across both data sets and all methods. We inferred that the most likely root, based on marginal likelihood comparisons, is to the east and south of the region of highest genetic diversity, which lies along the coast of the Gulf of California in Sonora, Mexico. Together with striking decreases in marginal likelihood found to the north, this supports hypotheses that saguaro's current range reflects post-glacial expansion from the refugia in the south of its range. We conclude with observations about practical and theoretical issues raised by phylogenomic data sets within species, in which SNP-based methods must be used rather than gene tree methods that are widely used when sequence divergence is higher. These include computational scalability, inference of gene flow, and proper assessment of statistical support in the presence of linkage effects.

Rapid radiation in a highly diverse marine environment

Rapid diversification is often observed when founding species invade isolated or newly formed habitats that provide ecological opportunity for adaptive radiation. However, most of the Earth's diversity arose in diverse environments where ecological opportunities appear to be more constrained. Here, we present a striking example of a rapid radiation in a highly diverse marine habitat. The hamlets, a group of reef fishes from the wider Caribbean, have radiated into a stunning diversity of color patterns but show low divergence across other ecological axes. Although the hamlet lineage is ∼26 My old, the radiation appears to have occurred within the last 10,000 generations in a burst of diversification that ranks among the fastest in fishes. As such, the hamlets provide a compelling backdrop to uncover the genomic elements associated with phenotypic diversification and an excellent opportunity to build a broader comparative framework for understanding the drivers of adaptive radiation. The analysis of 170 genomes suggests that color pattern diversity is generated by different combinations of alleles at a few large-effect loci. Such a modular genomic architecture of diversification has been documented before in Heliconius butterflies, capuchino finches, and munia finches, three other tropical radiations that took place in highly diverse and complex environments. The hamlet radiation also occurred in a context of high effective population size, which is typical of marine populations. This allows for the accumulation of new variants through mutation and the retention of ancestral genetic variation, both of which appear to be important in this radiation.

Nucleotide Usage Biases Distort Inferences of the Species Tree

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.

Taxonomy notes on Vandeae (Orchidaceae) from China: Five new species and two new records

Five new species (Gastrochilus yei, Gastrochilus minimus, Luisia simaoensis, Taeniophyllum xizangense, Tuberolabium subulatum) and two newly recorded species (Cleisostoma tricornutum, Luisia inconspicua) of Vandeae (Orchidaceae) from China are described and illustrated. Gastrochilus yei is similar to G. affinis and G. nepalensis, but differs from them by having an epichile not lobed, the apex of the hypochile not bilobed, and a tine on the apex of the leaf. Gastrochilus minimus is similar to G. acinacifolius, but can be distinguished from the latter by having a flabellate epichile that is densely hirsute on the adaxial surface and an inconspicuous central cushion; in addition, the hypochile of G. minimus has a keel that extends to the apex of the epichile. Taeniophyllum xizangense is similar to T. stella and T. radiatum, but it is distinguished from them by having much bigger flowers, inflorescences densely covered with short-bristly hairs, papillae on the external surface of sepals, and bigger triangular-ovate viscidium. Luisia simaoensis is similar to L. magniflora and L. ramosii, but can be easily distinguished from them by having lateral sepals longer than dorsal sepals and petals, lip with irregular and waved margins, and lip with bilobed apex. Luisia inconspicua is moved from Gastrochilus to Luisia based on phylogenetic analyses of plastid matK sequence data. Tuberolabium subulatum is similar to T. carnosum, but it can be easily distinguished from the latter by having an inflorescence much shorter than the leaves, yellow sepals and petals, and many small papillae outside the lip lobes.

Genomic skimming and nanopore sequencing uncover cryptic hybridization in one of world's most threatened primates

The Brazilian buffy-tufted-ear marmoset (Callithrix aurita), one of the world's most endangered primates, is threatened by anthropogenic hybridization with exotic, invasive marmoset species. As there are few genetic data available for C. aurita, we developed a PCR-free protocol with minimal technical requirements to rapidly generate genomic data with genomic skimming and portable nanopore sequencing. With this direct DNA sequencing approach, we successfully determined the complete mitogenome of a marmoset that we initially identified as C. aurita. The obtained nanopore-assembled sequence was highly concordant with a Sanger sequenced version of the same mitogenome. Phylogenetic analyses unexpectedly revealed that our specimen was a cryptic hybrid, with a C. aurita phenotype and C. penicillata mitogenome lineage. We also used publicly available mitogenome data to determine diversity estimates for C. aurita and three other marmoset species. Mitogenomics holds great potential to address deficiencies in genomic data for endangered, non-model species such as C. aurita. However, we discuss why mitogenomic approaches should be used in conjunction with other data for marmoset species identification. Finally, we discuss the utility and implications of our results and genomic skimming/nanopore approach for conservation and evolutionary studies of C. aurita and other marmosets.

IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era

IQ-TREE (http://www.iqtree.org, last accessed February 6, 2020) is a user-friendly and widely used software package for phylogenetic inference using maximum likelihood. Since the release of version 1 in 2014, we have continuously expanded IQ-TREE to integrate a plethora of new models of sequence evolution and efficient computational approaches of phylogenetic inference to deal with genomic data. Here, we describe notable features of IQ-TREE version 2 and highlight the key advantages over other software.

Pervasive duplication of tumor suppressors in Afrotherians during the evolution of large bodies and reduced cancer risk

The risk of developing cancer is correlated with body size and lifespan within species. Between species, however, there is no correlation between cancer and either body size or lifespan, indicating that large, long-lived species have evolved enhanced cancer protection mechanisms. Elephants and their relatives (Proboscideans) are a particularly interesting lineage for the exploration of mechanisms underlying the evolution of augmented cancer resistance because they evolved large bodies recently within a clade of smaller-bodied species (Afrotherians). Here, we explore the contribution of gene duplication to body size and cancer risk in Afrotherians. Unexpectedly, we found that tumor suppressor duplication was pervasive in Afrotherian genomes, rather than restricted to Proboscideans. Proboscideans, however, have duplicates in unique pathways that may underlie some aspects of their remarkable anti-cancer cell biology. These data suggest that duplication of tumor suppressor genes facilitated the evolution of increased body size by compensating for decreasing intrinsic cancer risk.

AmAMP1 from Acropora millepora and damicornin define a family of coral-specific antimicrobial peptides related to the Shk toxins of sea anemones

A candidate antimicrobial peptide (AmAMP1) was identified by searching the whole genome sequence of Acropora millepora for short (<125AA) cysteine-rich predicted proteins with an N-terminal signal peptide but lacking clear homologs in the SwissProt database. It resembled but was not closely related to damicornin, the only other known AMP from a coral, and was shown to be active against both Gram-negative and Gram-positive bacteria. These proteins define a family of AMPs present in corals and their close relatives, the Corallimorpharia, and are synthesised as preproproteins in which the C-terminal mature peptide contains a conserved arrangement of six cysteine residues. Consistent with the idea of a common origin for AMPs and toxins, this Cys motif is shared between the coral AMPs and the Shk neurotoxins of sea anemones. AmAMP1 is expressed at late stages of coral development, in ectodermal cells that resemble the "ganglion neurons" of Hydra, in which it has recently been demonstrated that a distinct AMP known as NDA-1 is expressed.

The amylase gene cluster in house mice (Mus musculus) was subject to repeated introgression including the rescue of a pseudogene

Background

Amylase gene clusters have been implicated in adaptive copy number changes in response to the amount of starch in the diet of humans and mammals. However, this interpretation has been questioned for humans and for mammals there is a paucity of information from natural populations.

Results

Using optical mapping and genome read information, we show here that the amylase cluster in natural house mouse populations is indeed copy-number variable for Amy2b paralogous gene copies (called Amy2a1 - Amy2a5), but a direct connection to starch diet is not evident. However, we find that the amylase cluster was subject to introgression of haplotypes between Mus musculus sub-species. A very recent introgression can be traced in the Western European populations and this leads also to the rescue of an Amy2b pseudogene. Some populations and inbred lines derived from the Western house mouse (Mus musculus domesticus) harbor a copy of the pancreatic amylase (Amy2b) with a stop codon in the first exon, making it non-functional. But populations in France harbor a haplotype introgressed from the Eastern house mouse (M. m. musculus) with an intact reading frame. Detailed analysis of phylogenetic patterns along the amylase cluster suggest an additional history of previous introgressions.

Conclusions

Our results show that the amylase gene cluster is a hotspot of introgression in the mouse genome, making it an evolutionary active region beyond the previously observed copy number changes.