ASU Scholarship Showcase

Background: Meiotic recombination has traditionally been explained based on the structural requirement to stabilize homologous chromosome pairs to ensure their proper meiotic segregation. Competing hypotheses seek to explain the emerging findings of significant heterogeneity in recombination rates within and between genomes, but intraspecific comparisons of genome-wide recombination patterns are rare. The honey bee (Apis mellifera) exhibits the highest rate of genomic recombination among multicellular animals with about five cross-over events per chromatid.

Results: Here, we present a comparative analysis of recombination rates across eight genetic linkage maps of the honey bee genome to investigate which genomic sequence features are correlated with recombination rate and with its variation across the eight data sets, ranging in average marker spacing ranging from 1 Mbp to 120 kbp. Overall, we found that GC content explained best the variation in local recombination rate along chromosomes at the analyzed 100 kbp scale. In contrast, variation among the different maps was correlated to the abundance of microsatellites and several specific tri- and tetra-nucleotides.

Conclusions: The combined evidence from eight medium-scale recombination maps of the honey bee genome suggests that recombination rate variation in this highly recombining genome might be due to the DNA configuration instead of distinct sequence motifs. However, more fine-scale analyses are needed. The empirical basis of eight differing genetic maps allowed for robust conclusions about the correlates of the local recombination rates and enabled the study of the relation between DNA features and variability in local recombination rates, which is particularly relevant in the honey bee genome with its exceptionally high recombination rate.

Background: The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats.

Results: We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits.

Conclusions: These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

Honeybee workers are essentially sterile female helpers that make up the majority of individuals in a colony. Workers display a marked change in physiology when they transition from in-nest tasks to foraging. Recent technological advances have made it possible to unravel the metabolic modifications associated with this transition. Previous studies have revealed extensive remodeling of brain, thorax, and hypopharyngeal gland biochemistry. However, data on changes in the abdomen is scarce. To narrow this gap we investigated the proteomic composition of abdominal tissue in the days typically preceding the onset of foraging in honeybee workers.

In order to get a broader representation of possible protein dynamics, we used workers of two genotypes with differences in the age at which they initiate foraging. This approach was combined with RNA interference-mediated downregulation of an insulin/insulin-like signaling component that is central to foraging behavior, the insulin receptor substrate (irs), and with measurements of glucose and lipid levels.
Our data provide new insight into the molecular underpinnings of phenotypic plasticity in the honeybee, invoke parallels with vertebrate metabolism, and support an integrated and irs-dependent association of carbohydrate and lipid metabolism with the transition from in-nest tasks to foraging.

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.

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.

Background: The discovery of genetic associations is an important factor in the understanding of human illness to derive disease pathways. Identifying multiple interacting genetic mutations associated with disease remains challenging in studying the etiology of complex diseases. And although recently new single nucleotide polymorphisms (SNPs) at genes implicated in immune response, cholesterol/lipid metabolism, and cell membrane processes have been confirmed by genome-wide association studies (GWAS) to be associated with late-onset Alzheimer's disease (LOAD), a percentage of AD heritability continues to be unexplained. We try to find other genetic variants that may influence LOAD risk utilizing data mining methods.

Methods: Two different approaches were devised to select SNPs associated with LOAD in a publicly available GWAS data set consisting of three cohorts. In both approaches, single-locus analysis (logistic regression) was conducted to filter the data with a less conservative p-value than the Bonferroni threshold; this resulted in a subset of SNPs used next in multi-locus analysis (random forest (RF)). In the second approach, we took into account prior biological knowledge, and performed sample stratification and linkage disequilibrium (LD) in addition to logistic regression analysis to preselect loci to input into the RF classifier construction step.

Results: The first approach gave 199 SNPs mostly associated with genes in calcium signaling, cell adhesion, endocytosis, immune response, and synaptic function. These SNPs together with APOE and GAB2 SNPs formed a predictive subset for LOAD status with an average error of 9.8% using 10-fold cross validation (CV) in RF modeling. Nineteen variants in LD with ST5, TRPC1, ATG10, ANO3, NDUFA12, and NISCH respectively, genes linked directly or indirectly with neurobiology, were identified with the second approach. These variants were part of a model that included APOE and GAB2 SNPs to predict LOAD risk which produced a 10-fold CV average error of 17.5% in the classification modeling.

Conclusions: With the two proposed approaches, we identified a large subset of SNPs in genes mostly clustered around specific pathways/functions and a smaller set of SNPs, within or in proximity to five genes not previously reported, that may be relevant for the prediction/understanding of AD.

Background: Glioblastoma is the most aggressive primary central nervous tumor and carries a very poor prognosis. Invasion precludes effective treatment and virtually assures tumor recurrence. In the current study, we applied analytical and bioinformatics approaches to identify a set of microRNAs (miRs) from several different human glioblastoma cell lines that exhibit significant differential expression between migratory (edge) and migration-restricted (core) cell populations. The hypothesis of the study is that differential expression of miRs provides an epigenetic mechanism to drive cell migration and invasion.

Results: Our research data comprise gene expression values for a set of 805 human miRs collected from matched pairs of migratory and migration-restricted cell populations from seven different glioblastoma cell lines. We identified 62 down-regulated and 2 up-regulated miRs that exhibit significant differential expression in the migratory (edge) cell population compared to matched migration-restricted (core) cells. We then conducted target prediction and pathway enrichment analysis with these miRs to investigate potential associated gene and pathway targets. Several miRs in the list appear to directly target apoptosis related genes. The analysis identifies a set of genes that are predicted by 3 different algorithms, further emphasizing the potential validity of these miRs to promote glioblastoma.

Conclusions: The results of this study identify a set of miRs with potential for decreased expression in invasive glioblastoma cells. The verification of these miRs and their associated targeted proteins provides new insights for further investigation into therapeutic interventions. The methodological approaches employed here could be applied to the study of other diseases to provide biomedical researchers and clinicians with increased opportunities for therapeutic interventions.

Urban environmental measurements and observational statistics should reflect the properties generated over an adjacent area of adequate length where homogeneity is usually assumed. The determination of this characteristic source area that gives sufficient representation of the horizontal coverage of a sensing instrument or the fetch of transported quantities is of critical importance to guide the design and implementation of urban landscape planning strategies. In this study, we aim to unify two different methods for estimating source areas, viz. the statistical correlation method commonly used by geographers for landscape fragmentation and the mechanistic footprint model by meteorologists for atmospheric measurements. Good agreement was found in the intercomparison of the estimate of source areas by the two methods, based on 2-m air temperature measurement collected using a network of weather stations. The results can be extended to shed new lights on urban planning strategies, such as the use of urban vegetation for heat mitigation. In general, a sizable patch of landscape is required in order to play an effective role in regulating the local environment, proportional to the height at which stakeholders’ interest is mainly concerned.

Nocturnal cooling of urban areas governs the evolution of thermal state and many thermal-driven environmental issues in cities, especially those suffer strong urban heat island (UHI) effect. Advances in the fundamental understanding of the underlying physics of nighttime UHI involve disentangling complex contributing effects and remains an open challenge. In this study, we develop new numerical algorithms to characterize the thermodynamics of urban nocturnal cooling based on solving the energy balance equations for both the landscape surface and the overlying atmosphere. Further, a scaling law is proposed to relate the UHI intensity to a range of governing mechanisms, including the vertical and horizontal transport of heat in the surface layer, the urban-rural breeze, and the possible urban expansion. The accuracy of proposed methods is evaluated against in-situ urban measurements collected in cities with different geographic and climatic conditions. It is found that the vertical and horizontal contributors modulate the nocturnal UHI at distinct elevation in the atmospheric boundary layer.

Background: Phosphatase and TENsin (PTEN) homolog is a negative regulator that takes part in IIS (insulin/insulin-like signaling) and Egfr (epidermal growth factor receptor) activation in Drosophila melanogaster. IIS and Egfr signaling events are also involved in the developmental process of queen and worker differentiation in honey bees (Apis mellifera). Here, we characterized the bee PTEN gene homologue for the first time and begin to explore its potential function during bee development and adult life.

Results: Honey bee PTEN is alternatively spliced, resulting in three splice variants. Next, we show that the expression of PTEN can be down-regulated by RNA interference (RNAi) in the larval stage, when female caste fate is determined. Relative to controls, we observed that RNAi efficacy is dependent on the amount of PTEN dsRNA that is delivered to larvae. For larvae fed queen or worker diets containing a high amount of PTEN dsRNA, PTEN knockdown was significant at a whole-body level but lethal. A lower dosage did not result in a significant gene down-regulation. Finally, we compared same-aged adult workers with different behavior: nursing vs. foraging. We show that between nurses and foragers, PTEN isoforms were differentially expressed within brain, ovary and fat body tissues. All isoforms were expressed at higher levels in the brain and ovaries of the foragers. In fat body, isoform B was expressed at higher level in the nurse bees.

Conclusion: Our results suggest that PTEN plays a central role during growth and development in queen- and worker-destined honey bees. In adult workers, moreover, tissue-specific patterns of PTEN isoform expression are correlated with differences in complex division of labor between same-aged individuals. Therefore, we propose that knowledge on the roles of IIS and Egfr activity in developmental and behavioral control may increase through studies of how PTEN functions can impact bee social phenotypes.