When One's Not Enough: Colony Pool-Seq Outperforms Individual-Based Methods for Assessing Introgression in Apis mellifera mellifera.

The human management of honey bees (Apis mellifera) has resulted in the widespread introduction of subspecies outside of their native ranges. One well known example of this is Apis mellifera mellifera, native to Northern Europe, which has now been significantly introgressed by the introduction of C lineage honey bees. Introgression has consequences for species in terms of future adaptive potential and long-term viability. However, estimating introgression in colony-living haplodiploid species is challenging. Previous studies have estimated introgression using individual workers, individual drones, multiple drones, and pooled workers. Here, we compare introgression estimates via three genetic approaches: SNP array, individual RAD-seq, and pooled colony RAD-seq. We also compare two statistical approaches: a maximum likelihood cluster program (ADMIXTURE) and an incomplete lineage sorting model (ABBA BABA). Overall, individual approaches resulted in lower introgression estimates than pooled colonies when using ADMIXTURE. However, the pooled colony ABBA BABA approach resulted in generally lower introgression estimates than all three ADMIXTURE estimates. These results highlight that sometimes one individual is not enough to assess colony-level introgression, and future studies that do use colony pools should not be solely dependent on clustering programs for introgression estimates.

Choice of molecular assay determines ranavirus detection probability and inferences about prevalence and occurrence.

Ranaviruses are emerging pathogens that can cause morbidity, mortality and population declines in ectothermic hosts; however, there is no standardized approach to diagnostics. Here, we compared the inter-assay variation and intra-assay precision among 2 commonly used quantitative PCRs (qPCRs), a conventional and a nested PCR assay (used as a gold standard), using laboratory-propagated ranavirus (FV3 and CMTV) and field-collected samples. A qPCR assay ('Leung') detected viral DNA in dilutions 2 orders of magnitude lower than other assays regardless of the viral lineage of the cultured isolate (FV3/CMTV). The second qPCR ('Brunner') was slightly more sensitive than the conventional PCR ('Mao' assay). For field samples, the Leung qPCR detected all known positives, while the Mao assay PCR only detected 2.5% of the positive samples. Amplicon sequences from the 2 conventional PCRs were shown to be useful for inferring viral lineage. Inaccurate results will bias estimates of the distribution and prevalence of ranaviruses, and together these findings emphasize that molecular assays should be chosen carefully in the context of study aims.

Climate structuring of Batrachochytrium dendrobatidis infection in the threatened amphibians of the northern Western Ghats, India.

Batrachochytrium dendrobatidis (Bd) is a pathogen killing amphibians worldwide. Its impact across much of Asia is poorly characterized. This study systematically surveyed amphibians for Bd across rocky plateaus in the northern section of the Western Ghats biodiversity hotspot, India, including the first surveys of the plateaus in the coastal region. These ecosystems offer an epidemiological model system since they are characterized by differing levels of connectivity, edaphic and climatic conditions, and anthropogenic stressors. One hundred and eighteen individuals of 21 species of Anura and Apoda on 13 plateaus ranging from 67 to 1179 m above sea level and 15.89 to 17.92° North latitude were sampled. Using qPCR protocols, 79% of species and 27% of individuals tested were positive for Bd. This is the first record of Bd in caecilians in India, the Critically Endangered Xanthophryne tigerina and Endangered Fejervarya cf. sahyadris. Mean site prevalence was 28.15%. Prevalence below the escarpment was 31.2% and 25.4% above. The intensity of infection (GE) showed the reverse pattern. Infection may be related to elevational temperature changes, thermal exclusion, inter-site connectivity and anthropogenic disturbance. Coastal plateaus may be thermal refuges from Bd. Infected amphibians represented a wide range of ecological traits posing interesting questions about transmission routes.

Micro-habitat distribution drives patch quality for sub-tropical rocky plateau amphibians in the northern Western Ghats, India.

The importance of patch quality for amphibians is frequently overlooked in distribution models. Here we demonstrate that it is highly important for the persistence of endemic and endangered amphibians found in the threatened and fragile ecosystems that are the rocky plateaus in Western Maharashtra, India. These plateaus are ferricretes of laterite and characterise the northern section of the Western Ghats/Sri Lanka Biodiversity Hotspot, the eighth most important global hotspot and one of the three most threatened by population growth. We present statistically supported habitat associations for endangered and data-deficient Indian amphibians, demonstrating significant relationships between individual species and their microhabitats. Data were collected during early monsoon across two seasons. Twenty-one amphibian taxa were identified from 14 lateritic plateaus between 67 and 1179m above sea level. Twelve of the study taxa had significant associations with microhabitats using a stepwise analysis of the AICc subroutine (distLM, Primer-e, v7). Generalist taxa were associated with increased numbers of microhabitat types. Non-significant associations are reported for the remaining 9 taxa. Microhabitat distribution was spatially structured and driven by climate and human activity. Woody plants were associated with 44% of high-elevation taxa. Of the 8 low-elevation taxa 63% related to water bodies and 60% of those were associated with pools. Rock size and abundance were important for 33% of high elevation specialists. Three of the 4 caecilians were associated with rocks in addition to soil and stream presence. We conclude the plateaus are individualistic patches whose habitat quality is defined by their microhabitats within climatic zones.

Ecological Variation in Response to Mass-Flowering Oilseed Rape and Surrounding Landscape Composition by Members of a Cryptic Bumblebee Complex.

The Bombus sensu stricto species complex is a widespread group of cryptic bumblebee species which are important pollinators of many crops and wild plants. These cryptic species have, until now, largely been grouped together in ecological studies, and so little is known about their individual colony densities, foraging ranges or habitat requirements, which can be influenced by land use at a landscape scale. We used mass-flowering oilseed rape fields as locations to sample bees of this complex, as well as the second most common visitor to oilseed rape B. lapidarius, and molecular RFLP methods to distinguish between the cryptic species. We then used microsatellite genotyping to identify sisters and estimate colony densities, and related both proportions of cryptic species and their colony densities to the composition of the landscape surrounding the fields. We found B. lucorum was the most common member of the complex present in oilseed rape followed by B. terrestris. B. cryptarum was also present in all but one site, with higher proportions found in the east of the study area. High numbers of bumblebee colonies were estimated to be using oilseed rape fields as a forage resource, with B. terrestris colony numbers higher than previous estimates from non-mass-flowering fields. We also found that the cryptic species responded differently to surrounding landscape composition: both relative proportions of B. cryptarum in samples and colony densities of B. lucorum were negatively associated with the amount of arable land in the landscape, while proportions and colony densities of other species did not respond to landscape variables at the scale measured. This suggests that the cryptic species have different ecological requirements (which may be scale-dependent) and that oilseed rape can be an important forage resource for many colonies of bumblebees. Given this, we recommend sustainable management of this crop to benefit bumblebees.

Lack of variation at phosphoglucose isomerase (pgi) in bumblebees: implications for conservation genetics studies.

Assessing genetic variation underlying ecologically important traits is increasingly of interest and importance in population and conservation genetics. For some groups generally useful markers exist for examining the relative role of selection and drift in shaping genetic diversity e.g. the major histocompatibility complex in vertebrates and self-incompatibility loci in plants. For invertebrates there is no such generally useful locus. However, phosphoglucose isomerase (Pgi) has been proposed as a useful functional marker in the conservation genetics of invertebrates. Where thermal microclimate varies, balanced polymorphisms may be maintained due to trade-offs between thermally stable and kinetically advantageous allelic forms. We here report very low levels of Pgi variation in bumblebees rendering this locus to be of little use as an adaptive marker in a conservation genetics context in this group. Potential explanations for this lack of variation are considered.

Effects of formalin preservation on stable carbon and nitrogen isotope signatures in Calanoid copepods: implications for the use of Continuous Plankton Recorder Survey samples in stable isotope analyses.

Preserved and archived organic material offers huge potential for the conduct of retrospective and long-term historical ecosystem reconstructions using stable isotope analyses, but because of isotopic exchange with preservatives the obtained values require validation. The Continuous Plankton Recorder (CPR) Survey is the most extensive long-term monitoring program for plankton communities worldwide and has utilised ships of opportunity to collect samples since 1931. To keep the samples intact for subsequent analysis, they are collected and preserved in formalin; however, previous studies have found that this may alter stable carbon and nitrogen isotope ratios in zooplankton. A maximum ~0.9‰ increase of δ(15) N and a time dependent maximum ~1.0‰ decrease of δ(13) C were observed when the copepod, Calanus helgolandicus, was experimentally exposed to two formalin preservatives for 12 months. Applying specific correction factors to δ(15) N and δ(13) C values for similarly preserved Calanoid species collected by the CPR Survey within 12 months of analysis may be appropriate to enable their use in stable isotope studies. The isotope values of samples stored frozen did not differ significantly from those of controls. Although the impact of formalin preservation was relatively small in this and other studies of marine zooplankton, changes in isotope signatures are not consistent across taxa, especially for δ(15) N, indicating that species-specific studies may be required.

Establishment and expansion of Lake Malawi rock fish populations after a dramatic Late Pleistocene lake level rise.

Major environmental events that fragment populations among multiple island habitats have potential to drive large-scale episodes of speciation and adaptive radiation. A recent palaeolimnological study of sediment cores indicated that Lake Malawi underwent major climate-driven desiccation events 75,000-135,000 years ago that lowered the water level to at least 580 m below the present state and severely reduced surface area. After this period, lake levels rose and stabilized, creating multiple discontinuous littoral rocky habitats. Here, we present evidence supporting the hypothesis that establishment and expansion of isolated philopatric rock cichlid populations occurred after this rise and stabilization of lake level. We studied the Pseudotropheus (Maylandia) species complex, a group with both allopatric and sympatric populations that differ in male nuptial colour traits and tend to mate assortatively. Using coalescent analyses based on mitochondrial DNA, we found evidence that populations throughout the lake started to expand and accumulate genetic diversity after the lake level rise. Moreover, most haplotypes were geographically restricted, and the greatest genetic similarities were typically among sympatric or neighbouring populations. This is indicative of limited dispersal and establishment of assortative mating among populations following the lake level rise. Together, this evidence is compatible with a single large-scale environmental event being central to evolution of spatial patterns of genetic and species diversity in P. (Maylandia) and perhaps other Lake Malawi rock cichlids. Equivalent climate-driven pulses of habitat formation and fragmentation may similarly have contributed to observed rapid and punctuated cladogenesis in other adaptive radiations.

Bumblebee flight distances in relation to the forage landscape.

1. Foraging range is a key aspect of the ecology of 'central place foragers'. Estimating how far bees fly under different circumstances is essential for predicting colony success, and for estimating bee-mediated gene flow between plant populations. It is likely to be strongly influenced by forage distribution, something that is hard to quantify in all but the simplest landscapes; and theories of foraging distance tend to assume a homogeneous forage distribution. 2. We quantified the distribution of bumblebee Bombus terrestris L. foragers away from experimentally positioned colonies, in an agricultural landscape, using two methods. We mass-marked foragers as they left the colony, and analysed pollen from foragers returning to the colonies. The data were set within the context of the 'forage landscape': a map of the spatial distribution of forage as determined from remote-sensed data. To our knowledge, this is the first time that empirical data on foraging distances and forage availability, at this resolution and scale, have been collected and combined for bumblebees. 3. The bees foraged at least 1.5 km from their colonies, and the proportion of foragers flying to one field declined, approximately linearly, with radial distance. In this landscape there was great variation in forage availability within 500 m of colonies but little variation beyond 1 km, regardless of colony location. 4. The scale of B. terrestris foraging was large enough to buffer against effects of forage patch and flowering crop heterogeneity, but bee species with shorter foraging ranges may experience highly variable colony success according to location.

Laboratory mating trials indicate incipient speciation by sexual selection among populations of the cichlid fish Pseudotropheus zebra from Lake Malawi.

It has been suggested that sexual selection may have played a major role in the rapid evolution of hundreds of species of cichlid fishes in Lake Malawi. We report the results of a laboratory test of assortative mating among Lake Malawi cichlid fishes from five closely related geographical populations differing in male courtship colour. Paternity of clutches was tested using microsatellite DNA typing of offspring. Out of 1955 offspring typed, 1296 (66.3%) were sired by the male from the same population as the female, which is more than three times the rate expected if females do not differentiate among males of the different populations (20%). This result indicates that mate preferences of geographical races are strongly differentiated, consistent with the races representing incipient geographical species diverging under sexual selection exerted by female preferences for different male courtship colours.

Divergent selection during speciation of Lake Malawi cichlid fishes inferred from parallel radiations in nuptial coloration.

Repeated evolution of the same phenotypic difference during independent episodes of speciation is strong evidence for selection during speciation. More than 1,000 species of cichlids, >10% of the world's freshwater fish species, have arisen within the past million years in Lakes Malawi and Victoria in eastern Africa. Many pairs of closely related sympatric species differ in their nuptial coloration in very similar ways. Nuptial coloration is important in their mate choice, and speciation by sexual selection on genetically or ecologically constrained variation in nuptial coloration had been proposed, which would repeatedly produce similar nuptial types in different populations, a prediction that was difficult to test in the absence of population-level phylogenies. We measured genetic similarity between individuals within and between populations, species, and lake regions by typing 59 individuals at >2,000 polymorphic genetic loci. From these data, we reconstructed, to our knowledge, the first larger species level phylogeny for the most diverse group of Lake Malawi cichlids. We used the genetic and phylogenetic data to test the divergent selection scenario against colonization, character displacement, and hybridization scenarios that could also explain diverse communities. Diversity has arisen by replicated radiations into the same color types, resulting in phenotypically very different, yet closely related, species within and phenotypically highly similar yet unrelated sets of species between regions, which is consistent with divergent selection during speciation and is inconsistent with colonization and character displacement models.

Nuclear markers reveal unexpected genetic variation and a Congolese-Nilotic origin of the Lake Victoria cichlid species flock.

Phylogenetic analyses based on mitochondrial (mt) DNA have indicated that the cichlid species flock of the Lake Victoria region is derived from a single ancestral species found in East African rivers, closely related to the ancestor of the Lake Malawi cichlid species flock. The Lake Victoria flock contains ten times less mtDNA variation than the Lake Malawi radiation, consistent with current estimates of the ages of the lakes. We present results of a phylogenetic investigation using nuclear (amplified fragment length polymorphism) markers and a wider coverage of riverine haplochromines. We demonstrate that the Lake Victoria-Edward flock is derived from the morphologically and ecologically diverse cichlid genus Thoracochromis from the Congo and Nile, rather than from the phenotypically conservative East African Astatotilapia. This implies that the ability to express much of the morphological diversity found in the species flock may by far pre-date the origin of the flock. Our data indicate that the nuclear diversity of the Lake Victoria-Edward species flock is similar to that of the Lake Malawi flock, indicating that the genetic diversity is considerably older than the 15 000 years that have passed since the lake began to refill. Most of this variation is manifested in trans-species polymorphisms, indicating very recent cladogenesis from a genetically very diverse founder stock. Our data do not confirm strict monophyly of either of the species flocks, but raise the possibility that these flocks have arisen from hybrid swarms.