Genetic dissection of heterosis using epistatic association mapping in a partial NCII mating design

Heterosis refers to the phenomenon in which an F1 hybrid exhibits enhanced growth or agronomic performance. However, previous theoretical studies on heterosis have been based on bi-parental segregating populations instead of F1 hybrids. To understand the genetic basis of heterosis, here we used a subset of F1 hybrids, named a partial North Carolina II design, to perform association mapping for dependent variables: original trait value, general combining ability (GCA), specific combining ability (SCA) and mid-parental heterosis (MPH). Our models jointly fitted all the additive, dominance and epistatic effects. The analyses resulted in several important findings: 1) Main components are additive and additive-by-additive effects for GCA and dominance-related effects for SCA and MPH, and additive-by-dominant effect for MPH was partly identified as additive effect; 2) the ranking of factors affecting heterosis was dominance > dominance-by-dominance > over-dominance > complete dominance; and 3) increasing the proportion of F1 hybrids in the population could significantly increase the power to detect dominance-related effects, and slightly reduce the power to detect additive and additive-by-additive effects. Analyses of cotton and rapeseed datasets showed that more additive-by-additive QTL were detected from GCA than from trait phenotype, and fewer QTL were from MPH than from other dependent variables.

Effect of surface restoring on subsurface variability in a climate model during 1949-2005

Initializing the ocean for decadal predictability studies is a challenge, as it requires reconstructing the little observed subsurface trajectory of ocean variability. In this study we explore to what extent surface nudging using well-observed sea surface temperature (SST) can reconstruct the deeper ocean variations for the 1949–2005 period. An ensemble made with a nudged version of the IPSLCM5A model and compared to ocean reanalyses and reconstructed datasets. The SST is restored to observations using a physically-based relaxation coefficient, in contrast to earlier studies, which use a much larger value. The assessment is restricted to the regions where the ocean reanalyses agree, i.e. in the upper 500 m of the ocean, although this can be latitude and basin dependent. Significant reconstruction of the subsurface is achieved in specific regions, namely region of subduction in the subtropical Atlantic, below the thermocline in the equatorial Pacific and, in some cases, in the North Atlantic deep convection regions. Beyond the mean correlations, ocean integrals are used to explore the time evolution of the correlation over 20-year windows. Classical fixed depth heat content diagnostics do not exhibit any significant reconstruction between the different existing observation-based references and can therefore not be used to assess global average time-varying correlations in the nudged simulations. Using the physically based average temperature above an isotherm (14 °C) alleviates this issue in the tropics and subtropics and shows significant reconstruction of these quantities in the nudged simulations for several decades. This skill is attributed to the wind stress reconstruction in the tropics, as already demonstrated in a perfect model study using the same model. Thus, we also show here the robustness of this result in an historical and observational context.

A comprehensive meta-analysis of common genetic variants in autism spectrum conditions

Background Autism spectrum conditions (ASC) are a group of neurodevelopmental conditions characterized by difficulties in social interaction and communication alongside repetitive and stereotyped behaviours. ASC are heritable, and common genetic variants contribute substantial phenotypic variability. More than 600 genes have been implicated in ASC to date. However, a comprehensive investigation of candidate gene association studies in ASC is lacking. Methods In this study, we systematically reviewed the literature for association studies for 552 genes associated with ASC. We identified 58 common genetic variants in 27 genes that have been investigated in three or more independent cohorts and conducted a meta-analysis for 55 of these variants. We investigated publication bias and sensitivity and performed stratified analyses for a subset of these variants. Results We identified 15 variants nominally significant for the mean effect size, 8 of which had P values below a threshold of significance of 0.01. Of these 15 variants, 11 were re-investigated for effect sizes and significance in the larger Psychiatric Genomics Consortium dataset, and none of them were significant. Effect direction for 8 of the 11 variants were concordant between both the datasets, although the correlation between the effect sizes from the two datasets was poor and non-significant. Conclusions This is the first study to comprehensively examine common variants in candidate genes for ASC through meta-analysis. While for majority of the variants, the total sample size was above 500 cases and 500 controls, the total sample size was not large enough to accurately identify common variants that contribute to the aetiology of ASC.

Functional identity and diversity of animals predict ecosystem functioning better than species-based indices

Drastic biodiversity declines have raised concerns about the deterioration of ecosystem functions and have motivated much recent research on the relationship between species diversity and ecosystem functioning. A functional trait framework has been proposed to improve the mechanistic understanding of this relationship, but this has rarely been tested for organisms other than plants. We analysed eight datasets, including five animal groups, to examine how well a trait-based approach, compared with a more traditional taxonomic approach, predicts seven ecosystem functions below- and above-ground. Trait-based indices consistently provided greater explanatory power than species richness or abundance. The frequency distributions of single or multiple traits in the community were the best predictors of ecosystem functioning. This implies that the ecosystem functions we investigated were underpinned by the combination of trait identities (i.e. single-trait indices) and trait complementarity (i.e. multi-trait indices) in the communities. Our study provides new insights into the general mechanisms that link biodiversity to ecosystem functioning in natural animal communities and suggests that the observed responses were due to the identity and dominance patterns of the trait composition rather than the number or abundance of species per se.

LOFZY 2005 - Field experiment on Cyclones over the Norwegian Sea: meteorological measurements of the research aircraft Falcon, 23 autonomous water buoys and radiosoundings at the research vessel Celtic Explorer

The field campaign LOFZY 2005 (LOFoten ZYklonen, engl.: Cyclones) was carried out in the frame of Collaborative Research Centre 512, which deals with low-pressure systems (cyclones) and the climate system of the North Atlantic. Cyclones are of special interest due to their influence on the interaction between atmosphere and ocean. Cyclone activity in the northern part of the Atlantic Ocean is notably high and is of particular importance for the entire Atlantic Ocean. An area of maximum precipitation exists in front of the Norwegian Lofoten islands. One aim of the LOFZY field campaign was to clarify the role cyclones play in the interaction of ocean and atmosphere. In order to obtain a comprehensive dataset of cyclone activity and ocean-atmosphere interaction a field experiment was carried out in the Lofoten region during March and April 2005. Employed platforms were the Irish research vessel RV Celtic Explorer which conducted a meteorological (radiosondes, standard parameters, observations) and an oceanographic (CTD) program. The German research aircraft Falcon accomplished eight flight missions (between 4-21 March) to observe synoptic conditions with high spatial and temporal resolution. In addition 23 autonomous marine buoys were deployed in advance of the campaign in the observed area to measure drift, air-temperature and -pressure and water-temperature. In addition to the published datasets several other measurements were performed during the experiment. Corresonding datasets will be published in the near future and are available on request. Details about all used platforms and sensors and all performed measurements are listed in the fieldreport. The following datasets are available on request: ground data at RV Celtic Explorer

abctools: an R package for tuning approximate Bayesian computation analyses

Approximate Bayesian computation (ABC) is a popular family of algorithms which perform approximate parameter inference when numerical evaluation of the likelihood function is not possible but data can be simulated from the model. They return a sample of parameter values which produce simulations close to the observed dataset. A standard approach is to reduce the simulated and observed datasets to vectors of summary statistics and accept when the difference between these is below a specified threshold. ABC can also be adapted to perform model choice. In this article, we present a new software package for R, abctools which provides methods for tuning ABC algorithms. This includes recent dimension reduction algorithms to tune the choice of summary statistics, and coverage methods to tune the choice of threshold. We provide several illustrations of these routines on applications taken from the ABC literature.

MOLES3: implementing an ISO standards driven data catalogue

ISO19156 Observations and Measurements (O&M) provides a standardised framework for organising information about the collection of information about the environment. Here we describe the implementation of a specialisation of O&M for environmental data, the Metadata Objects for Linking Environmental Sciences (MOLES3). MOLES3 provides support for organising information about data, and for user navigation around data holdings. The implementation described here, “CEDA-MOLES”, also supports data management functions for the Centre for Environmental Data Archival, CEDA. The previous iteration of MOLES (MOLES2) saw active use over five years, being replaced by CEDA-MOLES in late 2014. During that period important lessons were learnt both about the information needed, as well as how to design and maintain the necessary information systems. In this paper we review the problems encountered in MOLES2; how and why CEDA-MOLES was developed and engineered; the migration of information holdings from MOLES2 to CEDA-MOLES; and, finally, provide an early assessment of MOLES3 (as implemented in CEDA-MOLES) and its limitations. Key drivers for the MOLES3 development included the necessity for improved data provenance, for further structured information to support ISO19115 discovery metadata export (for EU INSPIRE compliance), and to provide appropriate fixed landing pages for Digital Object Identifiers (DOIs) in the presence of evolving datasets. Key lessons learned included the importance of minimising information structure in free text fields, and the necessity to support as much agility in the information infrastructure as possible without compromising on maintainability both by those using the systems internally and externally (e.g. citing in to the information infrastructure), and those responsible for the systems themselves. The migration itself needed to ensure continuity of service and traceability of archived assets.

Assessing species’ habitat associations from occurrence records, standardised monitoring data and expert opinion: a test with British butterflies

Accurate knowledge of species’ habitat associations is important for conservation planning and policy. Assessing habitat associations is a vital precursor to selecting appropriate indicator species for prioritising sites for conservation or assessing trends in habitat quality. However, much existing knowledge is based on qualitative expert opinion or local scale studies, and may not remain accurate across different spatial scales or geographic locations. Data from biological recording schemes have the potential to provide objective measures of habitat association, with the ability to account for spatial variation. We used data on 50 British butterfly species as a test case to investigate the correspondence of data-derived measures of habitat association with expert opinion, from two different butterfly recording schemes. One scheme collected large quantities of occurrence data (c. 3 million records) and the other, lower quantities of standardised monitoring data (c. 1400 sites). We used general linear mixed effects models to derive scores of association with broad-leaf woodland for both datasets and compared them with scores canvassed from experts. Scores derived from occurrence and abundance data both showed strongly positive correlations with expert opinion. However, only for occurrence data did these fell within the range of correlations between experts. Data-derived scores showed regional spatial variation in the strength of butterfly associations with broad-leaf woodland, with a significant latitudinal trend in 26% of species. Sub-sampling of the data suggested a mean sample size of 5000 occurrence records per species to gain an accurate estimation of habitat association, although habitat specialists are likely to be readily detected using several hundred records. Occurrence data from recording schemes can thus provide easily obtained, objective, quantitative measures of habitat association.

Agricultural innovations at a Late Iron Age oppidum: archaeobotanical evidence for flax, food and fodder from Calleva Atrebatum, UK

The development of oppida in the late first millennium BC across north-western Europe represents a major change in settlement form and social organisation. The construction of extensive earthwork systems, the presence of nucleated settlement areas, long-distance trade links and the development of hierarchical societies have been evidenced. These imply that changes in the style and organisation of agriculture would have been required to support these proto-urban population centres. Hypotheses of the subsistence bases of these settlements, ranging from a reliance on surplus arable production from local rural settlements, to an emphasis on pastoral activities, are here reviewed and grounded against a wider understanding of the expansion of agriculture in the Late Iron Age. These agricultural models have not been previously evaluated. This paper presents archaeobotanical data from six well fills from large-scale excavations at Late Iron Age and Early Roman Silchester, a Late Iron Age territorial oppidum and subsequent Roman civitas capital located in central-southern Britain. This is the first large-scale study of waterlogged plant macrofossils from within a settlement area of an oppidum. Waterlogged plant macrofossils were studied from a series of wells within the settlement. An assessment of taphonomy, considering stratigraphic and contextual information, is reported, followed by an analysis of the diverse assemblages of the plant remains through univariate analysis. Key results evidence animal stabling, flax cultivation, hay meadow management and the use of heathland resources. The staple crops cultivated and consumed at Late Iron Age and Early Roman Silchester are consistent with those cultivated in the wider region, whilst a range of imported fruits and flavourings were also present. The adoption of new oil crops and new grassland management shows that agricultural innovations were associated with foddering for animals rather than providing food for the proto-urban population. The evidence from Silchester is compared with other archaeobotanical datasets from oppida in Europe in order to identify key trends in agricultural change.

Tests of sunspot number sequences: 1. Using ionosonde data

More than 70 years ago it was recognised that ionospheric F2-layer critical frequencies [foF2] had a strong relationship to sunspot number. Using historic datasets from the Slough and Washington ionosondes, we evaluate the best statistical fits of foF2 to sunspot numbers (at each Universal Time [UT] separately) in order to search for drifts and abrupt changes in the fit residuals over Solar Cycles 17-21. This test is carried out for the original composite of the Wolf/Zürich/International sunspot number [R], the new “backbone” group sunspot number [RBB] and the proposed “corrected sunspot number” [RC]. Polynomial fits are made both with and without allowance for the white-light facular area, which has been reported as being associated with cycle-to-cycle changes in the sunspot number - foF2 relationship. Over the interval studied here, R, RBB, and RC largely differ in their allowance for the “Waldmeier discontinuity” around 1945 (the correction factor for which for R, RBB and RC is, respectively, zero, effectively over 20 %, and explicitly 11.6 %). It is shown that for Solar Cycles 18-21, all three sunspot data sequences perform well, but that the fit residuals are lowest and most uniform for RBB. We here use foF2 for those UTs for which R, RBB, and RC all give correlations exceeding 0.99 for intervals both before and after the Waldmeier discontinuity. The error introduced by the Waldmeier discontinuity causes R to underestimate the fitted values based on the foF2 data for 1932-1945 but RBB overestimates them by almost the same factor, implying that the correction for the Waldmeier discontinuity inherent in RBB is too large by a factor of two. Fit residuals are smallest and most uniform for RC and the ionospheric data support the optimum discontinuity multiplicative correction factor derived from the independent Royal Greenwich Observatory (RGO) sunspot group data for the same interval.

Positive phenotypic selection inferred from phylogenies

Rates of phenotypic evolution vary widely in nature and these rates may often reflect the intensity of natural selection. Here we outline an approach for detecting exceptional shifts in the rate of phenotypic evolution across phylogenies. We introduce a simple new branch-specific metric ∆V/∆B that divides observed phenotypic change along a branch into two components: (1) that attributable to the background rate (∆B), and (2) that attributable to departures from the background rate (∆V). Where the amount of expected change derived from variation in the rate of morphological evolution doubles that explained by to the background rate (∆V/∆B > 2), we identify this as positive phenotypic selection. We apply our approach to six datasets, finding multiple instances of positive selection in each. Our results support the growing appreciation that the traditional gradual view of phenotypic evolution is rarely upheld, with a more episodic view taking its place. This moves focus away from viewing phenotypic evolution as a simple homogeneous process and facilitates reconciliation with macroevolutionary interpretations from a genetic perspective, paving the way to novel insights into the link between genotype and phenotype. The ability to detect positive selection when genetic data are unavailable or unobtainable represents an attractive prospect for extant species, but when applied to fossil data it can reveal patterns of natural selection in deep time that would otherwise be impossible.

A cautionary note on the use of Ornstein Uhlenbeck models in macroevolutionary studies

Phylogenetic comparative methods are increasingly used to give new insights into the dynamics of trait evolution in deep time. For continuous traits the core of these methods is a suite of models that attempt to capture evolutionary patterns by extending the Brownian constant variance model. However, the properties of these models are often poorly understood, which can lead to the misinterpretation of results. Here we focus on one of these models – the Ornstein Uhlenbeck (OU) model. We show that the OU model is frequently incorrectly favoured over simpler models when using Likelihood ratio tests, and that many studies fitting this model use datasets that are small and prone to this problem. We also show that very small amounts of error in datasets can have profound effects on the inferences derived from OU models. Our results suggest that simulating fitted models and comparing with empirical results is critical when fitting OU and other extensions of the Brownian model. We conclude by making recommendations for best practice in fitting OU models in phylogenetic comparative analyses, and for interpreting the parameters of the OU model.

Preconditioning 2D integer data for fast convex hull computations

In order to accelerate computing the convex hull on a set of n points, a heuristic procedure is often applied to reduce the number of points to a set of s points, s ≤ n, which also contains the same hull. We present an algorithm to precondition 2D data with integer coordinates bounded by a box of size p × q before building a 2D convex hull, with three distinct advantages. First, we prove that under the condition min(p, q) ≤ n the algorithm executes in time within O(n); second, no explicit sorting of data is required; and third, the reduced set of s points forms a simple polygonal chain and thus can be directly pipelined into an O(n) time convex hull algorithm. This paper empirically evaluates and quantifies the speed up gained by preconditioning a set of points by a method based on the proposed algorithm before using common convex hull algorithms to build the final hull. A speedup factor of at least four is consistently found from experiments on various datasets when the condition min(p, q) ≤ n holds; the smaller the ratio min(p, q)/n is in the dataset, the greater the speedup factor achieved.

Abrupt transitions in the NAO control of explosive North Atlantic cyclone development

Explosive cyclones are intense extra-tropical low pressure systems featuring large deepening rates. In the Euro-Atlantic sector, they are a major source of life-threatening weather impacts due to their associated strong wind gusts, heavy precipitation and storm surges. The wintertime variability of the North Atlantic cyclonic activity is primarily modulated by the North Atlantic Oscillation (NAO). In this study, we investigate the interannual and multi-decadal variability of explosive North Atlantic cyclones using track density data from two reanalysis datasets (NCEP and ERA-40) and a control simulation of an atmosphere/ocean coupled General Circulation Model (GCM—ECHAM5/MPIOM1). The leading interannual and multi-decadal modes of variability of explosive cyclone track density are characterized by a strengthening/weakening pattern between Newfoundland and Iceland, which is mainly modulated by the NAO at both timescales. However, the NAO control of interannual cyclone variability is not stationary in time and abruptly fluctuates during periods of 20–25 years long both in NCEP and ECHAM5/MPIOM1. These transitions are accompanied by structural changes in the leading mode of explosive cyclone variability, and by decreased/enhanced baroclinicity over the sub-polar/sub-tropical North Atlantic. The influence of the ocean is apparently important for both the occurrence and persistence of such anomalous periods. In the GCM, the Atlantic Meridional Overturning Circulation appears to influence the large-scale baroclinicity and explosive cyclone development over the North Atlantic. These results permit a better understanding of explosive cyclogenesis variability at different climatic timescales and might help to improve predictions of these hazardous events.

Transcriptomics and proteomics show that selenium affects inflammation, ctoskeleton, and cancer pathways in human rectal biopsies

Epidemiologic studies highlight the potential role of dietary selenium (Se) in colorectal cancer prevention. Our goal was to elucidate whether expression of factors crucial for colorectal homoeostasis is affected by physiologic differences in Se status. Using transcriptomics and proteomics followed by pathway analysis, we identified pathways affected by Se status in rectal biopsies from 22 healthy adults, including 11 controls with optimal status (mean plasma Se = 1.43 μM) and 11 subjects with suboptimal status (mean plasma Se = 0.86 μM). We observed that 254 genes and 26 proteins implicated in cancer (80%), immune function and inflammatory response (40%), cell growth and proliferation (70%), cellular movement, and cell death (50%) were differentially expressed between the 2 groups. Expression of 69 genes, including selenoproteins W1 and K, which are genes involved in cytoskeleton remodelling and transcription factor NFκB signaling, correlated significantly with Se status. Integrating proteomics and transcriptomics datasets revealed reduced inflammatory and immune responses and cytoskeleton remodelling in the suboptimal Se status group. This is the first study combining omics technologies to describe the impact of differences in Se status on colorectal expression patterns, revealing that suboptimal Se status could alter inflammatory signaling and cytoskeleton in human rectal mucosa and so influence cancer risk.