Phylogeography of the sand dune ant Mycetophylax simplex along the Brazilian Atlantic Forest coast: remarkably low mtDNA diversity and shallow population structure

Coastal dunes as drivers of genetic differentiation in the honeypot ant Myrmecocystus baja (Formicidae: Formicinae)

Coastal dunes are unique habitats, threatened by human activities. In biogeographical terms, coastal dunes are habitat islands, being discrete and distinct patches of similar habitat among themselves, separated from each other by a different type of habitat. Furthermore, coastal dunes harbor endemic species, adapted to living solely in the habitats found on specific dune systems. For example, the honeypot ant Myrmecocystus baja is endemic and restricted to coastal dunes of Mexico's Baja California Pacific coast. This ecological and biogeographical scenario led to the questions whether their geographical isolation is reflected in their genetic diversity and structuring, and how their demographic history is related with the formation of the dune system habitats. To answer these questions, population genetic, isolation-with-migration, and phylogeographical analyses were carried out, based on mitochondrial and five nuclear intronic markers. Minimal gene flow was detected only between two of the dune systems sampled; otherwise, the M. baja populations were found to be isolated and genetically structured, and their divergence generally pre-dated the modern-day dune systems. It is therefore highly likely that these ants were already present in paleodunes and that each of the populations was established from founder populations as the dunes formed. These findings highlight the importance of coastal dunes for species such as the honeypot ant from Baja California, in promoting genetic differentiation.

De novo assembly of the complete mitochondrial genome of Mycetophylax simplex Emery, 1888 through organelle targeting revels no substantial expansion of gene spacers, but rather some slightly shorter genes

Mitochondria play a key role in cell biology and have their own genome, residing in a highly oxidative environment that induces faster changes than the nuclear genome. Because of this, mitochondrial markers have been exploited to reconstruct phylogenetic and phylogeographic relationships in studies of adaptation and molecular evolution. In this study, we determined the complete mitogenome of the fungus-farming ant Mycetophylax simplex (Hymenoptera, Formicidae) and conducted a comparative analysis among 29 myrmicine ant mitogenomes. Mycetophylax simplex is an endemic ant that inhabits sand dunes along the southern Atlantic coast. Specifically, the species occur in the ecosystem known as "restinga", within the Atlantic Forest biome. Due to habitat degradation, land use and decline of restinga habitats, the species is considered locally extinct in extremely urban beaches and is listed as vulnerable on the Brazilian Red List (ICMBio). We employed a mitochondrion-targeting approach to obtain the complete mitogenome through high-throughput DNA sequencing technology. This method allowed us to determine the mitogenome with high performance, coverage and low cost. The circular mitogenome has a length of 16,367 base pairs enclosing 37 genes (13 protein-coding genes, 22 tRNAs and 2 rRNAs) along with one control region (CR). All the protein-coding genes begin with a typical ATN codon and end with the canonical stop codons. All tRNAs formed the fully paired acceptor stems and fold into the typical cloverleaf-shaped secondary structures. The gene order is consistent with the shared Myrmicinae structure, and the A + T content of the majority strand is 81.51%. Long intergenic spacers were not found but some gene are slightly shorter. The phylogenetic relationships based on concatenated nucleotide and amino acid sequences of the 13 protein-coding genes, using Maximum Likelihood and Bayesian Inference methods, indicated that mitogenome sequences were useful in resolving higher-level relationship within Formicidae.

The influence of climate oscillations and geological events on population differentiation of Camponotus japonicus in the Chinese mainland

Camponotus japonicus (Hymenoptera: Formicidae) is an omnivorous social insect which builds sizable colonies in sparse woodlands or cropland and spreads across multiple climatic zones in the Chinese mainland. This study aims to reveal the role of climate changes and geological events in driving the genetic structure of social insect populations by investigating the phylogenetics and historical demography of C. japonicus in the Chinese mainland. Phylogenetic analyses were conducted based on the mitochondria DNA dataset using MrBayes and IQ-TREE. We constructed a haplotype network, calculated analyses of molecular variance, estimated the divergence time, and reconstructed the maximum clade credibility tree. Mismatch distribution and Bayesian skyline plots were used to infer historical population fluctuations. Additionally, ecological niche modeling was employed to predict the potential distribution of the species during the present, mid-holocene, and last glacial maximum periods in the Chinese mainland. The phylogenetic tree and median-joining network analyses support the presence of four distinct lineages in C. japonicus. These lineages exhibit significant genetic differentiation and limited gene flow. The divergence among the four lineages began in the early Pleistocene, approximately 1.41 million years ago (Ma). Subsequently, the central lineage diverged from both the northern and southern lineages around 1.16 Ma, while the northern and southern lineages diverged from each other at approximately 1.07 Ma. Population expansion was observed in the southern, central, and northern lineages prior to the last glacial maximum, while the Yunnan-Sichuan lineage experienced a slight increase in population size in more recent times. The predicted distribution of the species corresponds well with the actual distribution. Furthermore, the current suitable habitat areas in northern Xinjiang, southern Tibet, and the southeast coastal regions have significantly decreased compared to the last glacial maximum and the mid-holocene periods. Our results suggest that climate oscillations and geological events play an important role in driving genetic patterns and differentiation of C. japonicus. Mountain barriers isolate populations from each other, hinder the flow of genes, and effectively prevent the spread of this species. But at the same time, it also formed refugia at low altitudes areas such as Qinling-Bashan Mountains and Yanshan-Taihang Mountains and provide suitable habitats during glaciation. This study provides a good model for understanding how complex climate changes and geological events affect population genetic differentiation of social insects in the Chinese mainland.

Male-biased dispersal in a fungus-gardening ant symbiosis

For nearly all organisms, dispersal is a fundamental life-history trait that can shape their ecology and evolution. Variation in dispersal capabilities within a species exists and can influence population genetic structure and ecological interactions. In fungus-gardening (attine) ants, co-dispersal of ants and mutualistic fungi is crucial to the success of this obligate symbiosis. Female-biased dispersal (and gene flow) may be favored in attines because virgin queens carry the responsibility of dispersing the fungi, but a paucity of research has made this conclusion difficult. Here, we investigate dispersal of the fungus-gardening ant Trachymyrmex septentrionalis using a combination of maternally (mitochondrial DNA) and biparentally inherited (microsatellites) markers. We found three distinct, spatially isolated mitochondrial DNA haplotypes; two were found in the Florida panhandle and the other in the Florida peninsula. In contrast, biparental markers illustrated significant gene flow across this region and minimal spatial structure. The differential patterns uncovered from mitochondrial DNA and microsatellite markers suggest that most long-distance ant dispersal is male-biased and that females (and concomitantly the fungus) have more limited dispersal capabilities. Consequently, the limited female dispersal is likely an important bottleneck for the fungal symbiont. This bottleneck could slow fungal genetic diversification, which has significant implications for both ant hosts and fungal symbionts regarding population genetics, species distributions, adaptive responses to environmental change, and coevolutionary patterns.

Development, characterization, and cross-amplification of polymorphic microsatellite markers for North American Trachymyrmex and Mycetomoellerius ants

Objective

The objective of this study is to develop and identify polymorphic microsatellite markers for fungus-gardening (attine) ants in the genus Trachymyrmex sensu lato. These ants are important ecosystem engineers and have been a model group for understanding complex symbiotic systems, but very little is understood about the intraspecific genetic patterns across most North American attine species. These microsatellite markers will help to better study intraspecific population genetic structure, gene flow, mating habits, and phylogeographic patterns in these species and potentially other congeners.

Results

Using next-generation sequencing techniques, we identified 17 and 12 polymorphic microsatellite markers from T. septentrionalis and Mycetomoellerius (formerly Trachymyrmex) turrifex, respectively, and assessed the genetic diversity of each marker. We also analyzed the cross-amplification success of the T. septentrionalis markers in two other closely related Trachymyrmex species, and identified 10 and 12 polymorphic markers for T. arizonensis and T. pomonae, respectively.

Karyotype and putative chromosomal inversion suggested by integration of cytogenetic and molecular data of the fungus-farming ant Mycetomoellerius iheringi Emery, 1888

Comparative cytogenetic analyses are being increasingly used to collect information on species evolution, for example, diversification of closely related lineages and identification of morphologically indistinguishable species or lineages. Here, we have described the karyotype of the fungus-farming ant Mycetomoellerius iheringi Emery, 1888 and investigated its evolutionary relationships on the basis of molecular and cytogenetic data. The M. iheringi karyotype consists of 2n = 20 chromosomes (2K = 18M + 2SM). We also demonstrated that this species has the classical insect TTAGG telomere organization. Phylogenetic reconstruction showed that M. iheringi is phylogenetically closer to M. cirratus Mayhé-Nunes & Brandão, 2005 and M. kempfi Fowler, 1982. We compared M. iheringi with other congeneric species such as M. holmgreni Wheeler, 1925 and inferred that M. iheringi probably underwent a major pericentric inversion in one of its largest chromosomes, making it submetacentric. We discussed our results in the light of the phylogenetic relationships and chromosomal evolution.

Chromosomal dynamics in space and time: evolutionary history of Mycetophylax ants across past climatic changes in the Brazilian Atlantic coast

Fungus-farming ants of the genus Mycetophylax exhibit intra and interspecific chromosome variability, which makes them suitable for testing hypotheses about possible chromosomal rearrangements that endure lineage diversification. We combined cytogenetic and molecular data from Mycetophylax populations from coastal environments to trace the evolutionary history of the clade in light of chromosomal changes under a historical and geographic context. Our cytogenetic analyses revealed chromosomal differences within and among species. M. morschi exhibited three distinct karyotypes and considerable variability in the localization of 45S rDNA clusters. The molecular phylogeny was congruent with our cytogenetic findings. Biogeographical and divergence time dating analyses estimated that the most recent common ancestor of Mycetophylax would have originated at about 30 Ma in an area including the Amazon and Southern Grasslands, and several dispersion and vicariance events may have occurred before the colonization of the Brazilian Atlantic coast. Diversification of the psammophilous Mycetophylax first took place in the Middle Miocene (ca. 18-10 Ma) in the South Atlantic coast, while "M. morschi" lineages diversified during the Pliocene-Pleistocene transition (ca. 3-2 Ma) through founder-event dispersal for the Northern coastal regions. Psammophilous Mycetophylax diversification fits into the major global climatic events that have had a direct impact on the changes in sea level as well as deep ecological impact throughout South America. We assume therefore that putative chromosomal rearrangements correlated with increased ecological stress during the past climatic transitions could have intensified and/or accompanied the divergence of the psammophilous Mycetophylax. We further reiterate that "M. morschi" comprises a complex of at least three well-defined lineages, and we emphasize the role of this integrative approach for the identification and delimitation of evolutionary lineages.

Dual influence of terrestrial and marine historical processes on the phylogeography of the Brazilian intertidal red alga Gracilaria caudata

In this study, we explored how past terrestrial and marine climate changes have interacted to shape the phylogeographic patterns of the intertidal red seaweed Gracilaria caudata, an economically important species exploited for agar production in the Brazilian north-east. Seven sites were sampled along the north-east tropical and south-east sub-tropical Brazilian coast. The genetic diversity and structure of G. caudata was inferred using a combination of mitochondrial (COI and cox2-3), chloroplast (rbcL) and 15 nuclear microsatellite markers. A remarkable congruence between nuclear, mitochondrial and chloroplast data revealed clear separation between the north-east (from 03° S to 08° S) and the south-east (from 20° S to 23° S) coast of Brazil. These two clades differ in their demographic histories, with signatures of recent demographic expansions in the north-east and divergent populations in the south-east, suggesting the maintenance of several refugia during the last glacial maximum due to sea-level rise and fall. The Bahia region (around 12° S) occupies an intermediate position between both clades. Microsatellites and mtDNA markers showed additional levels of genetic structure within each sampled site located south of Bahia. The separation between the two main groups in G. caudata is likely recent, probably occurring during the Quaternary glacial cycles. The genetic breaks are concordant with (i) those separating terrestrial refugia, (ii) major river outflows and (iii) frontiers between tropical and subtropical regions. Taken together with previously published eco-physiological studies that showed differences in the physiological performance of the strains from distinct locations, these results suggest that the divergent clades in G. caudata correspond to distinct ecotypes in the process of incipient speciation and thus should be considered for the management policy of this commercially important species.

Hybridization and geographic distribution shapes the spatial genetic structure of two co-occurring orchid species

Multiple ecological and life-history traits shape the fine-scale spatial genetic structure (FSGS) of a given population. The occurrence in core versus peripheral populations, levels of outcrossing, pollen and seed dispersal, and hybridization are important biological properties that influence the kinship of individuals within populations. We examined spatial genetic structure within 15 populations of Epidendrum fulgens and E. puniceoluteum distributed along a linear gradient of Brazilian coastal vegetation, including both allopatric and sympatric populations where the two orchid species hybridize. We analyzed 581 mapped specimens using nine simple sequence repeat loci, aiming to investigate how geographic distribution and hybridization shape within-population FSGS. A significant increase in FSGS was found towards peripheral populations, compared to core populations. Analysis of short-distance and long-distance components of FSGS identified biparental inbreeding and higher levels of FSGS at peripheral populations, when compared to core populations. In contrast, the relatively high density of reproductive adults in core populations potentially leads to highly overlapping seed and pollen movement, decreasing FSGS. Hybridization was an important factor shaping within-population spatial genetic structure at sympatric sites, decreasing the FSGS observed in parental species. Our results indicate that different ecological forces act in concert to create a gradient of FSGS along species distribution ranges, shaped by extensive levels of intraspecific and interspecific gene exchange.

Neither barriers nor refugia explain genetic structure in a major biogeographic break: phylogeography of praying mantises in the Brazilian Atlantic Forest

The Atlantic Forest is one of the world's top biodiversity hotspots, but the diversification processes of its biota are still poorly known, with competing models attributing dominant roles to either Quaternary climatic changes or geographic barriers. Many studies identify the Doce river as a major phylogeographic break, but the reasons for this phenomenon are highly debated. Here we test the predictions of the refugial and barrier models for a common species of praying mantis, Miobantia fuscata, focusing in the areas immediately south and north of the Doce river. Our analyses show high intraspecific genetic diversity, deep coalescence times and no evidence for recent population expansion. Phylogeographic structure is inconsistent with a refugial hypothesis. Significant gene flow between northern and southern populations also conflicts with a strong role for geographic barriers. This study highlights the need for considering invertebrate taxa to infer recent landscape changes, and points towards a more complex picture of genetic diversification in the Atlantic Forest.

The xeric side of the Brazilian Atlantic Forest: The forces shaping phylogeographic structure of cacti

In order to investigate biogeographic influences on xeric biota in the Brazilian Atlantic Forest (BAF), a biodiversity hotspot, we used a monophyletic group including three cactus taxa as a model to perform a phylogeographic study: Cereus fernambucensis subsp. fernambucensis, C. fernambucensis subsp. sericifer, and C. insularis. These cacti are allopatric and grow in xeric habitats along BAF, including isolated granite and gneiss rock outcrops (Inselbergs), sand dune vegetation (Restinga forest), and the rocky shore of an oceanic archipelago (islands of Fernando de Noronha). The nucleotide information from nuclear gene phytochrome C and plastid intergenic spacer trnS-trnG was used to perform different approaches and statistical analyses, comprising population structure, demographic changes, phylogenetic relationships, and biogeographic reconstruction in both spatial and temporal scales. We recovered four allopatric population groups with highly supported branches in the phylogenetic tree with divergence initiated in the middle Pleistocene: southern distribution of C. fernambucensis subsp. fernambucensis, northern distribution of C. fernambucensis subsp. fernambucensis together with C. insularis, southern distribution of C. fernambucensis subsp. sericifer, and northern distribution of C. fernambucensis subsp. sericifer. Further, the results suggest that genetic diversity of population groups was strongly shaped by an initial colonization event from south to north followed by fragmentation. The phylogenetic pattern found for C. insularis is plausible with peripatric speciation in the archipelago of Fernando de Noronha. To explain the phylogeographic patterns, the putative effects of both climatic and sea level changes as well as neotectonic activity during the Pleistocene are discussed.

The roles of barriers, refugia, and chromosomal clines underlying diversification in Atlantic Forest social wasps

Phylogeographic studies have sought to explain the genetic imprints of historical climatic changes and geographic barriers within the Brazilian Atlantic Forest (AF) biota, and consequently two processes of diversification (refugia and barriers) have been proposed. Additionally, there is evidence that eustatic changes influenced the biogeographic history of the AF. Here we evaluate these contrasting diversification processes using two AF social wasp species – the mid-montane Synoeca cyanea and the lowland Synoeca aff. septentrionalis. We analyzed several sources of data including multilocus DNA sequence, climatic niche models and chromosomal features. We find support for idiosyncratic phylogeographic patterns between these wasps, involving different levels of population structure and genetic diversity, contrary suitable climatic conditions during the last glaciation, and contrasting historical movements along the AF. Our data indicate that neotectonics and refugia played distinct roles in shaping the genetic structure of these wasps. However, we argue that eustatic changes influenced the demographic expansion but not population structure in AF biota. Notably, these wasps exhibited chromosomal clines, involving chromosome number and decreasing of GC content, latitudinally oriented along the AF. Together, these results reinforce the need to consider individual organismal histories and indicate that barriers and refugia are significant factors in understanding AF evolution.

Integrating Paleodistribution Models and Phylogeography in the Grass-Cutting Ant Acromyrmex striatus (Hymenoptera: Formicidae) in Southern Lowlands of South America

Past climate changes often have influenced the present distribution and intraspecific genetic diversity of organisms. The objective of this study was to investigate the phylogeography and historical demography of populations of Acromyrmex striatus (Roger, 1863), a leaf-cutting ant species restricted to the open plains of South America. Additionally, we modeled the distribution of this species to predict its contemporary and historic habitat. From the partial sequences of the mitochondrial gene cytochrome oxidase I of 128 A. striatus workers from 38 locations we estimated genetic diversity and inferred historical demography, divergence time, and population structure. The potential distribution areas of A. striatus for current and quaternary weather conditions were modeled using the maximum entropy algorithm. We identified a total of 58 haplotypes, divided into five main haplogroups. The analysis of molecular variance (AMOVA) revealed that the largest proportion of genetic variation is found among the groups of populations. Paleodistribution models suggest that the potential habitat of A. striatus may have decreased during the Last Interglacial Period (LIG) and expanded during the Last Maximum Glacial (LGM). Overall, the past potential distribution recovered by the model comprises the current potential distribution of the species. The general structuring pattern observed was consistent with isolation by distance, suggesting a balance between gene flow and drift. Analysis of historical demography showed that populations of A. striatus had remained constant throughout its evolutionary history. Although fluctuations in the area of their potential historic habitat occurred during quaternary climate changes, populations of A. striatus are strongly structured geographically. However, explicit barriers to gene flow have not been identified. These findings closely match those in Mycetophylax simplex, another ant species that in some areas occurs in sympatry with A. striatus. Ecophysiological traits of this species and isolation by distance may together have shaped the phylogeographic pattern.