Genomic Mismatch at LIMS1 Locus and Kidney Allograft Rejection.

BACKGROUND: In the context of kidney transplantation, genomic incompatibilities between donor and recipient may lead to allosensitization against new antigens. We hypothesized that recessive inheritance of gene-disrupting variants may represent a risk factor for allograft rejection. METHODS: We performed a two-stage genetic association study of kidney allograft rejection. In the first stage, we performed a recessive association screen of 50 common gene-intersecting deletion polymorphisms in a cohort of kidney transplant recipients. In the second stage, we replicated our findings in three independent cohorts of donor-recipient pairs. We defined genomic collision as a specific donor-recipient genotype combination in which a recipient who was homozygous for a gene-intersecting deletion received a transplant from a nonhomozygous donor. Identification of alloantibodies was performed with the use of protein arrays, enzyme-linked immunosorbent assays, and Western blot analyses. RESULTS: In the discovery cohort, which included 705 recipients, we found a significant association with allograft rejection at the LIMS1 locus represented by rs893403 (hazard ratio with the risk genotype vs. nonrisk genotypes, 1.84; 95% confidence interval [CI], 1.35 to 2.50; P = 9.8×10-5). This effect was replicated under the genomic-collision model in three independent cohorts involving a total of 2004 donor-recipient pairs (hazard ratio, 1.55; 95% CI, 1.25 to 1.93; P = 6.5×10-5). In the combined analysis (discovery cohort plus replication cohorts), the risk genotype was associated with a higher risk of rejection than the nonrisk genotype (hazard ratio, 1.63; 95% CI, 1.37 to 1.95; P = 4.7×10-8). We identified a specific antibody response against LIMS1, a kidney-expressed protein encoded within the collision locus. The response involved predominantly IgG2 and IgG3 antibody subclasses. CONCLUSIONS: We found that the LIMS1 locus appeared to encode a minor histocompatibility antigen. Genomic collision at this locus was associated with rejection of the kidney allograft and with production of anti-LIMS1 IgG2 and IgG3. (Funded by the Columbia University Transplant Center and others.).

Mutations in DSTYK and dominant urinary tract malformations.

BACKGROUND: Congenital abnormalities of the kidney and the urinary tract are the most common cause of pediatric kidney failure. These disorders are highly heterogeneous, and the etiologic factors are poorly understood. METHODS: We performed genomewide linkage analysis and whole-exome sequencing in a family with an autosomal dominant form of congenital abnormalities of the kidney or urinary tract (seven affected family members). We also performed a sequence analysis in 311 unrelated patients, as well as histologic and functional studies. RESULTS: Linkage analysis identified five regions of the genome that were shared among all affected family members. Exome sequencing identified a single, rare, deleterious variant within these linkage intervals, a heterozygous splice-site mutation in the dual serine-threonine and tyrosine protein kinase gene (DSTYK). This variant, which resulted in aberrant splicing of messenger RNA, was present in all affected family members. Additional, independent DSTYK mutations, including nonsense and splice-site mutations, were detected in 7 of 311 unrelated patients. DSTYK is highly expressed in the maturing epithelia of all major organs, localizing to cell membranes. Knockdown in zebrafish resulted in developmental defects in multiple organs, which suggested loss of fibroblast growth factor (FGF) signaling. Consistent with this finding is the observation that DSTYK colocalizes with FGF receptors in the ureteric bud and metanephric mesenchyme. DSTYK knockdown in human embryonic kidney cells inhibited FGF-stimulated phosphorylation of extracellular-signal-regulated kinase (ERK), the principal signal downstream of receptor tyrosine kinases. CONCLUSIONS: We detected independent DSTYK mutations in 2.3% of patients with congenital abnormalities of the kidney or urinary tract, a finding that suggests that DSTYK is a major determinant of human urinary tract development, downstream of FGF signaling. (Funded by the National Institutes of Health and others.).

A retrotransposon insertion in the 5' regulatory domain of Ptf1a results in ectopic gene expression and multiple congenital defects in Danforth's short tail mouse.

Danforth's short tail mutant (Sd) mouse, first described in 1930, is a classic spontaneous mutant exhibiting defects of the axial skeleton, hindgut, and urogenital system. We used meiotic mapping in 1,497 segregants to localize the mutation to a 42.8-kb intergenic segment on chromosome 2. Resequencing of this region identified an 8.5-kb early retrotransposon (ETn) insertion within the highly conserved regulatory sequences upstream of Pancreas Specific Transcription Factor, 1a (Ptf1a). This mutation resulted in up to tenfold increased expression of Ptf1a as compared to wild-type embryos at E9.5 but no detectable changes in the expression levels of other neighboring genes. At E9.5, Sd mutants exhibit ectopic Ptf1a expression in embryonic progenitors of every organ that will manifest a developmental defect: the notochord, the hindgut, and the mesonephric ducts. Moreover, at E 8.5, Sd mutant mice exhibit ectopic Ptf1a expression in the lateral plate mesoderm, tail bud mesenchyme, and in the notochord, preceding the onset of visible defects such as notochord degeneration. The Sd heterozygote phenotype was not ameliorated by Ptf1a haploinsufficiency, further suggesting that the developmental defects result from ectopic expression of Ptf1a. These data identify disruption of the spatio-temporal pattern of Ptf1a expression as the unifying mechanism underlying the multiple congenital defects in Danforth's short tail mouse. This striking example of an enhancer mutation resulting in profound developmental defects suggests that disruption of conserved regulatory elements may also contribute to human malformation syndromes.

Copy-number disorders are a common cause of congenital kidney malformations.

We examined the burden of large, rare, copy-number variants (CNVs) in 192 individuals with renal hypodysplasia (RHD) and replicated findings in 330 RHD cases from two independent cohorts. CNV distribution was significantly skewed toward larger gene-disrupting events in RHD cases compared to 4,733 ethnicity-matched controls (p = 4.8 × 10(-11)). This excess was attributable to known and novel (i.e., not present in any database or in the literature) genomic disorders. All together, 55/522 (10.5%) RHD cases harbored 34 distinct known genomic disorders, which were detected in only 0.2% of 13,839 population controls (p = 1.2 × 10(-58)). Another 32 (6.1%) RHD cases harbored large gene-disrupting CNVs that were absent from or extremely rare in the 13,839 population controls, identifying 38 potential novel or rare genomic disorders for this trait. Deletions at the HNF1B locus and the DiGeorge/velocardiofacial locus were most frequent. However, the majority of disorders were detected in a single individual. Genomic disorders were detected in 22.5% of individuals with multiple malformations and 14.5% of individuals with isolated urinary-tract defects; 14 individuals harbored two or more diagnostic or rare CNVs. Strikingly, the majority of the known CNV disorders detected in the RHD cohort have previous associations with developmental delay or neuropsychiatric diseases. Up to 16.6% of individuals with kidney malformations had a molecular diagnosis attributable to a copy-number disorder, suggesting kidney malformations as a sentinel manifestation of pathogenic genomic imbalances. A search for pathogenic CNVs should be considered in this population for the diagnosis of their specific genomic disorders and for the evaluation of the potential for developmental delay.

Combined linkage and association mapping reveals CYCD5;1 as a quantitative trait gene for endoreduplication in Arabidopsis.

Endoreduplication is the process where a cell replicates its genome without mitosis and cytokinesis, often followed by cell differentiation. This alternative cell cycle results in various levels of endoploidy, reaching 4× or higher one haploid set of chromosomes. Endoreduplication is found in animals and is widespread in plants, where it plays a major role in cellular differentiation and plant development. Here, we show that variation in endoreduplication between Arabidopsis thaliana accessions Columbia-0 and Kashmir is controlled by two major quantitative trait loci, ENDO-1 and ENDO-2. A local candidate gene association analysis in a set of 87 accessions, combined with expression analysis, identified CYCD5;1 as the most likely candidate gene underlying ENDO-2, operating as a rate-determining factor of endoreduplication. In accordance, both the overexpression and silencing of CYCD5;1 were effective in changing DNA ploidy levels, confirming CYCD5;1 to be a previously undescribed quantitative trait gene underlying endoreduplication in Arabidopsis.

APOL1 variants increase risk for FSGS and HIVAN but not IgA nephropathy.

A chromosome 22q13 locus strongly associates with increased risk for idiopathic focal segmental glomerulosclerosis (FSGS), HIV-1-associated nephropathy (HIVAN), and hypertensive ESRD among individuals of African descent. Although initial studies implicated MYH9, more recent analyses localized the strongest association within the neighboring APOL1 gene. In this replication study, we examined the six top-most associated variants in APOL1 and MYH9 in an independent cohort of African Americans with various nephropathies (44 with FSGS, 21 with HIVAN, 32 with IgA nephropathy, and 74 healthy controls). All six variants associated with FSGS and HIVAN (additive ORs, 1.8 to 3.0; P values 3 × 10(-2) to 5 × 10(-5)) but not with IgA nephropathy. In conditional and haplotype analyses, two APOL1 haplotypes accounted for virtually all of the association with FSGS and HIVAN on chromosome 22q13 (haplotype P value = 5.6 × 10(-8)). To assess the role of MYH9 deficiency in nephropathy, we crossbred Myh9-haploinsufficient mice (Myh9(+/-)) with HIV-1 transgenic mice. Myh9(+/-) mice were healthy and did not demonstrate overt proteinuria or nephropathy, irrespective of the presence of the HIV-1 transgene. These data further support the strong association of genetic variants in APOL1 with susceptibility to FSGS and HIVAN among African Americans.

Identification of the nephropathy-susceptibility locus HIVAN4.

HIVAN1, HIVAN2, and HIVAN3 are nephropathy-susceptibility loci previously identified in the HIV-1 transgenic mouse, a model of collapsing glomerulopathy. The HIVAN1 and HIVAN2 loci modulate expression of Nphs2, which encodes podocin and several other podocyte-expressed genes. To identify additional loci predisposing to nephropathy, we performed a genome-wide scan in 165 backcross mice generated between the nephropathy-sensitive HIV-1-transgenic FVB/NJ (TgFVB) strain and the resistant Balb/cJ (BALB) strain. We identified a major susceptibility locus (HIVAN4) on chromosome 6 G3-F3, with BALB alleles conferring a twofold reduction in severity (peak LOD score = 4.0). Similar to HIVAN1 and HIVAN2, HIVAN4 modulated expression of Nphs2, indicating a common pathway underlying these loci. We independently confirmed the HIVAN4 locus in a sister TgFVB colony that experienced a dramatic loss of nephropathy subsequent to a breeding bottleneck. In this low-penetrance line, 3% of the genome was admixed with BALB alleles, suggesting a remote contamination event. The admixture localized to discrete segments on chromosome 2 and at the HIVAN4 locus. HIVAN4 candidate genes include killer lectin-like receptor genes as well as A2m and Ptpro, whose gene products are enriched in the glomerulus and interact with HIV-1 proteins. In summary, these data identify HIVAN4 as a major quantitative trait locus for nephropathy and a transregulator of Nphs2. Furthermore, similar selective breeding strategies may help identify further susceptibility loci.

The molecular pathogenesis of HIV-1 associated nephropathy: recent advances.

HIV-1-associated nephropathy (HIVAN) is a major complication of HIV-1 infection, frequently resulting in kidney failure. HIVAN arises due to HIV-1-induced dysregulation of podocytes, the glomerular epithelial cells that establish and maintain the kidney filtration barrier. Host genetic factors are important for the development of HIVAN. The risk of HIVAN is greatest in populations of African ancestry, and is attributable to a genetic variation at the APOL1 locus on chromosome 22. Mouse models of HIVAN enable delineation of dysregulated pathways underlying disease. Identification of HIVAN susceptibility loci in a mouse model, combined with expression quantitative trait locus mapping, has demonstrated that murine HIVAN loci transregulate podocyte gene expression. HIV-1 induces perturbations in podocyte expression response, suggesting that HIV-1 potentially interferes with compensatory pathways that normally restore cellular homeostasis in the face of genetic mutations. These findings present a framework for identification of podocyte transregulators and reconstruction of the molecular networks connecting susceptibility genes to the development of nephropathy.

The genetics of albuminuria: from haplotype association mapping in mice to genetic causation in humans.

Genome-wide haplotype association mapping (HAM) in inbred mouse strains emerged as an efficient method for identifying novel quantitative trait loci for disease-related phenotypes. In this issue of Kidney International, Tsaih et al. present the results of the first HAM for age-related kidney damage in mice and examine the detected loci in the context of the human genome-wide association study for diabetic nephropathy.

A population genomics study of the Arabidopsis core cell cycle genes shows the signature of natural selection.

Large-scale comparison of sequence polymorphism and divergence at numerous genomic loci within and between closely related species can reveal signatures of natural selection. Here, we present a population genomics study based on direct sequencing of 61 mitotic cell cycle genes from 30 Arabidopsis thaliana accessions and comparison of the resulting data to the close relative Arabidopsis lyrata. We found that the Arabidopsis core cell cycle (CCC) machinery is not highly constrained but is subject to different modes of selection. We found patterns of purifying selection for the cyclin-dependent kinase (CDK), CDK subunit, retinoblastoma, and WEE1 gene families. Other CCC gene families often showed a mix of one or two constrained genes and relaxed purifying selection on the other genes. We found several large effect mutations in CDKB1;2 that segregate in the species. We found a strong signature of adaptive protein evolution in the Kip-related protein KRP6 and departures from equilibrium at CDKD;1 and CYCA3;3 consistent with the operation of selection in these gene regions. Our data suggest that within Arabidopsis, the genetic robustness of cell cycle-related processes is more due to functional redundancy than high selective constraint.

Genome-wide screening for cis-regulatory variation using a classical diallel crossing scheme.

Large-scale screening studies carried out to date for genetic variants that affect gene regulation are generally limited to descriptions of differences in allele-specific expression (ASE) detected in vivo. Allele-specific differences in gene expression provide evidence for a model whereby cis-acting genetic variation results in differential expression between alleles. Such gene surveys for regulatory variation are a first step in identifying the specific nucleotide changes that govern gene expression differences, but they leave the underlying mechanisms unexplored. Here, we propose a quantitative genetics approach to perform a genome-wide analysis of ASE differences (GASED). The GASED approach is based on a diallel design that is often used in plant breeding programs to estimate general combining abilities (GCA) of specific inbred lines and to identify high-yielding hybrid combinations of parents based on their specific combining abilities (SCAs). In a context of gene expression, the values of GCA and SCA parameters allow cis- and trans-regulatory changes to be distinguished and imbalances in gene expression to be ascribed to cis-regulatory variation. With this approach, a total of 715 genes could be identified that are likely to carry allelic polymorphisms responsible for at least a 1.5-fold allelic expression difference in a total of 10 diploid Arabidopsis thaliana hybrids. The major strength of the GASED approach, compared to other ASE detection methods, is that it is not restricted to genes with allelic transcript variants. Although a false-positive rate of 9/41 was observed, the GASED approach is a valuable pre-screening method that can accelerate systematic surveys of naturally occurring cis-regulatory variation among inbred lines for laboratory species, such as Arabidopsis, mouse, rat and fruitfly, and economically important crop species, such as corn.