Abstract
Congenital anomalies of the kidney and urinary tract (CAKUT) occur in three to six of 1000 live births, represent about 20% of the prenatally detected anomalies, and constitute the main cause of CKD in children. These disorders are phenotypically and genetically heterogeneous. Monogenic causes of CAKUT in humans and mice have been identified. However, despite high-throughput sequencing studies, the cause of the disease remains unknown in most patients, and several studies support more complex inheritance and the role of environmental factors and/or epigenetics in the pathophysiology of CAKUT. Here, we report the targeted exome sequencing of 330 genes, including genes known to be involved in CAKUT and candidate genes, in a cohort of 204 unrelated patients with CAKUT; 45% of the patients were severe fetal cases. We identified pathogenic mutations in 36 of 204 (17.6%) patients. These mutations included five de novo heterozygous loss of function mutations/deletions in the PBX homeobox 1 gene (PBX1), a gene known to have a crucial role in kidney development. In contrast, the frequency of SOX17 and DSTYK variants recently reported as pathogenic in CAKUT did not indicate causality. These findings suggest that PBX1 is involved in monogenic CAKUT in humans and call into question the role of some gene variants recently reported as pathogenic in CAKUT. Targeted exome sequencing also proved to be an efficient and cost-effective strategy to identify pathogenic mutations and deletions in known CAKUT genes.
Congenital anomalies of the kidney and urinary tract (CAKUT) refer to a heterogeneous group of malformations of the kidney and the urinary tract due to defects in embryonic kidney development, a complex process that involves reciprocal interaction between the ureteric bud and the metanephric mesenchymal tissue. Prevalence of CAKUT has been estimated between three and six per 1000 live births. These kidney defects are frequently detected in utero by prenatal ultrasonography and account for 15%–20% of all of the prenatally detected congenital anomalies .1 Their severity is extremely variable, but CAKUT represent the leading cause of CKD in children.2 Genetic and environmental factors are proposed to be involved in CAKUT. Except for posterior urethral valves that usually appear as sporadic, familial clustering of CAKUT has been frequently described3; however, the heritability estimates remain unknown. In some families, CAKUT are inherited as a Mendelian monogenic disease, mostly transmitted as an autosomal dominant trait with variable expressivity. Although in these cases CAKUT are frequently associated with extrarenal anomalies (>100 syndromes including renal or urinary abnormalities, with a Mendelian transmission, have been described),4 they can also present with isolated kidney diseases, particularly at prenatal ultrasound screening. Identification of the molecular defect in these patients is important to provide genetic counseling to the families. A monogenic cause of CAKUT was reported to be identified in approximately 10%–15% of patients presenting with isolated CAKUT.5 Among the approximately 50 genes with reported mutations, PAX2 and HNF1B are the two genes most frequently mutated in syndromic or isolated CAKUT.6,7 Copy number variations (CNVs) have also been shown to be frequently associated with syndromic as well as isolated CAKUT; however, these CNVs are frequently attributable to known genomic disorders.8–10 Recently, the simultaneous analysis of multiple candidate genes in parallel by next generation sequencing (NGS) was reported not to lead to a substantial increase in the proportion of patients with CAKUT with identified mutations,11 suggesting that inheritance may frequently be more complex and/or involve noncoding regions, epigenetic modifications, or somatic mutations.
Here, we report targeted NGS of 330 genes known to be either involved in CAKUT or candidate genes in a large cohort of 204 unrelated patients with CAKUT, which led to the identification of a novel CAKUT gene, PBX1.
Results
Patient Cohort
Two hundred four unrelated patients with CAKUT were included (K1–K204 in Supplemental Table 1). Criteria for inclusion were CAKUT involving both kidneys or unilateral CAKUT associated with either extrarenal defects or familial history of CAKUT (Table 1). Patients with posterior urethral valves were not included in the study. Fifty of the 204 patients had been previously tested for mutations in HNF1B (involved in renal cysts and diabetes syndrome) and/or PAX2 (involved in papillorenal syndrome) and/or EYA1 (involved in branchio-oto-renal syndrome) by Sanger sequencing and were revealed to be negative, and nine fetuses had been tested for mutations in RET (involved in bilateral kidney agenesis12), with one identified as carrying a variant of unknown significance.13 In addition, we also tested 11 patients suspected to suffer from branchio-oto-renal syndrome because of branchial and/or ear abnormalities but without any renal phenotype (normal renal ultrasound) (BO1–BO11 in Supplemental Table 1). Sex ratio was 116 boys, 89 girls, and ten patients for which sex was unknown. Not including the 11 BO patients, 79 patients presented extrarenal anomalies. Fifty-five patients had familial history of renal disease affecting first-degree relative(s) (parents, children, or affected siblings), and eight additional patients had history of renal disease in more distant family. Moreover, seven patients had a parent with branchial defects and/or deafness (five), coloboma (one), or diabetes (one) without renal defect (Supplemental Table 1, Table 1). Ninety-three (45%) cases were fetuses affected with extremely severe renal anomalies at fetal ultrasound screening, which led to termination of pregnancy after parents’ request and case by case evaluation by a multidisciplinary prenatal committee in accordance with the French law. Autopsy was performed after informed consent in all cases.
Renal and extrarenal phenotypes of the cohort encompassing 204 patients with CAKUT and 11 patients with branchial and/or ear defect (215 patients)
Screening the Candidate Genes and Identification of PBX1 as a Novel Gene Involved in Monogenic CAKUT
Among the 330 genes, 275 were candidate genes with no mutation reported in patients with CAKUT so far (Supplemental Table 2). Analysis of these genes in the 204 individuals with renal involvement led to the identification of 125 heterozygous variants never reported in the Exome Aggregation Consortium (ExAC) database in a total of 88 genes in 92 patients (Supplemental Tables 1 and 3). In one of these genes, PBX1, we identified five heterozygous loss of function mutations/deletions (Figure 1, Table 2). The first mutation (family 1) was a frameshift mutation (p.Asn143Thrfs*37) in a 21-year-old adult affected with renal hypoplasia and deafness. She has an eGFR of 40 ml/min per 1.73 m2. The second one (family 2) was a nonsense mutation (p.Arg184*) in a child presenting with small hyperechogenic kidneys with cysts, developmental delay, growth retardation, and long and narrow face. She had eGFR of 51 ml/min per 1.73 m2 at the age of 11 years old. The third one (family 3) was a splice mutation (c.511–2A>G) in a fetus with renal hypoplasia and oligoamnios. The autopsy showed extremely severe oligonephronia. RNA was not available to study the effect of that mutation on splicing, but the mutation is predicted to lead to total loss of the 3′ splice site and skipping of exon 4, resulting in a frameshift. The two last patients presented with a heterozygous deletion removing the whole PBX1 gene indicated by a reduced number of reads for the nine exons of the gene (Figure 1). The first one (family 4) was a 39-year-old adult woman followed for small and dysplastic horseshoe kidney also affected with profound deafness. She was the oldest patient of our series and had stage 3B renal failure (eGFR=40 ml/min per 1.73 m2) at last evaluation. The second one (family 5) was an infant with a single small hyperechogenic kidney also presenting developmental delay, microcephaly, and facial dysmorphism (with long and narrow face as well as abnormal ear lobes). He had normal renal function at the age of 18 months old. Comparative Genomic Hybrization (CGH) array analysis (Agilent 60 K) performed on DNA from patients’ lymphocytes confirmed presence of a 2.46-Mb deletion at 1q23.3q24.1 including PBX1 and seven other genes in a family 4′ patient and a 9.1-Mb deletion on the long arm of a chromosome 1 including PBX1 together with 130 other genes (Supplemental Figure 1) in a family 5′ patient. De novo occurrence of the deletion in patients K136 (family 5) was confirmed by fluorescence in situ hybridization on parents’ lymphocytes, which was normal (Supplemental Figure 2). The pelviureteric junction obstruction in the patient’s mother (shown in gray in Figure 1) is thus an incidental association. Sanger sequencing (families 1–3) or analysis of four microsatellite markers located in PBX1 (family 4) allowed us to show that the mutations/deletion in these families also occurred de novo (Figure 1, Supplemental Table 4). In all five families, misclassification of the parents was excluded by analysis of at least 12 unlinked microsatellite markers located on various chromosomes (Supplemental Table 4) and the identification of a paternally inherited CNV on 1p (family 5) (Supplemental Figure 2). Moreover, by retrospective analysis of our cohort of patients with CAKUT, we identified another individual not included in the 204 patients screened by targeted NGS, in whom a de novo 6.2-Mb deletion, including PBX1, had been identified by CGH array analysis. This patient presented with small hyperechogenic kidneys, developmental delay, and facial dysmorphism, similar to the patient from family 5. He already had severe proteinuria associated with a rapidly progressing renal failure (eGFR=60 ml/min per 1.73 m2) at the age of 10 years old.
Identification of de novo PBX1 mutations in three families and deletions in two families. (A) Pedigrees of three families with one affected individual (in black) carrying a de novo point mutations in PBX1. The affected individual in family 1 (K175) was a 21-year-old woman presenting with renal hypoplasia and deafness. The affected individual in family 2 (K179) was a 12-year-old girl presenting with small hyperechogenic kidneys with cysts, developmental delay, growth retardation, and long and narrow face. The affected individual in family 3 (K186) was a male fetus with extremely severe renal hypoplasia. Loss of function mutations (deletion of one base leading to a frameshift, nonsense mutation, or nucleotide change in the consensus acceptor splice site located in 3′ of intron 3) identified in patients K175, K179, and K186, respectively, were validated by Sanger sequencing. Absence of the mutation in the parents showed de novo occurrence in all patients. NT, not tested. (B) Pedigree of family 4 with a 40-year-old woman (K181) affected with small and dysplastic horseshoe kidney and deafness. NGS analysis of patient DNA revealed a reduced number of reads for all of the exons of PBX1 (in red) compared with the mean number of reads for the 37 other DNAs analyzed in the same run (in green), suggesting deletion of one PBX1 allele. Haplotypes of the affected woman and her parents were generated using four known microsatellite markers upstream and downstream of PBX1 (markers A, B, G, and H) as well as four intragenic PBX1 markers (markers C–F). Genomic positions of the microsatellite markers (on Human Genome version GRCh38.p7 from Ensembl) are A (D1S2675): chr1:162,240,203–162,240,364; B (D1S2844): chr1:162,979,036–162,979,218; C (in PBX1 intron 1): chr1:164,562,253–164,562,694; D (in PBX1 intron 2): chr1:164,651,731–164,652,170; E (in PBX1 intron 2): chr1:164,707,219–164,707,656; F (in PBX1 intron 8): chr1:164,833,830–164,834,270; G (D1S2762) chr1:166,986,900–166,987,137; and H (D1S196): chr1:167,635,063–167,635,195. Deletion of a <4-Mb region encompassing PBX1 located on the maternally inherited chromosome was shown by the lack of maternal contribution for markers C–F (genotypes are in Supplemental Table 3). (C) Pedigree of family 5 with a 2-year-old infant boy (K136) presenting with a single hyperechogenic kidney associated with developmental delay, microcephaly, and facial dysmorphism (with long and narrow face as well as abnormal ear lobes). The mother was presenting pyeloureteric junction obstruction. NGS analysis of K136 DNA revealed a reduced number of reads for all of the exons of PBX1 (in red) compared with the mean number of reads for the 38 other DNA analyzed in the same run (in green), suggesting deletion of one PBX1 allele. Validation of this deletion and identification of its extent and de novo occurrence were performed by CGH analysis (Supplemental Figure 1).
De novo loss of function variants in PBX1 (NM_002585)
To precisely evaluate the strength of association between PBX1 mutations and CAKUT, we used a binomial probability test. The P value for identifying three de novo loss of function mutations in that gene among 204 patients was <0.001 (Supplemental Figure 3). The identification of three de novo loss of function mutations together with two de novo large deletions encompassing the gene among 204 patients with CAKUT strongly support the causal effect of the PBX1 mutations.
Screening the Known CAKUT Genes
In 31 of the 204 patients, we identified 24 pathogenic mutations in seven known genes associated with highly penetrant CAKUT, namely HNF1B, PAX2, EYA1, ANOS1, GATA3, CHD7, and KIF14 (Figure 2, Table 3). Three of the mutations/deletions, affecting GATA3 or KIF14, were identified in patients who had been previously screened for HNF1B, PAX2, and/or EYA1 by Sanger sequencing and remained unsolved. Renal phenotypes, extrarenal anomalies, and mutation inheritance are shown in Table 3 (Supplemental Table 1). As frequently observed for patients with mutations in these genes, 16 of the patients had extrarenal symptoms. Notably, the fetuses with biallelic KIF14 mutations presented severe microcephaly and agenesis of the corpus callosum as recently reported.14 Fourteen patients were sporadic, and de novo occurrence of the mutation was shown for the five families with available parental DNA. Somatic mosaicism in the mother was identified in one patient with a PAX2 mutation. Ten patients had a family history of CAKUT (Supplemental Table 1, Table 3), and presence of the mutation was shown in the affected related individual in the seven tested families. In addition, in six patients with a family history of branchial and/or ear defect (three patients), deafness (one patient), diabetes (one patient), or coloboma (one patient), mutations were inherited from a parent with branchial and/or ear defect (EYA1), diabetes (HNF1B), coloboma (PAX2), or deafness (GATA3) but no renal defect in agreement with incomplete penetrance of the renal defect associated with mutations in these genes.15–17 Altogether, 11 mutations were large deletions removing all exons of HNF1B (eight), GATA3 (two), or PAX2 (one); five were intragenic deletions removing few exons of PAX2 (two), EYA1 (one), ANOS1 (one), or KIF14 (one; homozygous); and 18 were point mutations (including small deletions up to 9 bp), all absent from the ExAC, except for one EYA1 mutation that was present with a frequency of one of 66,722 in the European population. The 15 heterozygous large deletions were confirmed by either CGH or multiplex ligation–dependent probe amplification, and the homozygous deletion of KIF14 was confirmed by sequencing of a cDNA synthesized from kidney RNA of the affected fetus (Supplemental Figure 4).
Identification of causative mutations/deletions in 36 of the 204 CAKUT cases. For each gene (ANOS1, CHD7, EYA1, GATA3, HNF1B, KIF14, PAX2, and PBX1), the proportions of patients with mutations and deletions are shown in dark and light colors, respectively. All of the mutations/deletions are heterozygous, except for KIF4, for which they are biallelic.
Causative mutations in known genes associated with highly penetrant CAKUT phenotype
We also identified 21 rare heterozygous damaging variants in 12 genes reported as involved in dominant forms of CAKUT in a total of 17 individuals with CAKUT (Table 4). Extrarenal symptoms usually associated with mutations in some of these genes were never observed in these patients. All of the variants were missense predicted to be damaging by at least three of the five prediction programs used, and 13 of them, affecting a total of ten genes (BICC1, CDC5L, CHD1L, NOTCH2, RET, SALL1, SALL4, TBC1D1, TBX18, and TNXB), were absent from the ExAC. This suggests a causative effect for these variants, which remains to be shown. Of note, two of these variants were present in individuals with a pathogenic mutation identified (one variant in CDC5L in a child with a causative mutation in CHD7 and one in TBX18 in a fetus with a causative de novo mutation in GATA3), thus questioning the effect of these additional variants. Other patients were carrying variants in several genes, which may suggest oligogenic inheritance but also, highlights the arduousness of the interpretation of the variants (Table 4).
Variants in other genes reported in dominant forms of CAKUT, including previously reported variants and novel rare variants
In addition, using less stringent criteria for variant filtering, we identified four previously reported mutations with frequency in the ExAC that was above the one to 5000 threshold (Table 4). The RET p.Asp567Gln variation in a fetus with bilateral kidney agenesis (that patient has already been reported13) was inherited from his unaffected father and not carried by his sister, who presented with a less severe renal phenotype. This variant was identified in 19 individuals without CAKUT (seven being unrelated) from our in-house exome database (>8000 exomes). The RET p.Arg694Gln variation detected in a fetus with bilateral hypoplasia has been reported in a patient with Hirschsprung disease without associated renal involvement.18 This mutation has low pathogenicity scores. The SOX17 p.Tyr259Asn variation, previously reported to be associated with vesicoureteral reflux,19 was identified in a patient also carrying missense variants in BICC1 and DSTYK. That SOX17 variant was inherited from his unaffected father. It is found with a frequency of one in 181 in the European population in the ExAC and in 82 individuals without CAKUT in our in-house exome database. Finally, the DSTYK splice variant (c.654+1), recently reported in a large family with CAKUT,20 was identified in a BO individual with normal renal ultrasound scan. That variant is found in the ExAC database with a frequency of one in 3563 in the European population and was identified in 11 individuals without CAKUT (ten being unrelated) in our in-house exome database.
Burden Analyses of Rare Variants in Patients with Unsolved Cases
Finally, for the 168 patients with no mutation identified in PBX1 or any of the highly penetrant CAKUT genes (HNF1B, PAX2, EYA1, ANOS1, GATA3, CHD7, or KIF14), we performed a burden test to look for genes with an enrichment in rare variants. When considering only loss of function variants, we did not identify any gene with nominal P value <0.05. When considering loss of function plus missense, we identified a few genes with nominal P value <0.05, but none of them reached significance with Bonferroni correction (P<1.5×10−4). These results are shown in Supplemental Table 5.
Discussion
Whole-exome sequencing is currently considered as the optimal approach for identifying new genes involved in human diseases. However, targeted exome sequencing, because it offers the advantage to analyze larger series of patients for the same cost, may prove more efficient for mutation identification in the context of highly genetically heterogeneous diseases, such as CAKUT, provided that the gene selection is relevant. Analysis of candidate genes allowed us to identify five de novo PBX1 mutations/deletions, including three loss of function point mutations, showing that this gene is a novel CAKUT gene. PBX1 encodes the Pre-B Cell Leukemia Transcription Factor 1, a TALE homeodomain transcription factor known to play a crucial role in several developmental processes and modify Hox gene activity.21 In the mouse, Pbx1 was notably shown to play a role during kidney development.22,23 It is expressed in the nephrogenic mesenchyme, weakly in the induced mesenchyme and strongly in the stroma after ureteric bud invasion. Its expression decreases when cells undergo epithelial differentiation, whereas it increases in glomerular cells at the stage of capillary loop extension. Except for Pdgfrb, which was recently reported as a direct target gene during vascular patterning,24 Pbx1 transcriptional targets in the kidney remain largely unknown. Although no defect was reported in Pbx1+/− kidneys, a majority of Pbx1−/− embryos exhibit hypoplastic kidneys (30% with unilateral kidney agenesis), poorly defined cortical and medullary regions (E13.0), and a reduced number of differentiating nephrons (E14.5). Mutant kidneys were also reported to be mispositioned caudally. Histologic analysis revealed defects in ureteric branching and abnormal expansion of induced mesenchyme, arising from a mesenchymal dysfunction that led to a reduced mesenchymal to epithelial differentiation.23 Pbx1−/− mouse kidney phenotypes are thus highly reminiscent of those of the patients with heterozygous PBX1 mutations (hypoplastic kidneys [all patients], unilateral agenesis [K136], oligonephronia [K186], and horseshoe kidney with absence of corticomedullar differentiation [K181]). A phenotype associated with haploinsufficiency in human but not in mouse has already been reported for other transcription factor–encoding genes involved in kidney development (e.g., HNF1B, SIX1, and SIX2).
In 2012, Sanna-Cherchi et al.9 reported a 0.5-Mb de novo deletion at 1q32, including PBX1, in a patient with renal hypoplasia. More recently, the study of overlapping deletions on the long arm of chromosome 1 in several patients presenting with renal malformation associated with other defects suggested that haploinsufficiency of PBX1 might be responsible for the renal phenotype.25 Our data confirm that hypothesis, because three PBX1 mutations were de novo point mutations leading to a null allele. As for CAKUT associated with heterozygous mutations in genes encoding other transcription factors (HNF1B, PAX2, EYA1, and SALL1), the severity of the renal disease associated with PBX1 defect seems variable: the oldest patient presents with stage 3 renal failure (eGFR of 40 ml/min per 1.73 m2) at the age of 39 years old, whereas the fetal case showed major renal hypoplasia leading to oligohydramnios. This could be due to variants in modifier genes. Associated deafness was present in two of the five patients with PBX1 mutation/deletion who were not carrying additional mutation or deletion of PAX2, GATA3, EYA1, or SIX1, four genes with mutations that may lead to deafness. In a mouse cochlear hair cell line, Pbx1 is coexpressed with Gata3,26 suggesting that PBX1 mutations might be responsible for deafness. However, the majority of patients with large deletions at 1q23.3q24.1 removing PBX1 do not have deafness.25 Microcephaly, facial dysmorphism, ear anomalies, and developmental delay in a family 5 patient were reported in other patients with large deletions at that locus25; however, developmental delay and facial dysmorphism were also present in the patient from family 2 carrying a PBX1 nonsense mutation. Screening of larger series of patients with CAKUT for mutations in PBX1 will be necessary to more precisely evaluate their association with extrarenal phenotypes.
Our strategy of targeted exome sequencing also allowed us to improve the efficiency of identification of causative mutations in patients with known CAKUT compared with our previous workflow (by Sanger sequencing and MLPA analysis), in which only HNF1B, PAX2, and/or EYA1 were screened. Altogether, causative mutations, including the five PBX1 mutations, were identified in 36 of the 204 (17.6%) patients with CAKUT. This rate is far above that reported recently by a similar NGS approach (six of 453),11 and this is likely due, at least in part, to the fact that the spectrum of phenotypes included was different. We only included CAKUT affecting both kidneys and/or familial and/or syndromic forms, excluding posterior urethral valves known to be most frequently sporadic, whereas Nicolaou et al.11 included unilateral forms and many posterior urethral valves. The rate of mutation that we observed here is also higher than that reported by Thomas et al.6 but close to that reported by Weber et al.5 or Madariaga et al.7 Of note, the proportion of patients with HNF1B, PAX2, or EYA1 mutations is probably slightly biased downward, because our cohort included 50 patients already known not to carry mutations in these genes. Interestingly, 40% of the identified mutations are deletions, removing either a few exons or the whole gene. Indeed, a major advantage of NGS is the identification of both point mutations and deletions (intragenic and whole-gene deletions) in a sole experiment. In addition to HNF1B deletions, heterozygous deletions, which would not have been diagnosed by our previous workflow, were identified in PAX2, EYA1, GATA3, and ANOS1. Also, one homozygous intragenic deletion of KIF14 was shown in a consanguineous family. No mutation was identified in ITGA8 or FGF20, two genes recently reported as responsible for autosomal recessive bilateral kidney agenesis.27,28
Among variants identified in this study, some are in other genes previously reported as responsible for dominant CAKUT (Table 4). However, although these variants met criteria supporting their pathogenicity (absence or low frequency in the ExAC and pathogenicity scores indicating deleterious effect), in the seven patients whose parents could be tested, the variants were inherited from an unaffected parent (Table 4), thus questioning their causality. In one patient with renal cystic hypodysplasia and deafness, compound heterozygous variants in NOTCH2 were shown in the index patient, the significance of which is unknown. Growing numbers of available data regarding the frequency of variants among large control populations can lead variants initially considered as pathogenic to be reclassified. The frequency of the SOX17 p.Tyr259Asn variation in the ExAC and our in-house database indicates that this variant is likely not causative, although it could be an at-risk allele as recently proposed.29 As for the DSTYK splice variant (c.654+1)20 that we identified in a BO individual without any renal involvement, the rate of this variant in the ExAC and our in-house database supports that, if involved in the development of a CAKUT phenotype, it has an incomplete penetrance. This highlights the fact that one must be cautious when interpreting variants reported as pathogenic in the literature. Indeed, some of them can reveal as low-penetrance alleles, modifiers, or even frequent polymorphisms instead of causative mutations.
Except for KIF14, no ballelic mutation was identified in any of the genes previously reported in recessive forms of CAKUT. As in the work by Nicolaou et al.,11 we did not identify any significant increase in rare variants in any of the tested genes, reflecting the complexity of genetic studies in CAKUT. Moreover, no pathogenic mutation or rare variant of unknown significance was identified in 40% of our series (Supplemental Table 1). Although mutations or CNVs in gene(s), noncoding regions, or microRNAs not targeted by our capture library could be involved in some of these patients, somatic events and/or environmental factors or epigenetic mechanisms likely explain at least part of this large fraction of patients. The familial aggregation of CAKUT as well as the wide spectrum of severity of the phenotypes suggest a complex genetic architecture as observed in neurodevelopmental disorders.30 This would fit with a model in which various combinations of variants as well as environmental factors will cause early kidney development insults and act together to alter renal and urinary tract formation. The timing, severity, location, and extent of the “deviation from an optimal renal developmental program” may determine the ultimate phenotype. Collaborative effort will be essential in the future to test a larger number of individuals by whole-exome sequencing, classify variants, and establish high-throughput functional assays to successfully decipher the complex genetic mechanisms of CAKUT and be able to propose appropriate genetic tests and counseling for the families.
During the review process of the paper, Le Tanno et al.31 reported a novel series of de novo 1q23.3-q24.1 microdeletions in patients with syndromic CAKUT. The patient carrying the smallest deletion, limited to PBX1, presented with mild developmental delay and hearing loss, which was observed in two patients of our series harboring PBX1 point mutations.
Concise Methods
Design of the Targeted SureSelect Library
The 330 genes included in our panel are shown in Supplemental Table 2 and include 55 genes with mutations reported in the literature in patients with isolated and/or syndromic CAKUT and 275 candidate genes (including 104 genes with knockout in mouse that led to kidney developmental defects, 84 genes involved in cellular processes/signaling pathways relevant for kidney development, 11 genes with a role in ureter/bladder development, 57 genes encoding reported targets of transcription factors WT1 or HNF1B, and 19 homologs of genes in the sublists above expressed during kidney/lower urinary tract development). The custom SureSelect gene panel was designed using the SureDesign software (Agilent). The target regions covered 1.38 Mb, including coding exons and splice junctions.
Targeted Exome Sequencing and Prioritization of the Variants
DNA was extracted from blood cells (living patients) or frozen tissues (fetuses). The study was approved by the Comité de Protection des Personnes pour la Recherche Biomédicale Ile de France 2, and informed consent was obtained from the patients or parents. Illumina-compatible precapture barcoded genomic DNA libraries were constructed, and a series of 16 or 36 barcoded libraries was pooled at equimolar concentrations. The capture process was performed according to the SureSelect protocol (Agilent) and sequencing on an Illumina HiSeq2500. Sequences were aligned to the reference human genome hg19 using the Burrows–Wheeler Aligner.
The mean depth of coverage ranged from 215 to 719, with ≥98% of the bases covered at least 30× (Supplemental Material). Variants were prioritized according to their frequency in the ExAC database32 and their predicted damaging effect. For known CAKUT genes, we selected previously published variants as well as nonsense, frameshift, splice variants, and missense predicted as damaging by at least three of five used prediction programs (Polyphen2, Sift, MutationTaster, Grantham score, and Align-GVGD) (Supplemental Material) with a minor allele frequency ≤1/5000 for genes involved in dominant forms and ≤1/1000 for those involved in recessive forms. For candidate genes, only variants that were absent from the ExAC were retained.
To evaluate duplication and large deletion events, for each individual, the relative read count for each targeted region was determined as the ratio of the read count for that region divided by the total absolute read counts of all targeted regions of the design. The ratio of the relative read count of a region in a given individual over the average relative read counts in other individuals of the run resulted in the estimated copy number for that region in that individual (method adapted from ref. 33).
Burden Test
To look for an excess of rare variants in the 330 genes in patients with CAKUT, we conducted a burden test using the Madsen and Browning34 Wilcoxon rank sum test and the SKAT-o test35 as implemented in EPACTS software (http://genome.sph.umich.edu/wiki/EPACTS). For the patients, we included the 168 patients with no mutation identified in PBX1 or any of the highly penetrant CAKUT genes (HNF1B, PAX2, EYA1, ANOS1, GATA3, CHD7, or KIF14). For the controls, we used 426 unrelated individuals who were unaffected parents of children presenting with immune deficiency that were sequenced by whole-exome sequencing on the same genomic platform. Variants in patients and controls were called together. Variants were filtered according to their frequency in public databases (the ExAC, dbSNP, 1000G, and EVS) using two different thresholds (<0.01% and <0.001%), and we then considered either loss of function variants (nonsense, frameshift, splice, start loss, and stop loss) or loss of function plus missense predicted as damaging by Polyphen2 (score >0.470) and Sift (score <0.05). Large deletions were not taken into account in this analysis.
Disclosures
None.
Acknowledgments
We thank patients and families for their participation. We also thank the foetopathologists of the SOciété Française de FOETopatholgie for providing us with fetal tissue samples and the pediatric nephrologists from the French Society of Pediatric Nephrology for referring the patients. We thank Marie-Claire Gubler for her help with histology analysis of fetal kidneys. We thank Olivier Gribouval for help with microsatellite analysis.
The research leading to these results has received funding from the European Union’s Seventh Framework Programme (FP7/2007–2013) grant 305608 (EURenOmics). C. Humbert was funded by EURenOmics. This work was also supported by Programme Hospitalier de la Recherche Clinique Assistance Publique grant AOM07129 (to R.S.), Investments for the Future Program grant ANR-10-IAHY-01 (to S.S. and C.A.), and Fondation pour la Recherche Médicale grant DEQ20130326532 (to S.S.).
Footnotes
Published online ahead of print. Publication date available at www.jasn.org.
This article contains supplemental material online at http://jasn.asnjournals.org/lookup/suppl/doi:10.1681/ASN.2017010043/-/DCSupplemental.
- Copyright © 2017 by the American Society of Nephrology