Normalized Affymetrix expression data are biased by G-quadruplex formation

Probes with runs of four or more guanines (G-stacks) in their sequences can exhibit a level of hybridization that is unrelated to the expression levels of the mRNA that they are intended to measure. This is most likely caused by the formation of G-quadruplexes, where inter-probe guanines form Hoogsteen hydrogen bonds, which probes with G-stacks are capable of forming. We demonstrate that for a specific microarray data set using the Human HG-U133A Affymetrix GeneChip and RMA normalization there is significant bias in the expression levels, the fold change and the correlations between expression levels. These effects grow more pronounced as the number of G-stack probes in a probe set increases. Approximately 14 of the probe sets are directly affected. The analysis was repeated for a number of other normalization pipelines and two, FARMS and PLIER, minimized the bias to some extent. We estimate that ∼15 of the data sets deposited in the GEO database are susceptible to the effect. The inclusion of G-stack probes in the affected data sets can bias key parameters used in the selection and clustering of genes. The elimination of these probes from any analysis in such affected data sets outweighs the increase of noise in the signal. © 2011 The Author(s)

BioTIME: A database of biodiversity time series for the Anthropocene.

MotivationThe BioTIME database contains raw data on species identities and abundances in ecological assemblages through time. These data enable users to calculate temporal trends in biodiversity within and amongst assemblages using a broad range of metrics. BioTIME is being developed as a community-led open-source database of biodiversity time series. Our goal is to accelerate and facilitate quantitative analysis of temporal patterns of biodiversity in the Anthropocene.Main types of variables includedThe database contains 8,777,413 species abundance records, from assemblages consistently sampled for a minimum of 2 years, which need not necessarily be consecutive. In addition, the database contains metadata relating to sampling methodology and contextual information about each record.Spatial location and grainBioTIME is a global database of 547,161 unique sampling locations spanning the marine, freshwater and terrestrial realms. Grain size varies across datasets from 0.0000000158 km2 (158 cm2) to 100 km2 (1,000,000,000,000 cm2).Time period and grainBioTIME records span from 1874 to 2016. The minimal temporal grain across all datasets in BioTIME is a year.Major taxa and level of measurementBioTIME includes data from 44,440 species across the plant and animal kingdoms, ranging from plants, plankton and terrestrial invertebrates to small and large vertebrates.Software format.csv and .SQL

Response to Comment on “Plant diversity increases with the strength of negative density dependence at the global scale”

Hülsmann and Hartig suggest that ecological mechanisms other than specialized natural enemies or intraspecific competition contribute to our estimates of conspecific negative density dependence (CNDD). To address their concern, we show that our results are not the result of a methodological artifact and present a null-model analysis that demonstrates that our original findings—(i) stronger CNDD at tropical relative to temperate latitudes and (ii) a latitudinal shift in the relationship between CNDD and species abundance—persist even after controlling for other processes that might influence spatial relationships between adults and recruits

Local Factors Determine Plant Community Structure on Closely Neighbored Islands

Despite the recent popularity of the metacommunity concept, ecologists have not evaluated the applicability of different metacommunity frameworks to insular organisms. We surveyed 50 closely spaced islands in the Thousand-Island Lake of China to examine the role of local (environmental) and regional (dispersal) factors in structuring woody plant assemblages (tree and shrub species) on these islands. By partitioning the variation in plant community structure into local and regional causes, we showed that local environmental conditions, specifically island morphometric characteristics, accounted for the majority of the variation in plant community structure among the studied islands. Spatial variables, representing the potential importance of species dispersal, explained little variation. We conclude that one metacommunity framework–species sorting–best characterizes these plant communities. This result reinforces the idea that the traditional approach of emphasizing the local perspective when studying ecological communities continues to hold its value

Patterns and drivers of tree Mortality in Iberian Forests: climatic effects are modified by competition

Tree mortality is a key process underlying forest dynamics and community assembly. Understanding how tree mortality is driven by simultaneous drivers is needed to evaluate potential effects of climate change on forest composition. Using repeat-measure information fromc.400,000 trees from the Spanish Forest Inventory, we quantified the relative importance of tree size, competition, climate and edaphic conditions on tree mortality of 11 species, and explored the combined effect of climate and competition. Tree mortality was affected by all of these multiple drivers, especially tree size and asymmetric competition, and strong interactions between climate and competition were found. All species showed L-shaped mortality patterns (i.e. showed decreasing mortality with tree size), but pines were more sensitive to asymmetric competition than broadleaved species. Among climatic variables, the negative effect of temperature on tree mortality was much larger than the effect of precipitation. Moreover, the effect of climate (mean annual temperature and annual precipitation) on tree mortality was aggravated at high competition levels for all species, but especially for broadleaved species. The significant interaction between climate and competition on tree mortality indicated that global change in Mediterranean regions, causing hotter and drier conditions and denser stands, could lead to profound effects on forest structure and composition. Therefore, to evaluate the potential effects of climatic change on tree mortality, forest structure must be considered, since two systems of similar composition but different structure could radically differ in their response to climatic conditions

A niche remedy for the dynamical problems of neutral theory

We demonstrate how niche theory and Hubbell's original formulation of neutral theory can be blended together into a general framework modeling the combined effects of selection, drift, speciation, and dispersal on community dynamics. This framework connects many seemingly unrelated ecological population models, and allows for quantitative predictions to be made about the impact of niche stabilizing and destabilizing forces on population extinction times and abundance distributions. In particular, the existence of niche stabilizing forces in our blended framework can simultaneously resolve two major problems with the dynamics of neutral theory, namely predictions of species lifetimes that are too short and species ages that are too long.Comment: 47 pages, 4 figures, Accepted to Theoretical Ecolog

Bedform migration in a mixed sand and cohesive clay intertidal environment and implications for bed material transport predictions

Many coastal and estuarine environments are dominated by mixtures of non-cohesive sand and cohesive mud. The migration rate of bedforms, such as ripples and dunes, in these environments is important in determining bed material transport rates to inform and assess numerical models of sediment transport and geomorphology. However, these models tend to ignore parameters describing the physical and biological cohesion (resulting from clay and extracellular polymeric substances, EPS) in natural mixed sediment, largely because of a scarcity of relevant laboratory and field data. To address this gap in knowledge, data were collected on intertidal flats over a spring-neap cycle to determine the bed material transport rates of bedforms in biologically-active mixed sand-mud. Bed cohesive composition changed from below 2 vol% up to 5.4 vol% cohesive clay, as the tide progressed from spring towards neap. The amount of EPS in the bed sediment was found to vary linearly with the clay content. Using multiple linear regression, the transport rate was found to depend on the Shields stress parameter and the bed cohesive clay content. The transport rates decreased with increasing cohesive clay and EPS content, when these contents were below 2.8 vol% and 0.05 wt%, respectively. Above these limits, bedform migration and bed material transport was not detectable by the instruments in the study area. These limits are consistent with recently conducted sand-clay and sand-EPS laboratory experiments on bedform development. This work has important implications for the circumstances under which existing sand-only bedform migration transport formulae may be applied in a mixed sand-clay environment, particularly as 2.8 vol% cohesive clay is well within the commonly adopted definition of “clean sand”

Direct and indirect effects of climate on richness drive the latitudinal diversity gradient in forest trees

Data accessibility statement: Full census data are available upon reasonable request from the ForestGEO data portal, http://ctfs.si.edu/datarequest/ We thank Margie Mayfield, three anonymous reviewers and Jacob Weiner for constructive comments on the manuscript. This study was financially supported by the National Key R&D Program of China (2017YFC0506100), the National Natural Science Foundation of China (31622014 and 31570426), and the Fundamental Research Funds for the Central Universities (17lgzd24) to CC. XW was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB3103). DS was supported by the Czech Science Foundation (grant no. 16-26369S). Yves Rosseel provided us valuable suggestions on using the lavaan package conducting SEM analyses. Funding and citation information for each forest plot is available in the Supplementary Information Text 1.Peer reviewedPostprin

Aboveground biomass density models for NASA's Global Ecosystem Dynamics Investigation (GEDI) lidar mission

NASA's Global Ecosystem Dynamics Investigation (GEDI) is collecting spaceborne full waveform lidar data with a primary science goal of producing accurate estimates of forest aboveground biomass density (AGBD). This paper presents the development of the models used to create GEDI's footprint-level (similar to 25 m) AGBD (GEDI04_A) product, including a description of the datasets used and the procedure for final model selection. The data used to fit our models are from a compilation of globally distributed spatially and temporally coincident field and airborne lidar datasets, whereby we simulated GEDI-like waveforms from airborne lidar to build a calibration database. We used this database to expand the geographic extent of past waveform lidar studies, and divided the globe into four broad strata by Plant Functional Type (PFT) and six geographic regions. GEDI's waveform-to-biomass models take the form of parametric Ordinary Least Squares (OLS) models with simulated Relative Height (RH) metrics as predictor variables. From an exhaustive set of candidate models, we selected the best input predictor variables, and data transformations for each geographic stratum in the GEDI domain to produce a set of comprehensive predictive footprint-level models. We found that model selection frequently favored combinations of RH metrics at the 98th, 90th, 50th, and 10th height above ground-level percentiles (RH98, RH90, RH50, and RH10, respectively), but that inclusion of lower RH metrics (e.g. RH10) did not markedly improve model performance. Second, forced inclusion of RH98 in all models was important and did not degrade model performance, and the best performing models were parsimonious, typically having only 1-3 predictors. Third, stratification by geographic domain (PFT, geographic region) improved model performance in comparison to global models without stratification. Fourth, for the vast majority of strata, the best performing models were fit using square root transformation of field AGBD and/or height metrics. There was considerable variability in model performance across geographic strata, and areas with sparse training data and/or high AGBD values had the poorest performance. These models are used to produce global predictions of AGBD, but will be improved in the future as more and better training data become available