A multifaceted ecological assessment reveals the invasion of the freshwater red macroalga Montagnia macrospora (Batrachospermales,Rhodophyta) in Taiwan

Invasive freshwater macroalgae are rarely described. Montagnia macrospora is a freshwater red alga introduced from South America to East Asia via the global aquarium trade. The earliest occurrence record of this alga in Taiwan is dated 2005. To determine whether M. macrospora has become invasive in Taiwan and to understand the traits that facilitated its invasion, we took a multifaceted approach that combines examination of ecological background and population genetic analysis. Our island‐wide survey showed that M. macrospora is widespread in the field across Taiwan, where the climate greatly differs from that of South America, and can self‐sustain for nearly a decade. Our population genetic analysis revealed a lack of genetic diversity of M. macrospora in Taiwan, consistent with the hypothesis that the alga expanded through asexual reproduction. Moreover, during our long‐term ecological assessments and field surveys, we observed that M. macrospora is an ecological generalist that can survive in a wide range of temperature, pH, illumination, and nutrient enrichment. Taken together, our data suggest that M. macrospora has successfully invaded the freshwater ecosystems of Taiwan, likely due to its ability to disperse asexually and to grow under broad environmental conditions. We hope that our study brings attention to invasive freshwater algae, which have been overlooked in conservation planning and management.     

Reduced genetic diversity associated with the northern expansion of an amphibian species with high habitat specialization,Ascaphus truei,resolved using two types of genetic markers

Reconstruction of historical relationships between geographic regions within a species’ range can indicate dispersal patterns and help predict future responses to shifts in climate. Ascaphus truei (coastal tailed frog) is an indicator species of the health of forests and perennial streams in the Coastal and Cascade Mountains of the Pacific Northwest of North America. We used two genetic techniques—microsatellite and genotype‐by‐sequencing (GBS)—to compare the within‐region genetic diversity of populations near the northern extent of the species’ range (British Columbia, Canada) to two geographic regions in British Columbia and two in Washington, USA, moving toward the core of the range. Allelic richness and heterozygosity declined substantially as latitude increased. The northernmost region had the lowest mean expected heterozygosities for both techniques (microsatellite, M = 0.20, SE = 0.080; GBS, M = 0.025, SE = 0.0010) and the southernmost region had the highest (microsatellite, M = 0.88, SE = 0.054; GBS, M = 0.20, SE = 0.0029). The northernmost regions (NC and MC) clustered together in population structure models for both genetic techniques. Our discovery of reduced diversity may have important conservation and management implications for population connectivity and the response of A. truei to climate change.     

Population genomic analyses reveal hybridization and marked differences in genetic structure of Scurria limpet sister species with parapatric distributions across the South Eastern Pacific

The study of sister species that occur in parapatry around biogeographic transition zones can help understand the evolutionary processes that underlie the changes in species composition across biogeographic transition zones. The South Eastern Pacific (SEP) coast is a highly productive coastal system that exhibits a broad biogeographic transition zone around 30–35°S. Here, we present a comparative genome‐wide analysis of the sister species Scurria viridula and Scurria zebrina, that occur in parapatry and whose poleward and equatorward range edges intersect in the 30–35°S SEP biogeographic transition zone. We sampled 118 specimens sourced from nine sites from Tocopilla (22°S) to Chiloé (41°S) including one site where both species overlap and analyzed over 8000 biallelic single nucleotide polymorphisms. We found evidence of hybridization between these species in the contact zone and found significant but contrasting population structures for both species. Our results indicate that the genetic structure in S. viridula, which is currently expanding its range poleward, follows a simple isolation by distance model with no traces of natural selection (no evidence of outlier loci). In contrast, S. zebrina, which has  its equatorward range edge at the transition zone, displayed a pronounced genetic break approximately at 32–34°S, along a region of marked environmental heterogeneity in association with a semi‐permanent coastal upwelling regime. For S. zebrina we also found 43 outlier loci associated with this genetic break, with a significant proportion of them clustering in a single linkage group. This marked difference in the presence of outlier loci between species suggests that they could be responding differently to local environmental challenges found at their overlapping geographic range edges, thus providing important new insights about genomic changes around biogeographic transition zones in sister species and the forces that shape genetic diversity in intertidal marine species.     

EST‐SSR‐based landscape genetics of Pseudotaxus chienii,a tertiary relict conifer endemic to China

Pseudotaxus chienii, belonging to the monotypic genus Pseudotaxus (Taxaceae), is a relict conifer endemic to China. Its populations are usually small and patchily distributed, having a low capacity of natural regeneration. To gain a clearer understanding of how landscape variables affect the local adaptation of P. chienii, we applied EST‐SSR markers in conjunction with landscape genetics methods: (a) to examine the population genetic pattern and spatial genetic structure; (b) to perform genome scan and selection scan to identify outlier loci and the associated landscape variables; and (c) to model the ecological niche under climate change. As a result, P. chienii was found to have a moderate level of genetic variation and a high level of genetic differentiation. Its populations displayed a significant positive relationship between the genetic and geographical distance (i.e., “isolation by distance” pattern) and a strong fine‐scale spatial genetic structure within 2 km. A putatively adaptive locus EMS6 (functionally annotated to cellulose synthase A catalytic subunit 7) was identified, which was found significantly associated with soil Cu, K, and Pb content and the combined effects of temperature and precipitation. Moreover, P. chienii was predicted to experience significant range contractions in future climate change scenarios. Our results highlight the potential of specific soil metal content and climate variables as the driving force of adaptive genetic differentiation in P. chienii. The data would also be useful to develop a conservation action plan for P. chienii.     

Major biogeographic barriers in eastern Australia have shaped the population structure of widely distributed Eucalyptus moluccana and its putative subspecies

We have investigated the impact of recognized biogeographic barriers on genetic differentiation of grey box (Eucalyptus moluccana), a common and widespread tree species of the family Myrtaceae in eastern Australian woodlands, and its previously proposed four subspecies moluccana, pedicellata, queenslandica, and crassifolia. A range of phylogeographic analyses were conducted to examine the population genetic differentiation and subspecies genetic structure in E. moluccana in relation to biogeographic barriers. Slow evolving markers uncovering long term processes (chloroplast DNA) were used to generate a haplotype network and infer phylogeographic barriers. Additionally, fast evolving, hypervariable markers (microsatellites) were used to estimate demographic processes and genetic structure among five geographic regions (29 populations) across the entire distribution of E. moluccana. Morphological features of seedlings, such as leaf and stem traits, were assessed to evaluate population clusters and test differentiation of the putative subspecies. Haplotype network analysis revealed twenty chloroplast haplotypes with a main haplotype in a central position shared by individuals belonging to the regions containing the four putative subspecies. Microsatellite analysis detected the genetic structure between Queensland (QLD) and New South Wales (NSW) populations, consistent with the McPherson Range barrier, an east‐west spur of the Great Dividing Range. The substructure was detected within QLD and NSW in line with other barriers in eastern Australia. The morphological analyses supported differentiation between QLD and NSW populations, with no difference within QLD, yet some differentiation within NSW populations. Our molecular and morphological analyses provide evidence that several geographic barriers in eastern Australia, including the Burdekin Gap and the McPherson Range have contributed to the genetic structure of E. moluccana. Genetic differentiation among E. moluccana populations supports the recognition of some but not all the four previously proposed subspecies, with crassifolia being the most differentiated.     

The dispersal between Amazonia and Atlantic Forest during the Early Neogene revealed by the biogeography of the treefrog tribe Sphaenorhynchini (Anura,Hylidae)

The Amazonia and the Atlantic Forest, separated by the diagonal of open formations, are two ecoregions that comprise the most diverse tropical forests in the world. The Sphaenorhynchini tribe is among the few tribes of anurans that occur in both rainforests, and their historical biogeographic have never been proposed. In this study, we infer a dated phylogeny for the species of the Sphaenorhynchini and we reconstructed the biogeographic history describing the diversification chronology, and possible patterns of dispersion and vicariance, providing information about how orogeny, forest dynamics and allopatric speciation affected their evolution in South America. We provided a dated phylogeny and biogeography study for the Sphaenorhynchini tribe using mitochondrial and nuclear genes. We analyzed 41 samples to estimate the ancestral areas using biogeographical analysis based on the estimated divergence times and the current geographical ranges of the species of Sphaenorhynchini. We recovered three characteristic clades that we recognize as groups of species (S. lacteus, S. planicola, and S. platycephalus groups), with S. carneus and G. pauloalvini being the sister taxa of all other species from the tribe. We found that the diversification of the tribe lineages coincided with the main climatic and geological factors that shaped the Neotropical landscape during the Cenozoic. The most recent common ancestor of the Sphaenorhynchini species emerged in the North of the Atlantic Forest and migrated to the Amazonia in different dispersion events that occurred during the connections between these ecoregions. This is the first large‐scale study to include an almost complete calibrated phylogeny of Sphaenorhynchini, presenting important information about the evolution and diversification of the tribe. Overall, we suggest that biogeographic historical of Sphaenorhynchini have resulted from a combination of repeated range expansion and contraction cycles concurrent with climate fluctuations and dispersal events between the Atlantic Forest and Amazonia.     

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.     

Babesia bigemina,Babesia argentina, and Anaplasma marginale: Coinfectious immunity in bovines

Forty-eight intact and eight splenectomized cattle were used to evaluate different systems of coinfectious immunization against Babesia bigemina, Babesia argentina, and Anaplasma marginale. Coinfectious immunity was induced by two methods: (1) blood of cattle acutely infected with B. bigemina, B. argentina and A. marginale was used as the source of inoculum and the post vaccination reactions were chemotherapeutically controlled with Imidocarb, Ganaseg, Gloxazone, and Liquamycin, and (2) by artificially inducing babesiosis with the blood of carrier cattle with chronic infections of B. bigemina and B. argentina without chemotherapy. The degree of resistance was determined by bloodborne and tick-borne challenges. Ticks were collected from cattle and identified as Boophilus microplus and Dermacentor nitens. Vaccinated cattle demonstrated a high degree of resistance to babesiosis and anaplasmosis; however, cattle without coinfectious immunity were treated chemotherapeutically to prevent death losses.     

Genetic Diversity and Population Structure of Miscanthus sinensis Germplasm in China

Miscanthus is a perennial rhizomatous C4 grass native to East Asia. Endowed with great biomass yield, high ligno-cellulose composition, efficient use of radiation, nutrient and water, as well as tolerance to stress, Miscanthus has great potential as an excellent bioenergy crop. Despite of the high potential for biomass production of the allotriploid hybrid M. ×giganteus, derived from M. sacchariflorus and M. sinensis, other options need to be explored to improve the narrow genetic base of M. ×giganteus, and also to exploit other Miscanthus species, including M. sinensis (2n = 2x = 38), as bioenergy crops. In the present study, a large number of 459 M. sinensis accessions, collected from the wide geographical distribution regions in China, were genotyped using 23 SSR markers transferable from Brachypodium distachyon. Genetic diversity and population structure were assessed. High genetic diversity and differentiation of the germplasm were observed, with 115 alleles in total, a polymorphic rate of 0.77, Nei’s genetic diversity index (He) of 0.32 and polymorphism information content (PIC) of 0.26. Clustering of germplasm accessions was primarily in agreement with the natural geographic distribution. AMOVA and genetic distance analyses confirmed the genetic differentiation in the M. sinensis germplasm and it was grouped into five clusters or subpopulations. Significant genetic variation among subpopulations indicated obvious genetic differentiation in the collections, but within-subpopulation variation (83%) was substantially greater than the between-subpopulation variation (17%). Considerable phenotypic variation was observed for multiple traits among 300 M. sinensis accessions. Nine SSR markers were found to be associated with heading date and biomass yield. The diverse Chinese M. sinensis germplasm and newly identified SSR markers were proved to be valuable for breeding Miscanthus varieties with desired bioenergy traits.     

Climate Warming Since the Holocene Accelerates West–East Communication for the Eurasian Temperate Water Strider Species Aquarius paludum

Holocene climate warming has dramatically altered biological diversity and distributions. Recent human-induced emissions of greenhouse gases will exacerbate global warming and thus induce threats to cold-adapted taxa. However, the impacts of this major climate change on transcontinental temperate species are still poorly understood. Here, we generated extensive genomic datasets for a water strider, Aquarius paludum, which was sampled across its entire distribution in Eurasia and used these datasets in combination with ecological niche modeling (ENM) to elucidate the influence of the Holocene and future climate warming on its population structure and demographic history. We found that A. paludum consisted of two phylogeographic lineages that diverged in the middle Pleistocene, which resulted in a “west–east component” genetic pattern that was probably triggered by Central Asia-Mongoxin aridification and Pleistocene glaciations. The diverged western and eastern lineages had a second contact in the Holocene, which shaped a temporary hybrid zone located at the boundary of the arid–semiarid regions of China. Future predictions detected a potentially novel northern corridor to connect the western and eastern populations, indicating west–east gene flow would possibly continue to intensify under future warming climate conditions. Further integrating phylogeographic and ENM analyses of multiple Eurasian temperate taxa based on published studies reinforced our findings on the “west–east component” genetic pattern and the predicted future northern corridor for A. paludum. Our study provided a detailed paradigm from a phylogeographic perspective of how transcontinental temperate species differ from cold-adapted taxa in their response to climate warming.     

The Influence of Pleistocene Climatic Changes and Ocean Currents on the Phylogeography of the Southern African Barnacle,Tetraclita serrata (Thoracica; Cirripedia)

The evolutionary effects of glacial periods are poorly understood for Southern Hemisphere marine intertidal species, particularly obligatory sessile organisms. We examined this by assessing the phylogeographic patterns of the southern African volcano barnacle, Tetraclita serrata, a dominant species on rocky intertidal shores. Restricted gene flow in some geographical areas was hypothesized based on oceanic circulation patterns and known biogeographic regions. Barnacle population genetic structure was investigated using the mitochondrial cytochrome oxidase subunit 1 (COI) region for 410 individuals sampled from 20 localities spanning the South African coast. The mtDNA data were augmented by generating nuclear internal transcribed spacer 1 (ITS1) sequences from a subset of samples. Phylogenetic and population genetic analyses of mitochondrial DNA data reveal two distinct clades with mostly sympatric distributions, whereas nuclear analyses reveal only a single lineage. Shallow, but significant structure (0.0041–0.0065, P<0.01) was detected for the mtDNA data set, with the south-west African region identified as harbouring the highest levels of genetic diversity. Gene flow analyses on the mtDNA data show that individuals sampled in south-western localities experience gene flow primarily in the direction of the Benguela Current, while south and eastern localities experience bi-directional gene flow, suggesting an influence of both the inshore currents and the offshore Agulhas Current in the larval distribution of T. serrata. The mtDNA haplotype network, Bayesian Skyline Plots, mismatch distributions and time since expansion indicate that T. serrata population numbers were not severely affected by the Last Glacial Maximum (LGM), unlike other southern African marine species. The processes resulting in the two morphologically cryptic mtDNA lineages may be the result of a recent historical allopatric event followed by secondary contact or could reflect selective pressures due to differing environmental conditions.     

Subspecies divergence and pronounced phylogenetic incongruence in the East-Asia-endemic shrub Magnolia sieboldii

Background and AimsThe biogeographic patterns of the East-Asia-endemic shrub Magnolia sieboldii, in which the range of the subsp. sieboldii is interposed with the disjunct distribution of subsp. japonica, implies a complex evolutionary history, involving rapid speciation and hybridization. Here, we aim to reveal the evolutionary and phylogeographic histories of the species with a particular focus on the time of subspecies divergence, the hypothesis of secondary hybridization and the Pleistocene survival of each subspecies, using a combination of genetic analyses and ecological niche modelling.MethodsGenetic variation, genetic structures and phylogenetic relationships were elucidated based on nuclear low-copy genes, chloroplast DNA, and nuclear simple sequence repeats (SSRs). A scenario selection analysis and divergence time estimation were performed using coalescent simulation in DIYABC and *BEAST. Ecological niche modelling and a test of niche differentiation were performed using Maxent and ENMTools.Key ResultsAll marker types showed deep, but pronouncedly incongruent, west–east genetic divergences, with the subspecies being delineated only by the nuclear low-copy genes. Phylogenetic tree topologies suggested that ancient hybridization and introgression were likely to have occurred; however, this scenario did not receive significant support in the DIYABC analysis. The subspecies differentiated their niches, but both showed a dependence on high humidity and were predicted to have persisted during the last glacial cycle by maintaining a stable latitudinal distribution via migration to lower altitudes.ConclusionsWe found a deep genetic divergence and a pronounced phylogenetic incongruence among the two subspecies of M. sieboldii, which may have been driven by major paleogeographic and paleoclimatic events that have occurred since the Neogene in East Asia, including global cooling, climate oscillations and the formation of land bridges. Both subspecies were, however, considered to persist in situ in stable climatic conditions during the late Pleistocene.     

Biogeographic,climatic and spatial drivers differentially affect α-, β- and γ-diversities on oceanic archipelagos

Island biogeographic studies traditionally treat single islands as units of analysis. This ignores the fact that most islands are spatially nested within archipelagos. Here, we took a fundamentally different approach and focused on entire archipelagos using species richness of vascular plants on 23 archipelagos worldwide and their 174 constituent islands. We assessed differential effects of biogeographic factors (area, isolation, age, elevation), current and past climate (temperature, precipitation, seasonality, climate change velocity) and intra-archipelagic spatial structure (archipelago area, number of islands, area range, connectivity, environmental volume, inter-island distance) on plant diversity. Species diversity of each archipelago (γ) was additively partitioned into α, β, nestedness and replacement β-components to investigate the relative importance of environmental and spatial drivers. Multiple regressions revealed strong effects of biogeography and climate on α and γ, whereas spatial factors, particularly number of islands, inter-island distance and area range, were key to explain β. Structural equation models additionally suggested that γ is predominantly determined by indirect abiotic effects via its components, particularly β. This highlights that β and the spatial arrangement of islands are essential to understand insular ecology and evolution. Our methodological framework can be applied more widely to other taxa and archipelago-like systems, allowing new insights into biodiversity origin and maintenance.     

Population genetic structure in the endemic cyprinid fish Microphysogobio alticorpus in Taiwan: Evidence for a new phylogeographical area

The Formosa Bank is a significant zoogeographical barrier for the freshwater fish in Southern Taiwan. Here, we present population genetic structure analyses on Microphysogobio alticorpus populations in Taiwan, and biogeographic pattern to infer the relationship between historic dispersal geological dynamics. A total of 24 haplotypes were genotyped for mitochondrial CR + cyt b region in 134 specimens collected from 9 populations. Relatively high levels of haplotype diversity (h = 0.896) and low levels of nucleotide diversity (π = 0.005) were detected in M. alticorpus. Two major phylogenetic haplotype groups, Clade N1+N2 and S, were revealed and estimated to diverge about 0.121 myr (million years) ago. Haplotype network and population genetic analyses indicated significant genetic structure (F_ST = 0.876), largely concordant with the geographical location of the populations. According to SAMOVA analysis, M. alticorpus populations were divided into five units that can be explained by episodes of dispersal and population expansion followed by population fragmentation and restricted gene flow. The results of the topological pattern and network suggested that the Kaoping River was the first habitat colonized by the ancestor of M. alticorpus, before the northward dispersal.     

Evolution of genes involved in the unusual genitals of the bear macaque,Macaca arctoides

Genital divergence is thought to contribute to reproductive barriers by establishing a “lock‐and‐key" mechanism for reproductive compatibility. One such example, Macaca arctoides, the bear macaque, has compensatory changes in both male and female genital morphology as compared to close relatives. M. arctoides also has a complex evolutionary history, having extensive introgression between the fascicularis and sinica macaque species groups. Here, phylogenetic relationships were analyzed via whole‐genome sequences from five species, including M. arctoides, and two species each from the putative parental species groups. This analysis revealed ~3x more genomic regions supported placement in the sinica species group as compared to the fascicularis species group. Additionally, introgression analysis of the M. arctoides genome revealed it is a mosaic of recent polymorphisms shared with both species groups. To examine the evolution of their unique genital morphology further, the prevalence of candidate genes involved in genital morphology was compared against genome‐wide outliers in various population genetic metrics of diversity, divergence, introgression, and selection, while accounting for background variation in recombination rate. This analysis identified 67 outlier genes, including several genes that influence baculum morphology in mice, which were of interest since the bear macaque has the longest primate baculum. The mean of four of the seven population genetic metrics was statistically different in the candidate genes as compared to the rest of the genome, suggesting that genes involved in genital morphology have increased divergence and decreased diversity beyond expectations. These results highlight specific genes that may have played a role in shaping the unique genital morphology in the bear macaque.     

The Tien Shan vole (Microtus ilaeus; Rodentia: Cricetidae) as a new species in the Late Pleistocene of Europe

Grey voles (subgenus Microtus) represent a complex of at least seven closely related and partly cryptic species. The range of these species extends from the Atlantic to the Altai Mountains, but most of them occur east of the Black Sea. Using ancient DNA analyses of the Late Pleistocene specimens, we identified a new mtDNA lineage of grey voles in Europe. Phylogenetic analysis of mitochondrial DNA cytochrome b sequences from 23 voles from three caves, namely, Emine–Bair–Khosar (Crimea, Ukraine), Cave 16 (Bulgaria), and Bacho Kiro (Bulgaria), showed that 14 specimens form a previously unrecognized lineage, sister to the Tien Shan vole. The average sequence divergence of this lineage and the extant Tien Shan vole was 4.8%, which is similar to the divergence of grey vole forms, which are considered distinct species or being on the verge of speciation; M. arvalis and M. obscurus or M. mystacinus and M. rossiaemeridionalis. We estimated the time to the most recent common ancestor of the grey voles to be 0.66 Ma, which is over twice the recent estimates, while the divergence of the extant Tien Shan vole and the new lineage to be 0.29 Ma. Our discovery suggests that grey voles may have been more diversified in the past and that their ranges may have differed substantially from current ones. It also underlines the utility of ancient DNA to decipher the evolutionary history of voles.     

Ancient and recent introgression shape the evolutionary history of pollinator adaptation and speciation in a model monkeyflower radiation (Mimulus section Erythranthe)

Inferences about past processes of adaptation and speciation require a gene-scale and genome-wide understanding of the evolutionary history of diverging taxa. In this study, we use genome-wide capture of nuclear gene sequences, plus skimming of organellar sequences, to investigate the phylogenomics of monkeyflowers in Mimulus section Erythranthe (27 accessions from seven species). Taxa within Erythranthe, particularly the parapatric and putatively sister species M. lewisii (bee-pollinated) and M. cardinalis (hummingbird-pollinated), have been a model system for investigating the ecological genetics of speciation and adaptation for over five decades. Across >8000 nuclear loci, multiple methods resolve a predominant species tree in which M. cardinalis groups with other hummingbird-pollinated taxa (37% of gene trees), rather than being sister to M. lewisii (32% of gene trees). We independently corroborate a single evolution of hummingbird pollination syndrome in Erythranthe by demonstrating functional redundancy in genetic complementation tests of floral traits in hybrids; together, these analyses overturn a textbook case of pollination-syndrome convergence. Strong asymmetries in allele sharing (Patterson’s D-statistic and related tests) indicate that gene tree discordance reflects ancient and recent introgression rather than incomplete lineage sorting. Consistent with abundant introgression blurring the history of divergence, low-recombination and adaptation-associated regions support the new species tree, while high-recombination regions generate phylogenetic evidence for sister status for M. lewisii and M. cardinalis. Population-level sampling of core taxa also revealed two instances of chloroplast capture, with Sierran M. lewisii and Southern Californian M. parishii each carrying organelle genomes nested within respective sympatric M. cardinalis clades. A recent organellar transfer from M. cardinalis, an outcrosser where selfish cytonuclear dynamics are more likely, may account for the unexpected cytoplasmic male sterility effects of selfer M. parishii organelles in hybrids with M. lewisii. Overall, our phylogenomic results reveal extensive reticulation throughout the evolutionary history of a classic monkeyflower radiation, suggesting that natural selection (re-)assembles and maintains species-diagnostic traits and barriers in the face of gene flow. Our findings further underline the challenges, even in reproductively isolated species, in distinguishing re-use of adaptive alleles from true convergence and emphasize the value of a phylogenomic framework for reconstructing the evolutionary genetics of adaptation and speciation.     

New karyotype for Mesomys stimulax (Rodentia,Echimyidae) from the Brazilian Amazon: A case for species complex?

Mesomys Wagner, 1845 (Rodentia, Echimyidae, Eumysopinae) currently has four recognized species, three of which occur in Brazil: Mesomys hispidus (probably a species complex), M. occultus, and M. stimulax. Mesomys leniceps is found in montane forests of northern Peru. Mesomys stimulax, the focus of the present study, has a distribution that is restricted to the central and eastern Amazonia south of the Amazon River, extending from the left bank of the Tapajós River to the right bank of the Tocantins River, and south to the southeast portion of Pará State. The genus presents karyotypes with diploid number 2n = 60 and Fundamental Number (FN) = 116 for M. hispidus and M. stimulax, and 2n = 42, FN = 54 for M. occultus. We studied the karyotype of a female specimen of M. stimulax collected from the Tapirapé‐Aquiri National Forest, Marabá, Pará, Brazil, in the Xingu/Tocantins interfluvium. The obtained karyotype (2n = 60 and FN = 110) differs from that described in the literature for both M. stimulax and M. hispidus by exhibiting more biarmed chromosomes, probably due to pericentric inversions and/or centromeric repositioning, and exhibiting differences in the amount and distribution of constitutive heterochromatin (CH). These results suggest that, similar to what has already been proposed for M. hispidus, M. stimulax may represent a species complex and/or cryptic species. The mechanisms of chromosomal diversification in Mesomys and the biogeographic implications are discussed reinforcing the need for broad systematic review for Mesomys.     

Combining phylogeography and landscape genetics to infer the evolutionary history of a short-range Mediterranean relict, <Emphasis Type="Italic">Salamandra salamandra longirostris</Emphasis>

Examining historical and contemporary processes underlying current patterns of genetic variation is key to reconstruct the evolutionary history of species and implement conservation measures promoting their long-term persistence. Combining phylogeographic and landscape genetic approaches can provide valuable insights, especially in regions harboring high levels of biodiversity that are currently threatened by climate and land cover changes, like southern Iberia. We used genetic (mtDNA and microsatellites) and spatial data (climate and land cover) to infer the evolutionary history and contemporary genetic connectivity in a short-range endemic salamander subspecies, Salamandra salamandra longirostris, using a combination of ecological niche modelling, phylogeographic, and landscape genetic analyses. Ecological-based analyses support a role of the Guadalquivir River Basin as a major vicariant agent in this taxon. The lower genetic diversity and greater differentiation of peripheral populations, together with analyses of climatically stable areas throughout time, suggest the persistence of a population in the central part of the current range since the Last Inter Glacial [LIG; ~?120,000–140,000 years BP], and a micro refugium in the eastern end of the range. Habitat heterogeneity plays a major role in shaping patterns of genetic differentiation in S. s. longirostris, with forests representing key areas for its long-term persistence under scenarios of environmental change. Our study stresses the importance of maintaining population genetic connectivity in low-dispersal organisms under rapidly changing environments, and will inform management plans for the long-term survival of this evolutionarily distinct Mediterranean endemic.