Proteomic networks and related genetic variants associated with smoking and chronic obstructive pulmonary disease.

BACKGROUND: Studies have identified individual blood biomarkers associated with chronic obstructive pulmonary disease (COPD) and related phenotypes. However, complex diseases such as COPD typically involve changes in multiple molecules with interconnections that may not be captured when considering single molecular features. METHODS: Leveraging proteomic data from 3,173 COPDGene Non-Hispanic White (NHW) and African American (AA) participants, we applied sparse multiple canonical correlation network analysis (SmCCNet) to 4,776 proteins assayed on the SomaScan v4.0 platform to derive sparse networks of proteins associated with current vs. former smoking status, airflow obstruction, and emphysema quantitated from high-resolution computed tomography scans. We then used NetSHy, a dimension reduction technique leveraging network topology, to produce summary scores of each proteomic network, referred to as NetSHy scores. We next performed a genome-wide association study (GWAS) to identify variants associated with the NetSHy scores, or network quantitative trait loci (nQTLs). Finally, we evaluated the replicability of the networks in an independent cohort, SPIROMICS. RESULTS: We identified networks of 13 to 104 proteins for each phenotype and exposure in NHW and AA, and the derived NetSHy scores significantly associated with the variable of interests. Networks included known (sRAGE, ALPP, MIP1) and novel molecules (CA10, CPB1, HIS3, PXDN) and interactions involved in COPD pathogenesis. We observed 7 nQTL loci associated with NetSHy scores, 4 of which remained after conditional analysis. Networks for smoking status and emphysema, but not airflow obstruction, demonstrated a high degree of replicability across race groups and cohorts. CONCLUSIONS: In this work, we apply state-of-the-art molecular network generation and summarization approaches to proteomic data from COPDGene participants to uncover protein networks associated with COPD phenotypes. We further identify genetic associations with networks. This work discovers protein networks containing known and novel proteins and protein interactions associated with clinically relevant COPD phenotypes across race groups and cohorts.

Type 1 Diabetes Genetic Risk in 109,954 Veterans With Adult-Onset Diabetes: The Million Veteran Program (MVP).

OBJECTIVE: To characterize high type 1 diabetes (T1D) genetic risk in a population where type 2 diabetes (T2D) predominates. RESEARCH DESIGN AND METHODS: Characteristics typically associated with T1D were assessed in 109,594 Million Veteran Program participants with adult-onset diabetes, 2011-2021, who had T1D genetic risk scores (GRS) defined as low (0 to <45%), medium (45 to <90%), high (90 to <95%), or highest (≥95%). RESULTS: T1D characteristics increased progressively with higher genetic risk (P < 0.001 for trend). A GRS ≥90% was more common with diabetes diagnoses before age 40 years, but 95% of those participants were diagnosed at age ≥40 years, and their characteristics resembled those of individuals with T2D in mean age (64.3 years) and BMI (32.3 kg/m2). Compared with the low-risk group, the highest-risk group was more likely to have diabetic ketoacidosis (low GRS 0.9% vs. highest GRS 3.7%), hypoglycemia prompting emergency visits (3.7% vs. 5.8%), outpatient plasma glucose <50 mg/dL (7.5% vs. 13.4%), a shorter median time to start insulin (3.5 vs. 1.4 years), use of a T1D diagnostic code (16.3% vs. 28.1%), low C-peptide levels if tested (1.8% vs. 32.4%), and glutamic acid decarboxylase antibodies (6.9% vs. 45.2%), all P < 0.001. CONCLUSIONS: Characteristics associated with T1D were increased with higher genetic risk, and especially with the top 10% of risk. However, the age and BMI of those participants resemble those of people with T2D, and a substantial proportion did not have diagnostic testing or use of T1D diagnostic codes. T1D genetic screening could be used to aid identification of adult-onset T1D in settings in which T2D predominates.

Proteomic Networks and Related Genetic Variants Associated with Smoking and Chronic Obstructive Pulmonary Disease.

Background: Studies have identified individual blood biomarkers associated with chronic obstructive pulmonary disease (COPD) and related phenotypes. However, complex diseases such as COPD typically involve changes in multiple molecules with interconnections that may not be captured when considering single molecular features. Methods: Leveraging proteomic data from 3,173 COPDGene Non-Hispanic White (NHW) and African American (AA) participants, we applied sparse multiple canonical correlation network analysis (SmCCNet) to 4,776 proteins assayed on the SomaScan v4.0 platform to derive sparse networks of proteins associated with current vs. former smoking status, airflow obstruction, and emphysema quantitated from high-resolution computed tomography scans. We then used NetSHy, a dimension reduction technique leveraging network topology, to produce summary scores of each proteomic network, referred to as NetSHy scores. We next performed genome-wide association study (GWAS) to identify variants associated with the NetSHy scores, or network quantitative trait loci (nQTLs). Finally, we evaluated the replicability of the networks in an independent cohort, SPIROMICS. Results: We identified networks of 13 to 104 proteins for each phenotype and exposure in NHW and AA, and the derived NetSHy scores significantly associated with the variable of interests. Networks included known (sRAGE, ALPP, MIP1) and novel molecules (CA10, CPB1, HIS3, PXDN) and interactions involved in COPD pathogenesis. We observed 7 nQTL loci associated with NetSHy scores, 4 of which remained after conditional analysis. Networks for smoking status and emphysema, but not airflow obstruction, demonstrated a high degree of replicability across race groups and cohorts. Conclusions: In this work, we apply state-of-the-art molecular network generation and summarization approaches to proteomic data from COPDGene participants to uncover protein networks associated with COPD phenotypes. We further identify genetic associations with networks. This work discovers protein networks containing known and novel proteins and protein interactions associated with clinically relevant COPD phenotypes across race groups and cohorts.

Hippocampus Glutathione S Reductase Potentially Confers Genetic Resilience to Cognitive Decline in the AD-BXD Mouse Population.

Alzheimer's disease (AD) is a prevalent and costly age-related dementia. Heritable factors account for 58-79% of variation in late-onset AD, but substantial variation remains in age-of- onset, disease severity, and whether those with high-risk genotypes acquire AD. To emulate the diversity of human populations, we utilized the AD-BXD mouse panel. This genetically diverse resource combines AD genotypes with multiple BXD strains to discover new genetic drivers of AD resilience. Comparing AD-BXD carriers to noncarrier littermates, we computed a novel quantitative metric for resilience to cognitive decline in the AD-BXDs. Our quantitative AD resilience trait was heritable and genetic mapping identified a locus on chr8 associated with resilience to AD mutations that resulted in amyloid brain pathology. Using a hippocampus proteomics dataset, we nominated the mitochondrial glutathione S reductase protein (GR or GSHR) as a resilience factor, finding that the DBA/2J genotype was associated with substantially higher GR abundance. By mapping protein QTLs (pQTLs), we identified synaptic organization and mitochondrial proteins coregulated in trans with a cis-pQTL for GR. We found four coexpression modules correlated with the quantitative resilience score in aged 5XFAD mice using paracliques, which were related to cell structure, protein folding, and postsynaptic densities. Finally, we found significant positive associations between human GSR transcript abundance in the brain and better outcomes on AD-related cognitive and pathology traits in the Religious Orders Study/Memory and Aging project (ROSMAP). Taken together, these data support a framework for resilience in which neuronal antioxidant pathway activity provides for stability of synapses within the hippocampus.

Heterogeneous effects on type 2 diabetes and cardiovascular outcomes of genetic variants and traits associated with fasting insulin.

Hyperinsulinemia is a complex and heterogeneous phenotype that characterizes molecular alterations that precede the development of type 2 diabetes (T2D). It results from a complex combination of molecular processes, including insulin secretion and insulin sensitivity, that differ between individuals. To better understand the physiology of hyperinsulinemia and ultimately T2D, we implemented a genetic approach grouping fasting insulin (FI)-associated genetic variants based on their molecular and phenotypic similarities. We identified seven distinctive genetic clusters representing different physiologic mechanisms leading to rising FI levels, ranging from clusters of variants with effects on increased FI, but without increased risk of T2D (non-diabetogenic hyperinsulinemia), to clusters of variants that increase FI and T2D risk with demonstrated strong effects on body fat distribution, liver, lipid, and inflammatory processes (diabetogenic hyperinsulinemia). We generated cluster-specific polygenic scores in 1,104,258 individuals from five multi-ancestry cohorts to show that the clusters differed in associations with cardiometabolic traits. Among clusters characterized by non-diabetogenic hyperinsulinemia, there was both increased and decreased risk of coronary artery disease despite the non-increased risk of T2D. Similarly, the clusters characterized by diabetogenic hyperinsulinemia were associated with an increased risk of T2D, yet had differing risks of cardiovascular conditions, including coronary artery disease, myocardial infarction, and stroke. The strongest cluster-T2D associations were observed with the same direction of effect in non-Hispanic Black, Hispanic, non-Hispanic White, and non-Hispanic East Asian populations. These genetic clusters provide important insights into granular metabolic processes underlying the physiology of hyperinsulinemia, notably highlighting specific processes that decouple increasing FI levels from T2D and cardiovascular risk. Our findings suggest that increasing FI levels are not invariably associated with adverse cardiometabolic outcomes.

The Microbiome, Epigenome, and Diet in Adults with Obesity during Behavioral Weight Loss.

Obesity has been linked to the gut microbiome, epigenome, and diet, yet these factors have not been studied together during obesity treatment. Our objective was to evaluate associations among gut microbiota (MB), DNA methylation (DNAme), and diet prior to and during a behavioral weight loss intervention. Adults (n = 47, age 40.9 ± 9.7 years, body mass index (BMI) 33.5 ± 4.5 kg/m2, 77% female) with data collected at baseline (BL) and 3 months (3 m) were included. Fecal MB was assessed via 16S sequencing and whole blood DNAme via the Infinium EPIC array. Food group and nutrient intakes and Healthy Eating Index (HEI) scores were calculated from 7-day diet records. Linear models were used to test for the effect of taxa relative abundance on DNAme and diet cross-sectionally at each time point, adjusting for confounders and a false discovery rate of 5%. Mean weight loss was 6.2 ± 3.9% at 3 m. At BL, one MB taxon, Ruminiclostridium, was associated with DNAme of the genes COL20A1 (r = 0.651, p = 0.029), COL18A1 (r = 0.578, p = 0.044), and NT5E (r = 0.365, p = 0.043). At 3 m, there were 14 unique MB:DNAme associations, such as Akkermansia with DNAme of GUSB (r = -0.585, p = 0.003), CRYL1 (r = -0.419, p = 0.007), C9 (r = -0.439, p = 0.019), and GMDS (r = -0.559, p = 0.046). Among taxa associated with DNAme, no significant relationships were seen with dietary intakes of relevant nutrients, food groups, or HEI scores. Our findings indicate that microbes linked to mucin degradation, short-chain fatty acid production, and body weight are associated with DNAme of phenotypically relevant genes. These relationships offer an initial understanding of the possible routes by which alterations in gut MB may influence metabolism during weight loss.

Mendelian randomization study of diabetes and dementia in the Million Veteran Program.

INTRODUCTION: Diabetes and dementia are diseases of high health-care burden worldwide. Individuals with diabetes have 1.4 to 2.2 times higher risk of dementia. Our objective was to evaluate evidence of causality between these two common diseases. METHODS: We conducted a one-sample Mendelian randomization (MR) analysis in the US Department of Veterans Affairs Million Veteran program. The study included 334,672 participants ≥65 years of age with type 2 diabetes and dementia case-control status and genotype data. RESULTS: For each standard deviation increase in genetically predicted diabetes, we found increased odds of three dementia diagnoses in non-Hispanic White participants (all-cause: odds ratio [OR] = 1.07 [1.05-1.08], P = 3.40E-18; vascular: OR = 1.11 [1.07-1.15], P = 3.63E-09, Alzheimer's disease [AD]: OR = 1.06 [1.02-1.09], P = 6.84E-04) and non-Hispanic Black participants (all-cause: OR = 1.06 [1.02-1.10], P = 3.66E-03, vascular: OR = 1.11 [1.04-1.19], P = 2.20E-03, AD: OR = 1.12 [1.02-1.23], P = 1.60E-02) but not in Hispanic participants (all P > 0.05). DISCUSSION: We found evidence of causality between diabetes and dementia using a one-sample MR study, with access to individual level data, overcoming limitations of prior studies using two-sample MR techniques.

Large scale proteomic studies create novel privacy considerations.

Privacy protection is a core principle of genomic but not proteomic research. We identified independent single nucleotide polymorphism (SNP) quantitative trait loci (pQTL) from COPDGene and Jackson Heart Study (JHS), calculated continuous protein level genotype probabilities, and then applied a naïve Bayesian approach to link SomaScan 1.3K proteomes to genomes for 2812 independent subjects from COPDGene, JHS, SubPopulations and InteRmediate Outcome Measures In COPD Study (SPIROMICS) and Multi-Ethnic Study of Atherosclerosis (MESA). We correctly linked 90-95% of proteomes to their correct genome and for 95-99% we identify the 1% most likely links. The linking accuracy in subjects with African ancestry was lower (~ 60%) unless training included diverse subjects. With larger profiling (SomaScan 5K) in the Atherosclerosis Risk Communities (ARIC) correct identification was > 99% even in mixed ancestry populations. We also linked proteomes-to-proteomes and used the proteome only to determine features such as sex, ancestry, and first-degree relatives. When serial proteomes are available, the linking algorithm can be used to identify and correct mislabeled samples. This work also demonstrates the importance of including diverse populations in omics research and that large proteomic datasets (> 1000 proteins) can be accurately linked to a specific genome through pQTL knowledge and should not be considered unidentifiable.

Propensity scores as a novel method to guide sample allocation and minimize batch effects during the design of high throughput experiments.

BACKGROUND: We developed a novel approach to minimize batch effects when assigning samples to batches. Our algorithm selects a batch allocation, among all possible ways of assigning samples to batches, that minimizes differences in average propensity score between batches. This strategy was compared to randomization and stratified randomization in a case-control study (30 per group) with a covariate (case vs control, represented as ß1, set to be null) and two biologically relevant confounding variables (age, represented as ß2, and hemoglobin A1c (HbA1c), represented as ß3). Gene expression values were obtained from a publicly available dataset of expression data obtained from pancreas islet cells. Batch effects were simulated as twice the median biological variation across the gene expression dataset and were added to the publicly available dataset to simulate a batch effect condition. Bias was calculated as the absolute difference between observed betas under the batch allocation strategies and the true beta (no batch effects). Bias was also evaluated after adjustment for batch effects using ComBat as well as a linear regression model. In order to understand performance of our optimal allocation strategy under the alternative hypothesis, we also evaluated bias at a single gene associated with both age and HbA1c levels in the 'true' dataset (CAPN13 gene). RESULTS: Pre-batch correction, under the null hypothesis (ß1), maximum absolute bias and root mean square (RMS) of maximum absolute bias, were minimized using the optimal allocation strategy. Under the alternative hypothesis (ß2 and ß3 for the CAPN13 gene), maximum absolute bias and RMS of maximum absolute bias were also consistently lower using the optimal allocation strategy. ComBat and the regression batch adjustment methods performed well as the bias estimates moved towards the true values in all conditions under both the null and alternative hypotheses. Although the differences between methods were less pronounced following batch correction, estimates of bias (average and RMS) were consistently lower using the optimal allocation strategy under both the null and alternative hypotheses. CONCLUSIONS: Our algorithm provides an extremely flexible and effective method for assigning samples to batches by exploiting knowledge of covariates prior to sample allocation.

Mendelian randomization study of diabetes and dementia in the Million Veteran Program.

INTRODUCTION: Diabetes and dementia are diseases of high healthcare burden worldwide. Individuals with diabetes have 1.4 to 2.2 times higher risk of dementia. Our objective was to evaluate evidence of causality between these two common diseases. METHODS: We conducted a one-sample Mendelian randomization (MR) analysis in the U.S. Department of Veterans Affairs Million Veteran program. The study included 334,672 participants ≥65 years of age with type 2 diabetes and dementia case-control status and genotype data. RESULTS: For each standard deviation increase in genetically-predicted diabetes, we found increased odds of three dementia diagnoses in non-Hispanic White participants (all-cause: OR=1.07[1.05-1.08], P =3.40E-18; vascular: OR=1.11[1.07-1.15], P =3.63E-09, Alzheimer's: OR=1.06[1.02-1.09], P =6.84E-04) and non-Hispanic Black participants (all-cause: OR=1.06[1.02-1.10], P =3.66E-03, vascular: OR=1.11[1.04-1.19], P =2.20E-03, Alzheimer's: OR=1.12 [1.02-1.23], P =1.60E-02) but not in Hispanic participants (all P >.05). DISCUSSION: We found evidence of causality between diabetes and dementia using a one-sample MR study, with access to individual level data, overcoming limitations of prior studies utilizing two-sample MR techniques.

NetSHy: network summarization via a hybrid approach leveraging topological properties.

MOTIVATION: Biological networks can provide a system-level understanding of underlying processes. In many contexts, networks have a high degree of modularity, i.e. they consist of subsets of nodes, often known as subnetworks or modules, which are highly interconnected and may perform separate functions. In order to perform subsequent analyses to investigate the association between the identified module and a variable of interest, a module summarization, that best explains the module's information and reduces dimensionality is often needed. Conventional approaches for obtaining network representation typically rely only on the profiles of the nodes within the network while disregarding the inherent network topological information. RESULTS: In this article, we propose NetSHy, a hybrid approach which is capable of reducing the dimension of a network while incorporating topological properties to aid the interpretation of the downstream analyses. In particular, NetSHy applies principal component analysis (PCA) on a combination of the node profiles and the well-known Laplacian matrix derived directly from the network similarity matrix to extract a summarization at a subject level. Simulation scenarios based on random and empirical networks at varying network sizes and sparsity levels show that NetSHy outperforms the conventional PCA approach applied directly on node profiles, in terms of recovering the true correlation with a phenotype of interest and maintaining a higher amount of explained variation in the data when networks are relatively sparse. The robustness of NetSHy is also demonstrated by a more consistent correlation with the observed phenotype as the sample size decreases. Lastly, a genome-wide association study is performed as an application of a downstream analysis, where NetSHy summarization scores on the biological networks identify more significant single nucleotide polymorphisms than the conventional network representation. AVAILABILITY AND IMPLEMENTATION: R code implementation of NetSHy is available at https://github.com/thaovu1/NetSHy. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

African ancestry GWAS of dementia in a large military cohort identifies significant risk loci.

While genome wide association studies (GWASs) of Alzheimer's Disease (AD) in European (EUR) ancestry cohorts have identified approximately 83 potentially independent AD risk loci, progress in non-European populations has lagged. In this study, data from the Million Veteran Program (MVP), a biobank which includes genetic data from more than 650,000 US Veteran participants, was used to examine dementia genetics in an African descent (AFR) cohort. A GWAS of Alzheimer's disease and related dementias (ADRD), an expanded AD phenotype including dementias such as vascular and non-specific dementia that included 4012 cases and 18,435 controls age 60+ in AFR MVP participants was performed. A proxy dementia GWAS based on survey-reported parental AD or dementia (n = 4385 maternal cases, 2256 paternal cases, and 45,970 controls) was also performed. These two GWASs were meta-analyzed, and then subsequently compared and meta-analyzed with the results from a previous AFR AD GWAS from the Alzheimer's Disease Genetics Consortium (ADGC). A meta-analysis of common variants across the MVP ADRD and proxy GWASs yielded GWAS significant associations in the region of APOE (p = 2.48 × 10-101), in ROBO1 (rs11919682, p = 1.63 × 10-8), and RNA RP11-340A13.2 (rs148433063, p = 8.56 × 10-9). The MVP/ADGC meta-analysis yielded additional significant SNPs near known AD risk genes TREM2 (rs73427293, p = 2.95 × 10-9), CD2AP (rs7738720, p = 1.14 × 10-9), and ABCA7 (rs73505251, p = 3.26 × 10-10), although the peak variants observed in these genes differed from those previously reported in EUR and AFR cohorts. Of the genes in or near suggestive or genome-wide significant associated variants, nine (CDA, SH2D5, DCBLD1, EML6, GOPC, ABCA7, ROS1, TMCO4, and TREM2) were differentially expressed in the brains of AD cases and controls. This represents the largest AFR GWAS of AD and dementia, finding non-APOE GWAS-significant common SNPs associated with dementia. Increasing representation of AFR participants is an important priority in genetic studies and may lead to increased insight into AD pathophysiology and reduce health disparities.

A Diabetes Genetic Risk Score Is Associated With All-Cause Dementia and Clinically Diagnosed Vascular Dementia in the Million Veteran Program.

OBJECTIVE: Diabetes and dementia are diseases of high health care burden worldwide, and studies have shown that diabetes is associated with an increased relative risk of dementia. We set out to examine whether type 2 diabetes-associated genetic variants were associated with dementia and whether they differed by race/ethnicity or clinical dementia diagnosis. RESEARCH DESIGN AND METHODS: We evaluated associations of two type 2 diabetes genetic risk scores (GRS and GRS-nonAPOE: a score without rs429358, a variant associated with Alzheimer disease [AD]) with three classifications of clinical dementia diagnoses in the Million Veteran Program (MVP): all-cause dementia, vascular dementia (VaD), and AD. We conducted our analysis stratified by European (EUR), African (AFR), and Hispanic (HIS) races/ethnicities. RESULTS: In EUR, we found associations of the GRS with all-cause dementia (odds ratio [OR] 1.06, P = 1.60e-07) and clinically diagnosed VaD (OR 1.12, P = 5.2e-05) but not with clinically diagnosed AD (OR 1.02, P = 0.43). The GRS was not associated with any dementia outcome in AFR or HIS. When testing with GRS-nonAPOE, we found that effect size estimates in EUR increased and P values decreased for all-cause dementia (OR 1.08, P = 2.6e-12), for VaD (OR 1.14, P = 7.2e-07), and for AD (OR 1.06, P = 0.018). For AFR, the association of GRS-nonAPOE and clinically diagnosed VaD (OR 1.15, P = 0.016) was statistically significant. There were no significant findings for HIS. CONCLUSIONS: We found evidence suggesting shared genetic pathogenesis of diabetes with all-cause dementia and clinically diagnosed VaD.

A multi-population phenome-wide association study of genetically-predicted height in the Million Veteran Program.

BACKGROUND: Height has been associated with many clinical traits but whether such associations are causal versus secondary to confounding remains unclear in many cases. To systematically examine this question, we performed a Mendelian Randomization-Phenome-wide association study (MR-PheWAS) using clinical and genetic data from a national healthcare system biobank. METHODS AND FINDINGS: Analyses were performed using data from the US Veterans Affairs (VA) Million Veteran Program in non-Hispanic White (EA, n = 222,300) and non-Hispanic Black (AA, n = 58,151) adults in the US. We estimated height genetic risk based on 3290 height-associated variants from a recent European-ancestry genome-wide meta-analysis. We compared associations of measured and genetically-predicted height with phenome-wide traits derived from the VA electronic health record, adjusting for age, sex, and genetic principal components. We found 345 clinical traits associated with measured height in EA and an additional 17 in AA. Of these, 127 were associated with genetically-predicted height at phenome-wide significance in EA and 2 in AA. These associations were largely independent from body mass index. We confirmed several previously described MR associations between height and cardiovascular disease traits such as hypertension, hyperlipidemia, coronary heart disease (CHD), and atrial fibrillation, and further uncovered MR associations with venous circulatory disorders and peripheral neuropathy in the presence and absence of diabetes. As a number of traits associated with genetically-predicted height frequently co-occur with CHD, we evaluated effect modification by CHD status of genetically-predicted height associations with risk factors for and complications of CHD. We found modification of effects of MR associations by CHD status for atrial fibrillation/flutter but not for hypertension, hyperlipidemia, or venous circulatory disorders. CONCLUSIONS: We conclude that height may be an unrecognized but biologically plausible risk factor for several common conditions in adults. However, more studies are needed to reliably exclude horizontal pleiotropy as a driving force behind at least some of the MR associations observed in this study.

Genetic Interactions with Intrauterine Diabetes Exposure in Relation to Obesity: The EPOCH and Project Viva Studies.

To examine whether BMI-associated genetic risk variants modify the association of intrauterine diabetes exposure with childhood BMI z-scores, we assessed the interaction between 95 BMI-associated genetic variants and in utero exposure to maternal diabetes among 459 children in the Exploring Perinatal Outcomes among Children historical prospective cohort study (n = 86 exposed; 373 unexposed) in relation to age- and sex-standardized childhood BMI z-scores (mean age = 10.3 years, standard deviation = 1.5 years). For the genetic variants showing a nominally significant interaction, we assessed the relationship in an additional 621 children in Project Viva, which is an independent longitudinal cohort study, and used meta-analysis to combine the results for the two studies. Seven of the ninety-five genetic variants tested exhibited a nominally significant interaction with in utero exposure to maternal diabetes in relation to the offspring BMI z-score in EPOCH. Five of the seven variants exhibited a consistent direction of interaction effect across both EPOCH and Project Viva. While none achieved statistical significance in the meta-analysis after accounting for multiple testing, three variants exhibited a nominally significant interaction with in utero exposure to maternal diabetes in relation to offspring BMI z-score: rs10733682 near LMX1B (interaction ß = 0.39; standard error (SE) = 0.17), rs17001654 near SCARB2 (ß = 0.53; SE = 0.22), and rs16951275 near MAP2K5 (ß = 0.37; SE = 0.17). BMI-associated genetic variants may enhance the association between exposure to in utero diabetes and higher childhood BMI, but larger studies of in utero exposures are necessary to confirm the observed nominally significant relationships.

Genetic Risk Score for Type 2 Diabetes and Traits Related to Glucose-Insulin Homeostasis in Youth: The Exploring Perinatal Outcomes Among Children (EPOCH) Study.

OBJECTIVE: The metabolic phenotype of youth-onset type 2 diabetes (T2D) differs from that of adult-onset T2D, but little is known about genetic contributions. We aimed to evaluate the association between a T2D genetic risk score (GRS) and traits related to glucose-insulin homeostasis among healthy youth. RESEARCH DESIGN AND METHODS: We used data from 356 youth (mean age 16.7 years; 50% female) in the Exploring Perinatal Outcomes Among Children (EPOCH) cohort to calculate a standardized weighted GRS based on 271 single nucleotide polymorphisms associated with T2D in adults. We used linear regression to assess associations of the GRS with log-transformed fasting glucose, 2-h glucose, HOMA of insulin resistance (HOMA-IR), oral disposition index, and insulinogenic index adjusted for age, sex, BMI z score, in utero exposure to maternal diabetes, and genetic principal components. We also evaluated effect modification by BMI z score, in utero exposure to maternal diabetes, and ethnicity. RESULTS: Higher weighted GRS was associated with lower oral disposition index (ß = -0.11; 95% CI -0.19, -0.02) and insulinogenic index (ß = -0.08; 95% CI -0.17, -0.001), but not with fasting glucose (ß = 0.01; 95% CI -0.01, 0.02), 2-h glucose (ß = 0.03; 95% CI -0.0004, 0.06), or HOMA-IR (ß = 0.02; 95% CI -0.04, 0.07). BMI z score and in utero exposure to maternal diabetes increased the effect of the GRS on glucose levels. CONCLUSIONS: Our results suggest that T2D genetic risk factors established in adults are relevant to glucose-insulin homeostasis in youth and that maintaining a healthy weight may be particularly important for youth with high genetic risk of T2D.

Comparing Statistical Tests for Differential Network Analysis of Gene Modules.

Genes often work together to perform complex biological processes, and "networks" provide a versatile framework for representing the interactions between multiple genes. Differential network analysis (DiNA) quantifies how this network structure differs between two or more groups/phenotypes (e.g., disease subjects and healthy controls), with the goal of determining whether differences in network structure can help explain differences between phenotypes. In this paper, we focus on gene co-expression networks, although in principle, the methods studied can be used for DiNA for other types of features (e.g., metabolome, epigenome, microbiome, proteome, etc.). Three common applications of DiNA involve (1) testing whether the connections to a single gene differ between groups, (2) testing whether the connection between a pair of genes differs between groups, or (3) testing whether the connections within a "module" (a subset of 3 or more genes) differs between groups. This article focuses on the latter, as there is a lack of studies comparing statistical methods for identifying differentially co-expressed modules (DCMs). Through extensive simulations, we compare several previously proposed test statistics and a new p-norm difference test (PND). We demonstrate that the true positive rate of the proposed PND test is competitive with and often higher than the other methods, while controlling the false positive rate. The R package discoMod (differentially co-expressed modules) implements the proposed method and provides a full pipeline for identifying DCMs: clustering tools to derive gene modules, tests to identify DCMs, and methods for visualizing the results.

Inverse probability weighting is an effective method to address selection bias during the analysis of high dimensional data.

Omics studies frequently use samples collected during cohort studies. Conditioning on sample availability can cause selection bias if sample availability is nonrandom. Inverse probability weighting (IPW) is purported to reduce this bias. We evaluated IPW in an epigenome-wide analysis testing the association between DNA methylation (261,435 probes) and age in healthy adolescent subjects (n = 114). We simulated age and sex to be correlated with sample selection and then evaluated four conditions: complete population/no selection bias (all subjects), naïve selection bias (no adjustment), and IPW selection bias (selection bias with IPW adjustment). Assuming the complete population condition represented the "truth," we compared each condition to the complete population condition. Bias or difference in associations between age and methylation was reduced in the IPW condition versus the naïve condition. However, genomic inflation and type 1 error were higher in the IPW condition relative to the naïve condition. Postadjustment using bacon, type 1 error and inflation were similar across all conditions. Power was higher under the IPW condition compared with the naïve condition before and after inflation adjustment. IPW methods can reduce bias in genome-wide analyses. Genomic inflation is a potential concern that can be minimized using methods that adjust for inflation.

Genetic Risk for Hepatic Fat among an Ethnically Diverse Cohort of Youth: The Exploring Perinatal Outcomes among Children Study.

OBJECTIVE: To assess the importance of genetic and nongenetic risk factors contributing to hepatic fat accumulation in a multiethnic population of youth. STUDY DESIGN: We investigated the relationship between genetic factors and hepatic fat fraction (HFF) in 347 children aged 12.5-19.5 years. We examined 5 single nucleotide polymorphisms previously associated with HFF and a weighted genetic risk score (GRS) and examined how these associations varied with ethnicity (Hispanic vs non-Hispanic white) and body mass index (BMI) category. We also compared how much variation in HFF was explained by genetic factors vs cardiometabolic factors (BMI z-score and the Homeostasis Model of Insulin Resistance) or diet. RESULTS: PNPLA3 rs738409 and the GRS were each associated with HFF among Hispanic (ß = 0.39; 95% CI, 0.16-0.62; P = .001; and ß = 0.20; 95% CI, 0.05-0.34; P = .007, respectively) but not non-Hispanic white (ß = 0.04; 95% CI, -0.18 to 0.26; P = .696; and ß = 0.03; 95% CI, -0.09 to 0.14; P = .651, respectively) youth. Cardiometabolic risk factors explained more of the variation in HFF than genetic risk factors among non-lean Hispanic individuals (27.2% for cardiometabolic markers vs 6.4% for rs738409 and 4.3% for the GRS), and genetic risk factors were more important among lean individuals (2.7% for cardiometabolic markers vs 12.6% for rs738409 and 4.4% for the GRS). CONCLUSIONS: Poor cardiometabolic health may be more important than genetic factors when predicting HFF in overweight and obese young populations. Genetic risk is an important contributor to pediatric HFF among lean Hispanics, but further studies are necessary to elucidate the strength of the association between genetic risk and HFF in non-Hispanic white youth.