Methods to account for spatial autocorrelation in the analysis of species distributional data: a review

Species distributional or trait data based on range map (extent‐of‐occurrence) or atlas survey data often display spatial autocorrelation, i.e. locations close to each other exhibit more similar values than those further apart. If this pattern remains present in the residuals of a statistical model based on such data, one of the key assumptions of standard statistical analyses, that residuals are independent and identically distributed (i.i.d), is violated. The violation of the assumption of i.i.d. residuals may bias parameter estimates and can increase type I error rates (falsely rejecting the null hypothesis of no effect). While this is increasingly recognised by researchers analysing species distribution data, there is, to our knowledge, no comprehensive overview of the many available spatial statistical methods to take spatial autocorrelation into account in tests of statistical significance. Here, we describe six different statistical approaches to infer correlates of species’ distributions, for both presence/absence (binary response) and species abundance data (poisson or normally distributed response), while accounting for spatial autocorrelation in model residuals: autocovariate regression; spatial eigenvector mapping; generalised least squares; (conditional and simultaneous) autoregressive models and generalised estimating equations. A comprehensive comparison of the relative merits of these methods is beyond the scope of this paper. To demonstrate each method's implementation, however, we undertook preliminary tests based on simulated data. These preliminary tests verified that most of the spatial modeling techniques we examined showed good type I error control and precise parameter estimates, at least when confronted with simplistic simulated data containing spatial autocorrelation in the errors. However, we found that for presence/absence data the results and conclusions were very variable between the different methods. This is likely due to the low information content of binary maps. Also, in contrast with previous studies, we found that autocovariate methods consistently underestimated the effects of environmental controls of species distributions. Given their widespread use, in particular for the modelling of species presence/absence data (e.g. climate envelope models), we argue that this warrants further study and caution in their use. To aid other ecologists in making use of the methods described, code to implement them in freely available software is provided in an electronic appendix.     

Striatal pre- and postsynaptic profile of adenosine A(2A) receptor antagonists

Striatal adenosine A(2A) receptors (A(2A)Rs) are highly expressed in medium spiny neurons (MSNs) of the indirect efferent pathway, where they heteromerize with dopamine D(2) receptors (D(2)Rs). A(2A)Rs are also localized presynaptically in cortico-striatal glutamatergic terminals contacting MSNs of the direct efferent pathway, where they heteromerize with adenosine A(1) receptors (A(1)Rs). It has been hypothesized that postsynaptic A(2A)R antagonists should be useful in Parkinson's disease, while presynaptic A(2A)R antagonists could be beneficial in dyskinetic disorders, such as Huntington's disease, obsessive-compulsive disorders and drug addiction. The aim or this work was to determine whether selective A(2A)R antagonists may be subdivided according to a preferential pre- versus postsynaptic mechanism of action. The potency at blocking the motor output and striatal glutamate release induced by cortical electrical stimulation and the potency at inducing locomotor activation were used as in vivo measures of pre- and postsynaptic activities, respectively. SCH-442416 and KW-6002 showed a significant preferential pre- and postsynaptic profile, respectively, while the other tested compounds (MSX-2, SCH-420814, ZM-241385 and SCH-58261) showed no clear preference. Radioligand-binding experiments were performed in cells expressing A(2A)R-D(2)R and A(1)R-A(2A)R heteromers to determine possible differences in the affinity of these compounds for different A(2A)R heteromers. Heteromerization played a key role in the presynaptic profile of SCH-442416, since it bound with much less affinity to A(2A)R when co-expressed with D(2)R than with A(1)R. KW-6002 showed the best relative affinity for A(2A)R co-expressed with D(2)R than co-expressed with A(1)R, which can at least partially explain the postsynaptic profile of this compound. Also, the in vitro pharmacological profile of MSX-2, SCH-420814, ZM-241385 and SCH-58261 was is in accordance with their mixed pre- and postsynaptic profile. On the basis of their preferential pre- versus postsynaptic actions, SCH-442416 and KW-6002 may be used as lead compounds to obtain more effective antidyskinetic and antiparkinsonian compounds, respectively.     

Fungal antagonists of the plant pathogen Rhizoctonia solani: selection, control efficacy and influence on the indigenous microbial community

A broad spectrum of fungal antagonists was evaluated as potential biocontrol agents (BCAs) against the soil-borne pathogen Rhizoctonia solani using a new combination of in vitro and in vivo assays. The in vitro characterisation of diverse parameters including the ability to parasitise mycelium and to inhibit the germination of Rhizoctonia sclerotia at different temperatures resulted in the selection of six potential fungal antagonists. These were genotypically characterised by their BOX-PCR fingerprints, and identified as Trichoderma reesei and T. viride by partial 18S rDNA sequencing. When potato sprouts were treated with Trichoderma, all isolates significantly reduced the incidence of Rhizoctonia symptoms. Evaluated under growth chamber conditions, the selected Trichoderma isolates either partly or completely controlled the dry mass loss of lettuce caused by R. solani. Furthermore, the antagonistic Trichoderma strains were active under field conditions. To analyse the effect of Trichoderma treatment on indigenous root-associated microbial communities, we performed a DNA-dependent SSCP (Single-Strand Conformation Polymorphism) analysis of 16S rDNA/ITS sequences. In this first assessment study for Trichoderma it was shown that the pathogen and the vegetation time had much more influence on the composition of the microbiota than the BCA treatment. After evaluation of all results, three Trichoderma strains originally isolated from Rhizoctonia sclerotia were selected as promising BCAs.     

Decimal Place Slope,A Fast and Precise Method for Quantifying 13C Incorporation Levels for Detecting the Metabolic Activity of Microbial Species

The metabolic incorporation of stable isotopes such as ¹³C or ¹⁵N into proteins has become a powerful tool for qualitative and quantitative proteome studies. We recently introduced a method that monitors heavy isotope incorporation into proteins and presented data revealing the metabolic activity of various species in a microbial consortium using this technique. To further develop our method using an liquid chromatography (LC)-mass spectrometry (MS)-based approach, we present here a novel approach for calculating the incorporation level of ¹³C into peptides by using the information given in the decimal places of peptide masses obtained by modern high-resolution MS. In the present study, the applicability of this approach is demonstrated using Pseudomonas putida ML2 proteins uniformly labeled via the consumption of [¹³C₆]benzene present in the medium at concentrations of 0, 10, 25, 50, and 100 atom %. The incorporation of ¹³C was calculated on the basis of several labeled peptides derived from one band on an SDS-PAGE gel. The accuracy of the calculated incorporation level depended upon the number of peptide masses included in the analysis, and it was observed that at least 100 peptide masses were required to reduce the deviation below 4 atom %. This accuracy was comparable with calculations of incorporation based on the isotope envelope. Furthermore, this method can be extended to the calculation of the labeling efficiency for a wide range of biomolecules, including RNA and DNA. The technique will therefore allow a highly accurate determination of the carbon flux in microbial consortia with a direct approach based solely on LC-MS.The metabolic incorporation of stable isotopes such as ¹³C or ¹⁵N into proteins has become a powerful component of qualitative and quantitative proteome studies (1). Incorporation of heavy isotopes can be used to analyze microbial processes such as turnover rates and also to help to establish structure-function relationships within microbial communities. Stable isotope probing (SIP¹) techniques based on DNA-SIP (2) and RNA-SIP (3) have been used for this purpose previously. With the introduction of protein-SIP (4), the need for an accurate alternative method for calculating label incorporation into biomolecules arose. Protein-SIP has several advantages compared with DNA/RNA-SIP, the most important being its capacity to detect dynamic levels of incorporation, whereas only labeled or unlabeled states can be categorized by means of DNA/RNA-SIP because of the need to separate ¹³C-DNA/RNA by density gradient centrifugation. Quantitative analysis of ¹³C incorporation levels is of the utmost importance, especially when unraveling carbon fluxes through either microbial communities or food webs with different trophic levels.In contrast to the incorporation of isotopically labeled amino acids, which is often used in quantitative proteomics (5), metabolic labeling by growth substrates and nutrients (e.g. salts) is often imperfect and makes the processing of mass spectrometry (MS) data difficult. For example, when the incorporation of ¹³C exceeds ∼2 atom %, common database search algorithms fail to identify peptides and proteins. The problem can only be managed successfully if a stable, known degree of ¹³C incorporation can be achieved during the experiment (6). Using a low labeling efficiency of roughly 5 atom %, Huttlin et al. (6) chose the altered envelope chain for calculating the incorporation and simultaneously used the signal intensity for a quantitative comparison with the sample that had a natural abundance of ¹³C. Database approaches for peptide identification can cope only with the natural abundance of carbon isotopes; they fail if the incorporation of ¹³C significantly exceeds the natural isotope abundance or if incorporation patterns occur in unpredictable ways (7).The simplest method for determining the incorporation level is to compare the unlabeled average mass of the monoisotopic peptide with the mass of the labeled protein, as estimated by matrix-assisted laser desorption/ionization or electrospray ionization MS (8, 9). A more advanced approach for determining the isotopic mass distribution of peptides is based on the isotopic distribution of the peaks of a peptide envelope (10, 11). Here, for a given isotopomer, the incorporation efficiency is defined as the percentage of incorporated ¹³C atoms with relation to the total number of carbon atoms with the natural isotope abundance (approximately 1.01 atom % ¹³C). As a reference, the theoretical isotopic distribution of a peptide is calculated based upon an algorithm described elsewhere (12). The isotope distribution of both unlabeled and labeled peptides can subsequently be used to calculate the incorporation level. For this method, an Excel spreadsheet (ProSIPQuant.xls) was developed (4). A similar approach, also based on the calculation of isotopic distributions, has been used in other studies (7). In these studies, however, the identification of the peptides is limited to those that have unlabeled counterparts; in addition, an exact calculation can be hampered by overlapping signals coming from additional peaks with similar masses.In the present study, we describe a new way of determining the isotope incorporation level. Our method makes use of characteristic patterns in the digits after the decimal point of the peptide masses generated by high-accuracy instruments such as the linear ion trap LTQ-Orbitrap (Thermo Fisher Scientific, Bremen, Germany). For tryptic peptides, typical regularities in the decimal places of the monoisotopic masses have been observed (13, 14). These observations have been explored in detail for theoretical and experimental data of proteins originating from Helicobacter pylori (15). As a result, a rule called the “half decimal place rule” (HDPR) was defined; it states that the decimal place is nearly half of the first digit for tryptic peptides with masses in the range of 500–1,000 Da. In other words, the exact mass of a peptide is equal to its nominal mass times ∼1.005. Because the difference between ¹²C and ¹³C is slightly greater than 1 Da, exactly 1.0033548378, the decimal places of a tryptic peptide''s mass are shifted in a regular manner by the incorporation level and lead to a significantly increased slope for the digits in the third and fourth place after the decimal point. This shift can be used to estimate the incorporation level of heavy isotopes into the protein. Detecting such shifts requires the highly accurate measurement possible with modern mass spectrometers such as the LTQ-Orbitrap, the Fourier transform ion cyclotron resonance, or the quadrupole time of flight. In this communication, we demonstrate the applicability of this approach using Pseudomonas putida ML2 proteins labeled uniformly via the consumption of [¹³C₆]benzene with five different substrate concentrations (0, 10, 25, 50, and 100 atom % of ¹³C). The ¹³C incorporation was calculated based on several labeled peptides derived from different proteins in one SDS-PAGE band. By these means, we have established a method that allows the determination of ¹³C incorporation into proteins and can be used to assess the metabolic activity of a given species within a mixed community.     

Genome mining reveals the evolutionary origin and biosynthetic potential of basidiomycete polyketide synthases

Numerous polyketides are known from bacteria, plants, and fungi. However, only a few have been isolated from basidiomycetes. Large scale genome sequencing projects now help anticipate the capacity of basidiomycetes to synthesize polyketides. In this study, we identified and annotated 111 type I and three type III polyketide synthase (PKS) genes from 35 sequenced basidiomycete genomes. Phylogenetic analysis of PKS genes suggests that all main types of fungal iterative PKS had already evolved before the Ascomycota and Basidiomycota diverged. A comparison of genomic and metabolomic data shows that the number of polyketide genes exceeds the number of known polyketide structures by far. Exploiting these results to design degenerate PCR primers, we amplified and cloned the complete sequence of armB, a PKS gene from the melleolide producer Armillaria mellea. We expect this study will serve as a guide for future genomic mining projects to discover structurally diverse mushroom-derived polyketides.     

Molecular characterization of phytoplasmas in lilies with fasciation in the Czech Republic

Lilium spp. with symptoms of severe fasciation were observed in Southern and central Bohemia during the period 1999-2003. Nucleic acids extracted from symptomatic and asymptomatic plants were used in nested-PCR assays with primers amplifying 16S-23S rRNA sequences specific for phytoplasmas. The subsequent nested-PCR with phytoplasma group-specific primers followed by RFLP analyses and the 16S ribosomal gene sequencing, allowed classification of the detected phytoplasmas in the aster yellows group, subgroups 16SrI-B and 16SrI-C alone, and in mixed infection. Samples infected by 16SrI-C phytoplasmas showed different overlapping RFLP profiles after TruI digestion of R16F2/R2 amplicons. Two of these amplicons were sequenced, one of them directly and the other after cloning; sequence analyses and blast alignment confirmed the presence of two different overlapping patterns in samples studied. The sequences obtained were closely related, respectively, to operon A and operon B ribosomal sequences of the clover phyllody phytoplasma. Direct PCR followed by RFLP analyses of the tuf gene with two restriction enzymes showed no differences from reference strain of subgroup 16SrI-C. Infection with aster yellows phytoplasmas of 16SrI-B subgroup in asymptomatic lilies cv. Sunray was also detected.     

Induced Variations in Brassinosteroid Genes Define Barley Height and Sturdiness,and Expand the Green Revolution Genetic Toolkit

Reduced plant height and culm robustness are quantitative characteristics important for assuring cereal crop yield and quality under adverse weather conditions. A very limited number of short-culm mutant alleles were introduced into commercial crop cultivars during the Green Revolution. We identified phenotypic traits, including sturdy culm, specific for deficiencies in brassinosteroid biosynthesis and signaling in semidwarf mutants of barley (Hordeum vulgare). This set of characteristic traits was explored to perform a phenotypic screen of near-isogenic short-culm mutant lines from the brachytic, breviaristatum, dense spike, erectoides, semibrachytic, semidwarf, and slender dwarf mutant groups. In silico mapping of brassinosteroid-related genes in the barley genome in combination with sequencing of barley mutant lines assigned more than 20 historic mutants to three brassinosteroid-biosynthesis genes (BRASSINOSTEROID-6-OXIDASE, CONSTITUTIVE PHOTOMORPHOGENIC DWARF, and DIMINUTO) and one brassinosteroid-signaling gene (BRASSINOSTEROID-INSENSITIVE1 [HvBRI1]). Analyses of F2 and M2 populations, allelic crosses, and modeling of nonsynonymous amino acid exchanges in protein crystal structures gave a further understanding of the control of barley plant architecture and sturdiness by brassinosteroid-related genes. Alternatives to the widely used but highly temperature-sensitive uzu1.a allele of HvBRI1 represent potential genetic building blocks for breeding strategies with sturdy and climate-tolerant barley cultivars.The introduction of dwarfing genes to increase culm sturdiness of cereal crops was crucial for the first Green Revolution (Hedden, 2003). The culms of tall cereal crops were not strong enough to support the heavy spikes of high-yielding cultivars, especially under high-nitrogen conditions. As a result, plants fell over, a process known as lodging. This caused losses in yield and grain-quality issues attributable to fungal infections, mycotoxin contamination, and preharvest germination (Rajkumara, 2008). Today, a second Green Revolution is on its way, to revolutionize the agricultural sector and to ensure food production for a growing world population. Concurrently, global climate change is expected to cause more frequent occurrences of extreme weather conditions, including thunderstorms with torrential rain and strong winds, thus promoting cereal culm breakage (Porter and Semenov, 2005; National Climate Assessment Development Advisory Committee, 2013). Accordingly, plant architectures that resist lodging remain a major crop-improvement goal and identification of genes that regulate culm length is required to enhance the genetic toolbox in order to facilitate efficient marker-assisted breeding. The mutations and the corresponding genes that enabled the Green Revolution in wheat (Triticum aestivum) and rice (Oryza sativa) have been identified (Hedden, 2003). They all relate to gibberellin metabolism and signal transduction. It is now known that other plant hormones such as brassinosteroids are also involved in the regulation of plant height. Knowledge of the molecular mechanisms underlying the effects of the two hormones on cell elongation and division has mainly come from studies in Arabidopsis (Arabidopsis thaliana; Bai et al., 2012). Mutant-based breeding strategies to fine-tune brassinosteroid metabolism and signaling pathways could improve lodging behavior in modern crops (Vriet et al., 2012) such as barley (Hordeum vulgare), which is the fourth most abundant cereal in both area and tonnage harvested (http://faostat.fao.org).A short-culm phenotype in crops is often accompanied by other phenotypic changes. Depending on the penetrance of such pleiotropic characters, but also the parental background and different scientific traditions and expertise, short-culmed barley mutants were historically divided into groups, such as brachytic (brh), breviaristatum (ari), dense spike (dsp), erectoides (ert), semibrachytic (uzu), semidwarf (sdw), or slender dwarf (sld; Franckowiak and Lundqvist, 2012). Subsequent mutant characterization was limited to intragroup screens and very few allelism tests between mutants from different groups have been reported (Franckowiak and Lundqvist, 2012). Although the total number of short-culm barley mutants exceeds 500 (Franckowiak and Lundqvist, 2012), very few have been characterized at the DNA level (Helliwell et al., 2001; Jia et al., 2009; Chandler and Harding, 2013; Houston et al., 2013). One of the first identified haplotypes was uzu barley (Chono et al., 2003). The Uzu1 gene encodes the brassinosteroid hormone receptor and is orthologous to the BRASSINOSTEROID-INSENSITIVE1 (BRI1) gene of Arabidopsis, a crucial promoter of plant growth (Li and Chory, 1997). The uzu1.a allele has been used in East Asia for over a century and is presently distributed in winter barley cultivars in Japan, the Korean peninsula, and China (Saisho et al., 2004). Its agronomic importance comes from the short and sturdy culm that provides lodging resistance, and an upright plant architecture that tolerates dense planting.Today, more than 50 different brassinosteroids have been identified in plants (Bajguz and Tretyn, 2003). Most are intermediates of the complex biosynthetic pathway (Shimada et al., 2001). Approximately nine genes code for the enzymes that participate in the biosynthetic pathway from episterol to brassinolide (Supplemental Fig. S1). Brassinosteroid deficiency is caused by down-regulation of these genes, but it can also be associated with brassinosteroid signaling. The first protein in the signaling network is the brassinosteroid receptor encoded by BRI1 (Li and Chory, 1997; Kim and Wang, 2010). In this work, we show how to visually identify brassinosteroid-mutant barley plants and we describe more than 20 relevant mutations in four genes of the brassinosteroid biosynthesis and signaling pathways that can be used in marker-assisted breeding strategies.     

What processes shape early-successional vegetation in fly ash and mine tailings?

Local abiotic filters and regional processes (i.e., regional pools of species that are dispersal-limited to varying degrees) interactively structure the development of vegetation in human-disturbed habitats, yet their relative contributions to this process are still to be determined. In this study conducted in the Czech Republic, we related plant species diversity and composition of 10 fly ash and 7 mine tailings to local edaphic conditions, and to vegetation from a 100-m perimeter (regional species pool). We found that the species richness and composition on the tailings were significantly associated with diversity and composition of vegetation in the surroundings, but not with the local edaphic conditions. Species from adjacent vegetation that were more abundant and those producing lighter seeds were more likely to establish on the tailings. The same characteristics also enhanced species abundance on the tailings, but the two predictors explained less than 10% in variation of establishment success or of species abundance. A non-significant relationship between species number and tailings size, but a significant association between diversity and time of vegetation development indicate that the study systems are still far from equilibrium. Our study provides evidence for a strong effect of regional processes, with a limited influence of measured edaphic conditions on plant communities developing de novo. It also highlights the necessity to consider the broader spatial context in the analysis of vegetation succession in human-disturbed habitats.     

Meiotic chromosomes in motion: a perspective from Mus musculus and Caenorhabditis elegans

Vigorous chromosome movement during the extended prophase of the first meiotic division is conserved in most eukaryotes. The movement is crucial for the faithful segregation of homologous chromosomes into daughter cells, and thus for fertility. A prerequisite for meiotic chromosome movement is the stable and functional attachment of telomeres or chromosome ends to the nuclear envelope and their cytoplasmic coupling to the cytoskeletal forces responsible for generating movement. Important advances in understanding the components, mechanisms, and regulation of chromosome end attachment and movement have recently been made. This review focuses on insights gained from experiments into two major metazoan model organisms: the mouse, Mus musculus, and the nematode, Caenorhabditis elegans.

     

Ovarian carcinoma CDK12 mutations misregulate expression of DNA repair genes via deficient formation and function of the Cdk12/CycK complex

The Cdk12/CycK complex promotes expression of a subset of RNA polymerase II genes, including those of the DNA damage response. CDK12 is among only nine genes with recurrent somatic mutations in high-grade serous ovarian carcinoma. However, the influence of these mutations on the Cdk12/CycK complex and their link to cancerogenesis remain ill-defined. Here, we show that most mutations prevent formation of the Cdk12/CycK complex, rendering the kinase inactive. By examining the mutations within the Cdk12/CycK structure, we find that they likely provoke structural rearrangements detrimental to Cdk12 activation. Our mRNA expression analysis of the patient samples containing the CDK12 mutations reveals coordinated downregulation of genes critical to the homologous recombination DNA repair pathway. Moreover, we establish that the Cdk12/CycK complex occupies these genes and promotes phosphorylation of RNA polymerase II at Ser2. Accordingly, we demonstrate that the mutant Cdk12 proteins fail to stimulate the faithful DNA double strand break repair via homologous recombination. Together, we provide the molecular basis of how mutated CDK12 ceases to function in ovarian carcinoma. We propose that CDK12 is a tumor suppressor of which the loss-of-function mutations may elicit defects in multiple DNA repair pathways, leading to genomic instability underlying the genesis of the cancer.