Christopher E. Gillies

Integrative Genomics Identifies Novel Associations with APOL1 Risk Genotypes in Black NEPTUNE Subjects

APOL1 variants have been associated with renal phenotypes in blacks. To refine clinical outcomes and discover mechanisms of APOL1-associated kidney injury, we analyzed clinical and genomic datasets derived from 90 black subjects in the Nephrotic Syndrome Study Network (NEPTUNE), stratified by APOL1 risk genotype. Ninety subjects with proteinuria ≥0.5 g/d were enrolled at first biopsy for primary nephrotic syndrome and followed. Clinical outcomes were determined, and renal histomorphometry and sequencing of Mendelian nephrotic syndrome genes were performed. APOL1 variants were genotyped, and glomerular and tubulointerstitial transcriptomes from protocol renal biopsy cores were analyzed for differential and correlative gene expression. Analyses were performed under the recessive model (high-risk genotype defined by two risk alleles). APOL1 high-risk genotype was significantly associated with a 17 ml/min per 1.73 m(2) lower eGFR and a 69% reduction in the probability of complete remission at any time, independent of histologic diagnosis. Neither APOL1 risk group was enriched for Mendelian mutations. On renal biopsy, high-risk genotype was associated with increased fractional interstitial area, interstitial fibrosis, and tubular atrophy. Risk genotype was not associated with intrarenal APOL1 mRNA expression levels. Differential expression analysis demonstrated an increased steady-state level of five genes associated with the high-risk genotype (CXCL9, CXCL11, and UBD in glomerulus; SNOR14B and MUC13 in tubulointerstitium). APOL1 tubulointerstitial coexpression analysis showed coexpression of APOL1 mRNA levels with a group of intrarenal transcripts that together were associated with increased interstitial fibrosis and tubular atrophy. These data indicate the high-risk APOL1 genotype confers renal risk across histopathologic diagnoses.

Genetic Drivers of Kidney Defects in the DiGeorge Syndrome

Background The DiGeorge syndrome, the most common of the microdeletion syndromes, affects multiple organs, including the heart, the nervous system, and the kidney. It is caused by deletions on chromosome 22q11.2; the genetic driver of the kidney defects is unknown. Methods We conducted a genomewide search for structural variants in two cohorts: 2080 patients with congenital kidney and urinary tract anomalies and 22,094 controls. We performed exome and targeted resequencing in samples obtained from 586 additional patients with congenital kidney anomalies. We also carried out functional studies using zebrafish and mice. Results We identified heterozygous deletions of 22q11.2 in 1.1% of the patients with congenital kidney anomalies and in 0.01% of population controls (odds ratio, 81.5; P=4.5×10‐14). We localized the main drivers of renal disease in the DiGeorge syndrome to a 370‐kb region containing nine genes. In zebrafish embryos, an induced loss of function in snap29, aifm3, and crkl resulted in renal defects; the loss of crkl alone was sufficient to induce defects. Five of 586 patients with congenital urinary anomalies had newly identified, heterozygous protein‐altering variants, including a premature termination codon, in CRKL. The inactivation of Crkl in the mouse model induced developmental defects similar to those observed in patients with congenital urinary anomalies. Conclusions We identified a recurrent 370‐kb deletion at the 22q11.2 locus as a driver of kidney defects in the DiGeorge syndrome and in sporadic congenital kidney and urinary tract anomalies. Of the nine genes at this locus, SNAP29, AIFM3, and CRKL appear to be critical to the phenotype, with haploinsufficiency of CRKL emerging as the main genetic driver. (Funded by the National Institutes of Health and others.)

Using Population Genetics to Interrogate the Monogenic Nephrotic Syndrome Diagnosis in a Case Cohort

To maximize clinical benefits of genetic screening of patients with nephrotic syndrome (NS) to diagnose monogenic causes, reliably distinguishing NS-causing variants from the background of rare, noncausal variants prevalent in all genomes is vital. To determine the prevalence of monogenic NS in a North American case cohort while accounting for background prevalence of genetic variation, we sequenced 21 implicated monogenic NS genes in 312 participants from the Nephrotic Syndrome Study Network and 61 putative controls from the 1000 Genomes Project (1000G). These analyses were extended to available sequence data from approximately 2500 subjects from the 1000G. A typical pathogenicity filter identified causal variants for NS in 4.2% of patients and 5.8% of subjects from the 1000G. We devised a more stringent pathogenicity filtering strategy, reducing background prevalence of causal variants to 1.5%. When applying this stringent filter to patients, prevalence of monogenic NS was 2.9%; of these patients, 67% were pediatric, and 44% had FSGS on biopsy. The rate of complete remission did not associate with monogenic classification. Thus, we identified factors contributing to inaccurate monogenic classification of NS and developed a more accurate variant filtering strategy. The prevalence and clinical correlates of monogenic NS in this sporadically affected cohort differ substantially from those reported for patients referred for genetic analysis. Particularly in unselected, population-based cases, considering putative causal variants in known NS genes from a probabilistic rather than a deterministic perspective may be more precise. We also introduce GeneVetter, a web tool for monogenic assessment of rare disease.

APOL1-Associated glomerular disease among African-American children: A collaboration of the chronic kidney disease in children (CKiD) and nephrotic syndrome study network (NEPTUNE) cohorts

Background Individuals of African ancestry harboring two variant alleles within apolipoprotein L1 ( APOL1 ) are classified with a high-risk (HR) genotype. Adults with an HR genotype have increased risk of focal segmental glomerulosclerosis and chronic kidney disease compared with those with a low-risk (LR) genotype (0 or 1 variants). The role of APOL1 risk genotypes in children with glomerular disease is less well known. Methods This study characterized 104 African-American children with a glomerular disease by APOL1 genotype in two cohorts: the Chronic Kidney Disease in Children (CKiD) and Nephrotic Syndrome Study Network (NEPTUNE). Results Among these subjects, 46% had an HR genotype with a similar age at cohort enrollment. For APOL1 HR children, the median age of disease onset was older (CKiD: 4.5 versus 11.5 years for LR versus HR; NEPTUNE: 11 versus 14 years for LR versus HR, respectively) and preterm birth was more common [CKiD: 27 versus 4%; NEPTUNE: 26 versus 12%; combined odds ratio 4.6 (95% confidence interval: 1.4, 15.5)]. Within studies, HR children had lower initial estimated glomerular filtration rate (eGFR) (CKiD: 53 versus 69 mL/min/1.73 m 2 ; NEPTUNE: 74 versus 94 mL/min/1.73 m 2 ). Longitudinal eGFR decline was faster among HR children versus LR (CKiD: -18 versus -8% per year; NEPTUNE: -13 versus -3% per year). Conclusions Children with an HR genotype in CKiD and NEPTUNE seem to have a more aggressive form of glomerular disease, in part due to a higher prevalence of focal segmental glomerulosclerosis. These consistent findings across independent cohorts suggest a common natural history for children with APOL1 -associated glomerular disease. Further study is needed to determine the generalizability of these findings.

Evaluating Mendelian nephrotic syndrome genes for evidence for risk alleles or oligogenicity that explain heritability

BackgroundMore than 30 genes can harbor rare exonic variants sufficient to cause nephrotic syndrome (NS), and the number of genes implicated in monogenic NS continues to grow. However, outside the first year of life, the majority of affected patients, particularly in ancestrally mixed populations, do not have a known monogenic form of NS. Even in those children classified with a monogenic form of NS, there is phenotypic heterogeneity. Thus, we have only discovered a fraction of the heritability of NS—the underlying genetic factors contributing to phenotypic variation. Part of the “missing heritability” for NS has been posited to be explained by patients harboring coding variants across one or more previously implicated NS genes, insufficient to cause NS in a classical Mendelian manner, but that nonetheless have a sufficient impact on protein function to cause disease. However, systematic evaluation in patients with NS for rare or low-frequency risk alleles within single genes, or in combination across genes (“oligogenicity”), has not been reported. To determine whether, compared with a reference population, patients with NS have either a significantly increased burden of protein-altering variants (“risk-alleles”), or a unique combination of them (“oligogenicity”), in a set of 21 genes implicated in Mendelian forms of NS.MethodsIn 303 patients with NS enrolled in the Nephrotic Syndrome Study Network (NEPTUNE), we performed targeted amplification paired with next-generation sequencing of 21 genes implicated in monogenic NS. We created a high-quality variant call set and compared it with a variant call set of the same genes in a reference population composed of 2,535 individuals from phase 3 of the 1000 Genomes Project. We created both a “stringent” and a “relaxed” pathogenicity-filtering pipeline, applied them to both cohorts, and computed the burden of variants in the entire gene set per cohort, the burden of variants in the entire gene set per individual, the burden of variants within a single gene per cohort, and unique combinations of variants across two or more genes per cohort.ResultsWith few exceptions when using the relaxed filter, and which are likely the result of confounding by population stratification, NS patients did not have a significantly increased burden of variants in Mendelian NS genes in comparison to a reference cohort, nor was there any evidence for oligogenicity. This was true when using both the relaxed and the stringent variant pathogenicity filter.ConclusionIn our study, there were no significant differences in the burden or particular combinations of low-frequency or rare protein-altering variants in a previously implicated Mendelian NS genes cohort between North American patients with NS and a reference population. Studies in larger independent cohorts or meta-analyses are needed to assess the generalizability of our discoveries and also address whether there is in fact small but significant enrichment of risk alleles or oligogenicity in NS cases that was undetectable with this current sample size. It is still possible that rare protein-altering variants in these genes, insufficient to cause Mendelian disease, still contribute to NS as risk alleles and/or via oligogenicity. However, we suggest that more accurate bioinformatic analyses and the incorporation of functional assays would be necessary to identify bona fide instances of this form of genetic architecture as a contributor to the heritability of NS.

Gene-level Integrated Metric of negative Selection (GIMS) Prioritizes Candidate Genes for Nephrotic Syndrome

Nephrotic syndrome (NS) gene discovery efforts are now occurring in small kindreds and cohorts of sporadic cases. Power to identify causal variants in these groups beyond a statistical significance threshold is challenging due to small sample size and/or lack of family information. There is a need to develop novel methods to identify NS-associated variants. One way to determine putative functional relevance of a gene is to measure its strength of negative selection, as variants in genes under strong negative selection are more likely to be deleterious. We created a gene-level, integrated metric of negative selection (GIMS) score for 20,079 genes by combining multiple comparative genomics and population genetics measures. To understand the utility of GIMS for NS gene discovery, we examined this score in a diverse set of NS-relevant gene sets. These included genes known to cause monogenic forms of NS in humans as well as genes expressed in the cells of the glomerulus and, particularly, the podocyte. We found strong negative selection in the following NS-relevant gene sets: (1) autosomal-dominant Mendelian focal segmental glomerulosclerosis (FSGS) genes (p= 0.03 compared to reference), (2) glomerular expressed genes (p = 4×10-23), and (3) predicted podocyte genes (p = 3×10-9). Eight genes causing autosomal dominant forms of FSGS had a stronger combined score of negative selection and podocyte enrichment as compared to all other genes (p=1 x 10-3). As a whole, recessive FSGS genes were not enriched for negative selection. Thus, we also created a transcript-level, integrated metric of negative selection (TIMS) to quantify negative selection on an isoform level. These revealed transcripts of known autosomal recessive disease-causing genes that were nonetheless under strong selection. We suggest that a filtering strategy that includes measuring negative selection on a gene or isoform level could aid in identifying NS-related genes. Our GIMS and TIMS scores are available at http://glom.sph.umich.edu/GIMS/.

An investigation of APOL1 risk genotypes and preterm birth in African American population cohorts.

Background Two genetic variants in apolipoprotein L1 (APOL1) are associated with increased risk of focal segmental glomerulosclerosis as well as other glomerular phenotypes. These risk variants are common in individuals of African ancestry but absent in other racial groups. Yet, the majority of individuals with two APOL1 risk alleles [high-risk (HR) genotype] do not have renal disease. It is critical to identify environmental and secondary genetic influences that, when combined with these alleles, lead to kidney disease. In a recent study of black children with glomerular disease enrolled in the Nephrotic Syndrome Study Network (NEPTUNE) and Chronic Kidney Disease in Children Study (n = 104), we found that subjects with an HR genotype had a 4.6-fold increase in the odds of preterm birth as compared to those with a low risk (LR) genotype [odds ratio 4.6 (CI 1.4-15.5)]. There are known racial disparities in preterm birth, which itself is a known risk factor for chronic kidney disease and focal segmental glomerulosclerosis. Thus, we questioned whether an HR APOL1 genotype is associated with prematurity in the general African American population. Methods We analyzed two publically available genetic datasets of preterm birth in African Americans, including 867 infants and 519 mothers from the Gene Environment Association Studies (GENEVA) study of preterm delivery and 960 mothers from the Boston Medical Center genome-wide association study of preterm birth. We performed multivariable analyses testing for association between HR APOL1 and birth outcomes. Results In both studies, there was no association between HR APOL1 in mothers and prematurity, gestational age or birthweight. Additionally, in the GENEVA study, we saw no association between infant HR APOL1 and prematurity, gestational age or birthweight. Conclusion From these data, we conclude that the previously observed association between HR APOL1 and prematurity is specific to those with glomerular disease, suggesting prematurity may act as an additional risk factor in APOL1-associated renal disease.

Whole Exome Sequencing Reveals Novel PHEX Splice Site Mutations in Patients with Hypophosphatemic Rickets

Objective Hypophosphatemic rickets (HR) is a heterogeneous genetic phosphate wasting disorder. The disease is most commonly caused by mutations in the PHEX gene located on the X-chromosome or by mutations in CLCN5, DMP1, ENPP1, FGF23, and SLC34A3. The aims of this study were to perform molecular diagnostics for four patients with HR of Indian origin (two independent families) and to describe their clinical features. Methods We performed whole exome sequencing (WES) for the affected mother of two boys who also displayed the typical features of HR, including bone malformations and phosphate wasting. B-lymphoblast cell lines were established by EBV transformation and subsequent RT-PCR to investigate an uncommon splice site variant found by WES. An in silico analysis was done to obtain accurate nucleotide frequency occurrences of consensus splice positions other than the canonical sites of all human exons. Additionally, we applied direct Sanger sequencing for all exons and exon/intron boundaries of the PHEX gene for an affected girl from an independent second Indian family. Results WES revealed a novel PHEX splice acceptor mutation in intron 9 (c.1080-3C>A) in a family with 3 affected individuals with HR. The effect on splicing of this mutation was further investigated by RT-PCR using RNA obtained from a patient’s EBV-transformed lymphoblast cell line. RT-PCR revealed an aberrant splice transcript skipping exons 10-14 which was not observed in control samples, confirming the diagnosis of X-linked dominant hypophosphatemia (XLH). The in silico analysis of all human splice sites adjacent to all 327,293 exons across 81,814 transcripts among 20,345 human genes revealed that cytosine is, with 64.3%, the most frequent nucleobase at the minus 3 splice acceptor position, followed by thymidine with 28.7%, adenine with 6.3%, and guanine with 0.8%. We generated frequency tables and pictograms for the extended donor and acceptor splice consensus regions by analyzing all human exons. Direct Sanger sequencing of all PHEX exons in a sporadic case with HR from the Indian subcontinent revealed an additional novel PHEX mutation (c.1211_1215delACAAAinsTTTACAT, p.Asp404Valfs*5, de novo) located in exon 11. Conclusions Mutation analyses revealed two novel mutations and helped to confirm the clinical diagnoses of XLH in two families from India. WES helped to analyze all genes implicated in the underlying disease complex. Mutations at splice positions other than the canonical key sites need further functional investigation to support the assertion of pathogenicity.

NPHS2 V260E Is a Frequent Cause of Steroid-Resistant Nephrotic Syndrome in Black South African Children

Introduction In South Africa (SA), steroid-resistant nephrotic syndrome (SRNS) is more frequent in black than in Indian children. Methods Seeking a genetic basis for this disparity, we enrolled 33 Indian and 31 black children with steroid-sensitive nephrotic syndrome (SSNS) and SRNS from KwaZulu-Natal, SA; SRNS children underwent kidney biopsy. We sequenced NPHS2 and genotyped APOL1 in 15 SSNS and 64 SRNS unrelated patients and 104 controls and replicated results in 18 black patients with steroid-resistant focal segmental glomerulosclerosis (SR-FSGS). Known FSGS genes (n = 21) were sequenced in a subset of patients. Results Homozygosity for NPHS2 V260E was found in 8 of 30 black children with SRNS (27%); all 260E/E carriers had SR-FSGS. Combining SR-FSGS patients from the 2 groups, 14 of 42 (33%) were homozygous for V260E. One black control was heterozygous for V260E; no Indian patients or controls were carriers. Haplotype analysis indicated that homozygosity for V260E was not explained by cryptic consanguinity. Children with NPHS2 260E/E developed SRNS at earlier age than noncarriers (34 vs. 78 months, P = 0.01), and none achieved partial or complete remission (0% vs. 47%, P = 0.002). APOL1 variants did not associate with NS. Sequencing FSGS genes identified a CD2AP predicted pathogenic variant in the heterozygous state in 1 Indian case with SR-FSGS. Conclusion NPHS2 260E/E was present in one-third of black FSGS patients, was absent in black controls and Indian patients, and affected patients were unresponsive to therapy. Genotyping V260E in black children from South Africa with NS will identify a substantial group with SR-FSGS, potentially sparing these children biopsy and ineffective steroid treatment.