Comparative phylogeography in the Atlantic forest and Brazilian savannas: pleistocene fluctuations and dispersal shape spatial patterns in two bumblebees

Phylogeographic and genetic insights into Sinonychia martensi: an endemic cave-dwelling harvestman in Beijing

Caves are one of the most exciting environments on earth, often considered an evolutionary laboratory due to the suite of convergent adaptive traits (troglomorphisms) of organisms inhabiting them. Sinonychia martensi Zhang & Derkarabetian, 2021, is the first and only Travunioidea species recorded in China and is endemic to Beijing, being known from multiple caves. However, nothing is known regarding its phylogeographic or evolutionary history. In this study, we assessed the species boundaries of S. martensi from nine caves using morphological and molecular methods to elucidate its phylogenetic position and genealogical relationships. We also investigated the genetic diversity, population genetic structure and demographic history of S. martensi to clarify the population-level relationships and make inferences about historical phylogeography. The results indicate that the species from different caves all belonged to S. martensi but represent different populations. These populations exhibit strong population structure and low genetic diversity. Cave populations may share a common ancestor and multiple independent invasions to different caves. The diversification within S. martensi was likely driven by climate change and subtropical evergreen broadleaf forests associated with the middle Miocene. This study highlights the need for further conservation efforts and exploration in Beijing caves.

Montane diversification as a mechanism of speciation in neotropical butterflies

The mountains in the Atlantic Forest domain are environments that harbor a high biodiversity, including species adapted to colder climates that were probably influenced by the climatic variations of the Pleistocene. To understand the phylogeographic pattern and assess the taxonomic boundaries between two sister montane species, a genomic study of the butterflies Actinote mantiqueira and A. alalia (Nymphalidae: Acraeini) was conducted. Analyses based on partial sequences of the mitochondrial gene COI (barcode region) failed to recover any phylogenetic or genetic structure discriminating the two species or sampling localities. However, single nucleotide polymorphisms gathered using Genotyping-by-Sequencing provided a strong isolation pattern in all analyses (genetic distance, phylogenetic hypothesis, clustering analyses, and F ST statistics) which is consistent with morphology, separating all individuals of A. alalia from all populations of A. mantiqueira. The three sampled mountain ranges where A. mantiqueira populations occur-Serra do Mar, Serra da Mantiqueira, and Poços de Caldas Plateau-were identified as three isolated clusters. Paleoclimate simulations indicate that both species' distributions changed according to climatic oscillations in the Pleistocene period, with the two species potentially occurring in areas of lower altitude during glacial periods when compared to the interglacial periods (as the present). Besides, a potential path between their distribution through the Serra do Mar Mountain range was inferred. Therefore, the Pleistocene climatic fluctuation had a significant impact on the speciation process between A. alalia and A. mantiqueira, which was brought on by isolation at different mountain summits during interglacial periods, as shown by the modeled historical distribution and the observed genetic structure.

The complete mitochondrial genome of Trigonisca nataliae (Hymenoptera, Apidae) assemblage reveals heteroplasmy in the control region

The evolution of mitochondrial genomes in the stingless bees is surprisingly dynamic, making them a model system to understand mitogenome structure, function, and evolution. Out of the seven mitogenomes available in this group, five exhibit atypical characteristics, including extreme rearrangements, rapid evolution and complete mitogenome duplication. To further explore the mitogenome diversity in these bees, we utilized isolated mtDNA and Illumina sequencing to assemble the complete mitogenome of Trigonisca nataliae, a species found in Northern Brazil. The mitogenome of T. nataliae was highly conserved in gene content and structure when compared to Melipona species but diverged in the control region (CR). Using PCR amplification, cloning and Sanger sequencing, six different CR haplotypes, varying in size and content, were recovery. These findings indicate that heteroplasmy, where different mitochondrial haplotypes coexist within individuals, occurs in T. nataliae. Consequently, we argue that heteroplasmy might indeed be a common phenomenon in bees that could be associated with variations in mitogenome size and challenges encountered during the assembly process.

Unparalleled mitochondrial heteroplasmy and Wolbachia co-infection in the non-model bee, Amphylaeus morosus

Mitochondrial heteroplasmy is the occurrence of more than one type of mitochondrial DNA within a single individual. Although generally reported to occur in a small subset of individuals within a species, there are some instances of widespread heteroplasmy across entire populations. Amphylaeus morosus is an Australian native bee species in the diverse and cosmopolitan bee family Colletidae. This species has an extensive geographical range along the eastern Australian coast, from southern Queensland to western Victoria, covering approximately 2,000 km. Seventy individuals were collected from five localities across this geographical range and sequenced using Sanger sequencing for the mitochondrial cytochrome c oxidase subunit I (COI) gene. These data indicate that every individual had the same consistent heteroplasmic sites but no other nucleotide variation, suggesting two conserved and widespread heteroplasmic mitogenomes. Ion Torrent shotgun sequencing revealed that heteroplasmy occurred across multiple mitochondrial protein-coding genes and is unlikely explained by transposition of mitochondrial genes into the nuclear genome (NUMTs). DNA sequence data also demonstrated a consistent co-infection of Wolbachia across the A. morosus distribution with every individual infected with both bacterial strains. Our data are consistent with the presence of two mitogenomes within all individuals examined in this species and suggest a major divergence from standard patterns of mitochondrial inheritance. Because the host's mitogenome and the Wolbachia genome are genetically linked through maternal inheritance, we propose three possible hypotheses that could explain maintenance of the widespread and conserved co-occurring bacterial and mitochondrial genomes in this species.

Phylogeography and genetic diversity of the widespread katydid Ducetia japonica (Thunberg, 1815) across China

Habitat fragmentation can lower migration rates and genetic connectivity among remaining populations of native species. Ducetia japonica is one of the most widespread katydids in China, but little is known about its genetic structure and phylogeographic distribution. We combined the five-prime region of cytochrome c oxidase subunit I (COI-5P), 11 newly developed microsatellite loci coupled with an ecological niche model (ENM) to examine the genetic diversity and population structure of D. japonica in China and beyond to Laos and Singapore. Both Bayesian inference (BI) and haplotype network methods revealed six mitochondrial COI-5P lineages. The distribution of COI-5P haplotypes may not demonstrate significant phylogeographic structure (N ST > G ST, p > .05). The STRUCTURE analysis based on microsatellite data also revealed six genetic clusters, but discordant with those obtained from COI-5P haplotypes. For both COI-5P and microsatellite data, Mantel tests revealed a significant positive correlation between geographic and genetic distances in mainland China. Bayesian skyline plot (BSP) analyses indicated that the population size of D. japonica's three major mitochondrial COI-5P lineages were seemingly not affected by last glacial maximum (LGM, 0.015-0.025 Mya). The ecological niche models showed that the current distribution of D. japonica was similar to the species' distribution during the LGM period and only slightly extended in northern China. Further phylogeographic studies based on more extensive sampling are needed to identify specific locations of glacial refugia in northern China.

North-south and climate-landscape-associated pattern of population structure for the Atlantic Forest White Morpho butterflies

Atlantic Forest White Morpho butterflies, currently classified as Morpho epistrophus and M. iphitus, are endemic to the Atlantic Forest, where they are widely distributed throughout heterogeneous environmental conditions. Studies with endemic butterflies allow to elucidate questions on both patterns of diversity distribution and current and past processes acting on insect groups in this biodiversity hotspot. In the present study, we characterized one mtDNA marker (COI sequences) and developed 11 polymorphic loci of microsatellite for 22 sampling locations distributed throughout the entire Atlantic Forest domain. We investigated both the taxonomic limits of taxa classified as White Morpho and the structure and distribution of the genetic diversity throughout their populations. Genetic markers and distribution data failed to identify species diversification, population structure, or isolation among subpopulations attributed to different taxa proposed for the White Morpho, suggesting that the current distinction between two species is unreasonable. The Bayesian coalescent tree based on COI sequences also failed to recover monophyletic clades for the putative species, and pointed instead to a north-south oriented pattern of genetic structure, with the northern clade coalescing later than the southern clade. Northern samples also showed more intragroup structure than southern samples based on mtDNA data. Clustering tests based on microsatellites indicated the existence of three genetic clusters, with turnover between the states of Paraná and São Paulo. The north-south pattern found for the White Morpho populations is showed for the first time to a endemic AF insect and coincides with the two different bioclimatic domains previously described for vertebrates and plants. Population structure observed for these butterflies is related to climate- and landscape-associated variables, mainly precipitation and elevation.

Connection, isolation and reconnection: Quaternary climatic oscillations and the Andes shaped the phylogeographical patterns of the Patagonian bee Centris cineraria (Apidae)

The joint effect of the Andes as a geographical barrier and the Quaternary glaciations as promoters of genetic divergence remains virtually unexplored in southern South America. To help fill this knowledge gap, in this study we investigated the demographic history of Centris cineraria, a solitary bee mainly distributed in Patagonia. We used mitochondrial and nuclear markers and performed phylogeographical and dating analyses, adjusted spatio-temporal diffusion and species distribution models, and used Approximate Bayesian Computation to identify likely historical demographic scenarios. Our results revealed that during glacial periods the Andes represented a barrier due to the extent of the ice-sheets and the occurrence of unsuitable habitats, while interglacials allowed for gene flow across the Andes. Secondary contact between previously isolated lineages was evident across at least two low-altitude Andean areas, the northern one being a putative glacial refugium. Our findings also suggest that C. cineraria has persisted in situ in four periglacial refugia located along a north–south transect, congruent with the maximum extent of the ice sheet during the Greatest Patagonian Glaciation. As the first phylogeographical study of Patagonian insects, our work reveals that the interaction between Quaternary climatic oscillations and the Andes as a barrier was the main driver of the spatial and demographic history of C. cineraria.

Mitochondrial DNA intra-individual variation in a bumblebee species: A challenge for evolutionary studies and molecular identification

Mitochondrial DNA (mtDNA) regions have been widely used as molecular markers in evolutionary studies and species identification. However, the presence of heteroplasmy and NUMTs may represent obstacles. Heteroplasmy is a state where an organism has different mitochondrial haplotypes. NUMTs are nuclear pseudogenes originating from mtDNA sequences transferred to nuclear DNA. Evidences of heteroplasmy were already verified in the bumblebee Bombus morio in an earlier study. The present work investigated in more detail the presence of intra-individual haplotypes variation in this species. Heteroplasmy was detected in individuals from all the ten sampled locations, with an average of six heteroplasmic haplotypes per individual. In addition, some of these heteroplasmic haplotypes were shared among individuals from different locations, suggesting the existence of stable heteroplasmy in B. morio. These results demonstrated that heteroplasmy is likely to affect inferences based on mtDNA analysis, especially in phylogenetic, phylogeographic and population genetics studies. In addition, NUMTs were also detected. These sequences showed divergence of 2.7% to 12% in relation to the mitochondrial haplotypes. These levels of divergence could mislead conclusions in evolutionary studies and affect species identification through DNA barcoding.

Climatic dynamics and topography control genetic variation in Atlantic Forest montane birds

Montane organisms responded to Quaternary climate change by tracking suitable habitat along elevational gradients. However, it is unclear whether these past climatic dynamics generated predictable patterns of genetic diversity in co-occurring montane taxa. To test if the genetic variation is associated with historical changes in the elevational distribution of montane habitats, we integrated paleoclimatic data and a model selection approach for testing the demographic history of five co-distributed bird species occurring in the southern Atlantic Forest sky islands. We found that changes in historical population sizes and current genetic diversity are attributable to habitat dynamics among time periods and the current elevational distribution of populations. Taxa with populations restricted to the more climatically dynamic southern mountain block (SMB) had, on average, a six-fold demographic expansion, whereas the populations from the northern mountain block (NMB) remained constant. In the current configuration of the southern Atlantic Forest montane habitats, populations in the SMB have more widespread elevational distributions, occur at lower elevations, and harbor higher levels of genetic diversity than NMB populations. Despite the apparent coupling of demographic and climatic oscillations, our data rejected simultaneous population structuring due to historical habitat fragmentation. Demographic modeling indicated that the species had different modes of differentiation, and varied in the timing of divergence and the degree of gene flow across mountain blocks. Our results suggest that the heterogeneous distribution of genetic variation in birds of the Atlantic Forest sky islands is associated with the interplay between topography and climate of distinct mountains, leading to predictable patterns of genetic diversity.

Fidelity of DNA polymerases in the detection of intraindividual variation of mitochondrial DNA

Here we investigated the consequences of PCR amplification errors in the identification of intraindividual mtDNA variation. The bumblebee Bombus morio was chosen as model for the COI gene amplification tests with two DNA polymerases (Taq and Q5) presenting different error rates. The amplifications using Taq resulted in a significant increase of singleton haplotypes per individual in comparison to Q5. The sequence characteristics indicated that Taq resulted haplotypes are mostly due to amplification errors. Studies focusing on intraindividual variability should address special attention to the DNA polymerase fidelity to avoid overestimation of heteroplasmic haplotypes.

Conserved numts mask a highly divergent mitochondrial-COI gene in a species complex of Australian stingless bees Tetragonula (Hymenoptera: Apidae)

Tetragonula carbonaria, Tetragonula davenporti, Tetragonula hockingsi and Tetragonula mellipes comprise a species complex of Australian stingless bee species known as the 'Carbonaria' group. The species are difficult to distinguish morphologically and the major species-defining characters relate to comb architecture and nest entrance ornamentation. The taxonomy of the group is further complicated by likely nuclear mitochondrial pseudogenes (numts) and inter-specific hybrids. Here we demonstrate the existence of COI numts and isolate and characterize the 'true' mt-COI gene in T. carbonaria and T. hockingsi. Numts were isolated from enriched-nuclear DNA extraction followed by PCR amplification and Sanger sequencing, and were recognized by the presence of deletions and/or premature stop codons in the translated sequences. The mt-COI sequences were obtained from NGS sequencing using purified mtDNA. In T. carbonaria, two numts (numt1 and numt2) were identified and a third (numt3) was identified in T. hockingsi. Numt2 and numt3 are similar (1.2% sequence divergence), indicating a recent common origin. The genetic distance between the mt-COI of the two Tetragonula species was higher than might be expected for closely related species, 16.5%, corroborating previous studies in which T. carbonaria and T. hockingsi were regarded as separate species. The three numts are more similar to the COI of other stingless bee species, including Australian Austroplebia australis and South American Melipona bicolor (81.7-83.9%) than to the mt-COI of their own species (70-71.4%). This is because the mt-COI of T. carbonaria and T. hockingsi differ greatly from other Meliponinae. Our findings explain some formerly puzzling aspects of Carbonaria biogeography, and misinterpreted amplifications.

Insights into the evolutionary dynamics of Neotropical biomes from the phylogeography and paleodistribution modeling of Bromelia balansae

PREMISE OF THE STUDY

Historical abiotic and biotic factors have strongly affected species diversification and speciation. Although pre-Pleistocene events have been linked to the divergence of several Neotropical organisms, studies have highlighted a more prominent role of Pleistocene climatic oscillations in shaping current patterns of genetic variation of plants.

METHODS

We performed phylogeographic analyses based on plastidial markers and modeled the current distribution and paleodistribution of Bromelia balansae (Bromeliaceae), an herbaceous species with a wide geographical distribution in South America, to infer the processes underlying its evolutionary history.

KEY RESULTS

Combined molecular and paleodistributional modeling analyses indicated retraction during the Last Glacial Maximum followed by interglacial expansion. Populations occurring in the semideciduous Atlantic Forest and the Cerrado formed two distinct genetic clusters, which have been historically or ecologically isolated since late Pliocene to early Pleistocene. Populations located in the transition zone had higher levels of genetic diversity, as expected by the long-term climatic stability in the region detected in our ecological niche models.

CONCLUSIONS

Our study adds important information on how herbaceous species have been affected by past climate in Central and Southeast Brazil, helping to disentangle the complex processes that have triggered the evolution of Neotropical biota.

Quaternary climate instability is correlated with patterns of population genetic variability in Bombus huntii

Climate oscillations have left a significant impact on the patterns of genetic diversity observed in numerous taxa. In this study, we examine the effect of Quaternary climate instability on population genetic variability of a bumble bee pollinator species, Bombus huntii in western North America. Pleistocene and contemporary B. huntii habitat suitability (HS) was estimated with an environmental niche model (ENM) by associating 1,035 locality records with 10 bioclimatic variables. To estimate genetic variability, we genotyped 380 individuals from 33 localities at 13 microsatellite loci. Bayesian inference was used to examine population structure with and without a priori specification of geographic locality. We compared isolation by distance (IBD) and isolation by resistance (IBR) models to examine population differentiation within and among the Bayesian inferred genetic clusters. Furthermore, we tested for the effect of environmental niche stability (ENS) on population genetic diversity with linear regression. As predicted, high-latitude B. huntii habitats exhibit low ENS when compared to low-latitude habitats. Two major genetic clusters of B. huntii inhabit western North America: (a) a north genetic cluster predominantly distributed north of 28°N and (b) a south genetic cluster distributed south of 28°N. In the south genetic cluser, both IBD and IBR models are significant. However, in the north genetic cluster, IBD is significant but not IBR. Furthermore, the IBR models suggest that low-latitude montane populations are surrounded by habitat with low HS, possibly limiting dispersal, and ultimately gene flow between populations. Finally, we detected high genetic diversity across populations in regions that have been climatically unstable since the last glacial maximum (LGM), and low genetic diversity across populations in regions that have been climatically stable since the LGM. Understanding how species have responded to climate change has the potential to inform management and conservation decisions of both ecological and economic concerns.

Environmental constraints on the compositional and phylogenetic beta-diversity of tropical forest snake assemblages

The ongoing biodiversity crisis increases the importance and urgency of studies addressing the role of environmental variation on the composition and evolutionary history of species assemblages, but especially the tropics and ectotherms remain understudied. In regions with rainy summers, coexistence of tropical ectothermic species may be determined by the partitioning of the climatic niche, as ectotherms can rely on water availability and thermoregulatory behaviour to buffer constraints along their climatic niche. Conversely, tropical ectotherms facing dry summers would have fewer opportunities to climatic niche partitioning and other processes rather than environmental filtering would mediate species coexistence. We used 218 snake assemblages to quantify the compositional (CBD) and phylogenetic (PBD) beta-diversity of snakes in the Atlantic Forest (AF) hotspot, South America. We identify two AF regions with distinct climatological regimes: dry summers in the northern-AF and rainy summers in the southern-AF. While accounting for the influence of multiscale spatial processes, we disentangle the relative contribution of thermal, water-related and topographic conditions in structuring the CBD and PBD of snake assemblages, and determine the extent in which snake assemblages under distinct climatological regimes are affected by environmental filtering. Thermal conditions best explain CBD and PBD of snakes for the whole AF, whereas water-related factors best explain the structure of snake assemblages within a same climatological regime. CBD and PBD patterns are similarly explained by spatial factors but snake assemblages facing dry summers are more affected by spatial processes operating at fine to intermediate spatial scale, whereas those assemblages in regions with rainy summers have a stronger signature of coarse-scale processes. As expected, environmental filtering plays a stronger role in southern-AF than northern-AF, and the synergism between thermal and water-related conditions is the key cause behind this difference. Differences in climatological regimes within the tropics affect processes mediating species coexistence. The influence of broad-scale gradients (e.g. temperature and precipitation) in structuring tropical ectothermic assemblages is greater in regions with rainy summers where climatic niche partitioning is more likely. Our findings highlight the potential stronger role of biotic interactions and neutral processes in structuring ectothermic assemblages facing changes towards warmer and dryer climates.