Tin Louie

Exome sequencing of extreme phenotypes identifies DCTN4 as a modifier of chronic Pseudomonas aeruginosa infection in cystic fibrosis

Exome sequencing has become a powerful and effective strategy for the discovery of genes underlying Mendelian disorders. However, use of exome sequencing to identify variants associated with complex traits has been more challenging, partly because the sample sizes needed for adequate power may be very large. One strategy to increase efficiency is to sequence individuals who are at both ends of a phenotype distribution (those with extreme phenotypes). Because the frequency of alleles that contribute to the trait are enriched in one or both phenotype extremes, a modest sample size can potentially be used to identify novel candidate genes and/or alleles. As part of the National Heart, Lung, and Blood Institute (NHLBI) Exome Sequencing Project (ESP), we used an extreme phenotype study design to discover that variants in DCTN4, encoding a dynactin protein, are associated with time to first P. aeruginosa airway infection, chronic P. aeruginosa infection and mucoid P. aeruginosa in individuals with cystic fibrosis.

Identification of Rare Variants in ATP8B4 as a Risk Factor for Systemic Sclerosis by Whole-Exome Sequencing.

To determine the contribution of rare variants as genetic modifiers of the expressivity, penetrance, and severity of systemic sclerosis (SSc).

Whole‐exome sequencing identifies rare variants in ATP8B4 as a risk factor for systemic sclerosis

To determine the contribution of rare variants as genetic modifiers of the expressivity, penetrance, and severity of systemic sclerosis (SSc).

Exome Sequencing of Phenotypic Extremes Identifies CAV2 and TMC6 as Interacting Modifiers of Chronic Pseudomonas aeruginosa Infection in Cystic Fibrosis

Discovery of rare or low frequency variants in exome or genome data that are associated with complex traits often will require use of very large sample sizes to achieve adequate statistical power. For a fixed sample size, sequencing of individuals sampled from the tails of a phenotype distribution (i.e., extreme phenotypes design) maximizes power and this approach was recently validated empirically with the discovery of variants in DCTN4 that influence the natural history of P. aeruginosa airway infection in persons with cystic fibrosis (CF; MIM219700). The increasing availability of large exome/genome sequence datasets that serve as proxies for population-based controls affords the opportunity to test an alternative, potentially more powerful and generalizable strategy, in which the frequency of rare variants in a single extreme phenotypic group is compared to a control group (i.e., extreme phenotype vs. control population design). As proof-of-principle, we applied this approach to search for variants associated with risk for age-of-onset of chronic P. aeruginosa airway infection among individuals with CF and identified variants in CAV2 and TMC6 that were significantly associated with group status. These results were validated using a large, prospective, longitudinal CF cohort and confirmed a significant association of a variant in CAV2 with increased age-of-onset of P. aeruginosa airway infection (hazard ratio = 0.48, 95% CI=[0.32, 0.88]) and variants in TMC6 with diminished age-of-onset of P. aeruginosa airway infection (HR = 5.4, 95% CI=[2.2, 13.5]) A strong interaction between CAV2 and TMC6 variants was observed (HR=12.1, 95% CI=[3.8, 39]) for children with the deleterious TMC6 variant and without the CAV2 protective variant. Neither gene showed a significant association using an extreme phenotypes design, and conditions for which the power of an extreme phenotype vs. control population design was greater than that for the extreme phenotypes design were explored.

Genome-wide association study of heavy smoking and daily/nondaily smoking in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL)

Introduction Genetic variants associated with nicotine dependence have previously been identified, primarily in European-ancestry populations. No genome-wide association studies (GWAS) have been reported for smoking behaviors in Hispanics/Latinos in the United States and Latin America, who are of mixed ancestry with European, African, and American Indigenous components. Methods We examined genetic associations with smoking behaviors in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL) (N = 12 741 with smoking data, 5119 ever-smokers), using ~2.3 million genotyped variants imputed to the 1000 Genomes Project phase 3. Mixed logistic regression models accounted for population structure, sampling, relatedness, sex, and age. Results The known region of CHRNA5, which encodes the α5 cholinergic nicotinic receptor subunit, was associated with heavy smoking at genome-wide significance (p ≤ 5 × 10-8) in a comparison of 1929 ever-smokers reporting cigarettes per day (CPD) > 10 versus 3156 reporting CPD ≤ 10. The functional variant rs16969968 in CHRNA5 had a p value of 2.20 × 10-7 and odds ratio (OR) of 1.32 for the minor allele (A); its minor allele frequency was 0.22 overall and similar across Hispanic/Latino background groups (Central American = 0.17; South American = 0.19; Mexican = 0.18; Puerto Rican = 0.22; Cuban = 0.29; Dominican = 0.19). CHRNA4 on chromosome 20 attained p < 10-4, supporting prior findings in non-Hispanics. For nondaily smoking, which is prevalent in Hispanic/Latino smokers, compared to daily smoking, loci on chromosomes 2 and 4 achieved genome-wide significance; replication attempts were limited by small Hispanic/Latino sample sizes. Conclusions Associations of nicotinic receptor gene variants with smoking, first reported in non-Hispanic European-ancestry populations, generalized to Hispanics/Latinos despite different patterns of smoking behavior. Implications We conducted the first large-scale genome-wide association study (GWAS) of smoking behavior in a US Hispanic/Latino cohort, and the first GWAS of daily/nondaily smoking in any population. Results show that the region of the nicotinic receptor subunit gene CHRNA5, which in non-Hispanic European-ancestry smokers has been associated with heavy smoking as well as cessation and treatment efficacy, is also significantly associated with heavy smoking in this Hispanic/Latino cohort. The results are an important addition to understanding the impact of genetic variants in understudied Hispanic/Latino smokers.

Variation in Cilia Protein Genes and Progression of Lung Disease in Cystic Fibrosis

Rationale: Cystic fibrosis, like primary ciliary dyskinesia, is an autosomal recessive disorder characterized by abnormal mucociliary clearance and obstructive lung disease. We hypothesized that genes underlying the development or function of cilia may modify lung disease severity in persons with cystic fibrosis. Objectives: To test this hypothesis, we compared variants in 93 candidate genes in both upper and lower tertiles of lung function in a large cohort of children and adults with cystic fibrosis with those of a population control dataset. Methods: Variants within candidate genes were tested for association using the SKAT‐O test, comparing cystic fibrosis cases defined by poor (n = 127) or preserved (n = 127) lung function with population controls (n = 3,269 or 3,148, respectively). Associated variants were then tested for association with related phenotypes in independent datasets. Results: Variants in DNAH14 and DNAAF3 were associated with poor lung function in cystic fibrosis, whereas variants in DNAH14 and DNAH6 were associated with preserved lung function in cystic fibrosis. Associations between DNAH14 and lung function were replicated in disease‐related phenotypes characterized by obstructive lung disease in adults. Conclusions: Genetic variants within DNAH6, DNAH14, and DNAAF3 are associated with variation in lung function among persons with cystic fibrosis.

The Epithelial Sodium Channel Is a Modifier of the Long-Term Nonprogressive Phenotype Associated with F508del CFTR Mutations

&NA; Cystic fibrosis (CF) remains the most lethal genetic disease in the Caucasian population. However, there is great variability in clinical phenotypes and survival times, even among patients harboring the same genotype. We identified five patients with CF and a homozygous F508del mutation in the CFTR gene who were in their fifth or sixth decade of life and had shown minimal changes in lung function over a longitudinal period of more than 20 years. Because of the rarity of this long‐term nonprogressive phenotype, we hypothesized these individuals may carry rare genetic variants in modifier genes that ameliorate disease severity. Individuals at the extremes of survival time and lung‐function trajectory underwent whole‐exome sequencing, and the sequencing data were filtered to include rare missense, stopgain, indel, and splicing variants present with a mean allele frequency of <0.2% in general population databases. Epithelial sodium channel (ENaC) mutants were generated via site‐directed mutagenesis and expressed for Xenopus oocyte assays. Four of the five individuals carried extremely rare or never reported variants in the SCNN1D and SCNN1B genes of the ENaC. Separately, an independently enriched rare variant in SCNN1D was identified in the Exome Variant Server database associated with a milder pulmonary disease phenotype. Functional analysis using Xenopus oocytes revealed that two of the three variants in &dgr;‐ENaC encoded by SCNN1D exhibited hypomorphic channel activity. Our data suggest a potential role for &dgr;‐ENaC in controlling sodium reabsorption in the airways, and advance the plausibility of ENaC as a therapeutic target in CF.

Genome-wide association study of iron traits and relation to diabetes in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL): potential genomic intersection of iron and glucose regulation?

&NA; Genetic variants contribute to normal variation of iron‐related traits and may also cause clinical syndromes of iron deficiency or excess. Iron overload and deficiency can adversely affect human health. For example, elevated iron storage is associated with increased diabetes risk, although mechanisms are still being investigated. We conducted the first genome‐wide association study of serum iron, total iron binding capacity (TIBC), transferrin saturation, and ferritin in a Hispanic/Latino cohort, the Hispanic Community Health Study/Study of Latinos (>12 000 participants) and also assessed the generalization of previously known loci to this population. We then evaluated whether iron‐associated variants were associated with diabetes and glycemic traits. We found evidence for a novel association between TIBC and a variant near the gene for protein phosphatase 1, regulatory subunit 3B (PPP1R3B; rs4841132, &bgr; = −0.116, P = 7.44 × 10−8). The effect strengthened when iron deficient individuals were excluded (&bgr; = −0.121, P = 4.78 × 10−9). Ten of sixteen variants previously associated with iron traits generalized to HCHS/SOL, including variants at the transferrin (TF), hemochromatosis (HFE), fatty acid desaturase 2 (FADS2)/myelin regulatory factor (MYRF), transmembrane protease, serine 6 (TMPRSS6), transferrin receptor (TFR2), N‐acetyltransferase 2 (arylamine N‐acetyltransferase) (NAT2), ABO blood group (ABO), and GRB2 associated binding protein 3 (GAB3) loci. In examining iron variant associations with glucose homeostasis, an iron‐raising variant of TMPRSS6 was associated with lower HbA1c levels (P = 8.66 × 10−10). This association was attenuated upon adjustment for iron measures. In contrast, the iron‐raising allele of PPP1R3B was associated with higher levels of fasting glucose (P = 7.70 × 10−7) and fasting insulin (P = 4.79 × 10−6), but these associations were not attenuated upon adjustment for TIBC—so iron is not likely a mediator. These results provide new genetic information on iron traits and their connection with glucose homeostasis.

Genetic variation near IRS1 is associated with adiposity and a favorable metabolic profile in U.S. Hispanics/Latinos.

Associations of IRS1 genetic variation with adiposity and metabolic profile in U.S. Hispanic/Latino individuals of diverse backgrounds were examined.

Exome sequencing in HFE C282Y homozygous men with extreme phenotypes identifies a GNPAT variant associated with severe iron overload

To identify polymorphisms associated with variability of iron overload severity in HFE‐associated hemochromatosis, we performed exome sequencing of DNA from 35 male HFE C282Y homozygotes with either markedly increased iron stores (n = 22; cases) or with normal or mildly increased iron stores (n = 13; controls). The 35 participants, residents of the United States, Canada, and Australia, reported no or light alcohol consumption. Sequencing data included 82,068 single‐nucleotide variants, and 10,337 genes were tested for a difference between cases and controls. A variant in the GNPAT gene showed the most significant association with severe iron overload (P = 3 × 10−6; P = 0.033 by the likelihood ratio test after correction for multiple comparisons). Sixteen of twenty‐two participants with severe iron overload had glyceronephosphate O‐acyltransferase (GNPAT) polymorphism p.D519G (rs11558492; 15 heterozygotes, one homozygote). No control participant had this polymorphism. To examine functional consequences of GNPAT deficiency, we performed small interfering RNA–based knockdown of GNPAT in the human liver‐derived cell line, HepG2/C3A. This knockdown resulted in a >17‐fold decrease in expression of the messenger RNA encoding the iron‐regulatory hormone, hepcidin. Conclusion: GNPAT p.D519G is associated with a high‐iron phenotype in HFE C282Y homozygotes and may participate in hepcidin regulation. (Hepatology 2015;62:429–439