Genetic diversity fuels gene discovery for tobacco and alcohol use.

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

Genome-wide association study identifies four pan-ancestry loci for suicidal ideation in the Million Veteran Program.

Suicidal ideation (SI) often precedes and predicts suicide attempt and death, is the most common suicidal phenotype and is over-represented in veterans. The genetic architecture of SI in the absence of suicide attempt (SA) is unknown, yet believed to have distinct and overlapping risk with other suicidal behaviors. We performed the first GWAS of SI without SA in the Million Veteran Program (MVP), identifying 99,814 SI cases from electronic health records without a history of SA or suicide death (SD) and 512,567 controls without SI, SA or SD. GWAS was performed separately in the four largest ancestry groups, controlling for sex, age and genetic substructure. Ancestry-specific results were combined via meta-analysis to identify pan-ancestry loci. Four genome-wide significant (GWS) loci were identified in the pan-ancestry meta-analysis with loci on chromosomes 6 and 9 associated with suicide attempt in an independent sample. Pan-ancestry gene-based analysis identified GWS associations with DRD2, DCC, FBXL19, BCL7C, CTF1, ANNK1, and EXD3. Gene-set analysis implicated synaptic and startle response pathways (q's<0.05). European ancestry (EA) analysis identified GWS loci on chromosomes 6 and 9, as well as GWS gene associations in EXD3, DRD2, and DCC. No other ancestry-specific GWS results were identified, underscoring the need to increase representation of diverse individuals. The genetic correlation of SI and SA within MVP was high (rG = 0.87; p = 1.09e-50), as well as with post-traumatic stress disorder (PTSD; rG = 0.78; p = 1.98e-95) and major depressive disorder (MDD; rG = 0.78; p = 8.33e-83). Conditional analysis on PTSD and MDD attenuated most pan-ancestry and EA GWS signals for SI without SA to nominal significance, with the exception of EXD3 which remained GWS. Our novel findings support a polygenic and complex architecture for SI without SA which is largely shared with SA and overlaps with psychiatric conditions frequently comorbid with suicidal behaviors.

Lysyl oxidase-like 1-antisense 1 (LOXL1-AS1) lncRNA differentially regulates gene and protein expression, signaling and morphology of human ocular cells.

Pseudoexfoliation glaucoma (PEXG) is characterized by dysregulated extracellular matrix (ECM) homeostasis that disrupts conventional outflow function and increases intraocular pressure (IOP). Prolonged IOP elevation results in optic nerve head damage and vision loss. Uniquely, PEXG is a form of open angle glaucoma that has variable penetrance, is difficult to treat and does not respond well to common IOP-lowering pharmaceuticals. Therefore, understanding modulators of disease severity will aid in targeted therapies for PEXG. Genome-wide association studies have identified polymorphisms in the long non-coding RNA lysyl oxidase-like 1-antisense 1 (LOXL1-AS1) as a risk factor for PEXG. Risk alleles, oxidative stress and mechanical stretch all alter LOXL1-AS1 expression. As a long non-coding RNA, LOXL1-AS1 binds hnRNPL and regulates global gene expression. In this study, we focus on the role of LOXL1-AS1 in the ocular cells (trabecular meshwork and Schlemm's canal) that regulate IOP. We show that selective knockdown of LOXL1-AS1 leads to cell-type-specific changes in gene expression, ECM homeostasis, signaling and morphology. These results implicate LOXL1-AS1 as a modulator of cellular homeostasis, altering cell contractility and ECM turnover, both of which are well-known contributors to PEXG. These findings support LOXL1-AS1 as a key target for modifying the disease.

Nontargeted Plasma Proteomic Analysis of Renal Disease and Pulmonary Hypertension in Patients with Sickle Cell Disease.

Sickle cell disease (SCD) is characterized by red blood cell sickling, vaso-occlusion, hemolytic anemia, damage to multiple organ systems, and, as a result, shortened life expectancy. Sickle cell disease nephropathy (SCDN) and pulmonary hypertension (pHTN) are common and frequently co-occurring complications of SCD; both are associated with markedly accelerated mortality. To identify candidate circulating biomarkers of SCDN and pHTN, we used mass spectrometry to quantify the relative abundance of >1000 proteins in plasma samples from 189 adults with SCD from the Outcome Modifying Genes in SCD (OMG-SCD) cohort (ProteomeXchange identifier PXD048716). Forty-four proteins were differentially abundant in SCDN, most significantly cystatin-C and collagen α-1(XVIII) chain (COIA1), and 55 proteins were dysregulated in patients with SCDN and pHTN, most significantly insulin-like growth factor-binding protein 6 (IBP6). Network analysis identified a module of 133 coregulated proteins significantly associated with SCDN, that was enriched for extracellular matrix proteins, insulin-like growth factor binding proteins, cell adhesion proteins, EGF-like calcium binding proteins, and several cadherin family members. Collectively, these data provide a comprehensive understanding of plasma protein changes in SCDN and pHTN which validate numerous studies of chronic kidney disease and suggest shared profiles of protein disruption in kidney dysfunction and pHTN among SCD patients.

Variation and impact of polygenic hematologic traits in monogenic sickle cell disease.

Several of the complications observed in sickle cell disease (SCD) are influenced by variation in hematologic traits (HT), such as fetal hemoglobin (HbF) level and neutrophil count. Previous large-scale genome-wide association studies carried out in largely healthy individuals have identified thousands of variants associated with HT, which have then been used to develop multi-ancestry polygenic trait scores (PTS). Here, we tested whether these PTS associate with HT in SCD patients and if they can improve statistical models associated with SCD-related complications. In 2,056 SCD patients, we found that the PTS predicted less HT variance than in non-SCD individuals of African ancestry. This was particularly striking at the Duffy/DARC locus, where we observed an epistatic interaction between the SCD genotype and the Duffy null variant (rs2814778) that led to a two-fold weaker effect on neutrophil count. PTS for these HT which are measured as part of routine practice were not associated with complications in SCD. In contrast, we found that a simple PTS for HbF that includes only six variants explained a large fraction of the phenotypic variation (20.5-27.1%), associated with acute chest syndrome and stroke risk, and improved the statistical modeling of the vaso-occlusive crisis rate. Using Mendelian randomization, we found that increasing HbF by 4.8% reduces stroke risk by 39% (P=0.0006). Taken together, our results highlight the importance of validating PTS in large diseased populations before proposing their implementation in the context of precision medicine initiatives.

A genome-wide association study of suicide attempts in the million veterans program identifies evidence of pan-ancestry and ancestry-specific risk loci.

To identify pan-ancestry and ancestry-specific loci associated with attempting suicide among veterans, we conducted a genome-wide association study (GWAS) of suicide attempts within a large, multi-ancestry cohort of U.S. veterans enrolled in the Million Veterans Program (MVP). Cases were defined as veterans with a documented history of suicide attempts in the electronic health record (EHR; N = 14,089) and controls were defined as veterans with no documented history of suicidal thoughts or behaviors in the EHR (N = 395,064). GWAS was performed separately in each ancestry group, controlling for sex, age and genetic substructure. Pan-ancestry risk loci were identified through meta-analysis and included two genome-wide significant loci on chromosomes 20 (p = 3.64 × 10-9) and 1 (p = 3.69 × 10-8). A strong pan-ancestry signal at the Dopamine Receptor D2 locus (p = 1.77 × 10-7) was also identified and subsequently replicated in a large, independent international civilian cohort (p = 7.97 × 10-4). Additionally, ancestry-specific genome-wide significant loci were also detected in African-Americans, European-Americans, Asian-Americans, and Hispanic-Americans. Pathway analyses suggested over-representation of many biological pathways with high clinical significance, including oxytocin signaling, glutamatergic synapse, cortisol synthesis and secretion, dopaminergic synapse, and circadian rhythm. These findings confirm that the genetic architecture underlying suicide attempt risk is complex and includes both pan-ancestry and ancestry-specific risk loci. Moreover, pathway analyses suggested many commonly impacted biological pathways that could inform development of improved therapeutics for suicide prevention.

A Genocentric Approach to Discovery of Mendelian Disorders.

The advent of inexpensive, clinical exome sequencing (ES) has led to the accumulation of genetic data from thousands of samples from individuals affected with a wide range of diseases, but for whom the underlying genetic and molecular etiology of their clinical phenotype remains unknown. In many cases, detailed phenotypes are unavailable or poorly recorded and there is little family history to guide study. To accelerate discovery, we integrated ES data from 18,696 individuals referred for suspected Mendelian disease, together with relatives, in an Apache Hadoop data lake (Hadoop Architecture Lake of Exomes [HARLEE]) and implemented a genocentric analysis that rapidly identified 154 genes harboring variants suspected to cause Mendelian disorders. The approach did not rely on case-specific phenotypic classifications but was driven by optimization of gene- and variant-level filter parameters utilizing historical Mendelian disease-gene association discovery data. Variants in 19 of the 154 candidate genes were subsequently reported as causative of a Mendelian trait and additional data support the association of all other candidate genes with disease endpoints.

Potential causal role of l-glutamine in sickle cell disease painful crises: A Mendelian randomization analysis.

In a recent clinical trial, the metabolite l-glutamine was shown to reduce painful crises in sickle cell disease (SCD) patients. To support this observation and identify other metabolites implicated in SCD clinical heterogeneity, we profiled 129 metabolites in the plasma of 705 SCD patients. We tested correlations between metabolite levels and six SCD-related complications (painful crises, cholecystectomy, retinopathy, leg ulcer, priapism, aseptic necrosis) or estimated glomerular filtration rate (eGFR), and used Mendelian randomization (MR) to assess causality. We found a potential causal relationship between l-glutamine levels and painful crises (N = 1278, odds ratio (OR) [95% confidence interval] = 0.68 [0.52-0.89], P = 0.0048). In two smaller SCD cohorts (N = 299 and 406), the protective effect of l-glutamine was observed (OR = 0.82 [0.50-1.34]), although the MR result was not significant (P = 0.44). We identified 66 significant correlations between the levels of other metabolites and SCD-related complications or eGFR. We tested these correlations for causality using MR analyses and found no significant causal relationship. The baseline levels of quinolinic acid were associated with prospectively ascertained survival in SCD patients, and this effect was dependent on eGFR. Metabolomics provide a promising approach to prioritize small molecules that may serve as biomarkers or drug targets in SCD.

Rapid decline in estimated glomerular filtration rate in sickle cell anemia: results of a multicenter pooled analysis.

Not available.

Genetic variation in dopamine neurotransmission and motor development of infants born extremely-low-birthweight.

AIM: To determine if genetic variation associated with decreased dopamine neurotransmission predicts a decrease in motor development in a convenience cohort study of infants born extremely-low-birthweight (ELBW). METHOD: Four hundred and ninety-eight infants born ELBW had genome-wide genotyping and a neurodevelopmental evaluation at 18 to 22 months of age, corrected for preterm birth. A polygenic risk score (PRS) was created to combine into one predictor variable the hypothesized influences on motor development of alleles at seven independent single nucleotide polymorphisms previously associated with relative decreases in both dopamine neurotransmission and motor learning, by summing the number of alleles present in each infant (range=0-14). The motor development outcome was the Psychomotor Development Index (PDI) of the Bayley Scales of Infant Development, Second Edition. The linear regression models were adjusted for seven clinical and four genetic ancestry covariates. The mean PRS of infants with cerebral palsy (CP) was compared to those without CP. RESULTS: PRS was inversely related to PDI (p=0.011). Each 1-point increase in PRS resulted in an average decrease in PDI of 1.37 points. Patients with CP did not have a greater mean PRS than those without (p=0.67), both with and without adjustment for covariates. INTERPRETATION: Genetic variation that favors a decrease in dopamine neurotransmission predisposes to a decrease in motor development in infants born ELBW, but not to the diagnosis of CP. WHAT THIS PAPER ADDS: Genetic variation in dopamine neurotransmission was associated with a decrease in motor development in infants born at an extremely-low-birthweight. It does not predispose to the diagnosis of cerebral palsy.

Clinical and metabolomic risk factors associated with rapid renal function decline in sickle cell disease.

Sickle cell disease (SCD) nephropathy and lower estimated glomerular filtration rate (eGFR) are risk factors for early mortality. Furthermore, rate of eGFR decline predicts progression to end-stage renal disease in many clinical settings. However, factors predicting renal function decline in SCD are poorly documented. Using clinical, laboratory, genetic, and metabolomic data, we evaluated predictors of renal function decline in a longitudinal cohort of 288 adults (mean age 33.0 years). In 193 subjects with 5-year follow-up data, mean rate of eGFR decline was 2.35 mL/min/1.73 m2 /year, nearly twice that of African American adults overall. Hyperfiltration was prevalent at baseline (61.1%), and 36.8% of subjects experienced rapid eGFR decline (≥3 mL/min/1.73 m2 /year). Severe Hb genotype; proteinuria; higher platelet and reticulocyte counts, and systolic BP; and lower Hb level and BMI were associated with rapid decline. A risk scoring system was created using these 7 variables and was highly predictive of rapid eGFR decline, with odds of rapid decline increasing 1.635-fold for every point increment (P < 0.0001). Rapid eGFR decline was also associated with higher organ system severity score and peak creatinine. Additionally, two metabolites (asymmetric dimethylarginine and quinolinic acid) were associated with rapid decline. Further investigation into longitudinal SCD nephropathy (SCDN) trajectory, early markers of SCDN, and tools for risk stratification should inform interventional studies targeted to slowing GFR decline and improving SCD outcomes.

In vivo Modeling Implicates APOL1 in Nephropathy: Evidence for Dominant Negative Effects and Epistasis under Anemic Stress.

African Americans have a disproportionate risk for developing nephropathy. This disparity has been attributed to coding variants (G1 and G2) in apolipoprotein L1 (APOL1); however, there is little functional evidence supporting the role of this protein in renal function. Here, we combined genetics and in vivo modeling to examine the role of apol1 in glomerular development and pronephric filtration and to test the pathogenic potential of APOL1 G1 and G2. Translational suppression or CRISPR/Cas9 genome editing of apol1 in zebrafish embryos results in podocyte loss and glomerular filtration defects. Complementation of apol1 morphants with wild-type human APOL1 mRNA rescues these defects. However, the APOL1 G1 risk allele does not ameliorate defects caused by apol1 suppression and the pathogenicity is conferred by the cis effect of both individual variants of the G1 risk haplotype (I384M/S342G). In vivo complementation studies of the G2 risk allele also indicate that the variant is deleterious to protein function. Moreover, APOL1 G2, but not G1, expression alone promotes developmental kidney defects, suggesting a possible dominant-negative effect of the altered protein. In sickle cell disease (SCD) patients, we reported previously a genetic interaction between APOL1 and MYH9. Testing this interaction in vivo by co-suppressing both transcripts yielded no additive effects. However, upon genetic or chemical induction of anemia, we observed a significantly exacerbated nephropathy phenotype. Furthermore, concordant with the genetic interaction observed in SCD patients, APOL1 G2 reduces myh9 expression in vivo, suggesting a possible interaction between the altered APOL1 and myh9. Our data indicate a critical role for APOL1 in renal function that is compromised by nephropathy-risk encoding variants. Moreover, our interaction studies indicate that the MYH9 locus is also relevant to the phenotype in a stressed microenvironment and suggest that consideration of the context-dependent functions of both proteins will be required to develop therapeutic paradigms.

Genetic variants and cellular stressors associated with exfoliation syndrome modulate promoter activity of a lncRNA within the LOXL1 locus.

Exfoliation syndrome (XFS) is a common, age-related, systemic fibrillinopathy. It greatly increases risk of exfoliation glaucoma (XFG), a major worldwide cause of irreversible blindness. Coding variants in the lysyl oxidase-like 1 (LOXL1) gene are strongly associated with XFS in all studied populations, but a functional role for these variants has not been established. To identify additional candidate functional variants, we sequenced the entire LOXL1 genomic locus (∼40 kb) in 50 indigenous, black South African XFS cases and 50 matched controls. The variants with the strongest evidence of association were located in a well-defined 7-kb region bounded by the 3'-end of exon 1 and the adjacent region of intron 1 of LOXL1. We replicated this finding in US Caucasian (91 cases/1031 controls), German (771 cases/1365 controls) and Japanese (1484 cases/1188 controls) populations. The region of peak association lies upstream of LOXL1-AS1, a long non-coding RNA (lncRNA) encoded on the opposite strand of LOXL1. We show that this region contains a promoter and, importantly, that the strongly associated XFS risk alleles in the South African population are functional variants that significantly modulate the activity of this promoter. LOXL1-AS1 expression is also significantly altered in response to oxidative stress in human lens epithelial cells and in response to cyclic mechanical stress in human Schlemm's canal endothelial cells. Taken together, these findings support a functional role for the LOXL1-AS1 lncRNA in cellular stress response and suggest that dysregulation of its expression by genetic risk variants plays a key role in XFS pathogenesis.

The association of single-nucleotide polymorphisms in the oxytocin receptor and G protein-coupled receptor kinase 6 (GRK6) genes with oxytocin dosing requirements and labor outcomes.

BACKGROUND: Oxytocin is a potent uterotonic agent that is widely used for induction and augmentation of labor. Oxytocin has a narrow therapeutic index and the optimal dosing for any individual woman varies widely. OBJECTIVE: The objective of this study was to determine whether genetic variation in the oxytocin receptor (OXTR) or in the gene encoding G protein-coupled receptor kinase 6 (GRK6), which regulates desensitization of the oxytocin receptor, could explain variation in oxytocin dosing and labor outcomes among women being induced near term. STUDY DESIGN: Pregnant women with a singleton gestation residing in Durham County, NC, were prospectively enrolled as part of the Healthy Pregnancy, Healthy Baby cohort study. Those women undergoing an induction of labor at 36 weeks or greater were genotyped for 18 haplotype-tagging single-nucleotide polymorphisms in OXTR and 7 haplotype-tagging single-nucleotide polymorphisms in GRK6 using TaqMan assays. Linear regression was used to examine the relationship between maternal genotype and maximal oxytocin infusion rate, total oxytocin dose received, and duration of labor. Logistic regression was used to test for the association of maternal genotype with mode of delivery. For each outcome, backward selection techniques were utilized to control for important confounding variables and additive genetic models were used. Race/ethnicity was included in all models because of differences in allele frequencies across populations, and Bonferroni correction for multiple testing was used. RESULTS: DNA was available from 482 women undergoing induction of labor at 36 weeks or greater. Eighteen haplotype-tagging single-nucleotide polymorphisms within OXTR and 7 haplotype-tagging single-nucleotide polymorphisms within GRK6 were examined. Five single-nucleotide polymorphisms in OXTR showed nominal significance with maximal infusion rate of oxytocin, and two single-nucleotide polymorphisms in OXTR were associated with total oxytocin dose received. One single-nucleotide polymorphism in OXTR and two single-nucleotide polymorphisms in GRK6 were associated with duration of labor, one of which met the multiple testing threshold (P = .0014, rs2731664 [GRK6], mean duration of labor, 17.7 hours vs 20.2 hours vs 23.5 hours for AA, AC, and CC genotypes, respectively). Three single-nucleotide polymorphisms, two in OXTR and one in GRK6, showed nominal significance with mode of delivery. CONCLUSION: Genetic variation in OXTR and GRK6 is associated with the amount of oxytocin required as well as the duration of labor and risk for cesarean delivery among women undergoing induction of labor near term. With further research, pharmacogenomic approaches may potentially be utilized to develop personalized treatment to improve safety and efficacy outcomes among women undergoing induction of labor.

Epigenome-wide association of PTSD from heterogeneous cohorts with a common multi-site analysis pipeline.

Compelling evidence suggests that epigenetic mechanisms such as DNA methylation play a role in stress regulation and in the etiologic basis of stress related disorders such as Post traumatic Stress Disorder (PTSD). Here we describe the purpose and methods of an international consortium that was developed to study the role of epigenetics in PTSD. Inspired by the approach used in the Psychiatric Genomics Consortium, we brought together investigators representing seven cohorts with a collective sample size of N = 1147 that included detailed information on trauma exposure, PTSD symptoms, and genome-wide DNA methylation data. The objective of this consortium is to increase the analytical sample size by pooling data and combining expertise so that DNA methylation patterns associated with PTSD can be identified. Several quality control and analytical pipelines were evaluated for their control of genomic inflation and technical artifacts with a joint analysis procedure established to derive comparable data over the cohorts for meta-analysis. We propose methods to deal with ancestry population stratification and type I error inflation and discuss the advantages and disadvantages of applying robust error estimates. To evaluate our pipeline, we report results from an epigenome-wide association study (EWAS) of age, which is a well-characterized phenotype with known epigenetic associations. Overall, while EWAS are highly complex and subject to similar challenges as genome-wide association studies (GWAS), we demonstrate that an epigenetic meta-analysis with a relatively modest sample size can be well-powered to identify epigenetic associations. Our pipeline can be used as a framework for consortium efforts for EWAS.

Linkage of familial essential tremor to chromosome 5q35.

BACKGROUND: Essential tremor is a neurological condition characterized by tremor during voluntary movement. To date, 3 loci linked to familial essential tremor have been identified. METHODS: We examined 48 essential tremor patients in 5 large essential tremor pedigrees in our data set for genetic linkage using an Affymetrix Axiom array. Linkage analysis was performed using an affecteds-only dominant model in SIMWALK2. To incorporate all genotype information, GERMLINE was used to identify genome segments shared identical-by-descent in pairs of affecteds. Exome sequencing was performed in pedigrees showing evidence of linkage. RESULTS: For one family, chromosomes 5 and 18 showed genome-wide significant linkage to essential tremor. Shared segment analysis excluded the 18p11 candidate region and reduced the 5q35 region by 1 megabase. Exome sequencing did not identify a potential causative variant in this region. CONCLUSION: A locus on chromosome 5 is linked to essential tremor. Further research is needed to identify a causative variant. © 2016 International Parkinson and Movement Disorder Society.

Joint eQTL assessment of whole blood and dura mater tissue from individuals with Chiari type I malformation.

BACKGROUND: Expression quantitative trait loci (eQTL) play an important role in the regulation of gene expression. Gene expression levels and eQTLs are expected to vary from tissue to tissue, and therefore multi-tissue analyses are necessary to fully understand complex genetic conditions in humans. Dura mater tissue likely interacts with cranial bone growth and thus may play a role in the etiology of Chiari Type I Malformation (CMI) and related conditions, but it is often inaccessible and its gene expression has not been well studied. A genetic basis to CMI has been established; however, the specific genetic risk factors are not well characterized. RESULTS: We present an assessment of eQTLs for whole blood and dura mater tissue from individuals with CMI. A joint-tissue analysis identified 239 eQTLs in either dura or blood, with 79% of these eQTLs shared by both tissues. Several identified eQTLs were novel and these implicate genes involved in bone development (IPO8, XYLT1, and PRKAR1A), and ribosomal pathways related to marrow and bone dysfunction, as potential candidates in the development of CMI. CONCLUSIONS: Despite strong overall heterogeneity in expression levels between blood and dura, the majority of cis-eQTLs are shared by both tissues. The power to detect shared eQTLs was improved by using an integrative statistical approach. The identified tissue-specific and shared eQTLs provide new insight into the genetic basis for CMI and related conditions.

Hepatic gene expression profiles differentiate presymptomatic patients with mild versus severe nonalcoholic fatty liver disease.

UNLABELLED: Clinicians rely upon the severity of liver fibrosis to segregate patients with well-compensated nonalcoholic fatty liver disease (NAFLD) into subpopulations at high- versus low-risk for eventual liver-related morbidity and mortality. We compared hepatic gene expression profiles in high- and low-risk NAFLD patients to identify processes that distinguish the two groups and hence might be novel biomarkers or treatment targets. Microarray analysis was used to characterize gene expression in percutaneous liver biopsies from low-risk, "mild" NAFLD patients (fibrosis stage 0-1; n = 40) and high-risk, "severe" NAFLD patients (fibrosis stage 3-4; n = 32). Findings were validated in a second, independent cohort and confirmed by real-time polymerase chain reaction and immunohistochemistry (IHC). As a group, patients at risk for bad NAFLD outcomes had significantly worse liver injury and more advanced fibrosis (severe NAFLD) than clinically indistinguishable NAFLD patients with a good prognosis (mild NAFLD). A 64-gene profile reproducibly differentiated severe NAFLD from mild NAFLD, and a 20-gene subset within this profile correlated with NAFLD severity, independent of other factors known to influence NAFLD progression. Multiple genes involved with tissue repair/regeneration and certain metabolism-related genes were induced in severe NAFLD. Ingenuity Pathway Analysis and IHC confirmed deregulation of metabolic and regenerative pathways in severe NAFLD and revealed overlap among the gene expression patterns of severe NAFLD, cardiovascular disease, and cancer. CONCLUSION: By demonstrating specific metabolic and repair pathways that are differentially activated in livers with severe NAFLD, gene profiling identified novel targets that can be exploited to improve diagnosis and treatment of patients who are at greatest risk for NAFLD-related morbidity and mortality.

Gene-centric association study of acute chest syndrome and painful crisis in sickle cell disease patients.

Patients with sickle cell disease (SCD) present with a wide range of clinical complications. Understanding this clinical heterogeneity offers the prospects to tailor the right treatments to the right patients and also guide the development of novel therapies. Several environmental (eg, nutrition) and nonenvironmental (eg, fetal hemoglobin levels, α-thalassemia status) factors are known to modify SCD severity. To find new genetic modifiers of SCD severity, we performed a gene-centric association study in 1514 African American participants from the Cooperative Study of Sickle Cell Disease (CSSCD) for acute chest syndrome (ACS) and painful crisis. From the initial results, we selected 36 single nucleotide polymorphism (SNPs) and genotyped them for replication in 387 independent patients from the CSSCD, 318 SCD patients recruited at Georgia Health Sciences University, and 449 patients from the Duke SCD cohort. In the combined analysis, an association between ACS and rs6141803 reached array-wide significance (P = 4.1 × 10(-7)). This SNP is located 8.2 kilobases upstream of COMMD7, a gene highly expressed in the lung that interacts with nuclear factor-κB signaling. Our results provide new leads to gaining a better understanding of clinical variability in SCD, a "simple" monogenic disease.