BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis

Martínez-Castilla L, García-Sandoval R, Ballester P, Ferrándiz C, Sanchez MDLP, García-Ponce B, Garay-Arroyo A, R. Álvarez-Buylla E. HDACs MADS-domain protein interaction: a case study of HDA15 and XAL1 in Arabidopsis thaliana

Cellular behavior, cell differentiation and ontogenetic development in eukaryotes result from complex interactions between epigenetic and classic molecular genetic mechanisms, with many of these interactions still to be elucidated. Histone deacetylase enzymes (HDACs) promote the interaction of histones with DNA by compacting the nucleosome, thus causing transcriptional repression. MADS-domain transcription factors are highly conserved in eukaryotes and participate in controlling diverse developmental processes in animals and plants, as well as regulating stress responses in plants. In this work, we focused on finding out putative interactions of Arabidopsis thaliana HDACs and MADS-domain proteins using an evolutionary perspective combined with bioinformatics analyses and testing the more promising predicted interactions through classic molecular biology tools. Through bioinformatic analyses, we found similarities between HDACs proteins from different organisms, which allowed us to predict a putative protein-protein interaction between the Arabidopsis thaliana deacetylase HDA15 and the MADS-domain protein XAANTAL1 (XAL1). The results of two-hybrid and Bimolecular Fluorescence Complementation analysis demonstrated in vitro and in vivo HDA15-XAL1 interaction in the nucleus. Likely, this interaction might regulate developmental processes in plants as is the case for this type of interaction in animals.

Streamlining large-scale oceanic biomonitoring using passive eDNA samplers integrated into vessel's continuous pump underway seawater systems

Passive samplers are enabling the scaling of environmental DNA (eDNA) biomonitoring in our oceans, by circumventing the time-consuming process of water filtration. Designing a novel passive sampler that does not require extensive sample handling time and can be connected to ocean-going vessels without impeding normal underway activities has potential to rapidly upscale global biomonitoring efforts onboard the world's oceanic fleet. Here, we demonstrate the utility of an artificial sponge sampler connected to the continuous pump underway seawater system as a means to enable oceanic biomonitoring. We compared the performance of this passive sampling protocol with standard water filtration at six locations during a research voyage from New Zealand to Antarctica in early 2023. Eukaryote metabarcoding of the mitochondrial COI gene revealed no significant difference in phylogenetic α-diversity between sampling methods and both methods delineated a progressive reduction in number of Zero-Radius Operational Taxonomic Units (ZOTUs) with increased latitudes. While both sampling methods revealed comparable trends in geographical community compositions, distinct clusters were identified for passive samplers and water filtration at each location. Additionally, greater variability between replicates was observed for passive samplers, resulting in an increased estimated level of replication needed to recover 90 % of the biodiversity. Furthermore, traditional water filtration failed to detect three phyla observed by passive samplers and extrapolation analysis estimated passive samplers recover a larger number of ZOTUs compared to water filtration for all six locations. Our results demonstrate the potential of this passive eDNA sampler protocol and highlight areas where this emerging technology could be improved, thereby enabling large-scale offshore marine eDNA biomonitoring by leveraging the world's oceanic fleet without interfering with onboard activities.

Utility of various genome lengths in diversity and evolution analyses of Hepatitis E virus

The aim of this study was to investigate to what extent fragments of the HEV genome could be used for accurate diagnostics and inference of viral population-scale processes. For this, we selected all the published whole genome sequences from the NCBI GenBank and trimmed them to various fragment lengths (ORF1,2,3, ORF1, ORF2, ORF3, 493 nt in ORF2 and 148 nt in ORF2). Each of the fragment lengths was used to infer the richness and diversity of the viral sequence types, typing accuracy, and potential use in phylodynamics. The results obtained from the different fragments were compared. We observed that, generally, the longer the nucleic acid fragment used in typing, the better the accuracy in predicting the viral subtype. However, the dominant HEV subtypes circulating in Europe were relatively well classified even by the 493 nt fragment, with false negative rates as low as 8 in 1000 typed sequences. Most fragments also give comparable results in analyses of population size, albeit with shorter fragments showing a broader 95 % highest posterior density interval and less obvious increase of the viral effective population size. The reconstructed phylogenies of a heterochronous subset indicated a good concordance between all the fragments, with the major clades following similar branching patterns. Furthermore, we have used the HEV sequence data from the Netherlands available in the HEVnet database as a case study for reconstruction of population size changes in the past decades. This data showed that molecular and epidemiological results are concordant and point to an increase in the viral effective population size underlying the observed increase in incidence of acute HEV infection cases. In the absence of whole genome sequencing data, the 493 bp fragment can be used for analyzing HEV strains currently circulating in Europe, as it is informative for describing short term population-scale processes.

Congruent patterns of cryptic cladogenesis revealed using RADseq and Sanger sequencing in a velvet worm species complex (Onychophora: Peripatopsidae: Peripatopsis sedgwicki)

In the present study, first generation DNA sequencing (mitochondrial cytochrome c oxidase subunit one, COI) and reduced-representative genomic RADseq data were used to understand the patterns and processes of diversification of the velvet worm, Peripatopsis sedgwicki species complex across its distribution range in South Africa. For the RADseq data, three datasets (two primary and one supplementary) were generated corresponding to 1,259-11,468 SNPs, in order to assess the diversity and phylogeography of the species complex. Tree topologies for the two primary datasets were inferred using maximum likelihood and Bayesian inferences methods. Phylogenetic analyses using the COI datasets retrieved four distinct, well-supported clades within the species complex. Five species delimitation methods applied to the COI data (ASAP, bPTP, bGMYC, STACEY and iBPP) all showed support for the distinction of the Fort Fordyce Nature Reserve specimens. In the main P. sedgwicki species complex, the species delimitation methods revealed a variable number of operational taxonomic units and overestimated the number of putative taxa. Divergence time estimates coupled with the geographic exclusivity of species and phylogeographic results suggest recent cladogenesis during the Plio/Pleistocene. The RADseq data were subjected to a principal components analysis and a discriminant analysis of principal components, under a maximum-likelihood framework. The latter results corroborate the four main clades observed using the COI data, however, applying additional filtering revealed additional diversity. The high overall congruence observed between the RADseq data and COI data suggest that first generation sequence data remain a cheap and effective method for evolutionary studies, although RADseq does provide a far greater resolution of contemporary temporo-spatial patterns.

Phylogenomics and biogeography of the small carpenter bees (Apidae: Xylocopinae: Ceratina)

Small carpenter bees in the genus Ceratina are behaviourally diverse, species-rich, and cosmopolitan, with over 370 species and a range including all continents except Antarctica. Here, we present the first comprehensive phylogeny of the genus based on ultraconserved element (UCE) phylogenomic data, covering a total of 185 ingroup specimens representing 22 of the 25 current subgenera. Our results support most recognized subgenera as natural groups, but we also highlight several groups in need of taxonomic revision - particularly the nominate subgenus Ceratina sensu stricto - and several clades that likely need to be described as new subgenera. In addition to phylogeny, we explore the evolutionary history of Ceratina through divergence time estimation and biogeographic reconstruction. Our findings suggest that Ceratinini split from its sister tribe Allodapini about 72 million years ago. The common ancestor of Ceratina emerged in the Afrotropical realm approximately 42 million years ago, near the Middle Eocene Climatic Optimum. Multiple subsequent dispersal events led to the present cosmopolitan distribution of Ceratina, with the majority of transitions occurring between the Afrotropics, Indomalaya, and the Palearctic. Additional movements also led to the arrival of Ceratina in Madagascar, Australasia, and a single colonization of the Americas. Dispersal events were asymmetrical overall, with temperate regions primarily acting as destinations for migrations from tropical source regions.

Plasmid-mediated acquisition and chromosomal integration of blaCTX-M-14 in a subclade of Escherichia coli ST131-H30 clade C1

Escherichia coli ST131 is a multidrug-resistant lineage associated with the global spread of extended-spectrum β-lactamase-producing organisms. Particularly, ST131 clade C1 is the most predominant clade in Japan, harboring blaCTX-M-14 at a high frequency. However, the process of resistance gene acquisition and spread remains unclear. Here, we performed whole-genome sequencing of 19 E. coli strains belonging to 12 STs and 12 fimH types collected between 1997 and 2016. Additionally, we analyzed the full-length genome sequences of 96 ST131-H30 clade C0 and C1 strains, including those obtained from this study and those registered in public databases, to understand how ST131 clade C1 acquired and spread blaCTX-M-14. We detected conjugative IncFII plasmids and IncB/O/K/Z plasmids carrying blaCTX-M-14 in diverse genetic lineages of E. coli strains from the 1990s to the 2010s, suggesting that these plasmids played an important role in the spread of blaCTX-M-14. Molecular phylogenetic and molecular clock analyses of the 96 ST131-H30 clade C0 and C1 strains identified 8 subclades. Strains harboring blaCTX-M-14 were clustered in subclades 4 and 5, and it was inferred that clade C1 acquired blaCTX-M-14 around 1993. All 34 strains belonging to subclade 5 possessed blaCTX-M-14 with ISEcp1 upstream at the same chromosomal position, indicating their common ancestor acquired blaCTX-M-14 in a single ISEcp1-mediated transposition event during the early formation of the subclade around 1999. Therefore, both the horizontal transfer of plasmids carrying blaCTX-M-14 to diverse genetic lineages and chromosomal integration in the predominant genetic lineage have contributed to the spread of blaCTX-M-14.

Phylogenetic relationships and historical biogeography of silkmoths (Lepidoptera: Bombycidae) suggest an origin in Southern Gondwana

Silkmoths (Bombycidae) have a disjunct distribution predominantly in the Southern Hemisphere and Asia. Here we reconstruct the phylogenetic history of the family to test competing hypotheses on their origin and assess how vicariance and long-distance dispersal shaped their current distribution. We sequenced up to 5,074 base pairs from six loci (COI, EF1-α, wgl, CAD, GAPDH, and RpS5) to infer the historical biogeography of Bombycidae. The multilocus dataset covering 20 genera (80 %) of the family, including 17 genera (94 %) of Bombycinae and 3 genera (43 %) of Epiinae, was used to estimate phylogenetic patterns, divergence times and biogeographic reconstruction. Dating estimates extrapolated from secondary calibration sources indicate the Bombycidae stem-group originated approximately 64 Mya. The subfamilies Epiinae (South America) and Bombycinae (Australia, Asia, East Palaearctic, and Africa) were reciprocally monophyletic, diverging at c. 56 Mya (95 % credibility interval: 66-46 Mya). The 'basal' lineage of Bombycinae - Gastridiota + Elachyophtalma - split from the rest of Bombycinae c. 53 Mya (95 % credibility interval: 63-43 Mya). Gastridiota is a monobasic genus with a relictual distribution in subtropical forests of eastern Australia. The Oriental and African genera comprised a monophyletic group: the Oriental region was inferred to have been colonized from a long-distance dispersal event from Australia to South-East Asia c. 53 Mya or possibly later (c. 36-26 Mya); Africa was subsequently colonized by dispersal from Asia c. 16 Mya (95 % credibility interval: 21-12 Mya). Based on the strongly supported phylogenetic relationships and estimates of divergence times, we conclude that Bombycidae had its origin in the fragment of Southern Gondwana consisting of Australia, Antarctica and South America during the Paleocene. The disjunction between South America (Epiinae) and Australia (Bombycinae) is best explained by vicariance in the Eocene, whereas the disjunct distribution in Asia and Africa is best explained by more recent dispersal events.

The last of their kind: Is the genus Scutiger (Anura: Megophryidae) a relict element of the paleo-transhimalaya biota?

The orographic evolution of the Himalaya-Tibet Mountain system continues to be a subject of controversy, leading to considerable uncertainty regarding the environment and surface elevation of the Tibetan Plateau during the Cenozoic era. As many geoscientific (but not paleontological) studies suggest, elevations close to modern heights exist in vast areas of Tibet since at least the late Paleogene, implicating the presence of large-scale alpine environments for more than 30 million years. To explore a recently proposed alternative model that assumes a warm temperate environment across paleo-Tibet, we carried out a phylogeographic survey using genomic analyses of samples covering the range of endemic lazy toads (Scutiger) across the Himalaya-Tibet orogen. We identified two main clades, with several, geographically distinct subclades. The long temporal gap between the stem and crown age of Scutiger may suggest high extinction rates. Diversification within the crown group, depending on the calibration, occurred either from the Mid-Miocene or Late-Miocene and continued until the Holocene. The present-day Himalayan Scutiger fauna could have evolved from lineages that existed on the southern edges of the paleo-Tibetan area (the Transhimalaya = Gangdese Shan), while extant species living on the eastern edge of the Plateau originated probably from the eastern edges of northern parts of the ancestral Tibetan area (Hoh Xil, Tanggula Shan). Based on the Mid-Miocene divergence time estimation and ancestral area reconstruction, we propose that uplift-associated aridification of a warm temperate Miocene-Tibet, coupled with high extirpation rates of ancestral populations, and species range shifts along drainage systems and epigenetic transverse valleys of the rising mountains, is a plausible scenario explaining the phylogenetic structure of Scutiger. This hypothesis aligns with the fossil record but conflicts with geoscientific concepts of high elevated Tibetan Plateau since the late Paleogene. Considering a Late-Miocene/Pliocene divergence time, an alternative scenario of dispersal from SE Asia into the East, Central, and West Himalaya cannot be excluded, although essential evolutionary and biogeographic aspects remain unresolved within this model.

Cretaceous amber inclusions illuminate the evolutionary origin of tardigrades

Tardigrades are a diverse phylum of microscopic invertebrates widely known for their extreme survival capabilities. Molecular clocks suggest that tardigrades diverged from other panarthropods before the Cambrian, but their fossil record is extremely sparse. Only the fossil tardigrades Milnesium swolenskyi (Late Cretaceous) and Paradoryphoribius chronocaribbeus (Miocene) have resolved taxonomic positions, restricting the availability of calibration points for estimating for the origin of this phylum. Here, we revise two crown-group tardigrades from Canadian Cretaceous-aged amber using confocal fluorescence microscopy, revealing critical morphological characters that resolve their taxonomic positions. Formal morphological redescription of Beorn leggi reveals that it features Hypsibius-type claws. We also describe Aerobius dactylus gen. et sp. nov. based on its unique combination of claw characters. Phylogenetic analyses indicate that Beo. leggi and Aer. dactylus belong to the eutardigrade superfamily Hypsibioidea, adding a critical fossil calibration point to investigate tardigrade origins. Our molecular clock estimates suggest an early Paleozoic diversification of crown-group Tardigrada and highlight the importance of Beo. leggi as a calibration point that directly impacts estimates of shallow nodes. Our results suggest that independent terrestrialization of eutardigrades and heterotardigrades occurred around the end-Carboniferous and Lower Jurassic, respectively. These estimates also provide minimum ages for convergent acquisition of cryptobiosis.

A targeted gene phylogenetic framework to investigate diversification in the highly diverse yet geographically restricted red devil spiders (Araneae, Dysderidae)

The family Dysderidae is a highly diverse group of nocturnal ground-dwelling and active-hunter spiders. Dysderids are mostly restricted to the Western Palearctic, and particularly rich and abundant around the Mediterranean region. Interestingly, the distribution of species richness among its 24 genera and three subfamilies is highly biased-80% of its 644 documented species belong to just two genera, Dysdera (326) and Harpactea (211). Dysderidae provides an excellent study case for evolutionary and ecological research. It includes cases of trophic specialization, which are uncommon among spiders, and exhibit other remarkable biological (e.g. holocentric chromosomes), behavioural (e.g. cryptic female choice), evolutionary (e.g. adaptive radiation) and ecological features (e.g. recurrent colonization of the subterranean environment). The lack of a quantitative hypothesis on its phylogenetic structure has hampered its potential as a testing ground for evolutionary, biogeographical and ecological hypotheses. Here, we present the results of a target, multi-locus phylogenetic analysis, using mitochondrial (cox1, 16s and 12s) and nuclear genes (h3, 28s and 18s), of the most exhaustive taxonomic sample within Dysderidae (104 spp.) to date and across related families (Synspermiata) (83 spp.). We estimate divergence times using a combination of fossil and biogeographic node calibrations and use this timeline to identify shifts in diversification rates. Our results support the monophyly of the Dysderidae subfamilies Rhodinae and Dysderinae but reject Harpacteinae as currently defined. Moreover, the clades recovered within Harpacteinae do not support its current taxonomy. The origin of the family most likely post-dated the break-up of Pangea, and cave colonization may be older than previously considered. After correcting for the taxonomic artefacts, we identified a significant shift in diversification rates at the base of the genus Dysdera. Although the unique coexistence of specialist and generalist diets within the lineage could be suggested as the potential driver for the rate acceleration, further quantitative analyses would be necessary to test this hypothesis.

Seek, and you will find: Cryptic diversity of the cardiopulmonary nematode Angiostrongylus vasorum in the Americas

Angiostrongylus vasorum is a metastrongylid parasite infecting wild canids and domestic dogs. Its patchy distribution, high pathogenicity and taxonomical classification makes the evolutionary history of A. vasorum intriguing and important to study. First larval stages of A. vasorum were recovered from feces of two grey foxes, Urocyon cinereoargenteus, from Costa Rica. Sequencing and phylogenetic and haplotypic analyses of the ITS2, 18S and cytochrome oxidase subunit 1 (cox1) fragments were performed. Then p- and Nei´s genetic distance, nucleotide substitution rates and species delimitation analyses were conducted with cox1 data of the specimens collected herein and other Angiostrongylus spp. Cophylogenetic congruence and coevolutionary events of Angiostrongylus spp. and their hosts were evaluated using patristic and phenetic distances and maximum parsimony reconciliations. Specimens from Costa Rica clustered in a separate branch from European and Brazilian A. vasorum sequences in the phylogenetic and haplotype network analyses using the ITS2 and cox1 data. In addition, cox1 p-distance of the sequences derived from Costa Rica were up to 8.6 % different to the ones from Europe and Brazil, a finding mirrored in Nei´s genetic distance PCoA. Species delimitation analysis supported a separate group with the sequences from Costa Rica, suggesting that these worms may represent cryptic variants of A. vasorum, a new undescribed taxon or Angiocaulus raillieti, a synonym species of A. vasorum described in Brazil. Moreover, nucleotide substitution rates in A. vasorum were up to six times higher than in the congener Angiostrongylus cantonensis. This finding and the long time elapsed since the last common ancestor between both species may explain the larger diversity in A. vasorum. Finally, cophylogenetic congruence was observed between Angiostrongylus spp. and their hosts, with cospeciation events occurring at deeper taxonomic branching of host order. Altogether, our data suggest that the diversity of the genus Angiostrongylus is larger than expected, since additional species may be circulating in wild canids from the Americas.

Genome analysis reveals three distinct lineages of the cosmopolitan white shark

The white shark (Carcharodon carcharias) (Linnaeus, 1758), an iconic apex predator occurring in all oceans,1,2 is classified as Vulnerable globally3-with global abundance having dropped to 63% of 1970s estimates,4-and as Critically Endangered in Europe.5 Identification of evolutionary significant units and their management are crucial for conservation,6 especially as the white shark is facing various but often region-specific anthropogenic threats.7,8,9,10,11 Assessing connectivity in a cosmopolitan marine species requires worldwide sampling and high-resolution genetic markers.12 Both are lacking for the white shark, with studies to date typified by numerous but geographically limited sampling, and analyses relying largely on relatively small numbers of nuclear microsatellites,13,14,15,16,17,18,19 which can be plagued by various genotyping artefacts and thus require cautious interpretation.20 Sequencing and computational advances are finally allowing genomes21,22,23 to be leveraged into population studies,24,25,26,27 with datasets comprising thousands of single-nucleotide polymorphisms (SNPs). Here, combining target gene capture (TGC)28 sequencing (89 individuals, 4,000 SNPs) and whole-genome re-sequencing (17 individuals, 391,000 SNPs) with worldwide sampling across most of the distributional range, we identify three genetically distinct allopatric lineages (North Atlantic, Indo-Pacific, and North Pacific). These diverged 100,000-200,000 years ago during the Penultimate Glaciation, when low sea levels, different ocean currents, and water temperatures produced significant biogeographic barriers. Our results show that without high-resolution genomic analyses of samples representative of a species' range,12 the true extent of diversity, presence of past and contemporary barriers to gene flow, subsequent speciation, and local evolutionary events will remain enigmatic.

Combining distribution modelling and phylogeography to understand present, past and future of an endangered spider

BACKGROUND

Understanding how endangered species respond to climatic changes is fundamental for their conservation. Due to its restricted geographic range, its sensitivity to the ongoing global warming and its continuing decline, the Southwestern-Alpine endemic wolf spider Vesubia jugorum is currently classified as Endangered in the IUCN Red List. Here, we combined species distribution modelling (SDM) and phylogeographic inference to describe the present, the past and the future of this species in light of the mtDNA genetic structure of extant populations.

RESULTS

Phylogenetic and network analyses show a high level of genetic differentiation and a strong genetic structure of the populations, likely explicable by a long history of isolation and survival in separate refugia. The SDM projection into past climatic conditions supports these results by showing a smaller distribution range compared to present, mostly restricted to the Maritime and Ligurian Alps, which possibly served as main refugium. Future forecast shows a significant shift in the bioclimatic range towards higher altitudes and latitudes, with a drastic decrease of habitat suitability in the central and south-eastern parts of the range, with consequent general loss of haplotype diversity.

CONCLUSION

SDM and phylogeographic inference support the hypothesis that the current distribution and the genetic structure of the extant populations mirror the survival in situ of Vesubia jugorum across repeated glacial and interglacial phases, in line with the 'long-term stability hypothesis'. Future predictions show a significant shift in the bioclimatic range that V. jugorum will be likely unable to track, with profound impact on its long-term survival and its genetic diversity. Our considerations have implication for conservation genetics, highlighting the pivotal role of the transboundary protected areas of the SW-Alps in promoting conservation efforts for this species.

Comparative mitochondrial genomics of endemic Mexican vesper yellow bats genus Rhogeessa (Chiroptera: Vespertilionidae) and insights into internal relationships in the family Vespertilionidae

In the species-rich family Vespertilionidae, vesper yellow bats in the genus Rhogeessa include eleven species, three of them endemic to Mexico. These insectivorous bats provide important ecosystem services, including pest control. Even though some aspects of their biology are well- known, only a few genomic resources are available for these species, which limits our understanding of their biology. In this study, we assembled and annotated the mitochondrial genome of four species: R. aenea, R. genowaysi, R. mira, and R. parvula. We generated a phylomitogenomic hypothesis based on translated protein-coding genes for a total of 52 species in the family Vespertilionidae and examined the phylogenetic position of the genus Rhogeessa and species within the family. The AT-rich mitogenomes of R. aenea, R. genowaysi, R. mira, and R. parvula are 16,763, 16,781, 16,807, and 16,794 pb in length, respectively. Each studied mitogenome encodes 13 Protein Coding Genes (PCGs), 22 transfer RNA genes, and 2 rRNA genes, and contains a putative control region (CR). All tRNAs exhibit a 'cloverleaf' secondary structure, except tRNA-Serine-1 that lacked the DHU arm in all studied mitogenomes. Selective pressure analyses indicated that all protein-coding genes are exposed to purifying selection. The phylomitogenomic analysis supported the monophyletic status of the family Vespertilionidae, confirmed the placement of Rhogeessa within the tribe Antrozoini, and clarified phylogenetic relationships within and among subfamilies and tribes in this family. Our results indicate that phylomitogenomics are useful to explore the evolutionary history of vesper bats. The assembly and comprehensive analysis of mitochondrial genomes offer the potential to generate molecular references and resources beneficial for genetic analyses aimed at understanding the ecology and evolution of these remarkable bats.

Phylogenomics of the pantropical Connaraceae: revised infrafamilial classification and the evolution of heterostyly

Connaraceae is a pantropical family of about 200 species containing lianas and small trees with remarkably diverse floral polymorphisms, including distyly, tristyly, homostyly, and dioecy. To date, relationships within the family have not been investigated using a targeted molecular phylogenetic treatment, severely limiting systematic understanding and reconstruction of trait evolution. Accordingly, their last infrafamilial classification was based only on morphological data. Here, we used phylogenomic data obtained using the Angiosperms353 nuclear target sequence capture probes, sampling all tribes and almost all genera, entirely from herbarium specimens, to revise infrafamilial classification and investigate the evolution of heterostyly. The backbone of the resulting molecular phylogenetic tree is almost entirely resolved. Connaraceae consists of two clades, one containing only the African genus Manotes (4 or 5 species), which we newly recognize at the subfamily level. Vegetative and reproductive synapomorphies are proposed for Manotoideae. Within Connaroideae, Connareae is expanded to include the former Jollydoreae. The backbone of Cnestideae, which contains more than half of the Connaraceae species, remains incompletely resolved. Reconstructions of reproductive system evolution are presented that tentatively support tristyly as the ancestral state for the family, with multiple parallel losses, in agreement with previous hypotheses, plus possible re-gains. However, the great diversity of stylar polymorphisms and their phylogenetic lability preclude a definitive answer. Overall, this study reinforces the usefulness of herbarium phylogenomics, and unlocks the reproductive diversity of Connaraceae as a model system for the evolution of complex biological phenomena.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00606-024-01909-y.

Mitogenomic recognition of incognito lineages in the mud spiny lobster Panulirus polyphagus (Herbst, 1793): A tale of unique genetic structuring and diversification

This study provides the first documentation of three deep conspecific lineages within Panulirus polyphagus in the Indian Ocean, bridging the gap in genetic research. Comparative mitogenomics between lineages (L) at both species and family levels, evolutionary relationships and heterogeneity of sequence divergence within Decapoda, and divergence time estimation were performed. The characterized mitogenomes ranged from 15,685-15,705 bp in size and exhibited a typical pancrustacean pattern. Among the three lineages, L1 predominated the Bay of Bengal, L2 the Arabian Sea, and L2.a, a less common lineage genetically closer to L2, was restricted to the latter region. A minor lineage L1.a, was observed in the Coral Triangle area. All PCGs displayed evidence of purifying selection across species and family levels. The largest genetic distance (K2P) between lineages was 9 %, notably between L1.a and L2.a. The phylogenetic tree subdivided the Achelates into Palinuridae and Scyllaridae, and the topology demonstrated a distinct pattern of lineage diversification within P. polyphagus. AliGROOVE analysis revealed no discernible divergence in Decapoda. The diversification of P. polyphagus appears to have occurred during Miocene, with further diversification in Pliocene. Furthermore, genetic stocks and population connectivity recognized here will provide valuable insight for spatial management planning of this dwindling resource.

Cryptic diversity and phylogeographic patterns of Mediodactylus species in the Eastern Mediterranean region

Cryptic diversity poses a great obstacle in our attempts to assess the current biodiversity crisis and may hamper conservation efforts. The gekkonid genus Mediodactylus, a well-known case of hidden species and genetic diversity, has been taxonomically reclassified several times during the last decade. Focusing on the Mediterranean populations, a recent study within the M. kotschyi species complex using classic mtDNA/nuDNA markers suggested the existence of five distinct species, some being endemic and some possibly threatened, yet their relationships have not been fully resolved. Here, we generated genome-wide SNPs (using ddRADseq) and applied molecular species delimitation approaches and population genomic analyses to further disentangle these relationships. Τhe most extensive nuclear dataset, so far, encompassing 2,360 loci and ∼ 699,000 bp from across the genome of Mediodactylus gecko, enabled us to resolve previously obscure phylogenetic relationships among the five, recently elevated, Mediodactylus species and to support the hypothesis that the taxon includes several new, undescribed species. Population genomic analyses within each of the proposed species showed strong genetic structure and high levels of genetic differentiation among populations.

Phylogenetic climatic niche evolution and diversification of the Neurergus species (Salamandridae) in the Irano-Anatolian biodiversity hotspot

This study explores how climate variables influenced the evolution and diversification of Neurergus newts within the Irano-Anatolian biodiversity hotspot. We use a dated phylogenetic tree and climatic niche models to analyze their evolutionary history and ecological preferences. Using genetic data from nuclear (KIAA) and mitochondrial (16s and 12s) genes, we estimate divergence times and identify four major Neurergus clades. The initial speciation event occurred approximately 11.3 million years ago, coinciding with the uplift of the Zagros and Anatolian mountains. This geological transformation isolated newt populations, likely triggering the first speciation event. By integrating potential geographic distribution with climate variables, we reconstruct ancestral niche occupancy profiles. This highlights the critical roles of temperature and precipitation in shaping Neurergus habitat preferences and distribution. We observe both phylogenetic niche conservatism and divergence, with niche divergence playing a dominant role in diversification. This research emphasizes the complex interplay of geography, climate, and ecology in speciation and the vulnerability of isolated mountain newt populations to environmental changes.

Application of molecular epidemiology in revealing HIV-1 transmission network and recombination patterns in Tianjin, China

Using a comprehensive molecular epidemiological approach, we characterized the transmission dynamics of HIV-1 among the MSM population in Tianjin, China. Our findings revealed that 38.56% (386/1001) of individuals clustered across 109 molecular transmission clusters (TCs), with MSM aged 50 and below being the group most commonly transmitting HIV-1. Among the identified TCs, CRF01_AE predominated, followed by CRF07_BC. Notably, CRF07_BC demonstrated a higher propensity for forming large clusters compared to CRF01_AE. Birth-death skyline analyses of the two largest clusters indicated that the HIV/AIDS transmission may be at a critical point, nearly all had Re approximately 1 by now. A retrospective analysis revealed that the rapid expansion of these large clusters was primarily driven by the introduction of viruses in 2021, highlighting the crucial importance of continuous molecular surveillance in identifying newly emerging high-risk transmission chains and adapting measures to address evolving epidemic dynamics. Furthermore, we detected the transmission of drug-resistant mutations (DRMs) within the TCs, particularly in the CRF07_BC clusters (K103N, Y181C, and K101E) and CRF01_AE clusters (P225H and K219R), emphasizing the importance of monitoring to support the continued efficacy of first-line therapies and pre-exposure prophylaxis (PrEP). Recombination analyses indicated that complex recombinant patterns, associated with increased amino acid variability, could confer adaptive traits to the viruses, potentially providing a competitive advantage in certain host populations or regions. Our study highlights the potential of integrating molecular epidemiological and phylodynamic approaches to inform targeted interventions.

Population structure of Dugong dugon across the Indo-Pacific revealed by historical mitogenomes

Sirenia, an iconic marine taxon with a tropical and subtropical worldwide distribution, face an uncertain future. All species are designated 'Vulnerable' to extinction by the IUCN. Nonetheless, a comprehensive understanding of geographic structuring across the global range is lacking, impeding our ability to highlight particularly vulnerable populations for conservation priority. Here, we use ancient DNA to investigate dugong (Dugong dugon) population structure, analysing 56 mitogenomes from specimens comprising the known historical range. Our results reveal geographically structured and distinct monophyletic clades characterized by contrasting evolutionary histories. We observe deep-rooted and divergent lineages in the East (Indo-Pacific) and obtain new evidence for the relatively recent dispersal of dugongs into the western Indian Ocean. All populations are significantly differentiated from each other with western populations having approximately 10-fold lower levels of genetic variation than eastern Indo-Pacific populations. Additionally, we find a significant temporal loss of genetic diversity in western Indian Ocean dugongs since the mid-twentieth century, as well as a decline in population size beginning approximately 1000 years ago. Our results add to the growing body of evidence that dugong populations are becoming ever more susceptible to ongoing human action and global climate change.

Molecular characterisation of Amblyomma integrum circulating in southern India

Amblyomma integrum is a large gooseberry sized longirostrate tick (when fully repleted) found in India and Sri Lanka. In Kerala (India), this tick is commonly found in the forest and its fringe areas frequently infesting deer and hence it is locally known as "maan chellu / maanunny" (deer tick). In the present study, molecular characterisation and phylogenetic analysis of A. integrum collected from the area grazed by the sambar deer (Rusa unicolor) of Kerala, south India was performed using three molecular markers viz., the mitochondrial cytochrome c oxidase subunit 1 (COI), mitochondrial 16S ribosomal RNA, and nuclear 18S ribosomal RNA genes. Cytochrome c oxidase subunit 1 (COI) gene showed better resolving ability for elucidating the evolutionary relationship of A. integrum and identified two distinct clades, viz., A and B. The Tamil Nadu isolates of south India and Marayoor isolate 1 (from Idukki district of Kerala bordering with Tamil Nadu) belonged to clade A. Majority of Wayanad isolates from Kerala, occupied clade B. The intraspecific genetic distance among the A. integrum species ranged from 0.00 to 13.34%. Between clades A and B, the genetic distance observed was 11.49%. The clade B isolates were genetically close to A. geoemydae (GD: 1.22%). Morphological variations between the clades included darker exoskeletal coloration in clade A and distinct differences in the shape of basis capitulum. Further analysis using Assemble Species by Automatic Partitioning (ASAP) and Generalized Mixed Yule Coalescent (GMYC) provided additional insights. Assemble Species by Automatic Partitioning (ASAP) identified 26 Molecular Operational Taxonomic Units (MOTUs) at a threshold distance of 5.38%, supporting the species partition of A. integrum clade B. Generalized Mixed Yule Coalescent (GMYC) analysis retained the same species complex (A. integrum-geoemydae Complex) inferred from the ASAP analyses. It could be inferred from the present study that the A. integrum clades A and B could be two different putative pseudocryptic species.

Exploring Paleogene Tibet's warm temperate environments through target enrichment and phylogenetic niche modelling of Himalayan spiny frogs (Paini, Dicroglossidae)

The Cenozoic topographic development of the Himalaya-Tibet orogen (HTO) substantially affected the paleoenvironment and biodiversity patterns of High Asia. However, concepts on the evolution and paleoenvironmental history of the HTO differ massively in timing, elevational increase and sequence of surface uplift of the different elements of the orogen. Using target enrichment of a large set of transcriptome-derived markers, ancestral range estimation and paleoclimatic niche modelling, we assess a recently proposed concept of a warm temperate paleo-Tibet in Asian spiny frogs of the tribe Paini and reconstruct their historical biogeography. That concept was previously developed in invertebrates. Because of their early evolutionary origin, low dispersal capacity, high degree of local endemism, and strict dependence on temperature and humidity, the cladogenesis of spiny frogs may echo the evolution of the HTO paleoenvironment. We show that diversification of main lineages occurred during the early to Mid-Miocene, while the evolution of alpine taxa started during the late Miocene/early Pliocene. Our distribution and niche modelling results indicate range shifts and niche stability that may explain the modern disjunct distributions of spiny frogs. They probably maintained their (sub)tropical or (warm)temperate preferences and moved out of the ancestral paleo-Tibetan area into the Himalaya as the climate shifted, as opposed to adapting in situ. Based on ancestral range estimation, we assume the existence of low-elevation, climatically suitable corridors across paleo-Tibet during the Miocene along the Kunlun, Qiangtang and/or Gangdese Shan. Our results contribute to a deeper understanding of the mechanisms and processes of faunal evolution in the HTO.

Diversity and biogeography of South American mud turtles elucidated by multilocus DNA sequencing (Testudines: Kinosternidae)

Kinosternon is the most speciose genus of extant turtles, with 22 currently recognized species, distributed across large parts of the Americas. Most species have small distributions, but K. leucostomum and K. scorpioides range from Mexico to South America. Previous studies have found discordance between mitochondrial and nuclear phylogenies in some kinosternid groups, with the current taxonomy following the nuclear-based results. Herein, based on extended molecular, geographic, and taxonomic sampling, we explore the phylogeographic structure and taxonomic limits for K. leucostomum and the K. scorpioides group and present a fossil-calibrated nuclear time tree for Kinosternon. Our results reveal contrasting differentiation patterns for the K. scorpioides group and K. leucostomum, despite overlapping distributions. Kinosternon leucostomum shows only shallow geographic divergence, whereas the K. scorpioides group is polyphyletic with up to 10 distinct taxa, some of them undescribed. We support the elevation of K. s. albogulare and K. s. cruentatum to species level. Given the deep divergence within the genus Kinosternon, we propose the recognition of three subgenera, Kinosternon, Cryptochelys and Thyrosternum, and the abandonment of the group-based classification, at least for the K. leucostomum and K. scorpioides groups. Our results show an initial split in Kinosternon that gave rise to two main radiations, one Nearctic and one mainly Neotropical. Most speciation events in Kinosternon occurred during the Quaternary and we hypothesize that they were mediated by both climatic and geological events. Additionally, our data imply that at least three South American colonizations occurred, two in the K. leucostomum group, and one in the K. scorpioides group. Additionally, we hypothesize that discordance between mitochondrial and nuclear phylogenetic signal is due to mitochondrial capture from an extinct kinosternine lineage.

Clinging on the brink: Whole genomes reveal human-induced population declines and severe inbreeding in the Critically Endangered Emirati Leaf-toed Gecko (Asaccus caudivolvulus)

Human-mediated habitat destruction has had a profound impact on increased species extinction rates and population declines worldwide. The coastal development in the United Arab Emirates (UAE) over the last two decades, serves as an example of how habitat transformation can alter the landscape of a country in just a few years. Here, we study the genomic implications of habitat transformation in the Critically Endangered Emirati Leaf-toed Gecko (Asaccus caudivolvulus), the only endemic vertebrate of the UAE. We generate a high-quality reference genome for this gecko, representing the first reference genome for the family Phyllodactylidae, and produce whole-genome resequencing data for 23 specimens from 10 different species of leaf-toed geckos. Our results show that A. caudivolvulus has consistently lower genetic diversity than any other Arabian species of Asaccus, suggesting a history of ancient population declines. However, high levels of recent inbreeding are recorded among populations in heavily developed areas, with a more than 50% increase in long runs of homozygosity within a 9-year period. Moreover, results suggest that this species does not effectively purge deleterious mutations, hence making it more vulnerable to future stochastic threats. Overall, results show that A. caudivolvulus is in urgent need of protection, and habitat preservation must be warranted to ensure the species' survival.

Niche breadth and divergence in sympatric cryptic coral species (Pocillopora spp.) across habitats within reefs and among algal symbionts

While the presence of morphologically cryptic species is increasingly recognized, we still lack a useful understanding of what causes and maintains co-occurring cryptic species and its consequences for the ecology, evolution, and conservation of communities. We sampled 724 Pocillopora corals from five habitat zones (the fringing reef, back reef, and fore reef at 5, 10, and 20 m) at four sites around the island of Moorea, French Polynesia. Using validated genetic markers, we identified six sympatric species of Pocillopora, most of which cannot be reliably identified based on morphology: P. meandrina (42.9%), P. tuahiniensis (25.1%), P. verrucosa (12.2%), P. acuta (10.4%), P. grandis (7.73%), and P. cf. effusa (2.76%). For 423 colonies (58% of the genetically identified hosts), we also used psbA ncr or ITS2 markers to identify symbiont species (Symbiodiniaceae). The relative abundance of Pocillopora species differed across habitats within the reef. Sister taxa P. verrucosa and P. tuahiniensis had similar niche breadths and hosted the same specialist symbiont species (mostly Cladocopium pacificum) but the former was more common in the back reef and the latter more common deeper on the fore reef. In contrast, sister taxa P. meandrina and P. grandis had the highest niche breadths and overlaps and tended to host the same specialist symbiont species (mostly C. latusorum). Pocillopora acuta had the narrowest niche breadth and hosted the generalist, and more thermally tolerant, Durusdinium gynnii. Overall, there was a positive correlation between reef habitat niche breadth and symbiont niche breadth-Pocillopora species with a broader habitat niche also had a broader symbiont niche. Our results show how fine-scale variation within reefs plays an important role in the generation and coexistence of cryptic species. The results also have important implications for how niche differences affect community resilience, and for the success of coral restoration practices, in ways not previously appreciated.

Phylogenomics of mulberries (Morus, Moraceae) inferred from plastomes and single copy nuclear genes

Mulberries (genus Morus), belonging to the order Rosales, family Moraceae, are important woody plants due to their economic values in sericulture, as well as for nutritional benefits and medicinal values. However, the taxonomy and phylogeny of Morus, especially for the Asian species, remains challenging due to its wide geographical distribution, morphological plasticity, and interspecific hybridization. To better understand the evolutionary history of Morus, we combined plastomes and a large-scale nuclear gene analyses to investigate their phylogenetic relationships. We assembled the plastomes and screened 211 single-copy nuclear genes from 13 Morus species and related taxa. The plastomes of Morus species were relatively conserved in terms of genome size, gene content, synteny, IR boundary and codon usage. Using nuclear data, our results elucidated identical topologies based on coalescent and concatenation methods. The genus Morus was supported as monophyletic, with M. notabilis as the first diverging lineage and the two North American Morus species, M. microphylla and M. rubra, as sister to the other Asian species. In the Asian Morus species, interspecific relationships were completely resolved. However, cyto-nuclear discordances and gene tree-species tree conflicts were detected in the phylogenies of Morus, with multiple evidences supporting hybridization/introgression as the main cause of discordances between nuclear and plastid phylogenies, while gene tree-species tree conflicts were mainly caused by ILS.

Substrate specialisation drives an unexpectedly diverse radiation in barking geckos (Ptenopus: Gekkonidae)

Barking geckos (genus Ptenopus) are terrestrial, burrowing lizards endemic to southern Africa, currently with three recognised species. Two species are range-restricted (P. kochi and P. carpi) and display clear differences in substrate preference (soft sand vs. hard gravel). The third and most widespread species, P. garrulus, occurs on a variety of substrates of differing hardness, across potential geographic barriers, and over a steep climatic gradient. Variations in morphology and advertisement calls indicates that P. garrulus may be a species complex. Two subspecies of P. garrulus are currently recognised: P. g. maculatus and P. g. garrulus. To investigate species boundaries, we produced the first comprehensive phylogeny for the genus. We used a novel application of multiple regression on matrices models to assess multiple environmental drivers of diversification, as contrasted to isolation by distance. We show that P. kochi, P. carpi, and P. g. garrulus are valid species, but that P. g. maculatus is a paraphyletic complex of five previously unrecognised taxa. Specialisation onto different substrates was likely the main driver of divergence, with parapatric occurrence of two to four clades occurring at each of the three substrate transition zones identified a priori. The region encompasses diverse bioclimatic regions and potential geographic barriers, and these likely played a role in some divergence events.

Phylogeography of Parnassius citrinarius Based on Mitochondrial DNA Reveals Large Differences in Genetic Structure between the Eastern and Western Japan

The Japanese Archipelago hosts a rich butterfly fauna, and elucidating the genetic structures of multiple species is necessary to clarify their formation processes. This study aimed to reveal the genetic structure and distribution formation process of Parnassius citrinarius, which is widely distributed across the Japanese Archipelago from Hokkaido to Shikoku, through phylogeographic analysis based on the mitochondrial cytochrome c oxidase subunit I (COI) gene sequence. Thirty haplotypes were revealed from 311 individuals from 47 sites, indicating significant differences in the genetic structures between the eastern and western parts of the Japanese Archipelago. In Eastern Japan, multiple genetic clusters were found, with some sites harboring two clusters. The divergence times among populations in Eastern Japan were relatively recent, and no genetic differentiation was observed between regions, including between Hokkaido and Honshu, which are separated by a narrow strait. In contrast, in Western Japan, including Shikoku, unique genetic clusters were observed in each region. The phylogenetic relationships among populations were regionally clustered, and the divergence times were relatively ancient. The distribution and genetic structure of P. citrinarius in the Japanese Archipelago have been significantly influenced by temperature fluctuations and the presence of geographical barriers during the Pleistocene glacial-interglacial cycles, including the potential formation of refugia in Western Japan.

Selection against domestication alleles in introduced rabbit populations

Humans have moved domestic animals around the globe for thousands of years. These have occasionally established feral populations in nature, often with devastating ecological consequences. To understand how natural selection shapes re-adaptation into the wild, we investigated one of the most successful colonizers in history, the European rabbit. By sequencing the genomes of 297 rabbits across three continents, we show that introduced populations exhibit a mixed wild-domestic ancestry. We show that alleles that increased in frequency during domestication were preferentially selected against in novel natural environments. Interestingly, causative mutations for common domestication traits sometimes segregate at considerable frequencies if associated with less drastic phenotypes (for example, coat colour dilution), whereas mutations that are probably strongly maladaptive in nature are absent. Whereas natural selection largely targeted different genomic regions in each introduced population, some of the strongest signals of parallelism overlap genes associated with neuronal or brain function. This limited parallelism is probably explained by extensive standing genetic variation resulting from domestication together with the complex mixed ancestry of introduced populations. Our findings shed light on the selective and molecular mechanisms that enable domestic animals to re-adapt to the wild and provide important insights for the mitigation and management of invasive populations.

One hundred years of influenza A evolution

Leveraging the simplicity of nucleotide mismatch distributions, we provide an intuitive window into the evolution of the human influenza A 'nonstructural' (NS) gene segment. In an analysis suggested by the eminent Danish biologist Freddy B. Christiansen, we illustrate the existence of a continuous genetic "backbone" of influenza A NS sequences, steadily increasing in nucleotide distance to the 1918 root over more than a century. The 2009 influenza A/H1N1 pandemic represents a clear departure from this enduring genetic backbone. Utilizing nucleotide distance maps and phylogenetic analyses, we illustrate remaining uncertainties regarding the origin of the 2009 pandemic, highlighting the complexity of influenza evolution. The NS segment is interesting precisely because it experiences less pervasive positive selection, and departs less strongly from neutral evolution than e.g. the HA antigen. Consequently, sudden deviations from neutral diversification can indicate changes in other genes via the hitchhiking effect. Our approach employs two measures based on nucleotide mismatch counts to analyze the evolutionary dynamics of the NS gene segment. The rooted Hamming map of distances between a reference sequence and all other sequences over time, and the unrooted temporal Hamming distribution which captures the distribution of genotypic distances between simultaneously circulating viruses, thereby revealing patterns of nucleotide diversity and epi-evolutionary dynamics.

Shared alleles and genetic structures in different Thai domestic cat breeds: the possible influence of common racial origins

Over hundreds of years, cats have been domesticated and selectively bred, resulting in numerous pedigreed breeds expedited by recent cat shows and breeding associations. Concerns have been raised about the limited breeding options and the genetic implications of inbreeding, indicating challenges in maintaining genetic diversity and accurate identification in purebred cats. In this study, genetic variability and structure were examined in 5 Thai domestic cat breeds using 15 microsatellite markers and mitochondrial DNA (mtDNA) D-loop sequencing. In total, 184 samples representing the Wichien Maat (WCM), Suphalak (SL), Khao-Manee (KM), Korat (KR), and Konja (KJ) breeds were analyzed. High genetic diversity (Ho and He > 0.5) was observed in all breeds, and mtDNA analysis revealed two primary haplogroups (A and B) that were shared among all domestic cat breeds in Thailand and globally. However, minor differences were observed between Thai domestic cat breeds based on clustering analyses, in which a distinct genetic structure was observed in the WCM breed. This suggests that allele fixation for distinctive morphological traits has occurred in Thai domestic cat breeds that emerged in isolated regions with shared racial origins. Analysis of relationships among individuals within the breed revealed high identification efficiency in Thai domestic cat breeds (P(ID)sibs < 10-4). Additionally, diverse and effective individual identification can be ensured by optimizing marker efficiency by using only nine loci. This comprehensive genetic characterization provides valuable insights into conservation strategies and breeding practices for Thai domestic cat breeds.

Investigating linguistic and genetic shifts in East Indian tribal groups

South Asia is home to almost a quarter of the world's total population and is home to significant ethnolinguistic diversity. Previous studies of linguistic and genetic affiliations of Indian populations suggest that the formation of these distinct groups was a protracted and complex phenomenon involving multiple waves of migration, cultural assimilation, and genetic admixture. The evolutionary processes of migration, mixing and merging of populations thus impact the culture and linguistic diversity of different groups, some of which may retain their linguistic affinities despite genetic admixture with other groups, or vice versa. Our study examines the relationship of genetic and linguistic affinities between Austroasiatic and Indo-European speakers in adjacent geographical regions of Eastern India. We analyzed 224 mitogenomes and 0.65 million SNP genotypes from 40 unrelated individuals belonging to the Bathudi, Bhumij, Ho, and Mahali ethnic groups from the Eastern Indian state of Odisha. These four groups are speakers of Austroasiatic languages who have adopted elements from Indo-European languages spoken in neighbouring regions. Our results suggest that these groups have the greatest maternal genetic affinity with other Austroasiatic-speaking groups in India. Allele frequency-based analyses, genome-wide SNPs, haplotype-based methods and IBD sharing further support the genetic similarity of these East Indian groups to Austroasiatic speakers of South Asia rather than regional populations speaking Indo-European and Dravidian languages. Our study shows that these populations experienced linguistic mixing, likely due to industrialization and modernization that brought them into close cultural contact with neighbouring Indo-European-speaking groups. However, linguistic change in these groups is not reflected in genetic mixing in these populations, as they appear to maintain strict genetic boundaries while simultaneously experiencing cultural mixing.

Phylogenetic insights into the Salicaceae: The evolution of willows and beyond

The Salicaceae includes approximately 54 genera and over 1,400 species with a cosmopolitan distribution. Members of the family are well-known for their diverse secondary plant metabolites, and they play crucial roles in tropical and temperate forest ecosystems. Phylogenetic reconstruction of the Salicaceae has been historically challenging due to the limitations of molecular markers and the extensive history of hybridization and polyploidy within the family. Our study employs whole-genome sequencing of 74 species to generate an extensive phylogeny of the Salicaceae. We generated two RAD-Seq enriched whole-genome sequence datasets and extracted two additional gene sets corresponding to the universal Angiosperms353 and Salicaceae-specific targeted-capture arrays. We reconstructed maximum likelihood-based molecular phylogenies using supermatrix and coalescent-based supertree approaches. Our fossil-calibrated phylogeny estimates that the Salicaceae originated around 128 million years ago and unravels the complex taxonomic relationships within the family. Our findings confirm the non-monophyly of the subgenus Salix s.l. and further support the merging of subgenera Chamaetia and Vetrix, both of which exhibit intricate patterns within and among different sections. Overall, our study not only enhances our understanding of the evolution of the Salicaceae, but also provides valuable insights into the complex relationships within the family.

MAST: Phylogenetic Inference with Mixtures Across Sites and Trees

Hundreds or thousands of loci are now routinely used in modern phylogenomic studies. Concatenation approaches to tree inference assume that there is a single topology for the entire dataset, but different loci may have different evolutionary histories due to incomplete lineage sorting (ILS), introgression, and/or horizontal gene transfer; even single loci may not be treelike due to recombination. To overcome this shortcoming, we introduce an implementation of a multi-tree mixture model that we call mixtures across sites and trees (MAST). This model extends a prior implementation by Boussau et al. (2009) by allowing users to estimate the weight of each of a set of pre-specified bifurcating trees in a single alignment. The MAST model allows each tree to have its own weight, topology, branch lengths, substitution model, nucleotide or amino acid frequencies, and model of rate heterogeneity across sites. We implemented the MAST model in a maximum-likelihood framework in the popular phylogenetic software, IQ-TREE. Simulations show that we can accurately recover the true model parameters, including branch lengths and tree weights for a given set of tree topologies, under a wide range of biologically realistic scenarios. We also show that we can use standard statistical inference approaches to reject a single-tree model when data are simulated under multiple trees (and vice versa). We applied the MAST model to multiple primate datasets and found that it can recover the signal of ILS in the Great Apes, as well as the asymmetry in minor trees caused by introgression among several macaque species. When applied to a dataset of 4 Platyrrhine species for which standard concatenated maximum likelihood (ML) and gene tree approaches disagree, we observe that MAST gives the highest weight (i.e., the largest proportion of sites) to the tree also supported by gene tree approaches. These results suggest that the MAST model is able to analyze a concatenated alignment using ML while avoiding some of the biases that come with assuming there is only a single tree. We discuss how the MAST model can be extended in the future.

Bayesian Phylogenetic Analysis on Multi-Core Compute Architectures: Implementation and Evaluation of BEAGLE in RevBayes With MPI

Phylogenies are central to many research areas in biology and commonly estimated using likelihood-based methods. Unfortunately, any likelihood-based method, including Bayesian inference, can be restrictively slow for large datasets-with many taxa and/or many sites in the sequence alignment-or complex substitutions models. The primary limiting factor when using large datasets and/or complex models in probabilistic phylogenetic analyses is the likelihood calculation, which dominates the total computation time. To address this bottleneck, we incorporated the high-performance phylogenetic library BEAGLE into RevBayes, which enables multi-threading on multi-core CPUs and GPUs, as well as hardware specific vectorized instructions for faster likelihood calculations. Our new implementation of RevBayes+BEAGLE retains the flexibility and dynamic nature that users expect from vanilla RevBayes. In addition, we implemented native parallelization within RevBayes without an external library using the message passing interface (MPI); RevBayes+MPI. We evaluated our new implementation of RevBayes+BEAGLE using multi-threading on CPUs and 2 different powerful GPUs (NVidia Titan V and NVIDIA A100) against our native implementation of RevBayes+MPI. We found good improvements in speedup when multiple cores were used, with up to 20-fold speedup when using multiple CPU cores and over 90-fold speedup when using multiple GPU cores. The improvement depended on the data type used, DNA or amino acids, and the size of the alignment, but less on the size of the tree. We additionally investigated the cost of rescaling partial likelihoods to avoid numerical underflow and showed that unnecessarily frequent and inefficient rescaling can increase runtimes up to 4-fold. Finally, we presented and compared a new approach to store partial likelihoods on branches instead of nodes that can speed up computations up to 1.7 times but comes at twice the memory requirements.

Phylogeny and disease associations of a widespread and ancient intestinal bacteriophage lineage

Viruses are core components of the human microbiome, impacting health through interactions with gut bacteria and the immune system. Most human microbiome viruses are bacteriophages, which exclusively infect bacteria. Until recently, most gut virome studies focused on low taxonomic resolution (e.g., viral operational taxonomic units), hampering population-level analyses. We previously identified an expansive and widespread bacteriophage lineage in inhabitants of Amsterdam, the Netherlands. Here, we study their biodiversity and evolution in various human populations. Based on a phylogeny using sequences from six viral genome databases, we propose the Candidatus order Heliusvirales. We identify heliusviruses in 82% of 5441 individuals across 39 studies, and in nine metagenomes from humans that lived in Europe and North America between 1000 and 5000 years ago. We show that a large lineage started to diversify when Homo sapiens first appeared some 300,000 years ago. Ancient peoples and modern hunter-gatherers have distinct Ca. Heliusvirales populations with lower richness than modern urbanized people. Urbanized people suffering from type 1 and type 2 diabetes, as well as inflammatory bowel disease, have higher Ca. Heliusvirales richness than healthy controls. We thus conclude that these ancient core members of the human gut virome have thrived with increasingly westernized lifestyles.

Convergent Adaptation of True Crabs (Decapoda: Brachyura) to a Gradient of Terrestrial Environments

For much of terrestrial biodiversity, the evolutionary pathways of adaptation from marine ancestors are poorly understood and have usually been viewed as a binary trait. True crabs, the decapod crustacean infraorder Brachyura, comprise over 7600 species representing a striking diversity of morphology and ecology, including repeated adaptation to non-marine habitats. Here, we reconstruct the evolutionary history of Brachyura using new and published sequences of 10 genes for 344 tips spanning 88 of 109 brachyuran families. Using 36 newly vetted fossil calibrations, we infer that brachyurans most likely diverged in the Triassic, with family-level splits in the late Cretaceous and early Paleogene. By contrast, the root age is underestimated with automated sampling of 328 fossil occurrences explicitly incorporated into the tree prior, suggesting such models are a poor fit under heterogeneous fossil preservation. We apply recently defined trait-by-environment associations to classify a gradient of transitions from marine to terrestrial lifestyles. We estimate that crabs left the marine environment at least 7 and up to 17 times convergently, and returned to the sea from non-marine environments at least twice. Although the most highly terrestrial- and many freshwater-adapted crabs are concentrated in Thoracotremata, Bayesian threshold models of ancestral state reconstruction fail to identify shifts to higher terrestrial grades due to the degree of underlying change required. Lineages throughout our tree inhabit intertidal and marginal marine environments, corroborating the inference that the early stages of terrestrial adaptation have a lower threshold to evolve. Our framework and extensive new fossil and natural history datasets will enable future comparisons of non-marine adaptation at the morphological and molecular level. Crabs provide an important window into the early processes of adaptation to novel environments, and different degrees of evolutionary constraint that might help predict these pathways. [Brachyura; convergent evolution; crustaceans; divergence times; fossil calibration; molecular phylogeny; terrestrialization; threshold model.].

Insights from the first phylogenomic analysis of flat wasps (Hymenoptera, Bethylidae) reveal two new subfamilies

Despite significant advances in alpha level taxonomy in the past few decades, the higher-level phylogeny of flat wasps (Hymenoptera, Bethylidae) remains poorly explored. Herein we provide the first phylogenomic assessment of the family based on data from ultraconserved elements for 96 species in 61 genera of the family, with material from 29 countries and all biogeographic regions. Dataset cutoffs including ultraconserved element loci recovered for 50 and 70% of terminals resulted in matrices with 1513 and 451 loci, which were analysed in both parsimony and maximum likelihood frameworks. We also provide the first analyses of divergence dating for the family based on the calibration of 12 nodes. All analyses recovered the Bethylidae as a monophyletic group and estimated the origin of the family at around 143 Ma. The results suggest that all extant subfamilies had already diversified by the Late Cretaceous. All topologies suggest that Glenosema and Chilepyris form a clade separate from other Scleroderminae; owing to the morphological distinctiveness of each genus, we propose that they are accommodated in two new subfamilies, Glenoseminae subf. nov. and Chilepyrinae subf. nov. The monotypic sclerodermine genus Galodoxa was consistently recovered within Epyrinae and is transferred to the latter subfamily.

Phylogenetic and Comparative Genomics Study of Papilionidae Based on Mitochondrial Genomes

Most species of Papilionidae are large and beautiful ornamental butterflies. They are recognized as model organisms in ecology, evolutionary biology, genetics, and conservation biology but present numerous unresolved phylogenetic problems. Complete mitochondrial genomes (mitogenomes) have been widely used in phylogenetic studies of butterflies, but mitogenome knowledge within the family Papilionidae is limited, and its phylogeny is far from resolved. In this study, we first report the mitogenome of Byasa confusa from the subfamily Papilioninae of Papilionidae. The mitogenome of B. confusa is 15,135 bp in length and contains 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and an AT-rich control region (CR), closely mirroring the genomic structure observed in related butterfly species. Comparative analysis of 77 Papilionidae mitogenomes shows gene composition and order to be identical to that of an ancestral insect, and the AT bias, Ka/Ks, and relative synonymous codon usage (RSCU) are all consistent with that of other reported butterfly mitogenomes. We conducted phylogenetic analyses using maximum-likelihood (ML) and Bayesian-inference (BI) methods, with 77 Papilionidae species as ingroups and two species of Nymphalidae and Lycaenidae as outgroups. The phylogenetic analysis indicated that B. confusa were clustered within Byasa. The phylogenetic trees show the monophyly of the subfamily Papilioninae and the tribes Leptocircini, Papilionini, and Troidini. The data supported the following relationships in tribe level on Papilioninae: (((Troidini + Papilionini) + Teinopalpini) + Leptocircini). The divergence time analysis suggests that Papilionidae originated in the late Creataceous. Overall, utilizing the largest number of Papilionidae mitogenomes sequenced to date, with the current first exploration in a phylogenetic analysis on Papilionidae (including four subfamilies), this study comprehensively reveals the mitogenome characteristics and mitogenome-based phylogeny, providing information for further studies on the mitogenome, phylogeny, evolution, and taxonomic revision of the Papilionidae family.

Global Update on Measles Molecular Epidemiology

Molecular surveillance of circulating measles variants serves as a line of evidence for the absence of endemic circulation and provides a means to track chains of transmission. Molecular surveillance for measles (genotyping) is based on the sequence of 450 nucleotides at the end of the nucleoprotein coding region (N450) of the measles genome. Genotyping was established in 1998 and, with over 50,000 sequence submissions to the Measles Nucleotide Surveillance database, has proven to be an effective resource for countries attempting to trace pathways of transmission. This review summarizes the tools used for the molecular surveillance of measles and describes the challenge posed by the decreased number of circulating measles genotypes. The Global Measles and Rubella Laboratory Network addressed this challenge through the development of new tools such as named strains and distinct sequence identifiers that analyze the diversity within the currently circulating genotypes. The advantages and limitations of these approaches are discussed, together with the need to generate additional sequence data including whole genome sequences to ensure the continued utility of strain surveillance for measles.

Waves of Colonization and Gene Flow in a Great Speciator

Secondary contact between previously allopatric lineages offers a test of reproductive isolating mechanisms that may have accrued in isolation. Such instances of contact can produce stable hybrid zones-where reproductive isolation can further develop via reinforcement or phenotypic displacement-or result in the lineages merging. Ongoing secondary contact is most visible in continental systems, where steady input from parental taxa can occur readily. In oceanic island systems, however, secondary contact between closely related species of birds is relatively rare. When observed on sufficiently small islands, relative to population size, secondary contact likely represents a recent phenomenon. Here, we examine the dynamics of a group of birds whose apparent widespread hybridization influenced Ernst Mayr's foundational work on allopatric speciation: the whistlers of Fiji (Aves: Pachycephala). We demonstrate two clear instances of secondary contact within the Fijian archipelago, one resulting in a hybrid zone on a larger island, and the other resulting in a wholly admixed population on a smaller, adjacent island. We leveraged low genome-wide divergence in the hybrid zone to pinpoint a single genomic region associated with observed phenotypic differences. We use genomic data to present a new hypothesis that emphasizes rapid plumage evolution and post-divergence gene flow.

Stomatal dynamics in Alloteropsis semialata arise from the evolving interplay between photosynthetic physiology, stomatal size and biochemistry

C4 plants are expected to have faster stomatal movements than C3 species because they tend to have smaller guard cells. However, little is known about how the evolution of C4 photosynthesis influences stomatal dynamics in relation to guard cell size and environmental factors. We studied photosynthetically diverse populations of the grass Alloteropsis semialata, showing that the origin of C4 photosynthesis in this species was associated with a shortening of stomatal guard and subsidiary cells. However, for a given cell size, C4 and C3-C4 intermediate individuals had similar or slower light-induced stomatal opening speeds than C3 individuals. Conversely, when exposed to decreasing light, stomata in C4 plants closed as fast as those in non-C4 plants. Polyploid formation in some C4 plants led to larger stomatal cells and was associated with slower stomatal opening. Conversely, diversification of C4 diploid plants into wetter environments was associated with an acceleration of stomatal opening. Overall, there was significant relationship between light-saturated photosynthesis and stomatal opening speed in the C4 plants, implying that photosynthetic energy production was limiting for stomatal opening. Stomatal dynamics in this wild grass therefore arise from the evolving interplay between photosynthetic physiology and the size and biochemical function of stomatal complexes.

An 8000 years old genome reveals the Neolithic origin of the zoonosis Brucella melitensis

Brucella melitensis is a major livestock bacterial pathogen and zoonosis, causing disease and infection-related abortions in small ruminants and humans. A considerable burden to animal-based economies today, the presence of Brucella in Neolithic pastoral communities has been hypothesised but we lack direct genomic evidence thus far. We report a 3.45X B. melitensis genome preserved in an ~8000 year old sheep specimen from Menteşe Höyük, Northwest Türkiye, demonstrating that the pathogen had evolved and was circulating in Neolithic livestock. The genome is basal with respect to all known B. melitensis and allows the calibration of the B. melitensis speciation time from the primarily cattle-infecting B. abortus to approximately 9800 years Before Present (BP), coinciding with a period of consolidation and dispersal of livestock economies. We use the basal genome to timestamp evolutionary events in B. melitensis, including pseudogenization events linked to erythritol response, the supposed determinant of the pathogen's placental tropism in goats and sheep. Our data suggest that the development of herd management and multi-species livestock economies in the 11th-9th millennium BP drove speciation and host adaptation of this zoonotic pathogen.

Genetic diversity of the largest African mole-rat genus, Bathyergus. one, two or four species?

Recent advances in sequencing technology and phylogenetic methods allow us to solve puzzling taxonomic questions using detailed analyses of genetic diversity of populations and gene flow between them. The genus of solitary-living dune mole-rat, Bathyergus, is quite unique among six genera of African mole-rats. The animals are by far the largest and the only scratch digging mole-rat genus possessing a skull less adapted to digging, grooved upper incisors, and more surface locomotor activity. Most authors recognize two species of dune mole-rats, B. suillus and B. janetta, but according to others, the genus is monotypic. In addition, recent molecular studies have revealed cryptic genetic diversity and suggested the existence of up to four species. In our study, we used mitochondrial and genome-wide nuclear data collected throughout the distribution of the genus to investigate the number of species. In agreement with previous studies, we found Bathyergus to be differentiated into several distinct lineages, but we also found evidence for a degree of gene flow between some of them. Furthermore, we confirmed that B. janetta is nested within B. suillus, making the latter paraphyletic and we documented an instance of local mitochondrial introgression between these two nominal species. Phylogeographic structure of the genus was found to be very shallow. Although traditionally dated to the Miocene, we found the first split within the genus to be much younger estimated to 0.82 Ma before present. Genealogical distinctiveness of some lineages was very low, and the coancestry matrix showed extensive sharing of closely related haplotypes throughout the genus. Accordingly, Infomap clustering on the matrix showed all populations to form a single cluster. Overall, our study tends to support the existence of only one species of Bathyergus namely, B. suillus. Environmental niche modelling confirmed its dependence on sandy soils and the preference for soils with relatively high carbon content. Bayesian skyline plots indicate recent population decline in the janetta lineage, probably related to global environmental change.

Naming the other cousin: A new goldie barb (Cyprinidae: Smiliogastrininae) from the northeast escarpment in South Africa, with proposed taxonomic rearrangement of the goldie barb group in southern Africa

A growing body of evidence indicates that the global diversity of freshwater fishes has not been fully documented. Studies of freshwater fishes that were previously thought to be morphologically variable have revealed the existence of deeply divergent lineages, with many distinct species. In southern Africa a number of Enteromius species exhibit either exceedingly wide or divided distribution patterns that should be rare for freshwater fishes with limited dispersal opportunities between river systems. One such species is the sidespot barb, Enteromius neefi. As currently defined, E. neefi has a disjunct distribution that is divided between rivers in the northeast escarpment in South Africa and Eswatini, and tributaries of the Upper Zambezi in Zambia and southern Congo in the Democratic Republic of Congo, with a large geographic gap between these two populations. With the use of molecular and morphological methods, the level of divergence between the two populations was examined, and a new species was described from the Steelpoort River in the Limpopo River system of South Africa. Findings from this study provide further evidence for a number of taxonomic problems within the goldie barbs of southern Africa, and some taxonomic rearrangements are proposed for this group.

Spatial and temporal transmission dynamics of respiratory syncytial virus in New Zealand before and after the COVID-19 pandemic

Human respiratory syncytial virus (RSV) is a major cause of acute respiratory infection. In 2020, RSV was effectively eliminated from the community in New Zealand due to non-pharmaceutical interventions (NPI) used to control the spread of COVID-19. However, in April 2021, following a brief quarantine-free travel agreement with Australia, there was a large-scale nationwide outbreak of RSV that led to reported cases more than five times higher, and hospitalisations more than three times higher, than the typical seasonal pattern. In this study, we generated 1,471 viral genomes of both RSV-A and RSV-B sampled between 2015 and 2022 from across New Zealand. Using a phylodynamics approach, we used these data to better understand RSV transmission patterns in New Zealand prior to 2020, and how RSV became re-established in the community following the relaxation of COVID-19 restrictions. We found that in 2021, there was a large epidemic of RSV in New Zealand that affected a broader age group range compared to the usual pattern of RSV infections. This epidemic was due to an increase in RSV importations, leading to several large genomic clusters of both RSV-A ON1 and RSV-B BA9 genotypes in New Zealand. However, while a number of importations were detected, there was also a major reduction in RSV genetic diversity compared to pre-pandemic seasonal outbreaks. These genomic clusters were temporally associated with the increase of migration in 2021 due to quarantine-free travel from Australia at the time. The closest genetic relatives to the New Zealand RSV genomes, when sampled, were viral genomes sampled in Australia during a large, off-season summer outbreak several months prior, rather than cryptic lineages that were sustained but not detected in New Zealand. These data reveal the impact of NPI used during the COVID-19 pandemic on other respiratory infections and highlight the important insights that can be gained from viral genomes.

Whole genome phylogenomics helps to resolve the phylogenetic position of the Zygothrica genus group (Diptera, Drosophilidae) and the causes of previous incongruences

Incomplete Lineage Sorting (ILS) and introgression are among the two main factors causing incongruence between gene and species trees. Advances in phylogenomic studies have allowed us to overcome most of these issues, providing reliable phylogenetic hypotheses while revealing the underlying evolutionary scenario. Across the last century, many incongruent phylogenetic reconstructions were recovered for Drosophilidae, employing a limited sampling of genetic markers or species. In these studies, the monophyly and the phylogenetic positioning of the Zygothrica genus group stood out as one of the most controversial questions. Thus, here, we addressed these issues using a phylogenomic approach, while accessing the influence of ILS and introgressions on the diversification of these species and addressing the spatio-temporal scenario associated with their evolution. For this task, the genomes of nine specimens from six Neotropical species belonging to the Zygothrica genus group were sequenced and evaluated in a phylogenetic framework encompassing other 39 species of Drosophilidae. Nucleotide and amino acid sequences recovered for a set of 2,534 single-copy genes by BUSCO were employed to reconstruct maximum likelihood (ML) concatenated and multi-species coalescent (MSC) trees. Likelihood mapping, quartet sampling, and reticulation tests were employed to infer the level and causes of incongruence. Lastly, a penalized-likelihood molecular clock strategy with fossil calibrations was performed to infer divergence times. Taken together, our results recovered the subdivision of Drosophila into six different lineages, one of which clusters species of the Zygothrica genus group (except for H. duncani). The divergence of this lineage was dated to Oligocene ∼ 31 Mya and seems to have occurred in the same timeframe as other key diversification within Drosophila. According to the concatenated and MSC strategies, this lineage is sister to the clade joining Drosophila (Siphlodora) with the Hawaiian Drosophila and Scaptomyza. Likelihood mapping, quartet sampling, reticulation reconstructions as well as introgression tests revealed that this lineage was the target of several hybridization events involving the ancestors of different Drosophila lineages. Thus, our results generally show introgression as a major source of previous incongruence. Nevertheless, the similar diversification times recovered for several of the Neotropical Drosophila lineages also support the scenario of multiple and simultaneous diversifications taking place at the base of Drosophilidae phylogeny, at least in the Neotropics.

Early mutational signatures and transmissibility of SARS-CoV-2 Gamma and Lambda variants in Chile

Genomic surveillance (GS) programmes were crucial in identifying and quantifying the mutating patterns of SARS-CoV-2 during the COVID-19 pandemic. In this work, we develop a Bayesian framework to quantify the relative transmissibility of different variants tailored for regions with limited GS. We use it to study the relative transmissibility of SARS-CoV-2 variants in Chile. Among the 3443 SARS-CoV-2 genomes collected between January and June 2021, where sampling was designed to be representative, the Gamma (P.1), Lambda (C.37), Alpha (B.1.1.7), B.1.1.348, and B.1.1 lineages were predominant. We found that Lambda and Gamma variants' reproduction numbers were 5% (95% CI: [1%, 14%]) and 16% (95% CI: [11%, 21%]) larger than Alpha's, respectively. Besides, we observed a systematic mutation enrichment in the Spike gene for all circulating variants, which strongly correlated with variants' transmissibility during the studied period (r = 0.93, p-value = 0.025). We also characterised the mutational signatures of local samples and their evolution over time and with the progress of vaccination, comparing them with those of samples collected in other regions worldwide. Altogether, our work provides a reliable method for quantifying variant transmissibility under subsampling and emphasises the importance of continuous genomic surveillance.

Characterizing the Phan Rang Sheep: A First Look at the Y Chromosome, Mitochondrial DNA, and Morphometrics

The Phan Rang sheep, considered the sole indigenous breed of Vietnam, are primarily concentrated in the two central provinces of Ninh Thuan and Binh Thuan, with Ninh Thuan accounting for more than 90% of the country's sheep population. These provinces are known for their high temperatures and frequent droughts. The long-standing presence of the Phan Rang sheep in these regions suggests their potential resilience to heat stress-a trait of increasing interest in the face of global climate change. Despite the breed's significance, a critical knowledge gap hinders conservation and breeding programs. To address this, our study employed a two-pronged approach. First, we collected body conformational data to aid in breed identification. Second, we analyzed mitochondrial DNA (D-loop) and Y chromosome markers (SRY and SRYM18) to elucidate the maternal and paternal lineages. Among the 68 Phan Rang sheep analyzed for their D-loop, 19 belonged to mitochondrial haplogroup A, while 49 belonged to haplogroup B. The haplogroups can be subdivided into 16 unique haplotypes. All 19 rams surveyed for their paternal lineages belonged to haplotypes H5 and H6. These findings strongly support the hypothesis of dual origins for the Phan Rang sheep. This study presents the first genetic data for the Phan Rang breed, providing crucial insights for future research and conservation efforts.

Phylogenomic perspectives on speciation and reproductive isolation in a North American biodiversity hotspot: an example using California sages (Salvia subgenus Audibertia: Lamiaceae)

BACKGROUND AND AIMS

The California Floristic Province (CA-FP) is the most species-rich region of North America north of Mexico. One of several proposed hypotheses explaining the exceptional diversity of the region is that the CA-FP harbours myriad recently diverged lineages with nascent reproductive barriers. Salvia subgenus Audibertia is a conspicuous element of the CA-FP, with multiple sympatric and compatible species.

METHODS

Using 305 nuclear loci and both organellar genomes, we reconstruct species trees, examine genomic discordance, conduct divergence-time estimation, and analyse contemporaneous patterns of gene flow and mechanical reproductive isolation.

KEY RESULTS

Despite strong genomic discordance, an underlying bifurcating tree is supported. Organellar genomes capture additional introgression events not detected in the nuclear genome. Most interfertility is found within clades, indicating that reproductive barriers arise with increasing genetic divergence. Species are generally not mechanically isolated, suggesting that it is unlikely to be the primary factor leading to reproductive isolation.

CONCLUSIONS

Rapid, recent speciation with some interspecific gene flow in conjunction with the onset of a Mediterranean-like climate is the underlying cause of extant diversity in Salvia subgenus Audibertia. Speciation has largely not been facilitated by gene flow. Its signal in the nuclear genome seems to mostly be erased by backcrossing, but organellar genomes each capture different instances of historical gene flow, probably characteristic of many CA-FP lineages. Mechanical reproductive isolation appears to be only part of a mosaic of factors limiting gene flow.