The association of cigarette smoking with DNA methylation and gene expression in human tissue samples.

Cigarette smoking adversely affects many aspects of human health, and epigenetic responses to smoking may reflect mechanisms that mediate or defend against these effects. Prior studies of smoking and DNA methylation (DNAm), typically measured in leukocytes, have identified numerous smoking-associated regions (e.g., AHRR). To identify smoking-associated DNAm features in typically inaccessible tissues, we generated array-based DNAm data for 916 tissue samples from the GTEx (Genotype-Tissue Expression) project representing 9 tissue types (lung, colon, ovary, prostate, blood, breast, testis, kidney, and muscle). We identified 6,350 smoking-associated CpGs in lung tissue (n = 212) and 2,735 in colon tissue (n = 210), most not reported previously. For all 7 other tissue types (sample sizes 38-153), no clear associations were observed (false discovery rate 0.05), but some tissues showed enrichment for smoking-associated CpGs reported previously. For 1,646 loci (in lung) and 22 (in colon), smoking was associated with both DNAm and local gene expression. For loci detected in both lung and colon (e.g., AHRR, CYP1B1, CYP1A1), top CpGs often differed between tissues, but similar clusters of hyper- or hypomethylated CpGs were observed, with hypomethylation at regulatory elements corresponding to increased expression. For lung tissue, 17 hallmark gene sets were enriched for smoking-associated CpGs, including xenobiotic- and cancer-related gene sets. At least four smoking-associated regions in lung were impacted by lung methylation quantitative trait loci (QTLs) that co-localize with genome-wide association study (GWAS) signals for lung function (FEV1/FVC), suggesting epigenetic alterations can mediate the effects of smoking on lung health. Our multi-tissue approach has identified smoking-associated regions in disease-relevant tissues, including effects that are shared across tissue types.

Sequencing-based fine-mapping and in silico functional characterization of the 10q24.32 arsenic metabolism efficiency locus across multiple arsenic-exposed populations.

Inorganic arsenic is highly toxic and carcinogenic to humans. Exposed individuals vary in their ability to metabolize arsenic, and variability in arsenic metabolism efficiency (AME) is associated with risks of arsenic-related toxicities. Inherited genetic variation in the 10q24.32 region, near the arsenic methyltransferase (AS3MT) gene, is associated with urine-based measures of AME in multiple arsenic-exposed populations. To identify potential causal variants in this region, we applied fine mapping approaches to targeted sequencing data generated for exposed individuals from Bangladeshi, American Indian, and European American populations (n = 2,357, 557, and 648 respectively). We identified three independent association signals for Bangladeshis, two for American Indians, and one for European Americans. The size of the confidence sets for each signal varied from 4 to 85 variants. There was one signal shared across all three populations, represented by the same SNP in American Indians and European Americans (rs191177668) and in strong linkage disequilibrium (LD) with a lead SNP in Bangladesh (rs145537350). Beyond this shared signal, differences in LD patterns, minor allele frequency (MAF) (e.g., rs12573221 ~13% in Bangladesh ~0.2% among American Indians), and/or heterogeneity in effect sizes across populations likely contributed to the apparent population specificity of the additional identified signals. One of our potential causal variants influences AS3MT expression and nearby DNA methylation in numerous GTEx tissue types (with rs4919690 as a likely causal variant). Several SNPs in our confidence sets overlap transcription factor binding sites and cis-regulatory elements (from ENCODE). Taken together, our analyses reveal multiple potential causal variants in the 10q24.32 region influencing AME, including a variant shared across populations, and elucidate potential biological mechanisms underlying the impact of genetic variation on AME.

Systematic evaluation of transcriptomics-based deconvolution methods and references using thousands of clinical samples.

Estimating cell type composition of blood and tissue samples is a biological challenge relevant in both laboratory studies and clinical care. In recent years, a number of computational tools have been developed to estimate cell type abundance using gene expression data. Although these tools use a variety of approaches, they all leverage expression profiles from purified cell types to evaluate the cell type composition within samples. In this study, we compare 12 cell type quantification tools and evaluate their performance while using each of 10 separate reference profiles. Specifically, we have run each tool on over 4000 samples with known cell type proportions, spanning both immune and stromal cell types. A total of 12 of these represent in vitro synthetic mixtures and 300 represent in silico synthetic mixtures prepared using single-cell data. A final 3728 clinical samples have been collected from the Framingham cohort, for which cell populations have been quantified using electrical impedance cell counting. When tools are applied to the Framingham dataset, the tool Estimating the Proportions of Immune and Cancer cells (EPIC) produces the highest correlation, whereas Gene Expression Deconvolution Interactive Tool (GEDIT) produces the lowest error. The best tool for other datasets is varied, but CIBERSORT and GEDIT most consistently produce accurate results. We find that optimal reference depends on the tool used, and report suggested references to be used with each tool. Most tools return results within minutes, but on large datasets runtimes for CIBERSORT can exceed hours or even days. We conclude that deconvolution methods are capable of returning high-quality results, but that proper reference selection is critical.

Recurrent somatic mutations reveal new insights into consequences of mutagenic processes in cancer.

The sheer size of the human genome makes it improbable that identical somatic mutations at the exact same position are observed in multiple tumours solely by chance. The scarcity of cancer driver mutations also precludes positive selection as the sole explanation. Therefore, recurrent mutations may be highly informative of characteristics of mutational processes. To explore the potential, we use recurrence as a starting point to cluster >2,500 whole genomes of a pan-cancer cohort. We describe each genome with 13 recurrence-based and 29 general mutational features. Using principal component analysis we reduce the dimensionality and create independent features. We apply hierarchical clustering to the first 18 principal components followed by k-means clustering. We show that the resulting 16 clusters capture clinically relevant cancer phenotypes. High levels of recurrent substitutions separate the clusters that we link to UV-light exposure and deregulated activity of POLE from the one representing defective mismatch repair, which shows high levels of recurrent insertions/deletions. Recurrence of both mutation types characterizes cancer genomes with somatic hypermutation of immunoglobulin genes and the cluster of genomes exposed to gastric acid. Low levels of recurrence are observed for the cluster where tobacco-smoke exposure induces mutagenesis and the one linked to increased activity of cytidine deaminases. Notably, the majority of substitutions are recurrent in a single tumour type, while recurrent insertions/deletions point to shared processes between tumour types. Recurrence also reveals susceptible sequence motifs, including TT[C>A]TTT and AAC[T>G]T for the POLE and 'gastric-acid exposure' clusters, respectively. Moreover, we refine knowledge of mutagenesis, including increased C/G deletion levels in general for lung tumours and specifically in midsize homopolymer sequence contexts for microsatellite instable tumours. Our findings are an important step towards the development of a generic cancer diagnostic test for clinical practice based on whole-genome sequencing that could replace multiple diagnostics currently in use.

Validation and genotyping of multiple human polymorphic inversions mediated by inverted repeats reveals a high degree of recurrence.

In recent years different types of structural variants (SVs) have been discovered in the human genome and their functional impact has become increasingly clear. Inversions, however, are poorly characterized and more difficult to study, especially those mediated by inverted repeats or segmental duplications. Here, we describe the results of a simple and fast inverse PCR (iPCR) protocol for high-throughput genotyping of a wide variety of inversions using a small amount of DNA. In particular, we analyzed 22 inversions predicted in humans ranging from 5.1 kb to 226 kb and mediated by inverted repeat sequences of 1.6-24 kb. First, we validated 17 of the 22 inversions in a panel of nine HapMap individuals from different populations, and we genotyped them in 68 additional individuals of European origin, with correct genetic transmission in ∼ 12 mother-father-child trios. Global inversion minor allele frequency varied between 1% and 49% and inversion genotypes were consistent with Hardy-Weinberg equilibrium. By analyzing the nucleotide variation and the haplotypes in these regions, we found that only four inversions have linked tag-SNPs and that in many cases there are multiple shared SNPs between standard and inverted chromosomes, suggesting an unexpected high degree of inversion recurrence during human evolution. iPCR was also used to check 16 of these inversions in four chimpanzees and two gorillas, and 10 showed both orientations either within or between species, providing additional support for their multiple origin. Finally, we have identified several inversions that include genes in the inverted or breakpoint regions, and at least one disrupts a potential coding gene. Thus, these results represent a significant advance in our understanding of inversion polymorphism in human populations and challenge the common view of a single origin of inversions, with important implications for inversion analysis in SNP-based studies.

Integrative cross-omics and cross-context analysis elucidates molecular links underlying genetic effects on complex traits.

Genetic effects on functionally related 'omic' traits often co-occur in relevant cellular contexts, such as tissues. Motivated by the multi-tissue methylation quantitative trait loci (mQTLs) and expression QTLs (eQTLs) analysis, we propose X-ING (Cross-INtegrative Genomics) for cross-omics and cross-context integrative analysis. X-ING takes as input multiple matrices of association statistics, each obtained from different omics data types across multiple cellular contexts. It models the latent binary association status of each statistic, captures the major association patterns among omics data types and contexts, and outputs the posterior mean and probability for each input statistic. X-ING enables the integration of effects from different omics data with varying effect distributions. In the multi-tissue cis-association analysis, X-ING shows improved detection and replication of mQTLs by integrating eQTL maps. In the trans-association analysis, X-ING reveals an enrichment of trans-associations in many disease/trait-relevant tissues.

DNA methylation correlates of chronological age in diverse human tissue types.

BACKGROUND: While the association of chronological age with DNA methylation (DNAm) in whole blood has been extensively studied, the tissue-specificity of age-related DNAm changes remains an active area of research. Studies investigating the association of age with DNAm in tissues such as brain, skin, immune cells, fat, and liver have identified tissue-specific and non-specific effects, thus, motivating additional studies of diverse human tissue and cell types. RESULTS: Here, we performed an epigenome-wide association study, leveraging DNAm data (Illumina EPIC array) from 961 tissue samples representing 9 tissue types (breast, lung, colon, ovary, prostate, skeletal muscle, testis, whole blood, and kidney) from the Genotype-Tissue Expression (GTEx) project. We identified age-associated CpG sites (false discovery rate < 0.05) in 8 tissues (all except skeletal muscle, n = 47). This included 162,002 unique hypermethylated and 90,626 hypomethylated CpG sites across all tissue types, with 130,137 (80%) hypermethylated CpGs and 74,703 (82%) hypomethylated CpG sites observed in a single tissue type. While the majority of age-associated CpG sites appeared tissue-specific, the patterns of enrichment among genomic features, such as chromatin states and CpG islands, were similar across most tissues, suggesting common mechanisms underlying cellular aging. Consistent with previous findings, we observed that hypermethylated CpG sites are enriched in regions with repressed polycomb signatures and CpG islands, while hypomethylated CpG sites preferentially occurred in non-CpG islands and enhancers. To gain insights into the functional effects of age-related DNAm changes, we assessed the correlation between DNAm and local gene expression changes to identify age-related expression quantitative trait methylation (age-eQTMs). We identified several age-eQTMs present in multiple tissue-types, including in the CDKN2A, HENMT1, and VCWE regions. CONCLUSION: Overall, our findings will aid future efforts to develop biomarkers of aging and understand mechanisms of aging in diverse human tissue types.

The Impact of Inherited Genetic Variation on DNA Methylation in Prostate Cancer and Benign Tissues of African American and European American Men.

BACKGROUND: American men of African ancestry (AA) have higher prostate cancer incidence and mortality rates compared with American men of European ancestry (EA). Differences in genetic susceptibility mechanisms may contribute to this disparity. METHODS: To gain insights into the regulatory mechanisms of prostate cancer susceptibility variants, we tested the association between SNPs and DNA methylation (DNAm) at nearby CpG sites across the genome in benign and cancer prostate tissue from 74 AA and 74 EA men. Genome-wide SNP data (from benign tissue) and DNAm were generated using Illumina arrays. RESULTS: Among AA men, we identified 6,298 and 2,641 cis-methylation QTLs (meQTL; FDR of 0.05) in benign and tumor tissue, respectively, with 6,960 and 1,700 detected in EA men. We leveraged genome-wide association study (GWAS) summary statistics to identify previously reported prostate cancer GWAS signals likely to share a common causal variant with a detected meQTL. We identified nine GWAS-meQTL pairs with strong evidence of colocalization (four in EA benign, three in EA tumor, two in AA benign, and three in AA tumor). Among these colocalized GWAS-meQTL pairs, we identified colocalizing expression quantitative trait loci (eQTL) impacting four eGenes with known roles in tumorigenesis. CONCLUSIONS: These findings highlight epigenetic regulatory mechanisms by which prostate cancer-risk SNPs can modify local DNAm and/or gene expression in prostate tissue. IMPACT: Overall, our findings showed general consistency in the meQTL landscape of AA and EA men, but meQTLs often differ by tissue type (normal vs. cancer). Ancestry-based linkage disequilibrium differences and lack of AA representation in GWAS decrease statistical power to detect colocalization for some regions.

Multi-ancestry Genome-Wide Association Meta-Analysis Identifies Novel Loci in Atopic Dermatitis.

Atopic dermatitis (AD) is a highly heritable and common inflammatory skin condition affecting children and adults worldwide. Multi-ancestry approaches to AD genetic association studies are poised to boost power to detect genetic signal and identify ancestry-specific loci contributing to AD risk. Here, we present a multi-ancestry GWAS meta-analysis of twelve AD cohorts from five ancestral populations totaling 56,146 cases and 602,280 controls. We report 101 genomic loci associated with AD, including 15 loci that have not been previously associated with AD or eczema. Fine-mapping, QTL colocalization, and cell-type enrichment analyses identified genes and cell types implicated in AD pathophysiology. Functional analyses in keratinocytes provide evidence for genes that could play a role in AD through epidermal barrier function. Our study provides new insights into the etiology of AD by harnessing multiple genetic and functional approaches to unveil the mechanisms by which AD-associated variants impact genes and cell types.

DNA methylation QTL mapping across diverse human tissues provides molecular links between genetic variation and complex traits.

Studies of DNA methylation (DNAm) in solid human tissues are relatively scarce; tissue-specific characterization of DNAm is needed to understand its role in gene regulation and its relevance to complex traits. We generated array-based DNAm profiles for 987 human samples from the Genotype-Tissue Expression (GTEx) project, representing 9 tissue types and 424 subjects. We characterized methylome and transcriptome correlations (eQTMs), genetic regulation in cis (mQTLs and eQTLs) across tissues and e/mQTLs links to complex traits. We identified mQTLs for 286,152 CpG sites, many of which (>5%) show tissue specificity, and mQTL colocalizations with 2,254 distinct GWAS hits across 83 traits. For 91% of these loci, a candidate gene link was identified by integration of functional maps, including eQTMs, and/or eQTL colocalization, but only 33% of loci involved an eQTL and mQTL present in the same tissue type. With this DNAm-focused integrative analysis, we contribute to the understanding of molecular regulatory mechanisms in human tissues and their impact on complex traits.

5-Hydroxymethylcytosine Profiles in Circulating Cell-Free DNA Associate with Disease Burden in Children with Neuroblastoma.

PURPOSE: 5-Hydroxymethylcytosine (5-hmC) is an epigenetic marker of open chromatin and active gene expression. We profiled 5-hmC with Nano-hmC-Seal technology using 10 ng of plasma-derived cell-free DNA (cfDNA) in blood samples from patients with neuroblastoma to determine its utility as a biomarker. EXPERIMENTAL DESIGN: For the Discovery cohort, 100 5-hmC profiles were generated from 34 well children and 32 patients (27 high-risk, 2 intermediate-risk, and 3 low-risk) at various time points during the course of their disease. An independent Validation cohort encompassed 5-hmC cfDNA profiles (n = 29) generated from 21 patients (20 high-risk and 1 intermediate-risk). Metastatic burden was classified as high, moderate, low, or none per Curie metaiodobenzylguanidine scores and percentage of tumor cells in bone marrow. Genes with differential 5-hmC levels between samples according to metastatic burden were identified using DESeq2. RESULTS: Hierarchical clustering using 5-hmC levels of 347 genes identified from the Discovery cohort defined four clusters of samples that were confirmed in the Validation cohort and corresponded to high, high-moderate, moderate, and low/no metastatic burden. Samples from patients with increased metastatic burden had increased 5-hmC deposition on genes in neuronal stem cell maintenance and epigenetic regulatory pathways. Further, 5-hmC cfDNA profiles generated with 1,242 neuronal pathway genes were associated with subsequent relapse in the cluster of patients with predominantly low or no metastatic burden (sensitivity 65%, specificity 75.6%). CONCLUSIONS: cfDNA 5-hmC profiles in children with neuroblastoma correlate with metastatic burden and warrants development as a biomarker of treatment response and outcome.

Determinants of telomere length across human tissues.

Telomere shortening is a hallmark of aging. Telomere length (TL) in blood cells has been studied extensively as a biomarker of human aging and disease; however, little is known regarding variability in TL in nonblood, disease-relevant tissue types. Here, we characterize variability in TLs from 6391 tissue samples, representing >20 tissue types and 952 individuals from the Genotype-Tissue Expression (GTEx) project. We describe differences across tissue types, positive correlation among tissue types, and associations with age and ancestry. We show that genetic variation affects TL in multiple tissue types and that TL may mediate the effect of age on gene expression. Our results provide the foundational knowledge regarding TL in healthy tissues that is needed to interpret epidemiological studies of TL and human health.

Cell type-specific genetic regulation of gene expression across human tissues.

The Genotype-Tissue Expression (GTEx) project has identified expression and splicing quantitative trait loci in cis (QTLs) for the majority of genes across a wide range of human tissues. However, the functional characterization of these QTLs has been limited by the heterogeneous cellular composition of GTEx tissue samples. We mapped interactions between computational estimates of cell type abundance and genotype to identify cell type-interaction QTLs for seven cell types and show that cell type-interaction expression QTLs (eQTLs) provide finer resolution to tissue specificity than bulk tissue cis-eQTLs. Analyses of genetic associations with 87 complex traits show a contribution from cell type-interaction QTLs and enables the discovery of hundreds of previously unidentified colocalized loci that are masked in bulk tissue.

The impact of sex on gene expression across human tissues.

Many complex human phenotypes exhibit sex-differentiated characteristics. However, the molecular mechanisms underlying these differences remain largely unknown. We generated a catalog of sex differences in gene expression and in the genetic regulation of gene expression across 44 human tissue sources surveyed by the Genotype-Tissue Expression project (GTEx, v8 release). We demonstrate that sex influences gene expression levels and cellular composition of tissue samples across the human body. A total of 37% of all genes exhibit sex-biased expression in at least one tissue. We identify cis expression quantitative trait loci (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 of gene expression in a single sex. These findings provide an extensive characterization of sex differences in the human transcriptome and its genetic regulation.

Evolutionary and functional impact of common polymorphic inversions in the human genome.

Inversions are one type of structural variants linked to phenotypic differences and adaptation in multiple organisms. However, there is still very little information about polymorphic inversions in the human genome due to the difficulty of their detection. Here, we develop a new high-throughput genotyping method based on probe hybridization and amplification, and we perform a complete study of 45 common human inversions of 0.1-415 kb. Most inversions promoted by homologous recombination occur recurrently in humans and great apes and they are not tagged by SNPs. Furthermore, there is an enrichment of inversions showing signatures of positive or balancing selection, diverse functional effects, such as gene disruption and gene-expression changes, or association with phenotypic traits. Therefore, our results indicate that the genome is more dynamic than previously thought and that human inversions have important functional and evolutionary consequences, making possible to determine for the first time their contribution to complex traits.

Evaluación del traspaso de información (Hand Off) en equipos de enfermería de urgencias / Evaluation of the transfer of information (Hand Off) in emergency nursing teams

RESUMEN Introducción: los traspasos (Hand Off) de enfermería en urgencias aseguran la transmisión de información crítica, la continuidad de los cuidados y el tratamiento del paciente. Las incidencias o deficiencias en los Traspasos pueden ocasionar problemas en la seguridad del paciente, por lo que es necesario investigar y desarrollar estrategias para reducirlos. Objetivo: identificar el contenido y las deficiencias en los traspasos de enfermería en urgencias. Métodos: estudio descriptivo observacional prospectivo en urgencias del Hospital Clínic de Barcelona durante 2º trimestre 2016 sobre los 95265 traspasos de enfermería anuales realizados. Muestra de conveniencia calculada n=384. Se revisaron los traspasos de equipos de enfermería según modelo de 24 ítems creado por investigadoras basado en evidencia. Resultados: sobre 24 ítems del modelo de referencia de traspaso, se transmitieron media 7,61 (IC 95 por ciento 7,41 7,81), el 31,7 por ciento del total. Los ítems no comunicados fueron de media 16,38 (IC 95 por ciento 16,18 16,58), el 68,3 por ciento. Se obtuvieron diferencias significativas (p<0,001) según el Área de Urgencias evaluada (perfil de pacientes y cuidados diferenciados). Los ítems no comunicados en valores medios 16,24 (IC95 por ciento 15,97 16,52) en Observación Medicina; 16,74 (IC95 por ciento 16,16 17,32) en Observación Cirugía/Trauma; 15,72 (IC95 por ciento 15,20 16,25) en Emergencias Nivel 2; 16,85 (IC95 por ciento 16,37 17,34) en Urgencias Nivel 3; 17,23 (IC95 por ciento 16,77 17,69) en Urgencias Traumatología. Conclusiones: los traspasos de enfermería presentaron deficiencias en su contenido, con diferencias por área asistencial. La estandarización del traspaso en equipos de enfermería de urgencias puede mejorar su calidad, su eficiencia y repercutir en una mayor seguridad para el paciente(AU)

ABSTRACT Introduction: Emergency Nursing Hand off ensures the transmission of critical information and continuity of care and treatment of patients. Incidents or deficiencies in Hand off can cause problems in patient safety so it is necessary to investigate and develop strategies to reduce them. Objective: To identify the content and deficiencies in the Emergency Nursing Hand off. Methods: Prospective descriptive study in Emergency Service, Hospital Clínic of Barcelona during 2nd quarter 2016 over 95265 annual transfers of nursing made. Calculated convenience sample n = 384. Emergency Nursing Teams Hand off were reviewed by 24-item model created by researchers evidence based. Results: About 24 items reference Hand off model, were transmitted average 7.61 (95 percent CI 7.41 7.81), 31.7 percent of the total. The items were not reported average 16.38 (95 percent 16.18 16.58), 68.3 percent. Significant differences (p <0.001) were obtained according to the Emergency Department assessed (differentiated profile of patients and care). Items not reported mean = 16.24 (95 percent CI 15.97 16.52) in Observation Medicine; mean = 16.85 (95 percent CI 16.37 17.34) in Level 3 Emergency Room; mean = 16.74 (95 percent CI 16.16 17.32) in Observation Surgery / Trauma; mean = 15.72 (95 percent CI 15,20 16,25) in Level 2 Emergency Room; mean = 17.23 (95 percent CI 17.69 to 16.77) in Traumatology Emergency Room. Conclusions: Emergency Nursing Hand Off showed deficiencies in its content, with differences by healthcare area. Hand off standardization in Emergency Nursing Teams could improve their quality, efficiency and get greater impact on patient safety(AU)