Sequencing of LRP2 reveals multiple rare variants associated with urinary trefoil factor-3

Novel biomarkers are being investigated to identify patients with kidney disease. We measured a panel of 13 urinary biomarkers in participants from the Offspring Cohort of the Framingham Heart Study. Using an Affymetrix chip with imputation to 2.5 M single-nucleotide polymorphisms (SNPs), we conducted a GWAS of these biomarkers (n=2640) followed by exonic sequencing and genotyping. Functional studies in zebrafish were used to investigate histologic correlation with renal function. Across all 13 biomarkers, there were 97 significant SNPs at three loci. Lead SNPs at each locus were rs6555820 (P=6.7×10(-49); minor allele frequency [MAF]=0.49) in HAVCR1 (associated with kidney injury molecule-1), rs7565788 (P=2.15×10(-16); MAF=0.22) in LRP2 (associated with trefoil factor 3 [TFF3]), and rs11048230 (P=4.77×10(-8); MAF=0.10) in an intergenic region near RASSF8 (associated with vascular endothelial growth factor). Validation in the CKDGen Consortium (n=67,093) showed that only rs7565788 at LRP2, which encodes megalin, was associated with eGFR (P=0.003). Sequencing of exons 16-72 of LRP2 in 200 unrelated individuals at extremes of urinary TFF3 levels identified 197 variants (152 rare; MAF<0.05), 31 of which (27 rare) were nonsynonymous. In aggregate testing, rare variants were associated with urinary TFF3 levels (P=0.003), and the lead GWAS signal was not explained by these variants. Knockdown of LRP2 in zebrafish did not alter the renal phenotype in static or kidney injury models. In conclusion, this study revealed common variants associated with urinary levels of TFF3, kidney injury molecule-1, and vascular endothelial growth factor and identified a cluster of rare variants independently associated with TFF3.