A genealogical approach to quantifying lineage divergence

Recombination rate variation shapes genomic variability of phylogeographic structure in a widespread North American songbird (Aves: Certhia americana)

The nonrandom distribution of chromosomal characteristics and functional elements-genomic architecture-impacts the relative strengths and impacts of population genetic processes across the genome. Due to this relationship, genomic architecture has the potential to shape variation in population genetic structure across the genome. Population genetic structure has been shown to vary across the genome in a variety of taxa, but this body of work has largely focused on pairwise population genomic comparisons between closely related taxa. Here, we used whole genome sequencing of seven phylogeographically structured populations of a North American songbird, the Brown Creeper (Certhia americana), to determine the impacts of genomic architecture on phylogeographic structure variation across the genome. Using multiple methods to infer phylogeographic structure-ordination, clustering, and phylogenetic methods-we found that recombination rate variation explained a large proportion of phylogeographic structure variation. Genomic regions with low recombination showed phylogeographic structure consistent with the genome-wide pattern. In regions with high recombination, we found strong phylogeographic structure, but with discordant patterns relative to the genome-wide pattern. In regions with high recombination rate, we found that populations with small effective population sizes evolve relatively more rapidly than larger populations, leading to discordant signatures of phylogeographic structure. These results suggest that the interplay between recombination rate variation and effective population sizes shape the relative impacts of selection and genetic drift in different parts of the genome. Overall, the combined interactions of population genetic processes, genomic architecture, and effective population sizes shape patterns of variability in phylogeographic structure across the genome of the Brown Creeper.

Impact of Analytical Treatment Interruption on Burden and Diversification of HIV Peripheral Reservoir: A Pilot Study

BACKGROUND

If analytical antiretroviral-treatment (ART) interruption (ATI) might significantly impact quantitative or qualitative peripheral-total HIV-DNA is still debated.

METHODS

Six chronically HIV-1 infected patients enrolled in APACHE-study were analysed for peripheral-total HIV-DNA and residual viremia, major-resistance-mutations (MRMs) and C2-V3-C3 evolution at pre-ATI (T1), during ATI (T2) and at achievement of virological success after ART-resumption (post-ATI, T3). These data were obtained at three comparable time-points in five chronically HIV-1 infected patients on suppressive ART for ≥1 year, enrolled in MODAt-study.

RESULTS

At T1, APACHE and MODAt individuals had similar peripheral-total HIV-DNA and residual viremia (p = 0.792 and 0.662, respectively), and no significant changes for these parameters were observed between T1 and T3 in both groups. At T1, 4/6 APACHE and 2/5 MODAt carried HIV-DNA MRMs. MRMs disappeared at T3 in 3/4 APACHE. All disappearing MRMs were characterized by T1 intra-patient prevalence <80%, and mainly occurred in APOBEC3-related sites. All MRMs persisted over-time in the 2 MODAt. C2-V3-C3 genetic-distance significantly changed from T1 to T3 in APACHE individuals (+0.36[0.11-0.41], p = 0.04), while no significant changes were found in MODAt. Accordingly, maximum likelihood trees (bootstrap > 70%) and genealogical sorting indices (GSI > 0.50 with p-value < 0.05) showed that T1 C2-V3-C3 DNA sequences were distinct from T2 and T3 viruses in 4/6 APACHE. Virus populations at all three time-points were highly interspersed in MODAt.

CONCLUSIONS

This pilot study indicates that short ATI does not alter peripheral-total HIV-DNA burden and residual viremia, but in some cases could cause a genetic diversification of peripheral viral reservoir in term of both MRMs rearrangement and viral evolution.

Incorporating the speciation process into species delimitation

The "multispecies" coalescent (MSC) model that underlies many genomic species-delimitation approaches is problematic because it does not distinguish between genetic structure associated with species versus that of populations within species. Consequently, as both the genomic and spatial resolution of data increases, a proliferation of artifactual species results as within-species population lineages, detected due to restrictions in gene flow, are identified as distinct species. The toll of this extends beyond systematic studies, getting magnified across the many disciplines that rely upon an accurate framework of identified species. Here we present the first of a new class of approaches that addresses this issue by incorporating an extended speciation process for species delimitation. We model the formation of population lineages and their subsequent development into independent species as separate processes and provide for a way to incorporate current understanding of the species boundaries in the system through specification of species identities of a subset of population lineages. As a result, species boundaries and within-species lineages boundaries can be discriminated across the entire system, and species identities can be assigned to the remaining lineages of unknown affinities with quantified probabilities. In addition to the identification of species units in nature, the primary goal of species delimitation, the incorporation of a speciation model also allows us insights into the links between population and species-level processes. By explicitly accounting for restrictions in gene flow not only between, but also within, species, we also address the limits of genetic data for delimiting species. Specifically, while genetic data alone is not sufficient for accurate delimitation, when considered in conjunction with other information we are able to not only learn about species boundaries, but also about the tempo of the speciation process itself.

Fine-Scale Haplotype Structure Reveals Strong Signatures of Positive Selection in a Recombining Bacterial Pathogen

Identifying genetic variation in bacteria that has been shaped by ecological differences remains an important challenge. For recombining bacteria, the sign and strength of linkage provide a unique lens into ongoing selection. We show that derived alleles <300 bp apart in Neisseria gonorrhoeae exhibit more coupling linkage than repulsion linkage, a pattern that cannot be explained by limited recombination or neutrality as these couplings are significantly stronger for nonsynonymous alleles than synonymous alleles. This general pattern is driven by a small fraction of highly diverse genes, many of which exhibit evidence of interspecies horizontal gene transfer and an excess of intermediate frequency alleles. Extensive simulations show that two distinct forms of positive selection can create these patterns of genetic variation: directional selection on horizontally transferred alleles or balancing selection that maintains distinct haplotypes in the presence of recombination. Our results establish a framework for identifying patterns of selection in fine-scale haplotype structure that indicate specific ecological processes in species that recombine with distantly related lineages or possess coexisting adaptive haplotypes.

Optimal markers for the identification of Colletotrichum species

Colletotrichum is among the most important genera of fungal plant pathogens. Molecular phylogenetic studies over the last decade have resulted in a much better understanding of the evolutionary relationships and species boundaries within the genus. There are now approximately 200 species accepted, most of which are distributed among 13 species complexes. Given their prominence on agricultural crops around the world, rapid identification of a large collection of Colletotrichum isolates is routinely needed by plant pathologists, regulatory officials, and fungal biologists. However, there is no agreement on the best molecular markers to discriminate species in each species complex. Here we calculate the barcode gap distance and intra/inter-specific distance overlap to evaluate each of the most commonly applied molecular markers for their utility as a barcode for species identification. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), histone-3 (HIS3), DNA lyase (APN2), intergenic spacer between DNA lyase and the mating-type locus MAT1-2-1 (APN2/MAT-IGS), and intergenic spacer between GAPDH and a hypothetical protein (GAP2-IGS) have the properties of good barcodes, whereas sequences of actin (ACT), chitin synthase (CHS-1) and nuclear rDNA internal transcribed spacers (nrITS) are not able to distinguish most species. Finally, we assessed the utility of these markers for phylogenetic studies using phylogenetic informativeness profiling, the genealogical sorting index (GSI), and Bayesian concordance analyses (BCA). Although GAPDH, HIS3 and β-tubulin (TUB2) were frequently among the best markers, there was not a single set of markers that were best for all species complexes. Eliminating markers with low phylogenetic signal tends to decrease uncertainty in the topology, regardless of species complex, and leads to a larger proportion of markers that support each lineage in the Bayesian concordance analyses. Finally, we reconstruct the phylogeny of each species complex using a minimal set of phylogenetic markers with the strongest phylogenetic signal and find the majority of species are strongly supported as monophyletic.

Molecular phylogeny and species delimitation of the genus Dicerapanorpa (Mecoptera: Panorpidae)

Given that species is the fundamental unit in systematic biology, rigorous species delimitation is crucial for taxonomic studies, yet routine species delimitation remains an ongoing challenge in the taxonomic practice of insects. The two-horned scorpionfly Dicerapanorpa is a small genus in Panorpidae (Mecoptera) endemic to the Qinling-Bashan and Hengduan mountains, a biodiversity hotspot. However, species of Dicerapanorpa are difficult to delineate owing to marked intraspecific variation and interspecific similarity. Here, we investigate the diversity and species boundaries of Dicerapanorpa using an integrative approach based on DNA barcoding, morphological, geometric morphometric and molecular phylogenetic analyses. This integrative analyses confirmed the 13 described species of Dicerapanorpa and revealed three new species: Dicerapanorpa lativalva sp. nov., Dicerapanorpa hualongshana sp. nov. and Dicerapanorpa minshana sp. nov. Most molecular operational taxonomic units are in congruence with morphological clusters. Possible reasons for several discordances in Dicerapanorpa are tentatively discussed.

Unexpected species diversity in electric eels with a description of the strongest living bioelectricity generator

Is there only one electric eel species? For two and a half centuries since its description by Linnaeus, Electrophorus electricus has captivated humankind by its capacity to generate strong electric discharges. Despite the importance of Electrophorus in multiple fields of science, the possibility of additional species-level diversity in the genus, which could also reveal a hidden variety of substances and bioelectrogenic functions, has hitherto not been explored. Here, based on overwhelming patterns of genetic, morphological, and ecological data, we reject the hypothesis of a single species broadly distributed throughout Greater Amazonia. Our analyses readily identify three major lineages that diverged during the Miocene and Pliocene-two of which warrant recognition as new species. For one of the new species, we recorded a discharge of 860 V, well above 650 V previously cited for Electrophorus, making it the strongest living bioelectricity generator.

Multiple convergent events created a nominal widespread species: Triplophysa stoliczkae (Steindachner, 1866) (Cobitoidea: Nemacheilidae)

BACKGROUND

Triplophysa stoliczkae is the most widespread species in the genus Triplophysa and may have originated from morphological convergence. To understand the evolutionary history of T. stoliczkae, we employed a multilocus approach to investigate the phylogenetics and the morphological evolution of T. stoliczkae on the Qinghai-Tibetan Plateau.

RESULTS

All phylogenetic analyses (two mitochondrial and five nuclear loci), a genealogical sorting index and species tree inferences suggested that T. stoliczkae consists of distinct lineages that were not closest relatives. The time estimation indicated that the divergence events between "T. stoliczkae" and other Triplophysa species occurred from approximately 0.10 to 4.51 Ma. The ancestral state analyses supported the independent evolution of T. stoliczkae morphology in distinct lineages. The morphometric analysis and convergence estimates demonstrated significant phenotypic convergence among "T. stoliczkae" lineages.

CONCLUSIONS

Triplophysa stoliczkae includes 4 different lineages with similar morphologies. The increasingly harsh environments that have occurred since the Pliocene have driven the occurrences of scrape-feeding fish in the genus Triplophysa. Morphological adaptations associated with scrape-feeding behavior resulted in convergences and the artificial lumping of four different species in the nominal taxon T. stoliczkae. A taxonomic revision for T. stoliczkae is needed.

'At the crossroads towards polyploidy': genomic divergence and extent of homoploid hybridization are drivers for the formation of the ox-eye daisy polyploid complex (Leucanthemum, Compositae-Anthemideae)

Polyploidy plays a paramount role in phytodiversity, but the causes of this evolutionary pathway require further study. Here, we use phylogenetic methods to examine possible polyploidy-promoting factors by comparing diploid representatives of the comprehensive European polyploid complex Leucanthemum with members of its strictly diploid North African counterpart Rhodanthemum. We investigate genetic divergence and gene flow among all diploid lineages of both genera to evaluate the role of genomic differentiation and hybridization for polyploid speciation. To test whether hybridization in Leucanthemum has been triggered by the geological conditions during its diversification, we additionally generate a time-calibrated phylogeny of 46 species of the subtribe Leucantheminae. Leucanthemum shows a significantly higher genetic divergence and hybridization signal among diploid lineages compared with Rhodanthemum, in spite of a similar crown age and diversification pattern during the Quaternary. Our study demonstrates the importance of genetic differentiation among diploid progenitors and their concurrent affinity for natural hybridization for the formation of a polyploid complex. Furthermore, the role of climate-induced range overlaps on hybridization and polyploid speciation during the Quaternary is discussed.

Unraveling patterns and processes of diversification in the South Andean-Patagonian Nassauvia subgenus Strongyloma (Asteraceae, Nassauvieae)

Congruence among different sources of data is highly desirable in phylogenetic analyses. However, plastid and nuclear DNA may record different evolutionary processes such that incongruence among results from these sources can help unravel complex evolutionary histories. That is the case of Nassauvia subgenus Strongyloma (Asteraceae), a taxon with five putative species distributed in the southern Andes and Patagonian steppe. Morphometric and phylogeographic information cast doubt on the integrity of its species, and previous molecular data even questioned the monophyly of the subgenus. We tested those questions using plastid and nuclear DNA sequences by the application of different methods such as phylogenetic trees, networks, a test of genealogical sorting, an analysis of population structure, calibration of the trees, and hybridization test, assembling non-synchronous incongruent results at subgenus and species levels in a single reconstruction. The integration of our molecular analyses and previous taxonomic, morphological, and molecular studies support subgenus Strongyloma as a monophyletic group. However, the topology of the nuclear trees and the evidence of polyploids within subgenus Nassauvia, suggest a hypothetical origin and initial radiation of Nassauvia related to an ancient hybridization event that occurred around 17-6.3 Myr ago near the Andes in west-central Patagonia. Plastid data suggest a recent diversification within subgenus Strongyloma, at most 9.8 Myr ago, towards the Patagonian steppe east of the Andes. These processes cause phylogenies to deviate from the species tree since each putative species lack exclusive ancestry. The non-monophyly of its species using both plastid and nuclear data is caused mainly by incomplete lineage sorting occurred since the Miocene. The final uplift of the Andes and Pliocene-Pleistocene glacial-interglacial and its consequences on the landscape and climate structured the genetic composition of this group of plants in the Patagonian steppe. The molecular data presented here agree with previous morphological studies, in that the five putative species typically accepted in this subgenus are not independent taxa. This study emphasizes that adding more than one sequence per species, not combining data with dissimilar inheritance patterns without first performed incongruence tests, exploring data through different methodologies, considering the timing of events, and searching for the causes of poorly resolved and/or incongruent phylogenies help to reveal complex biological underlying processes, which might otherwise remain hidden.

Multilocus phylogeography of the brown-spotted pitviper Protobothrops mucrosquamatus (Reptilia: Serpentes: Viperidae) sheds a new light on the diversification pattern in Asia

Understanding the influence of geographical events and climate changes on genetic diversity is essential in explaining current patterns of genetic structure and geographic distribution of organisms. We inferred phylogenetic relationships, investigated historical demography, explored the evolutionary history, and clarified intraspecific taxonomy of Protobothrops mucrosquamatus, which is one of the commonest and most wide-ranging Asian pitvipers. A total of 184 samples from 54 localities were sequenced and analyzed for two mitochondrial gene fragments and two nuclear genes. Phylogenetic reconstruction based on mtDNA sequences revealed the existence of a minimum of five geographically structured and well-supported lineages within P. mucrosquamatus. Based on the mtDNA gene tree, and the geographic relationship between populations allied by matrilineal lineages, a complex longitudinal and latitudinal diversification pattern was uncovered in P. mucrosquamatus. The estimated date of the origin of the species (about 5.3 Ma) and divergence of the intraspecific lineages match the rapid uplifting of Qinghai-Xizang Plateau, and is also consistent with those of some other co-distributed Asian pitvipers. Formation of the two island lineages (Taiwan and Hainan) was generally congruent with the first isolation of the islands, but the two lineages showed different relationships with the continental Asian populations in comparison with some other pitvipers. Population historical demographic analyses, based on three methods, showed that all lineages have experienced slight population expansion in and around the Dali Glacial. Tests of intraspecific taxonomy indicated that no cryptic taxon is present within this widely distributed snake.

Lichen chemistry is concordant with multilocus gene genealogy in the genus Cetrelia (Parmeliaceae, Ascomycota)

The lichen genus Cetrelia represents a taxonomically interesting case where morphologically almost uniform populations differ considerably from each other chemically. Similar variation is not uncommon among lichenized fungi, but it is disputable whether such populations should be considered entities at the species level. Species boundaries in Cetrelia are traditionally delimited either as solely based on morphology or as combinations of morpho- and chemotypes. A dataset of four nuclear markers (ITS, IGS, Mcm7, RPB1) from 62 specimens, representing ten Cetrelia species, was analysed within Bayesian and maximum likelihood frameworks. Analyses recovered a well-resolved phylogeny where the traditional species generally were monophyletic, with the exception of Cetrelia chicitae and Cetrelia pseudolivetorum. Species delimitation analyses supported the distinction of 15 groups within the studied Cetrelia taxa, dividing three traditionally identified species into some species candidates. Chemotypes, distinguished according to the major medullary substance, clearly correlated with clades recovered within Cetrelia, while samples with the same reproductive mode were dispersed throughout the phylogenetic tree. Consequently, delimiting Cetrelia species based only on reproductive morphology is not supported phylogenetically. Character analyses suggest that chemical characters have been more consistent compared to reproductive mode and indicate that metabolite evolution in Cetrelia towards more complex substances is probable.

Azithromycin Resistance through Interspecific Acquisition of an Epistasis-Dependent Efflux Pump Component and Transcriptional Regulator in Neisseria gonorrhoeae

Mosaic interspecifically acquired alleles of the multiple transferable resistance (mtr) efflux pump operon correlate with increased resistance to azithromycin in Neisseria gonorrhoeae in epidemiological studies. However, whether and how these alleles cause resistance is unclear. Here, we use population genomics, transformations, and transcriptional analyses to dissect the relationship between variant mtr alleles and azithromycin resistance. We find that the locus encompassing the mtrR transcriptional repressor and the mtrCDE pump is a hot spot of interspecific recombination introducing alleles from Neisseria meningitidis and Neisseria lactamica into N. gonorrhoeae, with multiple rare haplotypes in linkage disequilibrium at mtrD and the mtr promoter region. Transformations demonstrate that resistance to azithromycin, as well as to other antimicrobial compounds such as polymyxin B and crystal violet, is mediated through epistasis between these two loci and that the full-length mosaic mtrD allele is required. Gene expression profiling reveals the mechanism of resistance in mosaics couples novel mtrD alleles with promoter mutations that increase expression of the pump. Overall, our results demonstrate that epistatic interactions at mtr gained from multiple neisserial species has contributed to increased gonococcal resistance to diverse antimicrobial agents.IMPORTANCENeisseria gonorrhoeae is the sexually transmitted bacterial pathogen responsible for more than 100 million cases of gonorrhea worldwide each year. The incidence of resistance to the macrolide azithromycin has increased in the past decade; however, a large proportion of the genetic basis of resistance remains unexplained. This study is the first to conclusively demonstrate the acquisition of macrolide resistance through mtr alleles from other Neisseria species, demonstrating that commensal Neisseria bacteria are a reservoir for antibiotic resistance to macrolides, extending the role of interspecies mosaicism in resistance beyond what has been previously described for cephalosporins. Ultimately, our results emphasize that future fine-mapping of genome-wide interspecies mosaicism may be valuable in understanding the pathways to antimicrobial resistance. Our results also have implications for diagnostics and public health surveillance and control, as they can be used to inform the development of sequence-based tools to monitor and control the spread of antibiotic-resistant gonorrhea.

Are species-pairs diverging lineages? A nine-locus analysis uncovers speciation among species-pairs of the Lobaria meridionalis-group (Ascomycota)

In spite of considerable effort to verify the theory of species-pairs, uncertainty still exists about the relationship between sexually or vegetatively reproducing populations of morphologically indistinguishable, sympatric lichen species. The current paper studies putative species-pairs within the Asian Lobaria meridionalis-group, using a nine-locus and time calibrated species-tree approach. Analyses demonstrate that pairs of sexually or vegetatively reproducing lineages split into highly supported monophyletic clades-confirming molecularly the species-pair concept for the L. meridionalis-group. In the broader context of evolution and speciation dynamics in lichenized fungi, this paper attempts to synthesize molecular findings from the last two decades to promote a more modern perception of the species-pair concept. Taxonomically, eight species were found to currently conform to the L. meridionalis-group, which differentiated during the Pliocene and Pleistocene. The coincidence of paleoclimatic events with estimated dates of divergence support a bioclimatic hypothesis for the evolution of species in the L. meridionalis-group, which also explains their current eco-geographic distribution patterns. Greater recognition for species with a long and independent evolutionary history, which merit high conservation priority, will be especially critical for preserving geographically restricted endemics from Southeast Asia, where habitat loss is driving rapid declines.

HIV-1 latent reservoir size and diversity are stable following brief treatment interruption

BACKGROUND

The effect of a brief analytical treatment interruption (ATI) on the HIV-1 latent reservoir of individuals who initiate antiretroviral therapy (ART) during chronic infection is unknown.

METHODS

We evaluated the impact of transient viremia on the latent reservoir in participants who underwent an ATI and at least 6 months of subsequent viral suppression in a clinical trial testing the effect of passive infusion of the broadly neutralizing Ab VRC01 during ATI.

RESULTS

Measures of total HIV-1 DNA, cell-associated RNA, and infectious units per million cells (IUPM) (measured by quantitative viral outgrowth assay [QVOA]) were not statistically different before or after ATI. Phylogenetic analyses of HIV-1 env sequences from QVOA and proviral DNA demonstrated little change in the composition of the virus populations comprising the pre- and post-ATI reservoir. Expanded clones were common in both QVOA and proviral DNA sequences. The frequency of clonal populations differed significantly between QVOA viruses, proviral DNA sequences, and the viruses that reactivated in vivo.

CONCLUSIONS

The results indicate that transient viremia from ATI does not substantially alter measures of the latent reservoir, that clonal expansion is prevalent within the latent reservoir, and that characterization of latent viruses that can reactivate in vivo remains challenging.

TRIAL REGISTRATION

ClinicalTrials.gov NCT02463227FUNDING. Funding was provided by the NIH.

Genetic diversity of the Chinese goat in the littoral zone of the Yangtze River as assessed by microsatellite and mtDNA

The objective of this study was to assess the genetic diversity and population structure of goats in the Yangtze River region using microsatellite and mtDNA to better understand the current status of those goat genetic diversity and the effects of natural landscape in fashion of domestic animal genetic diversity. The genetic variability of 16 goat populations in the littoral zone of the Yangtze River was estimated using 21 autosomal microsatellites, which revealed high diversity and genetic population clustering with a dispersed geographical distribution. A phylogenetic analysis of the mitochondrial D-loop region (482 bp) was conducted in 494 goats from the Yangtze River region. In total, 117 SNPs were reconstructed, and 173 haplotypes were identified, 94.5% of which belonged to lineages A and B. Lineages C, D, and G had lower frequencies (5.2%), and lineage F haplotypes were undetected. Several high-frequency haplotypes were shared by different ecogeographically distributed populations, and the close phylogenetic relationships among certain low-frequency haplotypes indicated the historical exchange of genetic material among these populations. In particular, the lineage G haplotype suggests that some west Asian goat genetic material may have been transferred to China via Muslim migration.

The genomic landscape at a late stage of stickleback speciation: High genomic divergence interspersed by small localized regions of introgression

Speciation is a continuous process and analysis of species pairs at different stages of divergence provides insight into how it unfolds. Previous genomic studies on young species pairs have revealed peaks of divergence and heterogeneous genomic differentiation. Yet less known is how localised peaks of differentiation progress to genome-wide divergence during the later stages of speciation in the presence of persistent gene flow. Spanning the speciation continuum, stickleback species pairs are ideal for investigating how genomic divergence builds up during speciation. However, attention has largely focused on young postglacial species pairs, with little knowledge of the genomic signatures of divergence and introgression in older stickleback systems. The Japanese stickleback species pair, composed of the Pacific Ocean three-spined stickleback (Gasterosteus aculeatus) and the Japan Sea stickleback (G. nipponicus), which co-occur in the Japanese islands, is at a late stage of speciation. Divergence likely started well before the end of the last glacial period and crosses between Japan Sea females and Pacific Ocean males result in hybrid male sterility. Here we use coalescent analyses and Approximate Bayesian Computation to show that the two species split approximately 0.68-1 million years ago but that they have continued to exchange genes at a low rate throughout divergence. Population genomic data revealed that, despite gene flow, a high level of genomic differentiation is maintained across the majority of the genome. However, we identified multiple, small regions of introgression, occurring mainly in areas of low recombination rate. Our results demonstrate that a high level of genome-wide divergence can establish in the face of persistent introgression and that gene flow can be localized to small genomic regions at the later stages of speciation with gene flow.

The role of isolation on contrasting phylogeographic patterns in two cave crustaceans

BACKGROUND

The underlying mechanisms and processes that prompt the colonisation of extreme environments, such as caves, constitute major research themes of evolutionary biology and biospeleology. The special adaptations required to survive in subterranean environments (low food availability, hypoxic waters, permanent darkness), and the geographical isolation of caves, nominate cave biodiversity as ideal subjects to answer long-standing questions concerning the interplay amongst adaptation, biogeography, and evolution. The present project aims to examine the phylogeographic patterns exhibited by two sympatric species of surface and cave-dwelling peracarid crustaceans (Asellus aquaticus and Niphargus hrabei), and in doing so elucidate the possible roles of isolation and exaptation in the colonisation and successful adaptation to the cave environment.

RESULTS

Specimens of both species were sampled from freshwater hypogean (cave) and epigean (surface) habitats in Hungary, and additional data from neighbouring countries were sourced from Genbank. Sequencing of mitochondrial and nuclear loci revealed, through haplotype network reconstruction (TCS) and phylogenetic inference, the genetic structure, phylogeographic patterns, and divergence-time estimates of A. aquaticus and N. hrabei surface and cave populations. Contrasting phylogeographic patterns were found between species, with A. aquaticus showing strong genetic differentiation between cave and surface populations and N. hrabei lacking any evidence of genetic structure mediated by the cave environment. Furthermore, N. hrabei populations show very low levels of genetic differentiation throughout their range, which suggests the possibility of recent expansion events over the last few thousand years.

CONCLUSIONS

Isolation by cave environment, rather than distance, is likely to drive the genetic structuring observed between immediately adjacent cave and surface populations of A. aquaticus, a predominantly surface species with only moderate exaptations to subterranean life. For N. hrabei, in which populations exhibit a fully 'cave-adapted' (troglomorphic) phenotype, the lack of genetic structure suggests that subterranean environments do not pose a dispersal barrier for this surface-cave species.

Genomics of adaptive divergence with chromosome-scale heterogeneity in crossover rate

Genetic differentiation between divergent populations is often greater in chromosome centres than peripheries. Commonly overlooked, this broadscale differentiation pattern is sometimes ascribed to heterogeneity in crossover rate and hence linked selection within chromosomes, but the underlying mechanisms remain incompletely understood. A literature survey across 46 organisms reveals that most eukaryotes indeed exhibit a reduced crossover rate in chromosome centres relative to the peripheries. Using simulations of populations diverging into ecologically different habitats through sorting of standing genetic variation, we demonstrate that such chromosome-scale heterogeneity in crossover rate, combined with polygenic divergent selection, causes stronger hitchhiking and especially barriers to gene flow across chromosome centres. Without requiring selection on new mutations, this rapidly leads to elevated population differentiation in the low-crossover centres relative to the high-crossover peripheries of chromosomes ("Chromosome Centre-Biased Differentiation", CCBD). Using simulated and empirical data, we then show that strong CCBD between populations can provide evidence of polygenic adaptive divergence with a phase of gene flow. We further demonstrate that chromosome-scale heterogeneity in crossover rate impacts analyses beyond that of population differentiation, including the inference of phylogenies and parallel adaptive evolution among populations, the detection of genetic loci under selection, and the interpretation of the strength of selection on genomic regions. Overall, our results call for a greater appreciation of chromosome-scale heterogeneity in crossover rate in evolutionary genomics.

Phylogeny of xerophilic aspergilli (subgenus Aspergillus) and taxonomic revision of section Restricti

Aspergillus section Restricti together with sister section Aspergillus (formerly Eurotium) comprises xerophilic species, that are able to grow on substrates with low water activity and in extreme environments. We adressed the monophyly of both sections within subgenus Aspergillus and applied a multidisciplinary approach for definition of species boundaries in sect. Restricti. The monophyly of sections Aspergillus and Restricti was tested on a set of 102 isolates comprising all currently accepted species and was strongly supported by Maximum likelihood (ML) and Bayesian inferrence (BI) analysis based on β-tubulin (benA), calmodulin (CaM) and RNA polymerase II second largest subunit (RPB2) loci. More than 300 strains belonging to sect. Restricti from various isolation sources and four continents were characterized by DNA sequencing, and 193 isolates were selected for phylogenetic analyses and phenotypic studies. Species delimitation methods based on multispecies coalescent model were employed on DNA sequences from four loci, i.e., ID region of rDNA (ITS + 28S), CaM, benA and RPB2, and supported recognition of 21 species, including 14 new. All these species were also strongly supported in ML and BI analyses. All recognised species can be reliably identified by all four examined genetic loci. Phenotype analysis was performed to support the delimitation of new species and includes colony characteristics on seven cultivation media incubated at several temperatures, growth on an osmotic gradient (six media with NaCl concentration from 0 to 25 %) and analysis of morphology including scanning electron microscopy. The micromorphology of conidial heads, vesicle dimensions, temperature profiles and growth parameters in osmotic gradient were useful criteria for species identification.

The vast majority of species in sect. Restricti produce asperglaucide, asperphenamate or both in contrast to species in sect. Aspergillus. Mycophenolic acid was detected for the first time in at least six members of the section. The ascomata of A. halophilicus do not contain auroglaucin, epiheveadride or flavoglaucin which are common in sect. Aspergillus, but shares the echinulins with sect. Aspergillus.

Isolation and gene flow in a speciation continuum in newts

Because reproductive isolation often evolves gradually, differentiating lineages may retain the potential for genetic exchange for prolonged periods, providing an opportunity to quantify and to understand the fundamental role of gene flow during speciation. Here we delimit evolutionary lineages, reconstruct the phylogeny and infer gene flow in newts of the Lissotriton vulgaris species complex based on 74 nuclear markers sampled from 127 localities. We demonstrate that distinct lineages along the speciation continuum in newts exchange nontrivial amounts of genes, affecting their evolutionary trajectories. By integrating a wide array of methods, we delimit nine evolutionary lineages and show that two principal factors have driven their genetic differentiation: time since the last common ancestor determining levels of shared ancestral polymorphism, and shifts in geographic distributions determining the extent of secondary contact. Post-divergence gene flow, indicative of evolutionary non-independence, has been most extensive in Central Europe, while four southern European lineages have acquired the population-genetic hallmarks of independent species (L. graecus, L. kosswigi, L. lantzi, L. schmidtleri). We obtained strong statistical support for widespread mtDNA introgression following secondary contact, previously suggested by discordance between mtDNA phylogeny and morphology. Our study reveals long-term evolutionary persistence of evolutionary lineages that may periodically exchange genes with one another: although some of these lineages may become extinct or fuse, others will acquire complete reproductive isolation and will carry signatures of this complex history in their genomes.

Pan-genome and phylogeny of Bacillus cereus sensu lato

BACKGROUND

Bacillus cereus sensu lato (s. l.) is an ecologically diverse bacterial group of medical and agricultural significance. In this study, I use publicly available genomes and novel bioinformatic workflows to characterize the B. cereus s. l. pan-genome and perform the largest phylogenetic and population genetic analyses of this group to date in terms of the number of genes and taxa included. With these fundamental data in hand, I identify genes associated with particular phenotypic traits (i.e., "pan-GWAS" analysis), and quantify the degree to which taxa sharing common attributes are phylogenetically clustered.

METHODS

A rapid k-mer based approach (Mash) was used to create reduced representations of selected Bacillus genomes, and a fast distance-based phylogenetic analysis of this data (FastME) was performed to determine which species should be included in B. cereus s. l. The complete genomes of eight B. cereus s. l. species were annotated de novo with Prokka, and these annotations were used by Roary to produce the B. cereus s. l. pan-genome. Scoary was used to associate gene presence and absence patterns with various phenotypes. The orthologous protein sequence clusters produced by Roary were filtered and used to build HaMStR databases of gene models that were used in turn to construct phylogenetic data matrices. Phylogenetic analyses used RAxML, DendroPy, ClonalFrameML, PAUP*, and SplitsTree. Bayesian model-based population genetic analysis assigned taxa to clusters using hierBAPS. The genealogical sorting index was used to quantify the phylogenetic clustering of taxa sharing common attributes.

RESULTS

The B. cereus s. l. pan-genome currently consists of ≈60,000 genes, ≈600 of which are "core" (common to at least 99% of taxa sampled). Pan-GWAS analysis revealed genes associated with phenotypes such as isolation source, oxygen requirement, and ability to cause diseases such as anthrax or food poisoning. Extensive phylogenetic analyses using an unprecedented amount of data produced phylogenies that were largely concordant with each other and with previous studies. Phylogenetic support as measured by bootstrap probabilities increased markedly when all suitable pan-genome data was included in phylogenetic analyses, as opposed to when only core genes were used. Bayesian population genetic analysis recommended subdividing the three major clades of B. cereus s. l. into nine clusters. Taxa sharing common traits and species designations exhibited varying degrees of phylogenetic clustering.

CONCLUSIONS

All phylogenetic analyses recapitulated two previously used classification systems, and taxa were consistently assigned to the same major clade and group. By including accessory genes from the pan-genome in the phylogenetic analyses, I produced an exceptionally well-supported phylogeny of 114 complete B. cereus s. l. genomes. The best-performing methods were used to produce a phylogeny of all 498 publicly available B. cereus s. l. genomes, which was in turn used to compare three different classification systems and to test the monophyly status of various B. cereus s. l. species. The majority of the methodology used in this study is generic and could be leveraged to produce pan-genome estimates and similarly robust phylogenetic hypotheses for other bacterial groups.

DNA-barcoding of forensically important blow flies (Diptera: Calliphoridae) in the Caribbean Region

Correct identification of forensically important insects, such as flies in the family Calliphoridae, is a crucial step for them to be used as evidence in legal investigations. Traditional identification based on morphology has been effective, but has some limitations when it comes to identifying immature stages of certain species. DNA-barcoding, using COI, has demonstrated potential for rapid and accurate identification of Calliphoridae, however, this gene does not reliably distinguish among some recently diverged species, raising questions about its use for delimitation of species of forensic importance. To facilitate DNA based identification of Calliphoridae in the Caribbean we developed a vouchered reference collection from across the region, and a DNA sequence database, and further added the nuclear ITS2 as a second marker to increase accuracy of identification through barcoding. We morphologically identified freshly collected specimens, did phylogenetic analyses and employed several species delimitation methods for a total of 468 individuals representing 19 described species. Our results show that combination of COI + ITS2 genes yields more accurate identification and diagnoses, and better agreement with morphological data, than the mitochondrial barcodes alone. All of our results from independent and concatenated trees and most of the species delimitation methods yield considerably higher diversity estimates than the distance based approach and morphology. Molecular data support at least 24 distinct clades within Calliphoridae in this study, recovering substantial geographic variation for Lucilia eximia, Lucilia retroversa, Lucilia rica and Chloroprocta idioidea, probably indicating several cryptic species. In sum, our study demonstrates the importance of employing a second nuclear marker for barcoding analyses and species delimitation of calliphorids, and the power of molecular data in combination with a complete reference database to enable identification of taxonomically and geographically diverse insects of forensic importance.

Speciation history of a species complex of Primulina eburnea (Gesneriaceae) from limestone karsts of southern China, a biodiversity hot spot

Limestone karsts in southern China are characterized by high edaphic and topographic heterogeneity and host high levels of species richness and endemism. However, the evolutionary mechanisms for generating such biodiversity remain poorly understood. Here, we performed species delimitation, population genetic analyses, simulations of gene flow and analyses of floral morphological traits to infer the geographic history of speciation in a species complex of Primulina eburnea from limestone karsts of southern China. Using Bayesian species delimitation, we determined that there are seven distinct species that correspond well to the putative morphological species. Species tree reconstruction, Structure and Neighbor‐Net analyses all recovered four lineages in agreement with currently species geographic boundaries. High levels of genetic differentiation were observed both within and among species. Isolation–migration coalescent analysis provides evidence for significant but low gene flow among species. Approximate Bayesian computation (ABC) analysis supports a scenario of historical gene flow rather than recent contemporary gene flow for most species divergences. Finally, we found no evidence of divergent selection contributing to population differentiation of a suite of flower traits. These results support the prevalence of allopatric speciation and highlight the role of geographic isolation in the diversification process. At small geographic scales, limited hybridization occurred in the past between proximate populations but did not eliminate species boundaries. We conclude that limited gene flow might have been the predominant evolutionary force in promoting population differentiation and speciation.

The natural history of molecular functions inferred from an extensive phylogenomic analysis of gene ontology data

The origin and natural history of molecular functions hold the key to the emergence of cellular organization and modern biochemistry. Here we use a genomic census of Gene Ontology (GO) terms to reconstruct phylogenies at the three highest (1, 2 and 3) and the lowest (terminal) levels of the hierarchy of molecular functions, which reflect the broadest and the most specific GO definitions, respectively. These phylogenies define evolutionary timelines of functional innovation. We analyzed 249 free-living organisms comprising the three superkingdoms of life, Archaea, Bacteria, and Eukarya. Phylogenies indicate catalytic, binding and transport functions were the oldest, suggesting a 'metabolism-first' origin scenario for biochemistry. Metabolism made use of increasingly complicated organic chemistry. Primordial features of ancient molecular functions and functional recruitments were further distilled by studying the oldest child terms of the oldest level 1 GO definitions. Network analyses showed the existence of an hourglass pattern of enzyme recruitment in the molecular functions of the directed acyclic graph of molecular functions. Older high-level molecular functions were thoroughly recruited at younger lower levels, while very young high-level functions were used throughout the timeline. This pattern repeated in every one of the three mappings, which gave a criss-cross pattern. The timelines and their mappings were remarkable. They revealed the progressive evolutionary development of functional toolkits, starting with the early rise of metabolic activities, followed chronologically by the rise of macromolecular biosynthesis, the establishment of controlled interactions with the environment and self, adaptation to oxygen, and enzyme coordinated regulation, and ending with the rise of structural and cellular complexity. This historical account holds important clues for dissection of the emergence of biomcomplexity and life.

Transitions between phases of genomic differentiation during stick-insect speciation

Speciation can involve a transition from a few genetic loci that are resistant to gene flow to genome-wide differentiation. However, only limited data exist concerning this transition and the factors promoting it. Here, we study phases of speciation using data from >100 populations of 11 species of Timema stick insects. Consistent with early phases of genic speciation, adaptive colour-pattern loci reside in localized genetic regions of accentuated differentiation between populations experiencing gene flow. Transitions to genome-wide differentiation are also observed with gene flow, in association with differentiation in polygenic chemical traits affecting mate choice. Thus, intermediate phases of speciation are associated with genome-wide differentiation and mate choice, but not growth of a few genomic islands. We also find a gap in genomic differentiation between sympatric taxa that still exchange genes and those that do not, highlighting the association between differentiation and complete reproductive isolation. Our results suggest that substantial progress towards speciation may involve the alignment of multi-faceted aspects of differentiation.

Taxonomic position of Chinese voles of the tribe Arvicolini and the description of 2 new species from Xizang, China

China has 26 species in the tribe Arvicolini. The taxonomic status of these voles remains controversial despite much effort. Herein, we evaluate the taxonomic position of 22 species plus 2 unidentified taxa using mitochondrial DNA gene sequences (cytb + CO1). We also evaluate 18 species and 2 unidentified taxa using morphological data. Phylogenetic analyses of cytb resolve monophyly for the genera Alexandromys, Lasiopodomys, Microtus, Neodon, Proedromys, and Volemys with strong support. Stenocranius clusters with Chionomys but with very weak support. Analyses of concatenated cytb + CO1 resolve the same genera with strong support, but the topology of the tree differs from that of cytb in that Chionomys roots at the base of the tree independent of Stenocranius, which forms the sister-group of Lasiopodomys in a more terminal position. The matrilineal genealogy excludes the type species Arvicola amphibius from the rest of the Arvicolini. This species forms the sister-group of Ondatra with high support. Neodon includes N. irene, N. linzhiensis, N. fuscus, N. leucurus, N. sikimensis, Microtus clarkei, and 2 unidentified specimens. Alexandromys includes the former species Microtus oeconomus, M. kikuchii, M. limnophilus, M. fortis, and M. maximowiczii. Finally, Microtus has the subgenera Blanfordimys, Microtus, Mynomes, Pedomys, Pitymys, and Terricola, which includes the Chinese species M. agrestis, M. arvalis, and Blanfordimys juldaschi. General mixed Yule-coalescent species delimitation modeling demarcates 6 currently recognized species and 2 new species of Neodon. A principal component analysis of the morphological data among 7 matrilines shows that all variables have positive loadings of high magnitude on the 1st component. Canonical discriminant analysis for Neodon (including M. clarkei and 2 unidentified species) correctly classifies 93.0% of specimens. Overall, our analyses support the recognition of Alexandromys, Lasiopodomys, Microtus, Neodon, Proedromys, and Volemys as genera. Stenocranius includes Microtus gregalis, and the genealogical position of Stenocranius remains uncertain. The status of Arvicola requires further study. We assign M. clarkei to Neodon and describe 2 new species of Neodon.

Catchments catch all in South African coastal lowlands: topography and palaeoclimate restricted gene flow in Nymania capensis (Meliaceae)-a multilocus phylogeographic and distribution modelling approach

BACKGROUND

This study investigates orbitally-forced range dynamics at a regional scale by exploring the evolutionary history of Nymania capensis (Meliaceae) across the deeply incised landscapes of the subescarpment coastal lowlands of South Africa; a region that is home to three biodiversity hotspots (Succulent Karoo, Fynbos, and Maputaland-Pondoland-Albany hotspots).

METHODS

A range of methods are used including: multilocus phylogeography (chloroplast and high- and low-copy nuclear DNA), molecular dating and species distribution modelling (SDM).

RESULTS

The results support an 'evolutionarily distinct catchment' hypothesis where: (1) different catchments contain genetically distinct lineages, (2) limited genetic structuring was detected within basins whilst high structuring was detected between basins, and (3) within primary catchment populations display a high degree of genealogical lineage sorting. In addition, the results support a glacial refugia hypothesis as: (a) the timing of chloroplast lineage diversification is restricted to the Pleistocene in a landscape that has been relatively unchanged since the late Pliocene, and (b) the projected LGM distribution of suitable climate for N. capensis suggest fragmentation into refugia that correspond to the current phylogeographic populations.

DISCUSSION

This study highlights the interaction of topography and subtle Pleistocene climate variations as drivers limiting both seed and pollen flow along these lowlands. This lends support to the region's large-scale conservation planning efforts, which used catchments as foundational units for conservation as these are likely to be evolutionarily significant units.

Resistance to type 1 interferons is a major determinant of HIV-1 transmission fitness

Sexual transmission of HIV-1 is an inefficient process, with only one or few variants of the donor quasispecies establishing the new infection. A critical, and as yet unresolved, question is whether the mucosal bottleneck selects for viruses with increased transmission fitness. Here, we characterized 300 limiting dilution-derived virus isolates from the plasma, and in some instances genital secretions, of eight HIV-1 donor and recipient pairs. Although there were no differences in the amount of virion-associated envelope glycoprotein, recipient isolates were on average threefold more infectious (P = 0.0001), replicated to 1.4-fold higher titers (P = 0.004), were released from infected cells 4.2-fold more efficiently (P < 0.00001), and were significantly more resistant to type I IFNs than the corresponding donor isolates. Remarkably, transmitted viruses exhibited 7.8-fold higher IFNα2 (P < 0.00001) and 39-fold higher IFNβ (P < 0.00001) half-maximal inhibitory concentrations (IC50) than did donor isolates, and their odds of replicating in CD4+ T cells at the highest IFNα2 and IFNβ doses were 35-fold (P < 0.00001) and 250-fold (P < 0.00001) greater, respectively. Interestingly, pretreatment of CD4+ T cells with IFNβ, but not IFNα2, selected donor plasma isolates that exhibited a transmitted virus-like phenotype, and such viruses were also detected in the donor genital tract. These data indicate that transmitted viruses are phenotypically distinct, and that increased IFN resistance represents their most distinguishing property. Thus, the mucosal bottleneck selects for viruses that are able to replicate and spread efficiently in the face of a potent innate immune response.

Phylogeny and Haplotype Analysis of Fungi Within the Fusarium incarnatum-equiseti Species Complex

Fusarium spp. are ranked among the top 10 most economically and scientifically important plant-pathogenic fungi in the world and are associated with plant diseases that include fruit decay of a number of crops. Fusarium isolates infecting bell pepper in Trinidad were identified based on sequence comparisons of the translation elongation factor gene (EF-1a) with sequences of Fusarium incarnatum-equiseti species complex (FIESC) verified in the FUSARIUM-ID database. Eighty-two isolates were identified as belonging to one of four phylogenetic species within the subclades FIESC-1, FIESC-15, FIESC-16, and FIESC-26, with the majority of isolates belonging to FIESC-15. A comparison of the level of DNA polymorphism and phylogenetic inference for sequences of the internal transcribed spacer region (ITS1-5.8S-ITS2) and EF-1a sequences for Trinidad and FUSARIUM-ID type species was carried out. The ITS sequences were less informative, had lower haplotype diversity and restricted haplotype distribution, and resulted in poor resolution and taxa placement in the consensus maximum-likelihood tree. EF-1a sequences enabled strongly supported phylogenetic inference with highly resolved branching patterns of the 30 phylogenetic species within the FIESC and placement of representative Trinidad isolates. Therefore, global phylogeny was inferred from EF-1a sequences representing 11 countries, and separation into distinct Incarnatum and Equiseti clades was again evident. In total, 42 haplotypes were identified: 12 were shared and the remaining were unique haplotypes. The most diverse haplotype was represented by sequences from China, Indonesia, Malaysia, and Trinidad and consisted exclusively of F. incarnatum isolates. Spain had the highest haplotype diversity, perhaps because both F. equiseti and F. incarnatum sequences were represented; followed by the United States, which contributed both F. equiseti and F. incarnatum sequences to the data set; then by countries representing Southeast Asia (China, Indonesia, Malaysia, Thailand, and Philippines) and Trinidad; both of these regions were represented by only F. incarnatum sequences. Trinidad shared two haplotypes with China and one haplotype with the United States for only F. incarnatum isolates. The findings of this study are important for devising disease management strategies and for understanding the phylogenetic relationships among members of the FIESC.

Sharpening the species boundaries in the Cladonia mediterranea complex (Cladoniaceae, Ascomycota)

The complex Cladonia mediterranea belongs to the section Impexae and is formed by C. azorica, C. macaronesica and C.mediterranea. These species are basically distributed in the Mediterranean and Macaronesian Regions. In the present work the limits between the species of this complex are re-examined. To this end, the morphological characters were studied along with the secondary metabolites and the DNA sequences from three loci (ITS rDNA, IGS rDNA and rpb2). The morphological data were studied by principal component analysis (PCA), while the DNA sequences were analyzed using several approaches available to delimit species: genealogical concordance phylogenetic species recognition, species tree (BEAST* and spedeSTEM) and cohesion species recognition. In addition, the genealogical sorting index was used in order to assess the monophyly of the species. The different procedures used in our study turned out to be highly congruent with respect to the limits they establish, but these limits are not the ones separating the prior species. Either the morphological analysis or the different approaches to species delimitation indicate that C. mediterranea is a different species from C. macaronesica, while C. azorica and C. macaronesica, which are reduced to synonyms of C. portentosa, constitute a separate lineage.

Paired quantitative and qualitative assessment of the replication-competent HIV-1 reservoir and comparison with integrated proviral DNA

HIV-1-infected individuals harbor a latent reservoir of infected CD4⁺ T cells that is not eradicated by antiretroviral therapy (ART). This reservoir presents the greatest barrier to an HIV-1 cure and has remained difficult to characterize, in part, because the vast majority of integrated sequences are defective and incapable of reactivation. To characterize the replication-competent reservoir, we have combined two techniques, quantitative viral outgrowth and qualitative sequence analysis of clonal outgrowth viruses. Leukapheresis samples from four fully ART-suppressed, chronically infected individuals were assayed at two time points separated by a 4- to 6-mo interval. Overall, 54% of the viruses emerging from the latent reservoir showed gp160 env sequences that were identical to at least one other virus. Moreover, 43% of the env sequences from viruses emerging from the reservoir were part of identical groups at the two time points. Groups of identical expanded sequences made up 54% of proviral DNA, and, as might be expected, the sequences of replication-competent viruses in the active reservoir showed limited overlap with integrated proviral DNA, most of which is known to represent defective viruses. Finally, there was an inverse correlation between proviral DNA clone size and the probability of reactivation, suggesting that replication-competent viruses are less likely to be found among highly expanded provirus-containing cell clones.

Diversification in wild populations of the model organism Anolis carolinensis: A genome-wide phylogeographic investigation

The green anole (Anolis carolinensis) is a lizard widespread throughout the southeastern United States and is a model organism for the study of reproductive behavior, physiology, neural biology, and genomics. Previous phylogeographic studies of A. carolinensis using mitochondrial DNA and small numbers of nuclear loci identified conflicting and poorly supported relationships among geographically structured clades; these inconsistencies preclude confident use of A. carolinensis evolutionary history in association with morphological, physiological, or reproductive biology studies among sampling localities and necessitate increased effort to resolve evolutionary relationships among natural populations. Here, we used anchored hybrid enrichment of hundreds of genetic markers across the genome of A. carolinensis and identified five strongly supported phylogeographic groups. Using multiple analyses, we produced a fully resolved species tree, investigated relative support for each lineage across all gene trees, and identified mito‐nuclear discordance when comparing our results to previous studies. We found fixed differences in only one clade—southern Florida restricted to the Everglades region—while most polymorphisms were shared between lineages. The southern Florida group likely diverged from other populations during the Pliocene, with all other diversification during the Pleistocene. Multiple lines of support, including phylogenetic relationships, a latitudinal gradient in genetic diversity, and relatively more stable long‐term population sizes in southern phylogeographic groups, indicate that diversification in A. carolinensis occurred northward from southern Florida.

Genomic and Evolutionary Analysis of Two Salmonella enterica Serovar Kentucky Sequence Types Isolated from Bovine and Poultry Sources in North America

Salmonella enterica subsp. enterica serovar Kentucky is frequently isolated from healthy poultry and dairy cows and is occasionally isolated from people with clinical disease. A genomic analysis of 119 isolates collected in the United States from dairy cows, ground beef, poultry and poultry products, and human clinical cases was conducted. Results of the analysis demonstrated that the majority of poultry and bovine-associated S. Kentucky were sequence type (ST) 152. Several bovine-associated (n = 3) and food product isolates (n = 3) collected from the United States and the majority of human clinical isolates were ST198, a sequence type that is frequently isolated from poultry and occasionally from human clinical cases in Northern Africa, Europe and Southeast Asia. A phylogenetic analysis indicated that both STs are more closely related to other Salmonella serovars than they are to each other. Additionally, there was strong evidence of an evolutionary divergence between the poultry-associated and bovine-associated ST152 isolates that was due to polymorphisms in four core genome genes. The ST198 isolates recovered from dairy farms in the United States were phylogenetically distinct from those collected from human clinical cases with 66 core genome SNPs differentiating the two groups, but more isolates are needed to determine the significance of this distinction. Identification of S. Kentucky ST198 from dairy animals in the United States suggests that the presence of this pathogen should be monitored in food-producing animals.

High regional genetic diversity and lack of host-specificity in Ostrinia nubilalis (Lepidoptera: Crambidae) as revealed by mtDNA variation

The European corn borer (Ostrinia nubilalis) infests a wide array of host plants and is considered one of the most serious pests of maize in Europe. Recent studies suggest that individuals feeding on maize in Europe should be referred to O. nubilalis (sensu nov.), while those infesting dicots as Ostrinia scapulalis (sensu nov.). We test if the clear genetic distinctiveness among individuals of O. nubilalis living on maize vs. dicots is tracked by mitochondrial DNA (mtDNA). We used fragments of COI and COII genes of 32 individuals traditionally recognized as O. nubilalis collected on three host plants, maize, mugwort and hop, growing in different parts of Poland. In addition, we reconstructed the mtDNA phylogeny of Ostrinia species based on our data and sequences retrieved from GenBank to assess host and/or biogeographic patterns. We also compared haplotype variation found in Poland (east-central Europe) with other regions (Anatolia, Eastern Europe, Balkans, Far East, North America). Our study showed high mtDNA diversity of O. nubilalis in Poland in comparison with other regions and revealed rare haplotypes likely of Asian origin. We did not find distinct mtDNA haplotypes in larvae feeding on maize vs. dicotyledonous plants. Phylogenetic analyses showed an apparent lack of mtDNA divergence among putatively distinct lineages belonging to the O. nubilalis group as identical haplotypes are shared by Asian and European individuals. We argue that human-mediated dispersal, hybridization and sporadic host jumps are likely responsible for the lack of a geographic pattern in mtDNA variation.

Origin of a cryptic lineage in a threatened reptile through isolation and historical hybridization

Identifying phylogenetically distinct lineages and understanding the evolutionary processes by which they have arisen are important goals of phylogeography. This information can also help define conservation units in endangered species. Such analyses are being transformed by the availability of genomic-scale data sets and novel analytical approaches for statistically comparing different historical scenarios as causes of phylogeographic patterns. Here, we use genomic-scale restriction-site-associated DNA sequencing (RADseq) data to test for distinct lineages in the endangered Eastern Massasauga Rattlesnake (Sistrurus catenatus). We then use coalescent-based modeling techniques to identify the evolutionary mechanisms responsible for the origin of the lineages in this species. We find equivocal evidence for distinct phylogenetic lineages within S. catenatus east of the Mississippi River, but strong support for a previously unrecognized lineage on the western edge of the range of this snake, represented by populations from Iowa, USA. Snakes from these populations show patterns of genetic admixture with a nearby non-threatened sister species (Sistrurus tergeminus). Tests of historical demographic models support the hypothesis that the genetic distinctiveness of Iowa snakes is due to a combination of isolation and historical introgression between S. catenatus and S. tergeminus. Our work provides an example of how model-based analysis of genomic-scale data can help identify conservation units in rare species.

Caribbean Bulimulus revisited: physical moves and molecular traces (Mollusca, Gastropoda, Bulimulidae)

Twenty-five samples of Bulimulus species are studied, partly from localities within their known distribution range, partly based on interceptions where the material originates from localities where the species seem to be recently introduced and non-native. Molecular study of cytochrome oxidase 1 (CO1) reveals the origin of some of these introductions, but is less conclusive for others. Four different methods for species delimitation were applied, which did not result in unambiguous species hypotheses. For a rapid identification of morphologically indistinct species, a more comprehensive database of sequences is needed.