Copy Number Variant Analysis and Genome-wide Association Study Identify Loci with Large Effect for Vesicoureteral Reflux

BACKGROUND

Vesicoureteral reflux (VUR) is a common, familial genitourinary disorder, and a major cause of pediatric urinary tract infection (UTI) and kidney failure. The genetic basis of VUR is not well understood.

METHODS

A diagnostic analysis sought rare, pathogenic copy number variant (CNV) disorders among 1737 patients with VUR. A GWAS was performed in 1395 patients and 5366 controls, of European ancestry.

RESULTS

Altogether, 3% of VUR patients harbored an undiagnosed rare CNV disorder, such as the 1q21.1, 16p11.2, 22q11.21, and triple X syndromes ((OR, 3.12; 95% CI, 2.10 to 4.54; P=6.35×10^-8) The GWAS identified three study-wide significant and five suggestive loci with large effects (ORs, 1.41-6.9), containing canonical developmental genes expressed in the developing urinary tract (WDPCP, OTX1, BMP5, VANGL1, and WNT5A). In particular, 3.3% of VUR patients were homozygous for an intronic variant in WDPCP (rs13013890; OR, 3.65; 95% CI, 2.39 to 5.56; P=1.86×10^-9). This locus was associated with multiple genitourinary phenotypes in the UK Biobank and eMERGE studies. Analysis of Wnt5a mutant mice confirmed the role of Wnt5a signaling in bladder and ureteric morphogenesis.

CONCLUSIONS

These data demonstrate the genetic heterogeneity of VUR. Altogether, 6% of patients with VUR harbored a rare CNV or a common variant genotype conferring an OR >3. Identification of these genetic risk factors has multiple implications for clinical care and for analysis of outcomes in VUR.

Investigation of DNA variants specific to ROBO2 Isoform 'a' in Irish vesicoureteric reflux patients reveals marked CpG island variation

ROBO2 gene disruption causes vesicoureteric reflux (VUR) amongst other congenital anomalies. Several VUR patient cohorts have been screened for variants in the ubiquitously expressed transcript, ROBO2b, but, apart from low levels in a few adult tissues, ROBO2a expression is confined to the embryo, and might be more relevant to VUR, a developmental disorder. ROBO2a has an alternative promoter and two alternative exons which replace the first exon of ROBO2b. We screened probands from 251 Irish VUR families for DNA variants in these. The CpG island of ROBO2a, which includes the non-coding first exon, was found to contain a run of six variants abolishing/creating CpG dinucleotides, including a novel variant, present in the VUR cases in one family, that was not present in 592 healthy Irish controls. In three of these positions, the CpG was created by the non-reference allele, and the reference allele was not the nucleotide that would result from spontaneous deamination of methylcytosine to thymine, suggesting that there might have been selection for variability in number of CpGs in this island. This is in marked contrast to the CpG island at the start of ROBO2b, which only contained a single variant that abolishes a CpG.

Author Correction: The copy number variation landscape of congenital anomalies of the kidney and urinary tract

In the version of this article initially published, affiliation 38 incorrectly read "ICNU-Nephrology and Urology Department, Barcelona, Spain"; "Renal Division, Hospital Clinic, IDIBAPS, University of Barcelona, Barcelona, Spain" is the correct affiliation. The error has been corrected in the HTML and PDF versions of the article.

The copy number variation landscape of congenital anomalies of the kidney and urinary tract

Congenital anomalies of the kidney and urinary tract (CAKUT) are a major cause of pediatric kidney failure. We performed a genome-wide analysis of copy number variants (CNVs) in 2,824 cases and 21,498 controls. Affected individuals carried a significant burden of rare exonic (i.e. affecting coding regions) CNVs and were enriched for known genomic disorders (GD). Kidney anomaly (KA) cases were most enriched for exonic CNVs, encompassing GD-CNVs and novel deletions; obstructive uropathy (OU) had a lower CNV burden and an intermediate prevalence of GD-CNVs; vesicoureteral reflux (VUR) had the fewest GD-CNVs but was enriched for novel exonic CNVs, particularly duplications. Six loci (1q21, 4p16.1-p16.3, 16p11.2, 16p13.11, 17q12, and 22q11.2) accounted for 65% of patients with GD-CNVs. Deletions at 17q12, 4p16.1-p16.3, and 22q11.2 were specific for KA; the 16p11.2 locus showed extensive pleiotropy. Using a multidisciplinary approach, we identified TBX6 as a driver for the CAKUT subphenotypes in the 16p11.2 microdeletion syndrome.

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.).

Urinary tract effects of HPSE2 mutations

Urofacial syndrome (UFS) is an autosomal recessive congenital disease featuring grimacing and incomplete bladder emptying. Mutations of HPSE2, encoding heparanase 2, a heparanase 1 inhibitor, occur in UFS, but knowledge about the HPSE2 mutation spectrum is limited. Here, seven UFS kindreds with HPSE2 mutations are presented, including one with deleted asparagine 254, suggesting a role for this amino acid, which is conserved in vertebrate orthologs. HPSE2 mutations were absent in 23 non-neurogenic neurogenic bladder probands and, of 439 families with nonsyndromic vesicoureteric reflux, only one carried a putative pathogenic HPSE2 variant. Homozygous Hpse2 mutant mouse bladders contained urine more often than did wild-type organs, phenocopying human UFS. Pelvic ganglia neural cell bodies contained heparanase 1, heparanase 2, and leucine-rich repeats and immunoglobulin-like domains-2 (LRIG2), which is mutated in certain UFS families. In conclusion, heparanase 2 is an autonomic neural protein implicated in bladder emptying, but HPSE2 variants are uncommon in urinary diseases resembling UFS.

A new genome scan for primary nonsyndromic vesicoureteric reflux emphasizes high genetic heterogeneity and shows linkage and association with various genes already implicated in urinary tract development

Primary vesicoureteric reflux (VUR), the retrograde flow of urine from the bladder toward the kidneys, results from a developmental anomaly of the vesicoureteric valve mechanism, and is often associated with other urinary tract anomalies. It is the most common urological problem in children, with an estimated prevalence of 1–2%, and is a major cause of hypertension in childhood and of renal failure in childhood or adult life. We present the results of a genetic linkage and association scan using 900,000 markers. Our linkage results show a large number of suggestive linkage peaks, with different results in two groups of families, suggesting that VUR is even more genetically heterogeneous than previously imagined. The only marker achieving P < 0.02 for linkage in both groups of families is 270 kb from EMX2. In three sibships, we found recessive linkage to KHDRBS3, previously reported in a Somali family. In another family we discovered sex-reversal associated with VUR, implicating PRKX, for which there was weak support for dominant linkage in the overall data set. Several other candidate genes are suggested by our linkage or association results, and four of our linkage peaks are within copy-number variants recently found to be associated with renal hypodysplasia. Undoubtedly there are many genes related to VUR. Our study gives support to some loci suggested by earlier studies as well as suggesting new ones, and provides numerous indications for further investigations.

The increased incidence of the RET p.Gly691Ser variant in French-Canadian vesicoureteric reflux patients is not replicated by a larger study in Ireland

The p.Gly691Ser variant of the RET protein, resulting from the 'A' allele of the SNP rs1799939 in exon 11 of the RET gene, was recently found to be present in a high proportion of primary vesicoureteric reflux (pVUR) patients in Quebec. We have determined the genotype of this SNP in 221 unrelated index cases of pVUR from the Irish population, in 190 full siblings of 160 of the index cases, and in 592 healthy controls. We found no significant difference in genotype or allele frequencies in patients and controls, and no tendency of affected siblings to share the same genotype. We also found no difference in the presence of additional phenotypic features such as duplex kidneys, between patients with and without the 'A' allele, and no difference in grade of reflux. We find no evidence of any influence of RET SNP rs1799939 on pVUR phenotype.

A genome-wide scan for genes involved in primary vesicoureteric reflux

BACKGROUND

Vesicoureteric reflux (VUR) is the retrograde flow of urine from the bladder into the ureters. It is the most common urological anomaly in children, and a major cause of end-stage renal failure and hypertension in both children and adults. VUR is seen in approximately 1-2% of Caucasian newborns and is frequently familial.

OBJECTIVE AND METHODS

In order to search for genetic loci involved in VUR, we performed a genome-wide linkage scan using 4710 single-nucleotide polymorphisms (SNPs) in 609 individuals from 129 Irish families with >1 affected member.

RESULTS

Nonparametric linkage (NPL) analysis of the dataset yielded moderately suggestive linkage at chromosome 2q37 (NPL(max) = 2.67, p<0.001). Analysis of a subset without any additional features, such as duplex kidneys, yielded a maximum NPL score of 4.1 (p = 0.001), reaching levels of genome-wide statistical significance. Suggestive linkage was also seen at 10q26 and 6q27, and there were several smaller peaks.

CONCLUSION

Our results confirm the previous conclusion that VUR is genetically heterogeneous, and support the identification of several disease-associated regions indicated by smaller studies, as well as indicating new regions of interest for investigation.

Gene conversion in human rearranged immunoglobulin genes

Over the past 20 years, many DNA sequences have been published suggesting that all or part of the V(H) segment of a rearranged immunoglobulin gene may be replaced in vivo. Two different mechanisms appear to be operating. One of these is very similar to primary V(D)J recombination, involving the RAG proteins acting upon recombination signal sequences, and this has recently been proven to occur. Other sequences, many of which show partial V(H) replacements with no addition of untemplated nucleotides at the V(H)-V(H) joint, have been proposed to occur by an unusual RAG-mediated recombination with the formation of hybrid (coding-to-signal) joints. These appear to occur in cells already undergoing somatic hypermutation in which, some authors are convinced, RAG genes are silenced. We recently proposed that the latter type of V(H) replacement might occur by homologous recombination initiated by the activity of AID (activation-induced cytidine deaminase), which is essential for somatic hypermutation and gene conversion. The latter has been observed in other species, but not in human Ig genes, so far. In this paper, we present a new analysis of sequences published as examples of the second type of rearrangement. This not only shows that AID recognition motifs occur in recombination regions but also that some sequences show replacement of central sections by a sequence from another gene, similar to gene conversion in the immunoglobulin genes of other species. These observations support the proposal that this type of rearrangement is likely to be AID-mediated rather than RAG-mediated and is consistent with gene conversion.

V(H) replacement in rearranged immunoglobulin genes

Examples suggesting that all or part of the V(H) segment of a rearranged V(H)DJ(H) may be replaced by all or part of another V(H) have been appearing since the 1980s. Evidence has been presented of two rather different types of replacement. One of these has gained acceptance and has now been clearly demonstrated to occur. The other, proposed more recently, has not yet gained general acceptance because the same effect can be produced by polymerase chain reaction artefact. We review both types of replacement including a critical examination of evidence for the latter. The first type involves RAG proteins and recombination signal sequences (RSS) and occurs in immature B cells. The second was also thought to be brought about by RAG proteins and RSS. However, it has been reported in hypermutating cells which are not thought to express RAG proteins but in which activation-induced cytidine deaminase (AID) has recently been shown to initiate homologous recombination. Re-examination of the published sequences reveals AID target sites in V(H)-V(H) junction regions and examples that resemble gene conversion.

Non-functional immunoglobulin G transcripts in a case of hyper-immunoglobulin M syndrome similar to type 4

Summary 86% of immunoglobulin G (IgG) heavy-chain gene transcripts were found to be non-functional in the peripheral blood B cells of a patient initially diagnosed with common variable immunodeficiency, who later developed raised IgM, whereas no non-functionally rearranged transcripts were found in the cells of seven healthy control subjects. All the patient's IgM heavy-chain and kappa light-chain transcripts were functional, suggesting that either non-functional rearrangements were being selectively class-switched to IgG, or that receptor editing was rendering genes non-functional after class-switching. The functional gamma-chain sequences showed a normal rate of somatic hypermutation while non-functional sequences contained few somatic mutations, suggesting that most came from cells that had no functional gene and therefore were not receiving signals for hypermutation. However, apoptosis of peripheral blood lymphocytes was not impaired. No defects have been found in any of the genes currently known to be responsible for hyper-IgM syndrome but the phenotype fits best to type 4.

Secondary structures in d(CGG) and d(CCG) repeat tracts

Several studies have been made to elucidate the nature of secondary structures in the single strands of d(CGG).d(CCG) repeat tracts but with conflicting conclusions. Here, we review this work and attempt to come towards consensus. Some investigators find that the G-rich strand forms hairpins. Of these, some conclude that pairing is in the alignment d(GGC).d(GGC) with two Watson-Crick bonds and one G.G bond per duplex repeat, others conclude that the alignment is d(GCG).d(GCG) with two G.G bonds and one C.C bond per duplex repeat. Others find quadruplex formation and conclude that this is in the latter alignment with two G4-quartets per quadruplex repeat and C.C bonds. We investigate why these different results were obtained and conclude that quadruplexes are likely to form under physiological conditions. We argue that they are probably bonded in the alignment d(GGC).d(GGC) with one G4-quartet and two C.G.C.G. quartets per quadruplex repeat. The C-rich strand does not appear to form quadruplexes under physiological conditions but forms hairpins. Apparently, short hairpins adopt the alignment d(CCG).d(CCG) with mismatched cytosine residues stacked into the helix but with 15 or more repeat units, the dominant form is a distorted hairpin aligned as d(GCC).d(GCC) with unpaired cytosine residues possibly turned outwards and stacked in the minor groove.

Evidence for two preferred hairpin folding patterns in d(CGG).d(CCG) repeat tracts in vivo

Unusual DNA secondary structures have been implicated in the expansion of trinucleotide repeat tracts that has been found to be responsible for a growing number of human inherited disorders and folate-sensitive fragile chromosome sites. By inserting trinucleotide repeat sequences into a palindromic clamp in lambda phage we are able to investigate their tendencies to form hairpins in vivo in any particular alignment and with odd or even numbers of repeat units in the hairpin. We previously showed that with d(CAG).d(CTG) repeat tracts there was a markedly greater tendency to form hairpins with even numbers of repeat units than with odd numbers, whereas d(GAC).d(GTC) repeats showed no such alternation despite having the same base composition. We expected that d(CGG).d(CCG) repeats, might show the same pattern as d(CAG).(CTG) repeats since they are also involved in trinucleotide repeat expansion disorders. The pattern was not so clear and we wondered whether this might be because d(CGG).d(CCG) repeats have more than one possible alignment in which they could self-anneal. We now present results for all three alignments, which suggest that while even-membered hairpins are preferred in the frame d(CGG).d(CCG), hairpins with odd numbers of trinucleotides are more stable in the frame d(GGC).d(GCC). In both cases the base-pair predicted to close the terminal loop of unpaired bases is 5'C.3'G which has previously been found to be a favoured loop-closing pair.

The effects of trinucleotide repeats found in human inherited disorders on palindrome inviability in Escherichia coli suggest hairpin folding preferences in vivo

Unusual DNA secondary structures have been implicated in the expansion of trinucleotide repeat tracts that are associated with several human inherited disorders. We present evidence consistent with the folding of these trinucleotide repeats into hairpin loops at the center of a long DNA palindrome in vivo. Our assay utilizes a palindrome in bacteriophage lambda, the center of which determines its ability to inhibit plaque formation in a manner that is consistent with folding into a hairpin or cruciform structure. We show that central inserts of even numbers of d(CAG).d(CTG) repeats inhibit plaque formation more than do odd numbers. Both d(CAG)2.d(CTG)2 and d(CGG)2.d(CCG)2 central sequences behave like DNA sequences known to form two-base loops in vitro, suggesting that they may also form compact and stable loops. By contrast, repeats of d(GAC).d(GTC) do not show any evidence consistent with unusual loop stability. These results agree with in vitro evidence that the unstable repeats can form hairpin secondary structures and suggest a favored position of folding. We discuss the potential roles of secondary structures, DNA replication and recombination in models of repeat tract expansion.