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.     

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.     

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 .     

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².     

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.     

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.     

High‐throughput identification of protein mutant stability computed from a double mutant fitness landscape

The effect of a mutation on protein stability is traditionally measured by genetic construction, expression, purification, and physical analysis using low‐throughput methods. This process is tedious and limits the number of mutants able to be examined in a single study. In contrast, functional fitness effects can be measured in a high‐throughput manner by various deep mutational scanning tools. Using protein GB 1, we have recently demonstrated the feasibility of estimating the mutational stability effect ( G) of single‐substitution based on the functional fitness profile of all double‐substitutions. The principle is to identify genetic backgrounds that have an exhausted stability margin. The functional effect of an additional substitution on these genetic backgrounds can then be used to compute the mutational G based on the biophysical relationship between functional fitness and thermodynamic stability. However, to identify such genetic backgrounds, the approach described in our previous study required a benchmark dataset, which is a set of known mutational G. In this study, a benchmark‐independent approach is developed. The genetic backgrounds of interest are identified using k‐means clustering with the integration of structural information. We further demonstrated that a reasonable approximation of G can also be obtained without taking structural information into account. In summary, this study describes a novel method for computing G from double‐substitution functional fitness profiles alone, without relying on any known mutational G as a benchmark.     

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.     

Deep soil carbon dynamics are driven more by soil type than by climate: a worldwide meta‐analysis of radiocarbon profiles

The response of soil carbon dynamics to climate and land‐use change will affect both the future climate and the quality of ecosystems. Deep soil carbon (>20 cm) is the primary component of the soil carbon pool, but the dynamics of deep soil carbon remain poorly understood. Therefore, radiocarbon activity (C), which is a function of the age of carbon, may help to understand the rates of soil carbon biodegradation and stabilization. We analyzed the published C contents in 122 profiles of mineral soil that were well distributed in most of the large world biomes, except for the boreal zone. With a multivariate extension of a linear mixed‐effects model whose inference was based on the parallel combination of two algorithms, the expectation–maximization (EM) and the Metropolis–Hasting algorithms, we expressed soil C profiles as a four‐parameter function of depth. The four‐parameter model produced insightful predictions of soil C as dependent on depth, soil type, climate, vegetation, land‐use and date of sampling (). Further analysis with the model showed that the age of topsoil carbon was primarily affected by climate and cultivation. By contrast, the age of deep soil carbon was affected more by soil taxa than by climate and thus illustrated the strong dependence of soil carbon dynamics on other pedologic traits such as clay content and mineralogy.     

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.     

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.     

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 chronic elevated nitrate concentrations on the structure and function of river biofilms

相似文献   

Characterization of the superoxide anion radical scavenging activity by tetracycline antibiotics in aprotic media

The tetracycline family antibiotics are widely used as human and veterinary treatments. The drugs are effective as antibiotics and also show antimicrobial and non‐microbial action. However, the antioxidant properties of tetracyclines have not been characterized in aprotic media. To better understand their biological functions, the in vitro superoxide anion radical () scavenging activities of tetracycline, chlortetracycline, oxytetracycline, doxycycline and methacycline were characterized, along with a very efficient scavenger, tiron, in dimethyl sulphoxide (DMSO), using ultra‐weak chemiluminescence (CL). We found that tetracycline, chlortetracycline and doxycycline efficiently inhibited CL from the ‐generating system at concentration levels of 0.02–1.0 mmol/L. Methacycline and oxytetracycline were the scavengers at concentration levels of 0.01–0.1 mmol/L, whereas when their concentration was lowered the drugs were capable of generating , leading to CL enhancement. For all the data obtained in this study, the scavenging activity for the compounds tested decreased in the following order: tetracycline > doxycycline > chlortetracycline > tiron methacycline > oxytetracycline. These results indicate that the tetracycline drugs directly alter redox chemistry in aprotic media. Copyright © 2011 John Wiley & Sons, Ltd.     

The role of a cytosolic superoxide dismutase in barley–pathogen interactions

Reactive oxygen species (ROS), including superoxide ( / ) and hydrogen peroxide (H₂O₂), are differentially produced during resistance responses to biotrophic pathogens and during susceptible responses to necrotrophic and hemi‐biotrophic pathogens. Superoxide dismutase (SOD) is responsible for the catalysis of the dismutation of / to H₂O₂, regulating the redox status of plant cells. Increased SOD activity has been correlated previously with resistance in barley to the hemi‐biotrophic pathogen Pyrenophora teres f. teres (Ptt, the causal agent of the net form of net blotch disease), but the role of individual isoforms of SOD has not been studied. A cytosolic CuZnSOD, HvCSD1, was isolated from barley and characterized as being expressed in tissue from different developmental stages. HvCSD1 was up‐regulated during the interaction with Ptt and to a greater extent during the resistance response. Net blotch disease symptoms and fungal growth were not as pronounced in transgenic HvCSD1 knockdown lines in a susceptible background (cv. Golden Promise), when compared with wild‐type plants, suggesting that cytosolic / contributes to the signalling required to induce a defence response to Ptt. There was no effect of HvCSD1 knockdown on infection by the hemi‐biotrophic rice blast pathogen Magnaporthe oryzae or the biotrophic powdery mildew pathogen Blumeria graminis f. sp. hordei, but HvCSD1 also played a role in the regulation of lesion development by methyl viologen. Together, these results suggest that HvCSD1 could be important in the maintenance of the cytosolic redox status and in the differential regulation of responses to pathogens with different lifestyles.     

Nucleotide‐dependent structural fluctuations and regulation of microtubule‐binding affinity of KIF1A

Molecular motors such as kinesin regulate affinity to a rail protein during the ATP hydrolysis cycle. The regulation mechanism, however, is yet to be determined. To understand this mechanism, we investigated the structural fluctuations of the motor head of the single‐headed kinesin called KIF1A in different nucleotide states using molecular dynamics simulations of a Gō‐like model. We found that the helix at the microtubule (MT) binding site intermittently exhibits a large structural fluctuation when MT is absent. Frequency of this fluctuation changes systematically according to the nucleotide states and correlates strongly with the experimentally observed binding affinity to MT. We also showed that thermal fluctuation enhances the correlation and the interaction with the nucleotide suppresses the fluctuation of the helix . These results suggest that KIF1A regulates affinity to MT by changing the flexibility of the helix during the ATP hydrolysis process: the binding site becomes more flexible in the strong binding state than in the weak binding state. Proteins 2015; 83:809–819. © 2015 Wiley Periodicals, Inc.     

Redox regulation of plant stem cell fate

Despite the importance of stem cells in plant and animal development, the common mechanisms of stem cell maintenance in both systems have remained elusive. Recently, the importance of hydrogen peroxide (H₂O₂) signaling in priming stem cell differentiation has been extensively studied in animals. Here, we show that different forms of reactive oxygen species (ROS) have antagonistic roles in plant stem cell regulation, which were established by distinct spatiotemporal patterns of ROS‐metabolizing enzymes. The superoxide anion () is markedly enriched in stem cells to activate WUSCHEL and maintain stemness, whereas H₂O₂ is more abundant in the differentiating peripheral zone to promote stem cell differentiation. Moreover, H₂O₂ negatively regulates biosynthesis in stem cells, and increasing H₂O₂ levels or scavenging leads to the termination of stem cells. Our results provide a mechanistic framework for ROS‐mediated control of plant stem cell fate and demonstrate that the balance between and H₂O₂ is key to stem cell maintenance and differentiation.     

Superoxide anion radical scavenging property of catecholamines

The direct effect of the four catecholamines (adrenaline, noradrenaline, dopamine and isoproterenol) on superoxide anion radicals () was investigated. The reaction between 18‐crown‐6‐ether and potassium superoxide in dimethylsulfoxide was used as a source of . The reactivity of catecholamines with was examined using chemiluminescence, reduction of nitroblue tetrazolium and electron paramagnetic resonance spin‐trapping techniques. 5,5‐Dimethyl‐1‐pyrroline‐N‐oxide was included as the spin trap. The results showed that the four catecholamines were effective and efficient in inhibiting chemiluminescence accompanying the potassium superoxide/18‐crown‐6‐ether system in a dose‐dependent manner over the range 0.05–2 mm in the following order: adrenaline > noradrenaline > dopamine > isoproterenol, with, IC₅₀ = 0.15 ± 0.02 mm 0.21 ± 0.03 mm , 0.27 ± 0.03 mm and 0.50 ± 0.04 mm , respectively. The catecholamines examined also exhibited a strong scavenging effect towards when evaluated this property by the inhibition of nitroblue tetrazolium reduction (56–73% at 1 m concentration). A very similar capacity of scavenging was monitored in the 5,5‐dimethyl‐1‐pyrroline‐N‐oxide spin‐trapping assay. The results suggest that catecholamines tested may involve a direct effect on scavenging radicals. Copyright © 2013 John Wiley & Sons, Ltd.