Chuan Gao

HEFT: eQTL analysis of many thousands of expressed genes while simultaneously controlling for hidden factors

MOTIVATION Identification of expression Quantitative Trait Loci (eQTL), the genetic loci that contribute to heritable variation in gene expression, can be obstructed by factors that produce variation in expression profiles if these factors are unmeasured or hidden from direct analysis. METHODS We have developed a method for Hidden Expression Factor analysis (HEFT) that identifies individual and pleiotropic effects of eQTL in the presence of hidden factors. The HEFT model is a combined multivariate regression and factor analysis, where the complete likelihood of the model is used to derive a ridge estimator for simultaneous factor learning and detection of eQTL. HEFT requires no pre-estimation of hidden factor effects; it provides P-values and is extremely fast, requiring just a few hours to complete an eQTL analysis of thousands of expression variables when analyzing hundreds of thousands of single nucleotide polymorphisms on a standard 8 core 2.6 G desktop. RESULTS By analyzing simulated data, we demonstrate that HEFT can correct for an unknown number of hidden factors and significantly outperforms all related hidden factor methods for eQTL analysis when there are eQTL with univariate and multivariate (pleiotropic) effects. To demonstrate a real-world application, we applied HEFT to identify eQTL affecting gene expression in the human lung for a study that included presumptive hidden factors. HEFT identified all of the cis-eQTL found by other hidden factor methods and 91 additional cis-eQTL. HEFT also identified a number of eQTLs with direct relevance to lung disease that could not be found without a hidden factor analysis, including cis-eQTL for GTF2H1 and MTRR, genes that have been independently associated with lung cancer. AVAILABILITY Software is available at http://mezeylab.cb.bscb.cornell.edu/Software.aspx. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

A Comprehensive Analysis of Common and Rare Variants to Identify Adiposity Loci in Hispanic Americans: The IRAS Family Study (IRASFS)

Obesity is growing epidemic affecting 35% of adults in the United States. Previous genome-wide association studies (GWAS) have identified numerous loci associated with obesity. However, the majority of studies have been completed in Caucasians focusing on total body measures of adiposity. Here we report the results from genome-wide and exome chip association studies focusing on total body measures of adiposity including body mass index (BMI), percent body fat (PBF) and measures of fat deposition including waist circumference (WAIST), waist-hip ratio (WHR), subcutaneous adipose tissue (SAT), and visceral adipose tissue (VAT) in Hispanic Americans (nmax = 1263) from the Insulin Resistance Atherosclerosis Family Study (IRASFS). Five SNPs from two novel loci attained genome-wide significance (P<5.00x10-8) in IRASFS. A missense SNP in the isocitrate dehydrogenase 1 gene (IDH1) was associated with WAIST (rs34218846, MAF = 6.8%, PDOM = 1.62x10-8). This protein is postulated to play an important role in fat and cholesterol biosynthesis as demonstrated in cell and knock-out animal models. Four correlated intronic SNPs in the Zinc finger, GRF-type containing 1 gene (ZGRF1; SNP rs1471880, MAF = 48.1%, PDOM = 1.00x10-8) were strongly associated with WHR. The exact biological function of ZGRF1 and the connection with adiposity remains unclear. SNPs with p-values less than 5.00x10-6 from IRASFS were selected for replication. Meta-analysis was computed across seven independent Hispanic-American cohorts (nmax = 4156) and the strongest signal was rs1471880 (PDOM = 8.38x10-6) in ZGRF1 with WAIST. In conclusion, a genome-wide and exome chip association study was conducted that identified two novel loci (IDH1 and ZGRF1) associated with adiposity. While replication efforts were inconclusive, when taken together with the known biology, IDH1 and ZGRF1 warrant further evaluation.

Vitamin D-responsive SGPP2 variants associated with lung cell expression and lung function

BackgroundVitamin D is associated with lung health in epidemiologic studies, but mechanisms mediating observed associations are poorly understood. This study explores mechanisms for an effect of vitamin D in lung through an in vivo gene expression study, an expression quantitative trait loci (eQTL) analysis in lung tissue, and a population-based cohort study of sequence variants.MethodsMicroarray analysis investigated the association of gene expression in small airway epithelial cells with serum 25(OH)D in adult non-smokers. Sequence variants in candidate genes identified by the microarray were investigated in a lung tissue eQTL database, and also in relation to cross-sectional pulmonary function in the Health, Aging, and Body Composition (Health ABC) study, stratified by race, with replication in the Framingham Heart Study (FHS).Results13 candidate genes had significant differences in expression by serum 25(OH)D (nominal p < 0.05), and a genome-wide significant eQTL association was detected for SGPP2. In Health ABC, SGPP2 SNPs were associated with FEV1 in both European- and African-Americans, and the gene-level association was replicated in European-American FHS participants. SNPs in 5 additional candidate genes (DAPK1, FSTL1, KAL1, KCNS3, and RSAD2) were associated with FEV1 in Health ABC participants.ConclusionsSGPP2, a sphingosine-1-phosphate phosphatase, is a novel vitamin D-responsive gene associated with lung function. The identified associations will need to be followed up in further studies.

Association Analysis of the Cubilin (CUBN) and Megalin (LRP2) Genes with ESRD in African Americans

BACKGROUND AND OBJECTIVES Genetic variation in the cubilin (CUBN) gene is associated with albuminuria and CKD. Common and rare coding variants in CUBN and the gene encoding its transport partner megalin (LRP2) were assessed for association with ESRD in blacks. DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS Sixty-six CUBN and LRP2 single-nucleotide polymorphisms (SNPs) were selected and analyzed in this multistage study. Exome sequencing data from 529 blacks with type 2 diabetes (T2D) -associated ESRD and 535 controls lacking T2D or nephropathy (the Type 2 Diabetes Genes [T2D-GENES] Consortium) were first evaluated, focusing on coding variants in CUBN and LRP2; 15 potentially associated SNPs identified from the T2D-GENES Consortium as well as 51 other selected SNPs were then assessed in an independent T2D-ESRD sample set of blacks (the Affymetrix Axiom Biobank Genotyping Array [AXIOM]; 2041 patients with T2D-ESRD, 627 patients with T2D without nephropathy, and 1140 nondiabetic, non-nephropathy controls). A meta-analysis combining the T2D-GENES Consortium and the AXIOM data was performed for 18 overlapping SNPs. Additionally, all 66 SNPs were genotyped in the Wake Forest School of Medicine samples of blacks with nondiabetic ESRD (885 patients with nondiabetic ESRD and 721 controls). Association testing with ESRD was performed in models including age, sex, African ancestry proportion, and apolipoprotein L1 gene renal-risk variants. RESULTS CUBN SNP rs1801239 (I2984V), previously associated with albuminuria, was significantly associated with T2D-ESRD in blacks (the T2D-GENES Consortium and the AXIOM meta-analysis, P=0.03; odds ratio, 1.31; 95% confidence interval, 1.03 to 1.67; minor allele frequency =0.028). A novel LRP2 missense variant, rs17848169 (N2632D), was also significantly protective from T2D-ESRD (the T2D-GENES Consortium and the AXIOM, P<0.002; odds ratio, 0.47; 95% confidence interval, 0.29 to 0.75; meta-analysis minor allele frequency =0.007). Neither SNP was associated with T2D when contrasting patients with T2D with controls lacking diabetes. CUBN and LRP2 SNPs were not associated with nondiabetic etiologies of ESRD. CONCLUSIONS Evidence for genetic association exists between a cubilin and a rare megalin variant with diabetes-associated ESRD in populations with recent African ancestry.

Genome-Wide Study of Subcutaneous and Visceral Adipose Tissue Reveals Novel Sex-Specific Adiposity Loci in Mexican Americans

This study aimed to explore the genetic mechanisms of regional fat deposition, which is a strong risk factor for metabolic diseases beyond total adiposity.

Differential Expression of Vitamin E and Selenium-Responsive Genes by Disease Severity in Chronic Obstructive Pulmonary Disease

Abstract Antioxidant nutritional status is hypothesized to influence chronic obstructive pulmonary disease (COPD) susceptibility and progression. Although past studies relate antioxidants to gene expression, there are no data in patients with COPD. This study investigated the hypothesis that antioxidant status is compromised in patients with COPD, and antioxidant-responsive genes differentially express in a similar pattern. Lung tissue samples from patients with COPD were assayed for vitamin E and gene expression. Selenium and vitamin E were assayed in corresponding plasma samples. Discovery based genome-wide expression analysis compared moderate, severe, and very severe COPD (GOLD II-IV) patients to mild and at-risk/normal (GOLD 0-I). Hypotheses-driven analyses assessed differential gene expression by disease severity for vitamin E-responsive and selenium-responsive genes. GOLD II-IV COPD patients had 30% lower lung tissue vitamin E levels compared to GOLD 0-I participants (p = 0.0082). No statistically significant genome-wide differences in expression by disease severity were identified. Hypothesis-driven analyses of 109 genes found 16 genes differentially expressed (padjusted < 0.05) by disease severity including 6 selenium-responsive genes (range in fold-change -1.39 to 2.25), 6 vitamin E-responsive genes (fold-change -2.30 to 1.51), and 4 COPD-associated genes. Lung tissue vitamin E in patients with COPD was associated with disease severity and vitamin E-responsive genes were differentially expressed by disease severity. Although nutritional status is hypothesized to contribute to COPD risk, and is of therapeutic interest, evidence to date is mainly observational. The findings reported herein are novel, and support a role of vitamin E in COPD progression.

A genome-wide linkage and association analysis of imputed insertions and deletions with cardiometabolic phenotypes in Mexican Americans: The Insulin Resistance Atherosclerosis Family Study

Insertions and deletions (INDELs) represent a significant fraction of interindividual variation in the human genome yet their contribution to phenotypes is poorly understood. To confirm the quality of imputed INDELs and investigate their roles in mediating cardiometabolic phenotypes, genome‐wide association and linkage analyses were performed for 15 phenotypes with 1,273,952 imputed INDELs in 1,024 Mexican‐origin Americans. Imputation quality was validated using whole exome sequencing with an average kappa of 0.93 in common INDELs (minor allele frequencies [MAFs] ≥ 5%). Association analysis revealed one genome‐wide significant association signal for the cholesterylester transfer protein gene (CETP) with high‐density lipoprotein levels (rs36229491, P = 3.06 × 10−12); linkage analysis identified two peaks with logarithm of the odds (LOD) > 5 (rs60560566, LOD = 5.36 with insulin sensitivity (SI) and rs5825825, LOD = 5.11 with adiponectin levels). Suggestive overlapping signals between linkage and association were observed: rs59849892 in the WSC domain containing 2 gene (WSCD2) was associated and nominally linked with SI (P = 1.17 × 10−7, LOD = 1.99). This gene has been implicated in glucose metabolism in human islet cell expression studies. In addition, rs201606363 was linked and nominally associated with low‐density lipoprotein (P = 4.73 × 10−4, LOD = 3.67), apolipoprotein B (P = 1.39 × 10−3, LOD = 4.64), and total cholesterol (P = 1.35 × 10−2, LOD = 3.80) levels. rs201606363 is an intronic variant of the UBE2F‐SCLY (where UBE2F is ubiquitin‐conjugating enzyme E2F and SCLY is selenocysteine lyase) fusion gene that may regulate cholesterol through selenium metabolism. In conclusion, these results confirm the feasibility of imputing INDELs from array‐based single nucleotide polymorphism (SNP) genotypes. Analysis of these variants using association and linkage replicated previously identified SNP signals and identified multiple novel INDEL signals. These results support the inclusion of INDELs into genetic studies to more fully interrogate the spectrum of genetic variation.

Exome Sequencing Identifies Genetic Variants Associated with Circulating Lipid Levels in Mexican Americans: The Insulin Resistance Atherosclerosis Family Study (IRASFS)

Genome-wide association studies have identified numerous variants associated with lipid levels; yet, the majority are located in non-coding regions with unclear mechanisms. In the Insulin Resistance Atherosclerosis Family Study (IRASFS), heritability estimates suggest a strong genetic basis: low-density lipoprotein (LDL, h2 = 0.50), high-density lipoprotein (HDL, h2 = 0.57), total cholesterol (TC, h2 = 0.53), and triglyceride (TG, h2 = 0.42) levels. Exome sequencing of 1,205 Mexican Americans (90 pedigrees) from the IRASFS identified 548,889 variants and association and linkage analyses with lipid levels were performed. One genome-wide significant signal was detected in APOA5 with TG (rs651821, PTG = 3.67 × 10−10, LODTG = 2.36, MAF = 14.2%). In addition, two correlated SNPs (r2 = 1.0) rs189547099 (PTG = 6.31 × 10−08, LODTG = 3.13, MAF = 0.50%) and chr4:157997598 (PTG = 6.31 × 10−08, LODTG = 3.13, MAF = 0.50%) reached exome-wide significance (P < 9.11 × 10−08). rs189547099 is an intronic SNP in FNIP2 and SNP chr4:157997598 is intronic in GLRB. Linkage analysis revealed 46 SNPs with a LOD > 3 with the strongest signal at rs1141070 (LODLDL = 4.30, PLDL = 0.33, MAF = 21.6%) in DFFB. A total of 53 nominally associated variants (P < 5.00 × 10−05, MAF ≥ 1.0%) were selected for replication in six Mexican-American cohorts (N = 3,280). The strongest signal observed was a synonymous variant (rs1160983, PLDL = 4.44 × 10−17, MAF = 2.7%) in TOMM40. Beyond primary findings, previously reported lipid loci were fine-mapped using exome sequencing in IRASFS. These results support that exome sequencing complements and extends insights into the genetics of lipid levels.

Genome-wide interaction with the insulin secretion locus MTNR1B reveals CMIP as a novel type 2 diabetes susceptibility gene in African Americans

Although type 2 diabetes (T2D) results from metabolic defects in insulin secretion and insulin sensitivity, most of the genetic risk loci identified to date relates to insulin secretion. We reported that T2D loci influencing insulin sensitivity may be identified through interactions with insulin secretion loci, thereby leading to T2D. Here, we hypothesize that joint testing of variant main effects and interaction effects with an insulin secretion locus increases power to identify genetic interactions leading to T2D. We tested this hypothesis with an intronic MTNR1B SNP, rs10830963, which is associated with acute insulin response to glucose, a dynamic measure of insulin secretion. rs10830963 was tested for interaction and joint (main + interaction) effects with genome‐wide data in African Americans (2,452 cases and 3,772 controls) from five cohorts. Genome‐wide genotype data (Affymetrix Human Genome 6.0 array) was imputed to a 1000 Genomes Project reference panel. T2D risk was modeled using logistic regression with rs10830963 dosage, age, sex, and principal component as predictors. Joint effects were captured using the Kraft two degrees of freedom test. Genome‐wide significant (P < 5 × 10−8) interaction with MTNR1B and joint effects were detected for CMIP intronic SNP rs17197883 (Pinteraction = 1.43 × 10−8; Pjoint = 4.70 × 10−8). CMIP variants have been nominally associated with T2D, fasting glucose, and adiponectin in individuals of East Asian ancestry, with high‐density lipoprotein, and with waist‐to‐hip ratio adjusted for body mass index in Europeans. These data support the hypothesis that additional genetic factors contributing to T2D risk, including insulin sensitivity loci, can be identified through interactions with insulin secretion loci.

XX Disorder of Sex Development is associated with an insertion on chromosome 9 and downregulation of RSPO1 in dogs (Canis lupus familiaris)

Remarkable progress has been achieved in understanding the mechanisms controlling sex determination, yet the cause for many Disorders of Sex Development (DSD) remains unknown. Of particular interest is a rare XX DSD subtype in which individuals are negative for SRY, the testis determining factor on the Y chromosome, yet develop testes or ovotestes, and both of these phenotypes occur in the same family. This is a naturally occurring disorder in humans (Homo sapiens) and dogs (C. familiaris). Phenotypes in the canine XX DSD model are strikingly similar to those of the human XX DSD subtype. The purposes of this study were to identify 1) a variant associated with XX DSD in the canine model and 2) gene expression alterations in canine embryonic gonads that could be informative to causation. Using a genome wide association study (GWAS) and whole genome sequencing (WGS), we identified a variant on C. familiaris autosome 9 (CFA9) that is associated with XX DSD in the canine model and in affected purebred dogs. This is the first marker identified for inherited canine XX DSD. It lies upstream of SOX9 within the canine ortholog for the human disorder, which resides on 17q24. Inheritance of this variant indicates that XX DSD is a complex trait in which breed genetic background affects penetrance. Furthermore, the homozygous variant genotype is associated with embryonic lethality in at least one breed. Our analysis of gene expression studies (RNA-seq and PRO-seq) in embryonic gonads at risk of XX DSD from the canine model identified significant RSPO1 downregulation in comparison to XX controls, without significant upregulation of SOX9 or other known testis pathway genes. Based on these data, a novel mechanism is proposed in which molecular lesions acting upstream of RSPO1 induce epigenomic gonadal mosaicism.