Temporal variability of ecological niches: a study on intertidal macrobenthic fauna

The determination of temporal niche dynamics under field conditions is an important component of a species’ ecology. Recent developments in niche mapping, and the possibility to account for spatial autocorrelation in species distributions, hold promise for the statistical approach explored here. Using species counts from a landscape‐scale benthic monitoring programme in the western Dutch Wadden Sea during 1997–2005 in combination with sediment characteristics and tidal height as explanatory variables, we statistically derive realised niches for two bivalves, two crustaceans and three polychaetes, encompassing predators, suspension and bottom feeding functional groups. Unsurprisingly, realized niches varied considerably between species. Intraspecific temporal variation was assessed as overlap between the year‐specific niche and the overall mean niche, and this analysis revealed considerable variation between years. The main functional groups represented by these species showed idiosyncratic and wide variability through the study period. There were no strong associations between niche characteristics and mean abundance or body size. Our assessment of intraspecific niche variability has ramifications for species distribution models in general and offers advances from previous methods. 1) By assessing species’ realized niches in the multivariate environmental space, analyses are independent from the relative availability of particular environments. Predicted realized niches present differences between years, rather than annual differences in environmental conditions. 2) Using spatially explicit models to predict species habitat preferences provide more precise and unbiased estimates of species–environment relationships. 3) Current niche models assume constant niches, whereas we illustrate how much these can vary over only a few generations. This emphasizes the potentially limited scope of global change studies with forecasts based on single‐time species distribution snapshots.     

Combined monte carlo and molecular dynamics simulation of fully hydrated dioleyl and palmitoyl-oleyl phosphatidylcholine lipid bilayers

We have applied a new equilibration procedure for the atomic level simulation of a hydrated lipid bilayer to hydrated bilayers of dioleyl-phosphatidylcholine (DOPC) and palmitoyl-oleyl phosphatidylcholine (POPC). The procedure consists of alternating molecular dynamics trajectory calculations in a constant surface tension and temperature ensemble with configurational bias Monte Carlo moves to different regions of the configuration space of the bilayer in a constant volume and temperature ensemble. The procedure is applied to bilayers of 128 molecules of POPC with 4628 water molecules, and 128 molecules of DOPC with 4825 water molecules. Progress toward equilibration is almost three times as fast in central processing unit (CPU) time compared with a purely molecular dynamics (MD) simulation. Equilibration is complete, as judged by the lack of energy drift in 200-ps runs of continuous MD. After the equilibrium state was reached, as determined by agreement between the simulation volume per lipid molecule with experiment, continuous MD was run in an ensemble in which the lateral area was restrained to fluctuate about a mean value and a pressure of 1 atm applied normal to the bilayer surface. Three separate continuous MD runs, 200 ps in duration each, separated by 10,000 CBMC steps, were carried out for each system. Properties of the systems were calculated and averaged over the three separate runs. Results of the simulations are presented and compared with experimental data and with other recent simulations of POPC and DOPC. Analysis of the hydration environment in the headgroups supports a mechanism by which unsaturation contributes to reduced transition temperatures. In this view, the relatively horizontal orientation of the unsaturated bond increases the area per lipid, resulting in increased water penetration between the headgroups. As a result the headgroup-headgroup interactions are attenuated and shielded, and this contributes to the lowered transition temperature.     

Expansion of mass-flowering crops leads to transient pollinator dilution and reduced wild plant pollination

Agricultural land use results in direct biodiversity decline through loss of natural habitat, but may also cause indirect cross-habitat effects on conservation areas. We conducted three landscape-scale field studies on 67 sites to test the hypothesis that mass flowering of oilseed rape (Brassica napus) results in a transient dilution of bees in crop fields, and in increased competition between crop plants and grassland plants for pollinators. Abundances of bumble-bees, which are the main pollinators of the grassland plant Primula veris, but also pollinate oilseed rape (OSR), decreased with increasing amount of OSR. This landscape-scale dilution affected bumble-bee abundances strongly in OSR fields and marginally in grasslands, where bumble-bee abundances were generally low at the time of Primula flowering. Seed set of Primula veris, which flowers during OSR bloom, was reduced by 20 per cent when the amount of OSR within 1 km radius increased from 0 to 15 per cent. Hence, the current expansion of bee-attractive biofuel crops results in transient dilution of crop pollinators, which means an increased competition for pollinators between crops and wild plants. In conclusion, mass-flowering crops potentially threaten fitness of concurrently flowering wild plants in conservation areas, despite the fact that, in the long run, mass-flowering crops can enhance abundances of generalist pollinators and their pollination service.     

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.     

Squares of different sizes: effect of geographical projection on model parameter estimates in species distribution modeling

In species distribution analyses, environmental predictors and distribution data for large spatial extents are often available in long‐lat format, such as degree raster grids. Long‐lat projections suffer from unequal cell sizes, as a degree of longitude decreases in length from approximately 110 km at the equator to 0 km at the poles. Here we investigate whether long‐lat and equal‐area projections yield similar model parameter estimates, or result in a consistent bias. We analyzed the environmental effects on the distribution of 12 ungulate species with a northern distribution, as models for these species should display the strongest effect of projectional distortion. Additionally we choose four species with entirely continental distributions to investigate the effect of incomplete cell coverage at the coast. We expected that including model weights proportional to the actual cell area should compensate for the observed bias in model coefficients, and similarly that using land coverage of a cell should decrease bias in species with coastal distribution. As anticipated, model coefficients were different between long‐lat and equal‐area projections. Having progressively smaller and a higher number of cells with increasing latitude influenced the importance of parameters in models, increased the sample size for the northernmost parts of species ranges, and reduced the subcell variability of those areas. However, this bias could be largely removed by weighting long‐lat cells by the area they cover, and marginally by correcting for land coverage. Overall we found little effect of using long‐lat rather than equal‐area projections in our analysis. The fitted relationship between environmental parameters and occurrence probability differed only very little between the two projection types. We still recommend using equal‐area projections to avoid possible bias. More importantly, our results suggest that the cell area and the proportion of a cell covered by land should be used as a weight when analyzing distribution of terrestrial species.     

Proteolytic processing of TAR DNA binding protein-43 by caspases produces C-terminal fragments with disease defining properties independent of progranulin

Neuronal and glial deposition of misfolded, proteolytically processed, polyubiquitinated and abnormally phosphorylated C-terminal fragments (CTFs) of the TAR DNA binding protein-43 (TDP-43) is a pathological hallmark of frontotemporal lobar degeneration with ubiquitin positive inclusions (FTLD-U) and certain cases of amyotrophic lateral sclerosis. We demonstrate that TDP-43 can be proteolytically processed by caspases upon induction of apoptosis to a major 35 kDa and a minor 25 kDa CTF. These fragments are initially soluble, but over time they accumulate as insoluble and pathologically phosphorylated derivatives. However, proteolytic processing appears not to be absolutely required for the deposition of insoluble TDP-43 species, since a caspase resistant mutant of TDP-43 is also converted into insoluble species. Phosphorylation at S409/410 apparently occurs late during the conversion of soluble to insoluble TDP-43, suggesting that phosphorylation is not a prerequisite for aggregation. Loss of function of the progranulin (PGRN) gene causes FTLD-U with TDP-43 positive inclusions and has been suggested to lead to caspase activation and subsequent TDP-43 processing. However, siRNA-mediated knockdown of PGRN in cell culture as well as a PGRN gene knockout in mice failed to cause the formation of the disease characterizing CTFs of TDP-43. Our findings therefore suggest that caspase-mediated processing generates CTFs of similar biochemical properties as those occurring in nuclear and cytoplasmic deposits of FTLD-U patients independent of PGRN levels.     

TaqMan Real-Time PCR Assays To Assess Arbuscular Mycorrhizal Responses to Field Manipulation of Grassland Biodiversity: Effects of Soil Characteristics,Plant Species Richness,and Functional Traits

Large-scale (temporal and/or spatial) molecular investigations of the diversity and distribution of arbuscular mycorrhizal fungi (AMF) require considerable sampling efforts and high-throughput analysis. To facilitate such efforts, we have developed a TaqMan real-time PCR assay to detect and identify AMF in environmental samples. First, we screened the diversity in clone libraries, generated by nested PCR, of the nuclear ribosomal DNA internal transcribed spacer (ITS) of AMF in environmental samples. We then generated probes and forward primers based on the detected sequences, enabling AMF sequence type-specific detection in TaqMan multiplex real-time PCR assays. In comparisons to conventional clone library screening and Sanger sequencing, the TaqMan assay approach provided similar accuracy but higher sensitivity with cost and time savings. The TaqMan assays were applied to analyze the AMF community composition within plots of a large-scale plant biodiversity manipulation experiment, the Jena Experiment, primarily designed to investigate the interactive effects of plant biodiversity on element cycling and trophic interactions. The results show that environmental variables hierarchically shape AMF communities and that the sequence type spectrum is strongly affected by previous land use and disturbance, which appears to favor disturbance-tolerant members of the genus Glomus. The AMF species richness of disturbance-associated communities can be largely explained by richness of plant species and plant functional groups, while plant productivity and soil parameters appear to have only weak effects on the AMF community.Arbuscular mycorrhizae are mutualistic associations between roots of plants and fungi that have been present for more than 400 million years (54). Approximately 80% of examined land plants (71), and almost all fungi of the phylum Glomeromycota (60), are capable of forming such associations. The main benefit of this relationship for plants is that it facilitates their acquisition of nutrients (especially P and N), while the fungus receives photoassimilates (7, 62). About 200 Glomeromycota species have been described to date, based on spore morphology (http://www.lrz-muenchen.de/∼schuessler/amphylo/amphylogeny.html), but there is increasing molecular evidence of significantly higher diversity in arbuscular mycorrhizal fungi (AMF) (10, 72).Diverse AMF communities have been detected in a wide range of plant communities (inter alia grasslands, boreal forests, and tropical communities; for an overview, see reference 48). Hence, AMF have been considered to be tolerant of wide ranges of ecological conditions and capable of associating with diverse plant partners. Identifying the factors regulating their community assemblages is challenging, but AMF community composition has been shown to be influenced by plant species diversity (e.g., see references 10, 22, and 33), and conversely, significant effects of AMF species and communities on the diversity and productivity of plant communities have been described (25, 68). Soil physicochemical parameters like phosphorus, nitrogen, and carbon availability (e.g., see references 4, 9, and 31); pH (17); moisture content (53); and disturbance (30) also reportedly influence AMF distribution. Hence, there is some support for niche theory, which presumes that two species of the same trophic level cannot coexist in a limited system and, if two species are present in such circumstances, one should become extinct (21). As a corollary, two cooccurring species must occupy niches that differ in some dimensions, e.g., plant hosts and/or soil properties (28). However, there are also indications that neutral ecological processes, as well as niche-defining parameters, may influence AMF diversity and community composition (17, 39). In contrast to niche theory, neutral theory (27) postulates that all individuals of every species at a given trophic level in a food web have ecological equivalence, and thus, all species within trophically defined communities can be regarded as open nonequilibrium assemblages that are solely shaped by dispersal and distinctions in spatiotemporal dimensions. According to the work of Hubbell (27), neutrality is defined at the level of individual organisms with identical probabilities of birth, death, migration, and speciation and not at the species level. In order to explore AMF communities more thoroughly and to test competing hypotheses, such as those raised by the niche and neutral theories, robust methods for high-throughput analyses of the communities are required.Recent investigations of variables that affect the structure of AMF communities have considered relationships between niche-defining dimensions, such as soil types (39) and pH gradients (17), and spatial variations in AMF community structure but not the role of plant diversity or functional traits of host plants. There have been several plant diversity manipulation experiments designed for coanalyzing multiple sets of ecological variables (e.g., the BIODEPTH and Cedar Creek projects) that would have been ideal for detailed examinations of effects of ecological variables on AMF, but previously reported AMF analyses in these experiments have been limited to counts of spores in a single study (11). However, not all AMF species regularly sporulate, and when present, spores poorly reflect AMF diversity (69), since active AMF occur as mycelia in roots and soils (e.g., see references 12 and 26). PCR-based molecular techniques enable much more rigorous characterization of AMF communities in these compartments (e.g., see references 26, 36, and 72), but assessments of broad spatial (42) and/or temporal (52) variations in AMF communities require analysis of large numbers of samples, which is not feasible using conventional PCR amplification followed by cloning and sequencing. This challenge can be potentially met by real-time PCR-based approaches, in which the AMF sequence types present in compartments of interest are first identified and then sequence type-specific probes are used for large-scale screening in real-time PCR TaqMan assays.In the study presented here, we explored AMF diversity in plots used in the Jena Experiment, a grassland plant diversity manipulation of 60 plant species representing four functional groups in 81 plots of 400 m² (56). The overall AMF diversity and community structure were first assessed by PCR amplification, cloning, and sequencing (55) of internal transcribed spacer (ITS) ribosomal DNA (rDNA) gene sequences in soil samples from 23 representative plots. Using the acquired data, we then developed sequence type-specific probes, which were applied in high-throughput real-time PCR TaqMan assays of samples from all 81 experimental plots, and the effects of 15 plant and soil variables on the AMF community assemblage were investigated.     

Disease‐linked TDP‐43 hyperphosphorylation suppresses TDP‐43 condensation and aggregation

Post‐translational modifications (PTMs) have emerged as key modulators of protein phase separation and have been linked to protein aggregation in neurodegenerative disorders. The major aggregating protein in amyotrophic lateral sclerosis and frontotemporal dementia, the RNA‐binding protein TAR DNA‐binding protein (TDP‐43), is hyperphosphorylated in disease on several C‐terminal serine residues, a process generally believed to promote TDP‐43 aggregation. Here, we however find that Casein kinase 1δ‐mediated TDP‐43 hyperphosphorylation or C‐terminal phosphomimetic mutations reduce TDP‐43 phase separation and aggregation, and instead render TDP‐43 condensates more liquid‐like and dynamic. Multi‐scale molecular dynamics simulations reveal reduced homotypic interactions of TDP‐43 low‐complexity domains through enhanced solvation of phosphomimetic residues. Cellular experiments show that phosphomimetic substitutions do not affect nuclear import or RNA regulatory functions of TDP‐43, but suppress accumulation of TDP‐43 in membrane‐less organelles and promote its solubility in neurons. We speculate that TDP‐43 hyperphosphorylation may be a protective cellular response to counteract TDP‐43 aggregation.     

Identification and characterization of 10 polymorphic microsatellite loci in the German cockroach,Blattella germanica

Primer sequences and initial characterization are presented for 10 microsatellite loci isolated from the German cockroach, Blattella germanica. In a sample of 30 individuals from a single population sample, all loci were polymorphic with two to 12 alleles segregating per locus and levels of observed heterozygosity ranging from 0.27 to 0.92. One locus showed a deficit of heterozygotes. Experimental conditions are described for polymerase chain reaction multiplexing, which enables the genotyping of eight loci in three electrophoretic runs consisting of one set of three and two sets of two markers. Seven primer sets cross‐amplify in the related Blattella asahinai.