Transcriptome sequencing of Festulolium accessions under salt stress

peer-reviewedObjectives The objective of this study was to establish transcriptome assemblies of Festulolium hybrids under salt stress, and identify genes regulated across the hybrids in response to salt stress. The development of transcriptome assemblies for Festulolium hybrids and cataloguing of genes regulated under salt stress will facilitate further downstream studies. Results Plants were grown at three salt concentrations (0.5%, 1% and 1.5%) and phenotypic and transcriptomic data was collected. Salt stress was confirmed by progressive loss of green leaves as salt concentration increased from 0 to 1.5%. We generated de-novo transcriptome assemblies for two Festulolium pabulare festucoid genotypes, for a single Festulolium braunii genotype, and a single F. pabulare loloid genotype. We also identified 1555 transcripts that were up regulated and 1264 transcripts that were down regulated in response to salt stress in the Festulolium hybrids. Some of the identified transcripts showed significant sequence similarity with genes known to be regulated during salt and other abiotic stresses

BoostMe accurately predicts DNA methylation values in whole-genome bisulfite sequencing of multiple human tissues

Abstract Background Bisulfite sequencing is widely employed to study the role of DNA methylation in disease; however, the data suffer from biases due to coverage depth variability. Imputation of methylation values at low-coverage sites may mitigate these biases while also identifying important genomic features associated with predictive power. Results Here we describe BoostMe, a method for imputing low-quality DNA methylation estimates within whole-genome bisulfite sequencing (WGBS) data. BoostMe uses a gradient boosting algorithm, XGBoost, and leverages information from multiple samples for prediction. We find that BoostMe outperforms existing algorithms in speed and accuracy when applied to WGBS of human tissues. Furthermore, we show that imputation improves concordance between WGBS and the MethylationEPIC array at low WGBS depth, suggesting improved WGBS accuracy after imputation. Conclusions Our findings support the use of BoostMe as a preprocessing step for WGBS analysis.https://deepblue.lib.umich.edu/bitstream/2027.42/143848/1/12864_2018_Article_4766.pd

Single-cell ATAC-Seq in human pancreatic islets and deep learning upscaling of rare cells reveals cell-specific type 2 diabetes regulatory signatures

Objective: Type 2 diabetes (T2D) is a complex disease characterized by pancreatic islet dysfunction, insulin resistance, and disruption of blood glucose levels. Genome-wide association studies (GWAS) have identified > 400 independent signals that encode genetic predisposition. More than 90% of associated single-nucleotide polymorphisms (SNPs) localize to non-coding regions and are enriched in chromatin-defined islet enhancer elements, indicating a strong transcriptional regulatory component to disease susceptibility. Pancreatic islets are a mixture of cell types that express distinct hormonal programs, so each cell type may contribute differentially to the underlying regulatory processes that modulate T2D-associated transcriptional circuits. Existing chromatin profiling methods such as ATAC-seq and DNase-seq, applied to islets in bulk, produce aggregate profiles that mask important cellular and regulatory heterogeneity. Methods: We present genome-wide single-cell chromatin accessibility profiles in >1,600 cells derived from a human pancreatic islet sample using single-cell combinatorial indexing ATAC-seq (sci-ATAC-seq). We also developed a deep learning model based on U-Net architecture to accurately predict open chromatin peak calls in rare cell populations. Results: We show that sci-ATAC-seq profiles allow us to deconvolve alpha, beta, and delta cell populations and identify cell-type-specific regulatory signatures underlying T2D. Particularly, T2D GWAS SNPs are significantly enriched in beta cell-specific and across cell-type shared islet open chromatin, but not in alpha or delta cell-specific open chromatin. We also demonstrate, using less abundant delta cells, that deep learning models can improve signal recovery and feature reconstruction of rarer cell populations. Finally, we use co-accessibility measures to nominate the cell-specific target genes at 104 non-coding T2D GWAS signals. Conclusions: Collectively, we identify the islet cell type of action across genetic signals of T2D predisposition and provide higher-resolution mechanistic insights into genetically encoded risk pathways. Published by Elsevier GmbH.Open access journalThis item from the UA Faculty Publications collection is made available by the University of Arizona with support from the University of Arizona Libraries. If you have questions, please contact us at [email protected]

Modeling soil organic carbon dynamics in temperate forests with Yasso07

In a context of global changes, modeling and predicting the dynamics of soil carbon stocks (CSs) in forest ecosystems are vital but challenging. Yasso07 is considered to be one of the most promising models for such a purpose. We examine the accuracy of its prediction of soil carbon dynamics over the whole French metropolitan territory at a decennial timescale. We used data from 101 sites in the RENECOFOR network, which encompasses most of the French temperate forests. These data include (i) the quantity of above-ground litterfall from 1994 to 2008, measured yearly, and (ii) the soil CSs measured twice at an interval of approximately 15 years (once in the early 1990s and around 2010). We used Yasso07 to simulate the annual changes in carbon stocks (ACCs; in tC&thinsp;ha−1&thinsp;yr−1) for each site and then compared the estimates with actual recorded data. We carried out meta-analyses to reveal the variability in litter biochemistry in different tree organs for conifers and broadleaves. We also performed sensitivity analyses to explore Yasso07's sensitivity to annual litter inputs and model initialization settings. At the national level, the simulated ACCs (+0.00±0.07&thinsp;tC&thinsp;ha−1&thinsp;yr−1, mean&thinsp;±&thinsp;SE) were of the same order of magnitude as the observed ones (+0.34±0.06&thinsp;tC&thinsp;ha−1&thinsp;yr−1). However, the correlation between predicted and measured ACCs remained weak (R2&lt;0.1). There was significant overestimation for broadleaved stands and underestimation for coniferous sites. Sensitivity analyses showed that the final estimated CS was strongly affected by settings in the model initialization, including litter and soil carbon quantity and quality and also by simulation length. Carbon quality set with the partial steady-state assumption gave a better fit than the model with the complete steady-state assumption. With Yasso07 as the support model, we showed that there is currently a bottleneck in soil carbon modeling and prediction due to a lack of knowledge or data on soil carbon quality and fine-root quantity in the litter.</p

Integrative analysis of gene expression, DNA methylation, physiological traits, and genetic variation in human skeletal muscle

We integrate comeasured gene expression and DNA methylation (DNAme) in 265 human skeletal muscle biopsies from the FUSION study with >7 million genetic variants and eight physiological traits: height, waist, weight, waist-hip ratio, body mass index, fasting serum insulin, fasting plasma glucose, and type 2 diabetes. We find hundreds of genes and DNAme sites associated with fasting insulin, waist, and body mass index, as well as thousands of DNAme sites associated with gene expression (eQTM). We find that controlling for heterogeneity in tissue/muscle fiber type reduces the number of physiological trait associations, and that long-range eQTMs (>1 Mb) are reduced when controlling for tissue/muscle fiber type or latent factors. We map genetic regulators (quantitative trait loci; QTLs) of expression (eQTLs) and DNAme (mQTLs). Using Mendelian randomization (MR) and mediation techniques, we leverage these genetic maps to predict 213 causal relationships between expression and DNAme, approximately two-thirds of which predict methylation to causally influence expression. We use MR to integrate FUSION mQTLs, FUSION eQTLs, and GTEx eQTLs for 48 tissues with genetic associations for 534 diseases and quantitative traits. We identify hundreds of genes and thousands of DNAme sites that may drive the reported disease/quantitative trait genetic associations. We identify 300 gene expression MR associations that are present in both FUSION and GTEx skeletal muscle and that show stronger evidence of MR association in skeletal muscle than other tissues, which may partially reflect differences in power across tissues. As one example, we find that increased RXRA muscle expression may decrease lean tissue mass.Peer reviewe

Genetic regulatory signatures underlying islet gene expression and type 2 diabetes

The majority of genetic variants associated with type 2 diabetes (T2D) are located outside of genes in noncoding regions that may regulate gene expression in disease-relevant tissues, like pancreatic islets. Here, we present the largest integrated analysis to date of high-resolution, high-throughput human islet molecular profiling data to characterize the genome (DNA), epigenome (DNA packaging), and transcriptome (gene expression). We find that T2D genetic variants are enriched in regions of the genome where transcription Regulatory Factor X (RFX) is predicted to bind in an islet-specific manner. Genetic variants that increase T2D risk are predicted to disrupt RFX binding, providing a molecular mechanism to explain how the genome can influence the epigenome, modulating gene expression and ultimately T2D risk

Is There Such a Thing as Psychological Pain? and Why It Matters

Medicine regards pain as a signal of physical injury to the body despite evidence contradicting the linkage and despite the exclusion of vast numbers of sufferers who experience psychological pain. By broadening our concept of pain and making it more inclusive, we would not only better accommodate the basic science of pain but also would recognize what is already appreciated by the layperson—that pain from diverse sources, physical and psychological, share an underlying felt structure

Genetic and genomic analysis of hyperlipidemia, obesity and diabetes using (C57BL/6J × TALLYHO/JngJ) F2 mice

<p>Abstract</p> <p>Background</p> <p>Type 2 diabetes (T2D) is the most common form of diabetes in humans and is closely associated with dyslipidemia and obesity that magnifies the mortality and morbidity related to T2D. The genetic contribution to human T2D and related metabolic disorders is evident, and mostly follows polygenic inheritance. The TALLYHO/JngJ (TH) mice are a polygenic model for T2D characterized by obesity, hyperinsulinemia, impaired glucose uptake and tolerance, hyperlipidemia, and hyperglycemia.</p> <p>Results</p> <p>In order to determine the genetic factors that contribute to these T2D related characteristics in TH mice, we interbred TH mice with C57BL/6J (B6) mice. The parental, F1, and F2 mice were phenotyped at 8, 12, 16, 20, and 24 weeks of age for 4-hour fasting plasma triglyceride, cholesterol, insulin, and glucose levels and body, fat pad and carcass weights. The F2 mice were genotyped genome-wide and used for quantitative trait locus (QTL) mapping. We also applied a genetical genomic approach using a subset of the F2 mice to seek candidate genes underlying the QTLs. Major QTLs were detected on chromosomes (Chrs) 1, 11, 4, and 8 for hypertriglyceridemia, 1 and 3 for hypercholesterolemia, 4 for hyperglycemia, 11 and 1 for body weight, 1 for fat pad weight, and 11 and 14 for carcass weight. Most alleles, except for Chr 3 and 14 QTLs, increased phenotypic values when contributed by the TH strain. Fourteen pairs of interacting loci were detected, none of which overlapped the major QTLs. The QTL interval linked to hypercholesterolemia and hypertriglyceridemia on distal Chr 1 contains <it>Apoa2 </it>gene. Sequencing analysis revealed polymorphisms of <it>Apoa2 </it>in TH mice, suggesting <it>Apoa2 </it>as the candidate gene for the hyperlipidemia QTL. Gene expression analysis added novel information and aided in selection of candidates underlying the QTLs.</p> <p>Conclusions</p> <p>We identified several genetic loci that affect the quantitative variations of plasma lipid and glucose levels and obesity traits in a TH × B6 intercross. Polymorphisms in <it>Apoa2 </it>gene are suggested to be responsible for the Chr 1 QTL linked to hypercholesterolemia and hypertriglyceridemia. Further, genetical genomic analysis led to potential candidate genes for the QTLs.</p