Rany M. Salem

Heritability of Nonalcoholic Fatty Liver Disease

BACKGROUND & AIMS Nonalcoholic fatty liver disease (NAFLD) is the most common chronic liver disease in the United States. The etiology is believed to be multifactorial with a substantial genetic component; however, the heritability of NAFLD is undetermined. Therefore, a familial aggregation study was performed to test the hypothesis that NAFLD is highly heritable. METHODS Overweight children with biopsy-proven NAFLD and overweight children without NAFLD served as probands. Family members were studied, including the use of magnetic resonance imaging to quantify liver fat fraction. Fatty liver was defined as a liver fat fraction of 5% or higher. Etiologies for fatty liver other than NAFLD were excluded. Narrow-sense heritability estimates for fatty liver (dichotomous) and fat fraction (continuous) were calculated using variance components analysis adjusted for covariate effects. RESULTS Fatty liver was present in 17% of siblings and 37% of parents of overweight children without NAFLD. Fatty liver was significantly more common in siblings (59%) and parents (78%) of children with NAFLD. Liver fat fraction was correlated with body mass index, although the correlation was significantly stronger for families of children with NAFLD than those without NAFLD. Adjusted for age, sex, race, and body mass index, the heritability of fatty liver was 1.000 and of liver fat fraction was 0.386. CONCLUSIONS Family members of children with NAFLD should be considered at high risk for NAFLD. These data suggest that familial factors are a major determinant of whether an individual has NAFLD. Studies examining the complex relations between genes and environment in the development and progression of NAFLD are warranted.

Distribution and Medical Impact of Loss-of-Function Variants in the Finnish Founder Population.

Exome sequencing studies in complex diseases are challenged by the allelic heterogeneity, large number and modest effect sizes of associated variants on disease risk and the presence of large numbers of neutral variants, even in phenotypically relevant genes. Isolated populations with recent bottlenecks offer advantages for studying rare variants in complex diseases as they have deleterious variants that are present at higher frequencies as well as a substantial reduction in rare neutral variation. To explore the potential of the Finnish founder population for studying low-frequency (0.5–5%) variants in complex diseases, we compared exome sequence data on 3,000 Finns to the same number of non-Finnish Europeans and discovered that, despite having fewer variable sites overall, the average Finn has more low-frequency loss-of-function variants and complete gene knockouts. We then used several well-characterized Finnish population cohorts to study the phenotypic effects of 83 enriched loss-of-function variants across 60 phenotypes in 36,262 Finns. Using a deep set of quantitative traits collected on these cohorts, we show 5 associations (p<5×10−8) including splice variants in LPA that lowered plasma lipoprotein(a) levels (P = 1.5×10−117). Through accessing the national medical records of these participants, we evaluate the LPA finding via Mendelian randomization and confirm that these splice variants confer protection from cardiovascular disease (OR = 0.84, P = 3×10−4), demonstrating for the first time the correlation between very low levels of LPA in humans with potential therapeutic implications for cardiovascular diseases. More generally, this study articulates substantial advantages for studying the role of rare variation in complex phenotypes in founder populations like the Finns and by combining a unique population genetic history with data from large population cohorts and centralized research access to National Health Registers.

Efficient Bayesian mixed model analysis increases association power in large cohorts

Linear mixed models are a powerful statistical tool for identifying genetic associations and avoiding confounding. However, existing methods are computationally intractable in large cohorts and may not optimize power. All existing methods require time cost O(MN2) (where N is the number of samples and M is the number of SNPs) and implicitly assume an infinitesimal genetic architecture in which effect sizes are normally distributed, which can limit power. Here we present a far more efficient mixed-model association method, BOLT-LMM, which requires only a small number of O(MN) time iterations and increases power by modeling more realistic, non-infinitesimal genetic architectures via a Bayesian mixture prior on marker effect sizes. We applied BOLT-LMM to 9 quantitative traits in 23,294 samples from the Womens Genome Health Study (WGHS) and observed significant increases in power, consistent with simulations. Theory and simulations show that the boost in power increases with cohort size, making BOLT-LMM appealing for genome-wide association studies in large cohorts.

New susceptibility loci associated with kidney disease in Type 1 diabetes

Diabetic kidney disease, or diabetic nephropathy (DN), is a major complication of diabetes and the leading cause of end-stage renal disease (ESRD) that requires dialysis treatment or kidney transplantation. In addition to the decrease in the quality of life, DN accounts for a large proportion of the excess mortality associated with type 1 diabetes (T1D). Whereas the degree of glycemia plays a pivotal role in DN, a subset of individuals with poorly controlled T1D do not develop DN. Furthermore, strong familial aggregation supports genetic susceptibility to DN. However, the genes and the molecular mechanisms behind the disease remain poorly understood, and current therapeutic strategies rarely result in reversal of DN. In the GEnetics of Nephropathy: an International Effort (GENIE) consortium, we have undertaken a meta-analysis of genome-wide association studies (GWAS) of T1D DN comprising ∼2.4 million single nucleotide polymorphisms (SNPs) imputed in 6,691 individuals. After additional genotyping of 41 top ranked SNPs representing 24 independent signals in 5,873 individuals, combined meta-analysis revealed association of two SNPs with ESRD: rs7583877 in the AFF3 gene (P = 1.2×10−8) and an intergenic SNP on chromosome 15q26 between the genes RGMA and MCTP2, rs12437854 (P = 2.0×10−9). Functional data suggest that AFF3 influences renal tubule fibrosis via the transforming growth factor-beta (TGF-β1) pathway. The strongest association with DN as a primary phenotype was seen for an intronic SNP in the ERBB4 gene (rs7588550, P = 2.1×10−7), a gene with type 2 diabetes DN differential expression and in the same intron as a variant with cis-eQTL expression of ERBB4. All these detected associations represent new signals in the pathogenesis of DN.

Catecholamine release-inhibitory peptide catestatin (chromogranin A352-372): Naturally occurring amino acid variant Gly364Ser causes profound changes in human autonomic activity and alters risk for hypertension

Background— Chromogranin A, coreleased with catecholamines by exocytosis, is cleaved to the catecholamine release–inhibitory fragment catestatin. We identified a natural nonsynonymous variant of catestatin, Gly364Ser, that alters human autonomic function and blood pressure. Methods and Results— Gly364Ser heterozygotes and controls underwent physiological and biochemical phenotyping, including catecholamine production, chromogranin A precursor, and its catestatin product. Case-control studies replicated effects of the gene on blood pressure in the population. Gly364Ser displayed diminished inhibition of catecholamine secretion from cultured neurons. Gly/Ser heterozygotes displayed increased baroreceptor slope during upward deflections (by ≈47%) and downward deflections (by ≈44%), increased cardiac parasympathetic index (by ≈2.4-fold), and decreased cardiac sympathetic index (by ≈26%). Renal norepinephrine excretion was diminished by ≈26% and epinephrine excretion by ≈34% in Gly/Ser heterozygotes. The coalescent dated emergence of the variant to ≈70 000 years ago. Gly364Ser was in linkage disequilibrium with 1 major Chromogranin A promoter haplotype, although promoter haplotypes did not predict autonomic phenotypes. The 364Ser variant was associated with lower diastolic blood pressure in 2 independent/confirmatory groups of patients with hypertension; genotype groups differed by ≈5 to 6 mm Hg, and the polymorphism accounted for ≈1.8% of population diastolic blood pressure variance, although a significant gene-by-sex interaction existed, with an enhanced effect in men. Conclusions— The catestatin Gly364Ser variant causes profound changes in human autonomic activity, both parasympathetic and sympathetic, and seems to reduce risk of developing hypertension, especially in men. A model for catestatin action in the baroreceptor center of the nucleus of the tractus solitarius accounts for these actions.

Longitudinal genome-wide association of cardiovascular disease risk factors in the Bogalusa heart study

Cardiovascular disease (CVD) is the leading cause of death worldwide. Recent genome-wide association (GWA) studies have pinpointed many loci associated with CVD risk factors in adults. It is unclear, however, if these loci predict trait levels at all ages, if they are associated with how a trait develops over time, or if they could be used to screen individuals who are pre-symptomatic to provide the opportunity for preventive measures before disease onset. We completed a genome-wide association study on participants in the longitudinal Bogalusa Heart Study (BHS) and have characterized the association between genetic factors and the development of CVD risk factors from childhood to adulthood. We report 7 genome-wide significant associations involving CVD risk factors, two of which have been previously reported. Top regions were tested for replication in the Young Finns Study (YF) and two associations strongly replicated: rs247616 in CETP with HDL levels (combined P = 9.7×10−24), and rs445925 at APOE with LDL levels (combined P = 8.7×10−19). We show that SNPs previously identified in adult cross-sectional studies tend to show age-independent effects in the BHS with effect sizes consistent with previous reports. Previously identified variants were associated with adult trait levels above and beyond those seen in childhood; however, variants with time-dependent effects were also promising predictors. This is the first GWA study to evaluate the role of common genetic variants in the development of CVD risk factors in children as they advance through adulthood and highlights the utility of using longitudinal studies to identify genetic predictors of adult traits in children.

Tyrosine Hydroxylase, the Rate-Limiting Enzyme in Catecholamine Biosynthesis Discovery of Common Human Genetic Variants Governing Transcription, Autonomic Activity, and Blood Pressure In Vivo

Background— Tyrosine hydroxylase (TH) is the rate-limiting enzyme in catecholamine biosynthesis. Does common genetic variation at human TH alter autonomic activity and predispose to cardiovascular disease? We undertook systematic polymorphism discovery at the TH locus and then tested variants for contributions to sympathetic function and blood pressure. Methods and Results— We resequenced 80 ethnically diverse individuals across the TH locus. One hundred seventy-two twin pairs were evaluated for sympathetic traits, including catecholamine production, reflex control of the circulation, and environmental (cold) stress responses. To evaluate hypertension, we genotyped subjects selected from the most extreme diastolic blood pressure percentiles in the population. Human TH promoter haplotype/reporter plasmids were transfected into chromaffin cells. Forty-nine single-nucleotide polymorphisms were discovered, but coding region polymorphism did not account for common phenotypic variation. A block of linkage disequilibrium spanned 4 common variants in the proximal promoter. Catecholamine secretory traits were significantly heritable (h2), as were stress-induced blood pressure changes. In the TH promoter, significant associations were found for urinary catecholamine excretion and for blood pressure response to stress. TH promoter haplotype 2 (TGGG) showed pleiotropy, increasing both norepinephrine excretion and blood pressure during stress. Coalescent simulations suggest that TH haplotype 2 likely arose ≈380 000 years ago. In hypertension, 2 independent case-control studies (1266 subjects with 53% women and 927 subjects with 24% women) replicated the effect of C-824T in the determination of blood pressure. Conclusions— We conclude that human catecholamine secretory traits are heritable, displaying joint genetic determination (pleiotropy) with autonomic activity and finally with blood pressure in the population. Catecholamine secretion is influenced by genetic variation in the adrenergic pathway encoding catecholamine synthesis, especially at the classically rate-limiting step, TH. The results suggest novel pathophysiological links between a key adrenergic locus, catecholamine metabolism, and blood pressure and suggest new strategies to approach the mechanism, diagnosis, and treatment of systemic hypertension.

A comprehensive literature review of haplotyping software and methods for use with unrelated individuals

Interest in the assignment and frequency analysis of haplotypes in samples of unrelated individuals has increased immeasurably as a result of the emphasis placed on haplotype analyses by, for example, the International HapMap Project and related initiatives. Although there are many available computer programs for haplotype analysis applicable to samples of unrelated individuals, many of these programs have limitations and/or very specific uses. In this paper, the key features of available haplotype analysis software for use with unrelated individuals, as well as pooled DNA samples from unrelated individuals, are summarised. Programs for haplotype analysis were identified through keyword searches on PUBMED and various internet search engines, a review of citations from retrieved papers and personal communications, up to June 2004. Priority was given to functioning computer programs, rather than theoretical models and methods. The available software was considered in light of a number of factors: the algorithm(s) used, algorithm accuracy, assumptions, the accommodation of genotyping error, implementation of hypothesis testing, handling of missing data, software characteristics and web-based implementations. Review papers comparing specific methods and programs are also summarised. Forty-six haplotyping programs were identified and reviewed. The programs were divided into two groups: those designed for individual genotype data (a total of 43 programs) and those designed for use with pooled DNA samples (a total of three programs). The accuracy of programs using various criteria are assessed and the programs are categorised and discussed in light of: algorithm and method, accuracy, assumptions, genotyping error, hypothesis testing, missing data, software characteristics and web implementation. Many available programs have limitations (eg some cannot accommodate missing data) and/or are designed with specific tasks in mind (eg estimating haplotype frequencies rather than assigning most likely haplotypes to individuals). It is concluded that the selection of an appropriate haplotyping program for analysis purposes should be guided by what is known about the accuracy of estimation, as well as by the limitations and assumptions built into a program.

Direct Vasoactive Effects of the Chromogranin A (CHGA) Peptide Catestatin in Humans In Vivo

Catestatin is a bioactive peptide of chromogranin A (CHGA) that is co-released with catecholamines from secretory vesicles. Catestatin may function as a vasodilator and is diminished in hypertension. To evaluate this potential vasodilator in vivo without systemic counterregulation, we infused catestatin to target concentrations of ∼ 50, ∼ 500, ∼ 5000 nM into dorsal hand veins of 18 normotensive men and women, after pharmacologic venoconstriction with phenylephrine. Pancreastatin, another CHGA peptide, was infused as a negative control. After preconstriction to ∼ 69%, increasing concentrations of catestatin resulted in dose-dependent vasodilation (P = 0.019), in female subjects (to ∼ 44%) predominantly. The EC50 (∼ 30 nM) for vasodilation induced by catestatin was the same order of magnitude to circulating endogenous catestatin (4.4 nM). No vasodilation occurred during the control infusion with pancreastatin. Plasma CHGA, catestatin, and CHGA-to-catestatin processing were then determined in 622 healthy subjects without hypertension. Female subjects had higher plasma catestatin levels than males (P = 0.001), yet lower CHGA precursor concentrations (P = 0.006), reflecting increased processing of CHGA-to-catestatin (P < 0.001). Our results demonstrate that catestatin dilates human blood vessels in vivo, especially in females. Catestatin may contribute to sex differences in endogenous vascular tone, thereby influencing the complex predisposition to hypertension.

Association Testing of Previously Reported Variants in a Large Case–Control Meta-Analysis of Diabetic Nephropathy

We formed the GEnetics of Nephropathy–an International Effort (GENIE) consortium to examine previously reported genetic associations with diabetic nephropathy (DN) in type 1 diabetes. GENIE consists of 6,366 similarly ascertained participants of European ancestry with type 1 diabetes, with and without DN, from the All Ireland-Warren 3-Genetics of Kidneys in Diabetes U.K. and Republic of Ireland (U.K.-R.O.I.) collection and the Finnish Diabetic Nephropathy Study (FinnDiane), combined with reanalyzed data from the Genetics of Kidneys in Diabetes U.S. Study (U.S. GoKinD). We found little evidence for the association of the EPO promoter polymorphism, rs161740, with the combined phenotype of proliferative retinopathy and end-stage renal disease in U.K.-R.O.I. (odds ratio [OR] 1.14, P = 0.19) or FinnDiane (OR 1.06, P = 0.60). However, a fixed-effects meta-analysis that included the previously reported cohorts retained a genome-wide significant association with that phenotype (OR 1.31, P = 2 × 10−9). An expanded investigation of the ELMO1 locus and genetic regions reported to be associated with DN in the U.S. GoKinD yielded only nominal statistical significance for these loci. Finally, top candidates identified in a recent meta-analysis failed to reach genome-wide significance. In conclusion, we were unable to replicate most of the previously reported genetic associations for DN, and significance for the EPO promoter association was attenuated.