T-cell receptor diversity in minimal change disease in the NEPTUNE study.

BACKGROUND: Minimal change disease (MCD) is the major cause of childhood idiopathic nephrotic syndrome, which is characterized by massive proteinuria and debilitating edema. Proteinuria in MCD is typically rapidly reversible with corticosteroid therapy, but relapses are common, and children often have many adverse events from the repeated courses of immunosuppressive therapy. The pathobiology of MCD remains poorly understood. Prior clinical observations suggest that abnormal T-cell function may play a central role in MCD pathogenesis. Based on these observations, we hypothesized that T-cell responses to specific exposures or antigens lead to a clonal expansion of T-cell subsets, a restriction in the T-cell repertoire, and an elaboration of specific circulating factors that trigger disease onset and relapses. METHODS: To test these hypotheses, we sequenced T-cell receptors in fourteen MCD, four focal segmental glomerulosclerosis (FSGS), and four membranous nephropathy (MN) patients with clinical data and blood samples drawn during active disease and during remission collected by the Nephrotic Syndrome Study Network (NEPTUNE). We calculated several T-cell receptor diversity metrics to assess possible differences between active disease and remission states in paired samples. RESULTS: Median productive clonality did not differ between MCD active disease (0.0083; range: 0.0042, 0.0397) and remission (0.0088; range: 0.0038, 0.0369). We did not identify dominant clonotypes in MCD active disease, and few clonotypes were shared with FSGS and MN patients. CONCLUSIONS: While these data do not support an obvious role of the adaptive immune system T-cells in MCD pathogenesis, further study is warranted given the limited sample size. A higher resolution version of the Graphical abstract is available as Supplementary information.

Polygenic Risk Scores.

As genome-wide association studies have continued to identify loci associated with complex traits, the implications of and necessity for proper use of these findings, including prediction of disease risk, have become apparent. Many complex diseases have numerous associated loci with detectable effects implicating risk for or protection from disease. A common contemporary approach to using this information for disease prediction is through the application of genetic risk scores. These scores estimate an individual's liability for a specific outcome by aggregating the effects of associated loci into a single measure as described in the previous version of this article. Although genetic risk scores have traditionally included variants that meet criteria for genome-wide significance, an extension known as the polygenic risk score has been developed to include the effects of more variants across the entire genome. Here, we describe common methods and software packages for calculating and interpreting polygenic risk scores. In this revised version of the article, we detail information that is needed to perform a polygenic risk score analysis, considerations for planning the analysis and interpreting results, as well as discussion of the limitations based on the choices made. We also provide simulated sample data and a walkthrough for four different polygenic risk score software. © 2021 Wiley Periodicals LLC.

Author Correction: Cross-ancestry genome-wide association analysis of corneal thickness strengthens link between complex and Mendelian eye diseases.

Emmanuelle Souzeau, who contributed to analysis of data, was inadvertently omitted from the author list in the originally published version of this Article. This has now been corrected in both the PDF and HTML versions of the Article.

Somatic T-cell Receptor Diversity in a Chronic Kidney Disease PatientPopulation Linked to Electronic Health Records.

Germline and somatic genomic variation represent the bulk of 'omics data available for precision medicine research. These data, however, may fail to capture the dynamic biological processes that underlie disease development, particularly for chronic diseases of aging such as chronic kidney disease (CKD). To demonstrate the value of additional dynamic precision medicine data, we sequenced somatic T-cell receptor rearrangements, markers of the adaptive immune response, from genomic DNA collected during a clinical encounter from 15 participants with CKD and associated co-morbidities. Participants were consented as part of a larger precision medicine research project at the MetroHealth System, a large urban public hospital in Cleveland, Ohio. Despite the limited sample size, we observed reduced T-cell receptor diversity in relation to biomarkers (creatinine and BUN) of CKD status in this older and mostly African American sample. Overall, these data suggest a relationship between advanced CKD and premature aging of the adaptive immune system and highlight the potential of dynamic 'omic data to generate novel hypotheses about disease mechanisms and unique opportunities for precision medicine applications.

Cross-ancestry genome-wide association analysis of corneal thickness strengthens link between complex and Mendelian eye diseases.

Central corneal thickness (CCT) is a highly heritable trait associated with complex eye diseases such as keratoconus and glaucoma. We perform a genome-wide association meta-analysis of CCT and identify 19 novel regions. In addition to adding support for known connective tissue-related pathways, pathway analyses uncover previously unreported gene sets. Remarkably, >20% of the CCT-loci are near or within Mendelian disorder genes. These included FBN1, ADAMTS2 and TGFB2 which associate with connective tissue disorders (Marfan, Ehlers-Danlos and Loeys-Dietz syndromes), and the LUM-DCN-KERA gene complex involved in myopia, corneal dystrophies and cornea plana. Using index CCT-increasing variants, we find a significant inverse correlation in effect sizes between CCT and keratoconus (r = -0.62, P = 5.30 × 10-5) but not between CCT and primary open-angle glaucoma (r = -0.17, P = 0.2). Our findings provide evidence for shared genetic influences between CCT and keratoconus, and implicate candidate genes acting in collagen and extracellular matrix regulation.

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.

Systems genetics identifies a role for Cacna2d1 regulation in elevated intraocular pressure and glaucoma susceptibility.

Glaucoma is a multi-factorial blinding disease in which genetic factors play an important role. Elevated intraocular pressure is a highly heritable risk factor for primary open angle glaucoma and currently the only target for glaucoma therapy. Our study helps to better understand underlying genetic and molecular mechanisms that regulate intraocular pressure, and identifies a new candidate gene, Cacna2d1, that modulates intraocular pressure and a promising therapeutic, pregabalin, which binds to CACNA2D1 protein and lowers intraocular pressure significantly. Because our study utilizes a genetically diverse population of mice with known sequence variants, we are able to determine that the intraocular pressure-lowering effect of pregabalin is dependent on the Cacna2d1 haplotype. Using human genome-wide association study (GWAS) data, evidence for association of a CACNA2D1 single-nucleotide polymorphism and primary open angle glaucoma is found. Importantly, these results demonstrate that our systems genetics approach represents an efficient method to identify genetic variation that can guide the selection of therapeutic targets.

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.

New insights into the genetics of primary open-angle glaucoma based on meta-analyses of intraocular pressure and optic disc characteristics.

Primary open-angle glaucoma (POAG), the most common optic neuropathy, is a heritable disease. Siblings of POAG cases have a ten-fold increased risk of developing the disease. Intraocular pressure (IOP) and optic nerve head characteristics are used clinically to predict POAG risk. We conducted a genome-wide association meta-analysis of IOP and optic disc parameters and validated our findings in multiple sets of POAG cases and controls. Using imputation to the 1000 genomes (1000G) reference set, we identified 9 new genomic regions associated with vertical cup-disc ratio (VCDR) and 1 new region associated with IOP. Additionally, we found 5 novel loci for optic nerve cup area and 6 for disc area. Previously it was assumed that genetic variation influenced POAG either through IOP or via changes to the optic nerve head; here we present evidence that some genomic regions affect both IOP and the disc parameters. We characterized the effect of the novel loci through pathway analysis and found that pathways involved are not entirely distinct as assumed so far. Further, we identified a novel association between CDKN1A and POAG. Using a zebrafish model we show that six6b (associated with POAG and optic nerve head variation) alters the expression of cdkn1a. In summary, we have identified several novel genes influencing the major clinical risk predictors of POAG and showed that genetic variation in CDKN1A is important in POAG risk.

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.

Whole exome sequencing of extreme age-related macular degeneration phenotypes.

PURPOSE: Demographic, environmental, and genetic risk factors for age-related macular degeneration (AMD) have been identified; however, a substantial portion of the variance in AMD disease risk and heritability remains unexplained. To identify AMD risk variants and generate hypotheses for future studies, we performed whole exome sequencing for 75 individuals whose phenotype was not well predicted by their genotype at known risk loci. We hypothesized that these phenotypically extreme individuals were more likely to carry rare risk or protective variants with large effect sizes. METHODS: A genetic risk score was calculated in a case-control set of 864 individuals (467 AMD cases, 397 controls) based on 19 common (≥1% minor allele frequency, MAF) single nucleotide variants previously associated with the risk of advanced AMD in a large meta-analysis of advanced cases and controls. We then selected for sequencing 39 cases with bilateral choroidal neovascularization with the lowest genetic risk scores to detect risk variants and 36 unaffected controls with the highest genetic risk score to detect protective variants. After minimizing the influence of 19 common genetic risk loci on case-control status, we targeted single variants of large effect and the aggregate effect of weaker variants within genes and pathways. Single variant tests were conducted on all variants, while gene-based and pathway analyses were conducted on three subsets of data: 1) rare (≤1% MAF in the European population) stop, splice, or damaging missense variants, 2) all rare variants, and 3) all variants. All analyses controlled for the effects of age and sex. RESULTS: No variant, gene, or pathway outside regions known to be associated with risk for advanced AMD reached genome-wide significance. However, we identified several variants with substantial differences in allele frequency between cases and controls with strong additive effects on affection status after controlling for age and sex. Protective effects trending toward significance were detected at two loci identified in single-variant analyses: an intronic variant in FBLN7 (the gene encoding fibulin 7) and at three variants near pyridoxal (pyridoxine, vitamin B6) kinase (PDXK). Aggregate rare-variant analyses suggested evidence for association at ASRGL1, a gene previously linked to photoreceptor cell death, and at BSDC1. In known AMD loci we also identified 29 novel or rare damaging missense or stop/splice variants in our sample of cases and controls. CONCLUSIONS: Identified variants and genes may highlight regions important in the pathogenesis of AMD and are key targets for replication.

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 worldwide. To identify new susceptibility loci, we performed meta-analysis on genome-wide association study (GWAS) results from eight independent studies from the United States (3,853 cases and 33,480 controls) and investigated the most significantly associated SNPs in two Australian studies (1,252 cases and 2,592 controls), three European studies (875 cases and 4,107 controls) and a Singaporean Chinese study (1,037 cases and 2,543 controls). A meta-analysis of the top SNPs identified three new associated loci: rs35934224[T] in TXNRD2 (odds ratio (OR) = 0.78, P = 4.05 × 10(-11)) encoding a mitochondrial protein required for redox homeostasis; rs7137828[T] in ATXN2 (OR = 1.17, P = 8.73 × 10(-10)); and rs2745572[A] upstream of FOXC1 (OR = 1.17, P = 1.76 × 10(-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 new targets for preventative therapies.

A large genome-wide association study of age-related macular degeneration highlights contributions of rare and common variants.

Advanced age-related macular degeneration (AMD) is the leading cause of blindness in the elderly, with limited therapeutic options. Here we report on a study of >12 million variants, including 163,714 directly genotyped, mostly rare, protein-altering variants. Analyzing 16,144 patients and 17,832 controls, we identify 52 independently associated common and rare variants (P < 5 × 10(-8)) distributed across 34 loci. Although wet and dry AMD subtypes exhibit predominantly shared genetics, we identify the first genetic association signal specific to wet AMD, near MMP9 (difference P value = 4.1 × 10(-10)). Very rare coding variants (frequency <0.1%) in CFH, CFI and TIMP3 suggest causal roles for these genes, as does a splice variant in SLC16A8. Our results support the hypothesis that rare coding variants can pinpoint causal genes within known genetic loci and illustrate that applying the approach systematically to detect new loci requires extremely large sample sizes.

Meta-analysis of genome-wide association studies identifies novel loci that influence cupping and the glaucomatous process.

Glaucoma is characterized by irreversible optic nerve degeneration and is the most frequent cause of irreversible blindness worldwide. Here, the International Glaucoma Genetics Consortium conducts a meta-analysis of genome-wide association studies of vertical cup-disc ratio (VCDR), an important disease-related optic nerve parameter. In 21,094 individuals of European ancestry and 6,784 individuals of Asian ancestry, we identify 10 new loci associated with variation in VCDR. In a separate risk-score analysis of five case-control studies, Caucasians in the highest quintile have a 2.5-fold increased risk of primary open-angle glaucoma as compared with those in the lowest quintile. This study has more than doubled the known loci associated with optic disc cupping and will allow greater understanding of mechanisms involved in this common blinding condition.

Common variants near ABCA1, AFAP1 and GMDS confer risk of primary open-angle glaucoma.

Primary open-angle glaucoma (POAG) is a major cause of irreversible blindness worldwide. We performed a genome-wide association study in an Australian discovery cohort comprising 1,155 cases with advanced POAG and 1,992 controls. We investigated the association of the top SNPs from the discovery stage in two Australian replication cohorts (932 cases and 6,862 controls total) and two US replication cohorts (2,616 cases and 2,634 controls total). Meta-analysis of all cohorts identified three loci newly associated with development of POAG. These loci are located upstream of ABCA1 (rs2472493[G], odds ratio (OR) = 1.31, P = 2.1 × 10(-19)), within AFAP1 (rs4619890[G], OR = 1.20, P = 7.0 × 10(-10)) and within GMDS (rs11969985[G], OR = 1.31, P = 7.7 × 10(-10)). Using RT-PCR and immunolabeling, we show that these genes are expressed within human retina, optic nerve and trabecular meshwork and that ABCA1 and AFAP1 are also expressed in retinal ganglion cells.

Genome-wide analysis of multi-ancestry cohorts identifies new loci influencing intraocular pressure and susceptibility to glaucoma.

Elevated intraocular pressure (IOP) is an important risk factor in developing glaucoma, and variability in IOP might herald glaucomatous development or progression. We report the results of a genome-wide association study meta-analysis of 18 population cohorts from the International Glaucoma Genetics Consortium (IGGC), comprising 35,296 multi-ancestry participants for IOP. We confirm genetic association of known loci for IOP and primary open-angle glaucoma (POAG) and identify four new IOP-associated loci located on chromosome 3q25.31 within the FNDC3B gene (P = 4.19 × 10(-8) for rs6445055), two on chromosome 9 (P = 2.80 × 10(-11) for rs2472493 near ABCA1 and P = 6.39 × 10(-11) for rs8176693 within ABO) and one on chromosome 11p11.2 (best P = 1.04 × 10(-11) for rs747782). Separate meta-analyses of 4 independent POAG cohorts, totaling 4,284 cases and 95,560 controls, showed that 3 of these loci for IOP were also associated with POAG.

The impact of the human genome project on complex disease.

In the decade that has passed since the initial release of the Human Genome, numerous advancements in science and technology within and beyond genetics and genomics have been encouraged and enhanced by the availability of this vast and remarkable data resource. Progress in understanding three common, complex diseases: age-related macular degeneration (AMD), Alzheimer's disease (AD), and multiple sclerosis (MS), are three exemplars of the incredible impact on the elucidation of the genetic architecture of disease. The approaches used in these diseases have been successfully applied to numerous other complex diseases. For example, the heritability of AMD was confirmed upon the release of the first genome-wide association study (GWAS) along with confirmatory reports that supported the findings of that state-of-the art method, thus setting the foundation for future GWAS in other heritable diseases. Following this seminal discovery and applying it to other diseases including AD and MS, the genetic knowledge of AD expanded far beyond the well-known APOE locus and now includes more than 20 loci. MS genetics saw a similar increase beyond the HLA loci and now has more than 100 known risk loci. Ongoing and future efforts will seek to define the remaining heritability of these diseases; the next decade could very well hold the key to attaining this goal.

Hypothesis-independent pathway analysis implicates GABA and acetyl-CoA metabolism in primary open-angle glaucoma and normal-pressure glaucoma.

Primary open-angle glaucoma (POAG) is a leading cause of blindness worldwide. Using genome-wide association single-nucleotide polymorphism data from the Glaucoma Genes and Environment study and National Eye Institute Glaucoma Human Genetics Collaboration comprising 3,108 cases and 3,430 controls, we assessed biologic pathways as annotated in the KEGG database for association with risk of POAG. After correction for genic overlap among pathways, we found 4 pathways, butanoate metabolism (hsa00650), hematopoietic cell lineage (hsa04640), lysine degradation (hsa00310) and basal transcription factors (hsa03022) related to POAG with permuted p < 0.001. In addition, the human leukocyte antigen (HLA) gene family was significantly associated with POAG (p < 0.001). In the POAG subset with normal-pressure glaucoma (NPG), the butanoate metabolism pathway was also significantly associated (p < 0.001) as well as the MAPK and Hedgehog signaling pathways (hsa04010 and hsa04340), glycosaminoglycan biosynthesis-heparan sulfate pathway (hsa00534) and the phenylalanine, tyrosine and tryptophan biosynthesis pathway (hsa0400). The butanoate metabolism pathway overall, and specifically the aspects of the pathway that contribute to GABA and acetyl-CoA metabolism, was the only pathway significantly associated with both POAG and NPG. Collectively these results implicate GABA and acetyl-CoA metabolism in glaucoma pathogenesis, and suggest new potential therapeutic targets.

Analysis of coding variants identified from exome sequencing resources for association with diabetic and non-diabetic nephropathy in African Americans.

Prior studies have identified common genetic variants influencing diabetic and non-diabetic nephropathy, diseases which disproportionately affect African Americans. Recently, exome sequencing techniques have facilitated identification of coding variants on a genome-wide basis in large samples. Exonic variants in known or suspected end-stage kidney disease (ESKD) or nephropathy genes can be tested for their ability to identify association either singly or in combination with known associated common variants. Coding variants in genes with prior evidence for association with ESKD or nephropathy were identified in the NHLBI-ESP GO database and genotyped in 5,045 African Americans (3,324 cases with type 2 diabetes associated nephropathy [T2D-ESKD] or non-T2D ESKD, and 1,721 controls) and 1,465 European Americans (568 T2D-ESKD cases and 897 controls). Logistic regression analyses were performed to assess association, with admixture and APOL1 risk status incorporated as covariates. Ten of 31 SNPs were associated in African Americans; four replicated in European Americans. In African Americans, SNPs in OR2L8, OR2AK2, C6orf167 (MMS22L), LIMK2, APOL3, APOL2, and APOL1 were nominally associated (P = 1.8 × 10(-4)-0.044). Haplotype analysis of common and coding variants increased evidence of association at the OR2L13 and APOL1 loci (P = 6.2 × 10(-5) and 4.6 × 10(-5), respectively). SNPs replicating in European Americans were in OR2AK2, LIMK2, and APOL2 (P = 0.0010-0.037). Meta-analyses highlighted four SNPs associated in T2D-ESKD and all-cause ESKD. Results from this study suggest a role for coding variants in the development of diabetic, non-diabetic, and/or all-cause ESKD in African Americans and/or European Americans.

The Role of Copy Number Variation in African Americans with Type 2 Diabetes-Associated End Stage Renal Disease.

This study investigated the association of copy number variants (CNVs) in type 2 diabetes (T2D) and T2D-associated end-stage renal disease (ESRD) in African Americans. Using the Affymetrix 6.0 array, >900,000 CNV probes spanning the genome were interrogated in 965 African Americans with T2D-ESRD and 1029 non-diabetic African American controls. Previously identified and novel CNVs were separately analyzed and were evaluated for insertion/deletion status and then used as predictors in a logistic regression model to test for association. One common CNV insertion on chromosome 1 was significantly associated with T2D-ESRD (p=6.17×10-5, OR=1.63) after multiple comparison correction. This CNV region encompasses the genes AMY2A and AMY2B, which encode amylase isoenzymes produced by the pancreas. Additional common and novel CNVs approaching significance with disease were also detected. These exploratory results require further replication but suggest the involvement of the AMY2A/AMY2B CNV in T2D and/or T2D-ESRD, and indicate that CNVs may contribute to susceptibility for these diseases.