Identifying X-Chromosome Variants Associated with Age-Related Macular Degeneration

Purpose

In genome-wide association studies (GWAS), X chromosome (ChrX) variants are often not investigated. Sex-specific effects and ChrX-specific quality control (QC) are needed to examine these effects. Previous analyses identified 52 autosomal variants associated with age-related macular degeneration (AMD) via the International AMD Genomics Consortium (IAMDGC), but did not analyze ChrX. Therefore, our goal was to investigate ChrX variants for association with AMD.

Methods

We genotyped 29,629 non-Hispanic White (NHW) individuals (M/F:10,404/18,865; AMD12,087/14723) via a custom chip and imputed after ChrX-specific QC (XWAS 3.0) using the Michigan Imputation Server. Imputation generated 1,221,623 variants on ChrX. Age, informative PCs, and subphenotyeps were covariates for logistic association analyses with Fishers correction. Gene/pathway analyses were performed with VEGAS, GSEASNP, ICSNPathway, DAVID, and mirPath.

Results

Logistic association on NHW individuals with sex correction, identified variants in/near the genes SLITRK4, ARHGAP6, FGF13 and DMD associated with AMD (P<1x10 -6 ,Fishers combined-corrected). Via association testing of subphenotypes of choroidal neovascularization and geographic atrophy (GA), variants in DMD associated with GA (P<1x10 -6 , Fishers combined-corrected). Via gene-based analysis with VEGAS, several genes were associated with AMD (P<0.05, both truncated tail strength/truncated product P) including SLITRK4 and BHLHB9 . Pathway analysis using GSEASNP and DAVID showed genes associated with nervous system development (FDR: P:0.02), and blood coagulation (FDR: P:0.03). Variants in the region of a microRNA (miR) were associated with AMD (P<0.05, truncated tail strength/truncated product P). Via DIANA mirPath analysis, downstream targets of miRs show association with brain disorders and fatty acid elongation (P<0.05). A long-non coding RNA on ChrX near the DMD locus was also associated with AMD (P=4x10 -7 ). Epistatic analysis (t-statistic) for a quantitative trait of AMD vs control including covariates found a suggestive association in the XG gene (P=2x10^-5).

Conclusions

Analysis of ChrX variants demonstrates association with AMD and these variants may be linked to novel pathways. Further analysis is needed to confirm results and to understand their biological significance and relationship with AMD development in worldwide populations.

Genomic Disorders in CKD across the Lifespan

SIGNIFICANCE STATEMENT

Pathogenic structural genetic variants, also known as genomic disorders, have been associated with pediatric CKD. This study extends those results across the lifespan, with genomic disorders enriched in both pediatric and adult patients compared with controls. In the Chronic Renal Insufficiency Cohort study, genomic disorders were also associated with lower serum Mg, lower educational performance, and a higher risk of death. A phenome-wide association study confirmed the link between kidney disease and genomic disorders in an unbiased way. Systematic detection of genomic disorders can provide a molecular diagnosis and refine prediction of risk and prognosis.

BACKGROUND

Genomic disorders (GDs) are associated with many comorbid outcomes, including CKD. Identification of GDs has diagnostic utility.

METHODS

We examined the prevalence of GDs among participants in the Chronic Kidney Disease in Children (CKiD) cohort II ( n =248), Chronic Renal Insufficiency Cohort (CRIC) study ( n =3375), Columbia University CKD Biobank (CU-CKD; n =1986), and the Family Investigation of Nephropathy and Diabetes (FIND; n =1318) compared with 30,746 controls. We also performed a phenome-wide association analysis (PheWAS) of GDs in the electronic MEdical Records and GEnomics (eMERGE; n =11,146) cohort.

RESULTS

We found nine out of 248 (3.6%) CKiD II participants carried a GD, replicating prior findings in pediatric CKD. We also identified GDs in 72 out of 6679 (1.1%) adult patients with CKD in the CRIC, CU-CKD, and FIND cohorts, compared with 199 out of 30,746 (0.65%) GDs in controls (OR, 1.7; 95% CI, 1.3 to 2.2). Among adults with CKD, we found recurrent GDs at the 1q21.1, 16p11.2, 17q12, and 22q11.2 loci. The 17q12 GD (diagnostic of renal cyst and diabetes syndrome) was most frequent, present in 1:252 patients with CKD and diabetes. In the PheWAS, dialysis and neuropsychiatric phenotypes were the top associations with GDs. In CRIC participants, GDs were associated with lower serum magnesium, lower educational achievement, and higher mortality risk.

CONCLUSION

Undiagnosed GDs are detected both in children and adults with CKD. Identification of GDs in these patients can enable a precise genetic diagnosis, inform prognosis, and help stratify risk in clinical studies. GDs could also provide a molecular explanation for nephropathy and comorbidities, such as poorer neurocognition for a subset of patients.

Younger Age and Albuminuria are Associated with Proliferative Diabetic Retinopathy and Diabetic Macular Edema in the South Indian GeNetics of DiAbeTic Retinopathy (SIGNATR) Study

Purpose: The purpose of the South Indian GeNetics of DiAbeTic Retinopathy (SIGNATR) Study is to identify non-genetic and genetic risk factors associated with diabetic retinopathy (DR). This report examines the non-genetic risk factors for DR in South Indian patients.Methods: Participants with South Indian ancestry and type 2 diabetes (T2D) were included from two sources: the Sankara Nethralaya Diabetic Retinopathy and Molecular Genetics Study (SN-DREAMS) and prospective recruitment at Sankara Nethralaya affiliates. Fundus photography and optical coherence tomography (OCT) were obtained on participants. Fundus images were graded for DR severity and OCTs were graded for center-involved diabetic macular edema (ciDME). Multivariate analyses were performed using stepwise logistic regression to assess effects of the demographic and clinical factors on proliferative DR (PDR) and DME.Results: Among the 2941 participants with DR grading, participants with PDR were more likely to be younger [odds ratio (OR)=0.95], men (OR = 1.83), have a longer duration of diabetes (OR = 1.10), have a higher hemoglobin A1c (OR = 1.12), have albuminuria (OR = 5.83), have hypertension (OR = 1.69), have a higher HDL (OR = 1.02) and a lower total cholesterol (OR = 0.99) (all p < 0.05). Among the 483 participants with gradable OCT scans, participants who had ciDME were more likely to be younger (OR = 0.97), men (OR = 2.80), have a longer duration of diabetes (OR = 1.06), have lower triglycerides (OR = 0.99), and have albuminuria (OR = 3.12) (all p < 0.05).Conclusions: Younger age, male sex, longer duration of diabetes, higher HbA1c, and presence of albuminuria were identified as risk factors for PDR and DME in a South Indian population with T2D.

Gene Set Enrichment Analsyes Identify Pathways Involved in Genetic Risk for Diabetic Retinopathy

To identify functionally related genes associated with diabetic retinopathy (DR) risk using gene set enrichment analyses applied to genome-wide association study meta-analyses.

METHODS

We analyzed DR GWAS meta-analyses performed on 3246 Europeans and 2611 African Americans with type 2 diabetes. Gene sets relevant to 5 key DR pathophysiology processes were investigated: tissue injury, vascular events, metabolic events and glial dysregulation, neuronal dysfunction, and inflammation. Keywords relevant to these processes were queried in 4 pathway and ontology databases. Two GSEA methods, Meta-Analysis Gene set Enrichment of variaNT Associations (MAGENTA) and Multi-marker Analysis of GenoMic Annotation (MAGMA), were used. Gene sets were defined to be enriched for gene associations with DR if the P value corrected for multiple testing (Pcorr) was <.05.

RESULTS

Five gene sets were significantly enriched for numerous modest genetic associations with DR in one method (MAGENTA or MAGMA) and also at least nominally significant (uncorrected P < .05) in the other method. These pathways were regulation of the lipid catabolic process (2-fold enrichment, Pcorr = .014); nitric oxide biosynthesis (1.92-fold enrichment, Pcorr = .022); lipid digestion, mobilization, and transport (1.6-fold enrichment, P = .032); apoptosis (1.53-fold enrichment, P = .041); and retinal ganglion cell degeneration (2-fold enrichment, Pcorr = .049). The interferon gamma (IFNG) gene, previously implicated in DR by protein-protein interactions in our GWAS, was among the top ranked genes in the nitric oxide pathway (best variant P = .0001).

CONCLUSIONS

These GSEA indicate that variants in genes involved in oxidative stress, lipid transport and catabolism, and cell degeneration are enriched for genes associated with DR risk. NOTE: Publication of this article is sponsored by the American Ophthalmological Society.

Methylome-wide Analysis Reveals Epigenetic Marks Associated With Resistance to Tuberculosis in Human Immunodeficiency Virus-Infected Individuals From East Africa

BACKGROUND

Tuberculosis (TB) is the most deadly infectious disease globally and is highly prevalent in the developing world. For individuals infected with both Mycobacterium tuberculosis (Mtb) and human immunodeficiency virus (HIV), the risk of active TB is 10% or more annually. Previously, we identified in a genome-wide association study (GWAS) a region on chromosome 5 associated with resistance to TB, which included epigenetic marks that could influence gene regulation. We hypothesized that HIV-infected individuals exposed to Mtb who remain disease free carry epigenetic changes that strongly protect them from active TB.

METHODS

We conducted a methylome-wide study in HIV-infected, TB-exposed cohorts from Uganda and Tanzania and integrated data from our GWAS.

RESULTS

We identified 3 regions of interest that included markers that were differentially methylated between TB cases and controls with latent TB infection: chromosome 1 (RNF220, P = 4 × 10-5), chromosome 2 (between COPS8 and COL6A3, P = 2.7 × 10-5), and chromosome 5 (CEP72, P = 1.3 × 10-5). These methylation results co-localized with associated single-nucleotide polymorphisms (SNPs), methylation QTLs, and methylation × SNP interaction effects. These markers were in regions with regulatory markers for cells involved in TB immunity and/or lung.

CONCLUSIONS

Epigenetic regulation is a potential biologic factor underlying resistance to TB in immunocompromised individuals that can act in conjunction with genetic variants.

Association between genes regulating neural pathways for quantitative traits of speech and language disorders

Speech sound disorders (SSD) manifest as difficulties in phonological memory and awareness, oral motor function, language, vocabulary, reading, and spelling. Families enriched for SSD are rare, and typically display a cluster of deficits. We conducted a genome-wide association study (GWAS) in 435 children from 148 families in the Cleveland Family Speech and Reading study (CFSRS), examining 16 variables representing 6 domains. Replication was conducted using the Avon Longitudinal Study of Parents and Children (ALSPAC). We identified 18 significant loci (combined p < 10-8) that we pursued bioinformatically. We prioritized 5 novel gene regions with likely functional repercussions on neural pathways, including those which colocalized with differentially methylated regions in our sample. Polygenic risk scores for receptive language, expressive vocabulary, phonological awareness, phonological memory, spelling, and reading decoding associated with increasing clinical severity. In summary, neural-genetic influence on SSD is primarily multigenic and acts on genomic regulatory elements, similar to other neurodevelopmental disorders.

Genome-wide meta-analysis identifies 127 open-angle glaucoma loci with consistent effect across ancestries

Primary open-angle glaucoma (POAG), is a heritable common cause of blindness world-wide. To identify risk loci, we conduct a large multi-ethnic meta-analysis of genome-wide association studies on a total of 34,179 cases and 349,321 controls, identifying 44 previously unreported risk loci and confirming 83 loci that were previously known. The majority of loci have broadly consistent effects across European, Asian and African ancestries. Cross-ancestry data improve fine-mapping of causal variants for several loci. Integration of multiple lines of genetic evidence support the functional relevance of the identified POAG risk loci and highlight potential contributions of several genes to POAG pathogenesis, including SVEP1, RERE, VCAM1, ZNF638, CLIC5, SLC2A12, YAP1, MXRA5, and SMAD6. Several drug compounds targeting POAG risk genes may be potential glaucoma therapeutic candidates.

Association of Rare CYP39A1 Variants With Exfoliation Syndrome Involving the Anterior Chamber of the Eye

IMPORTANCE

Exfoliation syndrome is a systemic disorder characterized by progressive accumulation of abnormal fibrillar protein aggregates manifesting clinically in the anterior chamber of the eye. This disorder is the most commonly known cause of glaucoma and a major cause of irreversible blindness.

OBJECTIVE

To determine if exfoliation syndrome is associated with rare, protein-changing variants predicted to impair protein function.

DESIGN, SETTING, AND PARTICIPANTS

A 2-stage, case-control, whole-exome sequencing association study with a discovery cohort and 2 independently ascertained validation cohorts. Study participants from 14 countries were enrolled between February 1999 and December 2019. The date of last clinical follow-up was December 2019. Affected individuals had exfoliation material on anterior segment structures of at least 1 eye as visualized by slit lamp examination. Unaffected individuals had no signs of exfoliation syndrome.

EXPOSURES

Rare, coding-sequence genetic variants predicted to be damaging by bioinformatic algorithms trained to recognize alterations that impair protein function.

MAIN OUTCOMES AND MEASURES

The primary outcome was the presence of exfoliation syndrome. Exome-wide significance for detected variants was defined as P < 2.5 × 10-6. The secondary outcomes included biochemical enzymatic assays and gene expression analyses.

RESULTS

The discovery cohort included 4028 participants with exfoliation syndrome (median age, 78 years [interquartile range, 73-83 years]; 2377 [59.0%] women) and 5638 participants without exfoliation syndrome (median age, 72 years [interquartile range, 65-78 years]; 3159 [56.0%] women). In the discovery cohort, persons with exfoliation syndrome, compared with those without exfoliation syndrome, were significantly more likely to carry damaging CYP39A1 variants (1.3% vs 0.30%, respectively; odds ratio, 3.55 [95% CI, 2.07-6.10]; P = 6.1 × 10-7). This outcome was validated in 2 independent cohorts. The first validation cohort included 2337 individuals with exfoliation syndrome (median age, 74 years; 1132 women; n = 1934 with demographic data) and 2813 individuals without exfoliation syndrome (median age, 72 years; 1287 women; n = 2421 with demographic data). The second validation cohort included 1663 individuals with exfoliation syndrome (median age, 75 years; 587 women; n = 1064 with demographic data) and 3962 individuals without exfoliation syndrome (median age, 74 years; 951 women; n = 1555 with demographic data). Of the individuals from both validation cohorts, 5.2% with exfoliation syndrome carried CYP39A1 damaging alleles vs 3.1% without exfoliation syndrome (odds ratio, 1.82 [95% CI, 1.47-2.26]; P < .001). Biochemical assays classified 34 of 42 damaging CYP39A1 alleles as functionally deficient (median reduction in enzymatic activity compared with wild-type CYP39A1, 94.4% [interquartile range, 78.7%-98.2%] for the 34 deficient variants). CYP39A1 transcript expression was 47% lower (95% CI, 30%-64% lower; P < .001) in ciliary body tissues from individuals with exfoliation syndrome compared with individuals without exfoliation syndrome.

CONCLUSIONS AND RELEVANCE

In this whole-exome sequencing case-control study, presence of exfoliation syndrome was significantly associated with carriage of functionally deficient CYP39A1 sequence variants. Further research is needed to understand the clinical implications of these findings.

Association of APOE With Primary Open-Angle Glaucoma Suggests a Protective Effect for APOE ε4

Purpose

Prior studies have demonstrated that microglial activation is involved in the pathogenesis of primary open-angle glaucoma (POAG). Here we sought to identify genetic associations between POAG and variants in APOE and TREM2, genes associated with Alzheimer disease (AD) that critically regulate microglial neurodegeneration-associated molecular signature.

Methods

APOE genotypes were called using imputed data from the NEIGHBOR consortium (2120 POAG cases, 2262 controls) and a second cohort from the Massachusetts Eye and Ear Infirmary (MEEI; 486 cases, 344 controls). TREM2 coding variants were genotyped by means of the Illumina HumanExome BeadArray. The data set was analyzed for association with POAG overall, as well as the high-tension glaucoma (HTG) and normal-tension glaucoma (NTG) subgroups, using logistic regression adjusting for age and sex.

Results

In the combined NEIGHBOR-MEEI data set, significant association was observed for APOE ε4 in POAG overall (odds ratio [OR], 0.83; 95% confidence interval [CI], 0.74–0.94; P = 0.0022) and in both the HTG subgroup (OR, 0.81; 95% CI, 0.70–0.94; P = 0.0052) and NTG subgroup (OR, 0.71; 95% CI, 0.58–0.87; P = 0.0014). A rare TREM2 variant (A105V) was found only in HTG cases (3 of 2863 cases) and in none of the controls (P = 0.03). Three TREM2 rare variants associated with AD were not significantly associated with POAG (P > 0.05).

Conclusions

We have found that the APOE ε4 allele is associated with a reduced risk of POAG. Interestingly, the same allele is adversely associated with AD, suggesting a mechanistic difference between neurodegenerative diseases of the eye and the brain. TREM2 variants associated with AD did not significantly contribute to POAG risk.

Statistical driver genes as a means to uncover missing heritability for age-related macular degeneration

Background

Age-related macular degeneration (AMD) is a progressive retinal disease contributing to blindness worldwide. Multiple estimates for AMD heritability (h²) exist; however, a substantial proportion of h² is not attributable to known genomic loci. The International AMD Genomics Consortium (IAMDGC) gathered the largest dataset of advanced AMD (ADV) cases and controls available and identified 34 loci containing 52 independent risk variants defining known AMD h². To better define AMD heterogeneity, we used Pathway Analysis by Randomization Incorporating Structure (PARIS) on the IAMDGC data and identified 8 statistical driver genes (SDGs), including 2 novel SDGs not discovered by the IAMDGC. We chose to further investigate these pathway-based risk genes and determine their contribution to ADV h², as well as the differential ADV subtype h².

Methods

We performed genomic-relatedness-based restricted maximum-likelihood (GREML) analyses on ADV, geographic atrophy (GA), and choroidal neovascularization (CNV) subtypes to investigate the h² of genotyped variants on the full DNA array chip, 34 risk loci (n = 2758 common variants), 52 variants from the IAMDGC 2016 GWAS, and the 8 SDGs, specifically the novel 2 SDGs, PPARA and PLCG2.

Results

Via GREML, full chip h² was 44.05% for ADV, 46.37% for GA, and 62.03% for CNV. The lead 52 variants’ h² (ADV: 14.52%, GA: 8.02%, CNV: 13.62%) and 34 loci h² (ADV: 13.73%, GA: 8.81%, CNV: 12.89%) indicate that known variants contribute ~ 14% to ADV h². SDG variants account for a small percentage of ADV, GA, and CNV heritability, but estimates based on the combination of SDGs and the 34 known loci are similar to those calculated for known loci alone. We identified modest epistatic interactions among variants in the 2 SDGs and the 52 IAMDGC variants, including modest interactions between variants in PPARA and PLCG2.

Conclusions

Pathway analyses, which leverage biological relationships among genes in a pathway, may be useful in identifying additional loci that contribute to the heritability of complex disorders in a non-additive manner. Heritability analyses of these loci, especially amongst disease subtypes, may provide clues to the importance of specific genes to the genetic architecture of AMD.

Interaction between host genes and Mycobacterium tuberculosis lineage can affect tuberculosis severity: Evidence for coevolution?

Genetic studies of both the human host and Mycobacterium tuberculosis (MTB) demonstrate independent association with tuberculosis (TB) risk. However, neither explains a large portion of disease risk or severity. Based on studies in other infectious diseases and animal models of TB, we hypothesized that the genomes of the two interact to modulate risk of developing active TB or increasing the severity of disease, when present. We examined this hypothesis in our TB household contact study in Kampala, Uganda, in which there were 3 MTB lineages of which L4-Ugandan (L4.6) is the most recent. TB severity, measured using the Bandim TBscore, was modeled as a function of host SNP genotype, MTB lineage, and their interaction, within two independent cohorts of TB cases, N = 113 and 121. No association was found between lineage and severity, but association between multiple polymorphisms in IL12B and TBscore was replicated in two independent cohorts (most significant rs3212227, combined p = 0.0006), supporting previous associations of IL12B with TB susceptibility. We also observed significant interaction between a single nucleotide polymorphism (SNP) in SLC11A1 and the L4-Ugandan lineage in both cohorts (rs17235409, meta p = 0.0002). Interestingly, the presence of the L4-Uganda lineage in the presence of the ancestral human allele associated with more severe disease. These findings demonstrate that IL12B is associated with severity of TB in addition to susceptibility, and that the association between TB severity and human genetics can be due to an interaction between genes in the two species, consistent with host-pathogen coevolution in TB.

Myocilin Mutations in Patients With Normal-Tension Glaucoma

Importance

Mutations in the myocilin (MYOC) gene are the most common molecularly defined cause of primary open-angle glaucoma that typically occurs in patients with high intraocular pressures (IOP). One MYOC mutation, p.Gln368Ter, has been associated with as many as 1.6% of primary open-angle glaucoma cases that had a mean maximum recorded IOP of 30 mm Hg. However, to our knowledge, the role of the p.Gln368Ter mutation in patients with normal-tension glaucoma (NTG) with an IOP of 21 mm Hg or lower has not been investigated.

Objective

To evaluate the role of the p.Gln368Ter MYOC mutation in patients with NTG.

Design, Setting, and Participants

In this case-control study of the prevalence of the p.Gln368Ter mutation in patients with NTG, cohort 1 was composed of 772 patients with NTG and 2152 controls from the United States (Iowa, Minnesota, and New York) and England and cohort 2 was composed of 561 patients with NTG and 2606 controls from the Massachusetts Eye and Ear Infirmary and the NEIGHBORHOOD consortium. Genotyping was conducted using real-time polymerase chain reaction that was confirmed with Sanger sequencing, the imputation of genome-wide association study data, or an analysis of whole-exome sequence data. Data analysis occurred between April 2007 and April 2018.

Main Outcomes and Measures

Comparison of the frequency of the p.Gln368Ter MYOC mutation between NTG cases and controls with the Fisher exact test.

Results

Of 6091 total participants, 3346 (54.9%) were women and 5799 (95.2%) were white. We detected the p.Gln368Ter mutation in 7 of 772 patients with NTG (0.91%) and 7 of 2152 controls (0.33%) in cohort 1 (P = .03). In cohort 2, we detected the p.Gln368Ter mutation in 4 of 561 patients with NTG (0.71%) and 10 of 2606 controls (0.38%; P = .15). When the cohorts were analyzed as a group, the p.Gln368Ter mutation was associated with NTG (odds ratio, 2.3; 95% CI, 0.98-5.3; P = .04).

Conclusions and Relevance

In cohorts 1 and 2, the p.Gln368Ter mutation in MYOC was found in patients with IOPs that were 21 mm Hg or lower (NTG), although at a frequency that is lower than previously detected in patients with higher IOP. These data suggest that the p.Gln368Ter mutation may be associated with glaucoma in patients with normal IOPs as well as in patients with IOPs that are greater than 21 mm Hg.

HSD3B1 genotype identifies glucocorticoid responsiveness in severe asthma

Asthma resistance to glucocorticoid treatment is a major health problem with unclear etiology. Glucocorticoids inhibit adrenal androgen production. However, androgens have potential benefits in asthma. HSD3B1 encodes for 3β-hydroxysteroid dehydrogenase-1 (3β-HSD1), which catalyzes peripheral conversion from adrenal dehydroepiandrosterone (DHEA) to potent androgens and has a germline missense-encoding polymorphism. The adrenal restrictive HSD3B1(1245A) allele limits conversion, whereas the adrenal permissive HSD3B1(1245C) allele increases DHEA metabolism to potent androgens. In the Severe Asthma Research Program (SARP) III cohort, we determined the association between DHEA-sulfate and percentage predicted forced expiratory volume in 1 s (FEV1PP). HSD3B1(1245) genotypes were assessed, and association between adrenal restrictive and adrenal permissive alleles and FEV1PP in patients with (GC) and without (noGC) daily oral glucocorticoid treatment was determined (n = 318). Validation was performed in a second cohort (SARP I&II; n = 184). DHEA-sulfate is associated with FEV1PP and is suppressed with GC treatment. GC patients homozygous for the adrenal restrictive genotype have lower FEV1PP compared with noGC patients (54.3% vs. 75.1%; P < 0.001). In patients with the homozygous adrenal permissive genotype, there was no FEV1PP difference in GC vs. noGC patients (73.4% vs. 78.9%; P = 0.39). Results were independently confirmed: FEV1PP for homozygous adrenal restrictive genotype in GC vs. noGC is 49.8 vs. 63.4 (P < 0.001), and for homozygous adrenal permissive genotype, it is 66.7 vs. 67.7 (P = 0.92). The adrenal restrictive HSD3B1(1245) genotype is associated with GC resistance. This effect appears to be driven by GC suppression of 3β-HSD1 substrate. Our results suggest opportunities for prediction of GC resistance and pharmacologic intervention.

Genetic Risk Scores

Genome-wide variation data with millions of genetic markers have become commonplace. However, the potential for interpretation and application of these data for clinical assessment of outcomes of interest, and prediction of disease risk, is currently not fully realized. Many common complex diseases now have numerous, well-established risk loci and likely harbor many genetic determinants with effects too small to be detected at genome-wide levels of statistical significance. A simple and intuitive approach for converting genetic data to a predictive measure of disease susceptibility is to aggregate the effects of these loci into a single measure, the genetic risk score. Here, we describe some common methods and software packages for calculating genetic risk scores and polygenic risk scores, with focus on studies of common complex diseases. We review the basic information needed, as well as important considerations for constructing genetic risk scores, including specific requirements for phenotypic and genetic data, and limitations in their application. © 2019 by John Wiley & Sons, Inc.

Differential Long-Term Outcomes for Individuals With Histories of Preschool Speech Sound Disorders

Purpose The goal of this study was to determine whether adolescent outcomes for individuals with histories of early speech sound disorders (SSD) could be differentiated by speech and language skills at earlier ages (preschool, 4-6 years, and school age, 7-10 years). Method The study used a retrospective longitudinal design. Participants with and without histories of early SSD were classified in adolescence as having no SSD, resolved SSD, low multisyllabic word (MSW; difficulty with MSW repetition but no errors in conversational speech), or persistent speech disorders (errors in both conversational speech and MSW repetition). Analysis of variance was employed to determine whether early speech, language, and literacy skills distinguished these adolescent outcome groups. Results Preschool and school-age skills differed for adolescents whose SSD had resolved from those who had persistent speech errors. Adolescents with errors solely in production of MSWs (Low MSW) did not differ in early speech and language skills from adolescents who had difficulty with both MSWs and persistent errors in conversation. Conclusions Speech and language assessments earlier in childhood can help establish risks for persistent SSD and other language and literacy difficulties in adolescence. Early identification of these clinically relevant subgroups of SSD may allow for early targeted interventions. Supplemental Material https://doi.org/10.23641/asha.9932279.

Association of Genetic Variants With Primary Open-Angle Glaucoma Among Individuals With African Ancestry

IMPORTANCE

Primary open-angle glaucoma presents with increased prevalence and a higher degree of clinical severity in populations of African ancestry compared with European or Asian ancestry. Despite this, individuals of African ancestry remain understudied in genomic research for blinding disorders.

OBJECTIVES

To perform a genome-wide association study (GWAS) of African ancestry populations and evaluate potential mechanisms of pathogenesis for loci associated with primary open-angle glaucoma.

DESIGN, SETTINGS, AND PARTICIPANTS

A 2-stage GWAS with a discovery data set of 2320 individuals with primary open-angle glaucoma and 2121 control individuals without primary open-angle glaucoma. The validation stage included an additional 6937 affected individuals and 14 917 unaffected individuals using multicenter clinic- and population-based participant recruitment approaches. Study participants were recruited from Ghana, Nigeria, South Africa, the United States, Tanzania, Britain, Cameroon, Saudi Arabia, Brazil, the Democratic Republic of the Congo, Morocco, Peru, and Mali from 2003 to 2018. Individuals with primary open-angle glaucoma had open iridocorneal angles and displayed glaucomatous optic neuropathy with visual field defects. Elevated intraocular pressure was not included in the case definition. Control individuals had no elevated intraocular pressure and no signs of glaucoma.

EXPOSURES

Genetic variants associated with primary open-angle glaucoma.

MAIN OUTCOMES AND MEASURES

Presence of primary open-angle glaucoma. Genome-wide significance was defined as P < 5 × 10-8 in the discovery stage and in the meta-analysis of combined discovery and validation data.

RESULTS

A total of 2320 individuals with primary open-angle glaucoma (mean [interquartile range] age, 64.6 [56-74] years; 1055 [45.5%] women) and 2121 individuals without primary open-angle glaucoma (mean [interquartile range] age, 63.4 [55-71] years; 1025 [48.3%] women) were included in the discovery GWAS. The GWAS discovery meta-analysis demonstrated association of variants at amyloid-β A4 precursor protein-binding family B member 2 (APBB2; chromosome 4, rs59892895T>C) with primary open-angle glaucoma (odds ratio [OR], 1.32 [95% CI, 1.20-1.46]; P = 2 × 10-8). The association was validated in an analysis of an additional 6937 affected individuals and 14 917 unaffected individuals (OR, 1.15 [95% CI, 1.09-1.21]; P < .001). Each copy of the rs59892895*C risk allele was associated with increased risk of primary open-angle glaucoma when all data were included in a meta-analysis (OR, 1.19 [95% CI, 1.14-1.25]; P = 4 × 10-13). The rs59892895*C risk allele was present at appreciable frequency only in African ancestry populations. In contrast, the rs59892895*C risk allele had a frequency of less than 0.1% in individuals of European or Asian ancestry.

CONCLUSIONS AND RELEVANCE

In this genome-wide association study, variants at the APBB2 locus demonstrated differential association with primary open-angle glaucoma by ancestry. If validated in additional populations this finding may have implications for risk assessment and therapeutic strategies.

Genetic Correlations Between Diabetes and Glaucoma: An Analysis of Continuous and Dichotomous Phenotypes

PURPOSE

A genetic correlation is the proportion of phenotypic variance between traits that is shared on a genetic basis. Here we explore genetic correlations between diabetes- and glaucoma-related traits.

DESIGN

Cross-sectional study.

METHODS

We assembled genome-wide association study summary statistics from European-derived participants regarding diabetes-related traits like fasting blood sugar (FBS) and type 2 diabetes (T2D) and glaucoma-related traits (intraocular pressure [IOP], central corneal thickness [CCT], corneal hysteresis [CH], corneal resistance factor [CRF], cup-to-disc ratio [CDR], and primary open-angle glaucoma [POAG]). We included data from the National Eye Institute Glaucoma Human Genetics Collaboration Heritable Overall Operational Database, the UK Biobank, and the International Glaucoma Genetics Consortium. We calculated genetic correlation (r_g) between traits using linkage disequilibrium score regression. We also calculated genetic correlations between IOP, CCT, and select diabetes-related traits based on individual level phenotype data in 2 Northern European population-based samples using pedigree information and Sequential Oligogenic Linkage Analysis Routines.

RESULTS

Overall, there was little r_g between diabetes- and glaucoma-related traits. Specifically, we found a nonsignificant negative correlation between T2D and POAG (r_g = -0.14; P = .16). Using Sequential Oligogenic Linkage Analysis Routines, the genetic correlations between measured IOP, CCT, FBS, fasting insulin, and hemoglobin A1c were null. In contrast, genetic correlations between IOP and POAG (r_g ≥ 0.45; P ≤ 3.0 × 10^-4) and between CDR and POAG were high (r_g = 0.57; P = 2.8 × 10^-10). However, genetic correlations between corneal properties (CCT, CRF, and CH) and POAG were low (r_g range -0.18 to 0.11) and nonsignificant (P ≥ .07).

CONCLUSION

These analyses suggest that there is limited genetic correlation between diabetes- and glaucoma-related traits.

Pathway Analysis Integrating Genome-Wide and Functional Data Identifies PLCG2 as a Candidate Gene for Age-Related Macular Degeneration

Purpose

Age-related macular degeneration (AMD) is the worldwide leading cause of blindness among the elderly. Although genome-wide association studies (GWAS) have identified AMD risk variants, their roles in disease etiology are not well-characterized, and they only explain a portion of AMD heritability.

Methods

We performed pathway analyses using summary statistics from the International AMD Genomics Consortium's 2016 GWAS and multiple pathway databases to identify biological pathways wherein genetic association signals for AMD may be aggregating. We determined which genes contributed most to significant pathway signals across the databases. We characterized these genes by constructing protein-protein interaction networks and performing motif analysis.

Results

We determined that eight genes (C2, C3, LIPC, MICA, NOTCH4, PLCG2, PPARA, and RAD51B) "drive" the statistical signals observed across pathways curated in the Kyoto Encyclopedia of Genes and Genomes (KEGG), Reactome, and Gene Ontology (GO) databases. We further refined our definition of statistical driver gene to identify PLCG2 as a candidate gene for AMD due to its significant gene-level signals (P < 0.0001) across KEGG, Reactome, GO, and NetPath pathways.

Conclusions

We performed pathway analyses on the largest available collection of advanced AMD cases and controls in the world. Eight genes strongly contributed to significant pathways from the three larger databases, and one gene (PLCG2) was central to significant pathways from all four databases. This is, to our knowledge, the first study to identify PLCG2 as a candidate gene for AMD based solely on genetic burden. Our findings reinforce the utility of integrating in silico genetic and biological pathway data to investigate the genetic architecture of AMD.

Genetically-guided algorithm development and sample size optimization for age-related macular degeneration cases and controls in electronic health records from the VA Million Veteran Program

Electronic health records (EHRs) linked to extensive biorepositories and supplemented with lifestyle, behavioral, and environmental exposure data, have enormous potential to contribute to genomic discovery, a necessary step in the pathway towards translational or precision medicine. A major bottleneck in incorporating EHRs into genomic studies is the extraction of research-grade variables for analysis, particularly when gold-standard measurements are not available or accessible. Here we develop algorithms for age-related macular degeneration (AMD), a common cause of blindness among the elderly, and controls free of AMD. These computable phenotypes were developed using billing codes (ICD-9-CM and ICD-10-CM) and Current Procedural Terminology (CPT) codes and evaluated in two study sites of the Veterans Affairs Million Veteran Program: Louis Stokes Cleveland VA Medical Center and the Providence VA Medical Center. After establishing a high overall positive and negative predictive values (93% and 95%, respectively) through manual chart review, the candidate algorithm was deployed in the full VA MVP dataset of >500,000 participants. The algorithm was then optimized in a data cube using a variety of approaches including adjusting inclusion age thresholds by examining previously-reported genetic associations for CFH (rs10801555, a proxy for rs1061170) and ARMS2 (rs10490924). The algorithm with the smallest p-values for the known genetic associations was selected for downstream and on-going AMD genomic discovery efforts. This two-phase approach to developing research-grade case/control variables for AMD genomic studies capitalizes on established genetic associations resulting in high precision and optimized sample sizes, an approach that can be applied to other large-scale biobanks linked to EHRs for precision medicine research.

Multiethnic Genome-Wide Association Study of Diabetic Retinopathy Using Liability Threshold Modeling of Duration of Diabetes and Glycemic Control

To identify genetic variants associated with diabetic retinopathy (DR), we performed a large multiethnic genome-wide association study. Discovery included eight European cohorts (n = 3,246) and seven African American cohorts (n = 2,611). We meta-analyzed across cohorts using inverse-variance weighting, with and without liability threshold modeling of glycemic control and duration of diabetes. Variants with a P value <1 × 10-5 were investigated in replication cohorts that included 18,545 European, 16,453 Asian, and 2,710 Hispanic subjects. After correction for multiple testing, the C allele of rs142293996 in an intron of nuclear VCP-like (NVL) was associated with DR in European discovery cohorts (P = 2.1 × 10-9), but did not reach genome-wide significance after meta-analysis with replication cohorts. We applied the Disease Association Protein-Protein Link Evaluator (DAPPLE) to our discovery results to test for evidence of risk being spread across underlying molecular pathways. One protein-protein interaction network built from genes in regions associated with proliferative DR was found to have significant connectivity (P = 0.0009) and corroborated with gene set enrichment analyses. These findings suggest that genetic variation in NVL, as well as variation within a protein-protein interaction network that includes genes implicated in inflammation, may influence risk for DR.

Testosterone Pathway Genetic Polymorphisms in Relation to Primary Open-Angle Glaucoma: An Analysis in Two Large Datasets

Purpose

Sex hormones may be associated with primary open-angle glaucoma (POAG), although the mechanisms are unclear. We previously observed that gene variants involved with estrogen metabolism were collectively associated with POAG in women but not men; here we assessed gene variants related to testosterone metabolism collectively and POAG risk.

Methods

We used two datasets: one from the United States (3853 cases and 33,480 controls) and another from Australia (1155 cases and 1992 controls). Both datasets contained densely called genotypes imputed to the 1000 Genomes reference panel. We used pathway- and gene-based approaches with Pathway Analysis by Randomization Incorporating Structure (PARIS) software to assess the overall association between a panel of single nucleotide polymorphisms (SNPs) in testosterone metabolism genes and POAG. In sex-stratified analyses, we evaluated POAG overall and POAG subtypes defined by maximum IOP (high-tension [HTG] or normal tension glaucoma [NTG]).

Results

In the US dataset, the SNP panel was not associated with POAG (permuted P = 0.77), although there was an association in the Australian sample (permuted P = 0.018). In both datasets, the SNP panel was associated with POAG in men (permuted P ≤ 0.033) and not women (permuted P ≥ 0.42), but in gene-based analyses, there was no consistency on the main genes responsible for these findings. In both datasets, the testosterone pathway association with HTG was significant (permuted P ≤ 0.011), but again, gene-based analyses showed no consistent driver gene associations.

Conclusions

Collectively, testosterone metabolism pathway SNPs were consistently associated with the high-tension subtype of POAG in two datasets.

A Genome-Wide Association Study of Diabetic Kidney Disease in Subjects With Type 2 Diabetes

Identification of sequence variants robustly associated with predisposition to diabetic kidney disease (DKD) has the potential to provide insights into the pathophysiological mechanisms responsible. We conducted a genome-wide association study (GWAS) of DKD in type 2 diabetes (T2D) using eight complementary dichotomous and quantitative DKD phenotypes: the principal dichotomous analysis involved 5,717 T2D subjects, 3,345 with DKD. Promising association signals were evaluated in up to 26,827 subjects with T2D (12,710 with DKD). A combined T1D+T2D GWAS was performed using complementary data available for subjects with T1D, which, with replication samples, involved up to 40,340 subjects with diabetes (18,582 with DKD). Analysis of specific DKD phenotypes identified a novel signal near GABRR1 (rs9942471, P = 4.5 × 10^-8) associated with microalbuminuria in European T2D case subjects. However, no replication of this signal was observed in Asian subjects with T2D or in the equivalent T1D analysis. There was only limited support, in this substantially enlarged analysis, for association at previously reported DKD signals, except for those at UMOD and PRKAG2, both associated with estimated glomerular filtration rate. We conclude that, despite challenges in addressing phenotypic heterogeneity, access to increased sample sizes will continue to provide more robust inference regarding risk variant discovery for DKD.

Genome-wide association meta-analysis highlights light-induced signaling as a driver for refractive error

Refractive errors, including myopia, are the most frequent eye disorders worldwide and an increasingly common cause of blindness. This genome-wide association meta-analysis in 160,420 participants and replication in 95,505 participants increased the number of established independent signals from 37 to 161 and showed high genetic correlation between Europeans and Asians (>0.78). Expression experiments and comprehensive in silico analyses identified retinal cell physiology and light processing as prominent mechanisms, and also identified functional contributions to refractive-error development in all cell types of the neurosensory retina, retinal pigment epithelium, vascular endothelium and extracellular matrix. Newly identified genes implicate novel mechanisms such as rod-and-cone bipolar synaptic neurotransmission, anterior-segment morphology and angiogenesis. Thirty-one loci resided in or near regions transcribing small RNAs, thus suggesting a role for post-transcriptional regulation. Our results support the notion that refractive errors are caused by a light-dependent retina-to-sclera signaling cascade and delineate potential pathobiological molecular drivers.

Genome-wide analyses identify 68 new loci associated with intraocular pressure and improve risk prediction for primary open-angle glaucoma

Glaucoma is the leading cause of irreversible blindness globally 1 . Despite its gravity, the disease is frequently undiagnosed in the community 2 . Raised intraocular pressure (IOP) is the most important risk factor for primary open-angle glaucoma (POAG)3,4. Here we present a meta-analysis of 139,555 European participants, which identified 112 genomic loci associated with IOP, 68 of which are novel. These loci suggest a strong role for angiopoietin-receptor tyrosine kinase signaling, lipid metabolism, mitochondrial function and developmental processes underlying risk for elevated IOP. In addition, 48 of these loci were nominally associated with glaucoma in an independent cohort, 14 of which were significant at a Bonferroni-corrected threshold. Regression-based glaucoma-prediction models had an area under the receiver operating characteristic curve (AUROC) of 0.76 in US NEIGHBORHOOD study participants and 0.74 in independent glaucoma cases from the UK Biobank. Genetic-prediction models for POAG offer an opportunity to target screening and timely therapy to individuals most at risk.

AmpliSeq transcriptome analysis of human alveolar and monocyte-derived macrophages over time in response to Mycobacterium tuberculosis infection

Human alveolar macrophages (HAM) are primary bacterial niche and immune response cells during Mycobacterium tuberculosis (M.tb) infection, and human blood monocyte-derived macrophages (MDM) are a model for investigating M.tb-macrophage interactions. Here, we use a targeted RNA-Seq method to measure transcriptome-wide changes in RNA expression patterns of freshly obtained HAM (used within 6 h) and 6 day cultured MDM upon M.tb infection over time (2, 24 and 72 h), in both uninfected and infected cells from three donors each. The Ion AmpliSeq™ Transcriptome Human Gene Expression Kit (AmpliSeq) uses primers targeting 18,574 mRNAs and 2,228 non-coding RNAs (ncRNAs) for a total of 20,802 transcripts. AmpliSeqTM yields highly precise and reproducible gene expression profiles (R2 >0.99). Taking advantage of AmpliSeq's reproducibility, we establish well-defined quantitative RNA expression patterns of HAM versus MDM, including significant M.tb-inducible genes, in networks and pathways that differ in part between MDM and HAM. A similar number of expressed genes are detected at all time-points between uninfected MDM and HAM, in common pathways including inflammatory and immune functions, but canonical pathway differences also exist. In particular, at 2 h, multiple genes relevant to the immune response are preferentially expressed in either uninfected HAM or MDM, while the HAM RNA profiles approximate MDM profiles over time in culture, highlighting the unique RNA expression profile of freshly obtained HAM. MDM demonstrate a greater transcriptional response than HAM upon M.tb infection, with 2 to >10 times more genes up- or down-regulated. The results identify key genes involved in cellular responses to M.tb in two different human macrophage types. Follow-up bioinformatics analysis indicates that approximately 30% of response genes have expression quantitative trait loci (eQTLs in GTEx), common DNA variants that can influence host gene expression susceptibility or resistance to M.tb, illustrated with the TREM1 gene cluster and IL-10.

Genetic correlations between intraocular pressure, blood pressure and primary open-angle glaucoma: a multi-cohort analysis

Primary open-angle glaucoma (POAG) is the most common chronic optic neuropathy worldwide. Epidemiological studies show a robust positive relation between intraocular pressure (IOP) and POAG and modest positive association between IOP and blood pressure (BP), while the relation between BP and POAG is controversial. The International Glaucoma Genetics Consortium (n=27 558), the International Consortium on Blood Pressure (n=69 395), and the National Eye Institute Glaucoma Human Genetics Collaboration Heritable Overall Operational Database (n=37 333), represent genome-wide data sets for IOP, BP traits and POAG, respectively. We formed genome-wide significant variant panels for IOP and diastolic BP and found a strong relation with POAG (odds ratio and 95% confidence interval: 1.18 (1.14-1.21), P=1.8 × 10^-27) for the former trait but no association for the latter (P=0.93). Next, we used linkage disequilibrium (LD) score regression, to provide genome-wide estimates of correlation between traits without the need for additional phenotyping. We also compared our genome-wide estimate of heritability between IOP and BP to an estimate based solely on direct measures of these traits in the Erasmus Rucphen Family (ERF; n=2519) study using Sequential Oligogenic Linkage Analysis Routines (SOLAR). LD score regression revealed high genetic correlation between IOP and POAG (48.5%, P=2.1 × 10^-5); however, genetic correlation between IOP and diastolic BP (P=0.86) and between diastolic BP and POAG (P=0.42) were negligible. Using SOLAR in the ERF study, we confirmed the minimal heritability between IOP and diastolic BP (P=0.63). Overall, IOP shares genetic basis with POAG, whereas BP has limited shared genetic correlation with IOP or POAG.

A chromosome 5q31.1 locus associates with tuberculin skin test reactivity in HIV-positive individuals from tuberculosis hyper-endemic regions in east Africa

One in three people has been infected with Mycobacterium tuberculosis (MTB), and the risk for MTB infection in HIV-infected individuals is even higher. We hypothesized that HIV-positive individuals living in tuberculosis-endemic regions who do not get infected by Mycobacterium tuberculosis are genetically resistant. Using an “experiment of nature” design that proved successful in our previous work, we performed a genome-wide association study of tuberculin skin test positivity using 469 HIV-positive patients from prospective study cohorts of tuberculosis from Tanzania and Uganda to identify genetic loci associated with MTB infection in the context of HIV-infection. Among these individuals, 244 tested were tuberculin skin test (TST) positive either at enrollment or during the >8 year follow up, while 225 were not. We identified a genome-wide significant association between a dominant model of rs877356 and binary TST status in the combined cohort (Odds ratio = 0.2671, p = 1.22x10^-8). Association was replicated with similar significance when examining TST induration as a continuous trait. The variant lies in the 5q31.1 region, 57kb downstream from IL9. Two-locus analyses of association of variants near rs877356 showed a haplotype comprised of rs877356 and an IL9 missense variant, rs2069885, had the most significant association (p = 1.59x10^-12). We also replicated previously linked loci on chromosomes 2, 5, and 11. IL9 is a cytokine produced by mast cells and T_H2 cells during inflammatory responses, providing a possible link between airway inflammation and protection from MTB infection. Our results indicate that studying uninfected, HIV-positive participants with extensive exposure increases the power to detect associations in complex infectious disease.

Approximately one-third of the world’s population has been exposed to Mycobacterium tuberculosis, the bacterium that causes tuberculosis. A small number of those infected develop active disease; however, there is a substantial portion of exposed people who do not even show evidence of an immunological response. These people who appear to resist infection, as measured by a negative tuberculin skin test, represent a subpopulation from which we can learn about resistance. We used a genome-wide approach to study the genetic basis of this resistance in unique cohorts of hypervulnerable, HIV-positive individuals from Uganda and Tanzania, in which exposure was virtually assured. We identified one locus that was highly significantly associated and conferred more than 70% protection from infection. The most significant variant, rs8773656, was near IL9 and SLC25A48, and a haplotype including this variant and a missense mutation in IL9 was even more significantly associated with negative skin tests. Although it is impossible based solely on our data to determine the causal variant or genes, IL9 is an attractive candidate as its product has previously been associated with bronchial hyperresponsiveness, thereby providing a possible link between inflammation and protection from Mycobacterium tuberculosis infection.

Genetic association study of exfoliation syndrome identifies a protective rare variant at LOXL1 and five new susceptibility loci

Exfoliation syndrome (XFS) is the most common known risk factor for secondary glaucoma and a major cause of blindness worldwide. Variants in two genes, LOXL1 and CACNA1A, have previously been associated with XFS. To further elucidate the genetic basis of XFS, we collected a global sample of XFS cases to refine the association at LOXL1, which previously showed inconsistent results across populations, and to identify new variants associated with XFS. We identified a rare protective allele at LOXL1 (p.Phe407, odds ratio (OR) = 25, P = 2.9 × 10^-14) through deep resequencing of XFS cases and controls from nine countries. A genome-wide association study (GWAS) of XFS cases and controls from 24 countries followed by replication in 18 countries identified seven genome-wide significant loci (P < 5 × 10^-8). We identified association signals at 13q12 (POMP), 11q23.3 (TMEM136), 6p21 (AGPAT1), 3p24 (RBMS3) and 5q23 (near SEMA6A). These findings provide biological insights into the pathology of XFS and highlight a potential role for naturally occurring rare LOXL1 variants in disease biology.

A Common Variant in MIR182 Is Associated With Primary Open-Angle Glaucoma in the NEIGHBORHOOD Consortium

Purpose

Noncoding microRNAs (miRNAs) have been implicated in the pathogenesis of glaucoma. We aimed to identify common variants in miRNA coding genes (MIR) associated with primary open-angle glaucoma (POAG).

Methods

Using the NEIGHBORHOOD data set (3853 cases/33,480 controls with European ancestry), we first assessed the relation between 85 variants in 76 MIR genes and overall POAG. Subtype-specific analyses were performed in high-tension glaucoma (HTG) and normal-tension glaucoma subsets. Second, we examined the expression of miR-182, which was associated with POAG, in postmortem human ocular tissues (ciliary body, cornea, retina, and trabecular meshwork [TM]), using miRNA sequencing (miRNA-Seq) and droplet digital PCR (ddPCR). Third, miR-182 expression was also examined in human aqueous humor (AH) by using miRNA-Seq. Fourth, exosomes secreted from primary human TM cells were examined for miR-182 expression by using miRNA-Seq. Fifth, using ddPCR we compared miR-182 expression in AH between five HTG cases and five controls.

Results

Only rs76481776 in MIR182 gene was associated with POAG after adjustment for multiple comparisons (odds ratio [OR] = 1.23, 95% confidence interval [CI]: 1.11–1.42, P = 0.0002). Subtype analysis indicated that the association was primarily in the HTG subset (OR = 1.26, 95% CI: 1.08–1.47, P = 0.004). The risk allele T has been associated with elevated miR-182 expression in vitro. Data from ddPCR and miRNA-Seq confirmed miR-182 expression in all examined ocular tissues and TM-derived exosomes. Interestingly, miR-182 expression in AH was 2-fold higher in HTG patients than nonglaucoma controls (P = 0.03) without controlling for medication treatment.

Conclusions

Our integrative study is the first to associate rs76481776 with POAG via elevated miR-182 expression.

Genome-wide association study identifies three novel loci in Fuchs endothelial corneal dystrophy

The structure of the cornea is vital to its transparency, and dystrophies that disrupt corneal organization are highly heritable. To understand the genetic aetiology of Fuchs endothelial corneal dystrophy (FECD), the most prevalent corneal disorder requiring transplantation, we conducted a genome-wide association study (GWAS) on 1,404 FECD cases and 2,564 controls of European ancestry, followed by replication and meta-analysis, for a total of 2,075 cases and 3,342 controls. We identify three novel loci meeting genome-wide significance (P<5 × 10-8): KANK4 rs79742895, LAMC1 rs3768617 and LINC00970/ATP1B1 rs1200114. We also observe an overwhelming effect of the established TCF4 locus. Interestingly, we detect differential sex-specific association at LAMC1, with greater risk in women, and TCF4, with greater risk in men. Combining GWAS results with biological evidence we expand the knowledge of common FECD loci from one to four, and provide a deeper understanding of the underlying pathogenic basis of FECD.

Age at natural menopause genetic risk score in relation to age at natural menopause and primary open-angle glaucoma in a US-based sample

OBJECTIVE

Several attributes of female reproductive history, including age at natural menopause (ANM), have been related to primary open-angle glaucoma (POAG). We assembled 18 previously reported common genetic variants that predict ANM to determine their association with ANM or POAG.

METHODS

Using data from the Nurses' Health Study (7,143 women), we validated the ANM weighted genetic risk score in relation to self-reported ANM. Subsequently, to assess the relation with POAG, we used data from 2,160 female POAG cases and 29,110 controls in the National Eye Institute Glaucoma Human Genetics Collaboration Heritable Overall Operational Database (NEIGHBORHOOD), which consists of 8 datasets with imputed genotypes to 5.6+ million markers. Associations with POAG were assessed in each dataset, and site-specific results were meta-analyzed using the inverse weighted variance method.

RESULTS

The genetic risk score was associated with self-reported ANM (P = 2.2 × 10) and predicted 4.8% of the variance in ANM. The ANM genetic risk score was not associated with POAG (Odds Ratio (OR) = 1.002; 95% Confidence Interval (CI): 0.998, 1.007; P = 0.28). No single genetic variant in the panel achieved nominal association with POAG (P ≥0.20). Compared to the middle 80 percent, there was also no association with the lowest 10 percentile or highest 90 percentile of genetic risk score with POAG (OR = 0.75; 95% CI: 0.47, 1.21; P = 0.23 and OR = 1.10; 95% CI: 0.72, 1.69; P = 0.65, respectively).

CONCLUSIONS

A genetic risk score predicting 4.8% of ANM variation was not related to POAG; thus, genetic determinants of ANM are unlikely to explain the previously reported association between the two phenotypes.

Genetic Risk Scores

The generation of genome-wide variation data has become commonplace. However, the potential for interpretation and application of these data for clinical assessment of outcomes of interest, and prediction of disease risk, is currently not fully realized. Many common, complex diseases now have numerous, well-established "risk" loci, and likely harbor many genetic determinants with effects too small to be detected at genome-wide levels of statistical significance. A simple and intuitive approach for converting genetic data to a predictive measure of disease susceptibility is to aggregate the risk effects of these loci into a single genetic risk score. Here, some common methods and software packages for calculating genetic risk scores, with focus on studies of common, complex diseases, are described. The basic information needed as well as important considerations for constructing genetic risk scores, including specific requirements for phenotypic and genetic data, and limitations in their application is reviewed. © 2016 by John Wiley & Sons, Inc.

Quality Control for the Illumina HumanExome BeadChip

The Illumina HumanExome BeadChip and other exome-based genotyping arrays offer inexpensive genotyping of some 240,000 mostly nonsynonymous coding variants across the human genome. The HumanExome chip, with its highly non-uniform distribution of markers and emphasis on rare coding variants, presents some unique challenges for quality control (QC) and data cleaning. Here, we describe QC procedures for HumanExome data, with examples of challenges specific to exome arrays from our experience cleaning a data set of ∼7,500 samples from the NEIGHBORHOOD Consortium. We focus on standard procedures for QC of genome-wide array data including genotype calling, sex verification, sample identity verification, relationship checking, and population structure that are complicated by the HumanExome panel's enrichment in rare, exonic variation. © 2016 by John Wiley & Sons, Inc.

Selecting SNPs informative for African, American Indian and European Ancestry: application to the Family Investigation of Nephropathy and Diabetes (FIND)

Background

The presence of population structure in a sample may confound the search for important genetic loci associated with disease. Our four samples in the Family Investigation of Nephropathy and Diabetes (FIND), European Americans, Mexican Americans, African Americans, and American Indians are part of a genome- wide association study in which population structure might be particularly important. We therefore decided to study in detail one component of this, individual genetic ancestry (IGA). From SNPs present on the Affymetrix 6.0 Human SNP array, we identified 3 sets of ancestry informative markers (AIMs), each maximized for the information in one the three contrasts among ancestral populations: Europeans (HAPMAP, CEU), Africans (HAPMAP, YRI and LWK), and Native Americans (full heritage Pima Indians). We estimate IGA and present an algorithm for their standard errors, compare IGA to principal components, emphasize the importance of balancing information in the ancestry informative markers (AIMs), and test the association of IGA with diabetic nephropathy in the combined sample.

Results

A fixed parental allele maximum likelihood algorithm was applied to the FIND to estimate IGA in four samples: 869 American Indians; 1385 African Americans; 1451 Mexican Americans; and 826 European Americans. When the information in the AIMs is unbalanced, the estimates are incorrect with large error. Individual genetic admixture is highly correlated with principle components for capturing population structure. It takes ~700 SNPs to reduce the average standard error of individual admixture below 0.01. When the samples are combined, the resulting population structure creates associations between IGA and diabetic nephropathy.

Conclusions

The identified set of AIMs, which include American Indian parental allele frequencies, may be particularly useful for estimating genetic admixture in populations from the Americas. Failure to balance information in maximum likelihood, poly-ancestry models creates biased estimates of individual admixture with large error. This also occurs when estimating IGA using the Bayesian clustering method as implemented in the program STRUCTURE. Odds ratios for the associations of IGA with disease are consistent with what is known about the incidence and prevalence of diabetic nephropathy in these populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2654-x) contains supplementary material, which is available to authorized users.

Meta-analysis of gene-environment-wide association scans accounting for education level identifies additional loci for refractive error

Myopia is the most common human eye disorder and it results from complex genetic and environmental causes. The rapidly increasing prevalence of myopia poses a major public health challenge. Here, the CREAM consortium performs a joint meta-analysis to test single-nucleotide polymorphism (SNP) main effects and SNP × education interaction effects on refractive error in 40,036 adults from 25 studies of European ancestry and 10,315 adults from 9 studies of Asian ancestry. In European ancestry individuals, we identify six novel loci (FAM150B-ACP1, LINC00340, FBN1, DIS3L-MAP2K1, ARID2-SNAT1 and SLC14A2) associated with refractive error. In Asian populations, three genome-wide significant loci AREG, GABRR1 and PDE10A also exhibit strong interactions with education (P<8.5 × 10(-5)), whereas the interactions are less evident in Europeans. The discovery of these loci represents an important advance in understanding how gene and environment interactions contribute to the heterogeneity of myopia.

A Locus at 5q33.3 Confers Resistance to Tuberculosis in Highly Susceptible Individuals

Immunosuppression resulting from HIV infection increases the risk of progression to active tuberculosis (TB) both in individuals newly exposed to Mycobacterium tuberculosis (MTB) and in those with latent infections. We hypothesized that HIV-positive individuals who do not develop TB, despite living in areas where it is hyperendemic, provide a model of natural resistance. We performed a genome-wide association study of TB resistance by using 581 HIV-positive Ugandans and Tanzanians enrolled in prospective cohort studies of TB; 267 of these individuals developed active TB, and 314 did not. A common variant, rs4921437 at 5q33.3, was significantly associated with TB (odds ratio = 0.37, p = 2.11 × 10(-8)). This variant lies within a genomic region that includes IL12B and is embedded in an H3K27Ac histone mark. The locus also displays consistent patterns of linkage disequilibrium across African populations and has signals of strong selection in populations from equatorial Africa. Along with prior studies demonstrating that therapy with IL-12 (the cytokine encoded in part by IL12B, associated with longer survival following MTB infection in mice deficient in CD4 T cells), our results suggest that this pathway might be an excellent target for the development of new modalities for treating TB, especially for HIV-positive individuals. Our results also indicate that studying extreme disease resistance in the face of extensive exposure can increase the power to detect associations in complex infectious disease.

Genome-wide association analysis identifies TXNRD2, ATXN2 and FOXC1 as susceptibility loci for primary open-angle glaucoma

Primary open angle glaucoma (POAG) is a leading cause of blindness world-wide. To identify new susceptibility loci, we meta-analyzed GWAS results from 8 independent studies from the United States (3,853 cases and 33,480 controls) and investigated the most significant SNPs in two Australian studies (1,252 cases and 2,592 controls), 3 European studies (875 cases and 4,107 controls) and a Singaporean Chinese study (1,037 cases and 2,543 controls). A meta-analysis of top SNPs identified three novel loci: rs35934224[T] within TXNRD2 (odds ratio (OR) = 0.78, P = 4.05×10⁻¹¹ encoding a mitochondrial protein required for redox homeostasis; rs7137828[T] within ATXN2 (OR = 1.17, P = 8.73×10⁻¹⁰), and rs2745572[A] upstream of FOXC1 (OR = 1.17, P = 1.76×10⁻¹⁰). Using RT-PCR and immunohistochemistry, we show TXNRD2 and ATXN2 expression in retinal ganglion cells and the optic nerve head. These results identify new pathways underlying POAG susceptibility and suggest novel targets for preventative therapies.

Genome-Wide Association and Trans-ethnic Meta-Analysis for Advanced Diabetic Kidney Disease: Family Investigation of Nephropathy and Diabetes (FIND)

Diabetic kidney disease (DKD) is the most common etiology of chronic kidney disease (CKD) in the industrialized world and accounts for much of the excess mortality in patients with diabetes mellitus. Approximately 45% of U.S. patients with incident end-stage kidney disease (ESKD) have DKD. Independent of glycemic control, DKD aggregates in families and has higher incidence rates in African, Mexican, and American Indian ancestral groups relative to European populations. The Family Investigation of Nephropathy and Diabetes (FIND) performed a genome-wide association study (GWAS) contrasting 6,197 unrelated individuals with advanced DKD with healthy and diabetic individuals lacking nephropathy of European American, African American, Mexican American, or American Indian ancestry. A large-scale replication and trans-ethnic meta-analysis included 7,539 additional European American, African American and American Indian DKD cases and non-nephropathy controls. Within ethnic group meta-analysis of discovery GWAS and replication set results identified genome-wide significant evidence for association between DKD and rs12523822 on chromosome 6q25.2 in American Indians (P = 5.74x10^-9). The strongest signal of association in the trans-ethnic meta-analysis was with a SNP in strong linkage disequilibrium with rs12523822 (rs955333; P = 1.31x10^-8), with directionally consistent results across ethnic groups. These 6q25.2 SNPs are located between the SCAF8 and CNKSR3 genes, a region with DKD relevant changes in gene expression and an eQTL with IPCEF1, a gene co-translated with CNKSR3. Several other SNPs demonstrated suggestive evidence of association with DKD, within and across populations. These data identify a novel DKD susceptibility locus with consistent directions of effect across diverse ancestral groups and provide insight into the genetic architecture of DKD.

Type 2 diabetes is the most common cause of severe kidney disease worldwide and diabetic kidney disease (DKD) associates with premature death. Individuals of non-European ancestry have the highest burden of type 2 DKD; hence understanding the causes of DKD remains critical to reducing health disparities. Family studies demonstrate that genes regulate the onset and progression of DKD; however, identifying these genes has proven to be challenging. The Family Investigation of Diabetes and Nephropathy consortium (FIND) recruited a large multi-ethnic collection of individuals with type 2 diabetes with and without kidney disease in order to detect genes associated with DKD. FIND discovered and replicated a DKD-associated genetic locus on human chromosome 6q25.2 (rs955333) between the SCAF8 and CNKSR genes. Findings were supported by significantly different expression of genes in this region from kidney tissue of subjects with, versus without DKD. The present findings identify a novel kidney disease susceptibility locus in individuals with type 2 diabetes which is consistent across subjects of differing ancestries. In addition, FIND results provide a rich catalogue of genetic variation in DKD patients for future research on the genetic architecture regulating this common and devastating disease.

Autosomal dominant hereditary spastic paraplegia with axonal sensory motor polyneuropathy maps to chromosome 21q 22.3

AIM

Hereditary spastic paraplegia (HSP) are a genetically and clinically heterogeneous group of disorders. At present, 19 autosomal dominant loci for HSP have been mapped. We ascertained an American family of European descent segregating an autosomal dominant HSP associated with peripheral neuropathy.

METHODS

A genome wide scan was performed with 410 microsatellite repeat marker (Weber lab screening set 16) and following linkage and haplotype analysis, fine mapping was performed. Established genes or loci for HSP were excluded by direct sequencing or haplotype analysis.

RESULTS

All established loci for HSP were excluded. Fine mapping suggested a locus on chromosome 21q22.3 flanked by markers D21S1411 and D21S1446 with a maximum logarithm of odds score of 2.05 and was supported by haplotype analysis. A number of candidate genes in this region were analyzed and no disease-producing mutations were detected.

CONCLUSION

We present the clinical and genetic analysis of an American family with autosomal dominant HSP with axonal sensory motor polyneuropathy mapping to a novel locus on chromosome 21q22.3 designated SPG56.

Practical barriers and ethical challenges in genetic data sharing

The underlying ethos of dbGaP is that access to these data by secondary data analysts facilitates advancement of science. NIH has required that genome-wide association study data be deposited in the Database of Genotypes and Phenotypes (dbGaP) since 2003. In 2013, a proposed updated policy extended this requirement to next-generation sequencing data. However, recent literature and anecdotal reports suggest lingering logistical and ethical concerns about subject identifiability, informed consent, publication embargo enforcement, and difficulty in accessing dbGaP data. We surveyed the International Genetic Epidemiology Society (IGES) membership about their experiences. One hundred and seventy five (175) individuals completed the survey, a response rate of 27%. Of respondents who received data from dbGaP (43%), only 32% perceived the application process as easy but most (75%) received data within five months. Remaining challenges include difficulty in identifying an institutional signing official and an overlong application process. Only 24% of respondents had contributed data to dbGaP. Of these, 31% reported local IRB restrictions on data release; an additional 15% had to reconsent study participants before depositing data. The majority of respondents (56%) disagreed that the publication embargo period was sufficient. In response, we recommend longer embargo periods and use of varied data-sharing models rather than a one-size-fits-all approach.

Meta-analysis of genome-wide association studies in African Americans provides insights into the genetic architecture of type 2 diabetes

Type 2 diabetes (T2D) is more prevalent in African Americans than in Europeans. However, little is known about the genetic risk in African Americans despite the recent identification of more than 70 T2D loci primarily by genome-wide association studies (GWAS) in individuals of European ancestry. In order to investigate the genetic architecture of T2D in African Americans, the MEta-analysis of type 2 DIabetes in African Americans (MEDIA) Consortium examined 17 GWAS on T2D comprising 8,284 cases and 15,543 controls in African Americans in stage 1 analysis. Single nucleotide polymorphisms (SNPs) association analysis was conducted in each study under the additive model after adjustment for age, sex, study site, and principal components. Meta-analysis of approximately 2.6 million genotyped and imputed SNPs in all studies was conducted using an inverse variance-weighted fixed effect model. Replications were performed to follow up 21 loci in up to 6,061 cases and 5,483 controls in African Americans, and 8,130 cases and 38,987 controls of European ancestry. We identified three known loci (TCF7L2, HMGA2 and KCNQ1) and two novel loci (HLA-B and INS-IGF2) at genome-wide significance (4.15×10⁻⁹⁴<P<5×10⁻⁸, odds ratio (OR) = 1.09 to 1.36). Fine-mapping revealed that 88 of 158 previously identified T2D or glucose homeostasis loci demonstrated nominal to highly significant association (2.2×10⁻²³ < locus-wide P<0.05). These novel and previously identified loci yielded a sibling relative risk of 1.19, explaining 17.5% of the phenotypic variance of T2D on the liability scale in African Americans. Overall, this study identified two novel susceptibility loci for T2D in African Americans. A substantial number of previously reported loci are transferable to African Americans after accounting for linkage disequilibrium, enabling fine mapping of causal variants in trans-ethnic meta-analysis studies.

Despite the higher prevalence of type 2 diabetes (T2D) in African Americans than in Europeans, recent genome-wide association studies (GWAS) were examined primarily in individuals of European ancestry. In this study, we performed meta-analysis of 17 GWAS in 8,284 cases and 15,543 controls to explore the genetic architecture of T2D in African Americans. Following replication in additional 6,061 cases and 5,483 controls in African Americans, and 8,130 cases and 38,987 controls of European ancestry, we identified two novel and three previous reported T2D loci reaching genome-wide significance. We also examined 158 loci previously reported to be associated with T2D or regulating glucose homeostasis. While 56% of these loci were shared between African Americans and the other populations, the strongest associations in African Americans are often found in nearby single nucleotide polymorphisms (SNPs) instead of the original SNPs reported in other populations due to differential genetic architecture across populations. Our results highlight the importance of performing genetic studies in non-European populations to fine map the causal genetic variants.

Clinical and epidemiological characteristics of individuals resistant to M. tuberculosis infection in a longitudinal TB household contact study in Kampala, Uganda

Background

Despite sustained exposure to a person with pulmonary tuberculosis (TB), some M. tuberculosis (Mtb) exposed individuals maintain a negative tuberculin skin test (TST). Our objective was to characterize these persistently negative TST (PTST-) individuals and compare them to TST converters (TSTC) and individuals who are TST positive at study enrollment.

Methods

During a TB household contact study in Kampala, Uganda, PTST-, TSTC, and TST + individuals were identified. PTST- individuals maintained a negative TST over a 2 year observation period despite prolonged exposure to an infectious tuberculosis (TB) case. Epidemiological and clinical characteristics were compared, a risk score developed by another group to capture risk for Mtb infection was computed, and an ordinal regression was performed.

Results

When analyzed independently, epidemiological risk factors increased in prevalence from PTST- to TSTC to TST+. An ordinal regression model suggested age (p < 0.01), number of windows (p < 0.01) and people (p = 0.07) in the home, and sleeping in the same room (p < 0.01) were associated with PTST- and TSTC. As these factors do not exist in isolation, we examined a risk score, which reflects an accumulation of risk factors. This compound exposure score did not differ significantly between PTST-, TSTC, and TST+, except for the 5–15 age group (p = 0.009).

Conclusions

Though many individual factors differed across all three groups, an exposure risk score reflecting a collection of risk factors did not differ for PTST-, TSTC and TST + young children and adults. This is the first study to rigorously characterize the epidemiologic risk profile of individuals with persistently negative TSTs despite close exposure to a person with TB. Additional studies are needed to characterize possible epidemiologic and host factors associated with this phenotype.

Mitochondrial polymorphism A10398G and Haplogroup I are associated with Fuchs' endothelial corneal dystrophy

PURPOSE

We investigated whether mitochondrial DNA (mtDNA) variants affect the susceptibility of Fuchs endothelial corneal dystrophy (FECD).

METHODS

Ten mtDNA variants defining European haplogroups were genotyped in a discovery dataset consisting of 530 cases and 498 controls of European descent from the Duke FECD cohort. Association tests for mtDNA markers and haplogroups were performed using logistic regression models with adjustment of age and sex. Subset analyses included controlling for additional effects of either the TCF4 SNP rs613872 or cigarette smoking. Our replication dataset was derived from the genome-wide association study (GWAS) of the FECD Genetics Consortium, where genotypes for three of 10 mtDNA markers were available. Replication analyses were performed to compare non-Duke cases to all GWAS controls (GWAS1, N = 3200), and to non-Duke controls (GWAS2, N = 3043).

RESULTS

The variant A10398G was significantly associated with FECD (odds ratio [OR] = 0.72; 95% confidence interval [CI] = [0.53, 0.98]; P = 0.034), and remains significant after adjusting for smoking status (min P = 0.012). This variant was replicated in GWAS1 (P = 0.019) and GWAS2 (P = 0.036). Haplogroup I was significantly associated with FECD (OR = 0.46; 95% CI = [0.22, 0.97]; P = 0.041) and remains significant after adjusting for the effect of smoking (min P = 0.008) or rs613872 (P = 0.034).

CONCLUSIONS

The 10398G allele and Haplogroup I appear to confer significant protective effects for FECD. The effect of A10398G and Haplogroup I to FECD is likely independent of the known TCF4 variant. More data are needed to decipher the interaction between smoking and mtDNA haplogroups.

Rare and common variants in extracellular matrix gene Fibrillin 2 (FBN2) are associated with macular degeneration

Neurodegenerative diseases affecting the macula constitute a major cause of incurable vision loss and exhibit considerable clinical and genetic heterogeneity, from early-onset monogenic disease to multifactorial late-onset age-related macular degeneration (AMD). As part of our continued efforts to define genetic causes of macular degeneration, we performed whole exome sequencing in four individuals of a two-generation family with autosomal dominant maculopathy and identified a rare variant p.Glu1144Lys in Fibrillin 2 (FBN2), a glycoprotein of the elastin-rich extracellular matrix (ECM). Sanger sequencing validated the segregation of this variant in the complete pedigree, including two additional affected and one unaffected individual. Sequencing of 192 maculopathy patients revealed additional rare variants, predicted to disrupt FBN2 function. We then undertook additional studies to explore the relationship of FBN2 to macular disease. We show that FBN2 localizes to Bruch's membrane and its expression appears to be reduced in aging and AMD eyes, prompting us to examine its relationship with AMD. We detect suggestive association of a common FBN2 non-synonymous variant, rs154001 (p.Val965Ile) with AMD in 10 337 cases and 11 174 controls (OR = 1.10; P-value = 3.79 × 10(-5)). Thus, it appears that rare and common variants in a single gene--FBN2--can contribute to Mendelian and complex forms of macular degeneration. Our studies provide genetic evidence for a key role of elastin microfibers and Bruch's membrane in maintaining blood-retina homeostasis and establish the importance of studying orphan diseases for understanding more common clinical phenotypes.

Genetic evidence for role of carotenoids in age-related macular degeneration in the Carotenoids in Age-Related Eye Disease Study (CAREDS)

PURPOSE

We tested variants in genes related to lutein and zeaxanthin status for association with age-related macular degeneration (AMD) in the Carotenoids in Age-Related Eye Disease Study (CAREDS).

METHODS

Of 2005 CAREDS participants, 1663 were graded for AMD from fundus photography and genotyped for 424 single nucleotide polymorphisms (SNPs) from 24 candidate genes for carotenoid status. Of 337 AMD cases 91% had early or intermediate AMD. The SNPs were tested individually for association with AMD using logistic regression. A carotenoid-related genetic risk model was built using backward selection and compared to existing AMD risk factors using the area under the receiver operating characteristic curve (AUC).

RESULTS

A total of 24 variants from five genes (BCMO1, BCO2, NPCL1L1, ABCG8, and FADS2) not previously related to AMD and four genes related to AMD in previous studies (SCARB1, ABCA1, APOE, and ALDH3A2) were associated independently with AMD, after adjusting for age and ancestry. Variants in all genes (not always the identical SNPs) were associated with lutein and zeaxanthin in serum and/or macula, in this or other samples, except for BCO2 and FADS2. A genetic risk score including nine variants significantly (P = 0.002) discriminated between AMD cases and controls beyond age, smoking, CFH Y402H, and ARMS2 A69S. The odds ratio (95% confidence interval) for AMD among women in the highest versus lowest quintile for the risk score was 3.1 (2.0-4.9).

CONCLUSIONS

Variants in genes related to lutein and zeaxanthin status were associated with AMD in CAREDS, adding to the body of evidence supporting a protective role of lutein and zeaxanthin in risk of AMD.

A genome-wide search for linkage of estimated glomerular filtration rate (eGFR) in the Family Investigation of Nephropathy and Diabetes (FIND)

Objective

Estimated glomerular filtration rate (eGFR), a measure of kidney function, is heritable, suggesting that genes influence renal function. Genes that influence eGFR have been identified through genome-wide association studies. However, family-based linkage approaches may identify loci that explain a larger proportion of the heritability. This study used genome-wide linkage and association scans to identify quantitative trait loci (QTL) that influence eGFR.

Methods

Genome-wide linkage and sparse association scans of eGFR were performed in families ascertained by probands with advanced diabetic nephropathy (DN) from the multi-ethnic Family Investigation of Nephropathy and Diabetes (FIND) study. This study included 954 African Americans (AA), 781 American Indians (AI), 614 European Americans (EA) and 1,611 Mexican Americans (MA). A total of 3,960 FIND participants were genotyped for 6,000 single nucleotide polymorphisms (SNPs) using the Illumina Linkage IVb panel. GFR was estimated by the Modification of Diet in Renal Disease (MDRD) formula.

Results

The non-parametric linkage analysis, accounting for the effects of diabetes duration and BMI, identified the strongest evidence for linkage of eGFR on chromosome 20q11 (log of the odds [LOD] = 3.34; P = 4.4×10⁻⁵) in MA and chromosome 15q12 (LOD = 2.84; P = 1.5×10⁻⁴) in EA. In all subjects, the strongest linkage signal for eGFR was detected on chromosome 10p12 (P = 5.5×10⁻⁴) at 44 cM near marker rs1339048. A subsequent association scan in both ancestry-specific groups and the entire population identified several SNPs significantly associated with eGFR across the genome.

Conclusion

The present study describes the localization of QTL influencing eGFR on 20q11 in MA, 15q21 in EA and 10p12 in the combined ethnic groups participating in the FIND study. Identification of causal genes/variants influencing eGFR, within these linkage and association loci, will open new avenues for functional analyses and development of novel diagnostic markers for DN.

Association of smoking and other risk factors with Fuchs' endothelial corneal dystrophy severity and corneal thickness

PURPOSE

We investigated effects of smoking and other risk factors on the development of advanced Fuchs' endothelial corneal dystrophy (FECD) and on central corneal thickness (CCT).

METHODS

Eyes from Caucasian probands, affected and unaffected family members, and unrelated controls matched for age from the FECD Genetics Multi-Center Study (n = 2044 subjects) were examined. Univariate and multivariate models, adjusted for family correlations, were used to determine the effect of smoking, sex, diabetes, and age on FECD case/control status and CCT.

RESULTS

In a multivariate model, sex and smoking were associated significantly with advanced FECD (grades 4-6) development (P = 0.016 and P = 0.047, respectively). Female sex increased odds by 34%. Smoking increased odds by 30%. In a multivariate model, diabetes was associated with an increase of 9.1 μm in average CCT (P = 0.021). Female sex was associated significantly with a decrease in average CCT by 6.9 μm (P = 0.015). Smoking had no significant effect on CCT in any model. As shown previously, advanced FECD was associated with large increases in CCT (31.4-94.2 μm).

CONCLUSIONS

Smoking was associated with an increased risk of advanced FECD and self-reported diabetes was associated with increased CCT. Further study of the impact of smoking and diabetes on FECD development and changes in corneal thickness is warranted.