Family-Based Rare Variant Association Analysis: A Fast and Efficient Method of Multivariate Phenotype Association Analysis

Sex differences in the influence of type 2 diabetes (T2D)-related genes, parental history of T2D, and obesity on T2D development: a case-control study

BACKGROUND

This study investigated the effect of sex and age at type 2 diabetes (T2D) diagnosis on the influence of T2D-related genes, parental history of T2D, and obesity on T2D development.

METHODS

In this case-control study, 1012 T2D cases and 1008 healthy subjects were selected from the Diabetes in Mexico Study database. Participants were stratified by sex and age at T2D diagnosis (early, ≤ 45 years; late, ≥ 46 years). Sixty-nine T2D-associated single nucleotide polymorphisms were explored and the percentage contribution (R²) of T2D-related genes, parental history of T2D, and obesity (body mass index [BMI] and waist-hip ratio [WHR]) on T2D development was calculated using univariate and multivariate logistic regression models.

RESULTS

T2D-related genes influenced T2D development most in males who were diagnosed early (R² = 23.5%; females, R² = 13.5%; males and females diagnosed late, R² = 11.9% and R² = 7.3%, respectively). With an early diagnosis, insulin production-related genes were more influential in males (76.0% of R²) while peripheral insulin resistance-associated genes were more influential in females (52.3% of R²). With a late diagnosis, insulin production-related genes from chromosome region 11p15.5 notably influenced males while peripheral insulin resistance and genes associated with inflammation and other processes notably influenced females. Influence of parental history was higher among those diagnosed early (males, 19.9%; females, 17.5%) versus late (males, 6.4%; females, 5,3%). Unilateral maternal T2D history was more influential than paternal T2D history. BMI influenced T2D development for all, while WHR exclusively influenced males.

CONCLUSIONS

The influence of T2D-related genes, maternal T2D history, and fat distribution on T2D development was greater in males than females.

MF-TOWmuT: Testing an optimally weighted combination of common and rare variants with multiple traits using family data

With rapid advancements of sequencing technologies and accumulations of electronic health records, a large number of genetic variants and multiple correlated human complex traits have become available in many genetic association studies. Thus, it becomes necessary and important to develop new methods that can jointly analyze the association between multiple genetic variants and multiple traits. Compared with methods that only use a single marker or trait, the joint analysis of multiple genetic variants and multiple traits is more powerful since such an analysis can fully incorporate the correlation structure of genetic variants and/or traits and their mutual dependence patterns. However, most of existing methods that simultaneously analyze multiple genetic variants and multiple traits are only applicable to unrelated samples. We develop a new method called MF-TOWmuT to detect association of multiple phenotypes and multiple genetic variants in a genomic region with family samples. MF-TOWmuT is based on an optimally weighted combination of variants. Our method can be applied to both rare and common variants and both qualitative and quantitative traits. Our simulation results show that (1) the type I error of MF-TOWmuT is preserved; (2) MF-TOWmuT outperforms two existing methods such as Multiple Family-based Quasi-Likelihood Score Test and Multivariate Family-based Rare Variant Association Test in terms of power. We also illustrate the usefulness of MF-TOWmuT by analyzing genotypic and phenotipic data from the Genetics of Kidneys in Diabetes study. R program is available at https://github.com/gaochengPRC/MF-TOWmuT.

Associations of Mitochondrial and Nuclear Mitochondrial Variants and Genes with Seven Metabolic Traits

Mitochondria (MT), the major site of cellular energy production, are under dual genetic control by 37 mitochondrial DNA (mtDNA) genes and numerous nuclear genes (MT-nDNA). In the CHARGEmtDNA+ Consortium, we studied genetic associations of mtDNA and MT-nDNA associations with body mass index (BMI), waist-hip-ratio (WHR), glucose, insulin, HOMA-B, HOMA-IR, and HbA1c. This 45-cohort collaboration comprised 70,775 (insulin) to 170,202 (BMI) pan-ancestry individuals. Validation and imputation of mtDNA variants was followed by single-variant and gene-based association testing. We report two significant common variants, one in MT-ATP6 associated (p ≤ 5E-04) with WHR and one in the D-loop with glucose. Five rare variants in MT-ATP6, MT-ND5, and MT-ND6 associated with BMI, WHR, or insulin. Gene-based meta-analysis identified MT-ND3 associated with BMI (p ≤ 1E-03). We considered 2,282 MT-nDNA candidate gene associations compiled from online summary results for our traits (20 unique studies with 31 dataset consortia's genome-wide associations [GWASs]). Of these, 109 genes associated (p ≤ 1E-06) with at least 1 of our 7 traits. We assessed regulatory features of variants in the 109 genes, cis- and trans-gene expression regulation, and performed enrichment and protein-protein interactions analyses. Of the identified mtDNA and MT-nDNA genes, 79 associated with adipose measures, 49 with glucose/insulin, 13 with risk for type 2 diabetes, and 18 with cardiovascular disease, indicating for pleiotropic effects with health implications. Additionally, 21 genes related to cholesterol, suggesting additional important roles for the genes identified. Our results suggest that mtDNA and MT-nDNA genes and variants reported make important contributions to glucose and insulin metabolism, adipocyte regulation, diabetes, and cardiovascular disease.

Strategic Methods for Recruiting Grandparents: The Tohoku Medical Megabank Birth and Three-Generation Cohort Study

Involvement of family members, especially grandparents, in genome epidemiological research is important to investigate both genetic and environmental factors of common diseases. The aim of the present study was to establish strategies to obtain enough number of family recruitment, especially focusing on grandparents, for the Tohoku Medical Megabank Birth and Three-Generation Cohort Study. Our main strategies are summarized below. 1) We standardized informed consent process with reference materials to help people understand the consent form, 2) we created an invitation letter to contact family members, and 3) we recruited family members in several settings. To obtain informed consent, we were careful of explaining clearly the complex reasons as well as drawing people's attention. By the end of March 2017, the number of invitation letters distributed to family members through the pregnant women was 23,806, including 18,702 grandparents. Among the grandparents who received invitation letters, 2,935 (15.7%) responded to us. Furthermore, some grandparents were asked to provide informed consent with other family members by staff at maternal clinics or Community Support Centers, and others directly booked Community Support Centers without responding to the invitation letter. Grandparents joined the study anytime during mother's maternal check-ups or delivery. Overall, 8,054 grandparents participated in our birth cohort study. The setting in which most grandparents were recruited was our own facilities. Importantly, both paternal and maternal grandparents more frequently participated in the study if the father also participated. In conclusion, we are able to recruit not only pregnant women but also fathers and grandparents.

CpG-set association assessment of lipid concentration changes and DNA methylation

Epigenome association studies that test a large number of methylation sites suffer from stringent multiple-testing corrections. This study's goals were to investigate region-based associations between DNA methylation sites and lipid-level changes in response to the treatment with fenofibrate in the GAW20 data and to investigate whether improvements in power could be obtained by taking into account correlations between DNA methylation at neighboring cytosine-phosphate-guanine (CpG) sites. To this end, we applied both a recently developed block-based data-dimension-reduction approach and a region-based variance-component (VC) linear mixed model to GAW20 data. We compared analyses of unrelated individuals with familial data. The region-based VC approach using unrelated (independent) individuals identified the gene LGALS9C as significantly associated with changes in triglycerides. However, univariate tests of individual CpG sites yielded no valid statistically significant results.

Pathway-based approach using hierarchical components of rare variants to analyze multiple phenotypes

BACKGROUND

As one possible solution to the "missing heritability" problem, many methods have been proposed that apply pathway-based analyses, using rare variants that are detected by next generation sequencing technology. However, while a number of methods for pathway-based rare-variant analysis of multiple phenotypes have been proposed, no method considers a unified model that incorporate multiple pathways.

RESULTS

Simulation studies successfully demonstrated advantages of multivariate analysis, compared to univariate analysis, and comparison studies showed the proposed approach to outperform existing methods. Moreover, real data analysis of six type 2 diabetes-related traits, using large-scale whole exome sequencing data, identified significant pathways that were not found by univariate analysis. Furthermore, strong relationships between the identified pathways, and their associated metabolic disorder risk factors, were found via literature search, and one of the identified pathway, was successfully replicated by an analysis with an independent dataset.

CONCLUSIONS

Herein, we present a powerful, pathway-based approach to investigate associations between multiple pathways and multiple phenotypes. By reflecting the natural hierarchy of biological behavior, and considering correlation between pathways and phenotypes, the proposed method is capable of analyzing multiple phenotypes and multiple pathways simultaneously.

WISARD: workbench for integrated superfast association studies for related datasets

BACKGROUND

A Mendelian transmission produces phenotypic and genetic relatedness between family members, giving family-based analytical methods an important role in genetic epidemiological studies-from heritability estimations to genetic association analyses. With the advance in genotyping technologies, whole-genome sequence data can be utilized for genetic epidemiological studies, and family-based samples may become more useful for detecting de novo mutations. However, genetic analyses employing family-based samples usually suffer from the complexity of the computational/statistical algorithms, and certain types of family designs, such as incorporating data from extended families, have rarely been used.

RESULTS

We present a Workbench for Integrated Superfast Association studies for Related Data (WISARD) programmed in C/C++. WISARD enables the fast and a comprehensive analysis of SNP-chip and next-generation sequencing data on extended families, with applications from designing genetic studies to summarizing analysis results. In addition, WISARD can automatically be run in a fully multithreaded manner, and the integration of R software for visualization makes it more accessible to non-experts.

CONCLUSIONS

Comparison with existing toolsets showed that WISARD is computationally suitable for integrated analysis of related subjects, and demonstrated that WISARD outperforms existing toolsets. WISARD has also been successfully utilized to analyze the large-scale massive sequencing dataset of chronic obstructive pulmonary disease data (COPD), and we identified multiple genes associated with COPD, which demonstrates its practical value.

Identifying pleiotropic genes in genome-wide association studies from related subjects using the linear mixed model and Fisher combination function

Background

A multivariate genome-wide association test is proposed for analyzing data on multivariate quantitative phenotypes collected from related subjects. The proposed method is a two-step approach. The first step models the association between the genotype and marginal phenotype using a linear mixed model. The second step uses the correlation between residuals of the linear mixed model to estimate the null distribution of the Fisher combination test statistic.

Results

The simulation results show that the proposed method controls the type I error rate and is more powerful than the marginal tests across different population structures (admixed or non-admixed) and relatedness (related or independent). The statistical analysis on the database of the Study of Addiction: Genetics and Environment (SAGE) demonstrates that applying the multivariate association test may facilitate identification of the pleiotropic genes contributing to the risk for alcohol dependence commonly expressed by four correlated phenotypes.

Conclusions

This study proposes a multivariate method for identifying pleiotropic genes while adjusting for cryptic relatedness and population structure between subjects. The two-step approach is not only powerful but also computationally efficient even when the number of subjects and the number of phenotypes are both very large.