ALBERT

Description: Connections between vegetation and soil thermal dynamics are critical for estimating the vulnerability of permafrost to thaw with continued climate warming and vegetation changes. The interplay of complex biophysical processes results in a highly heterogeneous soil temperature distribution on small spatial scales. Moreover, the link between topsoil temperature and active layer thickness remains poorly constrained. Sixty-eight temperature loggers were installed at 1–3 cm depth to record the distribution of topsoil temperatures at the Trail Valley Creek study site in the northwestern Canadian Arctic. The measurements were distributed across six different vegetation types characteristic for this landscape. Two years of topsoil temperature data were analysed statistically to identify temporal and spatial characteristics and their relationship to vegetation, snow cover, and active layer thickness. The mean annual topsoil temperature varied between −3.7 and 0.1 ∘C within 0.5 km2. The observed variation can, to a large degree, be explained by variation in snow cover. Differences in snow depth are strongly related with vegetation type and show complex associations with late-summer thaw depth. While cold winter soil temperature is associated with deep active layers in the following summer for lichen and dwarf shrub tundra, we observed the opposite beneath tall shrubs and tussocks. In contrast to winter observations, summer topsoil temperature is similar below all vegetation types with an average summer topsoil temperature difference of less than 1 ∘C. Moreover, there is no significant relationship between summer soil temperature or cumulative positive degree days and active layer thickness. Altogether, our results demonstrate the high spatial variability of topsoil temperature and active layer thickness even within specific vegetation types. Given that vegetation type defines the direction of the relationship between topsoil temperature and active layer thickness in winter and summer, estimates of permafrost vulnerability based on remote sensing or model results will need to incorporate complex local feedback mechanisms of vegetation change and permafrost thaw.

Description: In this paper, we quantify the CO2 and N2O emissions from denitrification over the Amazonian wetlands. The study concerns the entire Amazonian wetland ecosystem with a specific focus on three floodplain (FP) locations: the Branco FP, the Madeira FP and the FP alongside the Amazon River. We adapted a simple denitrification model to the case of tropical wetlands and forced it by open water surface extent products from the Soil Moisture and Ocean Salinity (SMOS) satellite. A priori model parameters were provided by in situ observations and gauging stations from the HYBAM Observatory. Our results show that the denitrification and the trace gas emissions present a strong cyclic pattern linked to the inundation processes that can be divided into three distinct phases: activation, stabilization and deactivation. We quantify the average yearly denitrification and associated emissions of CO2 and N2O over the entire watershed at 17.8 kgN ha−1 yr−1, 0.37 gC-CO2 m−2 yr−1 and 0.18 gN-N2O m−2 yr−1 respectively for the period 2011–2015. When compared to local observations, it was found that the CO2 emissions accounted for 0.01 % of the integrated ecosystem, which emphasizes the fact that minor changes to the land cover may induce strong impacts on the Amazonian carbon budget. Our results are consistent with the state of the art of global nitrogen models with a positive bias of 28 %. When compared to other wetlands in different pedoclimatic environments we found that the Amazonian wetlands have similar emissions of N2O with the Congo tropical wetlands and lower emissions than the temperate and tropical anthropogenic wetlands of the Garonne (France), the Rhine (Europe) and south-eastern Asia rice paddies. In summary our paper shows that a data-model-based approach can be successfully applied to quantify N2O and CO2 fluxes associated with denitrification over the Amazon basin. In the future, the use of higher-resolution remote sensing products from sensor fusion or new sensors like the Surface Water and Ocean Topography (SWOT) mission will permit the transposition of the approach to other large-scale watersheds in tropical environments.

Description: Anaerobic nitrate-dependent Fe(II) oxidation (NDFeO) is widespread in various aquatic environments and plays a major role in iron and nitrogen redox dynamics. However, evidence for truly enzymatic, autotrophic NDFeO remains limited, with alternative explanations involving the coupling of heterotrophic denitrification with the abiotic oxidation of structurally bound or aqueous Fe(II) by reactive intermediate nitrogen (N) species (chemodenitrification). The extent to which chemodenitrification is caused (or enhanced) by ex vivo surface catalytic effects has not been directly tested to date. To determine whether the presence of either an Fe(II)-bearing mineral or dead biomass (DB) catalyses chemodenitrification, two different sets of anoxic batch experiments were conducted: 2 mM Fe(II) was added to a low-phosphate medium, resulting in the precipitation of vivianite (Fe3(PO4)2), to which 2 mM nitrite (NO2-) was later added, with or without an autoclaved cell suspension (∼1.96×108 cells mL−1) of Shewanella oneidensis MR-1. Concentrations of nitrite (NO2-), nitrous oxide (N2O), and iron (Fe2+, Fetot) were monitored over time in both set-ups to assess the impact of Fe(II) minerals and/or DB as catalysts of chemodenitrification. In addition, the natural-abundance isotope ratios of NO2- and N2O (δ15N and δ18O) were analysed to constrain the associated isotope effects. Up to 90 % of the Fe(II) was oxidized in the presence of DB, whereas only ∼65 % of the Fe(II) was oxidized under mineral-only conditions, suggesting an overall lower reactivity of the mineral-only set-up. Similarly, the average NO2- reduction rate in the mineral-only experiments (0.004±0.003 mmol L−1 d−1) was much lower than in the experiments with both mineral and DB (0.053±0.013 mmol L−1 d−1), as was N2O production (204.02±60.29 nmol L−1 d−1). The N2O yield per mole NO2- reduced was higher in the mineral-only set-ups (4 %) than in the experiments with DB (1 %), suggesting the catalysis-dependent differential formation of NO. N-NO2- isotope ratio measurements indicated a clear difference between both experimental conditions: in contrast to the marked 15N isotope enrichment during active NO2- reduction (15εNO2=+10.3 ‰) observed in the presence of DB, NO2- loss in the mineral-only experiments exhibited only a small N isotope effect (

Description: Temperature appears to be the best predictor of species composition of planktonic foraminifera communities, making it possible to use their fossil assemblages to reconstruct sea surface temperature (SST) variation in the past. However, the role of other environmental factors potentially modulating the spatial and vertical distribution of planktonic foraminifera species is poorly understood. This is especially relevant for environmental factors affecting the subsurface habitat. If such factors play a role, changes in the abundance of subsurface-dwelling species may not solely reflect SST variation. In order to constrain the effect of subsurface parameters on species composition, we here characterize the vertical distribution of living planktonic foraminifera community across an east–west transect through the subtropical South Atlantic Ocean, where SST variability was small, but the subsurface water mass structure changed dramatically. Four planktonic foraminifera communities could be identified across the top 700 m of the transect. Gyre and Agulhas Leakage surface faunas were predominantly composed of Globigerinoides ruber, Globigerinoides tenellus, Trilobatus sacculifer, Globoturborotalita rubescens, Globigerinella calida, Tenuitella iota, and Globigerinita glutinata, and these only differed in terms of relative abundances (community composition). Upwelling fauna was dominated by Neogloboquadrina pachyderma, Neogloboquadrina incompta, Globorotalia crassaformis, and Globorotalia inflata. Thermocline fauna was dominated by Tenuitella fleisheri, Globorotalia truncatulinoides, and Globorotalia scitula in the west and by G. scitula only in the east. The largest part of the standing stock was consistently found in the surface layer, but SST was not the main predictor of species composition either for the depth-integrated fauna across the stations or at individual depth layers. Instead, we identified a pattern of vertical stacking of different parameters controlling species composition, reflecting different aspects of the pelagic habitat. Whereas productivity appears to dominate in the mixed layer (0–60 m), physical properties (temperature, salinity) become important at intermediate depths and in the subsurface, a complex combination of factors including oxygen concentration is required to explain the assemblage composition. These results indicate that the seemingly straightforward relationship between assemblage composition and SST in sedimentary assemblages reflects vertically and seasonally integrated processes that are only indirectly linked to SST. It also implies that fossil assemblages of planktonic foraminifera should also contain a signature of subsurface processes, which could be used for paleoceanographic reconstructions.

Description: Light-use efficiency defines the ability of primary producers to convert sunlight energy to primary production and is computed as the ratio between the gross primary production and the intercepted photosynthetic active radiation. While this measure has been applied broadly within terrestrial ecology to investigate habitat resource-use efficiency, it remains underused within the aquatic realm. This report provides a conceptual framework to compute hourly and daily light-use efficiency using underwater O2 eddy covariance, a recent technological development that produces habitat-scale rates of primary production under unaltered in situ conditions. The analysis, tested on two benthic flux datasets, documents that hourly light-use efficiency may approach the theoretical limit of 0.125 O2 per photon under low-light conditions, but it decreases rapidly towards the middle of the day and is typically 10-fold lower on a 24 h basis. Overall, light-use efficiency provides a useful measure of habitat functioning and facilitates site comparison in time and space.

Description: Fishermen are known to try to avoid fishing in stormy weather, as storms pose a physical threat to fishers, their vessels, and their gear. In this article, a dataset and methods are developed to investigate the degree to which fishers avoid storms, estimate storm aversion parameters, and explore how this response varies across vessel characteristics and across regions of the United States. The data consist of vessel-level trip-taking decisions from six federal fisheries across the United States combined with marine storm warning data from the National Weather Service. The estimates of storm aversion can be used to parameterize predictive models. Fishers’ aversion to storms decreases with increasing vessel size and increases with the severity of the storm warning. This information contributes to our understanding of the risk-to-revenue trade-off that fishers evaluate every time they consider going to sea, and of the propensity of fishers to take adaptive actions to avoid facing additional physical risk.

Description: Sardine Sardinops sagax is an ecologically and economically important Clupeid found off the entire South African coast that includes both coastal upwelling and western boundary current systems. Although the management of the sardine fisheries historically assumed a single, panmictic population, the existence of three, semi-discrete subpopulations has recently been hypothesized. We conducted otolith δ18O and microstructure analyses to investigate nursery habitat temperatures and early life growth rates, respectively, of sardine collected from three biogeographic regions around South Africa’s coast to test that hypothesis. Analyses indicated that for both summer- and winter-captured adults and summer-captured juveniles, fishes from the west coast grew significantly slower in water that was several degrees cooler than those from the south and east coasts. This suggests that mixing of sardines between regions, particularly the west and other coasts, is relatively limited and supports the hypothesis of semi-discrete subpopulations. However, the west-south differences disappeared in the results for winter-captured juveniles, suggesting that differences in early life conditions between regions may change seasonally, and/or that all or most winter-captured juveniles originated from the west coast. Further elucidating the interactions between South African sardine subpopulations and the mechanisms thereof is important for sustainable harvesting of this species.

Description: The CloudRoots field experiment was designed to obtain a comprehensive observational dataset that includes soil, plant, and atmospheric variables to investigate the interaction between a heterogeneous land surface and its overlying atmospheric boundary layer at the sub-hourly and sub-kilometre scale. Our findings demonstrate the need to include measurements at leaf level to better understand the relations between stomatal aperture and evapotranspiration (ET) during the growing season at the diurnal scale. Based on these observations, we obtain accurate parameters for the mechanistic representation of photosynthesis and stomatal aperture. Once the new parameters are implemented, the model reproduces the stomatal leaf conductance and the leaf-level photosynthesis satisfactorily. At the canopy scale, we find a consistent diurnal pattern on the contributions of plant transpiration and soil evaporation using different measurement techniques. From highly resolved vertical profile measurements of carbon dioxide (CO2) and other state variables, we infer a profile of the CO2 assimilation in the canopy with non-linear variations with height. Observations taken with a laser scintillometer allow us to quantify the non-steadiness of the surface turbulent fluxes during the rapid changes driven by perturbation of photosynthetically active radiation by cloud flecks. More specifically, we find 2 min delays between the cloud radiation perturbation and ET. To study the relevance of advection and surface heterogeneity for the land–atmosphere interaction, we employ a coupled surface–atmospheric conceptual model that integrates the surface and upper-air observations made at different scales from leaf to the landscape. At the landscape scale, we calculate a composite sensible heat flux by weighting measured fluxes with two different land use categories, which is consistent with the diurnal evolution of the boundary layer depth. Using sun-induced fluorescence measurements, we also quantify the spatial variability of ET and find large variations at the sub-kilometre scale around the CloudRoots site. Our study shows that throughout the entire growing season, the wide variations in stomatal opening and photosynthesis lead to large diurnal variations of plant transpiration at the leaf, plant, canopy, and landscape scales. Integrating different advanced instrumental techniques with modelling also enables us to determine variations of ET that depend on the scale where the measurement were taken and on the plant growing stage.

Description: In the western Mediterranean Sea, Levantine intermediate waters (LIW), which circulate below the surface productive zone, progressively accumulate nutrients along their pathway from the Tyrrhenian Sea to the Algerian Basin. This study addresses the role played by diffusion in the nutrient enrichment of the LIW, a process particularly relevant inside step-layer structures extending down to deep waters – structures known as thermohaline staircases. Profiling float observations confirmed that staircases develop over epicentral regions confined in large-scale circulation features and maintained by saltier LIW inflows on the periphery. Thanks to a high profiling frequency over the 4-year period 2013–2017, float observations reveal the temporal continuity of the layering patterns encountered during the cruise PEACETIME and document the evolution of layer properties by about +0.06 ∘C in temperature and +0.02 in salinity. In the Algerian Basin, the analysis of in situ lateral density ratios untangled double-diffusive convection as a driver of thermohaline changes inside epicentral regions and isopycnal diffusion as a driver of heat and salt exchanges with the surrounding sources. In the Tyrrhenian Sea, the nitrate flux across thermohaline staircases, as opposed to the downward salt flux, contributes up to 25 % of the total nitrate pool supplied to the LIW by vertical transfer. Overall, however, the nutrient enrichment of the LIW is driven mostly by other sources, coastal or atmospheric, as well as by inputs advected from the Algerian Basin.

Description: Peatlands store ∼ 20 %–30 % of the global soil organic carbon stock and are an important source of dissolved organic carbon (DOC) for inland waters. Recent improvements for in situ optical monitoring revealed that the DOC concentration in streams draining peatlands is highly variable, showing seasonal variation and short and intense DOC concentration peaks. This study aimed to statistically determine the variables driving stream DOC concentration variations at seasonal and event scales. Two mountainous peatlands (one fen and one bog) were monitored in the French Pyrenees to capture their outlet DOC concentration variability at a high-frequency rate (30 min). Abiotic variables including precipitation, stream temperature and water level, water table depth, and peat water temperature were also monitored at high frequency and used as potential predictors to explain DOC concentration variability. Results show that at both sites DOC concentration time series can be decomposed into a seasonal baseline interrupted by many short and intense peaks of higher concentrations. The DOC concentration baseline is driven, at the seasonal scale, by peat water temperature. At the event scale, DOC concentration increases are mostly driven by a rise in the water table within the peat at both sites. Univariate linear models between DOC concentration and peat water temperature or water table increases show greater efficiency at the fen site. Water recession times were derived from water level time series using master recession curve coefficients. They vary greatly between the two sites but also within one peatland site. They partly explain the differences between DOC dynamics in the studied peatlands, including peat porewater DOC concentrations and the links between stream DOC concentration and water table rise within the peatlands. This highlights that peatland complexes are composed of a mosaic of heterogeneous peat units distinctively producing or transferring DOC to streams.

Description: CerealsDB (www.cerealsdb.uk.net) is an online repository of mainly hexaploid wheat (Triticum aestivum) single nucleotide polymorphisms (SNPs) and genotyping data. The CerealsDB website has been designed to enable wheat breeders and scientists to select the appropriate markers for research breeding tasks, such as marker-assisted selection. We report a large update of genotyping information for over 6000 wheat accessions and describe new webtools for exploring and visualizing the data. We also describe a new database of quantitative trait loci that links phenotypic traits to CerealsDB SNP markers and allelic scores for each of those markers. CerealsDB is an open-access website that hosts information on wheat SNPs considered useful for both plant breeders and research scientists. The latest CerealsDB database is available at https://www.cerealsdb.uk.net/cerealgenomics/CerealsDB/indexNEW.php.

Description: Anoxic sediments, as compared to oxic settings, encompass a much higher proportion of relatively labile and thus more reactive organic matter, naturally giving rise to structural changes of the organic molecules themselves, as well as cross-linking between them (e.g., through reactive sulfur species). Both processes transform the original biomolecules into geomolecules. For the oxic environment, these intermolecular and intramolecular transformations also operate, but cross-linking may be less important since the labile, reactive component is rapidly removed. As such, one may expect a structurally better preservation of the more refractory initial biomolecules in the oxic environment. To test this hypothesis, initially identical biomolecules need to be compared between different preservational environments. Here, we use the species-specific morphology of organic microfossils to assure a single initial biosynthetic product (the cysts of the fossil dinoflagellate species Thalassiphora pelagica) for comparison. We assess the macromolecular structures of cysts from the Eocene (∼40 Ma) sulfidic Rhine Graben and the oxic Kerguelen Plateau and compare them with each other and the structures of recent cysts. While between the sites the T. pelagica cysts are morphologically identical and show no signs of morphological modification, pyrolysis gas chromatography mass spectroscopy and micro Fourier transform infrared analyses show that their macromolecular characteristics are strongly different. Comparison with recent cysts shows that the cysts deposited in the sulfidic Rhine Graben show a strong additional contribution of long-chain aliphatic moieties and thus less diagenetic intermolecular cross-linking. The presence of organic sulfur identifies natural volcanization as one of the diagenetic processes. Furthermore, we observe a loss of bound oxygen and no trace of the original carbohydrate signature of the cyst wall biomacromolecule. The material deposited in the oxic sediments of the Kerguelen Plateau shows no traces of sulfurization. It shows a minor contribution of short carbon chains only and thus less diagenetic intermolecular cross-linking. Furthermore, a carbohydrate signature was still preserved evidencing a better molecular preservation of the initial biomacromolecule, supporting our initial hypothesis. This shows that excellent morphological preservation does not imply excellent chemical preservation. It also leads to the conclusion that the best preservation of molecular structure is not necessarily where most organic matter gets preserved, which, in turn, is important for understanding the nature and fate of sedimentary organic matter and its isotopic signature.

Description: Progress towards ecosystem-based fisheries management calls for useful tools to prioritize actions. To select suitable methods for local circumstances, evaluating approaches used in other jurisdictions can be a cost-effective first step. We tested Productivity Susceptibility Analysis (PSA) to assess the potential vulnerability of the marine fish community in the Skagerrak–Kattegat (Eastern North Sea) to possible interactions with all Swedish fisheries operating in the area. This analysis combines attributes for a species productivity with attributes related to the susceptibility to capture to quantify a single score for vulnerability: high, medium, or low risk. Results indicate that demersal trawl and gillnet fisheries were associated with the highest risk levels if interaction occurs, i.e. having the highest prevalence of species with potentially high vulnerability to the fisheries. Mixed results were seen when comparing the assessment results with available data. The main benefit of utilizing PSA in the area is the comprehensiveness of the assessment, including data-deficient fisheries and species. Drawbacks include potential overestimation of actual risks. Overall, together with available data, PSA in the studied area provides a comprehensive map of potential risks for further actions and may progress a science-based, precautionary management of the area.

Description: While fisheries provide food and employment for hundreds of millions of people, they also can have significant impact on biodiversity. We explore the potential of area-based fisheries management to simultaneously maintain biodiversity and high levels of sustainable food production. We used two illustrative examples of fisheries that have different gear types, areas, and species to evaluate the trade-off between biodiversity and harvest. We calculate the optimal effort by gear and area that maximizes a weighted objective function of biodiversity and harvest, ranging from 100% of the weight on harvest to 100% on biodiversity. We found for both case studies that the trade-off was highly convex, with win–win solutions allowing for high levels of both fishery harvest and conservation. This is achieved by reducing or eliminating fishing effort that negatively impacts high conservation value species while maintaining fishing effort with gears and in areas where there is low conservation impact. We suggest that, in most fisheries, such situations can be found and that effective area-based management can provide for high levels of biodiversity protection and food production.

Description: The most prolific acacias in southern Europe (Acacia dealbata, A. melanoxylon and A. longifolia) are rapidly spreading in its westernmost area: Portugal and NW Spain, where congeners with invasion potential are already established. We performed a bibliographic search of acacia invasions in southern Europe and used spatial data on acacia distribution and abiotic parameters in Iberia to check the influence of abiotic factors on acacia invasion. According to our results, in Iberia A. dealbata and A. melanoxylon seem limited by high soil pH (pHCaCl2〉5.5), frequent frosts (〉21 to 40 d yr−1) and low annual precipitation (

Description: Boron isotope systematics of planktonic foraminifera from core-top sediments and culture experiments have been studied to investigate the sensitivity of δ11B of calcite tests to seawater pH. However, our knowledge of the relationship between δ11B and pH remains incomplete for many taxa. Thus, to expand the potential scope of application of this proxy, we report δ11B data for seven different species of planktonic foraminifera from sediment core tops. We utilize a method for the measurement of small samples of foraminifera and calculate the δ11B-calcite sensitivity to pH for Globigerinoides ruber, Trilobus sacculifer (sacc or without sacc), Orbulina universa, Pulleniatina obliquiloculata, Neogloboquadrina dutertrei, Globorotalia menardii, and Globorotalia tumida, including for unstudied core tops and species. These taxa have diverse ecological preferences and are from sites that span a range of oceanographic regimes, including some that are in regions of air–sea equilibrium and others that are out of equilibrium with the atmosphere. The sensitivity of δ11Bcarbonate to δ11Bborate (e.g., Δδ11Bcarbonate∕Δδ11Bborate) in core tops is consistent with previous studies for T. sacculifer and G. ruber and close to unity for N. dutertrei, O. universa, and combined deep-dwelling species. Deep-dwelling species closely follow the core-top calibration for O. universa, which is attributed to respiration-driven microenvironments likely caused by light limitation and/or symbiont–host interactions. Our data support the premise that utilizing boron isotope measurements of multiple species within a sediment core can be utilized to constrain vertical profiles of pH and pCO2 at sites spanning different oceanic regimes, thereby constraining changes in vertical pH gradients and yielding insights into the past behavior of the oceanic carbon pumps.

Description: Carbon–nitrogen (C–N) interactions regulate N availability for plant growth and for emissions of nitrous oxide (N2O) and the uptake of carbon dioxide. Future projections of these terrestrial greenhouse gas fluxes are strikingly divergent, leading to major uncertainties in projected global warming. Here we analyse the large increase in terrestrial N2O emissions over the past 21 000 years as reconstructed from ice-core isotopic data and presented in part 1 of this study. Remarkably, the increase occurred in two steps, each realized over decades and within a maximum of 2 centuries, at the onsets of the major deglacial Northern Hemisphere warming events. The data suggest a highly dynamic and responsive global N cycle. The increase may be explained by an increase in the flux of reactive N entering and leaving ecosystems or by an increase in N2O yield per unit N converted. We applied the LPX-Bern dynamic global vegetation model in deglacial simulations forced with Earth system model climate data to investigate N2O emission patterns, mechanisms, and C–N coupling. The N2O emission changes are mainly attributed to changes in temperature and precipitation and the loss of land due to sea-level rise. LPX-Bern simulates a deglacial increase in N2O emissions but underestimates the reconstructed increase by 47 %. Assuming time-independent N sources in the model to mimic progressive N limitation of plant growth results in a decrease in N2O emissions in contrast to the reconstruction. Our results appear consistent with suggestions of (a) biological controls on ecosystem N acquisition and (b) flexibility in the coupling of the C and N cycles during periods of rapid environmental change. A dominant uncertainty in the explanation of the reconstructed N2O emissions is the poorly known N2O yield per N lost through gaseous pathways and its sensitivity to soil conditions. The deglacial N2O record provides a constraint for future studies.

Description: Seamounts represent ideal systems to study the influence and interdependency of environmental gradients at a single geographic location. These topographic features represent a prominent habitat for various forms of life, including microbiota and macrobiota, spanning benthic as well as pelagic organisms. While it is known that seamounts are globally abundant structures, it still remains unclear how and to which extent the complexity of the sea floor is intertwined with the local oceanographic mosaic, biogeochemistry, and microbiology of a seamount ecosystem. Along these lines, the present study aimed to explore whether and to what extent seamounts can have an imprint on the microbial community composition of seawater and of sessile benthic invertebrates, sponges. For our high-resolution sampling approach of microbial diversity (16S rRNA gene amplicon sequencing) along with measurements of inorganic nutrients and other biogeochemical parameters, we focused on the Schulz Bank seamount ecosystem, a sponge ground ecosystem which is located on the Arctic Mid-Ocean Ridge. Seawater samples were collected at two sampling depths (mid-water, MW, and near-bed water, BW) from a total of 19 sampling sites. With a clustering approach we defined microbial microhabitats within the pelagic realm at Schulz Bank, which were mapped onto the seamount's topography and related to various environmental parameters (such as suspended particulate matter, SPM; dissolved inorganic carbon, DIC; silicate, SiO4-; phosphate, PO43-; ammonia, NH4+; nitrate, NO32-; nitrite, NO2-; depth; and dissolved oxygen, O2). The results of our study reveal a “seamount effect” (sensu stricto) on the microbial mid-water pelagic community at least 200 m above the sea floor. Further, we observed a strong spatial heterogeneity in the pelagic microbial landscape across the seamount, with planktonic microbial communities reflecting oscillatory and circulatory water movements, as well as processes of bentho-pelagic coupling. Depth, NO32-, SiO4-, and O2 concentrations differed significantly between the determined pelagic microbial clusters close to the sea floor (BW), suggesting that these parameters were presumably linked to changes in microbial community structures. Secondly, we assessed the associated microbial community compositions of three sponge species along a depth gradient of the seamount. While sponge-associated microbial communities were found to be mainly species-specific, we also detected significant intra-specific differences between individuals, depending on the pelagic near-bed cluster they originated from. The variable microbial phyla (i.e. phyla which showed significant differences across varying depth, NO32-, SiO4-, O2 concentrations, and different from local seawater communities) were distinct for every sponge species when considering average abundances per species. Variable microbial phyla included representatives of both those taxa traditionally counted for the variable community fraction and taxa counted traditionally for the core community fraction. Microbial co-occurrence patterns for the three examined sponge species Geodia hentscheli, Lissodendoryx complicata, and Schaudinnia rosea were distinct from each other. Over all, this study shows that topographic structures such as the Schulz Bank seamount can have an imprint (seamount effect sensu lato) on both the microbial community composition of seawater and sessile benthic invertebrates such as sponges by an interplay between the geology, physical oceanography, biogeochemistry, and microbiology of seamounts.

Description: Repeatability of study setups and reproducibility of research results by underlying data are major requirements in science. Until now, abstract models for describing the structural logic of studies in environmental sciences are lacking and tools for data management are insufficient. Mandatory for repeatability and reproducibility is the use of sophisticated data management solutions going beyond data file sharing. Particularly, it implies maintenance of coherent data along workflows. Design data concern elements from elementary domains of operations being transformation, measurement and transaction. Operation design elements and method information are specified for each consecutive workflow segment from field to laboratory campaigns. The strict linkage of operation design element values, operation values and objects is essential. For enabling coherence of corresponding objects along consecutive workflow segments, the assignment of unique identifiers and the specification of their relations are mandatory. The abstract model presented here addresses these aspects, and the software DiversityDescriptions (DWB-DD) facilitates the management of thusly connected digital data objects and structures. DWB-DD allows for an individual specification of operation design elements and their linking to objects. Two workflow design use cases, one for DNA barcoding and another for cultivation of fungal isolates, are given. To publish those structured data, standard schema mapping and XML-provision of digital objects are essential. Schemas useful for this mapping include the Ecological Markup Language, the Schema for Meta-omics Data of Collection Objects and the Standard for Structured Descriptive Data. Data pipelines with DWB-DD include the mapping and conversion between schemas and functions for data publishing and archiving according to the Open Archival Information System standard. The setting allows for repeatability of study setups, reproducibility of study results and for supporting work groups to structure and maintain their data from the beginning of a study. The theory of ‘FAIR++’ digital objects is introduced.

Description: The annual carbon and water dynamics of two eastern North American temperate forests were compared over a 6-year period from 2012 to 2017. The geographic location, forest age, soil, and climate were similar between the two stands; however, stand composition varied in terms of tree leaf-retention and shape strategy: one stand was a deciduous broadleaf forest, while the other was an evergreen needleleaf forest. The 6-year mean annual net ecosystem productivity (NEP) of the coniferous forest was slightly higher and more variable (218±109 g C m−2 yr−1) compared to that of the deciduous forest NEP (200±83 g C m−2 yr−1). Similarly, the 6-year mean annual evapotranspiration (ET) of the coniferous forest was higher (442±33 mm yr−1) than that of the deciduous forest (388±34 mm yr−1), but with similar interannual variability. Summer meteorology greatly impacted the carbon and water fluxes in both stands; however, the degree of response varied among the two stands. In general, warm temperatures caused higher ecosystem respiration (RE), resulting in reduced annual NEP values – an impact that was more pronounced at the deciduous broadleaf forest compared to the evergreen needleleaf forest. However, during warm and dry years, the evergreen forest had largely reduced annual NEP values compared to the deciduous forest. Variability in annual ET at both forests was related most to the variability in annual air temperature (Ta), with the largest annual ET observed in the warmest years in the deciduous forest. Additionally, ET was sensitive to prolonged dry periods that reduced ET at both stands, although the reduction at the coniferous forest was relatively larger than that of the deciduous forest. If prolonged periods (weeks to months) of increased Ta and reduced precipitation are to be expected under future climates during summer months in the study region, our findings suggest that the deciduous broadleaf forest will likely remain an annual carbon sink, while the carbon sink–source status of the coniferous forest remains uncertain.

Description: Aside from many well-known sources, the greenhouse gas methane (CH4) was recently discovered entrapped in the sediments of Swiss Alpine glacier forefields derived from calcareous bedrock. A first study performed in one glacial catchment indicated that CH4 was ubiquitous in sediments and rocks and was largely of thermogenic origin. Here, we present the results of a follow-up study that aimed at (1) determining the occurrence and origin of sediment-entrapped CH4 in other calcareous glacier forefields across Switzerland and (2) providing an inventory of this sediment-entrapped CH4, i.e., determining the contents and total mass of CH4 present, and its spatial distribution within and between five different Swiss glacier forefields situated on calcareous formations of the Helvetic nappes in the Central Alps. Sediment and bedrock samples were collected at high spatial resolution from the forefields of Im Griess, Griessfirn, Griessen, Wildstrubel, and Tsanfleuron glaciers, representing different geographic and geologic regions of the Helvetic nappes. We performed geochemical analyses on gas extracted from sediments and rocks, including the determination of CH4 contents, stable carbon-isotope analyses (δ13CCH4), and the determination of gas-wetness ratios (ratio of CH4 to ethane and propane contents). To estimate the total mass of CH4 entrapped in glacier-forefield sediments, the total volume of sediment was determined based on the measured forefield area and either literature values of mean sediment thickness or direct depth measurements using electrical resistivity tomography. Methane was found in all sediments (0.08–73.81 µg CH4 g−1 dry weight) and most rocks (0.06–108.58 µg CH4 g−1) collected from the five glacier forefields, confirming that entrapped CH4 is ubiquitous in these calcareous formations. Geochemical analyses further confirmed a thermogenic origin of the entrapped CH4 (average δ13CCH4 of sediment of −28.23 (± 3.42) ‰; average gas-wetness ratio of 75.2 (± 48.4)). Whereas sediment-entrapped CH4 contents varied moderately within individual forefields, we noted a large, significant difference in the CH4 content and total CH4 mass (range of 200–3881 t CH4) between glacier forefields at the regional scale. The lithology and tectonic setting within the Helvetic nappes appeared to be dominant factors determining rock and sediment CH4 contents. Overall, a substantial quantity of CH4 was found to be entrapped in Swiss calcareous glacier forefields. Its potential release and subsequent fate in this environment is the subject of ongoing studies.

Description: The dynamics of the physicochemical and biological parameters were followed during the decline of a Cymodocea nodosa meadow in the northern Adriatic Sea from July 2017 to October 2018. During the regular growth of C. nodosa from July 2017 to March 2018, the species successfully adapted to the changes in environmental conditions and prevented H2S accumulation by its reoxidation, supplying the sediment with O2 from the water column and/or leaf photosynthesis. The C. nodosa decline was most likely triggered in April 2018 when light availability to the plant was drastically reduced due to increased seawater turbidity that resulted from increased terrigenous input, indicated by a decrease in salinity accompanied with a substantial increase in particulate matter concentration, combined with resuspension of sediment and elevated autotrophic biomass. Light reduction impaired photosynthesis of C. nodosa and the oxidation capability of belowground tissue. Simultaneously, a depletion of oxygen due to intense oxidation of H2S occurred in the sediment, thus creating anoxic conditions in most of the rooted areas. These linked negative effects on the plant performance caused an accumulation of H2S in the sediments of the C. nodosa meadow. During the decay of aboveground and belowground tissues, culminating in August 2018, high concentrations of H2S were reached and accumulated in the sediment as well as in bottom waters. The influx of oxygenated waters in September 2018 led to the re-establishment of H2S oxidation in the sediment and remainder of the belowground tissue. Our results indicate that if disturbances of environmental conditions, particularly those compromising the light availability, take place during the recruitment phase of plant growth when metabolic needs are at a maximum and stored reserves minimal, a sudden and drastic decline of the seagrass meadow occurs.

Description: Recent application of Fourier transform near infra-red spectroscopy (FT-NIRS) to predict age in fish otoliths has gained attention among fisheries managers as a potential alternative to costly production ageing of managed species. We assessed the age prediction capability of FT-NIRS scans in whole otoliths from red snapper, Lutjanus campechanus, collected from the US Gulf of Mexico and US Atlantic Ocean (South Atlantic). Otoliths were scanned with an FT-NIR spectrometer and resulting spectral signatures were regressed with traditionally estimated ages via partial least squares regression to produce calibration models, which were validated for predictive capability against test sets of otoliths. Calibration models successfully predicted age with R2 ranging 0.94–0.95, mean squared error ≤1.8 years, and bias 31 years were not well predicted, possibly due to light attenuation in the thickest otoliths. Our results suggest that FT-NIRS can improve efficiency in production ageing for fisheries management while maintaining data quality standards.

Description: This article investigates the diet of the snow crab (Chionoecetes opilio) and its feeding intensity in the Barents Sea. Data show that snow crab has a diverse diet that includes almost all types of benthic invertebrates living in the Barents Sea. There are differences between the diets of females and males and of juveniles and adults. Juveniles and females typically occupy shallow areas with communities of bivalve molluscs, while males typically live deeper on slopes and depressions where polychaetes and crustaceans are the most abundant groups. Stomach contents were analysed to determine the species composition and frequency of occurrence of various benthic taxa. Consumption of food was estimated and compared with data from the Russian seas of the Pacific region. The total annual consumption of macrozoobenthos by snow crab was calculated in accordance with its current distribution in the Barents Sea. Snow crab consumes at least 30 000 tonnes of benthos annually, which amounts to 0.1–0.2% of the total macrozoobenthic biomass in the investigated area. The population of snow crab causes the largest impact on the benthic communities in the northeastern part of the Barents Sea and near the south side of the Novaya Zemlya archipelago.

Description: Until the late 19th century, extensive beds of flat oyster Ostrea edulis populated the Central North Sea, which have vanished after intensive fisheries. At present, various initiatives are being carried out to investigate the potential to restore this former key species in the area. This historical ecological study contributes by delineating the former oyster bed area and through an assessment of its limits against known gradients in the North Sea. Extensive data from historical maps, texts, and ship-based surveys were used to synthesize our knowledge on the former beds. It was revealed that the area with oyster beds covered ∼6.2% of the total North Sea bottom, with a delineation that could partly be explained by hydrodynamic and temperature gradients. The position and extent of the area are notably different from the area that is used in recent feasibility studies on the restoration of North Sea oyster beds. The offshore oysters lived on muddy sand in relatively cold conditions, and there are several indications that their reproductive rate was low. The apparent disappearance of cold water adapted flat oysters will challenge restoration projects. This study provides indispensable information for the future restoration of flat oyster beds in the North Sea.

Description: Many wetlands have been drained due to urbanization, agriculture, forestry or other purposes, which has resulted in a loss of their ecosystem services. To protect receiving waters and to achieve services such as flood control and storm water quality mitigation, new wetlands are created in urbanized areas. However, our knowledge of greenhouse gas exchange in newly created wetlands in urban areas is currently limited. In this paper we present measurements carried out at a created urban wetland in Southern Finland in the boreal climate. We conducted measurements of ecosystem CO2 flux and CH4 flux (FCH4) at the created storm water wetland Gateway in Nummela, Vihti, Southern Finland, using the eddy covariance (EC) technique. The measurements were commenced the fourth year after construction and lasted for 1 full year and two subsequent growing seasons. Besides ecosystem-scale fluxes measured by the EC tower, the diffusive CO2 and CH4 fluxes from the open-water areas (FwCO2 and FwCH4, respectively) were modelled based on measurements of CO2 and CH4 concentration in the water. Fluxes from the vegetated areas were estimated by applying a simple mixing model using the above-mentioned fluxes and the footprint-weighted fractional area. The half-hourly footprint-weighted contribution of diffusive fluxes from open water ranged from 0 % to 25.5 % in 2013. The annual net ecosystem exchange (NEE) of the studied wetland was 8.0 g C-CO2 m−2 yr−1, with the 95 % confidence interval between −18.9 and 34.9 g C-CO2 m−2 yr−1, and FCH4 was 3.9 g C-CH4 m−2 yr−1, with the 95 % confidence interval between 3.75 and 4.07 g C-CH4 m−2 yr−1. The ecosystem sequestered CO2 during summer months (June–August), while the rest of the year it was a CO2 source. CH4 displayed strong seasonal dynamics, higher in summer and lower in winter, with a sporadic emission episode in the end of May 2013. Both CH4 and CO2 fluxes, especially those obtained from vegetated areas, exhibited strong diurnal cycles during summer with synchronized peaks around noon. The annual FwCO2 was 297.5 g C-CO2 m−2 yr−1 and FwCH4 was 1.73 g C-CH4 m−2 yr−1. The peak diffusive CH4 flux was 137.6 nmol C-CH4 m−2 s−1, which was synchronized with the FCH4. Overall, during the monitored time period, the established storm water wetland had a climate-warming effect with 0.263 kg CO2-eq m−2 yr−1 of which 89 % was contributed by CH4. The radiative forcing of the open-water areas exceeded that of the vegetation areas (1.194 and 0.111 kg CO2-eq m−2 yr−1, respectively), which implies that, when considering solely the climate impact of a created wetland over a 100-year horizon, it would be more beneficial to design and establish wetlands with large patches of emergent vegetation and to limit the areas of open water to the minimum necessitated by other desired ecosystem services.

Description: Looking beyond stratification: a model-based analysis of the biological drivers of oxygen depletion in the North Sea Biogeosciences Discussions, 12, 12543-12610, 2015 Author(s): F. Große, N. Greenwood, M. Kreus, H. J. Lenhart, D. Machoczek, J. Pätsch, L. A. Salt, and H. Thomas The problem of low oxygen conditions, often referred to as hypoxia, occurs regularly in the North Sea, a temperate European shelf sea. Stratification represents a major process regulating the seasonal dynamics of bottom oxygen. However, lowest oxygen conditions in the North Sea do not occur in the regions of strongest stratification. This suggests that stratification is an important prerequisite for hypoxia, but that the complex interaction between hydrodynamics and the biological processes drives its development. In this study we use the ecosystem model HAMSOM-ECOHAM5 to provide a general characteristic of the different North Sea oxygen regimes, and to quantify the impact of the different physical and biological factors driving the oxygen dynamics below the thermocline and in the bottom layer. We show that the North Sea can be subdivided into three different regimes in terms of oxygen dynamics: (1) a highly productive, non-stratified coastal regime, (2) a productive, seasonally stratified regime with a small sub-thermocline volume, and (3) a productive, seasonally stratified regime with a large sub-thermocline volume, with regime 2 being highly susceptible to hypoxic conditions. Our analysis of the different processes driving the oxygen development reveals that inter-annual variations in the oxygen conditions are caused by variations in primary production, while spatial differences can be attributed to differences in stratification and water depth. In addition, we show that benthic bacteria represent the main oxygen consumers in the bottom layer, consistently accounting for more than 50 % of the overall consumption. By providing these valuable insights, we show that ecosystem models can be a useful tool for the interpretation of observations and the estimation of the impact of anthropogenic drivers on the North Sea oxygen conditions.

Description: : As sequencing becomes cheaper and more widely available, there is a greater need to quickly and effectively analyze large-scale genomic data. While the functionality of AVIA v1.0, whose implementation was based on ANNOVAR, was comparable with other annotation web servers, AVIA v2.0 represents an enhanced web-based server that extends genomic annotations to cell-specific transcripts and protein-level functional annotations. With AVIA’s improved interface, users can better visualize their data, perform comprehensive searches and categorize both coding and non-coding variants. Availability and implementation : AVIA is freely available through the web at http://avia.abcc.ncifcrf.gov . Contact : Hue.Vuong@fnlcr.nih.gov Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : As new methods for multivariate analysis of genome wide association studies become available, it is important to be able to combine results from different cohorts in a meta-analysis. The R package MultiMeta provides an implementation of the inverse-variance-based method for meta-analysis, generalized to an n -dimensional setting. Availability and implementation: The R package MultiMeta can be downloaded from CRAN. Contact: dragana.vuckovic@burlo.trieste.it ; vi1@sanger.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Optimizing the impact of temperature on bio-hydrogen production from food waste and its derivatives under no pH control using statistical modelling Biogeosciences Discussions, 12, 12823-12850, 2015 Author(s): A. Sattar, C. Arslan, C. Ji, S. Sattar, K. Yousaf, and S. Hashim The effect of temperature on bio-hydrogen production by co-digestion of sewerage sludge with food waste and its two derivatives, i.e. noodle waste and rice waste, was investigated by statistical modelling. Experimental results showed that increasing temperature from mesophilic (37 °C) to thermophilic (55 °C) was an effective mean for increasing bio-hydrogen production from food waste and noodle waste, but it caused a negative impact on bio-hydrogen production from rice waste. The maximum cumulative bio-hydrogen production of 650 mL was obtained from noodle waste under mesophilic temperature condition. Most of the production was observed during 48 h of incubation that continued till 72 h of incubation, and a decline in pH during this interval was 4.3 and 4.4 from a starting value of 7 under mesophilic and thermophilic conditions, respectively. Most of glucose consumption was also observed during 72 h of incubation and the maximum consumption was observed during the first 24 h, which was the same duration where the maximum pH drop occurred. The maximum hydrogen yields of 82.47 mL VS −1 , 131.38 mL COD −1 , and 44.90 mL glucose −1 were obtained from mesophilic food waste, thermophilic noodle waste and mesophilic rice waste respectively. The production of volatile fatty acids increased with an increase in time and temperature from food waste and noodle waste reactors whereas it decreased with temperature in rice waste reactors. The statistical modelling returned good results with high values of coefficient of determination ( R 2 ) for each waste type when it was opted for the study of cumulative hydrogen production, glucose consumption and volatile fatty acid production. The 3-D response surface plots developed by the statistical models helped a lot in developing better understanding of the impact of temperature and incubation time.

Description: Some of the most dangerous pathogens such as Mycobacterium tuberculosis and Yersinia pestis evolve clonally . This means that little or no recombination occurs between strains belonging to these species. Paradoxically, although different members of these species show extreme sequence similarity of orthologous genes, some show considerable intraspecies phenotypic variation, the source of which remains elusive. To examine the possible sources of phenotypic variation within clonal pathogenic bacterial species, we carried out an extensive genomic and pan-genomic analysis of the sources of genetic variation available to a large collection of clonal and nonclonal pathogenic bacterial species. We show that while nonclonal species diversify through a combination of changes to gene sequences, gene loss and gene gain, gene loss completely dominates as a source of genetic variation within clonal species. Indeed, gene loss is so prevalent within clonal species as to lead to levels of gene content variation comparable to those found in some nonclonal species that are much more diverged in their gene sequences and that acquire a substantial number of genes horizontally. Gene loss therefore needs to be taken into account as a potential dominant source of phenotypic variation within clonal bacterial species.

Description: Obligate bacterial symbionts are widespread in many invertebrates, where they are often confined to specialized host cells and are transmitted directly from mother to progeny. Increasing numbers of these bacteria are being characterized but questions remain about their population structure and evolution. Here we take a comparative genomics approach to investigate two prominent bacterial symbionts (BFo1 and BFo2) isolated from geographically separated populations of western flower thrips, Frankliniella occidentalis. Our multifaceted approach to classifying these symbionts includes concatenated multilocus sequence analysis (MLSA) phylogenies, ribosomal multilocus sequence typing (rMLST), construction of whole-genome phylogenies, and in-depth genomic comparisons. We showed that the BFo1 genome clusters more closely to species in the genus Erwinia, and is a putative close relative to Erwinia aphidicola . BFo1 is also likely to have shared a common ancestor with Erwinia pyrifoliae/Erwinia amylovora and the nonpathogenic Erwinia tasmaniensis and genetic traits similar to Erwinia billingiae . The BFo1 genome contained virulence factors found in the genus Erwinia but represented a divergent lineage. In contrast, we showed that BFo2 belongs within the Enterobacteriales but does not group closely with any currently known bacterial species. Concatenated MLSA phylogenies indicate that it may have shared a common ancestor to the Erwinia and Pantoea genera, and based on the clustering of rMLST genes, it was most closely related to Pantoea ananatis but represented a divergent lineage. We reconstructed a core genome of a putative common ancestor of Erwinia and Pantoea and compared this with the genomes of BFo bacteria. BFo2 possessed none of the virulence determinants that were omnipresent in the Erwinia and Pantoea genera. Taken together, these data are consistent with BFo2 representing a highly novel species that maybe related to known Pantoea .

Description: Stable isotope study of a new chondrichthyan fauna (Kimmeridgian, Porrentruy, Swiss Jura): an unusual freshwater-influenced isotopic composition for the hybodont shark Asteracanthus Biogeosciences Discussions, 12, 12899-12921, 2015 Author(s): L. Leuzinger, L. Kocsis, J.-P. Billon-Bruyat, S. Spezzaferri, and T. Vennemann Chondrichthyan teeth (sharks, rays and chimaeras) are mineralised in isotopic equilibrium with the surrounding water, and parameters such as water temperature and salinity can be inferred from the oxygen isotopic composition (δ 18 O p ) of their bioapatite. We analysed a new chondrichthyan assemblage, as well as teeth from bony fish (Pycnodontiformes). All specimens are from Kimmeridgian coastal marine deposits of the Swiss Jura (vicinity of Porrentruy, Ajoie district, NW Switzerland). While the overall faunal composition and the isotopic composition of bony fish are consistent with marine conditions, unusually low δ 18 O p values were measured for the hybodont shark Asteracanthus . These values are also lower compared to previously published data from older European Jurassic localities. Additional analyses on material from Solothurn (Kimmeridgian, NW Switzerland) also have comparable, low- 18 O isotopic compositions for Asteracanthus . The data are hence interpreted to represent a so far unique, freshwater-influenced isotopic composition for this shark that is classically considered as a marine genus. While reproduction in freshwater or brackish realms is established for other hybodonts, a similar behaviour for Asteracanthus is proposed here. Regular excursions into lower salinity waters can be linked to the age of the deposits and correspond to an ecological adaptation, most likely driven by the Kimmeridgian transgression and by the competition of the primitive shark Asteracanthus with the rapidly diversifying neoselachians (modern sharks).

Description: Climate change impacts on net primary production (NPP) and export production (EP) regulated by increasing stratification and phytoplankton community structure in CMIP5 models Biogeosciences Discussions, 12, 12851-12897, 2015 Author(s): W. Fu, J. Randerson, and J. K. Moore We examine climate change impacts on net primary production (NPP) and export production (sinking particulate flux; EP) with simulations from nine Earth System Models (ESMs) performed in the framework of the fifth Coupled Model Inter-comparison Project (CMIP5). Global NPP and EP are reduced considerably by the end of the century for the intense warming scenario of Representative Concentration Pathway (RCP) 8.5. Relative to the 1990s, global NPP in the 2090s is reduced by 2.3–16 % and EP by 7–18 %. The models with the largest increases in stratification (and largest relative reductions in NPP and EP) also show the largest positive biases in stratification for the contemporary period, suggesting some potential overestimation of climate impacts on NPP and EP. All of the CMIP5 models show an increase in stratification in response to surface ocean warming and freshening that is accompanied by decreases in NPP, EP, and surface macronutrient concentrations. There is considerable variability across models in the absolute magnitude of these fluxes, surface nutrient concentrations, and their perturbations by climate change, indicating large model uncertainties. The negative response of NPP and EP to stratification increases reflects a bottom-up control, as nutrient flux to the euphotic zone declines. Models with dynamic phytoplankton community structure show larger declines in EP than in NPP. This is driven by phytoplankton community composition shifts, with a reduced percentage of NPP by large phytoplankton under RCP 8.5, as smaller phytoplankton are favored under the increasing nutrient stress. Thus, projections of the NPP response to climate change in the CMIP5 models are critically dependent on the simulated phytoplankton community structure, the efficiency of the biological pump, and the resulting (highly variable) levels of regenerated production. Community composition is represented relatively simply in the CMIP5 models, and should be expanded to better capture the spatial patterns and the changes in export efficiency that are necessary for predicting climate impacts on NPP.

Description: Technical Note: A simple calculation algorithm to separate high-resolution CH 4 flux measurements into ebullition and diffusion-derived components Biogeosciences Discussions, 12, 12923-12945, 2015 Author(s): M. Hoffmann, M. Schulz-Hanke, J. Garcia Alba, N. Jurisch, U. Hagemann, T. Sachs, M. Sommer, and J. Augustin Processes driving the production, transformation and transport of methane (CH 4 ) in wetland ecosystems are highly complex. Thus, serious challenges are constitutes in terms of the mechanistic process understanding, the identification of potential environmental drivers and the calculation of reliable CH 4 emission estimates. We present a simple calculation algorithm to separate open-water CH 4 fluxes measured with automatic chambers into diffusion- and ebullition-derived components, which helps facilitating the identification of underlying dynamics and potential environmental drivers. Flux separation is based on ebullition related sudden concentration changes during single measurements. A variable ebullition filter is applied, using the lower and upper quartile and the interquartile range (IQR). Automation of data processing is achieved by using an established R-script, adjusted for the purpose of CH 4 flux calculation. The algorithm was tested using flux measurement data (July to September 2013) from a former fen grassland site, converted into a shallow lake as a result of rewetting ebullition and diffusion contributed 46 and 55 %, respectively, to total CH 4 emissions, which is comparable to those previously reported by literature. Moreover, the separation algorithm revealed a concealed shift in the diurnal trend of diffusive fluxes throughout the measurement period.

Description: Gene expression evolution occurs through changes in cis - or trans -regulatory elements or both. Interactions between transcription factors (TFs) and their binding sites (TFBSs) constitute one of the most important points where these two regulatory components intersect. In this study, we investigated the evolution of TFBSs in the promoter regions of different Saccharomyces strains and species. We divided the promoter of a gene into the proximal region and the distal region, which are defined, respectively, as the 200-bp region upstream of the transcription starting site and as the 200-bp region upstream of the proximal region. We found that the predicted TFBSs in the proximal promoter regions tend to be evolutionarily more conserved than those in the distal promoter regions. Additionally, Saccharomyces cerevisiae strains used in the fermentation of alcoholic drinks have experienced more TFBS losses than gains compared with strains from other environments (wild strains, laboratory strains, and clinical strains). We also showed that differences in TFBSs correlate with the cis component of gene expression evolution between species (comparing S. cerevisiae and its sister species Saccharomyces paradoxus ) and within species (comparing two closely related S. cerevisiae strains).

Description: Gene duplication is a key factor contributing to phenotype diversity across and within species. Although the availability of complete genomes has led to the extensive study of genomic duplications, the dynamics and variability of gene duplications mediated by retrotransposition are not well understood. Here, we predict mRNA retrotransposition and use comparative genomics to investigate their origin and variability across primates. Analyzing seven anthropoid primate genomes, we found a similar number of mRNA retrotranspositions (~7,500 retrocopies) in Catarrhini (Old Word Monkeys, including humans), but a surprising large number of retrocopies (~10,000) in Platyrrhini (New World Monkeys), which may be a by-product of higher long interspersed nuclear element 1 activity in these genomes. By inferring retrocopy orthology, we dated most of the primate retrocopy origins, and estimated a decrease in the fixation rate in recent primate history, implying a smaller number of species-specific retrocopies. Moreover, using RNA-Seq data, we identified approximately 3,600 expressed retrocopies. As expected, most of these retrocopies are located near or within known genes, present tissue-specific and even species-specific expression patterns, and no expression correlation to their parental genes. Taken together, our results provide further evidence that mRNA retrotransposition is an active mechanism in primate evolution and suggest that retrocopies may not only introduce great genetic variability between lineages but also create a large reservoir of potentially functional new genomic loci in primate genomes.

Description: Nitrogen export from a boreal stream network following forest harvesting: seasonal nitrate removal and conservative export of organic forms Biogeosciences Discussions, 12, 12061-12089, 2015 Author(s): J. Schelker, R. Sponseller, E. Ring, L. Högbom, S. Löfgren, and H. Laudon Boreal streams are under pressure from large scale disturbance by forestry. Recent scenarios predict an increase in forest production in Scandinavia to meet market demands and to mitigate higher anthropogenic CO 2 emissions. Increased fertilization and shorter forest rotations are anticipated which will likely enhance the pressure on boreal streams in the near future. Among the major environmental impacts of forest harvesting is the increased mobilization of inorganic nitrogen (N), primarily as nitrate (NO 3 - ) into surface waters. But whereas NO 3 - inputs to first-order streams have been previously described, their downstream fate and impact is not well understood. We evaluated the downstream fate of N inputs in a boreal landscape that has been altered by forest harvests over a 10 year period to estimate the effects of multiple clear-cuts on aquatic N export in a boreal stream network. Small streams showed substantial leaching of NO 3 - in response to harvests with concentrations increasing by ~ 15 fold. NO 3 - concentrations at two sampling stations further downstream in the network were strongly seasonal and increased significantly in response to harvesting at the medium size, but not at the larger stream. Nitrate removal efficiency, E r , calculated as the percentage of "forestry derived" NO 3 - that was retained within the landscape using a mass balance model was highest during the snow melt season followed by the growing season, but declined continuously throughout the dormant season. In contrast, export of organic N from the landscape indicated little removal and was essentially conservative. Overall, net removal of NO 3 - between 2008 and 2011 accounted for ~ 70 % of the total NO 3 - mass exported from harvested patches distributed across the landscape. These results highlight the capacity and limitation of N-limited terrestrial and aquatic ecosystems to buffer inorganic N mobilization that arises from multiple clear-cuts within meso-scale boreal watersheds.

Description: Viruses rely completely on the hosts’ machinery for translation of viral transcripts. However, for most viruses infecting humans, codon usage preferences (CUPrefs) do not match those of the host. Human papillomaviruses (HPVs) are a showcase to tackle this paradox: they present a large genotypic diversity and a broad range of phenotypic presentations, from asymptomatic infections to productive lesions and cancer. By applying phylogenetic inference and dimensionality reduction methods, we demonstrate first that genes in HPVs are poorly adapted to the average human CUPrefs, the only exception being capsid genes in viruses causing productive lesions. Phylogenetic relationships between HPVs explained only a small proportion of CUPrefs variation. Instead, the most important explanatory factor for viral CUPrefs was infection phenotype, as orthologous genes in viruses with similar clinical presentation displayed similar CUPrefs. Moreover, viral genes with similar spatiotemporal expression patterns also showed similar CUPrefs. Our results suggest that CUPrefs in HPVs reflect either variations in the mutation bias or differential selection pressures depending on the clinical presentation and expression timing. We propose that poor viral CUPrefs may be central to a trade-off between strong viral gene expression and the potential for eliciting protective immune response.

Description: Huntington's disease (HD) is a hereditary neurodegenerative disorder characterized by brain atrophy particularly in striatum leading to personality changes, chorea and dementia. Glycogen synthase kinase-3 (GSK-3) is a serine/threonine kinase in the crossroad of many signaling pathways that is highly pleiotropic as it phosphorylates more than hundred substrates including structural, metabolic, and signaling proteins. Increased GSK-3 activity is believed to contribute to the pathogenesis of neurodegenerative diseases like Alzheimer's disease and GSK-3 inhibitors have been postulated as therapeutic agents for neurodegeneration. Regarding HD, GSK-3 inhibitors have shown beneficial effects in cell and invertebrate animal models but no evident efficacy in mouse models. Intriguingly, those studies were performed without interrogating GSK-3 level and activity in HD brain. Here we aim to explore the level and also the enzymatic activity of GSK-3 in the striatum and other less affected brain regions of HD patients and of the R6/1 mouse model to then elucidate the possible contribution of its alteration to HD pathogenesis by genetic manipulation in mice. We report a dramatic decrease in GSK-3 levels and activity in striatum and cortex of HD patients with similar results in the mouse model. Correction of the GSK-3 deficit in HD mice, by combining with transgenic mice with conditional GSK-3 expression, resulted in amelioration of their brain atrophy and behavioral motor and learning deficits. Thus, our results demonstrate that decreased brain GSK-3 contributes to HD neurological phenotype and open new therapeutic opportunities based on increasing GSK-3 activity or attenuating the harmful consequences of its decrease.

Description: Spinocerebellar ataxia type 6 (SCA6) is dominantly inherited neurodegenerative disease, caused by an expansion of CAG repeat encoding a polyglutamine (PolyQ) tract in the Ca v 2.1 voltage-gated calcium channel. Its key pathological features include selective degeneration of the cerebellar Purkinje cells (PCs), a common target for PolyQ-induced toxicity in various SCAs. Mutant Ca v 2.1 confers toxicity primarily through a toxic gain-of-function mechanism; however, its molecular basis remains elusive. Here, we studied the cerebellar gene expression patterns of young Sca6 -MPI 118Q/118Q knockin (KI) mice, which expressed mutant Ca v 2.1 from an endogenous locus and recapitulated many phenotypic features of human SCA6. Transcriptional signatures in the MPI 118Q/118Q mice were distinct from those in the Sca1 154Q/2Q mice, a faithful SCA1 KI mouse model. Temporal expression profiles of the candidate genes revealed that the up-regulation of genes associated with microglial activation was initiated before PC degeneration and was augmented as the disease progressed. Histological analysis of the MPI 118Q/118Q cerebellum showed the predominance of M1-like pro-inflammatory microglia and it was concomitant with elevated expression levels of tumor necrosis factor, interleukin-6, Toll-like receptor (TLR) 2 and 7. Genetic ablation of MyD88, a major adaptor protein conveying TLR signaling, altered expression patterns of M1/M2 microglial phenotypic markers in the MPI 118Q/118Q cerebellum. More importantly, it ameliorated PC loss and partially rescued motor impairments in the early disease phase. These results suggest that early neuroinflammatory response may play an important role in the pathogenesis of SCA6 and its modulation could pave the way for slowing the disease progression during the early stage of the disease.

Description: Motivation: Stem cell differentiation is largely guided by master transcriptional regulators, but it also depends on the expression of other types of genes, such as cell cycle genes, signaling genes, metabolic genes, trafficking genes, etc. Traditional approaches to understanding gene expression patterns across multiple conditions, such as principal components analysis or K-means clustering, can group cell types based on gene expression, but they do so without knowledge of the differentiation hierarchy. Hierarchical clustering can organize cell types into a tree, but in general this tree is different from the differentiation hierarchy itself. Methods: Given the differentiation hierarchy and gene expression data at each node, we construct a weighted Euclidean distance metric such that the minimum spanning tree with respect to that metric is precisely the given differentiation hierarchy. We provide a set of linear constraints that are provably sufficient for the desired construction and a linear programming approach to identify sparse sets of weights, effectively identifying genes that are most relevant for discriminating different parts of the tree. Results: We apply our method to microarray gene expression data describing 38 cell types in the hematopoiesis hierarchy, constructing a weighted Euclidean metric that uses just 175 genes. However, we find that there are many alternative sets of weights that satisfy the linear constraints. Thus, in the style of random-forest training, we also construct metrics based on random subsets of the genes and compare them to the metric of 175 genes. We then report on the selected genes and their biological functions. Our approach offers a new way to identify genes that may have important roles in stem cell differentiation. Contact: tperkins@ohri.ca Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Principal component analysis (PCA) is a basic tool often used in bioinformatics for visualization and dimension reduction. However, it is known that PCA may not consistently estimate the true direction of maximal variability in high-dimensional, low sample size settings, which are typical for molecular data. Assuming that the underlying signal is sparse, i.e. that only a fraction of features contribute to a principal component (PC), this estimation consistency can be retained. Most existing sparse PCA methods use L1-penalization, i.e. the lasso , to perform feature selection. But, the lasso is known to lack variable selection consistency in high dimensions and therefore a subsequent interpretation of selected features can give misleading results. Results: We present S4VDPCA, a sparse PCA method that incorporates a subsampling approach, namely stability selection. S4VDPCA can consistently select the truly relevant variables contributing to a sparse PC while also consistently estimate the direction of maximal variability. The performance of the S4VDPCA is assessed in a simulation study and compared to other PCA approaches, as well as to a hypothetical oracle PCA that ‘knows’ the truly relevant features in advance and thus finds optimal, unbiased sparse PCs. S4VDPCA is computationally efficient and performs best in simulations regarding parameter estimation consistency and feature selection consistency. Furthermore, S4VDPCA is applied to a publicly available gene expression data set of medulloblastoma brain tumors. Features contributing to the first two estimated sparse PCs represent genes significantly over-represented in pathways typically deregulated between molecular subgroups of medulloblastoma. Availability and implementation: Software is available at https://github.com/mwsill/s4vdpca . Contact: m.sill@dkfz.de Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Glycans play critical roles in many biological processes, and their structural diversity is key for specific protein-glycan recognition. Comparative structural studies of biological molecules provide useful insight into their biological relationships. However, most computational tools are designed for protein structure, and despite their importance, there is no currently available tool for comparing glycan structures in a sequence order- and size-independent manner. Results: A novel method, GS-align, is developed for glycan structure alignment and similarity measurement. GS-align generates possible alignments between two glycan structures through iterative maximum clique search and fragment superposition. The optimal alignment is then determined by the maximum structural similarity score, GS-score, which is size-independent. Benchmark tests against the Protein Data Bank (PDB) N -linked glycan library and PDB homologous/non-homologous N -glycoprotein sets indicate that GS-align is a robust computational tool to align glycan structures and quantify their structural similarity. GS-align is also applied to template-based glycan structure prediction and monosaccharide substitution matrix generation to illustrate its utility. Availability and implementation: http://www.glycanstructure.org/gsalign . Contact: wonpil@ku.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Impedance-based technologies are advancing methods for measuring proliferation of adherent cell cultures non-invasively and in real time. The analysis of the resulting data has so far been hampered by inappropriate computational methods and the lack of systematic data to evaluate the characteristics of the assay. Results: We used a commercially available system for impedance-based growth measurement (xCELLigence) and compared the reported cell index with data from microscopy. We found that the measured signal correlates linearly with the cell number throughout the time of an experiment with sufficient accuracy in subconfluent cell cultures. The resulting growth curves for various colon cancer cells could be well described with the empirical Richards growth model, which allows for extracting quantitative parameters (such as characteristic cycle times). We found that frequently used readouts like the cell index at a specific time or the area under the growth curve cannot be used to faithfully characterize growth inhibition. We propose to calculate the average growth rate of selected time intervals to accurately estimate time-dependent IC50 values of drugs from growth curves. Contact: nils.bluethgen@charite.de Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Hydroxy fatty acids in fresh snow samples from northern Japan: long-range atmospheric transport of Gram-negative bacteria by Asian winter monsoon Biogeosciences Discussions, 12, 13375-13397, 2015 Author(s): P. Tyagi, S. Yamamoto, and K. Kawamura Hydroxy fatty acids (FAs) in fresh snow from Sapporo, one of the heaviest snowfall regions in the world, have been studied to ascertain the airborne bacterial endotoxin concentrations and their biomass. The presence of β-hydroxy FAs (C 9 –C 28 ), constituents of Gram-negative bacteria (GNB), suggests long-range transport of soil microbes. Likewise, the occurrence of α- and ω-hydroxy FAs (C 9 –C 30 and C 9 –C 28 , respectively) in snow reveals their contribution from epicuticular waxes and soil microorganisms. Estimated endotoxin and GNB mass can aid in assessing their possible impacts on the diversity and functioning of aquatic and terrestrial ecosystems, as well as lethal effects on pedestrians through dispersal of microbes. Air mass back trajectories together with hydroxy FAs unveil their sources from Siberia, Russian Far East and North China by the Asian monsoon. This study highlights the role of fresh snow that reduces the human health risk of GNB and endotoxin by scavenging from the air.

Description: Spliceosomal introns are a hallmark of eukaryotic genes that are hypothesized to play important roles in genome evolution but have poorly understood origins. Although most introns lack sequence homology to each other, new families of spliceosomal introns that are repeated hundreds of times in individual genomes have recently been discovered in a few organisms. The prevalence and conservation of these introner elements (IEs) or introner-like elements in other taxa, as well as their evolutionary relationships to regular spliceosomal introns, are still unknown. Here, we systematically investigate introns in the widespread marine green alga Micromonas and report new families of IEs, numerous intron presence–absence polymorphisms, and potential intron insertion hot-spots. The new families enabled identification of conserved IE secondary structure features and establishment of a novel general model for repetitive intron proliferation across genomes. Despite shared secondary structure, the IE families from each Micromonas lineage bear no obvious sequence similarity to those in the other lineages, suggesting that their appearance is intimately linked with the process of speciation. Two of the new IE families come from an Arctic culture ( Micromonas Clade E2) isolated from a polar region where abundance of this alga is increasing due to climate induced changes. The same two families were detected in metagenomic data from Antarctica—a system where Micromonas has never before been reported. Strikingly high identity between the Arctic isolate and Antarctic coding sequences that flank the IEs suggests connectivity between populations in the two polar systems that we postulate occurs through deep-sea currents. Recovery of Clade E2 sequences in North Atlantic Deep Waters beneath the Gulf Stream supports this hypothesis. Our research illuminates the dynamic relationships between an unusual class of repetitive introns, genome evolution, speciation, and global distribution of this sentinel marine alga.

Description: Cellular Ca 2+ homeostasis is tightly regulated and is pivotal to life. Inositol 1,4,5-trisphosphate receptors (IP 3 Rs) and ryanodine receptors (RyRs) are the major ion channels that regulate Ca 2+ release from intracellular stores. Although these channels have been extensively investigated in multicellular organisms, an appreciation of their evolution and the biology of orthologs in unicellular organisms is largely lacking. Extensive phylogenetic analyses reveal that the IP 3 R gene superfamily is ancient and diverged into two subfamilies, IP 3 R-A and IP 3 R-B/RyR, at the dawn of Opisthokonta. IP 3 R-B/RyR further diversified into IP 3 R-B and RyR at the stem of Filozoa. Subsequent evolution and speciation of Holozoa is associated with duplication of IP 3 R-A and RyR genes, and loss of IP 3 R-B in the vertebrate lineages. To gain insight into the properties of IP 3 R important for the challenges of multicellularity, the IP 3 R-A and IP 3 R-B family orthologs were cloned from Capsaspora owczarzaki, a close unicellular relative to Metazoa (designated as CO.IP 3 R-A and CO.IP 3 R-B). Both proteins were targeted to the endoplasmic reticulum. However, CO.IP 3 R-A, but strikingly not CO.IP 3 R-B, bound IP 3 , exhibited robust Ca 2+ release activity and associated with mammalian IP 3 Rs. These data indicate strongly that CO.IP 3 R-A as an exemplar of ancestral IP 3 R-A orthologs forms bona fide IP 3 -gated channels. Notably, however, CO.IP 3 R-A appears not to be regulated by Ca 2+ , ATP or Protein kinase A-phosphorylation. Collectively, our findings explore the origin, conservation, and diversification of IP 3 R gene families and provide insight into the functionality of ancestral IP 3 Rs and the added specialization of these proteins in Metazoa.

Description: Marriage rules, the community prescriptions that dictate who an individual can or cannot marry, are extremely diverse and universally present in traditional societies. A major focus of research in the early decades of modern anthropology, marriage rules impose social and economic forces that help structure societies and forge connections between them. However, in those early anthropological studies, the biological benefits or disadvantages of marriage rules could not be determined. We revisit this question by applying a novel simulation framework and genome-wide data to explore the effects of Asymmetric Prescriptive Alliance, an elaborate set of marriage rules that has been a focus of research for many anthropologists. Simulations show that strict adherence to these marriage rules reduces genetic diversity on the autosomes, X chromosome and mitochondrial DNA, but relaxed compliance produces genetic diversity similar to random mating. Genome-wide data from the Indonesian community of Rindi, one of the early study populations for Asymmetric Prescriptive Alliance, are more consistent with relaxed compliance than strict adherence. We therefore suggest that, in practice, marriage rules are treated with sufficient flexibility to allow social connectivity without significant degradation of biological diversity.

Description: The environment has profound effects on the expression of many traits and reaction norms describe the expression dynamics of a trait across a broad range of environmental conditions. Here, we analyze gene expression in Drosophila melanogaster across four different developmental temperatures (13–29 °C). Gene expression is highly plastic with 83.3% of the genes being differentially expressed. We distinguished three components of plasticity: 1) Dynamics of gene expression intensity (sum of change), 2) direction of change, and 3) curvature of the reaction norm (linear vs. quadratic). Studying their regulatory architecture we found that all three plasticity components were most strongly affected by the number of different transcription factors (TFs) binding to the target gene. More TFs were found in genes with less expression changes across temperatures. Although the effect of microRNAs was weaker, we consistently noted a trend in the opposite direction. The most plastic genes were regulated by fewer TFs and more microRNAs than less plastic genes. Different patterns of plasticity were also reflected by their functional characterization based on gene ontology. Our results suggest that reaction norms provide an important key to understand the functional requirements of natural populations exposed to variable environmental conditions.

Description: Recent advances in paleogenomic technologies have enabled an increasingly detailed understanding of the evolutionary relationships of now-extinct mammalian taxa. However, a number of enigmatic Quaternary species have never been characterized with molecular data, often because available fossils are rare or are found in environments that are not optimal for DNA preservation. Here, we analyze paleogenomic data extracted from bones attributed to the late Pleistocene western camel, Camelops cf. hesternus, a species that was distributed across central and western North America until its extinction approximately 13,000 years ago. Despite a modal sequence length of only around 35 base pairs, we reconstructed high-coverage complete mitochondrial genomes and low-coverage partial nuclear genomes for each specimen. We find that Camelops is sister to African and Asian bactrian and dromedary camels, to the exclusion of South American camelids (llamas, guanacos, alpacas, and vicuñas). These results contradict previous morphology-based phylogenetic models for Camelops , which suggest instead a closer relationship between Camelops and the South American camelids. The molecular data imply a Late Miocene divergence of the Camelops clade from lineages that separately gave rise to the extant camels of Eurasia. Our results demonstrate the increasing capacity of modern paleogenomic methods to resolve evolutionary relationships among distantly related lineages.

Description: Recent developments in the analysis of amino acid covariation are leading to breakthroughs in protein structure prediction, protein design, and prediction of the interactome. It is assumed that observed patterns of covariation are caused by molecular coevolution, where substitutions at one site affect the evolutionary forces acting at neighboring sites. Our theoretical and empirical results cast doubt on this assumption. We demonstrate that the strongest coevolutionary signal is a decrease in evolutionary rate and that unfeasibly long times are required to produce coordinated substitutions. We find that covarying substitutions are mostly found on different branches of the phylogenetic tree, indicating that they are independent events that may or may not be attributable to coevolution. These observations undermine the hypothesis that molecular coevolution is the primary cause of the covariation signal. In contrast, we find that the pairs of residues with the strongest covariation signal tend to have low evolutionary rates, and that it is this low rate that gives rise to the covariation signal. Slowly evolving residue pairs are disproportionately located in the protein’s core, which explains covariation methods’ ability to detect pairs of residues that are close in three dimensions. These observations lead us to propose the "coevolution paradox": The strength of coevolution required to cause coordinated changes means the evolutionary rate is so low that such changes are highly unlikely to occur. As modern covariation methods may lead to breakthroughs in structural genomics, it is critical to recognize their biases and limitations.

Description: We present a modified GPU (graphics processing unit) version of MrBayes, called ta(MC) 3 (GPU MrBayes V3.1), for Bayesian phylogenetic inference on protein data sets. Our main contributions are 1) utilizing 64-bit variables, thereby enabling ta(MC) 3 to process larger data sets than MrBayes; and 2) to use Kahan summation to improve accuracy, convergence rates, and consequently runtime. Versus the current fastest software, we achieve a speedup of up to around 2.5 (and up to around 90 vs. serial MrBayes), and more on multi-GPU hardware. GPU MrBayes V3.1 is available from http://sourceforge.net/projects/mrbayes-gpu/ .

Description: Understanding the genetic and molecular bases of the ability to distinguish self from nonself (allorecognition) and mechanisms underlying evolution of allorecognition systems is an important endeavor for understanding cases where it becomes dysfunctional, such as in autoimmune disorders. In filamentous fungi, allorecognition can result in vegetative or heterokaryon incompatibility, which is a type of programmed cell death that occurs following fusion of genetically different cells. Allorecognition is genetically controlled by het loci, with coexpression of any combination of incompatible alleles triggering vegetative incompatibility. Herein, we identified, characterized, and inferred the evolutionary history of candidate het loci in the filamentous fungus Neurospora crassa . As characterized het loci encode proteins carrying an HET domain, we annotated HET domain genes in 25 isolates from a natural population along with the N. crassa reference genome using resequencing data. Because allorecognition systems can be affected by frequency-dependent selection favoring rare alleles (i.e., balancing selection), we mined resequencing data for HET domain loci whose alleles displayed elevated levels of variability, excess of intermediate frequency alleles, and deep gene genealogies. From these analyses, 34 HET domain loci were identified as likely to be under balancing selection. Using transformation, incompatibility assays and genetic analyses, we determined that one of these candidates functioned as a het locus ( het-e ). The het-e locus has three divergent allelic groups that showed signatures of positive selection, intra- and intergroup recombination, and trans-species polymorphism. Our findings represent a compelling case of balancing selection functioning on multiple alleles across multiple loci potentially involved in allorecognition.

Description: Scoring the impact of noncoding variation on the function of cis -regulatory regions, on their chromatin state, and on the qualitative and quantitative expression levels of target genes is a fundamental problem in evolutionary genomics. A particular challenge is how to model the divergence of quantitative traits and to identify relationships between the changes across the different levels of the genome, the chromatin activity landscape, and the transcriptome. Here, we examine the use of the Ornstein–Uhlenbeck (OU) model to infer selection at the level of predicted cis -regulatory modules (CRMs), and link these with changes in transcription factor binding and chromatin activity. Using publicly available cross-species ChIP-Seq and STARR-Seq data we show how OU can be applied genome-wide to identify candidate transcription factors for which binding site and CRM turnover is correlated with changes in regulatory activity. Next, we profile open chromatin in the developing eye across three Drosophila species. We identify the recognition motifs of the chromatin remodelers, Trithorax-like and Grainyhead as mostly correlating with species-specific changes in open chromatin. In conclusion, we show in this study that CRM scores can be used as quantitative traits and that motif discovery approaches can be extended towards more complex models of divergence.

Description: How populations that inhabit the same geographical area become genetically differentiated is not clear. To investigate this, we characterized phenotypic and genetic differences between two populations of Saccharomyces cerevisiae that in some cases inhabit the same environment but show relatively little gene flow. We profiled stress sensitivity in a group of vineyard isolates and a group of oak-soil strains and found several niche-related phenotypes that distinguish the populations. We performed bulk-segregant mapping on two of the distinguishing traits: The vineyard-specific ability to grow in grape juice and oak-specific tolerance to the cell wall damaging drug Congo red. To implicate causal genes, we also performed a chemical genomic screen in the lab-strain deletion collection and identified many important genes that fell under quantitative trait loci peaks. One gene important for growth in grape juice and identified by both the mapping and the screen was SSU1 , a sulfite-nitrite pump implicated in wine fermentations. The beneficial allele is generated by a known translocation that we reasoned may also serve as a genetic barrier. We found that the translocation is prevalent in vineyard strains, but absent in oak strains, and presents a postzygotic barrier to spore viability. Furthermore, the translocation was associated with a fitness cost to the rapid growth rate seen in oak-soil strains. Our results reveal the translocation as a dual-function locus that enforces ecological differentiation while producing a genetic barrier to gene flow in these sympatric populations.

Description: Species tree reconstruction has been a subject of substantial research due to its central role across biology and medicine. A species tree is often reconstructed using a set of gene trees or by directly using sequence data. In either of these cases, one of the main confounding phenomena is the discordance between a species tree and a gene tree due to evolutionary events such as duplications and losses. Probabilistic methods can resolve the discordance by coestimating gene trees and the species tree but this approach poses a scalability problem for larger data sets. We present MixTreEM-DLRS: A two-phase approach for reconstructing a species tree in the presence of gene duplications and losses. In the first phase, MixTreEM, a novel structural expectation maximization algorithm based on a mixture model is used to reconstruct a set of candidate species trees, given sequence data for monocopy gene families from the genomes under study. In the second phase, PrIME-DLRS, a method based on the DLRS model (Åkerborg O, Sennblad B, Arvestad L, Lagergren J. 2009. Simultaneous Bayesian gene tree reconstruction and reconciliation analysis. Proc Natl Acad Sci U S A. 106(14):5714–5719), is used for selecting the best species tree. PrIME-DLRS can handle multicopy gene families since DLRS, apart from modeling sequence evolution, models gene duplication and loss using a gene evolution model (Arvestad L, Lagergren J, Sennblad B. 2009. The gene evolution model and computing its associated probabilities. J ACM. 56(2):1–44). We evaluate MixTreEM-DLRS using synthetic and biological data, and compare its performance with a recent genome-scale species tree reconstruction method PHYLDOG ( Boussau B, Szöllősi GJ, Duret L, Gouy M, Tannier E, Daubin V. 2013 . Genome-scale coestimation of species and gene trees. Genome Res. 23(2):323–330) as well as with a fast parsimony-based algorithm Duptree (Wehe A, Bansal MS, Burleigh JG, Eulenstein O. 2008. Duptree: a program for large-scale phylogenetic analyses using gene tree parsimony. Bioinformatics 24(13):1540–1541). Our method is competitive with PHYLDOG in terms of accuracy and runs significantly faster and our method outperforms Duptree in accuracy. The analysis constituted by MixTreEM without DLRS may also be used for selecting the target species tree, yielding a fast and yet accurate algorithm for larger data sets. MixTreEM is freely available at http://prime.scilifelab.se/mixtreem/ .

Description: Contrasting with birds and mammals, poikilothermic vertebrates often have homomorphic sex chromosomes, possibly resulting from high rates of sex-chromosome turnovers and/or occasional X–Y recombination. Strong support for the latter mechanism was provided by four species of European tree frogs, which inherited from a common ancestor (~5 Ma) the same pair of homomorphic sex chromosomes (linkage group 1, LG1), harboring the candidate sex-determining gene Dmrt1. Here, we test sex linkage of LG1 across six additional species of the Eurasian Hyla radiation with divergence times ranging from 6 to 40 Ma. LG1 turns out to be sex linked in six of nine resolved cases. Mapping the patterns of sex linkage to the Hyla phylogeny reveals several transitions in sex-determination systems within the last 10 My, including one switch in heterogamety. Phylogenetic trees of DNA sequences along LG1 are consistent with occasional X–Y recombination in all species where LG1 is sex linked. These patterns argue against one of the main potential causes for turnovers, namely the accumulation of deleterious mutations on nonrecombining chromosomes. Sibship analyses show that LG1 recombination is strongly reduced in males from most species investigated, including some in which it is autosomal. Intrinsically low male recombination might facilitate the evolution of male heterogamety, and the presence of important genes from the sex-determination cascade might predispose LG1 to become a sex chromosome.

Description: It has been hypothesised that positive associations between age and levels of oxidative stress-generated damage to DNA may be related to an age-dependent decline in DNA repair activity. The objective of this study was to investigate the association between age and repair activity of oxidatively damaged DNA in peripheral blood mononuclear cells (PBMCs). We isolated PBMCs from subjects aged 18–83 years, as part of a health survey of the Danish population that focussed on lifestyle factors. The level of DNA repair activity was measured as incisions on potassium bromate-damaged DNA by the comet assay. There was an inverse association between age and DNA repair activity with a 0.65% decline in activity per year from age 18 to 83 (95% confidence interval: 0.16–1.14% per year). Univariate regression analysis also indicated inverse associations between DNA repair activity and waist-hip ratio ( P 〈 0.05) and plasma concentrations of glycosylated hemoglobin ( P = 0.07). However, multivariate regression analysis only showed an inverse association between age and DNA repair activity ( P 〈 0.05), indicating that the decline in repair activity was not mediated by metabolic risk factors. In summary, the results show an inverse association between age and DNA repair activity of oxidatively damaged DNA.

Description: Exposure to traffic-related particulate matter (PM) has been associated with increased risk of lung disease, cancer and cardiovascular disease especially in elderly and overweight subjects. The proposed mechanisms involve intracellular production of reactive oxygen species (ROS), inflammation and oxidation-induced DNA damage studied mainly in young normal-weight subjects. We performed a controlled cross-over, randomised, single-blinded, repeated-measure study where 60 healthy subjects (25 males and 35 females) with age 55–83 years and body mass index above 25kg/m 2 were exposed for 5h to either particle-filtered or sham-filtered air from a busy street with number of concentrations and PM 2.5 levels of 1800/cm 3 versus 23 000/cm 3 and 3 µg/m 3 versus 24 µg/m 3 , respectively. Peripheral blood mononuclear cells (PBMCs) were collected and assayed for production of ROS with and without ex vivo exposure to nanosized carbon black as well as expression of genes related to inflammation ( chemokine (C-C motif) ligand 2 , interleukin-8 and tumour necrosis factor ), oxidative stress response ( heme oxygenase (decycling)-1 ) and DNA repair ( oxoguanine DNA glycosylase ). DNA strand breaks and oxidised purines were assayed by the alkaline comet assay. No statistically significant differences were found for any biomarker immediately after exposure to PM from urban street air although strand breaks and oxidised purines combined were significantly associated with the particle number concentration during exposure. In conclusion, 5h of controlled exposure to PM from urban traffic did not change the gene expression related to inflammation, oxidative stress or DNA repair, ROS production or oxidatively damaged DNA in PBMCs from elderly overweight human subjects.

Description: Ionising radiation causes free radical–mediated damage in cellular DNA. This damage is manifested as chromosomal aberrations and micronuclei (MN) in proliferating cells. Sesamol, present in sesame seeds, has the potential to scavenge free radicals; therefore, it can reduce radiation-induced cytogenetic damage in cells. The aim of this study was to investigate the radioprotective potential of sesamol in bone marrow cells of mice and related haematopoietic system against radiation-induced genotoxicity. A comparative study with melatonin was designed for assessing the radioprotective potential of sesamol. C57BL/6 mice were administered intraperitoneally with either sesamol or melatonin (10 and 20mg/kg body weight) 30min prior to 2-Gy whole-body irradiation (WBI) and sacrificed after 24h. Total chromosomal aberrations (TCA), MN and cell cycle analyses were performed using bone marrow cells. The comet assay was performed on bone marrow cells, splenocytes and lymphocytes. Blood was drawn to study haematological parameters. Prophylactic doses of sesamol (10 and 20mg/kg) in irradiated mice reduced TCA and micronucleated polychromatic erythrocyte frequency in bone marrow cells by 57% and 50%, respectively, in comparison with radiation-only groups. Sesamol-reduced radiation-induced apoptosis and facilitated cell proliferation. In the comet assay, sesamol (20mg/kg) treatment reduced radiation-induced comets (% DNA in tail) compared with radiation only ( P 〈 0.05). Sesamol also increased granulocyte populations in peripheral blood similar to melatonin. Overall, the radioprotective efficacy of sesamol was found to be similar to that of melatonin. Sesamol treatment also showed recovery of relative spleen weight at 24h of WBI. The results strongly suggest the radioprotective efficacy of sesamol in the haematopoietic system of mice.

Description: Individualized treatment rules recommend treatments on the basis of individual patient characteristics. A high-quality treatment rule can produce better patient outcomes, lower costs and less treatment burden. If a treatment rule learned from data is to be used to inform clinical practice or provide scientific insight, it is crucial that it be interpretable; clinicians may be unwilling to implement models they do not understand, and black-box models may not be useful for guiding future research. The canonical example of an interpretable prediction model is a decision tree. We propose a method for estimating an optimal individualized treatment rule within the class of rules that are representable as decision trees. The class of rules we consider is interpretable but expressive. A novel feature of this problem is that the learning task is unsupervised, as the optimal treatment for each patient is unknown and must be estimated. The proposed method applies to both categorical and continuous treatments and produces favourable marginal mean outcomes in simulation experiments. We illustrate it using data from a study of major depressive disorder.

Description: Sufficient dimension reduction in regression aims to reduce the predictor dimension by replacing the original predictors with some set of linear combinations of them without loss of information. Numerous dimension reduction methods have been developed based on this paradigm. However, little effort has been devoted to diagnostic studies within the context of dimension reduction. In this paper we introduce methods to check goodness-of-fit for a given dimension reduction subspace. The key idea is to extend the so-called distance correlation to measure the conditional dependence relationship between the covariates and the response given a reduction subspace. Our methods require minimal assumptions, which are usually much less restrictive than the conditions needed to justify the original methods. Asymptotic properties of the test statistic are studied. Numerical examples demonstrate the effectiveness of the proposed approach.

Description: We propose a geometric framework to assess sensitivity of Bayesian procedures to modelling assumptions based on the nonparametric Fisher–Rao metric. While the framework is general, the focus of this article is on assessing local and global robustness in Bayesian procedures with respect to perturbations of the likelihood and prior, and on the identification of influential observations. The approach is based on a square-root representation of densities, which enables analytical computation of geodesic paths and distances, facilitating the definition of naturally calibrated local and global discrepancy measures. An important feature of our approach is the definition of a geometric $\epsilon$ -contamination class of sampling distributions and priors via intrinsic analysis on the space of probability density functions. We demonstrate the applicability of our framework to generalized mixed-effects models and to directional and shape data.

Description: With the discovery of an increasing number of causal genes for complex human disorders, it is crucial to assess the genetic risk of disease onset for individuals who are carriers of these causal mutations and to compare the distribution of the age-at-onset for such individuals with the distribution for noncarriers. In many genetic epidemiological studies that aim to estimate causal gene effect on disease, the age-at-onset of disease is subject to censoring. In addition, the mutation carrier or noncarrier status of some individuals may be unknown, due to the high cost of in-person ascertainment by collecting DNA samples or because of the death of older individuals. Instead, the probability of such individuals’ mutation status can be obtained from various other sources. When mutation status is missing, the available data take the form of censored mixture data. Recently, various methods have been proposed for risk estimation using such data, but none is efficient for estimating a nonparametric distribution. We propose a fully efficient sieve maximum likelihood estimation method, in which we estimate the logarithm of the hazard ratio between genetic mutation groups using B-splines, while applying nonparametric maximum likelihood estimation to the reference baseline hazard function. Our estimator can be calculated via an expectation-maximization algorithm which is much faster than existing methods. We show that our estimator is consistent and semiparametrically efficient and establish its asymptotic distribution. Simulation studies demonstrate the superior performance of the proposed method, which is used to estimate the distribution of the age-at-onset of Parkinson's disease for carriers of mutations in the leucine-rich repeat kinase 2, LRRK2, gene.

Description: Smoothing splines provide flexible nonparametric regression estimators. However, the high computational cost of smoothing splines for large datasets has hindered their wide application. In this article, we develop a new method, named adaptive basis sampling, for efficient computation of smoothing splines in super-large samples. Except for the univariate case where the Reinsch algorithm is applicable, a smoothing spline for a regression problem with sample size n can be expressed as a linear combination of n basis functions and its computational complexity is generally O ( n 3 ). We achieve a more scalable computation in the multivariate case by evaluating the smoothing spline using a smaller set of basis functions, obtained by an adaptive sampling scheme that uses values of the response variable. Our asymptotic analysis shows that smoothing splines computed via adaptive basis sampling converge to the true function at the same rate as full basis smoothing splines. Using simulation studies and a large-scale deep earth core-mantle boundary imaging study, we show that the proposed method outperforms a sampling method that does not use the values of response variables.

Description: We propose tests for nonlinear serial dependence in time series under the null hypothesis of general linear dependence, in contrast to the more widely studied null hypothesis of independence. The approach is based on combining an entropy dependence metric, which possesses many desirable properties and is used as a test statistic, with a suitable extension of surrogate data methods, a class of Monte Carlo distribution-free tests for nonlinearity, and a smoothed sieve bootstrap scheme. We show how, in the same way as the autocorrelation function is used for linear models, our tests can in principle be employed to detect the lags at which a significant nonlinear relationship is present. We prove the asymptotic validity of the proposed procedures and the corresponding inferences. The small-sample performance of the tests in terms of power and size is assessed through a simulation study. Applications to real datasets of different kinds are also presented.

Description: Contamination caused by outliers is inevitable in data analysis, and robust statistical methods are often needed. In this paper we develop a new approach for robust data analysis on the basis of scoring rules. A scoring rule is a discrepancy measure to assess the quality of probabilistic forecasts. We propose a simple method of estimating not only parameters in the statistical model but also the contamination ratio, i.e., the ratio of outliers. The outliers are detected based on the estimated contamination ratio. For this purpose, we use scoring rules with extended statistical models called unnormalized models. Regression problems are also considered. We study complex heterogeneous contamination wherein the contamination ratio in a response variable may depend on covariate variables, and propose a simple method to estimate a robust regression function and expected contamination ratio. Simulation studies demonstrate the effectiveness of our method.

Description: The sample covariance matrix, which is well known to be highly nonrobust, plays a central role in many classical multivariate statistical methods. A popular way of making such multivariate methods more robust is to replace the sample covariance matrix with some robust scatter matrix. The aim of this paper is to point out that multivariate methods often require that certain properties of the covariance matrix hold also for the robust scatter matrix in order for the corresponding robust plug-in method to be a valid approach, but that not all scatter matrices possess the desired properties. Plug-in methods for independent components analysis, observational regression and graphical modelling are considered in more detail. For each case, it is shown that replacing the sample covariance matrix with a symmetrized robust scatter matrix yields a valid robust multivariate procedure.

Description: So-called big data are likely to have complex structure, in particular implying that estimates of precision obtained by applying standard statistical procedures are likely to be misleading, even if the point estimates of parameters themselves may be reasonably satisfactory. While this possibility is best explored in the context of each special case, here we outline a fairly general representation of the accretion of error in large systems and explore the possible implications for the estimation of regression coefficients. The discussion raises issues broadly parallel to the distinction between short-range and long-range dependence in time series theory.

Description: This paper extends the classical two-regime threshold autoregressive model by introducing hysteresis to its regime-switching structure, which leads to a new model: the hysteretic autoregressive model. The proposed model enjoys the piecewise linear structure of a threshold model but has a more flexible regime switching mechanism. A sufficient condition is given for geometric ergodicity. Conditional least squares estimation is discussed, and the asymptotic distributions of its estimators and information criteria for model selection are derived. Simulation results and an example support the model.

Description: This paper points out an error in Davidov and Iliopoulos's ( Biometrika 100 , 778–80) proof of convergence of an iterative algorithm for the proportional likelihood ratio model. It is shown that the iterative algorithm increases the likelihood in each iteration and converges under mild additional conditions when the odds ratio function is bounded.

Description: Nitrogen cycling in the subsurface biosphere: nitrate isotopes in porewaters underlying the oligotrophic North Atlantic Biogeosciences Discussions, 12, 13545-13591, 2015 Author(s): S. D. Wankel, C. Buchwald, W. Ziebis, C. B. Wenk, and M. F. Lehmann Nitrogen (N) is a key component of fundamental biomolecules. Hence, the cycling and availability of N is a central factor governing the extent of ecosystems across the Earth. In the organic-lean sediment porewaters underlying the oligotrophic ocean, where low levels of microbial activity persist despite limited organic matter delivery from overlying water, the extent and modes of nitrogen transformations have not been widely investigated. Here we use the N and oxygen (O) isotopic composition of porewater nitrate (NO 3 − ) from a site in the oligotrophic North Atlantic (IODP) to determine the extent and magnitude of microbial nitrate production (via nitrification) and consumption (via denitrification). We find that NO 3 − accumulates far above bottom seawater concentrations (∼ 21 μM) throughout the sediment column (up to ∼ 50 μM) down to the oceanic basement as deep as 90 mbsf, reflecting the predominance of aerobic nitrification/remineralization within the deep marine sediments. Large changes in the δ 15 N and δ 18 O of nitrate, however, reveal variable influence of nitrate respiration across the three sites. We use an inverse porewater diffusion–reaction model, constrained by the N and O isotope systematics of nitrification and denitrification and the porewater NO 3 − isotopic composition, to estimate rates of nitrification and denitrification throughout the sediment column. Results indicate variability of reaction rates across and within the three boreholes that are generally consistent with the differential distribution of dissolved oxygen at this site, though not necessarily with the canonical view of how redox thresholds separate nitrate regeneration from dissimilative consumption spatially. That is, we provide isotope evidence for expanded zones of co-ocurring nitrification and denitrification. The isotope biogeochemical modeling also yielded estimates for the δ 15 N and δ 18 O of newly produced nitrate (δ 15 N NTR and δ 18 O NTR ), as well as the isotope effect for denitrification ( 15 ϵ DNF ), parameters with high relevance to global ocean models of N cycling. Estimated values of δ 15 N NTR were generally lower than previously reported δ 15 N values for sinking PON in this region. We suggest that these values can be related to sedimentary N-fixation and remineralization of the newly fixed organic N. Values of δ 18 O NTR generally ranged between −2.8 and 0.0 ‰, consistent with recent estimates based on lab cultures of nitrifying bacteria. Notably, some δ 18 O NTR values were elevated, suggesting incorporation of 18 O-enriched dissolved oxygen during nitrification, and possibly indicating a tight coupling of NH 4 + and NO 2 − oxidation in this metabolically sluggish environment. Our findings indicate that the production of organic matter by in situ autotrophy (e.g., nitrification, nitrogen fixation) supply a large fraction of the biomass and organic substrate for heterotrophy in these sediments, supplementing the small organic matter pool derived from the overlying euphotic zone. This work sheds new light on an active nitrogen cycle operating, despite exceedingly low carbon inputs, in the deep sedimentary biosphere.

Description: Spring bloom onset in the Nordic Seas Biogeosciences Discussions, 12, 13631-13673, 2015 Author(s): A. Mignot, R. Ferrari, and K. A. Mork The North Atlantic spring bloom is a massive annual growth event of marine phytoplankton, tiny free-floating algae that form the base of the ocean's food web and generates a large fraction of the global primary production of organic matter. The conditions that trigger the onset of the spring bloom in the Nordic Seas, at the northern edge of the North Atlantic, are studied using in-situ data from five bio-optical floats released above the Arctic Circle. It is often assumed that spring blooms start as soon as phytoplankton cells daily irradiance is sufficiently abundant that division rates exceed losses. The bio-optical float data instead suggest the tantalizing hypothesis that Nordic Seas blooms start when the photoperiod, the number of daily light hours experienced by phytoplankton, exceeds a critical value, independently of division rates. This bloom behavior may be explained by realizing that photosynthesis is impossible during polar nights and phytoplankton enters in a dormant stage in winter, only to be awaken by a photoperiodic trigger. While the first accumulation of biomass recorded by the bio-optical floats is consistent with the photoperiod hypothesis, it is possible that some biomass accumulation started before the critical photoperiod but at levels too low to be detected by the fluorometers. Thus more precise observations are needed to test the photoperiod hypothesis.

Description: Unusual biogenic calcite structures in two shallow lakes, James Ross Island, Antarctica Biogeosciences Discussions, 12, 13593-13629, 2015 Author(s): J. Elster, L. Nedbalová, R. Vodrážka, K. Láska, J. Haloda, and J. Komárek The floors of two shallow endorheic lakes, located on volcanic surfaces on James Ross Island, are covered with calcareous organosedimentary structures. Their biological and chemical composition, lake water characteristics, and seasonal variability of the thermal regime are introduced. The lakes are frozen down to the bottom eight-nine months per year and their water chemistry is characterized by low conductivity and neutral to slightly alkaline pH. The photosynthetic microbial mat is composed of filamentous cyanobacteria and microalgae that are considered to be Antarctic endemic species. The mucilaginous black biofilm is covered by green spots formed by a green microalga and the macroscopic structures are packed together with fine material. Thin sections consist of rock substrate, soft biofilm, calcite spicules and mineral grains originating from different sources. The morphology of the spicules is typical of calcium carbonate monocrystals having a layered structure and worn surface, which reflect growth and degradation processes. The spicules chemical composition and structure correspond to pure calcite. Lakes age, altitude, morphometry, geomorphological and hydrological stability, including low sedimentation rates, together with thermal regime predispose the existence of this community. We hypothesize that the precipitation of calcite is connected with the photosynthetic activity of the green microalgae that were not recorded in any other lake in the region. This study has shown that the unique community producing biogenic calcite spicules is quite different to any yet described.

Description: Motivation: The storage and transmission of high-throughput sequencing data consumes significant resources. As our capacity to produce such data continues to increase, this burden will only grow. One approach to reduce storage and transmission requirements is to compress this sequencing data. Results: We present a novel technique to boost the compression of sequencing that is based on the concept of bucketing similar reads so that they appear nearby in the file. We demonstrate that, by adopting a data-dependent bucketing scheme and employing a number of encoding ideas, we can achieve substantially better compression ratios than existing de novo sequence compression tools, including other bucketing and reordering schemes. Our method, Mince, achieves up to a 45% reduction in file sizes (28% on average) compared with existing state-of-the-art de novo compression schemes. Availability and implementation : Mince is written in C++11, is open source and has been made available under the GPLv3 license. It is available at http://www.cs.cmu.edu/~ckingsf/software/mince . Contact: carlk@cs.cmu.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : Current methods for motif discovery from chromatin immunoprecipitation followed by sequencing (ChIP-seq) data often identify non-targeted transcription factor (TF) motifs, and are even further limited when peak sequences are similar due to common ancestry rather than common binding factors. The latter aspect particularly affects a large number of proteins from the Cys 2 His 2 zinc finger (C2H2-ZF) class of TFs, as their binding sites are often dominated by endogenous retroelements that have highly similar sequences. Here, we present recognition code-assisted discovery of regulatory elements (RCADE) for motif discovery from C2H2-ZF ChIP-seq data. RCADE combines predictions from a DNA recognition code of C2H2-ZFs with ChIP-seq data to identify models that represent the genuine DNA binding preferences of C2H2-ZF proteins. We show that RCADE is able to identify generalizable binding models even from peaks that are exclusively located within the repeat regions of the genome, where state-of-the-art motif finding approaches largely fail. Availability and implementation: RCADE is available as a webserver and also for download at http://rcade.ccbr.utoronto.ca/ . Supplementary information: Supplementary data are available at Bioinformatics online. Contact: t.hughes@utoronto.ca

Description: Motivation: Phylogenetic estimates from published studies can be archived using general platforms like Dryad (Vision, 2010) or TreeBASE (Sanderson et al. , 1994). Such services fulfill a crucial role in ensuring transparency and reproducibility in phylogenetic research. However, digital tree data files often require some editing (e.g. rerooting) to improve the accuracy and reusability of the phylogenetic statements. Furthermore, establishing the mapping between tip labels used in a tree and taxa in a single common taxonomy dramatically improves the ability of other researchers to reuse phylogenetic estimates. As the process of curating a published phylogenetic estimate is not error-free, retaining a full record of the provenance of edits to a tree is crucial for openness, allowing editors to receive credit for their work and making errors introduced during curation easier to correct. Results : Here, we report the development of software infrastructure to support the open curation of phylogenetic data by the community of biologists. The backend of the system provides an interface for the standard database operations of creating, reading, updating and deleting records by making commits to a git repository. The record of the history of edits to a tree is preserved by git’s version control features. Hosting this data store on GitHub ( http://github.com/ ) provides open access to the data store using tools familiar to many developers. We have deployed a server running the ‘phylesystem-api’, which wraps the interactions with git and GitHub. The Open Tree of Life project has also developed and deployed a JavaScript application that uses the phylesystem-api and other web services to enable input and curation of published phylogenetic statements. Availability and implementation : Source code for the web service layer is available at https://github.com/OpenTreeOfLife/phylesystem-api . The data store can be cloned from: https://github.com/OpenTreeOfLife/phylesystem . A web application that uses the phylesystem web services is deployed at http://tree.opentreeoflife.org/curator . Code for that tool is available from https://github.com/OpenTreeOfLife/opentree . Contact : mtholder@gmail.com

Description: Motivation: Chromatin Immunoprecipitation followed by sequencing (ChIP-seq) detects genome-wide DNA–protein interactions and chromatin modifications, returning enriched regions (ERs), usually associated with a significance score. Moderately significant interactions can correspond to true, weak interactions, or to false positives; replicates of a ChIP-seq experiment can provide co-localised evidence to decide between the two cases. We designed a general methodological framework to rigorously combine the evidence of ERs in ChIP-seq replicates, with the option to set a significance threshold on the repeated evidence and a minimum number of samples bearing this evidence. Results : We applied our method to Myc transcription factor ChIP-seq datasets in K562 cells available in the ENCODE project. Using replicates, we could extend up to 3 times the ER number with respect to single-sample analysis with equivalent significance threshold. We validated the ‘rescued’ ERs by checking for the overlap with open chromatin regions and for the enrichment of the motif that Myc binds with strongest affinity; we compared our results with alternative methods (IDR and jMOSAiCS), obtaining more validated peaks than the former and less peaks than latter, but with a better validation. Availability and implementation : An implementation of the proposed method and its source code under GPLv3 license are freely available at http://www.bioinformatics.deib.polimi.it/MSPC/ and http://mspc.codeplex.com/ , respectively. Contact : marco.morelli@iit.it Supplementary information: Supplementary Material are available at Bioinformatics online.

Description: : We announce the release of kSNP3.0, a program for SNP identification and phylogenetic analysis without genome alignment or the requirement for reference genomes. kSNP3.0 is a significantly improved version of kSNP v2. Availability and implementation : kSNP3.0 is implemented as a package of stand-alone executables for Linux and Mac OS X under the open-source BSD license. The executable packages, source code and a full User Guide are freely available at https://sourceforge.net/projects/ksnp/files/ Contact: barryghall@gmail.com

Description: Motivation: We have created an R package named phylogeo that provides a set of geographic utilities for sequencing-based microbial ecology studies. Although the geographic location of samples is an important aspect of environmental microbiology, none of the major software packages used in processing microbiome data include utilities that allow users to map and explore the spatial dimension of their data. phylogeo solves this problem by providing a set of plotting and mapping functions that can be used to visualize the geographic distribution of samples, to look at the relatedness of microbiomes using ecological distance, and to map the geographic distribution of particular sequences. By extending the popular phyloseq package and using the same data structures and command formats, phylogeo allows users to easily map and explore the geographic dimensions of their data from the R programming language. Availability and Implementation: phylogeo is documented and freely available http://zachcp.github.io/phylogeo Contact : zcharlop@rockefeller.edu

Description: : Gener is a development module for programming chemical controllers based on DNA strand displacement. Gener is developed with the aim of providing a simple interface that minimizes the opportunities for programming errors: Gener allows the user to test the computations of the DNA programs based on a simple two-domain strand displacement algebra, the minimal available so far. The tool allows the user to perform stepwise computations with respect to the rules of the algebra as well as exhaustive search of the computation space with different options for exploration and visualization. Gener can be used in combination with existing tools, and in particular, its programs can be exported to Microsoft Research’s DSD tool as well as to LaTeX. Availability and implementation : Gener is available for download at the Cosbi website at http://www.cosbi.eu/research/prototypes/gener as a windows executable that can be run on Mac OS X and Linux by using Mono. Contact : ozan@cosbi.eu

Description: Motivation: Molecular dynamics simulations provide atomic insight into the physicochemical characteristics of lipid membranes and hence, a wide range of force field families capable of modelling various lipid types have been developed in recent years. To model membranes in a biologically realistic lipid composition, simulation systems containing multiple different lipids must be assembled. Results: We present a new web service called MemGen that is capable of setting up simulation systems of heterogenous lipid membranes. MemGen is not restricted to certain lipid force fields or lipid types, but instead builds membranes from uploaded structure files which may contain any kind of amphiphilic molecule. MemGen works with any all-atom or united-atom lipid representation. Availability and implementation: MemGen is freely available without registration at http://memgen.uni-goettingen.de . Contact: jhub@gwdg.de Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Environmental control of flowering allows plant reproduction to occur under optimal conditions and facilitates adaptation to different locations. At high latitude, flowering of many plants is controlled by seasonal changes in day length. The photoperiodic flowering pathway confers this response in the Brassicaceae, which colonized temperate latitudes after divergence from the Cleomaceae, their subtropical sister family. The CONSTANS (CO) transcription factor of Arabidopsis thaliana , a member of the Brassicaceae, is central to the photoperiodic flowering response and shows characteristic patterns of transcription required for day-length sensing. CO is believed to be widely conserved among flowering plants; however, we show that it arose after gene duplication at the root of the Brassicaceae followed by divergence of transcriptional regulation and protein function. CO has two close homologs, CONSTANS-LIKE1 ( COL1 ) and COL2 , which are related to CO by tandem duplication and whole-genome duplication, respectively. The single CO homolog present in the Cleomaceae shows transcriptional and functional features similar to those of COL1 and COL2 , suggesting that these were ancestral. We detect cis -regulatory and codon changes characteristic of CO and use transgenic assays to demonstrate their significance in the day-length-dependent activation of the CO target gene FLOWERING LOCUS T. Thus, the function of CO as a potent photoperiodic flowering switch evolved in the Brassicaceae after gene duplication. The origin of CO may have contributed to the range expansion of the Brassicaceae and suggests that in other families CO genes involved in photoperiodic flowering arose by convergent evolution.

Description: The high regulatory complexity of vertebrates has been related to two rounds of whole genome duplication (2R-WGD) that occurred before the divergence of the major vertebrate groups. Following these events, many developmental transcription factors (TFs) were retained in multiple copies and subsequently specialized in diverse functions, whereas others reverted to their singleton state. TFs are known to be generally rich in amino acid repeats or low-complexity regions (LCRs), such as polyalanine or polyglutamine runs, which can evolve rapidly and potentially influence the transcriptional activity of the protein. Here we test the hypothesis that LCRs have played a major role in the diversification of TF gene duplicates. We find that nearly half of the TF gene families originated during the 2R-WGD contains LCRs. The number of gene duplicates with LCRs is 155 out of 550 analyzed (28%), about twice as many as the number of single copy genes with LCRs (15 out of 115, 13%). In addition, duplicated TFs preferentially accumulate certain LCR types, the most prominent of which are alanine repeats. We experimentally test the role of alanine-rich LCRs in two different TF gene families, PHOX2A/PHOX2B and LHX2/LHX9. In both cases, the presence of the alanine-rich LCR in one of the copies (PHOX2B and LHX2) significantly increases the capacity of the TF to activate transcription. Taken together, the results provide strong evidence that LCRs are important driving forces of evolutionary change in duplicated genes.

Description: Studies of species with continental distributions continue to identify intraspecific lineages despite continuous habitat. Lineages may form due to isolation by distance, adaptation, divergence across barriers, or genetic drift following range expansion. We investigated lineage diversification and admixture within American black bears ( Ursus americanus ) across their range using 22 k single nucleotide polymorphisms and mitochondrial DNA sequences. We identified three subcontinental nuclear clusters which we further divided into nine geographic regions: Alaskan (Alaska-East), eastern (Central Interior Highlands, Great Lakes, Northeast, Southeast), and western (Alaska-West, West, Pacific Coast, Southwest). We estimated that the western cluster diverged 67 ka, before eastern and Alaskan divergence 31 ka; these divergence dates contrasted with those from the mitochondrial genome where clades A and B diverged 1.07 Ma, and clades A-east and A-west diverged 169 ka. We combined estimates of divergence timing with hindcast species distribution models to infer glacial refugia for the species in Beringia, Pacific Northwest, Southwest, and Southeast. Our results show a complex arrangement of admixture due to expansion out of multiple refugia. The delineation of the genomic population clusters was inconsistent with the ranges for 16 previously described subspecies. Ranges for U. a. pugnax and U. a. cinnamomum were concordant with admixed clusters, calling into question how to order taxa below the species level. Additionally, our finding that U. a. floridanus has not diverged from U. a. americanus also suggests that morphology and genetics should be reanalyzed to assess taxonomic designations relevant to the conservation management of the species.

Description: Deinococcus bacteria are extremely resistant to radiation, oxidation, and desiccation. Resilience to these factors has been suggested to be due to enhanced damage prevention and repair mechanisms, as well as highly efficient antioxidant protection systems. Here, using mutation-accumulation experiments, we find that the GC-rich Deinococcus radiodurans has an overall background genomic mutation rate similar to that of E. coli , but differs in mutation spectrum, with the A/T to G/C mutation rate (based on a total count of 88 A:T-〉G:C transitions and 82 A:T-〉C:G transversions) per site per generation higher than that in the other direction (based on a total count of 157 G:C-〉A:T transitions and 33 G:C-〉T:A transversions). We propose that this unique spectrum is shaped mainly by the abundant uracil DNA glycosylases reducing G:C-〉A:T transitions, adenine methylation elevating A:T-〉C:G transversions, and absence of cytosine methylation decreasing G:C-〉A:T transitions. As opposed to the greater than 100 x elevation of the mutation rate in MMR – (DNA Mismatch Repair deficient) strains of most other organisms, MMR – D. radiodurans only exhibits a 4-fold elevation, raising the possibility that other DNA repair mechanisms compensate for a relatively low-efficiency DNA MMR pathway. As D. radiodurans has plentiful insertion sequence (IS) elements in the genome and the activities of IS elements are rarely directly explored, we also estimated the insertion (transposition) rate of the IS elements to be 2.50 x 10 –3 per genome per generation in the wild-type strain; knocking out MMR did not elevate the IS element insertion rate in this organism.

Description: Because mating systems affect population genetics and ecology, they are expected to impact the molecular evolution of species. Self-fertilizing species experience reduced effective population size, recombination rates, and heterozygosity, which in turn should decrease the efficacy of natural selection, both adaptive and purifying, and the strength of meiotic drive processes such as GC-biased gene conversion. The empirical evidence is only partly congruent with these predictions, depending on the analyzed species, some, but not all, of the expected effects have been observed. One possible reason is that self-fertilization is an evolutionary dead-end, so that most current selfers recently evolved self-fertilization, and their genome has not yet been strongly impacted by selfing. Here, we investigate the molecular evolution of two groups of freshwater snails in which mating systems have likely been stable for several millions of years. Analyzing coding sequence polymorphism, divergence, and expression levels, we report a strongly reduced genetic diversity, decreased efficacy of purifying selection, slower rate of adaptive evolution, and weakened codon usage bias/GC-biased gene conversion in the selfer Galba compared with the outcrosser Physa , in full agreement with theoretical expectations. Our results demonstrate that self-fertilization, when effective in the long run, is a major driver of population genomic and molecular evolutionary processes. Despite the genomic effects of selfing, Galba truncatula seems to escape the demographic consequences of the genetic load. We suggest that the particular ecology of the species may buffer the negative consequences of selfing, shedding new light on the dead-end hypothesis.

Description: Choosing the number of components in a finite mixture model is a challenging task. In this article, we study the behaviour of information criteria for selecting the mixture order, based on either the observed likelihood or the complete likelihood including component labels. We propose a new observed likelihood criterion called aic mix , which is shown to be order consistent. We further show that when there is a nontrivial level of classification uncertainty in the true model, complete likelihood criteria asymptotically underestimate the true number of components. A simulation study illustrates the potentially poor finite-sample performance of complete likelihood criteria, while aic mix and the Bayesian information criterion perform strongly regardless of the level of classification uncertainty.

Description: Outlier detection is an integral component of statistical modelling and estimation. For high-dimensional data, classical methods based on the Mahalanobis distance are usually not applicable. We propose an outlier detection procedure that replaces the classical minimum covariance determinant estimator with a high-breakdown minimum diagonal product estimator. The cut-off value is obtained from the asymptotic distribution of the distance, which enables us to control the Type I error and deliver robust outlier detection. Simulation studies show that the proposed method behaves well for high-dimensional data.

Description: This paper adopts a nonparametric Bayesian approach to testing whether a function is monotone. Two new families of tests are constructed. The first uses constrained smoothing splines with a hierarchical stochastic-process prior that explicitly controls the prior probability of monotonicity. The second uses regression splines together with two proposals for the prior over the regression coefficients. Via simulation, the finite-sample performance of the tests is shown to improve upon existing frequentist and Bayesian methods. The asymptotic properties of the Bayes factor for comparing monotone versus nonmonotone regression functions in a Gaussian model are also studied. Our results significantly extend those currently available, which chiefly focus on determining the dimension of a parametric linear model.

Description: The paper develops hierarchical empirical Bayes and benchmarked hierarchical empirical Bayes estimators of positive small area means under multiplicative models. The usual benchmarking requirement is that the small area estimates, when aggregated, should equal the direct estimates for the larger geographical areas. However, while estimating positive small area parameters, the conventional squared error or weighted squared error loss subject to the usual benchmark constraint may not produce positive estimators, so it is necessary to seek other loss functions. We consider a multiplicative model for the original data for estimating positive small area means, and suggest a variant of the Kullback–Leibler divergence as a loss function. The prediction errors of the suggested hierarchical empirical Bayes estimators are investigated asymptotically, and their second-order unbiased estimators are provided. Bootstrapped estimators of these prediction errors for both hierarchical empirical Bayes and benchmarked hierarchical empirical Bayes estimators are also given. The performance of the suggested procedures is investigated through simulation as well as with an example.

Description: Odds ratios can be estimated in case-control studies using standard logistic regression, ignoring the outcome-dependent sampling. In this paper we discuss an analogous result for treatment effects on the treated in matched cohort studies. Specifically, in studies where a sample of treated subjects is observed along with a separate sample of possibly matched controls, we show that efficient and doubly robust estimators of effects on the treated are computationally equivalent to standard estimators, which ignore the matching and exposure-based sampling. This is not the case for general average effects. We also show that matched cohort studies are often more efficient than random sampling for estimating effects on the treated, and derive the optimal number of matches for a given set of matching variables. We illustrate our results via simulation and in a matched cohort study of the effect of hysterectomy on the risk of cardiovascular disease.

Description: Current status data occur in contexts including demographic studies and tumorigenicity experiments. In such cases, each subject is observed only once and the failure time of interest is either left- or right-censored (Kalbfleisch & Prentice, 2002). Many methods have been developed for the analysis of such data (Huang, 1996; Sun, 2006), most of which assume that the failure time and the observation time are independent completely or given covariates. In this paper, we present a sieve maximum likelihood approach for current status data when independence does not hold. A copula model and monotone I-splines are used and the asymptotic properties of the resulting estimators are established. In particular, the estimated regression parameters are shown to be semiparametrically efficient. An illustrative example is provided.

Description: Characterizing Leaf Area Index (LAI) and Vertical Foliage Profile (VFP) over the United States Biogeosciences Discussions, 12, 13675-13710, 2015 Author(s): H. Tang, S. Ganguly, G. Zhang, M. A. Hofton, R. F. Nelson, and R. Dubayah Leaf area index (LAI) and vertical foliage profile (VFP) are among the important canopy structural variables. Recent advances in lidar remote sensing technology have demonstrated the capability of accurately mapping LAI and VFP over large areas. The primary objective of this study was to derive and validate a LAI and VFP product over the contiguous United States using spaceborne waveform lidar data. This product was derived at the footprint level from the Geoscience Laser Altimeter System (GLAS) using a biophysical model. We validated GLAS derived LAI and VFP across major forest biomes using airborne waveform lidar. The comparison results showed that GLAS retrievals of total LAI were generally accurate with little bias ( r 2 = 0.67, bias = −0.13, RMSE = 0.75). The derivations of GLAS retrievals of VFP within layers was not as accurate overall ( r 2 = 0.36, bias = −0.04, RMSE = 0.26), and these varied as a function of height, increasing from understory to overstory −0 to 5 m layer: r 2 = 0.04, bias = 0.09, RMSE = 0.31; 10 to 15 m layer: r 2 = 0.53, bias = −0.08, RMSE = 0.22; and 15 to 20 m layer: r 2 = 0.66, bias =−0.05, RMSE = 0.20. Significant relationships were also found between GLAS LAI products and different environmental factors, in particular elevation and annual precipitation. In summary, our results provide a unique insight into vertical canopy structure distribution across North American ecosystems. This data set is a first step towards a baseline of canopy structure needed for evaluating climate and land use induced forest changes at continental scale in the future and should help deepen our understanding of the role of vertical canopy structure on terrestrial ecosystem processes across varying scales.

Description: Motivation: In practice, identifying and interpreting the functional impacts of the regulatory relationships between micro-RNA and messenger-RNA is non-trivial. The sheer scale of possible micro-RNA and messenger-RNA interactions can make the interpretation of results difficult. Results: We propose a supervised framework, pMim, built upon concepts of significance combination, for jointly ranking regulatory micro-RNA and their potential functional impacts with respect to a condition of interest. Here, pMim directly tests if a micro-RNA is differentially expressed and if its predicted targets, which lie in a common biological pathway, have changed in the opposite direction. We leverage the information within existing micro-RNA target and pathway databases to stabilize the estimation and annotation of micro-RNA regulation making our approach suitable for datasets with small sample sizes. In addition to outputting meaningful and interpretable results, we demonstrate in a variety of datasets that the micro-RNA identified by pMim, in comparison to simpler existing approaches, are also more concordant with what is described in the literature. Availability and implementation: This framework is implemented as an R function, pMim , in the package sydSeq available from http://www.ellispatrick.com/r-packages . Contact: jean.yang@sydney.edu.au Supplementary information: Supplementary data are available at Bioinformatics online.