Imputing Gene Expression in Uncollected Tissues Within and Beyond GTEx

Gene expression and its regulation can vary substantially across tissue types. In order to generate knowledge about gene expression in human tissues, the Genotype-Tissue Expression (GTEx) program has collected transcriptome data in a wide variety of tissue types from post-mortem donors. However, many tissue types are difficult to access and are not collected in every GTEx individual. Furthermore, in non-GTEx studies, the accessibility of certain tissue types greatly limits the feasibility and scale of studies of multi-tissue expression. In this work, we developed multi-tissue imputation methods to impute gene expression in uncollected or inaccessible tissues. Via simulation studies, we showed that the proposed methods outperform existing imputation methods in multi-tissue expression imputation and that incorporating imputed expression data can improve power to detect phenotype-expression correlations. By analyzing data from nine selected tissue types in the GTEx pilot project, we demonstrated that harnessing expression quantitative trait loci (eQTLs) and tissue-tissue expression-level correlations can aid imputation of transcriptome data from uncollected GTEx tissues. More importantly, we showed that by using GTEx data as a reference, one can impute expression levels in inaccessible tissues in non-GTEx expression studies

Genes Involved in the Evolution of Herbivory by a Leaf-Mining, Drosophilid Fly

Herbivorous insects are among the most successful radiations of life. However, we know little about the processes underpinning the evolution of herbivory. We examined the evolution of herbivory in the fly, Scaptomyza flava, whose larvae are leaf miners on species of Brassicaceae, including the widely studied reference plant, Arabidopsis thaliana (Arabidopsis). Scaptomyza flava is phylogenetically nested within the paraphyletic genus Drosophila, and the whole genome sequences available for 12 species of Drosophila facilitated phylogenetic analysis and assembly of a transcriptome for S. flava. A time-calibrated phylogeny indicated that leaf mining in Scaptomyza evolved between 6 and 16 million years ago. Feeding assays showed that biosynthesis of glucosinolates, the major class of antiherbivore chemical defense compounds in mustard leaves, was upregulated by S. flava larval feeding. The presence of glucosinolates in wild-type (WT) Arabidopsis plants reduced S. flava larval weight gain and increased egg–adult development time relative to flies reared in glucosinolate knockout (GKO) plants. An analysis of gene expression differences in 5-day-old larvae reared on WT versus GKO plants showed a total of 341 transcripts that were differentially regulated by glucosinolate uptake in larval S. flava. Of these, approximately a third corresponded to homologs of Drosophila melanogaster genes associated with starvation, dietary toxin-, heat-, oxidation-, and aging-related stress. The upregulated transcripts exhibited elevated rates of protein evolution compared with unregulated transcripts. The remaining differentially regulated transcripts also contained a higher proportion of novel genes than the unregulated transcripts. Thus, the transition to herbivory in Scaptomyza appears to be coupled with the evolution of novel genes and the co-option of conserved stress-related genes.Organismic and Evolutionary Biolog

Physical versus biological spring transition: 2005

In 2005, the onset of spring conditions in the physics of the coastal ocean (lowered sea level, spin-up of vertically-sheared equatorward coastal jet) came about 50 days later than average off Newport Oregon, on May 24. There was a further delay of 50 days before the subsurface upwelled water penetrated into the anomalously stratified surface layer, becoming most available for biological activity. The warm anomaly in sea surface temperature which provided the surface cap was observed at mid-shelf locations from Washington to central California, but it ended sooner south of Oregon. Biological impacts of these delays to several trophic levels have been reported

Systematic Comparison of Three Methods for Fragmentation of Long-Range PCR Products for Next Generation Sequencing

Next Generation Sequencing (NGS) technologies are gaining importance in the routine clinical diagnostic setting. It is thus desirable to simplify the workflow for high-throughput diagnostics. Fragmentation of DNA is a crucial step for preparation of template libraries and various methods are currently known. Here we evaluated the performance of nebulization, sonication and random enzymatic digestion of long-range PCR products on the results of NGS. All three methods produced high-quality sequencing libraries for the 454 platform. However, if long-range PCR products of different length were pooled equimolarly, sequence coverage drastically dropped for fragments below 3,000 bp. All three methods performed equally well with regard to overall sequence quality (PHRED) and read length. Enzymatic fragmentation showed highest consistency between three library preparations but performed slightly worse than sonication and nebulization with regard to insertions/deletions in the raw sequence reads. After filtering for homopolymer errors, enzymatic fragmentation performed best if compared to the results of classic Sanger sequencing. As the overall performance of all three methods was equal with only minor differences, a fragmentation method can be chosen solely according to lab facilities, feasibility and experimental design

Supersymmetric Scenarios with Dominant Radiative Neutralino Decay

The radiative decay of the next-to-lightest neutralino into a lightest neutralino and a photon is analyzed in the MSSM. We find that significant regions of the supersymmetric parameter space with large radiative BR's (up to about 100%) do exist. The radiative channel turns out to be enhanced when the neutralino tree-level decays are suppressed either "kinematically" or "dynamically". In general, in the regions allowed by LEP data and not characterized by asymptotic values of the SuSy parameters, the radiative enhancement requires tan beta ~= 1 and/or M_1 ~= M_2, and negative values of \mu. We present typical specific scenarios where these "necessary" conditions are fulfilled, relaxing the usual relation M_1=(5/3)*tan^2(th_W)*M_2. The influence of varying the stop masses and mixing angle when the radiative decay is enhanced is also considered. Some phenomenological consequences of the above picture are discussed.Comment: 32 pages, LaTeX file + 23 figures embedded with epsf.sty. In this revised version, Eq.(3) plus some related notations and text passages have been changed. Minor error corrected in Fig.12(a). The numerical analysis and the conclusions of the paper are not affected. (Includes the erratum to appear in Phys. Rev. D.) Source and ps files are also available at ftp://hpteo.roma1.infn.it/pub/preprints/ambr-mele/Rome1-1148/ or at http://feynman.physics.lsa.umich.edu/~ambros/Physics.html#1

Association of Accelerometry-Measured Physical Activity and Cardiovascular Events in Mobility-Limited Older Adults: The LIFE (Lifestyle Interventions and Independence for Elders) Study.

BACKGROUND:Data are sparse regarding the value of physical activity (PA) surveillance among older adults-particularly among those with mobility limitations. The objective of this study was to examine longitudinal associations between objectively measured daily PA and the incidence of cardiovascular events among older adults in the LIFE (Lifestyle Interventions and Independence for Elders) study. METHODS AND RESULTS:Cardiovascular events were adjudicated based on medical records review, and cardiovascular risk factors were controlled for in the analysis. Home-based activity data were collected by hip-worn accelerometers at baseline and at 6, 12, and 24 months postrandomization to either a physical activity or health education intervention. LIFE study participants (n=1590; age 78.9±5.2 [SD] years; 67.2% women) at baseline had an 11% lower incidence of experiencing a subsequent cardiovascular event per 500 steps taken per day based on activity data (hazard ratio, 0.89; 95% confidence interval, 0.84-0.96; P=0.001). At baseline, every 30 minutes spent performing activities ≥500 counts per minute (hazard ratio, 0.75; confidence interval, 0.65-0.89 [P=0.001]) were also associated with a lower incidence of cardiovascular events. Throughout follow-up (6, 12, and 24 months), both the number of steps per day (per 500 steps; hazard ratio, 0.90, confidence interval, 0.85-0.96 [P=0.001]) and duration of activity ≥500 counts per minute (per 30 minutes; hazard ratio, 0.76; confidence interval, 0.63-0.90 [P=0.002]) were significantly associated with lower cardiovascular event rates. CONCLUSIONS:Objective measurements of physical activity via accelerometry were associated with cardiovascular events among older adults with limited mobility (summary score >10 on the Short Physical Performance Battery) both using baseline and longitudinal data. CLINICAL TRIAL REGISTRATION:URL: http://www.clinicaltrials.gov. Unique identifier: NCT01072500

Reproductive Factors and Serum Uric Acid Levels in Females from the General Population: The KORA F4 Study

Hyperuricemia is associated with an increased risk of metabolic and cardiovascular diseases. There are pronounced sex differences in the levels of uric acid. It is largely unknown whether or not reproductive parameters which induce hormonal changes are responsible for this. We examined if there are associations between reproductive parameters and uric acid levels in a female population-based sample. In this cross-sectional analysis, data of 1530 women aged 32 to 81 years participating in the KORA F4 study, conducted between 2006 and 2008 in Southern Germany were used. Reproductive parameters were obtained by standardized interviews. Uric acid levels were tested by the uricase method. The whole study sample and stratified in pre- and postmenopausal women was analyzed. Menopausal status and earlier age at menarche were associated with higher serum uric acid levels (age-adjusted: p-values 0.003, <0.001 respectively; after multivariable adjustment, including BMI: p-values 0.002, 0.036). A history of oral contraceptive use showed an association with uric acid levels only after multivariable adjustment (p-value 0.009). Hot flushes showed an association with uric acid levels only after age-adjustment (p-value 0.038), but lost significance after adding other confounders. Other reproductive factors, including parity, current or ever use of hormone replacement therapy, current use of oral contraceptives, hysterectomy, bilateral oophorectomy, or depressive mood related to menopausal transition were not associated with uric acid levels. Postmenopausal status, earlier age at menarche and a history of oral contraceptive use were independently associated with higher serum uric acid concentrations in women from the general population. Further studies, especially longitudinal population-based studies investigating the relationship of female reproductive parameters with uric acid levels are necessary to confirm our findings

The impact of sex on gene expression across human tissues

Many complex human phenotypes exhibit sex-differentiated characteristics, however the underlying molecular mechanisms of these differences remain largely unknown. Here, we present an extensive catalog of both sex differences in gene expression and its genetic regulation across 44 human tissue sources surveyed by GTEx (v8 release). We demonstrate that sex strongly influences gene expression levels and cellular composition of tissue samples across the human body. The effect of sex on gene expression is widespread, with a total of 37% of all genes exhibiting sex-biased expression in at least one tissue. This suggests that many if not most biological processes, and thus complex traits and diseases, are impacted by sex effects on the transciptome. We expand the identification of cis-eQTLs with sex-differentiated effects and characterize their cellular origin. By integrating sex-biased eQTLs with genome-wide association study data, we identify 58 gene-trait associations that are driven by genetic regulation in a single sex, including novel associations not detected with sex-agnostic approaches. Altogether we provide the most comprehensive characterization of sex differences in the human transcriptome and its regulation to date.Peer ReviewedPreprin

Initial Evaluation of the Effects of Aerosolized Florida Red Tide Toxins (Brevetoxins) in Persons with Asthma

Florida red tides annually occur in the Gulf of Mexico, resulting from blooms of the marine dinoflagellate Karenia brevis. K. brevis produces highly potent natural polyether toxins, known as brevetoxins, that activate voltage-sensitive sodium channels. In experimental animals, brevetoxins cause significant bronchoconstriction. A study of persons who visited the beach recreationally found a significant increase in self-reported respiratory symptoms after exposure to aerosolized Florida red tides. Anecdotal reports indicate that persons with underlying respiratory diseases may be particularly susceptible to adverse health effects from these aerosolized toxins. Fifty-nine persons with physician-diagnosed asthma were evaluated for 1 hr before and after going to the beach on days with and without Florida red tide. Study participants were evaluated with a brief symptom questionnaire, nose and throat swabs, and spirometry approved by the National Institute for Occupational Safety and Health. Environmental monitoring, water and air sampling (i.e., K. brevis, brevetoxins, and particulate size distribution), and personal monitoring (for toxins) were performed. Brevetoxin concentrations were measured by liquid chromatography mass spectrometry, high-performance liquid chromatography, and a newly developed brevetoxin enzyme-linked immunosorbent assay. Participants were significantly more likely to report respiratory symptoms after Florida red tide exposure. Participants demonstrated small but statistically significant decreases in forced expiratory volume in 1 sec, forced expiratory flow between 25 and 75%, and peak expiratory flow after exposure, particularly those regularly using asthma medications. Similar evaluation during nonexposure periods did not significantly differ. This is the first study to show objectively measurable adverse health effects from exposure to aerosolized Florida red tide toxins in persons with asthma. Future studies will examine the possible chronic effects of these toxins among persons with asthma and other chronic respiratory impairment

A pervasive role for biomass burning in tropical high ozone/low water structures.

Air parcels with mixing ratios of high O3 and low H2O (HOLW) are common features in the tropical western Pacific (TWP) mid-troposphere (300-700 hPa). Here, using data collected during aircraft sampling of the TWP in winter 2014, we find strong, positive correlations of O3 with multiple biomass burning tracers in these HOLW structures. Ozone levels in these structures are about a factor of three larger than background. Models, satellite data and aircraft observations are used to show fires in tropical Africa and Southeast Asia are the dominant source of high O3 and that low H2O results from large-scale descent within the tropical troposphere. Previous explanations that attribute HOLW structures to transport from the stratosphere or mid-latitude troposphere are inconsistent with our observations. This study suggest a larger role for biomass burning in the radiative forcing of climate in the remote TWP than is commonly appreciated.We thank L. Pan for coordinating the CONTRAST flights and her constructive criticism of an early version of the manuscript; S. Schauffler, V. Donets and R. Lueb for collecting and analysing AWAS samples; T. Robinson and O. Shieh for providing meteorology forecasts in the field; and the pilots and crews of the CAST BAe-146 and CONTRAST Gulfstream V aircrafts for their dedication and professionalism. CAST was funded by the Natural Environment Research Council; CONTRAST was funded by the National Science Foundation. Research at the Jet Propulsion Laboratory, California Institute of Technology, is performed under contract with the National Aeronautics and Space Administration (NASA). A number of the US-based investigators also benefitted from the support of NASA as well as the National Oceanic and Atmospheric Administration. The views, opinions, and findings contained in this report are those of the author(s) and should not be construed as an official National Oceanic and Atmospheric Administration or US Government position, policy or decision. We would like to acknowledge high-performance computing support from Yellowstone (ark:/85065/d7wd3xhc) provided by NCAR's Computational and Information Systems Laboratory. NCAR is sponsored by the National Science Foundation.This is the final version of the article. It first appeared from Nature Publishing Group via http://dx.doi.org/10.1038/ncomms1026