Background: Gene bodies are the most evolutionarily conserved targets of DNA methylation in eukaryotes. However, the regulatory functions of gene body DNA methylation remain largely unknown. DNA methylation in insects appears to be primarily confined to exons. Two recent studies in Apis mellifera (honeybee) and Nasonia vitripennis (jewel wasp) analyzed transcription and DNA methylation data for one gene in each species to demonstrate that exon-specific DNA methylation may be associated with alternative splicing events. In this study we investigated the relationship between DNA methylation, alternative splicing, and cross-species gene conservation on a genome-wide scale using genome-wide transcription and DNA methylation data.

Results: We generated RNA deep sequencing data (RNA-seq) to measure genome-wide mRNA expression at the exon- and gene-level. We produced a de novo transcriptome from this RNA-seq data and computationally predicted splice variants for the honeybee genome. We found that exons that are included in transcription are higher methylated than exons that are skipped during transcription. We detected enrichment for alternative splicing among methylated genes compared to unmethylated genes using fisher’s exact test. We performed a statistical analysis to reveal that the presence of DNA methylation or alternative splicing are both factors associated with a longer gene length and a greater number of exons in genes. In concordance with this observation, a conservation analysis using BLAST revealed that each of these factors is also associated with higher cross-species gene conservation.

Conclusions: This study constitutes the first genome-wide analysis exhibiting a positive relationship between exon-level DNA methylation and mRNA expression in the honeybee. Our finding that methylated genes are enriched for alternative splicing suggests that, in invertebrates, exon-level DNA methylation may play a role in the construction of splice variants by positively influencing exon inclusion during transcription. The results from our cross-species homology analysis suggest that DNA methylation and alternative splicing are genetic mechanisms whose utilization could contribute to a longer gene length and a slower rate of gene evolution.

The sensitivity of Earth’s wetlands to observed shifts in global precipitation and temperature patterns and their ability to produce large quantities of methane gas are key global change questions. We present a microwave satellite-based approach for mapping fractional surface water (FW) globally at 25-km resolution. The approach employs a land cover-supported, atmospherically-corrected dynamic mixture model applied to 20+ years (1992–2013) of combined, daily, passive/active microwave remote sensing data. The resulting product, known as Surface Water Microwave Product Series (SWAMPS), shows strong microwave sensitivity to sub-grid scale open water and inundated wetlands comprising open plant canopies. SWAMPS’ FW compares favorably (R² = 91%–94%) with higher-resolution, global-scale maps of open water from MODIS and SRTM-MOD44W. Correspondence of SWAMPS with open water and wetland products from satellite SAR in Alaska and the Amazon deteriorates when exposed wetlands or inundated forests captured by the SAR products were added to the open water fraction reflecting SWAMPS’ inability to detect water underneath the soil surface or beneath closed forest canopies. Except for a brief period of drying during the first 4 years of observation, the inundation extent for the global domain excluding the coast was largely stable. Regionally, inundation in North America is advancing while inundation is on the retreat in Tropical Africa and North Eurasia. SWAMPS provides a consistent and long-term global record of daily FW dynamics, with documented accuracies suitable for hydrologic assessment and global change-related investigations.

Epigenetic changes enable genomes to respond to changes in the environment, such as altered nutrition, activity, or social setting. Epigenetic modifications, thereby, provide a source of phenotypic plasticity in many species. The honey bee (Apis mellifera) uses nutritionally sensitive epigenetic control mechanisms in the development of the royal caste (queens) and the workers. The workers are functionally sterile females that can take on a range of distinct physiological and/or behavioral phenotypes in response to environmental changes. Honey bees have a wide repertoire of epigenetic mechanisms which, as in mammals, include cytosine methylation, hydroxymethylated cytosines, together with the enzymatic machinery responsible for these cytosine modifications. Current data suggests that honey bees provide an excellent system for studying the “social repertoire” of the epigenome. In this review, we elucidate what is known so far about the honey bee epigenome and its mechanisms. Our discussion includes what may distinguish honey bees from other model animals, how the epigenome can influence worker behavioral task separation, and how future studies can answer central questions about the role of the epigenome in social behavior.

A warming climate is altering land-atmosphere exchanges of carbon, with a potential for increased vegetation productivity as well as the mobilization of permafrost soil carbon stores. Here we investigate land-atmosphere carbon dioxide (CO₂) cycling through analysis of net ecosystem productivity (NEP) and its component fluxes of gross primary productivity (GPP) and ecosystem respiration (ER) and soil carbon residence time, simulated by a set of land surface models (LSMs) over a region spanning the drainage basin of Northern Eurasia. The retrospective simulations cover the period 1960–2009 at 0.5° resolution, which is a scale common among many global carbon and climate model simulations. Model performance benchmarks were drawn from comparisons against both observed CO₂ fluxes derived from site-based eddy covariance measurements as well as regional-scale GPP estimates based on satellite remote-sensing data.

The site-based comparisons depict a tendency for overestimates in GPP and ER for several of the models, particularly at the two sites to the south. For several models the spatial pattern in GPP explains less than half the variance in the MODIS MOD17 GPP product. Across the models NEP increases by as little as 0.01 to as much as 0.79 g C m^-2 yr^-2, equivalent to 3 to 340 % of the respective model means, over the analysis period. For the multimodel average the increase is 135 % of the mean from the first to last 10 years of record (1960–1969 vs. 2000–2009), with a weakening CO₂ sink over the latter decades. Vegetation net primary productivity increased by 8 to 30 % from the first to last 10 years, contributing to soil carbon storage gains. The range in regional mean NEP among the group is twice the multimodel mean, indicative of the uncertainty in CO₂ sink strength.

The models simulate that inputs to the soil carbon pool exceeded losses, resulting in a net soil carbon gain amid a decrease in residence time. Our analysis points to improvements in model elements controlling vegetation productivity and soil respiration as being needed for reducing uncertainty in land-atmosphere CO₂ exchange. These advances will require collection of new field data on vegetation and soil dynamics, the development of benchmarking data sets from measurements and remote-sensing observations, and investments in future model development and intercomparison studies.

Soil temperature (T_s) change is a key indicator of the dynamics of permafrost. On seasonal and interannual timescales, the variability of T_s determines the active-layer depth, which regulates hydrological soil properties and biogeochemical processes. On the multi-decadal scale, increasing T_s not only drives permafrost thaw/retreat but can also trigger and accelerate the decomposition of soil organic carbon. The magnitude of permafrost carbon feedbacks is thus closely linked to the rate of change of soil thermal regimes. In this study, we used nine process-based ecosystem models with permafrost processes, all forced by different observation-based climate forcing during the period 1960–2000, to characterize the warming rate of T_s in permafrost regions. There is a large spread of T_s trends at 20 cm depth across the models, with trend values ranging from 0.010 ± 0.003 to 0.031 ± 0.005 °C yr^-1. Most models show smaller increase in T_s with increasing depth. Air temperature (Tsub>a) and longwave downward radiation (LWDR) are the main drivers of T_s trends, but their relative contributions differ amongst the models. Different trends of LWDR used in the forcing of models can explain 61 % of their differences in T_s trends, while trends of T_a only explain 5 % of the differences in T_s trends. Uncertain climate forcing contributes a larger uncertainty in T_s trends (0.021 ± 0.008 °C yr^-1, mean ± standard deviation) than the uncertainty of model structure (0.012 ± 0.001 °C yr^-1), diagnosed from the range of response between different models, normalized to the same forcing. In addition, the loss rate of near-surface permafrost area, defined as total area where the maximum seasonal active-layer thickness (ALT) is less than 3 m loss rate, is found to be significantly correlated with the magnitude of the trends of T_s at 1 m depth across the models (R = −0.85, P = 0.003), but not with the initial total near-surface permafrost area (R = −0.30, P = 0.438). The sensitivity of the total boreal near-surface permafrost area to T_s at 1 m is estimated to be of −2.80 ± 0.67 million km²°C^-1. Finally, by using two long-term LWDR data sets and relationships between trends of LWDR and T_s across models, we infer an observation-constrained total boreal near-surface permafrost area decrease comprising between 39 ± 14  ×  10³ and 75 ± 14  ×  10³km²yr^-1 from 1960 to 2000. This corresponds to 9–18 % degradation of the current permafrost area.

A realistic simulation of snow cover and its thermal properties are important for accurate modelling of permafrost. We analyse simulated relationships between air and near-surface (20  cm) soil temperatures in the Northern Hemisphere permafrost region during winter, with a particular focus on snow insulation effects in nine land surface models, and compare them with observations from 268 Russian stations. There are large cross-model differences in the simulated differences between near-surface soil and air temperatures (ΔT; 3 to 14 °C), in the sensitivity of soil-to-air temperature (0.13 to 0.96 °C °^C-1), and in the relationship between ΔT and snow depth. The observed relationship between ΔT and snow depth can be used as a metric to evaluate the effects of each model's representation of snow insulation, hence guide improvements to the model's conceptual structure and process parameterisations. Models with better performance apply multilayer snow schemes and consider complex snow processes. Some models show poor performance in representing snow insulation due to underestimation of snow depth and/or overestimation of snow conductivity. Generally, models identified as most acceptable with respect to snow insulation simulate reasonable areas of near-surface permafrost (13.19 to 15.77 million  km²). However, there is not a simple relationship between the sophistication of the snow insulation in the acceptable models and the simulated area of Northern Hemisphere near-surface permafrost, because several other factors, such as soil depth used in the models, the treatment of soil organic matter content, hydrology and vegetation cover, also affect the simulated permafrost distribution.

Honey bees as other insects rely on the innate immune system for protection against diseases. The innate immune system includes the circulating hemocytes (immune cells) that clear pathogens from hemolymph (blood) by phagocytosis, nodulation or encapsulation. Honey bee hemocyte numbers have been linked to hemolymph levels of vitellogenin. Vitellogenin is a multifunctional protein with immune-supportive functions identified in a range of species, including the honey bee. Hemocyte numbers can increase via mitosis, and this recruitment process can be important for immune system function and maintenance. Here, we tested if hemocyte mediated phagocytosis differs among the physiologically different honey bee worker castes (nurses, foragers and winter bees), and study possible interactions with vitellogenin and hemocyte recruitment. To this end, we adapted phagocytosis assays, which—together with confocal microscopy and flow cytometry—allow qualitative and quantitative assessment of hemocyte performance. We found that nurses are more efficient in phagocytic uptake than both foragers and winter bees. We detected vitellogenin within the hemocytes, and found that winter bees have the highest numbers of vitellogenin-positive hemocytes. Connections between phagocytosis, hemocyte-vitellogenin and mitosis were worker caste dependent. Our results demonstrate that the phagocytic performance of immune cells differs significantly between honey bee worker castes, and support increased immune competence in nurses as compared to forager bees. Our data, moreover, provides support for roles of vitellogenin in hemocyte activity.

Honey bees (Apis mellifera) provide a system for studying social and food-related behavior. A caste of workers performs age-related tasks: young bees (nurses) usually feed the brood and other adult bees inside the nest, while older bees (foragers) forage outside for pollen, a protein/lipid source, or nectar, a carbohydrate source. The workers' transition from nursing to foraging and their foraging preferences correlate with differences in gustatory perception, metabolic gene expression, and endocrine physiology including the endocrine factors vitellogenin (Vg) and juvenile hormone (JH). However, the understanding of connections among social behavior, energy metabolism, and endocrine factors is incomplete. We used RNA interference (RNAi) to perturb the gene network of Vg and JH to learn more about these connections through effects on gustation, gene transcripts, and physiology. The RNAi perturbation was achieved by single and double knockdown of the genes ultraspiracle (usp) and vg, which encode a putative JH receptor and Vg, respectively. The double knockdown enhanced gustatory perception and elevated hemolymph glucose, trehalose, and JH. We also observed transcriptional responses in insulin like peptide 1 (ilp1), the adipokinetic hormone receptor (AKHR), and cGMP-dependent protein kinase (PKG, or “foraging gene” Amfor). Our study demonstrates that the Vg–JH regulatory module controls changes in carbohydrate metabolism, but not lipid metabolism, when worker bees shift from nursing to foraging. The module is also placed upstream of ilp1, AKHR, and PKG for the first time. As insulin, adipokinetic hormone (AKH), and PKG pathways influence metabolism and gustation in many animals, we propose that honey bees have conserved pathways in carbohydrate metabolism and conserved connections between energy metabolism and gustatory perception. Thus, perhaps the bee can make general contributions to the understanding of food-related behavior and metabolic disorders.

Insect immune systems can recognize specific pathogens and prime offspring immunity. High specificity of immune priming can be achieved when insect females transfer immune elicitors into developing oocytes. The molecular mechanism behind this transfer has been a mystery. Here, we establish that the egg-yolk protein vitellogenin is the carrier of immune elicitors. Using the honey bee, Apis mellifera, model system, we demonstrate with microscopy and western blotting that vitellogenin binds to bacteria, both Paenibacillus larvae – the gram-positive bacterium causing American foulbrood disease – and to Escherichia coli that represents gram-negative bacteria. Next, we verify that vitellogenin binds to pathogen-associated molecular patterns; lipopolysaccharide, peptidoglycan and zymosan, using surface plasmon resonance. We document that vitellogenin is required for transport of cell-wall pieces of E. coli into eggs by imaging tissue sections. These experiments identify vitellogenin, which is distributed widely in oviparous species, as the carrier of immune-priming signals. This work reveals a molecular explanation for trans-generational immunity in insects and a previously undescribed role for vitellogenin.

A data set of observed daily precipitation, maximum and minimum temperature, gridded to a 1/16° (~6 km) resolution, is described that spans the entire country of Mexico, the conterminous U.S. (CONUS), and regions of Canada south of 53° N for the period 1950-2013. The dataset improves previous products in spatial extent, orographic precipitation adjustment over Mexico and parts of Canada, and reduction of transboundary discontinuities. The impacts of adjusting gridded precipitation for orographic effects are quantified by scaling precipitation to an elevation-aware 1981-2010 precipitation climatology in Mexico and Canada. Differences are evaluated in terms of total precipitation as well as by hydrologic quantities simulated with a land surface model. Overall, orographic correction impacts total precipitation by up to 50% in mountainous regions outside CONUS. Hydrologic fluxes show sensitivities of similar magnitude, with discharge more sensitive than evapotranspiration and soil moisture. Because of the consistent gridding methodology, the current product reduces transboundary discontinuities as compared with a commonly used reanalysis product, making it suitable for estimating large-scale hydrometeorologic phenomena.