Etiology Of Complex Diseases Research Articles

VariFunNet, an integrated multiscale modeling framework to study the effects of rare non-coding variants in Genome-Wide Association Studies: applied to Alzheimer's Disease.

It is a grand challenge to reveal the causal effects of DNA variants in complex phenotypes. Although statistical techniques can establish correlations between genotypes and phenotypes in Genome-Wide Association Studies (GWAS), they often fail when the variant is rare. The emerging Network-based Association Studies aim to address this shortcoming in statistical analysis, but are mainly applied to coding variations. Increasing evidences suggest that non-coding variants play critical roles in the etiology of complex diseases. However, few computational tools are available to study the effect of rare non-coding variants on phenotypes. Here we have developed a multiscale modeling variant-to-function-to-network framework VariFunNet to address these challenges. VariFunNet first predict the functional variations of molecular interactions, which result from the non-coding variants. Then we incorporate the genes associated with the functional variation into a tissue-specific gene network, and identify subnetworks that transmit the functional variation to molecular phenotypes. Finally, we quantify the functional implication of the subnetwork, and prioritize the association of the non-coding variants with the phenotype. We have applied VariFunNet to investigating the causal effect of rare non-coding variants on Alzheimer's disease (AD). Among top 21 ranked causal non-coding variants, 16 of them are directly supported by existing evidences. The remaining 5 novel variants dysregulate multiple downstream biological processes, all of which are associated with the pathology of AD. Furthermore, we propose potential new drug targets that may modulate diverse pathways responsible for AD. These findings may shed new light on discovering new biomarkers and therapies for the prevention, diagnosis, and treatment of AD. Our results suggest that multiscale modeling is a potentially powerful approach to studying causal genotype-phenotype associations.

An Efficient Test for Gene-Environment Interaction in Generalized Linear Mixed Models with Family Data.

Gene-environment (GE) interaction has important implications in the etiology of complex diseases that are caused by a combination of genetic factors and environment variables. Several authors have developed GE analysis in the context of independent subjects or longitudinal data using a gene-set. In this paper, we propose to analyze GE interaction for discrete and continuous phenotypes in family studies by incorporating the relatedness among the relatives for each family into a generalized linear mixed model (GLMM) and by using a gene-based variance component test. In addition, we deal with collinearity problems arising from linkage disequilibrium among single nucleotide polymorphisms (SNPs) by considering their coefficients as random effects under the null model estimation. We show that the best linear unbiased predictor (BLUP) of such random effects in the GLMM is equivalent to the ridge regression estimator. This equivalence provides a simple method to estimate the ridge penalty parameter in comparison to other computationally-demanding estimation approaches based on cross-validation schemes. We evaluated the proposed test using simulation studies and applied it to real data from the Baependi Heart Study consisting of 76 families. Using our approach, we identified an interaction between BMI and the Peroxisome Proliferator Activated Receptor Gamma (PPARG) gene associated with diabetes.

Embracing polygenicity: a review of methods and tools for psychiatric genetics research.

The availability of genome-wide genetic data on hundreds of thousands of people has led to an equally rapid growth in methodologies available to analyse these data. While the motivation for undertaking genome-wide association studies (GWAS) is identification of genetic markers associated with complex traits, once generated these data can be used for many other analyses. GWAS have demonstrated that complex traits exhibit a highly polygenic genetic architecture, often with shared genetic risk factors across traits. New methods to analyse data from GWAS are increasingly being used to address a diverse set of questions about the aetiology of complex traits and diseases, including psychiatric disorders. Here, we give an overview of some of these methods and present examples of how they have contributed to our understanding of psychiatric disorders. We consider: (i) estimation of the extent of genetic influence on traits, (ii) uncovering of shared genetic control between traits, (iii) predictions of genetic risk for individuals, (iv) uncovering of causal relationships between traits, (v) identifying causal single-nucleotide polymorphisms and genes or (vi) the detection of genetic heterogeneity. This classification helps organise the large number of recently developed methods, although some could be placed in more than one category. While some methods require GWAS data on individual people, others simply use GWAS summary statistics data, allowing novel well-powered analyses to be conducted at a low computational burden.

Multidimensional Genetic Analysis of Repeated Seizures in the Hybrid Mouse Diversity Panel Reveals a Novel Epileptogenesis Susceptibility Locus.

Epilepsy has many causes and comorbidities affecting as many as 4% of people in their lifetime. Both idiopathic and symptomatic epilepsies are highly heritable, but genetic factors are difficult to characterize among humans due to complex disease etiologies. Rodent genetic studies have been critical to the discovery of seizure susceptibility loci, including Kcnj10 mutations identified in both mouse and human cohorts. However, genetic analyses of epilepsy phenotypes in mice to date have been carried out as acute studies in seizure-naive animals or in Mendelian models of epilepsy, while humans with epilepsy have a history of recurrent seizures that also modify brain physiology. We have applied a repeated seizure model to a genetic reference population, following seizure susceptibility over a 36-d period. Initial differences in generalized seizure threshold among the Hybrid Mouse Diversity Panel (HMDP) were associated with a well-characterized seizure susceptibility locus found in mice: Seizure susceptibility 1. Remarkably, Szs1 influence diminished as subsequent induced seizures had diminishing latencies in certain HMDP strains. Administration of eight seizures, followed by an incubation period and an induced retest seizure, revealed novel associations within the calmodulin-binding transcription activator 1, Camta1. Using systems genetics, we have identified four candidate genes that are differentially expressed between seizure-sensitive and -resistant strains close to our novel Epileptogenesis susceptibility factor 1 (Esf1) locus that may act individually or as a coordinated response to the neuronal stress of seizures.

Systematic tissue-specific functional annotation of the human genome highlights immune-related DNA elements for late-onset Alzheimer's disease.

Continuing efforts from large international consortia have made genome-wide epigenomic and transcriptomic annotation data publicly available for a variety of cell and tissue types. However, synthesis of these datasets into effective summary metrics to characterize the functional non-coding genome remains a challenge. Here, we present GenoSkyline-Plus, an extension of our previous work through integration of an expanded set of epigenomic and transcriptomic annotations to produce high-resolution, single tissue annotations. After validating our annotations with a catalog of tissue-specific non-coding elements previously identified in the literature, we apply our method using data from 127 different cell and tissue types to present an atlas of heritability enrichment across 45 different GWAS traits. We show that broader organ system categories (e.g. immune system) increase statistical power in identifying biologically relevant tissue types for complex diseases while annotations of individual cell types (e.g. monocytes or B-cells) provide deeper insights into disease etiology. Additionally, we use our GenoSkyline-Plus annotations in an in-depth case study of late-onset Alzheimer’s disease (LOAD). Our analyses suggest a strong connection between LOAD heritability and genetic variants contained in regions of the genome functional in monocytes. Furthermore, we show that LOAD shares a similar localization of SNPs to monocyte-functional regions with Parkinson’s disease. Overall, we demonstrate that integrated genome annotations at the single tissue level provide a valuable tool for understanding the etiology of complex human diseases. Our GenoSkyline-Plus annotations are freely available at http://genocanyon.med.yale.edu/GenoSkyline.

SEQSpark: A Complete Analysis Tool for Large-Scale Rare Variant Association Studies Using Whole-Genome and Exome Sequence Data

Massively parallel sequencing technologies provide great opportunities for discovering rare susceptibility variants involved in complex disease etiology via large-scale imputation and exome and whole-genome sequence-based association studies. Due to modest effect sizes, large sample sizes of tens to hundreds of thousands of individuals are required for adequately powered studies. Current analytical tools are obsolete when it comes to handling these large datasets. To facilitate the analysis of large-scale sequence-based studies, we developed SEQSpark which implements parallel processing based on Spark to increase the speed and efficiency of performing data quality control, annotation, and association analysis. To demonstrate the versatility and speed of SEQSpark, we analyzed whole-genome sequence data from the UK10K, testing for associations with waist-to-hip ratios. The analysis, which was completed in 1.5hr, included loading data, annotation, principal component analysis, and single variant and rare variant aggregate association analysis of >9 million variants. For rare variant aggregate analysis, an exome-wide significant association (p < 2.5× 10-6) was observed with CCDC62 (SKAT-O [p = 6.89× 10-7], combined multivariate collapsing [p = 1.48× 10-6], and burden of rare variants [p = 1.48× 10-6]). SEQSpark was also used to analyze 50,000 simulated exomes and it required 1.75hr for the analysis of a quantitative trait using several rare variant aggregate association methods. Additionally, the performance of SEQSpark was compared to Variant Association Tools and PLINK/SEQ. SEQSpark was always faster and in some situations computation was reduced to a hundredth of the time. SEQSpark will empower large sequence-based epidemiological studies to quickly elucidate genetic variation involved in the etiology of complex traits.

Epigenome-wide association study of DNA methylation in panic disorder

BackgroundPanic disorder (PD) is considered to be a multifactorial disorder emerging from interactions among multiple genetic and environmental factors. To date, although genetic studies reported several susceptibility genes with PD, few of them were replicated and the pathogenesis of PD remains to be clarified. Epigenetics is considered to play an important role in etiology of complex traits and diseases, and DNA methylation is one of the major forms of epigenetic modifications. In this study, we performed an epigenome-wide association study of PD using DNA methylation arrays so as to investigate the possibility that different levels of DNA methylation might be associated with PD.MethodsThe DNA methylation levels of CpG sites across the genome were examined with genomic DNA samples (PD, N = 48, control, N = 48) extracted from peripheral blood. Methylation arrays were used for the analysis. β values, which represent the levels of DNA methylation, were normalized via an appropriate pipeline. Then, β values were converted to M values via the logit transformation for epigenome-wide association study. The relationship between each DNA methylation site and PD was assessed by linear regression analysis with adjustments for the effects of leukocyte subsets.ResultsForty CpG sites showed significant association with PD at 5% FDR correction, though the differences of the DNA methylation levels were relatively small. Most of the significant CpG sites (37/40 CpG sites) were located in or around CpG islands. Many of the significant CpG sites (27/40 CpG sites) were located upstream of genes, and all such CpG sites with the exception of two were hypomethylated in PD subjects. A pathway analysis on the genes annotated to the significant CpG sites identified several pathways, including “positive regulation of lymphocyte activation.”ConclusionsAlthough future studies with larger number of samples are necessary to confirm the small DNA methylation abnormalities associated with PD, there is a possibility that several CpG sites might be associated, together as a group, with PD.

Proceedings of the third international molecular pathological epidemiology (MPE) meeting.

Molecular pathological epidemiology (MPE) is a transdisciplinary and relatively new scientific discipline that integrates theory, methods, and resources from epidemiology, pathology, biostatistics, bioinformatics, and computational biology. The underlying objective of MPE research is to better understand the etiology and progression of complex and heterogeneous human diseases with the goal of informing prevention and treatment efforts in population health and clinical medicine. Although MPE research has been commonly applied to investigating breast, lung, and colorectal cancers, its methodology can be used to study most diseases. Recent successes in MPE studies include: (1) the development of new statistical methods to address etiologic heterogeneity; (2) the enhancement of causal inference; (3) the identification of previously unknown exposure-subtype disease associations; and (4) better understanding of the role of lifestyle/behavioral factors on modifying prognosis according to disease subtype. Central challenges to MPE include the relative lack of transdisciplinary experts, educational programs, and forums to discuss issues related to the advancement of the field. To address these challenges, highlight recent successes in the field, and identify new opportunities, a series of MPE meetings have been held at the Dana-Farber Cancer Institute in Boston, MA. Herein, we share the proceedings of the Third International MPE Meeting, held in May 2016 and attended by 150 scientists from 17 countries. Special topics included integration of MPE with immunology and health disparity research. This meeting series will continue to provide an impetus to foster further transdisciplinary integration of divergent scientific fields.

Properties of human disease genes and the role of genes linked to Mendelian disorders in complex disease aetiology.

Do genes presenting variation that has been linked to human disease have different biological properties than genes that have never been related to disease? What is the relationship between disease and fitness? Are the evolutionary pressures that affect genes linked to Mendelian diseases the same to those acting on genes whose variation contributes to complex disorders? The answers to these questions could shed light on the architecture of human genetic disorders and may have relevant implications when designing mapping strategies in future genetic studies. Here we show that, relative to non-disease genes, human disease (HD) genes have specific evolutionary profiles and protein network properties. Additionally, our results indicate that the mutation-selection balance renders an insufficient account of the evolutionary history of some HD genes and that adaptive selection could also contribute to shape their genetic architecture. Notably, several biological features of HD genes depend on the type of pathology (complex or Mendelian) with which they are related. For example, genes harbouring both causal variants for Mendelian disorders and risk factors for complex disease traits (Complex-Mendelian genes), tend to present higher functional relevance in the protein network and higher expression levels than genes associated only with complex disorders. Moreover, risk variants in Complex-Mendelian genes tend to present higher odds ratios than those on genes associated with the same complex disorders but with no link to Mendelian diseases. Taken together, our results suggest that genetic variation at genes linked to Mendelian disorders plays an important role in driving susceptibility to complex disease.

Generalized partial linear varying multi-index coefficient model for gene-environment interactions.

Epidemiological studies have suggested the joint effect of simultaneous exposures to multiple environments on disease risk. However, how environmental mixtures as a whole jointly modify genetic effect on disease risk is still largely unknown. Given the importance of gene-environment (G×E) interactions on many complex diseases, rigorously assessing the interaction effect between genes and environmental mixtures as a whole could shed novel insights into the etiology of complex diseases. For this purpose, we propose a generalized partial linear varying multi-index coefficient model (GPLVMICM) to capture the genetic effect on disease risk modulated by multiple environments as a whole. GPLVMICM is semiparametric in nature which allows different index loading parameters in different index functions. We estimate the parametric parameters by a profile procedure, and the nonparametric index functions by a B-spline backfitted kernel method. Under some regularity conditions, the proposed parametric and nonparametric estimators are shown to be consistent and asymptotically normal. We propose a generalized likelihood ratio (GLR) test to rigorously assess the linearity of the interaction effect between multiple environments and a gene, while apply a parametric likelihood test to detect linear G×E interaction effect. The finite sample performance of the proposed method is examined through simulation studies and is further illustrated through a real data analysis.

Strategies for Enriching Variant Coverage in Candidate Disease Loci on a Multiethnic Genotyping Array.

Investigating genetic architecture of complex traits in ancestrally diverse populations is imperative to understand the etiology of disease. However, the current paucity of genetic research in people of African and Latin American ancestry, Hispanic and indigenous peoples in the United States is likely to exacerbate existing health disparities for many common diseases. The Population Architecture using Genomics and Epidemiology, Phase II (PAGE II), Study was initiated in 2013 by the National Human Genome Research Institute to expand our understanding of complex trait loci in ethnically diverse and well characterized study populations. To meet this goal, the Multi-Ethnic Genotyping Array (MEGA) was designed to substantially improve fine-mapping and functional discovery by increasing variant coverage across multiple ethnicities at known loci for metabolic, cardiovascular, renal, inflammatory, anthropometric, and a variety of lifestyle traits. Studying the frequency distribution of clinically relevant mutations, putative risk alleles, and known functional variants across multiple populations will provide important insight into the genetic architecture of complex diseases and facilitate the discovery of novel, sometimes population-specific, disease associations. DNA samples from 51,650 self-identified African ancestry (17,328), Hispanic/Latino (22,379), Asian/Pacific Islander (8,640), and American Indian (653) and an additional 2,650 participants of either South Asian or European ancestry, and other reference panels have been genotyped on MEGA by PAGE II. MEGA was designed as a new resource for studying ancestrally diverse populations. Here, we describe the methodology for selecting trait-specific content for use in multi-ethnic populations and how enriching MEGA for this content may contribute to deeper biological understanding of the genetic etiology of complex disease.

Using Systems Genetics to Understanding the Etiology of Complex Disease.

Here, we discuss Systems Genetics applications for systematic evaluation of candidate causal genes together with follow-up bioinformatics pathway analysis. The aim of this chapter is to illustrate analytic procedures and we provide examples in the context of Type 1 diabetes (T1day), the risk of which is conferred by over 60 loci. We also describe the Type 1 Diabetes Systems Genetics website and provide a guide for its use and application to other diseases.

MicroRNAs in Type 1 Diabetes: Complex Interregulation of the Immune System, β Cell Function and Viral Infections.

Since the discovery of the first mammalian microRNA (miRNA) more than two decades ago, a plethora of miRNAs has been identified in humans, now amounting to more than 2500. Essential for post-transcriptional regulation of gene networks integral for developmental pathways and immune response, it is not surprising that dysregulation of miRNAs is often associated with the aetiology of complex diseases including cancer, diabetes and autoimmune disorders. Despite massive expansion of small RNA studies and extensive investigation in diverse disease contexts, the role of miRNAs in type 1 diabetes has only recently been explored. Key studies using human islets have recently implicated virus-induced miRNA dysregulation as a pivotal mechanism of β cell destruction, while the interplay between miRNAs, the immune system and β cell survival has been illustrated in studies using animal and cellular models of disease. The role of specific miRNAs as major players in immune system homeostasis highlights their exciting potential as therapeutics and prognostic biomarkers of type 1 diabetes.

Gene-environment interactions between ERCC2, ERCC3, XRCC1 and cadmium exposure in nasal polyposis disease.

Gene-environment interactions have long been known to play an important role in complex disease aetiology, such as nasal polyposis (NP). The present study supports the concept that DNA repair gene polymorphisms play critical roles in modifying individual susceptibility to environmental diseases. In fact, we investigated the role of polymorphisms in DNA repair genes and cadmium as risk factors for Tunisian patients with NP. To the best of our knowledge, this is the first report on the impact of combined effects of cadmium and ERCC3 7122 A>G (rs4150407), ERCC2 Lys751Gln (rs13181) and XRCC1 Arg399Gln (rs25487) genes in the susceptibility to NP disease. Significant associations between the risk of developing NP disease and ERCC2 [odds ratio (OR) = 2.0, 95% confidence interval (CI) = 1.1-3.7, p = 0.023] and ERCC3 (OR = 2.2, 95% CI = 1.2-4.1, p = 0.013) genotypes polymorphisms were observed. Blood concentrations of Cd in NP patients (2.2μg/L) were significantly higher than those of controls (0.5μg/L). A significant interaction between ERCC3 (7122 A>G) polymorphism and blood-Cd levels (for the median of blood-Cd levels: OR = 3.8, 95% CI = 1.3-10.8, p = 0.014 and for the 75th percentiles of blood-Cd levels: OR = 2.7, 95% CI = 1.1-7.2, p = 0.041) was found in association with the risk of NP disease. In addition, when we stratified ERCC2, ERCC3 and XRCC1 polymorphism genotypes by the median and 75th percentiles of blood-Cd levels, we found also significant interactions between ERCC2 (Lys751Gln) and ERCC3 (7122 A>G) genotypes polymorphism and this metal in association with NP disease. However, no interaction was found between XRCC1 (Arg399Gln) polymorphism genotypes and Cd in association with NP disease.

Identification of low frequency and rare variants for hypertension using sparse-data methods.

Availability of genomic sequence data provides opportunities to study the role of low-frequency and rare variants in the etiology of complex disease. In this study, we conduct association analyses of hypertension status in the cohort of 1943 unrelated Mexican Americans provided by Genetic Analysis Workshop 19, focusing on exonic variants in MAP4 on chromosome 3. Our primary interest is to compare the performance of standard and sparse-data approaches for single-variant tests and variant-collapsing tests for sets of rare and low-frequency variants. We analyze both the real and the simulated phenotypes.

Constrained multivariate association with longitudinal phenotypes.

BackgroundThe incorporation of longitudinal data into genetic epidemiological studies has the potential to provide valuable information regarding the effect of time on complex disease etiology. Yet, the majority of research focuses on variables collected from a single time point. This aim of this study was to test for main effects on a quantitative trait across time points using a constrained maximum-likelihood measured genotype approach. This method simultaneously accounts for all repeat measurements of a phenotype in families. We applied this method to systolic blood pressure (SBP) measurements from three time points using the Genetic Analysis Workshop 19 (GAW19) whole-genome sequence family simulated data set and 200 simulated replicates. Data consisted of 849 individuals from 20 extended Mexican American pedigrees. Comparisons were made among 3 statistical approaches: (a) constrained, where the effect of a variant or gene region on the mean trait value was constrained to be equal across all measurements; (b) unconstrained, where the variant or gene region effect was estimated separately for each time point; and (c) the average SBP measurement from three time points. These approaches were run for nine genetic variants with known effect sizes (>0.001) for SBP variability and a known gene-centric kernel (MAP4)-based test under the GAW19 simulation model across 200 replicates.ResultsWhen compared to results using two time points, the constrained method utilizing all 3 time points increased power to detect association. Averaging SBP was equally effective when the variant has a large effect on the phenotype, but less powerful for variants with lower effect sizes. However, averaging SBP was far more effective than either the constrained or unconstrained approaches when using a gene-centric kernel-based test.ConclusionWe determined that this constrained multivariate approach improves genetic signal over the bivariate method. However, this method is still only effective in those variants that explain a moderate to large proportion of the phenotypic variance but is not as effective for gene-centric tests.

Genetic complexity at expression quantitative trait loci.

BackgroundIdentifying variants that regulate gene expression and delineating their genetic architecture is a critical next step in our endeavors to better understand the genetic etiology of complex diseases. The appropriate genomic tools are in place, and preliminary analytic strategies have been developed.MethodsHere we used Genetic Analysis Workshop (GAW) 19 data to investigate the genetic complexity of expression quantitative trait loci (eQTL), chromosomal regions likely to harbor regulatory elements responsible for gene expression. For this investigation, we analyzed the lymphocyte expression profiles of 653 individuals in 20 pedigrees who were also genotyped by single nucleotide polymorphism (SNP) arrays, followed by sequencing and imputation. We used these data to examine the degree of allelic heterogeneity, a contributor to genetic complexity at eQTL, by sequentially conditioning on the most significantly associated SNPs.ResultSOLAR (Sequential Oligogenic Linkage Analysis Routines)-MGA (measured genotype approach) and FaST-LMM (Factored Spectrally Transformed Linear Mixed Model) software allowed us to analyze pedigree data. The power and Type 1 error rates for single SNP association testing and multiple SNP sequential association testing were consistent for these programs. Sequential conditioning of the real expression data revealed substantial levels of allelic heterogeneity at the 2 eQTL examined, illustrating this feature of genetic complexity.ConclusionseQTL exhibit substantial genetic complexity among and within pedigrees.

A novel copy number variants kernel association test with application to autism spectrum disorders studies.

Copy number variants (CNVs) have been implicated in a variety of neurodevelopmental disorders, including autism spectrum disorders, intellectual disability and schizophrenia. Recent advances in high-throughput genomic technologies have enabled rapid discovery of many genetic variants including CNVs. As a result, there is increasing interest in studying the role of CNVs in the etiology of many complex diseases. Despite the availability of an unprecedented wealth of CNV data, methods for testing association between CNVs and disease-related traits are still under-developed due to the low prevalence and complicated multi-scale features of CNVs. We propose a novel CNV kernel association test (CKAT) in this paper. To address the low prevalence, CNVs are first grouped into CNV regions (CNVR). Then, taking into account the multi-scale features of CNVs, we first design a single-CNV kernel which summarizes the similarity between two CNVs, and next aggregate the single-CNV kernel to a CNVR kernel which summarizes the similarity between two CNVRs. Finally, association between CNVR and disease-related traits is assessed by comparing the kernel-based similarity with the similarity in the trait using a score test for variance components in a random effect model. We illustrate the proposed CKAT using simulations and show that CKAT is more powerful than existing methods, while always being able to control the type I error. We also apply CKAT to a real dataset examining the association between CNV and autism spectrum disorders, which demonstrates the potential usefulness of the proposed method. A R package to implement the proposed CKAT method is available at http://works.bepress.com/debashis_ghosh/ CONTACTS: xzhan@fhcrc.org or debashis.ghosh@ucdenver.eduSupplementary information: Supplementary data are available at Bioinformatics online.

A Clustered Multiclass Likelihood-Ratio Ensemble Method for Family-Based Association Analysis Accounting for Phenotypic Heterogeneity.

Although compelling evidence suggests that the genetic etiology of complex diseases could be heterogeneous in subphenotype groups, little attention has been paid to phenotypic heterogeneity in genetic association analysis of complex diseases. Simply ignoring phenotypic heterogeneity in association analysis could result in attenuated estimates of genetic effects and low power of association tests if subphenotypes with similar clinical manifestations have heterogeneous underlying genetic etiologies. To facilitate the family-based association analysis allowing for phenotypic heterogeneity, we propose a clustered multiclass likelihood-ratio ensemble (CMLRE) method. The proposed method provides an alternative way to model the complex relationship between disease outcomes and genetic variants. It allows for heterogeneous genetic causes of disease subphenotypes and can be applied to various pedigree structures. Through simulations, we found CMLRE outperformed the commonly adopted strategies in a variety of underlying disease scenarios. We further applied CMLRE to a family-based dataset from the International Consortium to Identify Genes and Interactions Controlling Oral Clefts (ICOC) to investigate the genetic variants and interactions predisposing to subphenotypes of oral clefts. The analysis suggested that two subphenotypes, nonsyndromic cleft lip without palate (CL) and cleft lip with palate (CLP), shared similar genetic etiologies, while cleft palate only (CP) had its own genetic mechanism. The analysis further revealed that rs10863790 (IRF6), rs7017252 (8q24), and rs7078160 (VAX1) were jointly associated with CL/CLP, while rs7969932 (TBK1), rs227731 (17q22), and rs2141765 (TBK1) jointly contributed to CP.

A Nonlinear Model for Gene-Based Gene-Environment Interaction

A vast amount of literature has confirmed the role of gene-environment (G×E) interaction in the etiology of complex human diseases. Traditional methods are predominantly focused on the analysis of interaction between a single nucleotide polymorphism (SNP) and an environmental variable. Given that genes are the functional units, it is crucial to understand how gene effects (rather than single SNP effects) are influenced by an environmental variable to affect disease risk. Motivated by the increasing awareness of the power of gene-based association analysis over single variant based approach, in this work, we proposed a sparse principle component regression (sPCR) model to understand the gene-based G×E interaction effect on complex disease. We first extracted the sparse principal components for SNPs in a gene, then the effect of each principal component was modeled by a varying-coefficient (VC) model. The model can jointly model variants in a gene in which their effects are nonlinearly influenced by an environmental variable. In addition, the varying-coefficient sPCR (VC-sPCR) model has nice interpretation property since the sparsity on the principal component loadings can tell the relative importance of the corresponding SNPs in each component. We applied our method to a human birth weight dataset in Thai population. We analyzed 12,005 genes across 22 chromosomes and found one significant interaction effect using the Bonferroni correction method and one suggestive interaction. The model performance was further evaluated through simulation studies. Our model provides a system approach to evaluate gene-based G×E interaction.