The biogeography of mitochondrial and nuclear discordance in animals

Mitochondrial Involvement in Vertebrate Speciation? The Case of Mito-nuclear Genetic Divergence in Chameleons

Compatibility between the nuclear (nDNA) and mitochondrial (mtDNA) genomes is important for organismal health. However, its significance for major evolutionary processes such as speciation is unclear, especially in vertebrates. We previously identified a sharp mtDNA-specific sequence divergence between morphologically indistinguishable chameleon populations (Chamaeleo chamaeleon recticrista) across an ancient Israeli marine barrier (Jezreel Valley). Because mtDNA introgression and gender-based dispersal were ruled out, we hypothesized that mtDNA spatial division was maintained by mito-nuclear functional compensation. Here, we studied RNA-seq generated from each of ten chameleons representing the north and south populations and identified candidate nonsynonymous substitutions (NSSs) matching the mtDNA spatial distribution. The most prominent NSS occurred in 14 nDNA-encoded mitochondrial proteins. Increased chameleon sample size (N = 70) confirmed the geographic differentiation in POLRMT, NDUFA5, ACO1, LYRM4, MARS2, and ACAD9. Structural and functionality evaluation of these NSSs revealed high functionality. Mathematical modeling suggested that this mito-nuclear spatial divergence is consistent with hybrid breakdown. We conclude that our presented evidence and mathematical model underline mito-nuclear interactions as a likely role player in incipient speciation in vertebrates.

Molecular and Morphological Differentiation of Common Dolphins (Delphinus sp.) in the Southwestern Atlantic: Testing the Two Species Hypothesis in Sympatry

The taxonomy of common dolphins (Delphinus sp.) has always been controversial, with over twenty described species since the original description of the type species of the genus (Delphinus delphis Linnaeus, 1758). Two species and four subspecies are currently accepted, but recent molecular data have challenged this view. In this study we investigated the molecular taxonomy of common dolphins through analyses of cytochrome b sequences of 297 individuals from most of their distribution. We included 37 novel sequences from the Southwestern Atlantic Ocean, a region where the short- and long-beaked morphotypes occur in sympatry, but which had not been well sampled before. Skulls of individuals from the Southwestern Atlantic were measured to test the validity of the rostral index as a diagnostic character and confirmed the presence of the two morphotypes in our genetic sample. Our genetic results show that all common dolphins in the Atlantic Ocean belong to a single species, Delphinus delphis. According to genetic data, the species Delphinus capensis is invalid. Long-beaked common dolphins from the Northeastern Pacific Ocean may constitute a different species. Our conclusions prompt the need for revision of currently accepted common dolphin species and subspecies and of Delphinus delphis distribution.

Molecular data indicate a cryptic species in Neotoma albigula (Cricetidae: Neotominae) from northwestern México

DNA sequences from the mitochondrial cytochrome- b gene were obtained from 41 specimens of Neotoma albigula from the southwestern United States and northwestern México. Phylogenetic analyses depicted that samples of N. a. melanura from southern Sonora and northern Sinaloa formed a clade separate from representatives of the other sampled subspecies of N. albigula ( albigula , laplataensis , mearnsi , seri , sheldoni , and venusta ). Genetic distances detected between these clades (7.41%) approached divergence levels reported for other sister species of woodrats and indicated that N. a. melanura is presumably a cryptic and genetically differentiated species relative to N . albigula . Analyses of DNA sequences from a nuclear gene (intron 2 of the alcohol dehydrogenase gene, Adh 1-I2) indicated that samples of N. a. melanura formed a separate, monophyletic clade relative to the remainder of N. albigula . Further, habitat, geographic distinctions, and morphological differences were apparent between members of the 2 clades. Together, those data support the elevation of N. a. melanura to species status.

Secuencias de ADN del gen mitocondrial citocromo- b se obtuvieron de 41 especímenes de Neotoma albigula provenientes del suroeste de los Estados Unidos y del noroeste de México. Análisis filogenéticos revelaron que las muestras de N. a. melanura del sureste de Sonora y del norte de Sinaloa formaron un clado separado de especímenes que representan a las otras subspecies de N. albigula (albigula, laplataensis, mearnsi, seri, sheldoni , y venusta ). Las distancias genéticas detectadas entre estos clados (7.41%) se aproximan a los valores de divergencia reportados para otras especies de ratas magueyeras, indicando que N. a. melanura es presumiblemente una especie críptica y relativamente diferenciada genéticamente de N. albigula . Los análisis de secuencias de ADN de un gen nuclear (intron 2 del gen alcohol deshidrogenasa, Adh -I2) indicaron que las muestras de N. a. melanura formaron un clado separado relativamente a especímenes de N. albigula . Adicionalmente, hábitat, distinciones geográficas, y diferencias morfológicas fueron aparentes entre miembros de los dos clados. Todos estos datos soportan la elevación de N. a. melanura al estatus de especie.

Genome sequencing of the extinct Eurasian wild aurochs, Bos primigenius, illuminates the phylogeography and evolution of cattle

Background

Domestication of the now-extinct wild aurochs, Bos primigenius, gave rise to the two major domestic extant cattle taxa, B. taurus and B. indicus. While previous genetic studies have shed some light on the evolutionary relationships between European aurochs and modern cattle, important questions remain unanswered, including the phylogenetic status of aurochs, whether gene flow from aurochs into early domestic populations occurred, and which genomic regions were subject to selection processes during and after domestication. Here, we address these questions using whole-genome sequencing data generated from an approximately 6,750-year-old British aurochs bone and genome sequence data from 81 additional cattle plus genome-wide single nucleotide polymorphism data from a diverse panel of 1,225 modern animals.

Results

Phylogenomic analyses place the aurochs as a distinct outgroup to the domestic B. taurus lineage, supporting the predominant Near Eastern origin of European cattle. Conversely, traditional British and Irish breeds share more genetic variants with this aurochs specimen than other European populations, supporting localized gene flow from aurochs into the ancestors of modern British and Irish cattle, perhaps through purposeful restocking by early herders in Britain. Finally, the functions of genes showing evidence for positive selection in B. taurus are enriched for neurobiology, growth, metabolism and immunobiology, suggesting that these biological processes have been important in the domestication of cattle.

Conclusions

This work provides important new information regarding the origins and functional evolution of modern cattle, revealing that the interface between early European domestic populations and wild aurochs was significantly more complex than previously thought.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0790-2) contains supplementary material, which is available to authorized users.

Mitochondrial introgression via ancient hybridization, and systematics of the Australian endemic pygopodid gecko genus Delma

Of the more than 1500 species of geckos found across six continents, few remain as unfamiliar as the pygopodids - Family Pygopodidae (Gray, 1845). These gekkotans are limited to Australia (44 species) and New Guinea (2 species), but have diverged extensively into the most ecologically diverse limbless radiation save Serpentes. Current phylogenetic understanding of the family has relied almost exclusively on two works, which have produced and synthesized an immense amount of morphological, geographical, and molecular data. However, current interspecific relationships within the largest genus Delma Gray 1831 are based chiefly upon data from two mitochondrial loci (16s, ND2). Here, we reevaluate the interspecific relationships within the genus Delma using two mitochondrial and four nuclear loci (RAG1, MXRA5, MOS, DYNLL1), and identify points of strong conflict between nuclear and mitochondrial genomic data. We address mito-nuclear discordance, and remedy this conflict by recognizing several points of mitochondrial introgression as the result of ancient hybridization events. Owing to the legacy value and intraspecific informativeness, we suggest the continued use of ND2 as a phylogenetic marker. Results identify strong support for species groups, but relationships among these clades, and the placement of several enigmatic taxa remain uncertain. We suggest a more careful review of Delma australis and the 'northwest Australia' clade. Accurately assessing and addressing species richness and relationships within this endemic Australian Gekkotan genus is relevant for understanding patterns of squamate speciation across the region.

Speciation processes in putative island endemic sister bat species: false impressions from mitochondrial DNA and microsatellite data

Cases of geographically restricted co-occurring sister taxa are rare and may point to potential divergence with gene flow. The two bat species Murina gracilis and Murina recondita are both endemic to Taiwan and are putative sister species. To test for nonallopatric divergence and gene flow in these taxa, we generated sequences using Sanger and next-generation sequencing, and combined these with microsatellite data for coalescent-based analyses. MtDNA phylogenies supported the reciprocally monophyletic sister relationship between M. gracilis and M. recondita; however, clustering of microsatellite genotypes revealed several cases of species admixture suggesting possible introgression. Sequencing of microsatellite flanking regions revealed that admixture signatures stemmed from microsatellite allele homoplasy rather than recent introgressive hybridization, and also uncovered an unexpected sister relationship between M. recondita and the continental species Murina eleryi, to the exclusion of M. gracilis. To dissect the basis of these conflicts between ncDNA and mtDNA, we analysed sequences from 10 anonymous ncDNA loci with *beast and isolation-with-migration and found two distinct clades of M. eleryi, one of which was sister to M. recondita. We conclude that Taiwan was colonized by the ancestor of M. gracilis first, followed by the ancestor of M. recondita after a period of allopatric divergence. After colonization, the mitochondrial genome of M. recondita was replaced by that of the resident M. gracilis. This study illustrates how apparent signatures of sympatric divergence can arise from complex histories of allopatric divergence, colonization and hybridization, thus highlighting the need for rigorous analyses to distinguish between such scenarios.

Introgression of mitochondrial DNA among lineages in a hybridogenetic ant

We report a remarkable pattern of incongruence between nuclear and mitochondrial variations in a social insect, the desert ant Cataglyphis hispanica. This species reproduces by social hybridogenesis. In all populations, two distinct genetic lineages coexist; non-reproductive workers develop from hybrid crosses between the lineages, whereas reproductive offspring (males and new queens) are typically produced asexually by parthenogenesis. Genetic analyses based on nuclear markers revealed that the two lineages remain highly differentiated despite constant hybridization for worker production. Here, we show that, in contrast with nuclear DNA, mitochondrial DNA (mtDNA) does not recover the two lineages as monophyletic. Rather, mitochondrial haplotypes cluster according to their geographical origin. We argue that this cytonuclear incongruence stems from introgression of mtDNA among lineages, and review the mechanisms likely to explain this pattern under social hybridogenesis.

Speciation Progress: A Case Study on the Bushcricket Poecilimon veluchianus

Different mechanisms such as selection or genetic drift permitted e.g. by geographical isolation can lead to differentiation of populations and could cause subsequent speciation. The two subspecies of Poecilimon veluchianus, a bushcricket endemic to central Greece, show a parapatric distribution and are partially reproductively isolated. Therefore, P. veluchianus is suitable to investigate an ongoing speciation process. We based our analysis on sequences of the internal transcribed spacer (ITS) and the mitochondrial control region (CR). The population genetic analysis based on the nuclear marker ITS revealed a barrier to gene flow within the range of Poecilimon veluchianus, which corresponds well to the described subspecies. In contrast to the results based on the nuclear ITS marker, the mitochondrial CR marker does not clearly support the separation into two subspecies with restricted gene flow and a clear contact zone. Furthermore, we could identify isolation by distance (IBD) as one important mechanism responsible for the observed genetic structure (based on the ITS marker). The population genetic analysis based on the nuclear marker ITS also suggests the existence of hybrids in the wild. Furthermore, the simultaneous lack of strong prezygotic barriers and the presence of postzygotic mating barriers, observed in previous laboratory experiments, suggest that a secondary contact after an allopatric phase is more likely than parapatric speciation.

Deeply divergent sympatric mitochondrial lineages of the earthworm Lumbricus rubellus are not reproductively isolated

Background

The accurate delimitation of species is essential to numerous areas of biological research. An unbiased assessment of the diversity, including the cryptic diversity, is of particular importance for the below ground fauna, a major component of global biodiversity. On the British Isles, the epigeic earthworm Lumbricus rubellus, which is a sentinel species in soil ecotoxicology, consists of two cryptic taxa that are differentiated in both the nuclear and the mitochondrial (mtDNA) genomes. Recently, several deeply divergent mtDNA lineages were detected in mainland Europe, but whether these earthworms also constitute cryptic species remains unclear. This information is important from an evolutionary perspective, but it is also essential for the interpretation and the design of ecotoxicological projects. In this study, we used genome-wide RADseq data to assess the reproductive isolation of the divergent mitochondrial lineages of L. rubellus that occur in sympatry in multiple localities in Central Europe.

Results

We identified five divergent (up to 16 % net p-distance) mitochondrial lineages of L. rubellus in sympatry. Because the clustering of the RADseq data was according to the population of origin and not the mtDNA lineage, reproductive isolation among the mtDNA lineages was not likely. Although each population contained multiple mtDNA lineages, subdivisions within the populations were not observed for the nuclear genome. The lack of fixed differences and sharing of the overwhelming majority of nuclear polymorphisms between localities, indicated that the populations did not constitute allopatric species. The nucleotide diversity within the populations was high, 0.7–0.8 %.

Conclusions

The deeply divergent mtDNA sympatric lineages of L. rubellus in Central Europe were not reproductively isolated groups. The earthworm L. rubellus, which is represented by several mtDNA lineages in continental Europe, apparently is a single highly polymorphic species rather than a complex of several cryptic species. This study demonstrated the critical importance of the use of multilocus nuclear data for the unbiased assessment of cryptic diversity and for the delimitation of species in soil invertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0488-9) contains supplementary material, which is available to authorized users.

Contrasting genetic structure between mitochondrial and nuclear markers in the dengue fever mosquito from Rio de Janeiro: implications for vector control

Dengue is the most prevalent global arboviral disease that affects over 300 million people every year. Brazil has the highest number of dengue cases in the world, with the most severe epidemics in the city of Rio de Janeiro (Rio). The effective control of dengue is critically dependent on the knowledge of population genetic structuring in the primary dengue vector, the mosquito Aedes aegypti. We analyzed mitochondrial and nuclear genomewide single nucleotide polymorphism markers generated via Restriction-site Associated DNA sequencing, as well as traditional microsatellite markers in Ae. aegypti from Rio. We found four divergent mitochondrial lineages and a strong spatial structuring of mitochondrial variation, in contrast to the overall nuclear homogeneity across Rio. Despite a low overall differentiation in the nuclear genome, we detected strong spatial structure for variation in over 20 genes that have a significantly altered expression in response to insecticides, xenobiotics, and pathogens, including the novel biocontrol agent Wolbachia. Our results indicate that high genetic diversity, spatially unconstrained admixing likely mediated by male dispersal, along with locally heterogeneous genetic variation that could affect insecticide resistance and mosquito vectorial capacity, set limits to the effectiveness of measures to control dengue fever in Rio.

Mechanisms Underlying Mammalian Hybrid Sterility in Two Feline Interspecies Models

The phenomenon of male sterility in interspecies hybrids has been observed for over a century, however, few genes influencing this recurrent phenotype have been identified. Genetic investigations have been primarily limited to a small number of model organisms, thus limiting our understanding of the underlying molecular basis of this well-documented "rule of speciation." We utilized two interspecies hybrid cat breeds in a genome-wide association study employing the Illumina 63 K single-nucleotide polymorphism array. Collectively, we identified eight autosomal genes/gene regions underlying associations with hybrid male sterility (HMS) involved in the function of the blood-testis barrier, gamete structural development, and transcriptional regulation. We also identified several candidate hybrid sterility regions on the X chromosome, with most residing in close proximity to complex duplicated regions. Differential gene expression analyses revealed significant chromosome-wide upregulation of X chromosome transcripts in testes of sterile hybrids, which were enriched for genes involved in chromatin regulation of gene expression. Our expression results parallel those reported in Mus hybrids, supporting the "Large X-Effect" in mammalian HMS and the potential epigenetic basis for this phenomenon. These results support the value of the interspecies feline model as a powerful tool for comparison to rodent models of HMS, demonstrating unique aspects and potential commonalities that underpin mammalian reproductive isolation.

In-solution hybridization for mammalian mitogenome enrichment: pros, cons and challenges associated with multiplexing degraded DNA

Here, we present a set of RNA-based probes for whole mitochondrial genome in-solution enrichment, targeting a diversity of mammalian mitogenomes. This probes set was designed from seven mammalian orders and tested to determine the utility for enriching degraded DNA. We generated 63 mitogenomes representing five orders and 22 genera of mammals that yielded varying coverage ranging from 0 to >5400X. Based on a threshold of 70% mitogenome recovery and at least 10× average coverage, 32 individuals or 51% of samples were considered successful. The estimated sequence divergence of samples from the probe sequences used to construct the array ranged up to nearly 20%. Sample type was more predictive of mitogenome recovery than sample age. The proportion of reads from each individual in multiplexed enrichments was highly skewed, with each pool having one sample that yielded a majority of the reads. Recovery across each mitochondrial gene varied with most samples exhibiting regions with gaps or ambiguous sites. We estimated the ability of the probes to capture mitogenomes from a diversity of mammalian taxa not included here by performing a clustering analysis of published sequences for 100 taxa representing most mammalian orders. Our study demonstrates that a general array can be cost and time effective when there is a need to screen a modest number of individuals from a variety of taxa. We also address the practical concerns for using such a tool, with regard to pooling samples, generating high quality mitogenomes and detail a pipeline to remove chimeric molecules.

Interactions between nuclear genes and a foreign mitochondrial genome in the redbelly dace Chrosomus eos

Given the coevolution process occurring between nuclear and mitochondrial genomes, the effects of introgressive hybridization remain puzzling. In this study, we take advantage of the natural co-occurrence of two biotypes bearing a similar nuclear genome (Chrosomus eos) but harbouring mitochondria from different species (wild type: C. eos; cybrids: Chrosomus neogaeus) to determine the extent of phenotype changes linked to divergence in the mitochondrial genome. Changes were assessed through differences in gene expression, enzymatic activity, proteomic and swimming activity. Our data demonstrate that complex IV activity was significantly higher in cybrids compared to wild type. This difference could result from one variable amino acid on the COX3 mitochondrial subunit and/or from a tremendous change in the proteome. We also show that cybrids present a higher swimming performance than wild type. Ultimately, our results demonstrate that the absence of coevolution for a period of almost ten million years between nuclear and mitochondrial genomes does not appear to be necessarily deleterious but could even have beneficial effects. Indeed, the capture of foreign mitochondria could be an efficient way to circumvent the selection process of genomic coevolution, allowing the rapid accumulation of new mutations in C. eos cybrids.

Timeframe of speciation inferred from secondary contact zones in the European tree frog radiation (Hyla arborea group)

Background

Hybridization between incipient species is expected to become progressively limited as their genetic divergence increases and reproductive isolation proceeds. Amphibian radiations and their secondary contact zones are useful models to infer the timeframes of speciation, but empirical data from natural systems remains extremely scarce. Here we follow this approach in the European radiation of tree frogs (Hyla arborea group). We investigated a natural hybrid zone between two lineages (Hyla arborea and Hyla orientalis) of Mio-Pliocene divergence (~5 My) for comparison with other hybrid systems from this group.

Results

We found concordant geographic distributions of nuclear and mitochondrial gene pools, and replicated narrow transitions (~30 km) across two independent transects, indicating an advanced state of reproductive isolation and potential local barriers to dispersal. This result parallels the situation between H. arborea and H. intermedia, which share the same amount of divergence with H. orientalis. In contrast, younger lineages show much stronger admixture at secondary contacts.

Conclusions

Our findings corroborate the negative relationship between hybridizability and divergence time in European tree frogs, where 5 My are necessary to achieve almost complete reproductive isolation. Speciation seems to progress homogeneously in this radiation, and might thus be driven by gradual genome-wide changes rather than single speciation genes. However, the timescale differs greatly from that of other well-studied amphibians. General assumptions on the time necessary for speciation based on evidence from unrelated taxa may thus be unreliable. In contrast, comparative hybrid zone analyses within single radiations such as our case study are useful to appreciate the advance of speciation in space and time.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0385-2) contains supplementary material, which is available to authorized users.

Biased introgression of mitochondrial and nuclear genes: a comparison of diploid and haplodiploid systems

Hybridization between recently diverged species, even if infrequent, can lead to the introgression of genes from one species into another. The rates of mitochondrial and nuclear introgression often differ, with some taxa showing biases for mitochondrial introgression and others for nuclear introgression. Several hypotheses exist to explain such biases, including adaptive introgression, sex differences in dispersal rates, sex-specific prezygotic isolation and sex-specific fitness of hybrids (e.g. Haldane's rule). We derive a simple population genetic model that permits an analysis of sex-specific demographic and fitness parameters and measures the relative rates of mitochondrial and nuclear introgression between hybridizing pairs. We do this separately for diploid and haplodiploid species. For diploid taxa, we recover results consistent with previous hypotheses: an excess of one sex among the hybridizing migrants or sex-specific prezygotic isolation causes a bias for one type of marker or the other; when Haldane's rule is obeyed, we find a mitochondrial bias in XY systems and a nuclear bias in ZW systems. For haplodiploid taxa, the model reveals that owing to their unique transmission genetics, they are seemingly assured of strong mitochondrial biases in introgression rates, unlike diploid taxa, where the relative fitness of male and female hybrids can tip the bias in either direction. This heretofore overlooked aspect of hybridization in haplodiploids provides what is perhaps the most likely explanation for differential introgression of mitochondrial and nuclear markers and raises concerns about the use of mitochondrial DNA barcodes for species delimitation in these taxa.

Biogeography and evolution of a widespread Central American lizard species complex: Norops humilis, (Squamata: Dactyloidae)

Background

Caribbean anole lizards (Dactyloidae) have frequently been used as models to study questions regarding biogeography and adaptive radiations, but the evolutionary history of Central American anoles (particularly those of the genus Norops) has not been well studied. Previous work has hypothesized a north-to-south dispersal pattern of Central American Norops, but no studies have examined dispersal within any Norops lineages. Here we test two major hypotheses for the dispersal of the N. humilis/quaggulus complex (defined herein, forming a subset within Savage and Guyer’s N. humilis group).

Results

Specimens of the N. humilis group were collected in Central America, from eastern Mexico to the Canal Zone of Panama. Major nodes were dated for comparison to the geologic history of Central America, and ancestral ranges were estimated for the N. humilis/quaggulus complex to test hypothesized dispersal patterns. These lineages displayed a northward dispersal pattern. We also demonstrate that the N. humilis/quaggulus complex consists of a series of highly differentiated mitochondrial lineages, with more conserved nuclear evolution. The paraphyly of the N. humilis species group is confirmed. A spatial analysis of molecular variance suggests that current populations are genetically distinct from one another, with limited mitochondrial gene flow occurring among sites.

Conclusions

The observed south-to-north colonization route within the Norops humilis/quaggulus complex represents the first evidence of a Norops lineage colonizing in a south-to-north pattern, (opposite to the previously held hypothesis for mainland Norops). One previously described taxon (N. quaggulus) was nested within N. humilis, demonstrating the paraphyly of this species; while our analyses also reject the monophyly of the Norops humilis species group (sensu Savage and Guyer), with N. tropidonotus, N. uniformis, and N. marsupialis being distantly related to/highly divergent from the N. humilis/quaggulus complex. Our work sheds light on mainland anole biogeography and past dispersal events, providing a pattern to test against other groups of mainland anoles.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0391-4) contains supplementary material, which is available to authorized users.

Genome-wide analysis of signatures of selection in populations of African honey bees (Apis mellifera) using new web-based tools

Background

With the development of inexpensive, high-throughput sequencing technologies, it has become feasible to examine questions related to population genetics and molecular evolution of non-model species in their ecological contexts on a genome-wide scale. Here, we employed a newly developed suite of integrated, web-based programs to examine population dynamics and signatures of selection across the genome using several well-established tests, including F_ST, pN/pS, and McDonald-Kreitman. We applied these techniques to study populations of honey bees (Apis mellifera) in East Africa. In Kenya, there are several described A. mellifera subspecies, which are thought to be localized to distinct ecological regions.

Results

We performed whole genome sequencing of 11 worker honey bees from apiaries distributed throughout Kenya and identified 3.6 million putative single-nucleotide polymorphisms. The dense coverage allowed us to apply several computational procedures to study population structure and the evolutionary relationships among the populations, and to detect signs of adaptive evolution across the genome. While there is considerable gene flow among the sampled populations, there are clear distinctions between populations from the northern desert region and those from the temperate, savannah region. We identified several genes showing population genetic patterns consistent with positive selection within African bee populations, and between these populations and European A. mellifera or Asian Apis florea.

Conclusions

These results lay the groundwork for future studies of adaptive ecological evolution in honey bees, and demonstrate the use of new, freely available web-based tools and workflows (http://usegalaxy.org/r/kenyanbee) that can be applied to any model system with genomic information.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1712-0) contains supplementary material, which is available to authorized users.

Regional differentiation and extensive hybridization between mitochondrial clades of the Southern Ocean giant sea spider Colossendeis megalonyx

Assessing the enormous diversity of Southern Ocean benthic species and their evolutionary histories is a central task in the era of global climate change. Based on mitochondrial markers, it was recently suggested that the circumpolar giant sea spider Colossendeis megalonyx comprises a complex of at least six cryptic species with mostly small and non-overlapping distribution ranges. Here, we expand the sampling to include over 500 mitochondrial COI sequences of specimens from around the Antarctic. Using multiple species delimitation approaches, the number of distinct mitochondrial OTUs increased from six to 15-20 with our larger dataset. In contrast to earlier studies, many of these clades show almost circumpolar distributions. Additionally, analysis of the nuclear internal transcribed spacer region for a subset of these specimens showed incongruence between nuclear and mitochondrial results. These mito-nuclear discordances suggest that several of the divergent mitochondrial lineages can hybridize and should not be interpreted as cryptic species. Our results suggest survival of C. megalonyx during Pleistocene glaciations in multiple refugia, some of them probably located on the Antarctic shelf, and emphasize the importance of multi-gene datasets to detect the presence of cryptic species, rather than their inference based on mitochondrial data alone.

Persistence, isolation and diversification of a naturally fragmented species in local refugia: the case of Hydromantes strinatii

The study of the European plethodontid salamander Hydromantes strinatii using allozyme and mitochondrial markers showed a strong geographical genetic structure. This was likely the outcome of different evolutionary mechanisms leaving their signature despite the effects of the genetic drift due to the low population size typical of this species. Two highly divergent clades were identified in the eastern and central-western part of the range, with further geographic sub-structure. Nuclear and mitochondrial markers substantially recovered the same population groups but were conflicting in reconstructing their relationships. This apparent incongruence highlighted the action of different mechanisms such as secondary contacts and incomplete lineage sorting in originating the observed genetic variation. The troglophilic habit of this species provided the opportunity to show the importance of caves as local refugia in maintaining the genetic diversity through the persistence of local populations. Accordingly, high nucleotide and haplotype diversity, strong geographic genetic structuring and lack of expansion were evidenced. This signature was found in the populations from the Ligurian and Maritime Alps, in agreement with the complex orography and paleoclimatic history of this Mediterranean hotspot.

Causations of phylogeographic barrier of some rocky shore species along the Chinese coastline

BACKGROUND

Substrate, ocean current and freshwater discharge are recognized as important factors that control the larval dispersal and recruitment of intertidal species. Life history traits of individual species will determine the differential responses to these physical factors, and hence resulting in contrasting phylogeography across the same biogeographic barrier. To determine how these factors affect genetic structure of rocky shore species along the China coast, a comparative phylogeographic study of four intertidal and subtidal species was conducted using mitochondrial and nuclear DNA by combining new sequences from Siphonaria japonica with previously published sequences from three species (Cellana toreuma, Sargassum horneri and Atrina pectinata).

RESULTS

Analysis of molecular variance and pairwise ΦST revealed significant genetic differences between the Yellow Sea (YS) and the other two marginal seas (East China Sea, ECS and South China Sea, SCS) for rocky-shore species (S. japonica, C. toreuma, S. horneri), but not for muddy-shore species Atrina pectinata. Demographic history analysis proved that the population size of all these four species were persistent though the Last Glacial Maximum (LGM, ~20 ka BP). Migration analysis revealed that gene flow differentiated northward and southward migration for these four species. However, the inferred direction of gene flow using alternatively mitochondrial or nuclear markers was contradictory in S. japonica.

CONCLUSIONS

It is concluded that there is a phylogeographical break at the Yangtze River estuary for the rocky shore species and the causation of the barrier is mainly due to the unsuitable substratum and freshwater discharge. All four intertidal and subtidal species appear to have persisted through the LGM in China, indicating the lower impact of LGM on intertidal and subtidal species than generally anticipated. The imbalanced gene flow between YS and ESCS groups for these four species could be explained by historical refugia. The discordance between mitochondrial and nuclear markers in the MIGRATE analysis of S. japonica prove the importance of employing multi-locus data in biogeographic study. Climate change, land reclamation and dam construction, which are changing substrate and hydrological conditions around Yangtze River estuary, will consequently affect the biogeographic pattern of intertidal species.

Morphological and genomic comparisons of Hawaiian and Japanese Black-footed Albatrosses (Phoebastria nigripes) using double digest RADseq: implications for conservation

Evaluating the genetic and demographic independence of populations of threatened species is important for determining appropriate conservation measures, but different technologies can yield different conclusions. Despite multiple studies, the taxonomic status and extent of gene flow between the main breeding populations of Black-footed Albatross (Phoebastria nigripes), a Near-Threatened philopatric seabird, are still controversial. Here, we employ double digest RADseq to quantify the extent of genomewide divergence and gene flow in this species. Our genomewide data set of 9760 loci containing 3455 single nucleotide polymorphisms yielded estimates of genetic diversity and gene flow that were generally robust across seven different filtering and sampling protocols and suggest a low level of genomic variation (θ per site = ∼0.00002-0.00028), with estimates of effective population size (N e = ∼500-15 881) falling far below current census size. Genetic differentiation was small but detectable between Japan and Hawaii (F ST ≈ 0.038-0.049), with no F ST outliers. Additionally, using museum specimens, we found that effect sizes of morphological differences by sex or population rarely exceeded 4%. These patterns suggest that the Hawaiian and Japanese populations exhibit small but significant differences and should be considered separate management units, although the evolutionary and adaptive consequences of this differentiation remain to be identified.

Estimating the temporal and spatial extent of gene flow among sympatric lizard populations (genus Sceloporus) in the southern Mexican highlands

Interspecific gene flow is pervasive throughout the tree of life. Although detecting gene flow between populations has been facilitated by new analytical approaches, determining the timing and geography of hybridization has remained difficult, particularly for historical gene flow. A geographically explicit phylogenetic approach is needed to determine the overlap of ancestral populations. In this study, we performed population genetic analyses, species delimitation, simulations and a recently developed approach of species tree diffusion to infer the phylogeographic history, timing and geographic extent of gene flow in lizards of the Sceloporus spinosus group. The two species in this group, S. spinosus and S. horridus, are distributed in eastern and western portions of Mexico, respectively, but populations of these species are sympatric in the southern Mexican highlands. We generated data consisting of three mitochondrial genes and eight nuclear loci for 148 and 68 individuals, respectively. We delimited six lineages in this group, but found strong evidence of mito-nuclear discordance in sympatric populations of S. spinosus and S. horridus owing to mitochondrial introgression. We used coalescent simulations to differentiate ancestral gene flow from secondary contact, but found mixed support for these two models. Bayesian phylogeography indicated more than 60% range overlap between ancestral S. spinosus and S. horridus populations since the time of their divergence. Isolation-migration analyses, however, revealed near-zero levels of gene flow between these ancestral populations. Interpreting results from both simulations and empirical data indicate that despite a long history of sympatry among these two species, gene flow in this group has only recently occurred.

Mitonuclear Ecology

Eukaryotes were born of a chimeric union between two prokaryotes--the progenitors of the mitochondrial and nuclear genomes. Early in eukaryote evolution, most mitochondrial genes were lost or transferred to the nucleus, but a core set of genes that code exclusively for products associated with the electron transport system remained in the mitochondrion. The products of these mitochondrial genes work in intimate association with the products of nuclear genes to enable oxidative phosphorylation and core energy production. The need for coadaptation, the challenge of cotransmission, and the possibility of genomic conflict between mitochondrial and nuclear genes have profound consequences for the ecology and evolution of eukaryotic life. An emerging interdisciplinary field that I call "mitonuclear ecology" is reassessing core concepts in evolutionary ecology including sexual reproduction, two sexes, sexual selection, adaptation, and speciation in light of the interactions of mitochondrial and nuclear genomes.

Phylogeny, hybridization, and life history evolution of Rhinogobius gobies in Japan, inferred from multiple nuclear gene sequences

Rhinogobius fishes (Gobiidae) are distributed widely in East and Southeast Asia, and represent the most species-rich group of freshwater gobies with diversified life histories (i.e., amphidromous, fluvial, and lentic). To reveal their phylogenetic relationships and life history evolution patterns, we sequenced six nuclear and three mitochondrial DNA (mtDNA) loci from 18 species, mainly from the mainland of Japan and the Ryukyu Archipelago. Our phylogenetic tree based on nuclear genes resolved three major clades, including several distinct subclades. The mtDNA and nuclear DNA phylogenies showed large discordance, which strongly suggested mitochondrial introgression through large-scale interspecific hybridization in these regions. On the basis of the molecular dating using geological data as calibration points, the hybridization occurred in the early to middle Pleistocene. Reconstruction of the ancestral states of life history traits based on nuclear DNA phylogeny suggests that the evolutionary change from amphidromous to freshwater life, accompanied by egg size change, occurred independently in at least three lineages. One of these lineages showed two life history alterations, i.e., from amphidromous (small egg) to fluvial (large egg) to lentic (small egg). Although more inclusive analysis using species outside Japan should be further conducted, the present results suggest the importance of the life history evolution associated with high adaptability to freshwater environments in the remarkable species diversification in this group. Such life history divergences may have contributed to the development of reproductive isolation.

Genetic evaluation of the evolutionary distinctness of a federally endangered butterfly, Lange's Metalmark

Background

The Mormon Metalmark (Apodemia mormo) species complex occurs as isolated and phenotypically variable colonies in dryland areas across western North America. Lange’s Metalmark, A. m. langei, one of the 17 subspecies taxonomically recognized in the complex, is federally listed under the U.S. Endangered Species Act of 1973. Metalmark taxa have traditionally been described based on phenotypic and ecological characteristics, and it is unknown how well this nomenclature reflects their genetic and evolutionary distinctiveness. Genetic variation in six microsatellite loci and mitochondrial cytochrome oxidase subunit I sequence was used to assess the population structure of the A. mormo species complex across 69 localities, and to evaluate A. m. langei’s qualifications as an Evolutionarily Significant Unit.

Results

We discovered substantial genetic divergence within the species complex, especially across the Continental Divide, with population genetic structure corresponding more closely with geographic proximity and local isolation than with taxonomic divisions originally based on wing color and pattern characters. Lange’s Metalmark was as genetically divergent as several other locally isolated populations in California, and even the unique phenotype that warranted subspecific and conservation status is reminiscent of the morphological variation found in some other populations.

Conclusions

This study is the first genetic treatment of the A. mormo complex across western North America and potentially provides a foundation for reassessing the taxonomy of the group. Furthermore, these results illustrate the utility of molecular markers to aid in demarcation of biological units below the species level. From a conservation point of view, Apodemia mormo langei’s diagnostic taxonomic characteristics may, by themselves, not support its evolutionary significance, which has implications for its formal listing as an Endangered Species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0354-9) contains supplementary material, which is available to authorized users.

Endosymbiosis and its implications for evolutionary theory

Historically, conceptualizations of symbiosis and endosymbiosis have been pitted against Darwinian or neo-Darwinian evolutionary theory. In more recent times, Lynn Margulis has argued vigorously along these lines. However, there are only shallow grounds for finding Darwinian concepts or population genetic theory incompatible with endosymbiosis. But is population genetics sufficiently explanatory of endosymbiosis and its role in evolution? Population genetics "follows" genes, is replication-centric, and is concerned with vertically consistent genetic lineages. It may also have explanatory limitations with regard to macroevolution. Even so, asking whether population genetics explains endosymbiosis may have the question the wrong way around. We should instead be asking how explanatory of evolution endosymbiosis is, and exactly which features of evolution it might be explaining. This paper will discuss how metabolic innovations associated with endosymbioses can drive evolution and thus provide an explanatory account of important episodes in the history of life. Metabolic explanations are both proximate and ultimate, in the same way genetic explanations are. Endosymbioses, therefore, point evolutionary biology toward an important dimension of evolutionary explanation.

Pleistocene niche stability and lineage diversification in the subtropical spider Araneus omnicolor (Araneidae)

The influence of Quaternary climate oscillations on the diversification of the South American fauna is being increasingly explored. However, most of these studies have focused on taxa that are endemic to tropical environments, and relatively few have treated organisms restricted to subtropical biomes. Here we used an integrative phylogeographical framework to investigate the effects of these climate events on the ecological niche and genetic patterns of the subtropical orb-weaver spider Araneus omnicolor (Araneidae). We analyzed the mitochondrial (Cytochrome Oxidase I, COI) and nuclear (Internal Transcribed Subunit II, ITS2) DNA of 130 individuals throughout the species' range, and generated distribution models in three different climate scenarios [present, Last Glacial Maximum (LGM), and Last Interglacial Maximum (LIG)]. Additionally, we used an Approximate Bayesian Computation (ABC) approach to compare possible demographic scenarios and select the hypothesis that better explains the genetic patterns of A. omnicolor. We obtained high haplotype diversity but low nucleotide variation among sequences. The population structure and demographic analyses showed discrepancies between markers, suggesting male-biased dispersal in the species. The time-calibrated COI phylogenetic inference showed a recent diversification of lineages (Middle/Late Pleistocene), while the paleoclimate modeling indicated niche stability since ~120 Kya. The ABC results agreed with the niche models, supporting a panmictic population as the most likely historical scenario for the species. These results indicate that A. omnicolor experienced no niche or population reductions during the Late Pleistocene, despite the intense landscape modifications that occurred in the subtropical region, and that other factors beside LGM and LIG climate oscillations might have contributed to the demographic history of this species. This pattern may be related to the high dispersal ability and wide environmental tolerance of A. omnicolor, highlighting the need for more phylogeographical studies with invertebrates and other generalist taxa, in order to understand the effects of Quaternary climate changes on Neotropical biodiversity.

Limited connectivity and a phylogeographic break characterize populations of the pink anemonefish, Amphiprion perideraion, in the Indo-Malay Archipelago: inferences from a mitochondrial and microsatellite loci

To enhance the understanding of larval dispersal in marine organisms, species with a sedentary adult stage and a pelagic larval phase of known duration constitute ideal candidates, because inferences can be made about the role of larval dispersal in population connectivity. Members of the immensely diverse marine fauna of the Indo-Malay Archipelago are of particular importance in this respect, as biodiversity conservation is becoming a large concern in this region. In this study, the genetic population structure of the pink anemonefish, Amphiprion perideraion, is analyzed by applying 10 microsatellite loci as well as sequences of the mitochondrial control region to also allow for a direct comparison of marker-derived results. Both marker systems detected a strong overall genetic structure (ΦST = 0.096, P < 0.0001; mean D est = 0.17; F ST = 0.015, P < 0.0001) and best supported regional groupings (ΦCT = 0.199 P < 0.0001; F CT = 0.018, P < 0.001) that suggested a differentiation of the Java Sea population from the rest of the archipelago. Differentiation of a New Guinea group was confirmed by both markers, but disagreed over the affinity of populations from west New Guinea. Mitochondrial data suggest higher connectivity among populations with fewer signals of regional substructure than microsatellite data. Considering the homogenizing effect of only a few migrants per generation on genetic differentiation between populations, marker-specific results have important implications for conservation efforts concerning this and similar species.

New insights into samango monkey speciation in South Africa

The samango monkey is South Africa's only exclusively forest dwelling primate and represents the southernmost extent of the range of arboreal guenons in Africa. The main threats to South Africa's forests and thus to the samango are linked to increasing land-use pressure and increasing demands for forest resources, resulting in deforestation, degradation and further fragmentation of irreplaceable habitats. The species belongs to the highly polytypic Cercopithecus nictitans group which is sometimes divided into two species C. mitis and C. albogularis. The number of subspecies of C. albogularis is also under debate and is based only on differences in pelage colouration and thus far no genetic research has been undertaken on South African samango monkey populations. In this study we aim to further clarify the number of samango monkey subspecies, as well as their respective distributions in South Africa by combining molecular, morphometric and pelage data. Overall, our study provides the most comprehensive view to date into the taxonomic description of samango monkeys in South Africa. Our data supports the identification of three distinct genetic entities namely; C. a. labiatus, C. a. erythrarchus and C. a. schwarzi and argues for separate conservation management of the distinct genetic entities defined by this study.