A Common RET Variant Is Associated with Reduced Newborn Kidney Size and Function
Zhao Zhang*,
Jackie Quinlan*,
Wendy Hoy,
Michael D. Hughson,
Mathieu Lemire,
Thomas Hudson,
Pierre-Alain Hueber*,
Alice Benjamin||,
Anne Roy¶,
Elena Pascuet*,
Meigan Goodyer*,
Chandhana Raju*,
Fiona Houghton*,
John Bertram** and
Paul Goodyer*
* McGill University Health Centre Research Institute, McGill University and Génome Québec Innovation Centre, and Departments of || Obstetrics and Gynecology and ¶ Radiology, McGill University, Montreal, Québec, Canada; Centre for Chronic Diseases, University of Queensland, Brisbane, and ** Department of Anatomy and Developmental Biology, Monash University, Victoria, Australia; and University of Mississippi Medical Center, Jackson, Mississippi
Correspondence: Dr. Paul Goodyer, 4060 Saint Catherine West, PT-413, Montreal Quebec, Canada H3Z 2Z3. Phone: 514-412-4461; Fax: 514-412-4478; E-mail: paul.goodyer{at}mcgill.ca
Received for publication October 12, 2007.
Accepted for publication April 18, 2008.
Congenital nephron number varies five-fold among normal humans,and individuals at the lower end of this range may have an increasedlifetime risk for essential hypertension or renal insufficiency;however, the mechanisms that determine nephron number are unknown.This study tested the hypothesis that common hypomorphic variantsof the RET gene, which encodes a tyrosine kinase receptor criticalfor renal branching morphogenesis, might account for subtlerenal hypoplasia in some normal newborns. A common single-nucleotidepolymorphism (rs1800860 G/A) was identified within an exonicsplicing enhancer in exon 7. The adenosine variant at mRNA position1476 reduced affinity for spliceosome proteins, enhanced thelikelihood of aberrant mRNA splicing, and diminished the levelof functional transcript in human cells. In vivo, normal whitenewborns with an rs1800860(1476A) allele had kidney volumes10% smaller and cord blood cystatin C levels 9% higher thanthose with the rs1800860(1476G) allele. These findings suggestthat the RET(1476A) allele, in combination with other commonpolymorphic developmental genes, may account for subtle renalhypoplasia in a significant proportion of the white population.Whether this gene variant affects clinical outcomes requiresfurther study.
Human kidney development begins at approximately 5 wk gestation,when RET-bearing cells of the descending nephric duct encounterglial-derived neurotrophic factor (GDNF) released by metanephricmesenchyme at somite 24.1 In response to the trophic effectsof GDNF, a ureteric bud sprouts from each nephric duct and arborizeswithin the lateral mesenchyme. Signals from each ureteric budbranch tip induce adjacent metanephric stem cells to form individualnephrons, which fuse to the tree-like collecting system. Becausenew nephron formation ends at approximately 36 wk gestation,the extent of branching nephrogenesis by this time determinesnephron endowment for life.
Interestingly, nephron number varies widely (0.3 to 1.3 millionnephrons per kidney) among normal humans.2 Although once dismissedas a benign reflection of human diversity, Brenner et al.3 proposedthat individuals born at the low end of the nephron endowmentspectrum may have increased risk for developing "essential"hypertension and renal insufficiency later in life. They hypothesizedthat signals driving compensatory hypertrophy of overworkednephrons cause glomerulosclerosis and a cycle of subtle, slowlyprogressive renal dysfunction.4
Recent evidence supports this theory; an autopsy study by Kelleret al.5 showed that German adults with essential hypertensionhad 47% fewer nephrons per kidney than well-matched normotensivecontrol subjects. As predicted, hypertensive patients had hypertrophicglomeruli (glomerular volume 233% of control subjects) and increasedglomerulosclerosis (5.5% of glomeruli versus 0% in control subjects).Other evidence suggests that racial differences in congenitalnephron number might also explain the relatively high incidenceof ESRD in Aboriginal versus white populations.6,7 At autopsy,Aboriginal individuals have 23% fewer glomeruli (683,174 perkidney) than white individuals (885,318 per kidney; P < 0.04).
Little is known about the factors that set nephron number, butthe GDNF/RET signaling pathway seems to play a central role.The RET gene (NM_020630) encodes a 1072 amino acid transmembranetyrosine kinase receptor (NP_065681) expressed at the tips ofthe branching ureteric buds during fetal kidney development8;homozygous Ret knockout mice are anephric,9 and heterozygoteshave a 22% reduction in nephron number at 15 d of age.10 Inthis study, we hypothesized that hypomorphic variants of theRET gene might be prevalent in "normal" humans, contributingto suboptimal nephron number and subtle renal hypoplasia ina significant portion of the population.
We scanned the NCBI dbSNP database (see Appendix) and identifiedtwo common variants (minor allele frequency >10%) of theRET gene coding region that might alter its tyrosine kinasereceptor function during kidney development. In exon 11, anA/G substitution at mRNA position 2251 (rs1799939) occurs inapproximately 16% of white alleles, changing gly691 to ser691.In exon 7, an A/G substitution at mRNA position 1476 (rs1800860)occurs in 25% of white alleles and lies within an exonic spliceenhancer (ESE) sequence (Figure 1A). The critical position ofthis substitution suggests that it might modify pre-mRNA splicing.We used ESEfinder11 to calculate the effect of the rs1800860polymorphism on binding to proteins (SF2/ASF, SC35) of the mRNAspliceosome (Table 1). In each case, the 1476(A) allele reducesRET affinity score below the predicted threshold for effectivebinding to these components of the splicing machinery.
Figure 1. 1476A SNP in human RET mRNA is associated with an aberrant 190-bp transcript. (A) Genomic structure of the human RET gene, indicating the 1476(G/A) SNP within an ESE of exon 7. (B) Nested primers spanning exon 6 to exon 17 were used to amplify RET transcripts in reverse-transcribed RNA extracted from various human cell lines. Only the predicted wild-type 1252-bp transcript is evident in the homozygous (1476 G/G) Wit49 Wilms tumor cell line (lane 3); however, an aberrant 190-bp transcript is noted in heterozygous (1476G/A) SH-SY5Y neuroblastoma cells (lane 2) and in homozygous 1476(A/A) SK-N-BE(2) neuroblastoma cells (lane 1). The aberrant 190-bp transcript was not seen in four other 1476G/G Wilms tumor, renal cell carcinoma, and ovarian carcinoma cell lines (data not shown). (C) Sequence analysis of the aberrant 190-bp band demonstrates fusion of RET exon 7a to exon 14b with a heterogeneous 18-bp intervening sequence which consisted of nucleotides continuing beyond 1476(A) into exon 7b (upper sequence) or nucleotides from exon 14b upstream of mRNA position 2500 through 2518 (down sequence). This suggests aberrant splicing between sites after the ESE in exon 7 and two alternative cryptic splice sites (position 2500 and 2518) in exon 14b.
Table 1. Predicted effect of SNP 1476G/A on affinity score for binding of SR protein
To confirm that the minor RET1476(A) allele alters normal mRNAsplicing, we studied three RET-expressing human cell lines:(1) SK-N-BE(2) neuroblastoma cells (1476A/A); (2) SH-SY5Y neuroblastomacells (1476G/A); and (3) Wit49 Wilms tumor cells (1476G/G).RET transcripts were amplified with nested primers spanningexon 6 through exon 15. Only the predicted wild-type 1252-bptranscript is evident in Wit49 (1476G/G) cells, whereas an aberrant190-bp transcript is noted in the two cell lines bearing oneor more 1476A alleles (Figure 1B). The 190-bp transcript wasnot seen in any of four other 1476(G/G) cell lines (data notshown). Sequencing of the 1252-bp band in the heterozygous (rs1800860)SH-SY5Y cell line demonstrated that normal splicing of exon7 can occur with either nucleotide at mRNA position 1476 (Figure 2);however, when the minor transcript in SK-N-BE(2) or SH-SY5Ycells was sequenced, we also noted aberrant splicing from 1476Ato two cryptic splice sites in exon 14 (Figure 1C). This del(7b-14a)transcript eliminates part of the cadherin-like domain 4 (CLD4),the entire cysteine-rich domain (CRD), transmembrane domain,and the first tyrosine kinase domain, undoubtedly renderingthe protein product nonfunctional (Figure 3). The CLD4 and CRDare the extracellular sites of interaction with RET ligand (GDNF)and co-receptor (GFR1).12,13 To confirm that this aberrant splicingactivity occurred in other cell lines, we also studied the heterozygous(1476G/A) Wilms tumor G401 cell line and found the samedel(7b-14a) transcript (data not shown).
Figure 2. Sequence of wild-type human RET exon 7 (NM_020630). After reverse transcription of mRNA from heterozygous (rs1800860) SH-SY5Y cells, exon 7 was amplified by PCR and the normal 1252-bp amplicon was sequenced. The presence of either G or A at position 1476 demonstrates that normal splicing of RET exon 7 can occur with both alleles.
Figure 3. The aberrant del(7b-14a) RET transcript lacks key functional domains. CLD1 through 4, CRD, transmembrane segment (TM), juxtamembrane segment (JM), and two cytoplasmic tyrosine kinase domains (TK1 and TK2) are displayed adjacent to the corresponding exon structure of the wild-type RET transcript (WT). To the left is the aberrant del(7b-14a) transcript associated with 1476(A) alleles; deletion of extracellular domains for ligand (GDNF) and co-receptor (GFR1), transmembrane segments, and a large portion of TK1 is predicted, rendering the transcript nonfunctional.
Because both the 1476A and 2251A single-nucleotide polymorphisms(SNP) are potentially hypomorphic variants of the RET gene,we hypothesized that either one could compromise nephrogenesisduring kidney development, resulting in subtle renal hypoplasiain a significant fraction of newborns. To validate newborn kidneyvolume as a surrogate for congenital nephron number, we measuredglomerular number in kidneys from 15 infants who died of congenitalanomalies or sudden infant death syndrome at or before the ageof 3 mo. Nephron number ranged from 246,181 to 1,106,062. Astrong direct relationship (P = 0.019) between kidney mass andtotal nephron number in children aged 3 mo was identified (Figure 4).Regression analysis predicted an additional 23,459 (95% confidenceinterval 4590 to 42,238) glomeruli per gram of kidney mass inthe range.
Figure 4. Relationship between Nglom (with 95% confidence interval) and mass of right kidney at autopsy in infants aged 3 mo. Nglom was measured by the disector/fractionator method as described previously in autopsied kidneys from infants who died from various congenital anomalies or sudden infant death syndrome within the first 3 mo of life.
We then recruited 136 normal term newborn white infants bornat the Royal Victoria Hospital in Montreal between 2004 and2006, excluding any with low birth weight (<2500 g), renalmalformation, maternal diabetes, or twin pregnancy. Total (combined)renal volume was estimated by ultrasonography within 48 h ofbirth and normalized for body surface; cord blood was obtainedfor assay of cystatin C as a surrogate of newborn GFR.14–16In previous studies, we found that newborn kidney volume correlateddirectly with body surface area and inversely with cord bloodcystatin C.14 Clinical characteristics of the cohort are summarizedin Figure 5. The distribution of cystatin and kidney size inour cohort approached normality (skewness <2 for all parameters).
Figure 5. Clinical characteristics of the newborns of white cohort (n = 136). (A) Total KidVol/BSA in the white newborn cohort (mean ± SD 132.16 ± 29.34 ml/m2). (B) Cord blood cystatin C (mean ± SD 1.93 ± 0.32 mg/L). (C) Clinical characteristics of the newborn cohort.
DNA extracted from cord blood was used to genotype each infantfor the rs1800860 (1476G/A) and rs1799939 (2251G/A) SNP. Genotypedistribution for each SNP conformed to the Hardy-Weinberg equilibrium(P = 0.39 and 0.47, respectively); SNP frequencies in our cohort(Table 2) were similar to those reported in the CAUC1 and CEUpopulations (NCBI dbSNP database; see Appendix). RET2251(G/A)SNP was associated with neither renal volume nor cord cystatinC; however, total kidney volume factored for body surface area(KidVol/BSA) in newborns bearing one or more 1476(A) alleleswas 9.7% smaller than that in newborns with the homozygous 1476(G/G)genotype (P = 0.009). The 1476(A) allele was also associatedwith 9.2% increase in cord blood cystatin C concentration (P= 0.002), suggesting a comparable reduction in functional renalmass15 (Table 2).
Table 2. Association between SNP and newborn kidney volume or cord blood cystatin C
The 1476(A) allele was widely distributed in our population;combined renal volume in infants with a 1476(A) allele rangedfrom 16.06 to 44.48 ml versus 15.01 to 50.33 ml in infants whowere homozygous for the more common 1476(G) allele. Thus, theassociation between 1476(A) allele and renal volume cannot beattributed to a few infants with very small kidneys. Similarly,the distribution of left/right kidney volume ratios was similarin newborns with one or more 1476(A) alleles (0.93 ±0.22 SD) versus newborns homozygous for the 1476(G) allele (0.96± 0.18 SD) (P > 0.05). Thus, the association between1476(A) and renal volume was not due to a few infants with unilateralrenal hypoplasia.
To confirm that the RET1476(A) allele compromises expressionof wild-type receptor mRNA, we examined allele-specific mRNAexpression in the heterozygous SH-SY5Y and G401 cell lines.RET exon 7 was amplified from genomic DNA and cellular mRNA,using high-sensitivity sequencing technology to quantify eachallele-specific amplicon as described by others.17 mRNA expressionfrom the 1476(A) allele was substantially reduced in both SH-SY5Y(by 25%) and G401 (by 50%) cells, compared with the expressionlevel from the 1476(G) allele (Figure 6).
Figure 6. Comparison of allele-specific mRNA expression in rs1800860(G/A) heterozygous G401 and SH-SY5Y cells. Peak area ratio of G/A alleles was measured in cDNA and genomic DNA from two cell lines heterozygous for the rs1800860 SNP. The G/A allelic ratio (2.09 ± 0.12) in cDNA was 2.1 times greater than the ratio (1.01 ± 0.02) in genomic DNA from G401 cells (*P = 0.004). Similarly, G/A ratio was 1.4 times greater in cDNA (1.41 ± 0.14) than in genomic DNA (0.98 ± 0.03) for SH-SY5Y cells (**P = 0.03).
Our previous studies demonstrated an association between a common(18.5% of white individuals) human PAX2 haplotype (AAA) anda 10% reduction in newborn renal size.16 Interestingly, RETtranscription is regulated by PAX2 gene dosage.10 When we analyzedour cohort for both genes, renal volume in the 17 of 136 newbornscarrying both the RET1476(A) and PAX2(AAA) minor alleles was23% lower than that in infants with the RET1476(G/G), PAX2(GGG)haplotype (Figure 7).
Figure 7. Comparison of KidVol/BSA in newborns with various combinations of hypomorphic RET1476(A) and PAX2AAA alleles. All 136 newborns in our cohort were genotyped for both the RET(1476G/A) SNP and the hypomorphic PAX2AAA haplotype previously described by Quinlan et al.15 KidVol/BSA among newborns bearing one or more hypomorphic RET1476(A) alleles (n = 46) was 90.3% of that for infants with the major RET1476(G)/PAX2GGG alleles (n = 41; P = 0.01). Similarly, KidVol/BSA among newborns with one or more hypomorphic PAX2AAA alleles (n = 14) was 89.5% of RET1476(G)/PAX2GGG newborns (P = 0.04). In the subset of newborns bearing both hypomorphic RET1476(A) and PAX2AAA alleles (n = 17), KidVol/BSA was only 77% of that in wild-type RET1476(G)/PAX2GGG infants (P = 0.00067).
During renal development, the branching ureteric bud expresseshigh levels of the transcription factor PAX2.1 Brophy and colleagues10,18recently showed that, among its many functions, PAX2 directlyactivates transcription of genes for both RET, a tyrosine kinasereceptor, and its ligand, GDNF. Epithelial cells expressingRET receptors cluster at the tip of each ureteric bud branchas it undergoes branching morphogenesis.8 When activated byGDNF from nearby mesenchyme, the RET receptor heterodimerizeswith GFR1, stimulating cell proliferation, migration, and survivalvia several intracellular signals, including the RAS/mitogen-activatedprotein kinase, phosphatidylinositol 3-kinase/AKT, and RAC1/JUNNH(2) terminal kinase pathways.19–21 Loss of RET, GDNF,or GFR1 results in renal agenesis, because of inhibition ofureteric bud growth and branching.9,22,23
Because RET integrity is critical for branching nephrogenesis,we hypothesized that heterozygous RET mutations might partiallycompromise the extent of ureteric bud branching during development,leading to suboptimal nephron number. In mutant mice, heterozygousnull alleles reduce nephron number by approximately 22% at postnatalday 15.10 Thus, a heterozygous hypomorphic RET allele such asthe 1476(A) variant should fit at the milder end of this spectrumand might plausibly produce the observed 10% decrease in newbornnephron number.
Mouse studies suggest a graded relationship between kidney sizeand the level of RET function; Ret knockout mice are anephric,whereas severely hypomorphic Ret alleles such as Tyr1062Phe24or RetDN25 produce severe renal hypoplasia. Thus, only a subtleeffect on kidney size would be expected for a RET polymorphismpredicted to have only a modest effect on total RET mRNA level;however, the relationship between kidney volume and nephronnumber is difficult to establish in mice; nephrogenesis continuesfor at least 2 wk after birth and potentially overlaps witha period of postnatal compensatory hypertrophy. For example,Clarke et al.10 noted that heterozygous Ret(+/–) postnatalday 15 mice had a 22% reduction in nephron number but only a10% reduction in kidney volume (NS). To establish that newbornkidney size was a valid surrogate for nephron number in humans,we measured glomerular number and kidney weight in newborn infantswho died before 3 mo of age. This showed a strong correlationbetween renal mass and nephron number in the newborn period.
In this study, we identified a common SNP within the ESE ofexon 7, which modifies the fidelity of mRNA splicing to thenormal splice site in exon 8. An adenine nucleotide at position1476 of this ESE increases the risk for aberrant splicing toalternative sites in exon 14. When this occurs, exons containingthe crucial GFR1 binding site, transmembrane, and first tyrosinekinase domains all are deleted, undoubtedly rendering the RETreceptor dysfunctional. On the basis of our studies in culturedheterozygous human cells, the presence of a 1476A allele reduceswild-type mRNA expression by approximately 38%; in a heterozygouscell, this would amount to an overall reduction of functionalRET transcript by 19%. Arguably, branching morphogenesis ofthe ureteric bud during fetal kidney development might be reducedin proportion to this reduction of RET expression. In our cohortof normal white infants from Montreal, the 1476(A) minor RETallele, infants with one or more 1476(A) alleles exhibited a9.7% reduction in newborn kidney volume (normalized for bodysurface area) and 9.2% reduction in newborn renal function.Thus, the degree of renal hypoplasia observed in infants (9to 10%) with one or more 1476(A) alleles is roughly commensuratewith the reduction in functional RET mRNA expression (19%) measuredin vitro.
To put this in a clinical perspective, the report of Kelleret al.5 suggested that adults with essential hypertension areborn with 47% fewer nephrons compared with those who remainnormotensive. Thus, the effect of the 1476(A) allele could accountfor only one fifth of this clinically relevant congenital nephrondeficit. Clearly, additional genes are involved in setting nephronnumber during development.
Heterozygous null mutations of RET have been reported in humanswith Hirschprung disease. If arborization of the ureteric budis proportional to the level of RET expression during renaldevelopment, then one might expect that congenital nephron numberin Hirschprung disease should be reduced. There are reportsof unilateral renal aplasia and various renal malformationsin Hirschprung disease,26 but nephron number has not been carefullyassessed. By the time patients with Hirschprung disease areidentified, there has been ample time for compensatory renalhypertrophy to erase the initial relationship between kidneysize and congenital nephron number. Thus, in the report of Kelleret al.,5 describing 50% reduction in nephron number among peoplewith essential hypertension, adult kidney mass was similar tothat in normal control subjects. In contrast, we measured renalvolume and function at birth before the relationship can bemasked by postnatal compensatory hypertrophy.10 Indeed, ourfindings from autopsied newborns suggested that kidney weightcorrelates with congenital nephron number (adjusted for age)for up to 3 mo.
Ureteric bud cells express high levels of PAX2 during renaldevelopment; homozygous Pax2 mutant mice are anephric, and heterozygousPax2 null mutants exhibit significant renal hypoplasia.27–29We recently identified a fairly common (allele frequency 0.2)hypomorphic human PAX2AAA allele, which reduces PAX2 expressionby 40% compared with the major wild-type PAX2GGG allele; thus,total PAX2 mRNA in a heterozygous cell is approximately 80%of normal.10 Total newborn kidney volume among infants withone or more PAX2AAA alleles was approximately 10% smaller thanin PAX2GGG/ PAX2GGG infants.16 Fifteen percent of the newborncohort in this study were compound PAX2AAA/RET1476A heterozygotes.In this subgroup, newborn KidVol/BSA was 23% smaller than ininfants with homozygous wild-type RET1476G/PAX2GGG alleles.This corresponds very nicely with the observations of Clarkeet al.,10 who found additive effects of compound heterozygosityfor Pax2 and Ret alleles on nephron number in mutant mice. Thus,together, the two hypomorphic RET and PAX2 alleles could accountfor up to half of the nephron deficit (48% of control) associatedwith essential hypertension reported by Keller et al.5
On the basis of HapMap linkage data in the white (CEPH) population,the 1476(A) SNP is not in high genetic linkage disequilibriumwith most other portions of the RET gene; however, to rule outthe unlikely possibility that some other site within the RETgene is responsible for reduced kidney size, we screened forassociation with 22 other haplotype-tagging SNPs spanning allmajor linkage blocks and 10 kb to either side of the codingsequence. No other haplotype-tagging SNPs were significantlyassociated with either kidney size or cystatin C level (datanot shown).
Interestingly, we identified a second RET SNP (rs1799939) thatalters a glycine at amino acid position 691 to a serine (RET2251G/A);however, this SNP was not associated with newborn kidney sizeor cystatin C. Because the G/A substitution lies within a linkerregion between the RET transmembrane and tyrosine kinase domains,the 2251A isoform may retain sufficient signaling activity topermit normal nephrogenesis.
In conclusion, we identified a polymorphic variant (RET1476A)of the human RET receptor that increases the risk for aberrantmRNA splicing and causes decreased expression of the functionalwild-type allele. The RET1476A is associated with a reduction(approximately 10%) of newborn kidney volume and an increase(approximately 9%) in umbilical cord cystatin C. Because kidneysize and function were assessed at birth, before the periodof postnatal hypertrophy, these measurements likely reflectcongenital nephron endowment. In mice, homozygous RET mutationsblock ureteric bud outgrowth, and heterozygous RET null mutationshave been shown to interfere with optimal nephrogenesis.10 Similarly,we previously showed that heterozygous PAX2 mutations interferewith ureteric bud branching and reduce congenital nephron numberin mice; a common polymorphic variant of PAX2 causes 10% reductionin human newborn kidney volume.16 Among the 15% of normal newbornswho inherit both a hypomorphic RET1476(A) and a hypomorphicPAX2AAA allele, kidney volume is reduced by 23% of wild-typecontrols. Clarke et al.10 also noted a synergistic effect ofPax2 and Ret mutations on ureteric bud branching and nephrogenesisin embryonic mouse kidney explants. Our observations suggesta model in which common polymorphic variants of genes involvedin renal branching morphogenesis account for subtle renal hypoplasiain the normal human population.
Cell Culture
Human neuroblastoma cells [SK-N-BE(2) SH-SY5Y] and Wilmstumor cells (G401, Wit49) were obtained from ATCC (Manassas,VA). The cells were grown and maintained using standard mediumand conditions according to ATCC protocols.
Reverse Transcriptase–PCR Analysis
Total RNA was isolated from cells using Qiagen RNeasy Mini-plusKit with gDNA eliminator column (Qiagen, Mississauga, ON, Canada).Two-step reverse transcriptase–PCR was performed; first-strandcDNA was primed with random hexamers and TaqMan MultiScribeReverse Transcriptase according to the manufacturer's instructions(Applied Biosystems, Foster City, CA). Nested PCR was performedwith the cDNA; the PCR primers are listed in Supplemental Table1.
Study Populations
For analysis of glomerular number, we studied autopsied rightkidneys from 15 infants who died of congenital anomalies withinthe first 3 mo of life at the University of Mississippi MedicalCenter. Infants with morphologic abnormalities of the kidneyswere excluded. These studies were approved by the institutionalreview board of the University of Mississippi Medical Center.All autopsies were performed and kidney tissues were used withthe permission and informed consent of county coroners and next-of-kin.
For association studies, healthy white infants (n = 136) bornto women with uncomplicated pregnancies were recruited withinformed parental consent at the final prenatal clinical visitto the Royal Victoria Hospital (Montreal, QC, Canada). The study(PED-04-016) was approved by the Montreal Children's HospitalResearch Ethics Board. Mothers with twins, diabetes, intrauterinegrowth restriction, genetic abnormalities, renal malformations,hydronephrosis, or delivery at <36 wk and newborns with lowbirth weight (<2500 g) or low serum albumin were excluded.
Estimation of Total Glomerular (Nephron) Number
The right kidney was perfusion-fixed with 10% buffered formalinand then weighed. Kidneys were excluded from the study whenthe two kidneys in the one subject were unequal in size or showedmacroscopic or microscopic evidence of pathology. Kidneys from15 infants aged 3 mo were analyzed. Subject age, gender, race,nephron number, and kidney weight are shown in SupplementalTable 3.
After perfusion, kidneys were immersion-fixed in formalin andsent to Monash University for stereologic analysis of totalnephron number (Nglom) using the physical disector/fractionatorcombination. This is an unbiased stereologic counting methodwith which all glomeruli are sampled and, thereby, counted withequal probability. Important with this method, glomeruli arecounted irrespective of their size, shape, and location. Fulldetails of this technique have been previously described indetail.30–32 The association between Nglom and kidneyweight was analyzed using Pearson product moment correlation.
Kidney Volume Measurement
Left and right kidney volumes were measured by ultrasonographyin newborns within the first 48 h of life using the formulakidney volume = 4/3II (length/2) (height/2) (width/2). Bodysurface area was calculated as the square root of [length (cm)and weight (kg)/3600] according to Mosteller.33
Renal Function Determination
Serum cystatin C was used as a surrogate of GFR.34 Cord bloodcystatin C was measured by nephelometry (normal newborn range1.17 to 3.06 ± 0.26 mg/L [SD]).15
Coding SNP Collection
The NCBI dbSNP database was screened for common coding SNP inthe human RET gene with minor allele frequency of >10% inwhite populations. One common nonsynonymous SNP changing anamino acid was identified, and one common synonymous SNP affectingan ESE was found using the RESCU-ESE program (see Appendix).
SNP Genotyping
Genomic DNA was isolated from cord blood with the FlexiGeneDNA kit (Qiagen) according to the manufacturer's protocol. Foreach infant, 15 ng of genomic DNA was used for multiplex genotyping,using Sequenom iPLEX PCR technology (Sequenom, San Diego, CA).This system involves extension of the PCR amplicon with modifiednucleotides to distinguish SNP alleles by matrix-assisted laserdesorption ionization–time of flight technology. Primersfor SNP detection were designed using MassARRAY AssayDesignsoftware (Sequenom, San Diego, CA).
Quantitative Allele Ratio Analysis
Two heterozygous cell lines (SH-SY5Y and G401) were selectedfor allelic expression analysis. Total RNA and genomic DNA wereisolated in triplicate from the cells. PCR and reverse transcriptase–PCRamplicons were sequenced in duplicate as described by Pastinenet al.17 Sequencing primer sequences did not contain any knownSNP and were used to amplify each RNA and DNA sample in duplicate(the PCR primers and sequencing primer are summarized in SupplementalTable 1). PeakPicker software (see Appendix) was used for quantitativeallele ratio analysis. This program normalizes nucleotide peakamplitude for the effect of surrounding bases, and normalizedratio values were calculated in genomic DNA and mRNA (cDNA).
Statistical Analysis
Data are presented as means ± SD. Deviation from Hardy-Weinbergequilibrium was calculated by the 2 test. Normality of datadistributions for kidney volume and serum cystatin C were confirmedby tests of skewness (values of 0.501 and 0.482, respectively)and kurtosis (values of 0.216 and 0.133, respectively). Associationbetween SNP genotypes and KidVol/BSA or cord blood cystatinC was assessed by two-tailed, independent-samples t test. Thecomparison of rs1800860 G/A allele expression ratios in DNAand RNA was analyzed with two-tailed t test. All data were analyzedwith SPSS for Windows 11.0 (SPSS, Chicago, IL) and MicrosoftExcel. The association between Nglom and kidney weight was analyzedusing Pearson product moment correlation.
This work was supported by operating grants from the CanadianInstitutes of Health Research (MOP 12954) and the McGill UniversityHealth Centre Research Institute, funded in part by the Fondsde Recherches en Santé du Québec. Measurementsof glomerular number were funded by grants from National Institutesof Health grant 1 R01 DK065970-01, National Institutes of HealthCenter of Excellence in Minority Health 5P20M000534-02, andthe Colonial Foundation of Australia.
Studies were conducted with informed consent from subjects andthe approval of the Montreal Children's Hospital institutionalreview board (PED 04-016) and the institutional review boardof the University of Mississippi Medical Center.
Dr. Goodyer is the recipient of a James McGill Research Chair.
Footnotes
Published online ahead of print. Publication date availableat www.jasn.org.
Dziarmaga A, Quinlan J, Goodyer P: Renal hypoplasia: Lessons from Pax2.
Pediatr Nephrol 21
: 26
–31, 2005[CrossRef][Medline]
Clark AT, Bertram JF: Molecular regulation of nephron endowment.
Am J Physiol 276
: F485
–F497, 1999[Medline]
Brenner BM, Garcia DL, Anderson S: Glomeruli and blood pressure: Less of one, more the other?
Am J Hypertens 1
: 335
–347, 1988[Medline]
Mackenzie HS, Lawler EV, Brenner BM: Congenital oligonephropathy: The fetal flaw in essential hypertension?
Kidney Int Suppl 55
: S30
–S34, 1996[Medline]
Keller G, Zimmer G, Mall G, Ritz E, Amann K: Nephron number in patients with primary hypertension.
N Engl J Med 348
: 101
–108, 2003[Abstract/Free Full Text]
Hoy WE, Wang Z, VanBuynder P, Baker PR, McDonald SM, Mathews JD: The natural history of renal disease in Australian Aborigines. Part 2. Albuminuria predicts natural death and renal failure.
Kidney Int 60
: 249
–256, 2001[CrossRef][Medline]
Cass A, Cunningham J, Wang Z, Hoy W: Regional variation in the incidence of end-stage renal disease in Indigenous Australians.
Med J Aust 175
: 24
–27, 2001[Medline]
Costantini F, Shakya R: GDNF/Ret signaling and the development of the kidney.
Bioessays 28
: 117
–127, 2006[CrossRef][Medline]
Schuchardt A, D'Agati V, Pachnis V, Costantini F: Renal agenesis and hypodysplasia in ret-k- mutant mice result from defects in ureteric bud development.
Development 122
: 1919
–1929, 1996[Abstract]
Clarke JC, Patel SR, Raymond RM Jr, Andrew S, Robinson BG, Dressler GR, Brophy PD: Regulation of c-Ret in the developing kidney is responsive to Pax2 gene dosage.
Hum Mol Genet 15
: 3420
–3428, 2006[Abstract/Free Full Text]
Cartegni L, Wang J, Zhu Z, Zhang MQ, Krainer AR: ESEfinder: A web resource to identify exonic splicing enhancers.
Nucleic Acids Res 31
: 3568
–3571, 2003[Abstract/Free Full Text]
Anders J, Kjar S, Ibanez CF: Molecular modeling of the extracellular domain of the RET receptor tyrosine kinase reveals multiple cadherin-like domains and a calcium-binding site.
J Biol Chem 276
: 35808
–35817, 2001[Abstract/Free Full Text]
Amoresano A, Incoronato M, Monti G, Pucci P, de Franciscis V, Cerchia L: Direct interactions among Ret, GDNF and GFRalpha1 molecules reveal new insights into the assembly of a functional three-protein complex.
Cell Signal 17
: 717
–727, 2005[CrossRef][Medline]
Goodyer P, Kurpad A, Rekha S, Muthayya S, Dwarkanath P, Iyengar A, Philip B, Mhaskar A, Benjamin A, Maharaj S, Laforte D, Raju C, Phadke K: Effects of maternal vitamin A status on kidney development: A pilot study.
Pediatr Nephrol 22
: 209
–214, 2007[CrossRef][Medline]
Filler G, Bokenkamp A, Hofmann W, Le Bricon T, Martinez-Bru C, Grubb A: Cystatin C as a marker of GFR: History, indications, and future research.
Clin Biochem 38
: 1
–8, 2005[CrossRef][Medline]
Quinlan J, Lemire M, Hudson T, Qu H, Benjamin A, Roy A, Pascuet E, Goodyer M, Raju C, Zhang Z, Houghton F, Goodyer P: A common variant of the PAX2 gene is associated with reduced newborn kidney size.
J Am Soc Nephrol 18
: 1915
–1921, 2007[Abstract/Free Full Text]
Pastinen T, Ge B, Gurd S, Gaudin T, Dore C, Lemire M, Lepage P, Harmsen E, Hudson TJ: Mapping common regulatory variants to human haplotypes.
Hum Mol Genet 14
: 3963
–3971, 2005[Abstract/Free Full Text]
Brophy PD, Ostrom L, Lang KM, Dressler GR: Regulation of ureteric bud outgrowth by Pax2-dependent activation of the glial derived neurotrophic factor gene.
Development 128
: 4747
–4756, 2001[Abstract/Free Full Text]
Kim D, Dressler GR: PTEN modulates GDNF/RET mediated chemotaxis and branching morphogenesis in the developing kidney.
Dev Biol 307
: 290
–299, 2007[CrossRef][Medline]
Plaza-Menacho I, van der Sluis T, Hollema H, Gimm O, Buys CH, Magee AI, Isacke CM, Hofstra RM, Eggen BJ: Ras/ERK1/2-mediated STAT3 Ser727 phosphorylation by familial medullary thyroid carcinoma-associated RET mutants induces full activation of STAT3 and is required for c-fos promoter activation, cell mitogenicity, and transformation.
J Biol Chem 282
: 6415
–6424, 2007[Abstract/Free Full Text]
Asai N, Fukuda T, Wu Z, Enomoto A, Pachnis V, Takahashi M, Costantini F: Targeted mutation of serine 697 in the Ret tyrosine kinase causes migration defect of enteric neural crest cells.
Development 133
: 4507
–4516, 2006[Abstract/Free Full Text]
Sanchez MP, Silos-Santiago I, Frisen J, He B, Lira SA, Barbacid M: Renal agenesis and the absence of enteric neurons in mice lacking GDNF.
Nature 382
: 70
–73, 1996[CrossRef][Medline]
Cacalano G, Farinas I, Wang LC, Hagler K, Forgie A, Moore M, Armanini M, Phillips H, Ryan AM, Reichardt LF, Hynes M, Davies A, Rosenthal A: GFRalpha1 is an essential receptor component for GDNF in the developing nervous system and kidney.
Neuron 21
: 53
–62, 1998[CrossRef][Medline]
Jijiwa M, Fukuda T, Kawai K, Nakamura A, Kurokawa K, Murakumo Y, Ichihara M, Takahashi M: A targeting mutation of tyrosine 1062 in Ret causes a marked decrease of enteric neurons and renal hypoplasia.
Mol Cell Biol 24
: 8026
–8036, 2004[Abstract/Free Full Text]
Jain S, Naughton CK, Yang M, Strickland A, Vij K, Encinas M, Golden J, Gupta A, Heuckeroth R, Johnson EM Jr, Milbrandt J: Mice expressing a dominant-negative Ret mutation phenocopy human Hirschsprung disease and delineate a direct role of Ret in spermatogenesis.
Development 131
: 5503
–5513, 2004[Abstract/Free Full Text]
Sinnassamy P, Yazbeck S, Brochu P, O'Regan S: Renal anomalies and agenesis associated with total intestinal aganglionosis.
Int J Pediatr Nephrol 7
: 1
–2, 1986[Medline]
Dressler GR, Deutsch U, Chowdhury K, Nornes HO, Gruss P: Pax2, a new murine paired-box-containing gene and its expression in the developing excretory system.
Development 109
: 787
–795, 1990[Abstract/Free Full Text]
Porteous S, Torban E, Cho NP, Cunliffe H, Chua L, McNoe L, Ward T, Souza C, Gus P, Giugliani R, Sato T, Yun K, Favor J, Sicotte M, Goodyer P, Eccles M: Primary renal hypoplasia in humans and mice with PAX2 mutations: Evidence of increased apoptosis in fetal kidneys of Pax2(1Neu) +/- mutant mice.
Hum Mol Genet 9
: 1
–11, 2000[Abstract/Free Full Text]
Dziarmaga A, Eccles M, Goodyer P: Suppression of ureteric bud apoptosis rescues nephron endowment and adult renal function in Pax2 mutant mice.
J Am Soc Nephrol 17
: 1568
–1575, 2006[Abstract/Free Full Text]
Bertram JF: Analyzing renal glomeruli with the new stereology.
Int Rev Cytol 161
: 111
–172, 1995[Medline]
Hughson M, Farris AB 3rd, Douglas-Denton R, Hoy WE, Bertram JF: Glomerular number and size in autopsy kidneys: The relationship to birth weight.
Kidney Int 63
: 2113
–2122, 2003[CrossRef][Medline]
Nyengaard JR: Stereologic methods and their application in kidney research.
J Am Soc Nephrol 10
: 1100
–1123, 1999[Abstract/Free Full Text]
Mosteller RD: Simplified calculation of body-surface area.
N Engl J Med 317
: 1098
, 1987[Medline]
Harmoinen A, Ylinen E, Ala-Houhala M, Janas M, Kaila M, Kouri T: Reference intervals for cystatin C in pre- and full-term infants and children.
Pediatr Nephrol 15
: 105
–108, 2000[CrossRef][Medline]
Related Article
This Month's Highlights
J. Am. Soc. Nephrol. 2008 19: A6.
[Full Text][PDF]