Mating behaviour in the endangered Sonoran topminnow: speciation in action

Two species of the endangered Sonoran topminnow, Poeciliopsis occidentalis and P. sonoriensis, occur in two isolated drainage systems in southeast Arizona, U.S.A., and are allopatric throughout their range. Although these two taxa are morphologically very similar, and have been previously described as the same species or subspecies, several molecular studies have since indicated that they differ more than their morphology suggests. To determine whether the behaviours of the two species function as premating barriers to reproduction, we investigated their mating preferences and behavioural patterns in a laboratory setting. Results from no-choice mating observations showed that the mating behaviours of the two species differ. Observations conducted during multiple-choice mating trials provided evidence of assortative mating, suggesting an early stage of premating reproductive isolation.     

Initial stages of reproductive isolation in two species of the endangered Sonoran topminnow

Long-term geographic isolation can result in reproductive incompatibilities due to forces such as mutation, genetic drift, and differential selection. In the Sonoran topminnow, molecular genetic studies of mtDNA, microsatellites, and MHC genes have shown that the endangered Gila and Yaqui topminnows are substantially different, suggesting that divergence took place approximately two million years ago. Here we examined hybrid crosses and backcrosses between these two allopatric taxa to evaluate the accumulation of postmating barriers to reproduction. These results are then compared with results from a previous study where male topminnows were shown to mate assortatively with conspecific females. Despite their preference for conspecific mates, both types of interspecific crosses successfully produced offspring. There was evidence of reduced hybrid fitness, including smaller mean brood size and male-biased sex ratio, for some classes of backcrosses. Brood sizes and interbrood intervals varied significantly when hybrids were subdivided into different cross categories. Our results illustrate the importance of distinctly defining hybrid classes in studies of reproductive isolation. To our knowledge, this is the first such detailed evolutionary analysis in endangered fish taxa.     

Using microsatellite and MHC variation to identify species, ESUs, and MUs in the endangered Sonoran topminnow

Highly variable loci can provide insight into the recognition of species, evolutionarily significant units (ESUs) and management units (MUs). In general, the ESU and MU categories are thought to be reflective of adaptive differences between them. Here we examine this premise by presenting a comprehensive examination of genetic variation for both microsatellite loci and a major histocompatibility complex (MHC) locus, thought to be of adaptive significance, in the endangered Sonoran topminnow. The extent of variation for the microsatellite loci and the MHC gene within the 13 populations of the Gila topminnow is highly correlated, suggesting that nonselective factors have played an important role in influencing variation within and between populations for the MHC locus. Therefore, using all of these loci, we found that the eight natural populations of the Gila topminnow fell into two different ESUs, one of which had four different MUs. The source of the Boyce Thompson sample, a population that was used extensively for restocking, appeared to be Monkey Spring. The source of the Watson Wash population also appeared to be Monkey Spring (or Boyce Thompson). The newly colonized Santa Cruz River population, which had the most genetic variation of any Gila topminnow population, appeared to descend primarily from Sonoita Creek populations. The Yaqui topminnow, presently considered another subspecies of the Sonoran topminnow, was very distinct for both microsatellite (only two of 25 alleles found in the Yaqui were in any of the Gila topminnow samples) and MHC alleles (nonoverlapping sets of alleles for the two groups). As a result, it appeared that the taxonomic status of the two subspecies should be re-evaluated and that full species status for Gila and Yaqui topminnows was appropriate. There was evidence for the importance of long-term selection at the MHC locus in the higher rate of nonsynonymous than synonymous substitution. In addition, there appeared to have been a duplication of the MHC locus that was present in most of the fish in six of the natural populations of the Gila topminnow.     

Sexual conflict does not drive reproductive isolation in experimental populations of Drosophila pseudoobscura

Sexual conflict has been predicted to drive reproductive isolation by generating arbitrary but rapid coevolutionary changes in reproductive traits among allopatric populations. A testable prediction of this proposal is that allopatric populations experiencing different levels of sexual conflict should exhibit different levels of reproductive isolation. We tested this prediction using experimentally evolved populations of the promiscuous Drosophila pseudoobscura. We manipulated sexual conflict by enforcing either monogamy, maintaining natural levels of promiscuity, or elevating promiscuity. Within each treatment, we carried out sympatric and allopatric crosses using replicated populations and examined pre-zygotic (number of mating pairs, mating speed and copulation duration) and post-zygotic (hybrid inviability and sterility) indicators of reproductive isolation. After 50 generations of selection, none of the measures conformed to predictions of sexual conflict driving reproductive isolation. Our results cannot be explained by lack of genetic variation or weak selection and suggest that sexual conflict may not be a widespread engine of speciation.     

Patch quality and context, but not patch number, drive multi-scale colonization dynamics in experimental aquatic landscapes

Colonization and extinction are primary drivers of local population dynamics, community structure, and spatial patterns of biological diversity. Existing paradigms of island biogeography, metapopulation biology, and metacommunity ecology, as well as habitat management and conservation biology based on those paradigms, emphasize patch size, number, and isolation as primary characteristics influencing colonization and extinction. Habitat selection theory suggests that patch quality could rival size, number, and isolation in determining rates of colonization and resulting community structure. We used naturally colonized experimental landscapes to address four issues: (a) how do colonizing aquatic beetles respond to variation in patch number, (b) how do they respond to variation in patch quality, (c) does patch context affect colonization dynamics, and (d) at what spatial scales do beetles respond to habitat variation? Increasing patch number had no effect on per patch colonization rates, while patch quality and context were critical in determining colonization rates and resulting patterns of abundance and species richness at multiple spatial scales. We graphically illustrate how variation in immigration rates driven by perceived predation risk (habitat quality) can further modify dynamics of the equilibrium theory of island biogeography beyond predator-driven effects on extinction rates. Our data support the importance of patch quality and context as primary determinants of colonization rate, occupancy, abundance, and resulting patterns of species richness, and reinforce the idea that management of metapopulations for species preservation, and metacommunities for local and regional diversity, should incorporate habitat quality into the predictive equation.     

Growth and demography of a declining,endangered cactus in the Sonoran Desert

Increases in the incidence and severity of drought threaten the viability of rare plants in arid regions. The endangered Nichol's Turk's head cactus (Echinocactus horizonthalonius Lemaire var. nicholii L. Benson) occurs only in four small, isolated populations in the Sonoran Desert of North America. Since 1995 we have monitored a population in southeastern Arizona (USA). Here we report 23 years of observations on abundance, growth, mortality, flowering and recruitment. Abundance of plants decreased from 132 in 1996 to 40 in 2017, with 100 individuals recruited and 203 dying during the study. Individual plants grew slowly, increasing annually by an average of 0.22 cm (95% confidence interval, 0.18–0.26 cm) in diameter and 0.27 cm (0.20–0.33 cm) in height. Growth was slowest when drought was most severe and slowed as plants reached maximum size. Annual mortality increased markedly across the study period and did not vary with plant size. Annual probability of flowering increased as plants increased in diameter but not in height, and varied with precipitation and drought but not with mean annual temperature. Recruitment was higher when average temperature was higher and the number of recruits per capita increased across the study period. The annual rate of change in abundance averaged −6%, but shifted markedly from −1% during 1995–2008 to −11% during 2008–2017. Our results indicate that the population's decline was not a consequence of failed recruitment but of increased mortality, which we discuss in the context of climate and herbivory.     

Fitness, genetic load and purging in experimental populations of the housefly

The relative effects of purging of the genetic load versus thefixation of deleterious alleles, under inbreeding, will influencea population's probability of extinction. The relative contributionof these two phenomena is expected to depend upon the rate ofinbreeding. A further complication is due to the fact that a purgingof the genetic load in one environment does not necessarily implya purging of the genetic load in other environments. To addressthese two issues, we compare fitness and genetic load in populationsexperiencing similar levels of inbreeding, but occurring as either ashort-term bottleneck or as a consequence of long-term reducedpopulation size, over a range of environments. Inbred populationshave consistently lower fitness than outbred populations acrossall environments tested. However, the bottlenecked populationssuffer less inbreeding depression for a given level of inbreeding,whether or not challenged by novel environments, than populationskept at a constant small size. The results of this study demonstratethat populations initiated from a small number of founders are ableto recover fitness and survive novel environmental challenges,provided that habitat is available for rapid population growth.     

An experimental test of character displacement's role in promoting postmating isolation between conspecific populations in contrasting competitive environments

Ecological character displacement takes place when two closely related species co-occur in only part of their geographical range, and selection to minimize competition between them promotes divergence in resource-use traits in sympatry but not in allopatry. Because populations sympatric with the heterospecific competitor will experience a different competitive environment than conspecific populations in allopatry, conspecific populations from these two competitive environments will also diverge in resource traits as an indirect consequence of interspecific ecological character displacement. Ultimately, ecologically dependent postmating isolation may arise between conspecific populations from these divergent competitive environments if offspring produced by matings between them are competitively inferior in either type of competitive environment. Yet, there are no direct tests of character displacement's role in initiating such postmating isolation. Here, we present a test by comparing the phenotypes and performances of spadefoot toad tadpoles produced from between-competitive-environment (BCE) matings versus those produced from within-competitive-environment (WCE) matings. When raised with naturally occurring competitors, BCE offspring grew significantly less and were significantly smaller than WCE offspring. BCE offspring generally performed worse even when raised alone, suggesting that they may have harbored intrinsic genetic incompatibilities. Moreover, the difference in growth and body size of BCE versus WCE offspring was significantly greater when each was raised with competitors than when each was raised alone, suggesting that BCE tadpoles were competitively inferior to WCE tadpoles. Presumably, this enhanced difference arose because BCE tadpoles produced an intermediate resource-use phenotype that is less well adapted to either competitive environment. Because larval size is under strong, positive, directional selection, reduced growth and size of BCE offspring may diminish gene flow between populations in divergent competitive environments, thereby generating postmating isolation. Thus, postmating isolation between conspecific populations, and possibly even speciation, may arise as a by-product of interactions between species.     

Selective recovery of founder genetic diversity in aquacultural broodstocks and captive, endangered fish populations

Hatchery broodstocks used for genetic conservation or aquaculture may represent their ancestral gene pools rather poorly. This is especially likely when the fish that found a broodstock are close relatives of each other. We re-analysed microsatellite data from a breeding experiment on red sea bream to demonstrate how lost genetic variation might be recovered when gene frequencies have been distorted by consanguineous founders in a hatchery. A minimal-kinship criterion based on a relatedness estimator was used to select subsets of breeders which represented the maximum number of founder lineages (i.e., carried the fewest identical copies of ancestral genes). UPGMA clustering of Nei's genetic distances grouped these selected subsets with the parental gene pool, rather than with the entire, highly drifted offspring generation. The selected subsets also captured much of the expected heterozygosity and allelic diversity of the parental gene pool. Independent pedigree data on the same fish showed that the selected subsets had more contributing parents and more founder equivalents than random subsets of the same size. The estimated mean coancestry was lower in the selected subsets, meaning that inbreeding in subsequent generations would be lower if they were used as breeders. The procedure appears suitable for reducing the genetic distortion due to consanguineous and over-represented founders of a hatchery gene pool.     

The influence of natural barriers in shaping the genetic structure of Maharashtra populations

Background

The geographical position of Maharashtra state makes it rather essential to study the dispersal of modern humans in South Asia. Several hypotheses have been proposed to explain the cultural, linguistic and geographical affinity of the populations living in Maharashtra state with other South Asian populations. The genetic origin of populations living in this state is poorly understood and hitherto been described at low molecular resolution level.

Methodology/Principal Findings

To address this issue, we have analyzed the mitochondrial DNA (mtDNA) of 185 individuals and NRY (non-recombining region of Y chromosome) of 98 individuals belonging to two major tribal populations of Maharashtra, and compared their molecular variations with that of 54 South Asian contemporary populations of adjacent states. Inter and intra population comparisons reveal that the maternal gene pool of Maharashtra state populations is composed of mainly South Asian haplogroups with traces of east and west Eurasian haplogroups, while the paternal haplogroups comprise the South Asian as well as signature of near eastern specific haplogroup J2a.

Conclusions/Significance

Our analysis suggests that Indian populations, including Maharashtra state, are largely derived from Paleolithic ancient settlers; however, a more recent (∼10 Ky older) detectable paternal gene flow from west Asia is well reflected in the present study. These findings reveal movement of populations to Maharashtra through the western coast rather than mainland where Western Ghats-Vindhya Mountains and Narmada-Tapti rivers might have acted as a natural barrier. Comparing the Maharastrian populations with other South Asian populations reveals that they have a closer affinity with the South Indian than with the Central Indian populations.     

Patterns of genetic variation in isolated Danish populations of the endangered butterfly Euphydryas aurinia

In the present study, we report an investigation on molecular variation in the endangered univoltine butterfly Euphydryas aurinia (Rottemburg, 1775), a species heavily affected by habitat degradation and fragmentation in Denmark. Levels of genetic variation in extant populations were estimated using six variable number tandem repeat loci and were found to be low compared to other butterfly species with low migration rates. An analysis of genetic structure, based on both allele frequencies and genotype distributions, divided the entire sample into four distinct clusters. This was partially concordant with the a priori subdivision based on collection areas. An overall F_ST value of 0.16 (pairwise values ranging from 0.087–0.276) indicated restrictions of gene flow. Especially two populations had higher F_ST values than the others, suggesting their isolation, and showed signs of bottlenecks/founder events. One population deviated significantly from Hardy–Weinberg equilibrium, suggesting a possible Wahlund effect or the presence of null alleles. The results suggest habitat fragmentation, resulting in genetic drift and possibly inbreeding. Future management is therefore recommended to increase gene flow between the remaining populations while habitats are restored in order to increase carrying capacity. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 677–687.     

Mapping the genetic architecture of complex traits in experimental populations

SUMMARY: Understanding how interactions among set of genes affect diverse phenotypes is having a greater impact on biomedical research, agriculture and evolutionary biology. Mapping and characterizing the isolated effects of single quantitative trait locus (QTL) is a first step, but we also need to assemble networks of QTLs and define non-additive interactions (epistasis) together with a host of potential environmental modulators. In this article, we present a full-QTL model with which to explore the genetic architecture of complex trait in multiple environments. Our model includes the effects of multiple QTLs, epistasis, QTL-by-environment interactions and epistasis-by-environment interactions. A new mapping strategy, including marker interval selection, detection of marker interval interactions and genome scans, is used to evaluate putative locations of multiple QTLs and their interactions. All the mapping procedures are performed in the framework of mixed linear model that are flexible to model environmental factors regardless of fix or random effects being assumed. An F-statistic based on Henderson method III is used for hypothesis tests. This method is less computationally greedy than corresponding likelihood ratio test. In each of the mapping procedures, permutation testing is exploited to control for genome-wide false positive rate, and model selection is used to reduce ghost peaks in F-statistic profile. Parameters of the full-QTL model are estimated using a Bayesian method via Gibbs sampling. Monte Carlo simulations help define the reliability and efficiency of the method. Two real-world phenotypes (BXD mouse olfactory bulb weight data and rice yield data) are used as exemplars to demonstrate our methods. AVAILABILITY: A software package is freely available at http://ibi.zju.edu.cn/software/qtlnetwork     

Prediction of offspring fitness based on parental genetic diversity in endangered salmonid populations

By simulating all possible offspring genotypic combinations based on pair-wise matings in three vulnerable or endangered populations of salmonids (two Atlantic salmon Salmo salar and one Arctic charr Salvelinus alpinus ), family level estimates of potential offspring genetic diversity were calculated using several commonly used genetic diversity estimators (     ,     and H _EST). By using this approach, the genetic variability of offspring could be predicted accurately based solely on parental genotype data. In addition, significant associations between offspring genetic variability and fitness-related traits were observed in the two Atlantic salmon populations and these associations were consistent between fitness measures and either estimated or observed genetic diversity. This indicates that associations with offspring fitness measures can be predicted without genetic analysis of the offspring themselves.     

Phylogeography of the endangered Otago skink, Oligosoma otagense: population structure, hybridisation and genetic diversity in captive populations

Climatic cooling and substantial tectonic activity since the late Miocene have had a pronounced influence on the evolutionary history of the fauna of New Zealand's South Island. However, many species have recently experienced dramatic range reductions due to habitat fragmentation and the introduction of mammalian predators and competitors. These anthropogenic impacts have been particularly severe in the tussock grasslands of the Otago region. The Otago skink (Oligosoma otagense), endemic to the region, is one of the most critically endangered vertebrates in New Zealand. We use mitochondrial DNA sequence data to investigate the evolutionary history of the Otago skink, examine its population genetic structure, and assess the level of genetic diversity in the individuals in the captive breeding program. Our data indicate that the Otago skink diverged from its closest relatives in the Miocene, consistent with the commencement of tectonic uplift of the Southern Alps. However, there is evidence for past introgression with the scree skink (O. waimatense) in the northern Otago-southern Canterbury region. The remnant populations in eastern Otago and western Otago are estimated to have diverged in the mid-Pliocene, with no haplotypes shared between these two regions. This divergence accounts for 95% of the genetic diversity in the species. Within both regions there is strong genetic structure among populations, although shared haplotypes are generally evident between adjacent localities. Although substantial genetic diversity is present in the captive population, all individuals originate from the eastern region and the majority had haplotypes that were not evident in the intensively managed populations at Macraes Flat. Our data indicate that eastern and western populations should continue to be regarded as separate management units. Knowledge of the genetic diversity of the breeding stock will act to inform the captive management of the Otago skink and contribute to a key recovery action for the species.     

Inbreeding rate modifies the dynamics of genetic load in small populations

The negative fitness consequences of close inbreeding are widely recognized, but predicting the long-term effects of inbreeding and genetic drift due to limited population size is not straightforward. As the frequency and homozygosity of recessive deleterious alleles increase, selection can remove (purge) them from a population, reducing the genetic load. At the same time, small population size relaxes selection against mildly harmful mutations, which may lead to accumulation of genetic load. The efficiency of purging and the accumulation of mutations both depend on the rate of inbreeding (i.e., population size) and on the nature of mutations. We studied how increasing levels of inbreeding affect offspring production and extinction in experimental Drosophila littoralis populations replicated in two sizes, N = 10 and N = 40. Offspring production and extinction were measured over 25 generations concurrently with a large control population. In the N = 10 populations, offspring production decreased strongly at low levels of inbreeding, then recovered only to show a consistent subsequent decline, suggesting early expression and purging of recessive highly deleterious alleles and subsequent accumulation of mildly harmful mutations. In the N = 40 populations, offspring production declined only after inbreeding reached higher levels, suggesting that inbreeding and genetic drift pose a smaller threat to population fitness when inbreeding is slow. Our results suggest that highly deleterious alleles can be purged in small populations already at low levels of inbreeding, but that purging does not protect the small populations from eventual genetic deterioration and extinction.     

The architecture of river networks can drive the evolutionary dynamics of aquatic populations

It is widely recognized that physical landscapes can shape genetic variation within and between populations. However, it is not well understood how riverscapes, with their complex architectures, affect patterns of neutral genetic diversity. Using a spatially explicit agent‐based modeling (ABM) approach, we evaluate the genetic consequences of dendritic river shapes on local population structure. We disentangle the relative contribution of specific river properties to observed patterns of genetic variation by evaluating how different branching architectures and downstream flow regimes affect the genetic structure of populations situated within river networks. Irrespective of the river length, our results illustrate that the extent of river branching, confluence position, and levels of asymmetric downstream migration dictate patterns of genetic variation in riverine populations. Comparisons between simple and highly branched rivers show a 20‐fold increase in the overall genetic diversity and a sevenfold increase in the genetic differentiation between local populations. Given that most rivers have complex architectures, these results highlight the importance of incorporating riverscape information into evolutionary models of aquatic species and could help explain why riverine fishes represent a disproportionately large amount of global vertebrate diversity per unit of habitable area.     

High levels of variation despite genetic fragmentation in populations of the endangered mountain pygmy-possum, Burramys parvus, in alpine Australia

In endangered mammals, levels of genetic variation are often low and this is accompanied by genetic divergence among populations. The mountain pygmy-possum (Burramys parvus) is an endangered marsupial restricted to the alpine region of Victoria and New South Wales, Australia. By scoring variation at eight microsatellite loci, we found that B. parvus populations exhibit high levels of genetic divergence and fall into three distinct groups from the northern, central and southern areas of the distribution of this species, consistent with previous assessments of mitochondrial DNA variation. F(ST) values between populations from these regions ranged from 0.19 to 0.54. Within the central area, there was further genetic fragmentation, and a linear association between genetic and geographical distance. This pattern is likely to reflect limited dispersal across barriers despite the fact that individual B. parvus can move several kilometres. Levels of genetic variation within populations were high with the exception of a southern population where there was evidence of inbreeding. From a conservation perspective, all three areas where B. parvus are found should be considered as separate gene pools; management of populations within these areas needs to take into account the low gene flow between populations, as well as threats posed by roads, resorts and other developments in the alpine region. The low genetic variability and inbreeding in the southern population is of particular concern given the high levels of variability in other B. parvus populations.     

Across the great divide: genetic forensics reveals misidentification of endangered cutthroat trout populations

Accurate assessment of species identity is fundamental for conservation biology. Using molecular markers from the mitochondrial and nuclear genomes, we discovered that many putatively native populations of greenback cutthroat trout (Oncorhynchus clarkii stomias) comprised another subspecies of cutthroat trout, Colorado River cutthroat trout (Oncorhynchus clarkii pleuriticus). The error can be explained by the introduction of Colorado River cutthroat trout throughout the native range of greenback cutthroat trout in the late 19th and early 20th centuries by fish stocking activities. Our results suggest greenback cutthroat trout within its native range is at a higher risk of extinction than ever before despite conservation activities spanning more than two decades.