Sensitivity analysis of effective population size to demographic parameters in house sparrow populations

The ratio between the effective and the census population size, , is an important measure of the long‐term viability and sustainability of a population. Understanding which demographic processes that affect most will improve our understanding of how genetic drift and the probability of fixation of alleles is affected by demography. This knowledge may also be of vital importance in management of endangered populations and species. Here, we use data from 13 natural populations of house sparrow (Passer domesticus) in Norway to calculate the demographic parameters that determine . Using the global variance‐based Sobol’ method for the sensitivity analyses, we found that was most sensitive to demographic variance, especially among older individuals. Furthermore, the individual reproductive values (that determine the demographic variance) were most sensitive to variation in fecundity. Our results draw attention to the applicability of sensitivity analyses in population management and conservation. For population management aiming to reduce the loss of genetic variation, a sensitivity analysis may indicate the demographic parameters towards which resources should be focused. The result of such an analysis may depend on the life history and mating system of the population or species under consideration, because the vital rates and sex–age classes that is most sensitive to may change accordingly.     

Effective number of breeders provides a link between interannual variation in stream flow and individual reproductive contribution in a stream salmonid

The effective number of breeders that give rise to a cohort (N_b) is a promising metric for genetic monitoring of species with overlapping generations; however, more work is needed to understand factors that contribute to variation in this measure in natural populations. We tested hypotheses related to interannual variation in N_b in two long‐term studies of brook trout populations. We found no supporting evidence for our initial hypothesis that reflects (defined as the number of adults in a population at the time of reproduction). was stable relative to and did not follow trends in abundance (one stream negative, the other positive). We used stream flow estimates to test the alternative hypothesis that environmental factors constrain N_b. We observed an intermediate optimum autumn stream flow for both (R² = 0.73, P = 0.02) and full‐sibling family evenness (R² = 0.77, P = 0.01) in one population and a negative correlation between autumn stream flow and full‐sib family evenness in the other population (r = ?0.95, P = 0.02). Evidence for greater reproductive skew at the lowest and highest autumn flow was consistent with suboptimal conditions at flow extremes. A series of additional tests provided no supporting evidence for a related hypothesis that density‐dependent reproductive success was responsible for the lack of relationship between N_b and N_C (so‐called genetic compensation). This work provides evidence that N_b is a useful metric of population‐specific individual reproductive contribution for genetic monitoring across populations and the link we provide between stream flow and N_b could be used to help predict population resilience to environmental change.     

Purging putative siblings from population genetic data sets: a cautionary view

Interest has surged recently in removing siblings from population genetic data sets before conducting downstream analyses. However, even if the pedigree is inferred correctly, this has the potential to do more harm than good. We used computer simulations and empirical samples of coho salmon to evaluate strategies for adjusting samples to account for family structure. We compared performance in full samples and sibling‐reduced samples of estimators of allele frequency (), population differentiation () and effective population size (). Results: (i) unless simulated samples included large family groups together with a component of unrelated individuals, removing siblings generally reduced precision of and ; (ii) based on the linkage disequilibrium method was largely unbiased using full random samples but became increasingly upwardly biased under aggressive purging of siblings. Under nonrandom sampling (some families over‐represented), using full samples was downwardly biased; removing just the right ‘Goldilocks’ fraction of siblings could produce an unbiased estimate, but this sweet spot varied widely among scenarios; (iii) weighting individuals based on the inferred pedigree (to produce a best linear unbiased estimator, BLUE) maximized precision of when the inferred pedigree was correct but performed poorly when the pedigree was wrong; (iv) a variant of sibling removal that leaves intact small sibling groups appears to be more robust to errors in inferences about family structure. Our results illustrate the complex challenges posed by presence of family structure, suggest that no single optimal solution exists and argue for caution in adjusting population genetic data sets for the presence of putative siblings without fully understanding the consequences.     

Species turnover (β‐diversity) in ectomycorrhizal fungi linked to uptake capacity

Ectomycorrhizal (EcM) fungal communities may be shaped by both deterministic and stochastic processes, potentially influencing ecosystem development and function. We evaluated community assembly processes for EcM fungi of Pseudotsuga menziesii among 12 sites up to 400 km apart in southwest British Columbia (Canada) by investigating species turnover (β‐diversity) in relation to soil nitrogen (N) availability and physical distance. We then examined functional traits for an N‐related niche by quantifying net fluxes of , and protons on excised root tips from three contrasting sites using a microelectrode ion flux measurement system. EcM fungal communities were well aligned with soil N availability and pH, with no effect of site proximity (distance–decay curve) on species assemblages. Species turnover was significant (β_1/2 = 1.48) along soil N gradients, with many more Tomentella species on high N than low N soils, in contrast to Cortinarius species. Ammonium uptake was greatest in the spring on the medium and rich sites and averaged over 190 nmol/m²/s for Tomentella species. The lowest uptake rates of were by nonmycorrhizal roots of axenically grown seedlings (10 nmol/m²/s), followed by Cortinarius species (60 nmol/m²/s). EcM roots from all sites displayed only marginal uptake of nitrate (8.3 nmol/m²/s). These results suggest uptake capacity is an important functional trait influencing the assembly of EcM fungal communities. The diversity of EcM fungal species across the region arguably provides critical belowground adaptations to organic and inorganic N supply that are integral to temperate rainforest ecology.     

Crossing to safety: dispersal,colonization and mate choice in evolutionarily distinct populations of Steller sea lions,Eumetopias jubatus

Population growth typically involves range expansion and establishment of new breeding sites, while the opposite occurs during declines. Although density dependence is widely invoked in theoretical studies of emigration and colonization in expanding populations, few empirical studies have documented the mechanisms. Still fewer have documented the direction and mechanisms of individual transfer in declining populations. Here, we screen large numbers of pups sampled on their natal rookeries for variation in mtDNA (n = 1106) and 16 microsatellite loci (n = 588) and show that new Steller sea lion breeding sites did not follow the typical paradigm and were instead colonized by sea lions from both a declining (Endangered) population and an increasing population. Dispersing individuals colonized rookeries in the distributional hiatus between two evolutionarily distinct ( = 0.222,  = 0.053, K = 2) metapopulations recently described as separate subspecies. Hardy–Weinberg, mixed‐stock and relatedness analysis revealed levels of interbreeding on the new rookeries that exclude (i) assortative mating among eastern and western forms, and (ii) inbreeding avoidance as primary motivations for dispersal. Positive and negative density dependence is implicated in both cases of individual transfer. Migration distance limits, and conspecific attraction and performance likely influenced the sequence of rookery colonizations. This study demonstrates that resource limitation may trigger an exodus of breeding animals from declining populations, with substantial impacts on distribution and patterns of genetic variation. It also revealed that this event is rare because colonists dispersed across an evolutionary boundary, suggesting that the causative factors behind recent declines are unusual or of larger magnitude than normally occur.     

The association between polyploidy and clonal reproduction in diploid and tetraploid Chamerion angustifolium

Clonal reproduction is associated with the incidence of polyploidy in flowering plants. This pattern may arise through selection for increased clonality in polyploids compared to diploids to reduce mixed‐ploidy mating. Here, we test whether clonal reproduction is greater in tetraploid than diploid populations of the mixed‐ploidy plant, Chamerion angustifolium, through an analysis of the size and spatial distribution of clones in natural populations using AFLP genotyping and a comparison of root bud production in a greenhouse study. Natural tetraploid populations (N = 5) had significantly more AFLP genotypes ( = 10.8) than diploid populations ( = 6.0). Tetraploid populations tended to have fewer ramets per genotype and fewer genotypes with >1 ramet. In a spatial autocorrelation analysis, ramets within genotypes were more spatially aggregated in diploid populations than in tetraploid populations. In the greenhouse, tetraploids allocated 90.4% more dry mass to root buds than diploids, but tetraploids produced no more root buds and 44% fewer root buds per unit root mass than diploids. Our results indicate that clonal reproduction is significant in most populations, but tetraploid populations are not more clonal than diploids, nor are their clones more spatially aggregated. As a result, tetraploids may be less sheltered from mixed‐ploidy mating and diploids more exposed to inbreeding, the balance of which could influence the establishment of tetraploids in diploid populations.     

When genes move farther than offspring: gene flow by male gamete dispersal in the highly philopatric bat species Thyroptera tricolor

For species characterized by philopatry of both sexes, mate selection represents an important behaviour for inbreeding avoidance, yet the implications for gene flow are rarely quantified. Here, we present evidence of male gamete‐mediated gene flow resulting from extra‐group mating in Spix's disc‐winged bat, Thyroptera tricolor, a species which demonstrates all‐offspring philopatry. We used microsatellite and capture–recapture data to characterize social group structure and the distribution of mated pairs at two sites in southwestern Costa Rica over four breeding seasons. Relatedness and genetic spatial autocorrelation analyses indicated strong kinship within groups and over short distances (<50 m), resulting from matrilineal group structure and small roosting home ranges (~0.2 ha). Despite high relatedness among‐group members, observed inbreeding coefficients were low (F_IS = 0.010 and 0.037). Parentage analysis indicated mothers and offspring belonged to the same social group, while fathers belonged to different groups, separated by large distances (~500 m) when compared to roosting home ranges. Simulated random mating indicated mate choice was not based on intermediate levels of relatedness, and mated pairs were less related than adults within social groups on average. Isolation‐by‐distance (IBD) models of genetic neighbourhood area based on father–offspring distances provided direct estimates of mean gamete dispersal distances () > 10 roosting home range equivalents. Indirect estimates based on genetic distance provided even larger estimates of , indicating direct estimates were biased low. These results suggest extra‐group mating reduces the incidence of inbreeding in T. tricolor, and male gamete dispersal facilitates gene flow in lieu of natal dispersal of young.     

No gene flow across the Eastern Pacific Barrier in the reef‐building coral Porites lobata

The expanse of deep water between the central Pacific islands and the continental shelf of the Eastern Tropical Pacific is regarded as the world's most potent marine biogeographic barrier. During recurrent climatic fluctuations (ENSO, El Niño Southern Oscillation), however, changes in water temperature and the speed and direction of currents become favourable for trans‐oceanic dispersal of larvae from central Pacific to marginal eastern Pacific reefs. Here, we investigate the population connectivity of the reef‐building coral Porites lobata across the Eastern Pacific Barrier (EPB). Patterns of recent gene flow in samples (n = 1173) from the central Pacific and the Eastern Tropical Pacific (ETP) were analysed with 12 microsatellite loci. Results indicated that P. lobata from the ETP are strongly isolated from those in the central Pacific and Hawaii (  = 0.509; P < 0.001). However, samples from Clipperton Atoll, an oceanic island on the eastern side of the EPB, grouped with the central Pacific. Within the central Pacific, Hawaiian populations were strongly isolated from three co‐occurring clusters found throughout the remainder of the central Pacific. No further substructure was evident in the ETP. Changes in oceanographic conditions during ENSO over the past several thousand years thus appear insufficient to support larval deliveries from the central Pacific to the ETP or strong postsettlement selection acts on ETP settlers from the central Pacific. Recovery of P. lobata populations in the frequently disturbed ETP thus must depend on local larval sources.     

The genomic consequences of adaptive divergence and reproductive isolation between species of manakins

The processes of adaptation and speciation are expected to shape genomic variation within and between diverging species. Here we analyze genomic heterogeneity of genetic differentiation and introgression in a hybrid zone between two bird species (Manacus candei and M. vitellinus) using 59 100 SNPs, a whole genome assembly, and Bayesian models. Measures of genetic differentiation () and introgression (genomic cline center [α] and rate [β]) were highly heterogeneous among loci. We identified thousands of loci with elevated parameter estimates, some of which are likely to be associated with variation in fitness in Manacus populations. To analyze the genomic organization of differentiation and introgression, we mapped SNPs onto a draft assembly of the M. vitellinus genome. Estimates of , α, and β were autocorrelated at very short physical distances (< 100 bp), but much less so beyond this. In addition, average statistical associations (linkage disequilibrium) between SNPs were generally low and were not higher in admixed populations than in populations of the parental species. Although they did not occur with a constant probability across the genome, loci with elevated , α, and β were not strongly co-localized in the genome. Contrary to verbal models that predict clustering of loci involved in adaptation and isolation in discrete genomic regions, these results are consistent with the hypothesis that genetic regions involved in adaptive divergence and reproductive isolation are scattered throughout the genome. We also found that many loci were characterized by both exceptional genetic differentiation and introgression, consistent with the hypothesis that loci involved in isolation are also often characterized by a history of divergent selection. However, the concordance between isolation and differentiation was only partial, indicating a complex architecture and history of loci involved in isolation.     

Population genomics from pool sequencing

Next generation sequencing of pooled samples is an effective approach for studies of variability and differentiation in populations. In this paper we provide a comprehensive set of estimators of the most common statistics in population genetics based on the frequency spectrum, namely the Watterson estimator , nucleotide pairwise diversity Π, Tajima's D, Fu and Li's D and F, Fay and Wu's H, McDonald‐Kreitman and HKA tests and , corrected for sequencing errors and ascertainment bias. In a simulation study, we show that pool and individual θ estimates are highly correlated and discuss how the performance of the statistics vary with read depth and sample size in different evolutionary scenarios. As an application, we reanalyse sequences from Drosophila mauritiana and from an evolution experiment in Drosophila melanogaster. These methods are useful for population genetic projects with limited budget, study of communities of individuals that are hard to isolate, or autopolyploid species.     

, Jost's D,and FST are similarly constrained by allele frequencies: A mathematical,simulation, and empirical study

Statistics and Jost's D have been proposed for replacing F_ST as measures of genetic differentiation. A principal argument in favour of these statistics is the independence of their maximal values with respect to the subpopulation heterozygosity H_S, a property not shared by F_ST. Nevertheless, it has been unclear if these alternative differentiation measures are constrained by other aspects of the allele frequencies. Here, for biallelic markers, we study the mathematical properties of the maximal values of and D, comparing them to those of F_ST. We show that and D exhibit the same peculiar frequency‐dependence phenomena as F_ST, including a maximal value as a function of the frequency of the most frequent allele that lies well below one. Although the functions describing , D, and F_ST in terms of the frequency of the most frequent allele are different, the allele frequencies that maximize them are identical. Moreover, we show using coalescent simulations that when taking into account the specific maximal values of the three statistics, their behaviours become similar across a large range of migration rates. We use our results to explain two empirical patterns: the similar values of the three statistics among North American wolves, and the low D values compared to and F_ST in Atlantic salmon. The results suggest that the three statistics are often predictably similar, so that they can make quite similar contributions to data analysis. When they are not similar, the difference can be understood in relation to features of genetic diversity.     

Multispecies genetic structure and hybridization in the Betula genus across Eurasia

Boreal and cool temperate forests are the major land cover of northern Eurasia, and information about continental‐scale genetic structure and past demographic history of forest species is important from an evolutionary perspective and has conservation implications. However, although many population genetic studies of forest tree species have been conducted in Europe or Eastern Asia, continental‐scale genetic structure and past demographic history remain poorly known. Here, we focus on the birch genus Betula, which is commonly distributed in boreal and cool temperate forests, and examine 129 populations of two tetraploid and four diploid species collected from Iceland to Japan. All individuals were genotyped at seven to 18 nuclear simple sequence repeats (nSSRs). Pairwise among the six species ranged from 0.285 to 0.903, and genetic differentiation among them was clear. structure analysis suggested that Betula pubescens is an allotetraploid and one of the parental species was Betula pendula. In both species pairs of B. pendula and B. plathyphylla, and B. pubescens and B. ermanii, genetic diversity was highest in central Siberia. A hybrid zone was detected around Lake Baikal for eastern and western species pairs regardless of ploidy level. Approximate Bayesian computation suggested that the divergence of B. pendula and B. platyphylla occurred around the beginning of the last ice age (36 300 years BP, 95% CI: 15 330–92 700) and hybridization between them was inferred to have occurred after the last glacial maximum (1614 years BP, 95% CI: 561–4710), with B. pendula providing a higher contribution to hybrids.     

Genetic structure of the northwestern Russian wolf populations and gene flow between Russia and Finland

We examined the genetic diversity and structure of wolf populations in northwestern Russia. Populations in Republic of Karelia and Arkhangelsk Oblast were sampled during 1995–2000, and 43 individuals were genotyped with 10 microsatellite markers. Moreover, 118 previously genotyped wolves from the neighbouring Finnish population were used as a reference population. A relatively large amount of genetic variation was found in the Russian populations, and the Karelian wolf population tended to be slightly more polymorphic than the Arkhangelsk population. We found significant inbreeding (F = 0.094) in the Karelian, but not in the Arkhangelsk population. The effective size estimates of the Karelian wolf population based on the approximate Bayesian computation and linkage disequilibrium methods were 39.9 and 46.7 individuals, respectively. AMOVA-analysis and exact test of population differentiation suggested clear differentiation between the Karelian, Arkhangelsk and Finnish wolf populations. Indirect estimates of gene flow based on the level of population differentiation (ϕ _ST  = 0.152) and frequency of private alleles (0.029) both suggested a low level of gene flow between the populations (Nm = 1.4 and Nm = 3.7, respectively). Assignment analysis of Karelian and Finnish populations suggested an even lower number of recent migrants (less than 0.03) between populations, with a larger amount of migration from Finland to Karelia than vice versa. Our findings emphasise the role of physical obstacles and territorial behaviour in creating barriers to gene flow between populations in relatively limited geographical areas, even in large-bodied mammalian species with long-distance dispersal capabilities and an apparently continuous population structure.     

Long‐distance dispersal and genetic structure of natural populations: an assessment of the inverse isolation hypothesis in peat mosses

It is well accepted that the shape of the dispersal kernel, especially its tail, has a substantial effect on the genetic structure of species. Theory predicts that dispersal by fat‐tailed kernels reshuffles genetic material, and thus, preserves genetic diversity during colonization. Moreover, if efficient long‐distance dispersal is coupled with random colonization, an inverse isolation effect is predicted to develop in which increasing genetic diversity per colonizer is expected with increasing distance from a genetically variable source. By contrast, increasing isolation leads to decreasing genetic diversity when dispersal is via thin‐tailed kernels. Here, we use a well‐established model group for dispersal biology (peat mosses: genus Sphagnum) with a fat‐tailed dispersal kernel, and the natural laboratory of the Stockholm archipelago to study the validity of the inverse isolation hypothesis in spore‐dispersed plants in island colonization. Population genetic structure of three species (Sphagnum fallax, Sphagnum fimbriatum and Sphagnum palustre) with contrasting life histories and ploidy levels were investigated on a set of islands using microsatellites. Our data show (, amova , IBD) that dispersal of the two most abundant species can be well approximated by a random colonization model. We find that genetic diversity per colonizer on islands increases with distance from the mainland for S. fallax and S. fimbriatum. By contrast, S. palustre deviates from this pattern, owing to its restricted distribution in the region, affecting its source pool strength. Therefore, the inverse isolation effect appears to hold in natural populations of peat mosses and, likely, in other organisms with small diaspores.     

Relationship between effective population size,inbreeding and adult fitness‐related traits in a steelhead (Oncorhynchus mykiss) population released in the wild

Inbreeding is of concern in supportive breeding programmes in Pacific salmonids, Oncorhynchus spp, where the number of breeding adults is limited by rearing space or poor survival to adulthood, and large numbers are released to supplement wild stocks and fisheries. We reconstructed the pedigree of 6602 migratory hatchery steelhead (Oncorhynchus mykiss) over four generations, to determine the incidence and fitness consequences of inbreeding in a northwest USA programme. The hatchery maintained an effective population size,  = 107.9 from F₀ to F₂, despite an increasing census size (N), which resulted in a decreasing N_e/N ratio (0.35 in F₀ to 0.08 in F₂). The reduced ratio was attributed to a small broodstock size, nonrandom transfers and high variance in reproductive success (particularly in males). We observed accumulation of inbreeding from the founder generation (in F₄, percentage individuals with inbreeding coefficients Δf > 0 = 15.7%). Generalized linear mixed models showed that body length and weight decreased significantly with increasing Δf, and inbred fish returned later to spawn in a model that included father identity. However, there was no significant correlation between Δf and age at return, female fecundity or gonad weight. Similarly, there was no relationship between Δf and reproductive success of F₂ and F₃ individuals, which might be explained by the fact that reproductive success is partially controlled by hatchery mating protocols. This study is one of the first to show that small changes in inbreeding coefficient can affect some fitness‐related traits in a monitored population propagated and released to the wild.     

Environment not dispersal limitation drives clonal composition of Arctic Daphnia in a recently deglaciated area

One of the most prominent manifestations of the ongoing climate warming is the retreat of glaciers and ice sheets around the world. Retreating glaciers result in the formation of new ponds and lakes, which are available for colonization. The gradual appearance of these new habitat patches allows us to determine to what extent the composition of asexual Daphnia (water flea) populations is affected by environmental drivers vs. dispersal limitation. Here, we used a landscape genetics approach to assess the processes structuring the clonal composition of species in the D. pulex species complex that have colonized periglacial habitats created by ice‐sheet retreat in western Greenland. We analysed 61 populations from a young (<50 years) and an old cluster (>150 years) of lakes and ponds. We identified 42 asexual clones that varied widely in spatial distribution. Beta‐diversity was higher among older than among younger systems. Lineage sorting by the environment explained 14% of the variation in clonal composition whereas the pure effect of geographical distance was very small and statistically insignificant ( = 0.010, P = 0.085). Dispersal limitation did not seem important, even among young habitat patches. The observation of several tens of clones colonizing the area combined with environmentally driven clonal composition of populations illustrates that population assembly of asexual species in the Arctic is structured by environmental gradients reflecting differences in the ecology of clones.     

Association mapping in Salix viminalis L. (Salicaceae) – identification of candidate genes associated with growth and phenology

Willow species (Salix) are important as short‐rotation biomass crops for bioenergy, which creates a demand for faster genetic improvement and breeding through deployment of molecular marker‐assisted selection (MAS). To find markers associated with important adaptive traits, such as growth and phenology, for use in MAS, we genetically dissected the trait variation of a Salix viminalis (L.) population of 323 accessions. The accessions were sampled throughout northern Europe and were established at two field sites in Pustnäs, Sweden, and at Woburn, UK, offering the opportunity to assess the impact of genotype‐by‐environment interactions (G × E) on trait–marker associations. Field measurements were recorded for growth and phenology traits. The accessions were genotyped using 1536 SNP markers developed from phenology candidate genes and from genes previously observed to be differentially expressed in contrasting environments. Association mapping between 1233 of these SNPs and the measured traits was performed taking into account population structure and threshold selection bias. At a false discovery rate (FDR) of 0.2, 29 SNPs were associated with bud burst, leaf senescence, number of shoots or shoot diameter. The percentage of accession variation () explained by these associations ranged from 0.3% to 4.4%, suggesting that the studied traits are controlled by many loci of limited individual impact. Despite this, a SNP in the EARLY FLOWERING 3 gene was repeatedly associated (FDR < 0.2) with bud burst. The rare homozygous genotype exhibited 0.4–1.0 lower bud burst scores than the other genotype classes on a five‐grade scale. Consequently, this marker could be promising for use in MAS and the gene deserves further study. Otherwise, associations were less consistent across sites, likely due to their small estimates and to considerable G × E interactions indicated by multivariate association analyses and modest trait accession correlations across sites (0.32–0.61).     

Exploring a Pool‐seq‐only approach for gaining population genomic insights in nonmodel species

Developing genomic insights is challenging in nonmodel species for which resources are often scarce and prohibitively costly. Here, we explore the potential of a recently established approach using Pool‐seq data to generate a de novo genome assembly for mining exons, upon which Pool‐seq data are used to estimate population divergence and diversity. We do this for two pairs of sympatric populations of brown trout (Salmo trutta): one naturally sympatric set of populations and another pair of populations introduced to a common environment. We validate our approach by comparing the results to those from markers previously used to describe the populations (allozymes and individual‐based single nucleotide polymorphisms [SNPs]) and from mapping the Pool‐seq data to a reference genome of the closely related Atlantic salmon (Salmo salar). We find that genomic differentiation (F_ST) between the two introduced populations exceeds that of the naturally sympatric populations (F_ST = 0.13 and 0.03 between the introduced and the naturally sympatric populations, respectively), in concordance with estimates from the previously used SNPs. The same level of population divergence is found for the two genome assemblies, but estimates of average nucleotide diversity differ ( ≈ 0.002 and  ≈ 0.001 when mapping to S. trutta and S. salar, respectively), although the relationships between population values are largely consistent. This discrepancy might be attributed to biases when mapping to a haploid condensed assembly made of highly fragmented read data compared to using a high‐quality reference assembly from a divergent species. We conclude that the Pool‐seq‐only approach can be suitable for detecting and quantifying genome‐wide population differentiation, and for comparing genomic diversity in populations of nonmodel species where reference genomes are lacking.     

The effect of O2 and pressure on thiosulfate oxidation by Thiomicrospira thermophila

Microbial sulfur cycling in marine sediments often occurs in environments characterized by transient chemical gradients that affect both the availability of nutrients and the activity of microbes. High turnover rates of intermediate valence sulfur compounds and the intermittent availability of oxygen in these systems greatly impact the activity of sulfur‐oxidizing micro‐organisms in particular. In this study, the thiosulfate‐oxidizing hydrothermal vent bacterium Thiomicrospira thermophila strain EPR85 was grown in continuous culture at a range of dissolved oxygen concentrations (0.04–1.9 mM) and high pressure (5–10 MPa) in medium buffered at pH 8. Thiosulfate oxidation under these conditions produced tetrathionate, sulfate, and elemental sulfur, in contrast to previous closed‐system experiments at ambient pressure during which thiosulfate was quantitatively oxidized to sulfate. The maximum observed specific growth rate at 5 MPa pressure under unlimited O₂ was 0.25 hr^?1. This is comparable to the μ_max (0.28 hr^?1) observed at low pH (<6) at ambient pressure when T. thermophila produces the same mix of sulfur species. The half‐saturation constant for O₂ () estimated from this study was 0.2 mM (at a cell density of 10⁵ cells/ml) and was robust at all pressures tested (0.4–10 MPa), consistent with piezotolerant behavior of this strain. The cell‐specific was determined to be 1.5 pmol O₂/cell. The concentrations of products formed were correlated with oxygen availability, with tetrathionate production in excess of sulfate production at all pressure conditions tested. This study provides evidence for transient sulfur storage during times when substrate concentration exceeds cell‐specific and subsequent consumption when oxygen dropped below that threshold. These results may be common among sulfur oxidizers in a variety of environments (e.g., deep marine sediments to photosynthetic microbial mats).     

Fine‐scale frequency differentiation along a herbivory gradient in the trichome dimorphism of a wild Arabidopsis

Geographic variation is commonly observed in plant resistance traits, where plant species might experience different selection pressure across a heterogeneous landscape. Arabidopsis halleri subsp. gemmifera is dimorphic for trichome production, generating two morphs, trichome‐producing (hairy) and trichomeless (glabrous) plants. Trichomes of A. halleri are known to confer resistance against the white butterfly, cabbage sawfly, and brassica leaf beetle, but not against flea beetles. We combined leaf damage, microclimate, and microsatellite loci data of 26 A. halleri populations in central Japan, to explore factors responsible for fine‐scale geographic variation in the morph frequency. We found that hairy plants were less damaged than glabrous plants within populations, but the among‐site variation was the most significant source of variation in the individual‐level damage. Fixation index () of a putative trichome locus exhibited a significant divergence along population‐level damage with an exception of an outlier population, inferring the local adaptation to herbivory. Notably, this outlier was a population wherein our previous study reported a balancing role of the brassica leaf beetle Phaedon brassicae on the morph frequency. This differentiation of the trichome locus was unrelated to neutral genetic differentiation (evaluated by of microsatellite loci) and meteorological factors (including temperature and solar radiation). The present findings, combined with those of our previous work, provide suggestive evidence that herbivore‐driven divergence and occasional outbreak of a specific herbivore have jointly contributed to the ecogeographic pattern in the frequency of two morphs.