Contemporary effective population and metapopulation size (Ne and meta‐Ne): comparison among three salmonids inhabiting a fragmented system and differing in gene flow and its asymmetries

We estimated local and metapopulation effective sizes ( and meta‐) for three coexisting salmonid species (Salmo salar, Salvelinus fontinalis, Salvelinus alpinus) inhabiting a freshwater system comprising seven interconnected lakes. First, we hypothesized that might be inversely related to within‐species population divergence as reported in an earlier study (i.e., F_ST: S. salar> S. fontinalis> S. alpinus). Using the approximate Bayesian computation method implemented in ONeSAMP, we found significant differences in () between species, consistent with a hierarchy of adult population sizes (). Using another method based on a measure of linkage disequilibrium (LDNE: ), we found more finite values for S. salar than for the other two salmonids, in line with the results above that indicate that S. salar exhibits the lowest among the three species. Considering subpopulations as open to migration (i.e., removing putative immigrants) led to only marginal and non‐significant changes in , suggesting that migration may be at equilibrium between genetically similar sources. Second, we hypothesized that meta‐ might be significantly smaller than the sum of local s (null model) if gene flow is asymmetric, varies among subpopulations, and is driven by common landscape features such as waterfalls. One ‘bottom‐up’ or numerical approach that explicitly incorporates variable and asymmetric migration rates showed this very pattern, while a number of analytical models provided meta‐ estimates that were not significantly different from the null model or from each other. Our study of three species inhabiting a shared environment highlights the importance and utility of differentiating species‐specific and landscape effects, not only on dispersal but also in the demography of wild populations as assessed through local s and meta‐s and their relevance in ecology, evolution and conservation.     

Estimators for QST and coalescence times

Comparisons of to can provide insights into the evolutionary processes that lead to differentiation, or lack thereof, among the phenotypes of different groups (e.g., populations, species), and these comparisons have been performed on a variety of taxa, including humans. Here, I show that for neutrally evolving (i.e., by genetic drift, mutation, and gene flow alone) quantitative characters, the two commonly used estimators have somewhat different interpretations in terms of coalescence times, particularly when the number of groups that have been sampled is small. A similar situation occurs for estimators. Consequently, when observations come from only a small number of groups, which is not an unusual situation, it is important to match estimators appropriately when comparing to .     

Soil moisture variations affect short‐term plant‐microbial competition for ammonium,glycine, and glutamate

We tested whether the presence of plant roots would impair the uptake of ammonium (), glycine, and glutamate by microorganisms in a deciduous forest soil exposed to constant or variable moisture in a short‐term (24‐h) experiment. The uptake of ¹⁵NH₄ and dual labeled amino acids by the grass Festuca gigantea L. and soil microorganisms was determined in planted and unplanted soils maintained at 60% WHC (water holding capacity) or subject to drying and rewetting. The experiment used a design by which competition was tested in soils that were primed by plant roots to the same extent in the planted and unplanted treatments. Festuca gigantea had no effect on microbial N uptake in the constant moist soil, but its presence doubled the microbial uptake in the dried and rewetted soil compared with the constant moist. The drying and rewetting reduced by half or more the uptake by F. gigantea, despite more than 60% increase in the soil concentration of . At the same time, the amino acid and ‐ N became equally valued in the plant uptake, suggesting that plants used amino acids to compensate for the lower acquisition. Our results demonstrate the flexibility in plant‐microbial use of different N sources in response to soil moisture fluctuations and emphasize the importance of including transient soil conditions in experiments on resource competition between plants and soil microorganisms. Competition between plants and microorganisms for N is demonstrated by a combination of removal of one of the potential competitors, the plant, and subsequent observations of the uptake of N in the organisms in soils that differ only in the physical presence and absence of the plant during a short assay. Those conditions are necessary to unequivocally test for competition.     

Identifying important conservation areas for the clouded leopard Neofelis nebulosa in a mountainous landscape: Inference from spatial modeling techniques

The survival of large carnivores is increasingly precarious due to extensive human development that causes the habitat loss and fragmentation. Habitat selection is influenced by anthropogenic as well as environmental factors, and understanding these relationships is important for conservation management. We assessed the environmental and anthropogenic variables that influence site use of clouded leopard Neofelis nebulosa in Bhutan, estimated their population density, and used the results to predict the species’ site use across Bhutan. We used a large camera‐trap dataset from the national tiger survey to estimate for clouded leopards, for the first time in Bhutan, (1) population density using spatially explicit capture–recapture models and (2) site‐use probability using occupancy models accounting for spatial autocorrelation. Population density was estimated at (0.10 SD) and (0.12 SE) per 100 km². Clouded leopard site use was positively associated with forest cover and distance to river while negatively associated with elevation. Mean site‐use probability (from the Bayesian spatial model) was (0.076 SD). When spatial autocorrelation was ignored, the probability of site use was overestimated, (0.066 SD). Predictive mapping allowed us to identify important conservation areas and priority habitats to sustain the future of these elusive, ambassador felids and associated guilds. Multiple sites in the south, many of them outside of protected areas, were identified as habitats suitable for this species, adding evidence to conservation planning for clouded leopards in continental South Asia.     

Effect of initial soil properties on six‐year growth of 15 tree species in tropical restoration plantings

In restoration plantings in degraded pastures, initial soil nutrient status may lead to differential growth of tropical tree species with diverse life history attributes and capacity for N₂ fixation. In 2006, we planted 1,440 seedlings of 15 native tree species in 16 fenced plots (30 × 30 m) in a 60‐year‐old pasture in Los Tuxtlas, Veracruz, Mexico, in two planting combinations. In the first year, we evaluated bulk density, pH, the concentration of organic carbon (C), total nitrogen (N), ammonia (), nitrate (), and total phosphorus (P) in the upper soil profile (0–20 cm in depth) of all plots. The first two axes of two principal component analyses explained more than 60% of the variation in soil variables: The axes were related to increasing bulk density, , , total N concentration, and pH. Average relative growth rates in diameter at the stem base of the juvenile trees after 6 years were higher for pioneer (45.7%) and N₂‐fixing species (47.6%) than for nonpioneer (34.7%) and nonfixing species (36.2%). Most N₂‐fixing species and those with the slowest growth rates did not respond to soil attributes. Tree species benefited from higher pH levels and existing litter biomass. The pioneers Ficus yoponensis, Cecropia obtusifolia, and Heliocarpus appendiculatus, and the N₂‐fixing nonpioneers Cojoba arborea, Inga sinacae, and Platymiscium dimorphandrum were promising for forest restoration on our site, given their high growth rates.     

Nitrate or ammonium: Influences of nitrogen source on the physiology of a green alga

In freshwaters, algal species are exposed to different inorganic nitrogen (N_i) sources whose incorporation varies in biochemical energy demand. We hypothesized that due to the lesser energy requirement of ammonium ()‐use, in contrast to nitrate ()‐use, more energy remains for other metabolic processes, especially under CO₂‐ and phosphorus (P_i) limiting conditions. Therefore, we tested differences in cell characteristics of the green alga Chlamydomonas acidophila grown on or under covariation of CO₂ and P_i‐supply in order to determine limitations, in a full‐factorial design. As expected, results revealed higher carbon fixation rates for ‐grown cells compared to growth with under low CO₂ conditions. ‐grown cells accumulated more of the nine analyzed amino acids, especially under P_i‐limited conditions, compared to cells provided with . This is probably due to a slower protein synthesis in cells provided with . In contrast to our expectations, compared to ‐grown cells ‐grown cells had higher photosynthetic efficiency under P_i‐limitation. In conclusion, growth on the N_i‐source did not result in a clearly enhanced C_i‐assimilation, as it was highly dependent on P_i and CO₂ conditions (replete or limited). Results are potentially connected to the fact that C. acidophila is able to use only CO₂ as its inorganic carbon (C_i) source.     

Oxygen and carbon isoscapes for the Baltic Sea: Testing their applicability in fish migration studies

Conventional tags applied to individuals have been used to investigate animal movement, but these methods require tagged individuals be recaptured. Maps of regional isotopic variability known as “isoscapes” offer potential for various applications in migration research without tagging wherein isotope values of tissues are compared to environmental isotope values. In this study, we present the spatial variability in oxygen () and dissolved inorganic carbon (δ¹³C_DIC) isotope values of Baltic Sea water. We also provide an example of how these isoscapes can reveal locations of individual animal via spatial probability surface maps, using the high‐resolution salmon otolith isotope data from salmon during their sea‐feeding phase in the Baltic Sea. A clear latitudinal and vertical gradient was found for both and δ¹³C_DIC values. The difference between summer and winter in the Baltic Sea values was only slight, whereas δ¹³C_DIC values exhibited substantial seasonal variability related to algal productivity. Salmon otolith δ¹⁸O_oto and δ¹³C_oto values showed clear differences between feeding areas and seasons. Our example demonstrates that dual isotope approach offers great potential for estimating probable fish habitats once issues in model parameterization have been resolved.     

Effective size of density‐dependent populations in fluctuating environments

Reliable estimates of effective population size are of central importance in population genetics and evolutionary biology. For populations that fluctuate in size, harmonic mean population size is commonly used as a proxy for (multi‐) generational effective size. This assumes no effects of density dependence on the ratio between effective and actual population size, which limits its potential application. Here, we introduce density dependence on vital rates in a demographic model of variance effective size. We derive an expression for the ratio in a density‐regulated population in a fluctuating environment. We show by simulations that yearly genetic drift is accurately predicted by our model, and not proportional to as assumed by the harmonic mean model, where N is the total population size of mature individuals. We find a negative relationship between and N. For a given N, the ratio depends on variance in reproductive success and the degree of resource limitation acting on the population growth rate. Finally, our model indicate that environmental stochasticity may affect not only through fluctuations in N, but also for a given N at a given time. Our results show that estimates of effective population size must include effects of density dependence and environmental stochasticity.     

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.     

Pyric‐carnivory: Raptor use of prescribed fires

Fire is a process that shaped and maintained most terrestrial ecosystems worldwide. Changes in land use and patterns of human settlement have altered fire regimes and led to fire suppression resulting in numerous undesirable consequences spanning individual species and entire ecosystems. Many obvious and direct consequences of fire suppression have been well studied, but several, albeit less obvious, costs of alteration to fire regimes on wildlife are unknown. One such phenomenon is the response of carnivores to fire events—something we refer to as pyric‐carnivory. To investigate the prevalence of pyric‐carnivory in raptors, we monitored 25 prescribed fires occurring during two different seasons and across two different locations in tallgrass prairie of the central United States. We used paired point counts occurring before and during prescribed fires to quantify the use of fires by raptors. We found a strong attraction to fires with average maximum abundance nearly seven times greater during fires than prior to ignitions (before:  = 2.90, SE = 0.42; during:  = 20.20; SE = 3.29) and an average difference between fire events and immediately before fires of 15.2 (±2.69) raptors. This result was driven by Swainson's hawks (Buteo swainsoni), which were the most abundant (n = 346) of the nine species we observed using fires. Our results illustrate the importance of fire as integral disturbance process that effects wildlife behavior through multiple mechanisms that are often overshadowed by the predominant view of fire as a tool used for vegetation management.     

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.     

speed‐ne: Software to simulate and estimate genetic effective population size (Ne) from linkage disequilibrium observed in single samples

The genetic effective population size, N_e, can be estimated from the average gametic disequilibrium () between pairs of loci, but such estimates require evaluation of assumptions and currently have few methods to estimate confidence intervals. speed‐ne is a suite of matlab computer code functions to estimate from with a graphical user interface and a rich set of outputs that aid in understanding data patterns and comparing multiple estimators. speed‐ne includes functions to either generate or input simulated genotype data to facilitate comparative studies of estimators under various population genetic scenarios. speed‐ne was validated with data simulated under both time‐forward and time‐backward coalescent models of genetic drift. Three classes of estimators were compared with simulated data to examine several general questions: what are the impacts of microsatellite null alleles on , how should missing data be treated, and does disequilibrium contributed by reduced recombination among some loci in a sample impact . Estimators differed greatly in precision in the scenarios examined, and a widely employed estimator exhibited the largest variances among replicate data sets. speed‐ne implements several jackknife approaches to estimate confidence intervals, and simulated data showed that jackknifing over loci and jackknifing over individuals provided ~95% confidence interval coverage for some estimators and should be useful for empirical studies. speed‐ne provides an open‐source extensible tool for estimation of from empirical genotype data and to conduct simulations of both microsatellite and single nucleotide polymorphism (SNP) data types to develop expectations and to compare estimators.     

Effects of distance on detectability of Arctic waterfowl using double‐observer sampling during helicopter surveys

Aerial survey is an important, widely employed approach for estimating free‐ranging wildlife over large or inaccessible study areas. We studied how a distance covariate influenced probability of double‐observer detections for birds counted during a helicopter survey in Canada’s central Arctic. Two observers, one behind the other but visually obscured from each other, counted birds in an incompletely shared field of view to a distance of 200 m. Each observer assigned detections to one of five 40‐m distance bins, guided by semi‐transparent marks on aircraft windows. Detections were recorded with distance bin, taxonomic group, wing‐flapping behavior, and group size. We compared two general model‐based estimation approaches pertinent to sampling wildlife under such situations. One was based on double‐observer methods without distance information, that provide sampling analogous to that required for mark–recapture (MR) estimation of detection probability, , and group abundance, , along a fixed‐width strip transect. The other method incorporated double‐observer MR with a categorical distance covariate (MRD). A priori, we were concerned that estimators from MR models were compromised by heterogeneity in due to un‐modeled distance information; that is, more distant birds are less likely to be detected by both observers, with the predicted effect that would be biased high, and biased low. We found that, despite increased complexity, MRD models (ΔAICc range: 0–16) fit data far better than MR models (ΔAICc range: 204–258). However, contrary to expectation, the more naïve MR estimators of were biased low in all cases, but only by 2%–5% in most cases. We suspect that this apparently anomalous finding was the result of specific limitations to, and trade‐offs in, visibility by observers on the survey platform used. While MR models provided acceptable point estimates of group abundance, their far higher stranded errors (0%–40%) compared to MRD estimates would compromise ability to detect temporal or spatial differences in abundance. Given improved precision of MRD models relative to MR models, and the possibility of bias when using MR methods from other survey platforms, we recommend avian ecologists use MRD protocols and estimation procedures when surveying Arctic bird populations.     

Evaluating indices of body condition in two cricket species

Body mass components (dry mass, lean dry mass, water mass, fat mass) in each sex correlate strongly with body mass and pronotum length in Gryllus texensis and Acheta domesticus. Ordinary least squares (OLS) regression underestimates the scaling relationship between body mass and structural size (i.e., pronotum length) in both cricket species compared with standard major axis (SMA) regression. Standardized mass components correlate more strongly with scaled mass index () than with residual body mass (R_i). R_i represents the residuals from an OLS regression of log body mass against log pronotum length. Neither condition index predicts energy stores (i.e., fat content) in G. texensis. R_i is not correlated with energy stores in A. domesticus whereas is negatively correlated. A comparison of condition index methods using published data showed that neither sex nor diet quality affected body condition at adulthood in G. texensis when using the scaled mass index. However, the residual index suggested that sex had a significant effect on body condition. Further, analysis of covariance (ANCOVA) suggested that diet quality significantly affects body mass while statistically controlling for body size (i.e., body condition). We conclude that the statistical assumptions of condition index methods must be met prior to use and urge caution when using methods that are based on least squares in the y‐plane (i.e., residual index ANCOVA).     

Summary goodness‐of‐fit statistics for binary generalized linear models with noncanonical link functions

Generalized linear models (GLM) with a canonical logit link function are the primary modeling technique used to relate a binary outcome to predictor variables. However, noncanonical links can offer more flexibility, producing convenient analytical quantities (e.g., probit GLMs in toxicology) and desired measures of effect (e.g., relative risk from log GLMs). Many summary goodness‐of‐fit (GOF) statistics exist for logistic GLM. Their properties make the development of GOF statistics relatively straightforward, but it can be more difficult under noncanonical links. Although GOF tests for logistic GLM with continuous covariates (GLMCC) have been applied to GLMCCs with log links, we know of no GOF tests in the literature specifically developed for GLMCCs that can be applied regardless of link function chosen. We generalize the Tsiatis GOF statistic originally developed for logistic GLMCCs, (), so that it can be applied under any link function. Further, we show that the algebraically related Hosmer–Lemeshow () and Pigeon–Heyse (J²) statistics can be applied directly. In a simulation study, , , and J² were used to evaluate the fit of probit, log–log, complementary log–log, and log models, all calculated with a common grouping method. The statistic consistently maintained Type I error rates, while those of and J² were often lower than expected if terms with little influence were included. Generally, the statistics had similar power to detect an incorrect model. An exception occurred when a log GLMCC was incorrectly fit to data generated from a logistic GLMCC. In this case, had more power than or J².     

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.     

Different fates of deposited and in a temperate forest in northeast China: a 15N tracer study

Increasing atmospheric reactive nitrogen (N) deposition due to human activities could change N cycling in terrestrial ecosystems. However, the differences between the fates of deposited and are still not fully understood. Here, we investigated the fates of deposited and , respectively, via the application of ¹⁵NH₄NO₃ and NH₄¹⁵NO₃ in a temperate forest ecosystem. Results showed that at 410 days after tracer application, most was immobilized in litter layer (50 ± 2%), while a considerable amount of penetrated into 0–5 cm mineral soil (42 ± 2%), indicating that litter layer and 0–5 cm mineral soil were the major N sinks of and , respectively. Broad‐leaved trees assimilated more ¹⁵N under NH₄¹⁵NO₃ treatment compared to under ¹⁵NH₄NO₃ treatment, indicating their preference for –N. At 410 days after tracer application, 16 ± 4% added ¹⁵N was found in aboveground biomass under treatment, which was twice more than that under treatment (6 ± 1%). At the same time, approximately 80% added ¹⁵N was recovered in soil and plants under both treatments, which suggested that this forest had high potential for retention of deposited N. These results provided evidence that there were great differences between the fates of deposited and , which could help us better understand the mechanisms and capability of forest ecosystems as a sink of reactive nitrogen.     

pKa predictions for proteins,RNAs, and DNAs with the Gaussian dielectric function using DelPhi pKa

We developed a Poisson‐Boltzmann based approach to calculate the values of protein ionizable residues (Glu, Asp, His, Lys and Arg), nucleotides of RNA and single stranded DNA. Two novel features were utilized: the dielectric properties of the macromolecules and water phase were modeled via the smooth Gaussian‐based dielectric function in DelPhi and the corresponding electrostatic energies were calculated without defining the molecular surface. We tested the algorithm by calculating values for more than 300 residues from 32 proteins from the PPD dataset and achieved an overall RMSD of 0.77. Particularly, the RMSD of 0.55 was achieved for surface residues, while the RMSD of 1.1 for buried residues. The approach was also found capable of capturing the large shifts of various single point mutations in staphylococcal nuclease (SNase) from ‐cooperative dataset, resulting in an overall RMSD of 1.6 for this set of pKa's. Investigations showed that predictions for most of buried mutant residues of SNase could be improved by using higher dielectric constant values. Furthermore, an option to generate different hydrogen positions also improves predictions for buried carboxyl residues. Finally, the calculations on two RNAs demonstrated the capability of this approach for other types of biomolecules. Proteins 2015; 83:2186–2197. © 2015 Wiley Periodicals, Inc.     

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.