Soil carbon loss with warming: New evidence from carbon‐degrading enzymes

Climate warming affects soil carbon (C) dynamics, with possible serious consequences for soil C stocks and atmospheric CO₂ concentrations. However, the mechanisms underlying changes in soil C storage are not well understood, hampering long‐term predictions of climate C‐feedbacks. The activity of the extracellular enzymes ligninase and cellulase can be used to track changes in the predominant C sources of soil microbes and can thus provide mechanistic insights into soil C loss pathways. Here we show, using meta‐analysis, that reductions in soil C stocks with warming are associated with increased ratios of ligninase to cellulase activity. Furthermore, whereas long‐term (≥5 years) warming reduced the soil recalcitrant C pool by 14%, short‐term warming had no significant effect. Together, these results suggest that warming stimulates microbial utilization of recalcitrant C pools, possibly exacerbating long‐term climate‐C feedbacks.     

Evolutionary force of AT‐rich repeats to trap genomic and episomal DNAs into the rice genome: lessons from endogenous pararetrovirus

In plant genomes, the incorporation of DNA segments is not a common method of artificial gene transfer. Nevertheless, various segments of pararetroviruses have been found in plant genomes in recent decades. The rice genome contains a number of segments of endogenous rice tungro bacilliform virus‐like sequences (ERTBVs), many of which are present between AT dinucleotide repeats (ATrs). Comparison of genomic sequences between two closely related rice subspecies, japonica and indica, allowed us to verify the preferential insertion of ERTBVs into ATrs. In addition to ERTBVs, the comparative analyses showed that ATrs occasionally incorporate repeat sequences including transposable elements, and a wide range of other sequences. Besides the known genomic sequences, the insertion sequences also represented DNAs of unclear origins together with ERTBVs, suggesting that ATrs have integrated episomal DNAs that would have been suspended in the nucleus. Such insertion DNAs might be trapped by ATrs in the genome in a host‐dependent manner. Conversely, other simple mono‐ and dinucleotide sequence repeats (SSR) were less frequently involved in insertion events relative to ATrs. Therefore, ATrs could be regarded as hot spots of double‐strand breaks that induce non‐homologous end joining. The insertions within ATrs occasionally generated new gene‐related sequences or involved structural modifications of existing genes. Likewise, in a comparison between Arabidopsis thaliana and Arabidopsis lyrata, the insertions preferred ATrs to other SSRs. Therefore ATrs in plant genomes could be considered as genomic dumping sites that have trapped various DNA molecules and may have exerted a powerful evolutionary force.     

Identification and sequence analysis of two thaumatin‐like protein (TmTLP) genes from Tenebrio molitor

Thaumatin‐like proteins (TLPs) share structural similarity with the sweet‐tasting thaumatin protein but exhibit antifungal activity by inhibiting growth of fungal pathogens. In a Tenebrio model, two TLP genes were identified by RNA‐sequencing analysis and genome sequencing. Both TmTLP1 and TmTLP2 genes contain 729 nucleotide sequences encoding 242 amino acid residues, including an initiation codon (ATG) and a termination codon (TAA). Interestingly, TmTLPs are proteins with 14 central cysteine residues that may have an important role in structure formation. These data will be used to characterize the innate immune functions of TmTLPs in Tenebrio molitor.     

ABC inference of multi‐population divergence with admixture from unphased population genomic data

Rapidly developing sequencing technologies and declining costs have made it possible to collect genome‐scale data from population‐level samples in nonmodel systems. Inferential tools for historical demography given these data sets are, at present, underdeveloped. In particular, approximate Bayesian computation (ABC) has yet to be widely embraced by researchers generating these data. Here, we demonstrate the promise of ABC for analysis of the large data sets that are now attainable from nonmodel taxa through current genomic sequencing technologies. We develop and test an ABC framework for model selection and parameter estimation, given histories of three‐population divergence with admixture. We then explore different sampling regimes to illustrate how sampling more loci, longer loci or more individuals affects the quality of model selection and parameter estimation in this ABC framework. Our results show that inferences improved substantially with increases in the number and/or length of sequenced loci, while less benefit was gained by sampling large numbers of individuals. Optimal sampling strategies given our inferential models included at least 2000 loci, each approximately 2 kb in length, sampled from five diploid individuals per population, although specific strategies are model and question dependent. We tested our ABC approach through simulation‐based cross‐validations and illustrate its application using previously analysed data from the oak gall wasp, Biorhiza pallida.     

Genetic barcoding of dark‐spored myxomycetes (Amoebozoa)—Identification,evaluation and application of a sequence similarity threshold for species differentiation in NGS studies

Unicellular, eukaryotic organisms (protists) play a key role in soil food webs as major predators of microorganisms. However, due to the polyphyletic nature of protists, no single universal barcode can be established for this group, and the structure of many protistean communities remains unresolved. Plasmodial slime moulds (Myxogastria or Myxomycetes) stand out among protists by their formation of fruit bodies, which allow for a morphological species concept. By Sanger sequencing of a large collection of morphospecies, this study presents the largest database to date of dark‐spored myxomycetes and evaluate a partial 18S SSU gene marker for species annotation. We identify and discuss the use of an intraspecific sequence similarity threshold of 99.1% for species differentiation (OTU picking) in environmental PCR studies (ePCR) and estimate a hidden diversity of putative species, exceeding those of described morphospecies by 99%. When applying the identified threshold to an ePCR data set (including sequences from both NGS and cloning), we find 64 OTUs of which 21.9% had a direct match (>99.1% similarity) to the database and the remaining had on average 90.2 ± 0.8% similarity to their best match, thus thought to represent undiscovered diversity of dark‐spored myxomycetes.     

A pragmatic approach to the analysis of diets of generalist predators: the use of next‐generation sequencing with no blocking probes

Predicting whether a predator is capable of affecting the dynamics of a prey species in the field implies the analysis of the complete diet of the predator, not simply rates of predation on a target taxon. Here, we employed the Ion Torrent next‐generation sequencing technology to investigate the diet of a generalist arthropod predator. A complete dietary analysis requires the use of general primers, but these will also amplify the predator unless suppressed using a blocking probe. However, blocking probes can potentially block other species, particularly if they are phylogenetically close. Here, we aimed to demonstrate that enough prey sequence could be obtained without blocking probes. In communities with many predators, this approach obviates the need to design and test numerous blocking primers, thus making analysis of complex community food webs a viable proposition. We applied this approach to the analysis of predation by the linyphiid spider Oedothorax fuscus in an arable field. We obtained over two million raw reads. After discarding the low‐quality and predator reads, the libraries still contained over 61 000 prey reads (3% of the raw reads; 6% of reads passing quality control). The libraries were rich in Collembola, Lepidoptera, Diptera and Nematoda. They also contained sequences derived from several spider species and from horticultural pests (aphids). Oedothorax fuscus is common in UK cereal fields, and the results showed that it is exploiting a wide range of prey. Next‐generation sequencing using general primers but without blocking probes provided ample sequences for analysis of the prey range of this spider and proved to be a simple and inexpensive approach.     

The association of six non‐synonymous variants in three DNA repair genes with hepatocellular carcinoma risk: a meta‐analysis

Hepatocellular carcinoma is a complex polygenic disease. Despite the huge advances in genetic epidemiology, it still remains a challenge to unveil the genetic architecture of hepatocellular carcinoma. We, therefore, decided to meta‐analytically assess the association of six non‐synonymous coding variants from XRCC1, XRCC3 and XPD genes with hepatocellular carcinoma risk by pooling the results of 20 English articles. This meta‐analysis was conducted according to the PRISMA statement, and data collection was independently completed in duplicate. In overall analyses, the minor alleles of four variants, Arg280His (odds ratio, 95% confidence interval, P: 1.37, 1.13–1.66, 0.001), Thr241Met (1.93, 1.17–3.20, 0.011), Asp312Asn (1.22, 1.08–1.38, 0.001) and Lys751Gln (1.42, 1.02–1.97, 0.038), were associated with the significant risk for hepatocellular carcinoma. There were low probabilities of publication bias for all variants. Subgroup analyses revealed significant association of XRCC1 gene Arg399Gln with hepatocellular carcinoma in Chinese especially from south China (odds ratio, 95% confidence interval, P: 1.57, 1.16–2.14, 0.004), in larger studies (1.48, 1.11–1.98, 0.007) and in studies with population‐based controls (1.33, 1.06–1.68, 0.016). Taken together, our findings demonstrated that XPD gene Asp312Asn and XRCC1 gene Arg399Gln might be candidate susceptibility loci for hepatocellular carcinoma. Considering the ubiquity of genetic heterogeneity, further validation in a broad range of ethnic populations is warranted.     

Linking genotype to phenotype in a changing ocean: inferring the genomic architecture of a blue mussel stress response with genome‐wide association

A key component to understanding the evolutionary response to a changing climate is linking underlying genetic variation to phenotypic variation in stress response. Here, we use a genome‐wide association approach (GWAS) to understand the genetic architecture of calcification rates under simulated climate stress. We take advantage of the genomic gradient across the blue mussel hybrid zone (Mytilus edulis and Mytilus trossulus) in the Gulf of Maine (GOM) to link genetic variation with variance in calcification rates in response to simulated climate change. Falling calcium carbonate saturation states are predicted to negatively impact many marine organisms that build calcium carbonate shells – like blue mussels. We sampled wild mussels and measured net calcification phenotypes after exposing mussels to a ‘climate change’ common garden, where we raised temperature by 3°C, decreased pH by 0.2 units and limited food supply by filtering out planktonic particles >5 μm, compared to ambient GOM conditions in the summer. This climate change exposure greatly increased phenotypic variation in net calcification rates compared to ambient conditions. We then used regression models to link the phenotypic variation with over 170 000 single nucleotide polymorphism loci (SNPs) generated by genotype by sequencing to identify genomic locations associated with calcification phenotype, and estimate heritability and architecture of the trait. We identified at least one of potentially 2–10 genomic regions responsible for 30% of the phenotypic variation in calcification rates that are potential targets of natural selection by climate change. Our simulations suggest a power of 13.7% with our study's average effective sample size of 118 individuals and rare alleles, but a power of >90% when effective sample size is 900.     

Cofactor‐binding sites in proteins of deviating sequence: Comparative analysis and clustering in torsion angle,cavity, and fold space

Small molecules are recognized in protein‐binding pockets through surface‐exposed physicochemical properties. To optimize binding, they have to adopt a conformation corresponding to a local energy minimum within the formed protein–ligand complex. However, their conformational flexibility makes them competent to bind not only to homologous proteins of the same family but also to proteins of remote similarity with respect to the shape of the binding pockets and folding pattern. Considering drug action, such observations can give rise tounexpected and undesired cross reactivity. In this study, datasets of six different cofactors (ADP, ATP, NAD(P)(H), FAD, and acetyl CoA, sharing an adenosine diphosphate moiety as common substructure), observed in multiple crystal structures of protein–cofactor complexes exhibiting sequence identity below 25%, have been analyzed for the conformational properties of the bound ligands, the distribution of physicochemical properties in the accommodating protein‐binding pockets, and the local folding patterns next to the cofactor‐binding site. State‐of‐the‐art clustering techniques have been applied to group the different protein–cofactor complexes in the different spaces. Interestingly, clustering in cavity (Cavbase) and fold space (DALI) reveals virtually the same data structuring. Remarkable relationships can be found among the different spaces. They provide information on how conformations are conserved across the host proteins and which distinct local cavity and fold motifs recognize the different portions of the cofactors. In those cases, where different cofactors are found to be accommodated in a similar fashion to the same fold motifs, only a commonly shared substructure of the cofactors is used for the recognition process. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Carbohydrate recognition by the rhamnose‐binding lectin SUL‐I with a novel three‐domain structure isolated from the venom of globiferous pedicellariae of the flower sea urchin Toxopneustes pileolus

The globiferous pedicellariae of the venomous sea urchin Toxopneustes pileolus contains several biologically active proteins. We have cloned the cDNA of one of the toxin components, SUL‐I, which is a rhamnose‐binding lectin (RBL) that acts as a mitogen through binding to carbohydrate chains on target cells. Recombinant SUL‐I (rSUL‐I) was produced in Escherichia coli cells, and its carbohydrate‐binding specificity was examined with the glycoconjugate microarray analysis, which suggested that potential target carbohydrate structures are galactose‐terminated N‐glycans. rSUL‐I exhibited mitogenic activity for murine splenocyte cells and toxicity against Vero cells. The three‐dimensional structure of the rSUL‐I/l ‐rhamnose complex was determined by X‐ray crystallographic analysis at a 1.8 Å resolution. The overall structure of rSUL‐I is composed of three distinctive domains with a folding structure similar to those of CSL3, a RBL from chum salmon (Oncorhynchus keta) eggs. The bound l ‐rhamnose molecules are mainly recognized by rSUL‐I through hydrogen bonds between its 2‐, 3‐, and 4‐hydroxy groups and Asp, Asn, and Glu residues in the binding sites, while Tyr and Ser residues participate in the recognition mechanism. It was also inferred that SUL‐I may form a dimer in solution based on the molecular size estimated via dynamic light scattering as well as possible contact regions in its crystal structure.     

An evaluation of behavior inferences from Bayesian state‐space models: A case study with the Pacific walrus

State‐space models offer researchers an objective approach to modeling complex animal location data sets, and state‐space model behavior classifications are often assumed to have a link to animal behavior. In this study, we evaluated the behavioral classification accuracy of a Bayesian state‐space model in Pacific walruses using Argos satellite tags with sensors to detect animal behavior in real time. We fit a two‐state discrete‐time continuous‐space Bayesian state‐space model to data from 306 Pacific walruses tagged in the Chukchi Sea. We matched predicted locations and behaviors from the state‐space model (resident, transient behavior) to true animal behavior (foraging, swimming, hauled out) and evaluated classification accuracy with kappa statistics (κ) and root mean square error (RMSE). In addition, we compared biased random bridge utilization distributions generated with resident behavior locations to true foraging behavior locations to evaluate differences in space use patterns. Results indicated that the two‐state model fairly classified true animal behavior (0.06 ≤ κ ≤ 0.26, 0.49 ≤ RMSE ≤ 0.59). Kernel overlap metrics indicated utilization distributions generated with resident behavior locations were generally smaller than utilization distributions generated with true foraging behavior locations. Consequently, we encourage researchers to carefully examine parameters and priors associated with behaviors in state‐space models, and reconcile these parameters with the study species and its expected behaviors.     

Increased salt tolerance with overexpression of cation/proton antiporter 1 genes: a meta‐analysis

Cation/proton antiporter 1 (CPA1) genes encode cellular Na⁺/H⁺ exchanger proteins, which act to adjust ionic balance. Overexpression of CPA1s can improve plant performance under salt stress. However, the diversified roles of the CPA1 family and the various parameters used in evaluating transgenic plants over‐expressing CPA1s make it challenging to assess the complex functions of CPA1s and their physiological mechanisms in salt tolerance. Using meta‐analysis, we determined how overexpression of CPA1s has influenced several plant characteristics involved in response and resilience to NaCl stress. We also evaluated experimental variables that favour or reduce CPA1 effects in transgenic plants. Viewed across studies, overexpression of CPA1s has increased the magnitude of 10 of the 19 plant characteristics examined, by 25% or more. Among the ten moderating variables, several had substantial impacts on the extent of CPA1 influence: type of culture media, donor and recipient type and genus, and gene family. Genes from monocotyledonous plants stimulated root K⁺, root K⁺/Na⁺, total chlorophyll, total dry weight and root length much more than genes from dicotyledonous species. Genes transformed to or from Arabidopsis have led to smaller CPA1‐induced increases in plant characteristics than genes transferred to or from other genera. Heterogeneous expression of CPA1s led to greater increases in leaf chlorophyll and root length than homologous expression. These findings should help guide future investigations into the function of CPA1s in plant salt tolerance and the use of genetic engineering for breeding of resistance.     

Variants of the IL‐10 gene associate with muscle strength in elderly from rural Africa: a candidate gene study

Recently, it has been shown that the capacity of the innate immune system to produce cytokines relates to skeletal muscle mass and strength in older persons. The interleukin‐10 (IL‐10) gene regulates the production capacities of IL‐10 and tumour necrosis factor‐α (TNF‐α). In rural Ghana, IL‐10 gene variants associated with different production capacities of IL‐10 and TNF‐α are enriched compared with Caucasian populations. In this setting, we explored the association between these gene variants and muscle strength. Among 554 Ghanaians aged 50 years and older, we determined 20 single nucleotide polymorphisms in the IL‐10 gene, production capacities of IL‐10 and TNF‐α in whole blood upon stimulation with lipopolysaccharide (LPS) and handgrip strength as a proxy for skeletal muscle strength. We distinguished pro‐inflammatory haplotypes associated with low IL‐10 production capacity and anti‐inflammatory haplotypes with high IL‐10 production capacity. We found that distinct haplotypes of the IL‐10 gene associated with handgrip strength. A pro‐inflammatory haplotype with a population frequency of 43.2% was associated with higher handgrip strength (P = 0.015). An anti‐inflammatory haplotype with a population frequency of 7.9% was associated with lower handgrip strength (P = 0.006). In conclusion, variants of the IL‐10 gene contributing to a pro‐inflammatory cytokine response associate with higher muscle strength, whereas those with anti‐inflammatory response associate with lower muscle strength. Future research needs to elucidate whether these effects of variation in the IL‐10 gene are exerted directly through its role in the repair of muscle tissue or indirectly through its role in the defence against infectious diseases.     

Inferring the geographic origin of a range expansion: Latitudinal and longitudinal coordinates inferred from genomic data in an ABC framework with the program x‐origin

Climatic or environmental change is not only driving distributional shifts in species today, but it has also caused distributions to expand and contract in the past. Inferences about the geographic locations of past populations especially regions that served as refugia (i.e., source populations) and migratory routes are a challenging endeavour. Refugial areas may be evidenced from fossil records or regions of temporal stability inferred from ecological niche models. Genomic data offer an alternative and broadly applicable source of information about the locality of refugial areas, especially relative to fossil data, which are either unavailable or incomplete for most species. Here, we present a pipeline we developed (called x ‐origin ) for statistically inferring the geographic origin of range expansion using a spatially explicit coalescent model and an approximate Bayesian computation testing framework. In addition to assessing the probability of specific latitudinal and longitudinal coordinates of refugial or source populations, such inferences can also be made accounting for the effects of temporal and spatial environmental heterogeneity, which may impact migration routes. We demonstrate x ‐origin with an analysis of genomic data collected in the Collared pika that underwent postglacial expansion across Alaska, as well as present an assessment of its accuracy under a known model of expansion to validate the approach.