Discerning structure versus speciation in phylogeographic analysis of Seepage Salamanders (Desmognathus aeneus) using demography, environment, geography, and phenotype

Numerous mechanisms can drive speciation, including isolation by adaptation, distance, and environment. These forces can promote genetic and phenotypic differentiation of local populations, the formation of phylogeographic lineages, and ultimately, completed speciation. However, conceptually similar mechanisms may also result in stabilizing rather than diversifying selection, leading to lineage integration and the long-term persistence of population structure within genetically cohesive species. Processes that drive the formation and maintenance of geographic genetic diversity while facilitating high rates of migration and limiting phenotypic differentiation may thereby result in population genetic structure that is not accompanied by reproductive isolation. We suggest that this framework can be applied more broadly to address the classic dilemma of "structure" versus "species" when evaluating phylogeographic diversity, unifying population genetics, species delimitation, and the underlying study of speciation. We demonstrate one such instance in the Seepage Salamander (Desmognathus aeneus) from the southeastern United States. Recent studies estimated up to 6.3% mitochondrial divergence and four phylogenomic lineages with broad admixture across geographic hybrid zones, which could potentially represent distinct species supported by our species-delimitation analyses. However, while limited dispersal promotes substantial isolation by distance, microhabitat specificity appears to yield stabilizing selection on a single, uniform, ecologically mediated phenotype. As a result, climatic cycles promote recurrent contact between lineages and repeated instances of high migration through time. Subsequent hybridization is apparently not counteracted by adaptive differentiation limiting introgression, leaving a single unified species with deeply divergent phylogeographic lineages that nonetheless do not appear to represent incipient species.

Phylogeographic reconstructions can be biased by ancestral shared alleles: The case of the polymorphic lichen Bryoria fuscescens in Europe and North Africa

Large phylogeographic studies on lichens are scarce, and none involves a single species within which different lineages show fixed alternative dispersal strategies. We investigated Bryoria fuscescens (including B. capillaris) in Europe and western North Africa by phenotypically characterizing 1400 specimens from 64 populations and genotyping them with 14 microsatellites. We studied population structure and genetic diversity at the local and continental scales, discussed the post-glacial phylogeography, and compared dispersal capacities of phenotypes with and without soralia. Our main hypothesis is that the estimated phylogeography, migration routes, and dispersal capacities may be strongly biased by ancestral shared alleles. Scandinavia is genetically the richest area, followed by the Iberian Peninsula, the Carpathians, and the Alps. Three gene pools were detected: two partially linked to phenotypic characteristics, and the third one genetically related to the American sister species B. pseudofuscescens. The comparison of one gene pool producing soredia and one not, suggested both as panmictic, with similar levels of isolation by distance (IBD). The migration routes were estimated to span from north to south, in disagreement with the assessed glacial refugia. The presence of ancestral shared alleles in distant populations can explain the similar IBD levels found in both gene pools while producing a false signal of panmixia, and also biasing the phylogeographic reconstruction. The incomplete lineage sorting recorded for DNA sequence loci also supports this hypothesis. Consequently, the high diversity in Scandinavia may rather come from recent immigration into northern populations than from an in situ diversification. Similar patterns of ancestral shared polymorphism may bias the phylogeographical reconstruction of other lichen species.

Population connectivity and larval dispersal of the exploited mangrove crab Ucides cordatus along the Brazilian coast

BACKGROUND

The mangrove crab Ucides cordatus is considered a key species for the ecological balance of mangrove forests and a major source of employment and income for traditional crab collectors in Brazil. Several studies evidenced weak genetic variation among populations due to an efficient larval transport. However, gene flow patterns of the species is poorly understood, with no information about migration rates. The influence of the two main Brazilian currents in larval dispersion is also not clear. In order to provide baseline information for conservation, planning and management of this important fishery resource, the present study aimed to estimate and evaluate spatial distribution of genetic diversity, migration rates and gene flow directivity among populations of U. cordatus in Brazil.

METHODS

Nine microsatellites were used to resolve population structure of 319 crabs collected from six sites located along the Brazilian coast. The degree of geographical differentiation included estimates of genetic diversity, population structure and gene flow models, with spatial analysis of shared alleles (SAShA), isolation by distance tests, AMOVA, discriminant analysis of principal components (DAPC) and Bayesian clustering. We estimated the amount of ongoing gene flow between clusters using the coalescent-based method implemented in Migrate-N.

RESULTS

Loci were highly polymorphic (average of 12.4 alleles per locus) evidencing high genetic variability. There was significant differentiation among localities, despite of the low value of FST (= 0.019; P < 0.001). FST and Jost's D indexes were also estimated in pairwise comparisons and showed significant differences between most of the surveyed site pairs (P < 0.05). Structure evidenced a single genetic group among samples, however SAShA pointed to a non-panmictic condition (P = 0.011). AMOVA detected four statistical significant clusters with low level of differentiation (FCT = 0.037; P = 0.023). The gene flow model that best described the population connectivity was the island model, with ∼24 crabs being exchanged among localities per generation.

DISCUSSION

The high migration rates found among localities seem to be the main force acting to sustain the distribution of the genetic diversity of U. cordatus. Despite the high gene flow and the weak population structure among samples, the significant genetic differences found suggest that gene flow alone does not bypass the effects of genetic drift, natural selection and/or human exploitation. These findings are vital for the establishment of a database to be used in the development of conservation programs.

Evidence for rangewide panmixia despite multiple barriers to dispersal in a marine mussel

Oceanographic features shape the distributional and genetic patterns of marine species by interrupting or promoting connections among populations. Although general patterns commonly arise, distributional ranges and genetic structure are species-specific and do not always comply with the expected trends. By applying a multimarker genetic approach combined with Lagrangian particle simulations (LPS) we tested the hypothesis that oceanographic features along northeastern Atlantic and Mediterranean shores influence dispersal potential and genetic structure of the intertidal mussel Perna perna. Additionally, by performing environmental niche modelling we assessed the potential and realized niche of P. perna along its entire native distributional range and the environmental factors that best explain its realized distribution. Perna perna showed evidence of panmixia across >4,000 km despite several oceanographic breaking points detected by LPS. This is probably the result of a combination of life history traits, continuous habitat availability and stepping-stone dynamics. Moreover, the niche modelling framework depicted minimum sea surface temperatures (SST) as the major factor shaping P. perna distributional range limits along its native areas. Forthcoming warming SST is expected to further change these limits and allow the species to expand its range polewards though this may be accompanied by retreat from warmer areas.

Asymmetrical dispersal and putative isolation-by-distance of an intertidal blenniid across the Atlantic-Mediterranean divide

Transition zones are of high evolutionary interest because unique patterns of spatial variation are often retained. Here, we investigated the phylogeography of the peacock blenny, Salaria pavo, a small marine intertidal fish that inhabits rocky habitats of the Mediterranean and the adjacent Atlantic Ocean. We screened 170 individuals using mitochondrial and nuclear sequence data from eight locations. Four models of genetic structure were tested: panmixia, isolation-by-distance, secondary contact and phylogeographic break. Results indicated clear asymmetric migration from the Mediterranean to the Atlantic but only marginally supported the isolation-by-distance model. Additionally, the species displays an imprint of demographic expansion compatible with the last glacial maximum. Although the existence of a refugium in the Mediterranean cannot be discarded, the ancestral lineage most likely originated in the Atlantic, where most of the genetic diversity occurs.

Panmixia and dispersal from the Mediterranean Basin to Macaronesian Islands of a macrolichen species

The Mediterranean region, comprising the Mediterranean Basin and the Macaronesian Islands, represents a center of diversification for many organisms. The genetic structure and connectivity of mainland and island microbial populations has been poorly explored, in particular in the case of symbiotic fungi. Here we investigated genetic diversity and spatial structure of the obligate outcrossing lichen-forming fungus Parmelina carporrhizans in the Mediterranean region. Using eight microsatellite and mating-type markers we showed that fungal populations are highly diverse but lack spatial structure. This is likely due to high connectivity and long distance dispersal of fungal spores. Consistent with low levels of linkage disequilibrium and lack of clonality, we detected both mating-type idiomorphs in all populations. Furthermore we showed that the Macaronesian Islands are the result of colonization from the Mediterranean Basin. The unidirectional gene flow, though, seemed not to be sufficient to counterbalance the effects of drift, resulting in comparatively allelic poor peripheral populations. Our study is the first to shed light on the high connectivity and lack of population structure in natural populations of a strictly sexual lichen fungus. Our data further support the view of the Macaronesian Islands as the end of the colonization road for this symbiotic ascomycete.

Genomic evidence for ecological divergence against a background of population homogeneity in the marine snail Chlorostoma funebralis

The balance between natural selection, gene flow and genetic drift is difficult to resolve in marine invertebrates with extensive dispersal and fluctuating population sizes. The intertidal snail Chlorostoma funebralis has planktonic larvae and previous work using mtDNA polymorphism reported no genetic population structure. Nevertheless, recent studies have documented differences in thermal tolerance and transcriptomic responses to heat stress between northern and southern California, USA, populations. To gain insight into the dynamics influencing adaptive divergence, we used double-digest restriction site-associated DNA (ddRAD) sequencing to identify 1861 genomewide, quality-filtered single-nucleotide polymorphism (SNP) loci for C. funebralis collected from three northern and three southern California sites (15 individuals per population). Considering all SNPs, there was no evidence for genetic differentiation among populations or regions (average FST  = 0.0042). However, outlier tests revealed 34 loci putatively under divergent selection between northern and southern populations, and structure and SNP tree analyses based on these outliers show clear genetic differentiation between geographic regions. Three of these outliers are known or hypothesized to be involved in stress granule formation, a response to environmental stress such as heat. Combined with previous work that found thermally tolerant southern populations show high baseline expression of stress response genes, these results further suggest that thermal stress is a strong selective pressure across C. funebralis populations. Overall, this study increases our understanding of the factors constraining local adaptation in marine organisms, while suggesting that ecologically driven, strong differentiation can occur at relevant loci in a species with planktonic larvae.

Stochasticity in space, persistence in time: genetic heterogeneity in harbour populations of the introduced ascidian Styela plicata

Spatio-temporal changes in genetic structure among populations provide crucial information on the dynamics of secondary spread for introduced marine species. However, temporal components have rarely been taken into consideration when studying the population genetics of non-indigenous species. This study analysed the genetic structure of Styela plicata, a solitary ascidian introduced in harbours and marinas of tropical and temperate waters, across spatial and temporal scales. A fragment of the mitochondrial gene Cytochrome Oxidase subunit I (COI) was sequenced from 395 individuals collected at 9 harbours along the NW Mediterranean coast and adjacent Atlantic waters (> 1,200 km range) at two time points 5 years apart (2009 and 2014). The levels of gene diversity were relatively low for all 9 locations in both years. Analyses of genetic differentiation and distribution of molecular variance revealed strong genetic structure, with significant differences among many populations, but no significant differences among years. A weak and marginally significant correlation between geographic distance and gene differentiation was found. Our results revealed spatial structure and temporal genetic homogeneity in S. plicata, suggesting a limited role of recurrent, vessel-mediated transport of organisms among small to medium-size harbours. Our study area is representative of many highly urbanized coasts with dense harbours. In these environments, the episodic chance arrival of colonisers appears to determine the genetic structure of harbour populations and the genetic composition of these early colonising individuals persists in the respective harbours, at least over moderate time frames (five years) that encompass ca. 20 generations of S. plicata.

Three in One--Multiple Faunal Elements within an Endangered European Butterfly Species

Ice ages within Europe forced many species to retreat to refugia, of which three major biogeographic basic types can be distinguished: "Mediterranean", "Continental" and "Alpine / Arctic" species. However, this classification often fails to explain the complex phylogeography of European species with a wide range of latitudinal and altitudinal distribution. Hence, we tested for the possibility that all three mentioned faunal elements are represented within one species. Our data was obtained by scoring 1,307 Euphydryas aurinia individuals (46 European locations) for 17 allozyme loci, and sequencing a subset of 492 individuals (21 sites) for a 626 base pairs COI fragment. Genetic diversity indices, F statistics, hierarchical analyses of molecular variance, individual-based clustering, and networks were used to explore the phylogeographic patterns. The COI fragment represented 18 haplotypes showing a strong geographic structure. All but one allozyme loci analysed were polymorphic with a mean F_ST of 0.20, supporting a pronounced among population structure. Interpretation of both genetic marker systems, using several analytical tools, calls for the recognition of twelve genetic groups. These analyses consistently distinguished different groups in Iberia (2), Italy, Provence, Alps (3), Slovenia, Carpathian Basin, the lowlands of West and Central Europe as well as Estonia, often with considerable additional substructures. The genetic data strongly support the hypothesis that E. aurinia survived the last glaciation in Mediterranean, extra-Mediterranean and perialpine refugia. It is thus a rare example of a model organism that combines attributes of faunal elements from all three of these sources. The observed differences between allozymes and mtDNA most likely result from recent introgression of mtDNA into nuclear allozyme groups. Our results indicate discrepancies with the morphologically-based subspecies models, underlining the need to revise the current taxonomy.

Three-dimensional post-glacial expansion and diversification of an exploited oceanic fish

Vertical divergence in marine organisms is being increasingly documented, yet much remains to be carried out to understand the role of depth in the context of phylogeographic reconstruction and the identification of management units. An ideal study system to address this issue is the beaked redfish, Sebastes mentella – one of four species of ‘redfish’ occurring in the North Atlantic – which is known for a widely distributed ‘shallow‐pelagic’ oceanic type inhabiting waters between 250 and 550 m, and a more localized ‘deep‐pelagic’ population dwelling between 550 and 800 m, in the oceanic habitat of the Irminger Sea. Here, we investigate the extent of population structure in relation to both depth and geographic spread of oceanic beaked redfish throughout most of its distribution range. By sequencing the mitochondrial control region of 261 redfish collected over a decadal interval, and combining 160 rhodopsin coding nuclear sequences and previously genotyped microsatellite data, we map the existence of two strongly divergent evolutionary lineages with significantly different distribution patterns and historical demography, and whose genetic variance is mostly explained by depth. Combined genetic data, analysed via independent approaches, are consistent with a Late Pleistocene lineage split, where segregation by depth probably resulted from the interplay of climatic and oceanographic processes with life history and behavioural traits. The ongoing process of diversification in North Atlantic S. mentella may serve as an ‘hourglass’ to understand speciation and adaptive radiation in Sebastes and in other marine taxa distributed across a depth gradient.

Spatio-temporal Genetic Structure of a Tropical Bee Species Suggests High Dispersal Over a Fragmented Landscape

Habitat destruction threatens biodiversity by reducing the amount of available resources and connectivity among geographic areas. For organisms living in fragmented habitats, population persistence may depend on dispersal, which maintains gene flow among fragments and can prevent inbreeding within them. It is centrally important to understand patterns of dispersal for bees living in fragmented areas given the importance of pollination systems and recently documented declines in bee populations. We used population and landscape genetic techniques to characterize patterns of dispersal over a large fragmented area in southern Costa Rica for the orchid bee species Euglossa championi. First, we estimated levels of genetic differentiation among forest fragments as φ_pt, an analog to the traditional summary statistic F_st, as well as two statistics that may more adequately represent levels of differentiation, G'_st and D_est . Second, we used a Bayesian approach to determine the number and composition of genetic groups in our sample. Third we investigated how genetic differentiation changes with distance. Fourth, we determined the extent to which deforested areas restrict dispersal. Finally, we estimated the extent to which there were temporal differences in allele frequencies within the same forest fragments. Within years we found low levels of differentiation even over 80 km, and no effect of land use type on level of genetic differentiation. However, we found significant genetic differentiation between years. Taken together our results suggest that there are high levels of gene flow over this geographic area, and that individuals show low site fidelity over time.

Is population structure in the European white stork determined by flyway permeability rather than translocation history?

European white stork are long considered to diverge to eastern and western migration pools as a result of independent overwintering flyways. In relatively recent times, the western and northern distribution has been subject to dramatic population declines and country-specific extirpations. A number of independent reintroduction programs were started in the mid 1950s to bring storks back to historical ranges. Founder individuals were sourced opportunistically from the Eastern and Western European distributions and Algeria, leading to significant artificial mixing between eastern and western flyways. Here we use mitochondrial and microsatellite DNA to test the contention that prior to translocation, eastern and western flyways were genetically distinct. The data show a surprising lack of structure at any spatial or temporal scale suggesting that even though birds were moved between flyways, there is evidence of natural mixing prior to the onset of translocation activities. Overall a high retention of genetic diversity, high N_ef, and an apparent absence of recent genetic bottleneck associated with early 20th century declines suggest that the species is well equipped to respond to future environmental pressures.

Phylogeography of the California sheephead, Semicossyphus pulcher: the role of deep reefs as stepping stones and pathways to antitropicality

In the past decade, the study of dispersal of marine organisms has shifted from focusing predominantly on the larval stage to a recent interest in adult movement. Antitropical distributions provide a unique system to assess vagility and dispersal. In this study, we have focused on an antitropical wrasse genus, Semicossyphus, which includes the California sheephead, S. pulcher, and Darwin's sheephead, S. darwini. Using a phylogenetic approach based on mitochondrial and nuclear markers, and a population genetic approach based on mitochondrial control region sequences and 10 microsatellite loci, we compared the phylogenetic relationships of these two species, as well as the population genetic characteristics within S. pulcher. While S. pulcher and S. darwini are found in the temperate eastern Pacific regions of the northern and southern hemispheres, respectively, their genetic divergence was very small (estimated to have occurred between 200 and 600 kya). Within S. pulcher, genetic structuring was generally weak, especially along mainland California, but showed weak differentiation between Sea of Cortez and California, and between mainland California and Channel Islands. We highlight the congruence of weak genetic differentiation both within and between species and discuss possible causes for maintenance of high gene flow. In particular, we argue that deep and cooler water refugia are used as stepping stones to connect distant populations, resulting in low levels of genetic differentiation.

Optimizing the trade-off between spatial and genetic sampling efforts in patchy populations: towards a better assessment of functional connectivity using an individual-based sampling scheme

Genetic data are increasingly used in landscape ecology for the indirect assessment of functional connectivity, that is, the permeability of landscape to movements of organisms. Among available tools, matrix correlation analyses (e.g. Mantel tests or mixed models) are commonly used to test for the relationship between pairwise genetic distances and movement costs incurred by dispersing individuals. When organisms are spatially clustered, a population-based sampling scheme (PSS) is usually performed, so that a large number of genotypes can be used to compute pairwise genetic distances on the basis of allelic frequencies. Because of financial constraints, this kind of sampling scheme implies a drastic reduction in the number of sampled aggregates, thereby reducing sampling coverage at the landscape level. We used matrix correlation analyses on simulated and empirical genetic data sets to investigate the efficiency of an individual-based sampling scheme (ISS) in detecting isolation-by-distance and isolation-by-barrier patterns. Provided that pseudo-replication issues are taken into account (e.g. through restricted permutations in Mantel tests), we showed that the use of interindividual measures of genotypic dissimilarity may efficiently replace interpopulation measures of genetic differentiation: the sampling of only three or four individuals per aggregate may be sufficient to efficiently detect specific genetic patterns in most situations. The ISS proved to be a promising methodological alternative to the more conventional PSS, offering much flexibility in the spatial design of sampling schemes and ensuring an optimal representativeness of landscape heterogeneity in data, with few aggregates left unsampled. Each strategy offering specific advantages, a combined use of both sampling schemes is discussed.

Dispersal at a snail's pace: historical processes affect contemporary genetic structure in the exploited wavy top snail (Megastraea undosa)

We used population genetics to assess historical and modern demography of the exploited wavy top snail, Megastraea undosa, which has a 5-10 day pelagic larval duration. Foot tissue was sampled from an average of 51 individuals at 17 sites across the range of M. undosa. Genetic structure at the mtDNA locus is strikingly high (ΦST of 0.19 across 1000 km), and a major cline occurs in northern Baja California (ΦCT of 0.29 between northern and southern populations). Genetic data indicate that the northern region is highly connected through larval dispersal, whereas the southern region exhibits low genetic structure. However, additional analyses based on patterns of haplotype diversity and relationships among haplotypes indicate that M. undosa has likely recently expanded into the Southern California Bight or expanded from a small refugial population, and analysis using isolation by distance to calculate dispersal distance indicates surprisingly short estimates of dispersal from 30 m to 3 km. This scenario of a northward expansion and limited larval dispersal is supported by coalescent-based simulations of genetic data. The different patterns of genetic variation between northern and southern populations are likely artifacts of evolutionary history rather than differences in larval dispersal and this may have applications to management of this species. Specifically, these data can help to inform the scale at which this species should be managed, and given the potentially very small dispersal distances, this species should be managed at local scales. Consideration of the evolutionary history of target species allows for a more accurate interpretation of genetic data for management.

Chiton phylogeny (Mollusca : Polyplacophora) and the placement of the enigmatic species Choriplax grayi (H. Adams & Angas)

Shallow marine chitons (Mollusca : Polyplacophora : Chitonida) are widespread and well described from established morphoanatomical characters, yet key aspects of polyplacophoran phylogeny have remained unresolved. Several species, including Hemiarthrum setulosum Carpenter in Dall, 1876, and especially the rare and enigmatic Choriplax grayi (Adams & Angas, 1864), defy systematic placement. Choriplax is known from only a handful of specimens and its morphology is a mosaic of key taxonomic features from two different clades. Here, new molecular evidence provides robust support for its correct association with a third different clade: Choriplax is placed in the superfamily Mopalioidea. Hemiarthrum is included in Cryptoplacoidea, as predicted from morphological evidence. Our multigene analysis of standard nuclear and mitochondrial markers demonstrates that the topology of the order Chitonida is divided into four clades, which have also been recovered in previous studies: Mopalioidea is sister to Cryptoplacoidea, forming a clade Acanthochitonina. The family Callochitonidae is sister to Acanthochitonina. Chitonoidea is resolved as the earliest diverging group within Chitonida. Consideration of this unexpected result for Choriplax and our well-supported phylogeny has revealed differing patterns of shell reduction separating the two superfamilies within Acanthochitonina. As in many molluscs, shell reduction as well as the de novo development of key shell features has occurred using different mechanisms, in multiple lineages of chitons.

Comparison of population-genetic structuring in congeneric kelp- versus rock-associated snails: a test of a dispersal-by-rafting hypothesis

Phylogeographic studies indicate that many marine invertebrates lacking autonomous dispersal ability are able to achieve trans-oceanic colonization by rafting on buoyant macroalgae. However, less is known about the impact of rafting on on-going population-genetic connectivity of intertidal species associated with buoyant macroalgae. We hypothesize that such species will have higher levels of population-genetic connectivity than those exploiting nonbuoyant substrates such as rock. We tested this hypothesis by comparing nuclear multilocus population-genetic structuring in two sister topshell species, which both have a planktonic larval phase but are fairly well segregated by their habitat preference of low-tidal bull-kelp holdfasts versus mid-to-low tidal bare rock. We analyzed population samples from four sympatric sites spanning 372 km of the east coast of southern New Zealand. The sampled region encompasses a 180 km wide habitat discontinuity and is influenced by a stable, northward coastal current. The level of connectivity was high in both species, and neither of them showed significant correlation between genetic and geographic distances. However, a significant negative partial correlation between genetic distance and habitat discontinuity was found in the rock-associated species, and estimates of migrant movement between sites were somewhat different between the two species, with the kelp-associated species more often yielding higher estimates across the habitat discontinuity, whereas the rock-associated species more often exhibited higher estimates between sites interspersed by rock habitats. We conclude that for species with substantial means of autonomous dispersal, the most conspicuous consequence of kelp dwelling may be enhanced long-distance dispersal across habitat discontinuities rather than a general increase of gene flow.

Delineating genetic groupings in continuously distributed species across largely homogeneous landscapes: a study of American black bears (Ursus americanus) in Ontario, Canada

There is a crucial need to understand the genetic consequences of landscape modifications on continuous populations that could become fragmented, and to evaluate the degree of differentiation of isolated populations that were historically part of the core. Using 15 microsatellite loci, we evaluated the genetic structure of American black bears ( Ursus americanus Pallas, 1780) across a vast, contiguous Ontario landscape (>1 × 106 km2) that largely represents their pre-European settlement distribution. Because geographic barriers are absent, we predicted that isolation by distance would drive genetic structure. We identified three genetic clusters (Northwest, Southeast, and Bruce Peninsula) that were less differentiated than when assessed with mtDNA, suggesting the influence of male-biased dispersal on large-scale genetic differentiation. Isolation by distance (r = 0.552, P = 0.001) was supported by a weak, clinal variation between Northwest and Southeast, illustrating the challenges to delineate populations in wide-ranging taxa. The Bruce Peninsula cluster, confined to a small area under strong anthropogenic pressures, was more differentiated from neighbouring clusters (FST > 0.13, P < 0.0001), with a genetic diversity corresponding to disjunct populations of black bears. Our results could be used in landscape genetics models to project the evolution of population differentiation based on upcoming landscape modifications in northern regions of North America.

Genetic population structure of the malaria vector Anopheles baimaii in north-east India using mitochondrial DNA

BACKGROUND

Anopheles baimaii is a primary vector of human malaria in the forest settings of Southeast Asia including the north-eastern region of India. Here, the genetic population structure and the basic population genetic parameters of An. baimaii in north-east India were estimated using DNA sequences of the mitochondrial cytochrome oxidase sub unit II (COII) gene.

METHODS

Anopheles baimaii were collected from 26 geo-referenced locations across the seven north-east Indian states and the COII gene was sequenced from 176 individuals across these sites. Fifty-seven COII sequences of An. baimaii from six locations in Bangladesh, Myanmar and Thailand from a previous study were added to this dataset. Altogether, 233 sequences were grouped into eight population groups, to facilitate analyses of genetic diversity, population structure and population history.

RESULTS

A star-shaped median joining haplotype network, unimodal mismatch distribution and significantly negative neutrality tests indicated population expansion in An. baimaii with the start of expansion estimated to be ~0.243 million years before present (MYBP) in north-east India. The populations of An. baimaii from north-east India had the highest haplotype and nucleotide diversity with all other populations having a subset of this diversity, likely as the result of range expansion from north-east India. The north-east Indian populations were genetically distinct from those in Bangladesh, Myanmar and Thailand, indicating that mountains, such as the Arakan mountain range between north-east India and Myanmar, are a significant barrier to gene flow. Within north-east India, there was no genetic differentiation among populations with the exception of the Central 2 population in the Barail hills area that was significantly differentiated from other populations.

CONCLUSIONS

The high genetic distinctiveness of the Central 2 population in the Barail hills area of the north-east India should be confirmed and its epidemiological significance further investigated. The lack of genetic population structure in the other north-east Indian populations likely reflects large population sizes of An. baimaii that, historically, were able to disperse through continuous forest habitats in the north-east India. Additional markers and analytical approaches are required to determine if recent deforestation is now preventing ongoing gene flow. Until such information is acquired, An. baimaii in north-east India should be treated as a single unit for the implementation of vector control measures.

Indoor fungal composition is geographically patterned and more diverse in temperate zones than in the tropics

Fungi are ubiquitous components of indoor human environments, where most contact between humans and microbes occurs. The majority of these organisms apparently play a neutral role, but some are detrimental to human lifestyles and health. Recent studies that used culture-independent sampling methods demonstrated a high diversity of indoor fungi distinct from that of outdoor environments. Others have shown temporal fluctuations of fungal assemblages in human environments and modest correlations with human activity, but global-scale patterns have not been examined, despite the manifest significance of biogeography in other microbial systems. Here we present a global survey of fungi from indoor environments (n = 72), using both taxonomic and phylogeny-informative molecular markers to determine whether global or local indoor factors determine indoor fungal composition. Contrary to common ecological patterns, we show that fungal diversity is significantly higher in temperate zones than in the tropics, with distance from the equator being the best predictor of phylogenetic community similarity. Fungal composition is significantly auto-correlated at the national and hemispheric spatial scales. Remarkably, building function has no significant effect on indoor fungal composition, despite stark contrasts between architecture and materials of some buildings in close proximity. Distribution of individual taxa is significantly range- and latitude-limited compared with a null model of randomized distribution. Our results suggest that factors driving fungal composition are primarily global rather than mediated by building design or function.