METRO: Multi-ancestry transcriptome-wide association studies for powerful gene-trait association detection.

Integrative analysis of genome-wide association studies (GWASs) and gene expression studies in the form of a transcriptome-wide association study (TWAS) has the potential to better elucidate the molecular mechanisms underlying disease etiology. Here we present a method, METRO, that can leverage gene expression data collected from multiple genetic ancestries to enhance TWASs. METRO incorporates expression prediction models constructed in different genetic ancestries through a likelihood-based inference framework, producing calibrated p values with substantially improved TWAS power. We illustrate the benefits of METRO in both simulations and applications to seven complex traits and diseases obtained from four GWASs. These GWASs include two of primarily European ancestry (n = 188,577 and 339,226) and two of primarily African ancestry (n = 42,752 and 23,827). In the real data applications, we leverage gene expression data measured on 1,032 African Americans and 801 European Americans from the Genetic Epidemiology Network of Arteriopathy (GENOA) study to identify a substantially larger number of gene-trait associations as compared to existing TWAS approaches. The benefits of METRO are most prominent in applications to GWASs of African ancestry where the sample size is much smaller than GWASs of European ancestry and where a more powerful TWAS method is crucial. Among the identified associations are high-density lipoprotein-associated genes including PLTP and PPARG that are critical for maintaining lipid homeostasis and the type II diabetes-associated gene MAPT that supports microtubule-associated protein tau as a key component underlying impaired insulin secretion.

Genetic Architecture of Gene Expression in European and African Americans: An eQTL Mapping Study in GENOA.

Most existing expression quantitative trait locus (eQTL) mapping studies have been focused on individuals of European ancestry and are underrepresented in other populations including populations with African ancestry. Lack of large-scale well-powered eQTL mapping studies in populations with African ancestry can both impede the dissemination of eQTL mapping results that would otherwise benefit individuals with African ancestry and hinder the comparable analysis for understanding how gene regulation is shaped through evolution. We fill this critical knowledge gap by performing a large-scale in-depth eQTL mapping study on 1,032 African Americans (AA) and 801 European Americans (EA) in the GENOA cohort. We identified a total of 354,931 eSNPs in AA and 371,309 eSNPs in EA, with 112,316 eSNPs overlapped between the two. We found that eQTL harboring genes (eGenes) are enriched in metabolic pathways and tend to have higher SNP heritability compared to non-eGenes. We found that eGenes that are common in the two populations tend to be less conserved than eGenes that are unique to one population, which are less conserved than non-eGenes. Through conditional analysis, we found that eGenes in AA tend to harbor more independent eQTLs than eGenes in EA, suggesting potentially diverse genetic architecture underlying expression variation in the two populations. Finally, the large sample sizes in GENOA allow us to construct accurate expression prediction models in both AA and EA, facilitating powerful transcriptome-wide association studies. Overall, our results represent an important step toward revealing the genetic architecture underlying expression variation in African Americans.

Leveraging gene co-expression patterns to infer trait-relevant tissues in genome-wide association studies.

Genome-wide association studies (GWASs) have identified many SNPs associated with various common diseases. Understanding the biological functions of these identified SNP associations requires identifying disease/trait relevant tissues or cell types. Here, we develop a network method, CoCoNet, to facilitate the identification of trait-relevant tissues or cell types. Different from existing approaches, CoCoNet incorporates tissue-specific gene co-expression networks constructed from either bulk or single cell RNA sequencing (RNAseq) studies with GWAS data for trait-tissue inference. In particular, CoCoNet relies on a covariance regression network model to express gene-level effect measurements for the given GWAS trait as a function of the tissue-specific co-expression adjacency matrix. With a composite likelihood-based inference algorithm, CoCoNet is scalable to tens of thousands of genes. We validate the performance of CoCoNet through extensive simulations. We apply CoCoNet for an in-depth analysis of four neurological disorders and four autoimmune diseases, where we integrate the corresponding GWASs with bulk RNAseq data from 38 tissues and single cell RNAseq data from 10 cell types. In the real data applications, we show how CoCoNet can help identify specific glial cell types relevant for neurological disorders and identify disease-targeted colon tissues as relevant for autoimmune diseases.

Identification of functionally important miRNA targeted genes associated with child obesity trait in genome-wide association studies.

BACKGROUND: Genome-wide association studies (GWAS) have uncovered thousands of genetic variants that are associated with complex human traits and diseases. miRNAs are single-stranded non-coding RNAs. In particular, genetic variants located in the 3'UTR region of mRNAs may play an important role in gene regulation through their interaction with miRNAs. Existing studies have not been thoroughly conducted to elucidate 3'UTR variants discovered through GWAS. The goal of this study is to analyze patterns of GWAS functional variants located in 3'UTRs about their relevance in the network between hosting genes and targeting miRNAs, and elucidate the association between the genes harboring these variants and genetic traits. METHODS: We employed MIGWAS, ANNOVAR, MEME, and DAVID software packages to annotate the variants obtained from GWAS for 31 traits and elucidate the association between their harboring genes and their related traits. We identified variants that occurred in the motif regions that may be functionally important in affecting miRNA binding. We also conducted pathway analysis and functional annotation on miRNA targeted genes harboring 3'UTR variants for a trait with the highest percentage of 3'UTR variants occurring. RESULTS: The Child Obesity trait has the highest percentage of 3'UTR variants (75%). Of the 16 genes related to the Child Obesity trait, 5 genes (ETV7, GMEB1, NFIX, ZNF566, ZBTB40) had a significant association with the term DNA-Binding (p < 0.05). EQTL analysis revealed 2 relevant tissues and 10 targeted genes associated with the Child Obesity trait. In addition, Red Blood Cells (RBC), Hemoglobin (HB), and Package Cell Volume (PCV) have overlapping variants. In particular, the PIM1 variant occurred inside the HB Motif region 37,174,641-37,174,660, and LUC7L3 variant occurred inside RBC Motif region 50,753,918-50,753,937. CONCLUSION: Variants located in 3'UTR can alter the binding affinity of miRNA and impact gene regulation, thus warranting further annotation and analysis. We have developed a bioinformatics bash pipeline to automatically annotate variants, determine the number of variants in different categories for each given trait, and check common variants across different traits. This is a valuable tool to annotate a large number of GWAS result files.

Genome-Wide Association Meta-Analysis of Individuals of European Ancestry Identifies Suggestive Loci for Sodium Intake, Potassium Intake, and Their Ratio Measured from 24-Hour or Half-Day Urine Samples.

BACKGROUND: Excess sodium intake and insufficient potassium intake are risk factors for hypertension, but there is limited knowledge regarding genetic factors that influence intake. Twenty-hour or half-day urine samples provide robust estimates of sodium and potassium intake, outperforming other measures such as spot urine samples and dietary self-reporting. OBJECTIVE: The aim of this study was to investigate genomic regions associated with sodium intake, potassium intake, and sodium-to-potassium ratio measured from 24-h or half-day urine samples. METHODS: Using samples of European ancestry (mean age: 54.2 y; 52.3% women), we conducted a meta-analysis of genome-wide association studies in 4 cohorts with 24-h or half-day urine samples (n = 6,519), followed by gene-based analysis. Suggestive loci (P < 10-6) were examined in additional European (n = 844), African (n = 1,246), and Asian (n = 2,475) ancestry samples. RESULTS: We found suggestive loci (P < 10-6) for all 3 traits, including 7 for 24-h sodium excretion, 4 for 24-h potassium excretion, and 4 for sodium-to-potassium ratio. The most significant locus was rs77958157 near cocaine- and amphetamine-regulated transcript prepropeptide (CARTPT) , a gene involved in eating behavior and appetite regulation (P = 2.3 × 10-8 with sodium-to-potassium ratio). Two suggestive loci were replicated in additional samples: for sodium excretion, rs12094702 near zinc finger SWIM-type containing 5 (ZSWIM5) was replicated in the Asian ancestry sample reaching Bonferroni-corrected significance (P = 0.007), and for potassium excretion rs34473523 near sodium leak channel (NALCN) was associated at a nominal P value with potassium excretion both in European (P = 0.043) and African (P = 0.043) ancestry cohorts. Gene-based tests identified 1 significant gene for sodium excretion, CDC42 small effector 1 (CDC42SE1), which is associated with blood pressure regulation. CONCLUSIONS: We identified multiple suggestive loci for sodium and potassium intake near genes associated with eating behavior, nervous system development and function, and blood pressure regulation in individuals of European ancestry. Further research is needed to replicate these findings and to provide insight into the underlying genetic mechanisms by which these genomic regions influence sodium and potassium intake.

SRTsim: spatial pattern preserving simulations for spatially resolved transcriptomics.

Spatially resolved transcriptomics (SRT)-specific computational methods are often developed, tested, validated, and evaluated in silico using simulated data. Unfortunately, existing simulated SRT data are often poorly documented, hard to reproduce, or unrealistic. Single-cell simulators are not directly applicable for SRT simulation as they cannot incorporate spatial information. We present SRTsim, an SRT-specific simulator for scalable, reproducible, and realistic SRT simulations. SRTsim not only maintains various expression characteristics of SRT data but also preserves spatial patterns. We illustrate the benefits of SRTsim in benchmarking methods for spatial clustering, spatial expression pattern detection, and cell-cell communication identification.

meQTL mapping in the GENOA study reveals genetic determinants of DNA methylation in African Americans.

Identifying genetic variants that are associated with variation in DNA methylation, an analysis commonly referred to as methylation quantitative trait locus (meQTL) mapping, is an important first step towards understanding the genetic architecture underlying epigenetic variation. Most existing meQTL mapping studies have focused on individuals of European ancestry and are underrepresented in other populations, with a particular absence of large studies in populations with African ancestry. We fill this critical knowledge gap by performing a large-scale cis-meQTL mapping study in 961 African Americans from the Genetic Epidemiology Network of Arteriopathy (GENOA) study. We identify a total of 4,565,687 cis-acting meQTLs in 320,965 meCpGs. We find that 45% of meCpGs harbor multiple independent meQTLs, suggesting potential polygenic genetic architecture underlying methylation variation. A large percentage of the cis-meQTLs also colocalize with cis-expression QTLs (eQTLs) in the same population. Importantly, the identified cis-meQTLs explain a substantial proportion (median = 24.6%) of methylation variation. In addition, the cis-meQTL associated CpG sites mediate a substantial proportion (median = 24.9%) of SNP effects underlying gene expression. Overall, our results represent an important step toward revealing the co-regulation of methylation and gene expression, facilitating the functional interpretation of epigenetic and gene regulation underlying common diseases in African Americans.

Spatially aware dimension reduction for spatial transcriptomics.

Spatial transcriptomics are a collection of genomic technologies that have enabled transcriptomic profiling on tissues with spatial localization information. Analyzing spatial transcriptomic data is computationally challenging, as the data collected from various spatial transcriptomic technologies are often noisy and display substantial spatial correlation across tissue locations. Here, we develop a spatially-aware dimension reduction method, SpatialPCA, that can extract a low dimensional representation of the spatial transcriptomics data with biological signal and preserved spatial correlation structure, thus unlocking many existing computational tools previously developed in single-cell RNAseq studies for tailored analysis of spatial transcriptomics. We illustrate the benefits of SpatialPCA for spatial domain detection and explores its utility for trajectory inference on the tissue and for high-resolution spatial map construction. In the real data applications, SpatialPCA identifies key molecular and immunological signatures in a detected tumor surrounding microenvironment, including a tertiary lymphoid structure that shapes the gradual transcriptomic transition during tumorigenesis and metastasis. In addition, SpatialPCA detects the past neuronal developmental history that underlies the current transcriptomic landscape across tissue locations in the cortex.

Spatial transcriptomics.

Spatial transcriptomics, with other spatial technologies, has enabled scientists to dissect the organization and interaction of different cell types within the tumor microenvironment. We asked experts to discuss some aspects of this technology from revealing the tumor microenvironment and heterogeneity, to tracking tumor evolution, to guiding tumor therapy, to current technical challenges.

DNA Methylation Mediates the Association Between Individual and Neighborhood Social Disadvantage and Cardiovascular Risk Factors.

Low socioeconomic status (SES) and living in a disadvantaged neighborhood are associated with poor cardiovascular health. Multiple lines of evidence have linked DNA methylation to both cardiovascular risk factors and social disadvantage indicators. However, limited research has investigated the role of DNA methylation in mediating the associations of individual- and neighborhood-level disadvantage with multiple cardiovascular risk factors in large, multi-ethnic, population-based cohorts. We examined whether disadvantage at the individual level (childhood and adult SES) and neighborhood level (summary neighborhood SES as assessed by Census data and social environment as assessed by perceptions of aesthetic quality, safety, and social cohesion) were associated with 11 cardiovascular risk factors including measures of obesity, diabetes, lipids, and hypertension in 1,154 participants from the Multi-Ethnic Study of Atherosclerosis (MESA). For significant associations, we conducted epigenome-wide mediation analysis to identify methylation sites mediating the relationship between individual/neighborhood disadvantage and cardiovascular risk factors using the JT-Comp method that assesses sparse mediation effects under a composite null hypothesis. In models adjusting for age, sex, race/ethnicity, smoking, medication use, and genetic principal components of ancestry, epigenetic mediation was detected for the associations of adult SES with body mass index (BMI), insulin, and high-density lipoprotein cholesterol (HDL-C), as well as for the association between neighborhood socioeconomic disadvantage and HDL-C at FDR q < 0.05. The 410 CpG mediators identified for the SES-BMI association were enriched for CpGs associated with gene expression (expression quantitative trait methylation loci, or eQTMs), and corresponding genes were enriched in antigen processing and presentation pathways. For cardiovascular risk factors other than BMI, most of the epigenetic mediators lost significance after controlling for BMI. However, 43 methylation sites showed evidence of mediating the neighborhood socioeconomic disadvantage and HDL-C association after BMI adjustment. The identified mediators were enriched for eQTMs, and corresponding genes were enriched in inflammatory and apoptotic pathways. Our findings support the hypothesis that DNA methylation acts as a mediator between individual- and neighborhood-level disadvantage and cardiovascular risk factors, and shed light on the potential underlying epigenetic pathways. Future studies are needed to fully elucidate the biological mechanisms that link social disadvantage to poor cardiovascular health.

Epigenome-wide association study identifies DNA methylation sites associated with target organ damage in older African Americans.

Target organ damage (TOD) manifests as vascular injuries in the body organ systems associated with long-standing hypertension. DNA methylation in peripheral blood leukocytes can capture inflammatory processes and gene expression changes underlying TOD. We investigated the association between epigenome-wide DNA methylation and five measures of TOD (estimated glomerular filtration rate (eGFR), urinary albumin-creatinine ratio (UACR), left ventricular mass index (LVMI), relative wall thickness (RWT), and white matter hyperintensity (WMH)) in 961 African Americans from hypertensive sibships. A multivariate (multi-trait) model of eGFR, UACR, LVMI, and RWT identified seven CpGs associated with at least one of the traits (cg21134922, cg04816311 near C7orf50, cg09155024, cg10254690 near OAT, cg07660512, cg12661888 near IFT43, and cg02264946 near CATSPERD) at FDR q < 0.1. Adjusting for blood pressure, body mass index, and type 2 diabetes attenuated the association for four CpGs. DNA methylation was associated with cis-gene expression for some CpGs, but no significant mediation by gene expression was detected. Mendelian randomization analyses suggested causality between three CpGs and eGFR (cg04816311, cg10254690, and cg07660512). We also assessed whether the identified CpGs were associated with TOD in 614 African Americans in the Hypertension Genetic Epidemiology Network (HyperGEN) study. Out of three CpGs available for replication, cg04816311 was significantly associated with eGFR (p = 0.0003), LVMI (p = 0.0003), and RWT (p = 0.002). This study found evidence of an association between DNA methylation and TOD in African Americans and highlights the utility of using a multivariate-based model that leverages information across related traits in epigenome-wide association studies.

Transcriptome-Wide Association Study of Blood Cell Traits in African Ancestry and Hispanic/Latino Populations.

BACKGROUND: Thousands of genetic variants have been associated with hematological traits, though target genes remain unknown at most loci. Moreover, limited analyses have been conducted in African ancestry and Hispanic/Latino populations; hematological trait associated variants more common in these populations have likely been missed. METHODS: To derive gene expression prediction models, we used ancestry-stratified datasets from the Multi-Ethnic Study of Atherosclerosis (MESA, including n = 229 African American and n = 381 Hispanic/Latino participants, monocytes) and the Depression Genes and Networks study (DGN, n = 922 European ancestry participants, whole blood). We then performed a transcriptome-wide association study (TWAS) for platelet count, hemoglobin, hematocrit, and white blood cell count in African (n = 27,955) and Hispanic/Latino (n = 28,324) ancestry participants. RESULTS: Our results revealed 24 suggestive signals (p < 1 × 10-4) that were conditionally distinct from known GWAS identified variants and successfully replicated these signals in European ancestry subjects from UK Biobank. We found modestly improved correlation of predicted and measured gene expression in an independent African American cohort (the Genetic Epidemiology Network of Arteriopathy (GENOA) study (n = 802), lymphoblastoid cell lines) using the larger DGN reference panel; however, some genes were well predicted using MESA but not DGN. CONCLUSIONS: These analyses demonstrate the importance of performing TWAS and other genetic analyses across diverse populations and of balancing sample size and ancestry background matching when selecting a TWAS reference panel.

A Review of Statistical Methods for Identifying Trait-Relevant Tissues and Cell Types.

Genome-wide association studies (GWASs) have identified and replicated many genetic variants that are associated with diseases and disease-related complex traits. However, the biological mechanisms underlying these identified associations remain largely elusive. Exploring the biological mechanisms underlying these associations requires identifying trait-relevant tissues and cell types, as genetic variants likely influence complex traits in a tissue- and cell type-specific manner. Recently, several statistical methods have been developed to integrate genomic data with GWASs for identifying trait-relevant tissues and cell types. These methods often rely on different genomic information and use different statistical models for trait-tissue relevance inference. Here, we present a comprehensive technical review to summarize ten existing methods for trait-tissue relevance inference. These methods make use of different genomic information that include functional annotation information, expression quantitative trait loci information, genetically regulated gene expression information, as well as gene co-expression network information. These methods also use different statistical models that range from linear mixed models to covariance network models. We hope that this review can serve as a useful reference both for methodologists who develop methods and for applied analysts who apply these methods for identifying trait relevant tissues and cell types.

Pinpointing miRNA and genes enrichment over trait-relevant tissue network in Genome-Wide Association Studies.

BACKGROUND: Understanding gene regulation is important but difficult. Elucidating tissue-specific gene regulation mechanism is even more challenging and requires gene co-expression network assembled from protein-protein interaction, transcription factor and gene binding, and post-transcriptional regulation (e.g., miRNA targeting) information. The miRNA binding affinity could therefore be changed by SNP(s) located at the 3' untranslated regions (3'UTR) of the target messenger RNA (mRNA) which miRNA(s) interacts with. Genome-wide association study (GWAS) has reported significant numbers of loci hosting SNPs associated with many traits. The goal of this study is to pinpoint GWAS functional variants located in 3'UTRs and elucidate if the genes harboring these variants along with their targeting miRNAs are associated with genetic traits relevant to certain tissues. METHODS: By applying MIGWAS, CoCoNet, ANNOVAR, and DAVID bioinformatics software and utilizing the gene expression database (e.g. GTEx data) to study GWAS summary statistics for 43 traits from 28 GWAS studies, we have identified a list of miRNAs and targeted genes harboring 3'UTR variants, which could contribute to trait-relevant tissue over miRNA-target gene network. RESULTS: Our result demonstrated that strong association between traits and tissues exists, and in particular, the Primary Biliary Cirrhosis (PBC) trait has the most significant p-value for all 180 tissues among all 43 traits used for this study. We reported SNPs located in 3'UTR regions of genes (SFMBT2, ZC3HAV1, and UGT3A1) targeted by miRNAs for PBC trait and its tissue association network. After employing Gene Ontology (GO) analysis for PBC trait, we have also identified a very important miRNA targeted gene over miRNA-target gene network, PFKL, which encodes the liver subunit of an enzyme. CONCLUSIONS: The non-coding variants identified from GWAS studies are casually assumed to be not critical to translated protein product. However, 3' untranslated regions (3'UTRs) of genes harbor variants can often change the binding affinity of targeting miRNAs playing important roles in protein translation degree. Our study has shown that GWAS variants could play important roles on miRNA-target gene networks by contributing the association between traits and tissues. Our analysis expands our knowledge on trait-relevant tissue network and paves way for future human disease studies.

Epigenetic loci for blood pressure are associated with hypertensive target organ damage in older African Americans from the genetic epidemiology network of Arteriopathy (GENOA) study.

BACKGROUND: Hypertension is a major modifiable risk factor for arteriosclerosis that can lead to target organ damage (TOD) of heart, kidneys, and peripheral arteries. A recent epigenome-wide association study for blood pressure (BP) identified 13 CpG sites, but it is not known whether DNA methylation at these sites is also associated with TOD. METHODS: In 1218 African Americans from the Genetic Epidemiology Network of Arteriopathy (GENOA) study, a cohort of hypertensive sibships, we evaluated the associations between methylation at these 13 CpG sites measured in peripheral blood leukocytes and five TOD traits assessed approximately 5 years later. RESULTS: Ten significant associations were found after adjustment for age, sex, blood cell counts, time difference between CpG and TOD measurement, and 10 genetic principal components (FDR q < 0.1): two with estimated glomerular filtration rate (eGFR, cg06690548, cg10601624), six with urinary albumin-to-creatinine ratio (UACR, cg16246545, cg14476101, cg19693031, cg06690548, cg00574958, cg22304262), and two with left ventricular mass indexed to height (LVMI, cg19693031, cg00574958). All associations with eGFR and four associations with UACR remained significant after further adjustment for body mass index (BMI), smoking status, and diabetes. We also found significant interactions between cg06690548 and BMI on UACR, and between 3 CpG sites (cg19693031, cg14476101, and cg06690548) and diabetes on UACR (FDR q < 0.1). Mediation analysis showed that 4.7% to 38.1% of the relationship between two CpG sites (cg19693031 and cg00574958) and two TOD measures (UACR and LVMI) was mediated by blood pressure (Bonferroni-corrected P < 0.05). Mendelian randomization analysis suggests that methylation at two sites (cg16246545 and cg14476101) in PHGDH may causally influence UACR. CONCLUSIONS: In conclusion, we found compelling evidence for associations between arteriosclerotic traits of kidney and heart and previously identified blood pressure-associated DNA methylation sites. This study may lend insight into the role of DNA methylation in pathological mechanisms underlying target organ damage from hypertension.

Ranking Cancer Proteins by Integrating PPI Network and Protein Expression Profiles.

Proteomics, the large-scale analysis of proteins, is contributing greatly to understanding gene function in the postgenomic era. However, disease protein ranking using shotgun proteomics data has not been fully evaluated. In this study, we prioritized disease-related proteins by integrating the protein-protein interaction (PPI) network and protein differential expression profiles from colon and rectal cancer (CRC) or breast cancer (BC) proteomics. We applied Local Ranking (LR) and Global Ranking (GR) methods in network with three kinds of protein sets as a priori knowledge, which were known disease proteins (KDPs) that were collected from the Online Mendelian Inheritance in Man (OMIM) database, differentially expressed proteins (DEPs), and the collection of KDPs and their direct neighborhood with differential expression (eKDPs). The cross-validations showed that GR method outperformed LR method while using eKDPs as the initial training showed significantly higher accuracy compared to using the other two a priori sets. And then we validated the top ranked proteins using RNAi-based loss-of-function screens in the DepMap database. The results showed that 75% of top 20 proteins in CRC are necessary for tumor survival. In summary, the network-based Global Ranking with protein differential expression can efficiently prioritize cancer-related proteins and discover new candidate cancer genes or proteins.