Gain and loss of TASK3 channel function and its regulation by novel variation cause KCNK9 imprinting syndrome.

BACKGROUND: Genomics enables individualized diagnosis and treatment, but large challenges remain to functionally interpret rare variants. To date, only one causative variant has been described for KCNK9 imprinting syndrome (KIS). The genotypic and phenotypic spectrum of KIS has yet to be described and the precise mechanism of disease fully understood. METHODS: This study discovers mechanisms underlying KCNK9 imprinting syndrome (KIS) by describing 15 novel KCNK9 alterations from 47 KIS-affected individuals. We use clinical genetics and computer-assisted facial phenotyping to describe the phenotypic spectrum of KIS. We then interrogate the functional effects of the variants in the encoded TASK3 channel using sequence-based analysis, 3D molecular mechanic and dynamic protein modeling, and in vitro electrophysiological and functional methodologies. RESULTS: We describe the broader genetic and phenotypic variability for KIS in a cohort of individuals identifying an additional mutational hotspot at p.Arg131 and demonstrating the common features of this neurodevelopmental disorder to include motor and speech delay, intellectual disability, early feeding difficulties, muscular hypotonia, behavioral abnormalities, and dysmorphic features. The computational protein modeling and in vitro electrophysiological studies discover variability of the impact of KCNK9 variants on TASK3 channel function identifying variants causing gain and others causing loss of conductance. The most consistent functional impact of KCNK9 genetic variants, however, was altered channel regulation. CONCLUSIONS: This study extends our understanding of KIS mechanisms demonstrating its complex etiology including gain and loss of channel function and consistent loss of channel regulation. These data are rapidly applicable to diagnostic strategies, as KIS is not identifiable from clinical features alone and thus should be molecularly diagnosed. Furthermore, our data suggests unique therapeutic strategies may be needed to address the specific functional consequences of KCNK9 variation on channel function and regulation.

CSNK2B: A broad spectrum of neurodevelopmental disability and epilepsy severity.

CSNK2B has recently been implicated as a disease gene for neurodevelopmental disability (NDD) and epilepsy. Information about developmental outcomes has been limited by the young age and short follow-up for many of the previously reported cases, and further delineation of the spectrum of associated phenotypes is needed. We present 25 new patients with variants in CSNK2B and refine the associated NDD and epilepsy phenotypes. CSNK2B variants were identified by research or clinical exome sequencing, and investigators from different centers were connected via GeneMatcher. Most individuals had developmental delay and generalized epilepsy with onset in the first 2 years. However, we found a broad spectrum of phenotypic severity, ranging from early normal development with pharmacoresponsive seizures to profound intellectual disability with intractable epilepsy and recurrent refractory status epilepticus. These findings suggest that CSNK2B should be considered in the diagnostic evaluation of patients with a broad range of NDD with treatable or intractable seizures.

Atlas-CNV: a validated approach to call single-exon CNVs in the eMERGESeq gene panel.

PURPOSE: To provide a validated method to confidently identify exon-containing copy-number variants (CNVs), with a low false discovery rate (FDR), in targeted sequencing data from a clinical laboratory with particular focus on single-exon CNVs. METHODS: DNA sequence coverage data are normalized within each sample and subsequently exonic CNVs are identified in a batch of samples, when the target log2 ratio of the sample to the batch median exceeds defined thresholds. The quality of exonic CNV calls is assessed by C-scores (Z-like scores) using thresholds derived from gold standard samples and simulation studies. We integrate an ExonQC threshold to lower FDR and compare performance with alternate software (VisCap). RESULTS: Thirteen CNVs were used as a truth set to validate Atlas-CNV and compared with VisCap. We demonstrated FDR reduction in validation, simulation, and 10,926 eMERGESeq samples without sensitivity loss. Sixty-four multiexon and 29 single-exon CNVs with high C-scores were assessed by Multiplex Ligation-dependent Probe Amplification (MLPA). CONCLUSION: Atlas-CNV is validated as a method to identify exonic CNVs in targeted sequencing data generated in the clinical laboratory. The ExonQC and C-score assignment can reduce FDR (identification of targets with high variance) and improve calling accuracy of single-exon CNVs respectively. We propose guidelines and criteria to identify high confidence single-exon CNVs.

ZMIZ1 Variants Cause a Syndromic Neurodevelopmental Disorder.

ZMIZ1 is a coactivator of several transcription factors, including p53, the androgen receptor, and NOTCH1. Here, we report 19 subjects with intellectual disability and developmental delay carrying variants in ZMIZ1. The associated features include growth failure, feeding difficulties, microcephaly, facial dysmorphism, and various other congenital malformations. Of these 19, 14 unrelated subjects carried de novo heterozygous single-nucleotide variants (SNVs) or single-base insertions/deletions, 3 siblings harbored a heterozygous single-base insertion, and 2 subjects had a balanced translocation disrupting ZMIZ1 or involving a regulatory region of ZMIZ1. In total, we identified 13 point mutations that affect key protein regions, including a SUMO acceptor site, a central disordered alanine-rich motif, a proline-rich domain, and a transactivation domain. All identified variants were absent from all available exome and genome databases. In vitro, ZMIZ1 showed impaired coactivation of the androgen receptor. In vivo, overexpression of ZMIZ1 mutant alleles in developing mouse brains using in utero electroporation resulted in abnormal pyramidal neuron morphology, polarization, and positioning, underscoring the importance of ZMIZ1 in neural development and supporting mutations in ZMIZ1 as the cause of a rare neurodevelopmental syndrome.

Phenotypic expansion in DDX3X - a common cause of intellectual disability in females.

De novo variants in DDX3X account for 1-3% of unexplained intellectual disability (ID) cases and are amongst the most common causes of ID especially in females. Forty-seven patients (44 females, 3 males) have been described. We identified 31 additional individuals carrying 29 unique DDX3X variants, including 30 postnatal individuals with complex clinical presentations of developmental delay or ID, and one fetus with abnormal ultrasound findings. Rare or novel phenotypes observed include respiratory problems, congenital heart disease, skeletal muscle mitochondrial DNA depletion, and late-onset neurologic decline. Our findings expand the spectrum of DNA variants and phenotypes associated with DDX3X disorders.

De novo mutations in the SET nuclear proto-oncogene, encoding a component of the inhibitor of histone acetyltransferases (INHAT) complex in patients with nonsyndromic intellectual disability.

The role of disturbed chromatin remodeling in the pathogenesis of intellectual disability (ID) is well established and illustrated by de novo mutations found in a plethora of genes encoding for proteins of the epigenetic regulatory machinery. We describe mutations in the "SET nuclear proto-oncogene" (SET), encoding a component of the "inhibitor of histone acetyltransferases" (INHAT) complex, involved in transcriptional silencing. Using whole exome sequencing, four patients were identified with de novo mutations in the SET gene. Additionally, an affected mother and child were detected who carried a frameshift variant in SET. Four patients were found in literature. The de novo mutations in patients affected all four known SET mRNA transcripts. LoF mutations in SET are exceedingly rare in the normal population and, if present, affect only one transcript. The pivotal role of SET in neurogenesis is evident from in vitro and animal models. SET interacts with numerous proteins involved in histone modification, including proteins encoded by known autosomal dominant ID genes, that is, EP300, CREBBP, SETBP1, KMT2A, RAC1, and CTCF. Our study identifies SET as a new component of epigenetic regulatory modules underlying human cognitive disorders, and as a first member of the Nucleosome Assembly Protein (NAP) family implicated in ID.

De novo apparent loss-of-function mutations in PRR12 in three patients with intellectual disability and iris abnormalities.

PRR12 encodes a proline-rich protein nuclear factor suspected to be involved in neural development. Its nuclear expression in fetal brains and in the vision system supports its role in brain and eye development more specifically. However, its function and potential role in human disease has not been determined. Recently, a de novo t(10;19) (q22.3;q13.33) translocation disrupting the PRR12 gene was detected in a girl with intellectual disability and neuropsychiatric alterations. Here we report on three unrelated patients with heterozygous de novo apparent loss-of-function mutations in PRR12 detected by clinical whole exome sequencing: c.1918G>T (p.Glu640*), c.4502_4505delTGCC (p.Leu1501Argfs*146) and c.903_909dup (p.Pro304Thrfs*46). All three patients had global developmental delay, intellectual disability, eye and vision abnormalities, dysmorphic features, and neuropsychiatric problems. Eye abnormalities were consistent among the three patients and consisted of stellate iris pattern and iris coloboma. Additional variable clinical features included hypotonia, skeletal abnormalities, sleeping problems, and behavioral issues such as autism and anxiety. In summary, we propose that haploinsufficiency of PRR12 is associated with this novel multisystem neurodevelopmental disorder.

Mutations in signal recognition particle SRP54 cause syndromic neutropenia with Shwachman-Diamond-like features.

Shwachman-Diamond syndrome (SDS) (OMIM #260400) is a rare inherited bone marrow failure syndrome (IBMFS) that is primarily characterized by neutropenia and exocrine pancreatic insufficiency. Seventy-five to ninety percent of patients have compound heterozygous loss-of-function mutations in the Shwachman-Bodian-Diamond syndrome (sbds) gene. Using trio whole-exome sequencing (WES) in an sbds-negative SDS family and candidate gene sequencing in additional SBDS-negative SDS cases or molecularly undiagnosed IBMFS cases, we identified 3 independent patients, each of whom carried a de novo missense variant in srp54 (encoding signal recognition particle 54 kDa). These 3 patients shared congenital neutropenia linked with various other SDS phenotypes. 3D protein modeling revealed that the 3 variants affect highly conserved amino acids within the GTPase domain of the protein that are critical for GTP and receptor binding. Indeed, we observed that the GTPase activity of the mutated proteins was impaired. The level of SRP54 mRNA in the bone marrow was 3.6-fold lower in patients with SRP54-mutations than in healthy controls. Profound reductions in neutrophil counts and chemotaxis as well as a diminished exocrine pancreas size in a SRP54-knockdown zebrafish model faithfully recapitulated the human phenotype. In conclusion, autosomal dominant mutations in SRP54, a key member of the cotranslation protein-targeting pathway, lead to syndromic neutropenia with a Shwachman-Diamond-like phenotype.

Use of Exome Sequencing for Infants in Intensive Care Units: Ascertainment of Severe Single-Gene Disorders and Effect on Medical Management.

Importance: While congenital malformations and genetic diseases are a leading cause of early infant death, to our knowledge, the contribution of single-gene disorders in this group is undetermined. Objective: To determine the diagnostic yield and use of clinical exome sequencing in critically ill infants. Design, Setting, and Participants: Clinical exome sequencing was performed for 278 unrelated infants within the first 100 days of life who were admitted to Texas Children's Hospital in Houston, Texas, during a 5-year period between December 2011 and January 2017. Exome sequencing types included proband exome, trio exome, and critical trio exome, a rapid genomic assay for seriously ill infants. Main Outcomes and Measures: Indications for testing, diagnostic yield of clinical exome sequencing, turnaround time, molecular findings, patient age at diagnosis, and effect on medical management among a group of critically ill infants who were suspected to have genetic disorders. Results: The mean (SEM) age for infants participating in the study was 28.5 (1.7) days; of these, the mean (SEM) age was 29.0 (2.2) days for infants undergoing proband exome sequencing, 31.5 (3.9) days for trio exome, and 22.7 (3.9) days for critical trio exome. Clinical indications for exome sequencing included a range of medical concerns. Overall, a molecular diagnosis was achieved in 102 infants (36.7%) by clinical exome sequencing, with relatively low yield for cardiovascular abnormalities. The diagnosis affected medical management for 53 infants (52.0%) and had a substantial effect on informed redirection of care, initiation of new subspecialist care, medication/dietary modifications, and furthering life-saving procedures in select patients. Critical trio exome sequencing revealed a molecular diagnosis in 32 of 63 infants (50.8%) at a mean (SEM) of 33.1 (5.6) days of life with a mean (SEM) turnaround time of 13.0 (0.4) days. Clinical care was altered by the diagnosis in 23 of 32 patients (71.9%). The diagnostic yield, patient age at diagnosis, and medical effect in the group that underwent critical trio exome sequencing were significantly different compared with the group who underwent regular exome testing. For deceased infants (n = 81), genetic disorders were molecularly diagnosed in 39 (48.1%) by exome sequencing, with implications for recurrence risk counseling. Conclusions and Relevance: Exome sequencing is a powerful tool for the diagnostic evaluation of critically ill infants with suspected monogenic disorders in the neonatal and pediatric intensive care units and its use has a notable effect on clinical decision making.

Clinical and molecular characterization of de novo loss of function variants in HNRNPU.

DNA alterations in the 1q43-q44 region are associated with syndromic neurodevelopmental disorders characterized by global developmental delay, intellectual disability, dysmorphic features, microcephaly, seizures, and agenesis of the corpus callosum. HNRNPU is located within the 1q43-q44 region and mutations in the gene have been reported in patients with early infantile epileptic encephalopathy. Here, we report on the clinical presentation of four patients with de novo heterozygous HNRNPU loss-of-function mutations detected by clinical whole exome sequencing: c.651_660del (p.Gly218Alafs*118), c.1089G>A (p.Trp363*), c.1714C>T (p.Arg572*), and c.2270_2271del (p.Pro757Argfs*7). All patients shared similar clinical features as previously reported including seizures, global developmental delay, intellectual disability, variable neurologic regression, behavior issues, and dysmorphic facial features. Features including heart defects and kidney abnormalities were not reported in our patients. These findings expands the clinical spectrum of HNRNPU-related disorder and shows that HNRNPU contributes to a subset of the clinical phenotypes associated with the contiguous 1q43-q44 deletion syndrome.

CRIPT exonic deletion and a novel missense mutation in a female with short stature, dysmorphic features, microcephaly, and pigmentary abnormalities.

Mutations in CRIPT encoding cysteine-rich PDZ domain-binding protein are rare, and to date have been reported in only two patients with autosomal recessive primordial dwarfism and distinctive facies. Here, we describe a female with biallelic mutations in CRIPT presenting with postnatal growth retardation, global developmental delay, and dysmorphic features including frontal bossing, high forehead, and sparse hair and eyebrows. Additional clinical features included high myopia, admixed hyper- and hypopigmented macules primarily on the face, arms, and legs, and syndactyly of 4-5 toes bilaterally. Using whole exome sequencing (WES) and chromosomal microarray analysis (CMA), we detected a c.8G>A (p.C3Y) missense variant in exon 1 of the CRIPT gene inherited from the mother and a 1,331 bp deletion encompassing exon 1, inherited from the father. The c.8G>A (p.C3Y) missense variant in CRIPT was apparently homozygous in the proband due to the exon 1 deletion. Our findings illustrate the clinical utility of combining WES with copy number variant (CNV) analysis to provide a molecular diagnosis to patients with rare Mendelian disorders. Our findings also illustrate the clinical spectrum of CRIPT related mutations. © 2016 Wiley Periodicals, Inc.

Recurrent Muscle Weakness with Rhabdomyolysis, Metabolic Crises, and Cardiac Arrhythmia Due to Bi-allelic TANGO2 Mutations.

The underlying genetic etiology of rhabdomyolysis remains elusive in a significant fraction of individuals presenting with recurrent metabolic crises and muscle weakness. Using exome sequencing, we identified bi-allelic mutations in TANGO2 encoding transport and Golgi organization 2 homolog (Drosophila) in 12 subjects with episodic rhabdomyolysis, hypoglycemia, hyperammonemia, and susceptibility to life-threatening cardiac tachyarrhythmias. A recurrent homozygous c.460G>A (p.Gly154Arg) mutation was found in four unrelated individuals of Hispanic/Latino origin, and a homozygous ∼34 kb deletion affecting exons 3-9 was observed in two families of European ancestry. One individual of mixed Hispanic/European descent was found to be compound heterozygous for c.460G>A (p.Gly154Arg) and the deletion of exons 3-9. Additionally, a homozygous exons 4-6 deletion was identified in a consanguineous Middle Eastern Arab family. No homozygotes have been reported for these changes in control databases. Fibroblasts derived from a subject with the recurrent c.460G>A (p.Gly154Arg) mutation showed evidence of increased endoplasmic reticulum stress and a reduction in Golgi volume density in comparison to control. Our results show that the c.460G>A (p.Gly154Arg) mutation and the exons 3-9 heterozygous deletion in TANGO2 are recurrent pathogenic alleles present in the Latino/Hispanic and European populations, respectively, causing considerable morbidity in the homozygotes in these populations.

De Novo GMNN Mutations Cause Autosomal-Dominant Primordial Dwarfism Associated with Meier-Gorlin Syndrome.

Meier-Gorlin syndrome (MGS) is a genetically heterogeneous primordial dwarfism syndrome known to be caused by biallelic loss-of-function mutations in one of five genes encoding pre-replication complex proteins: ORC1, ORC4, ORC6, CDT1, and CDC6. Mutations in these genes cause disruption of the origin of DNA replication initiation. To date, only an autosomal-recessive inheritance pattern has been described in individuals with this disorder, with a molecular etiology established in about three-fourths of cases. Here, we report three subjects with MGS and de novo heterozygous mutations in the 5' end of GMNN, encoding the DNA replication inhibitor geminin. We identified two truncating mutations in exon 2 (the 1(st) coding exon), c.16A>T (p.Lys6(∗)) and c.35_38delTCAA (p.Ile12Lysfs(∗)4), and one missense mutation, c.50A>G (p.Lys17Arg), affecting the second-to-last nucleotide of exon 2 and possibly RNA splicing. Geminin is present during the S, G2, and M phases of the cell cycle and is degraded during the metaphase-anaphase transition by the anaphase-promoting complex (APC), which recognizes the destruction box sequence near the 5' end of the geminin protein. All three GMNN mutations identified alter sites 5' to residue Met28 of the protein, which is located within the destruction box. We present data supporting a gain-of-function mechanism, in which the GMNN mutations result in proteins lacking the destruction box and hence increased protein stability and prolonged inhibition of replication leading to autosomal-dominant MGS.

Clinical whole-exome sequencing for the diagnosis of mendelian disorders.

BACKGROUND: Whole-exome sequencing is a diagnostic approach for the identification of molecular defects in patients with suspected genetic disorders. METHODS: We developed technical, bioinformatic, interpretive, and validation pipelines for whole-exome sequencing in a certified clinical laboratory to identify sequence variants underlying disease phenotypes in patients. RESULTS: We present data on the first 250 probands for whom referring physicians ordered whole-exome sequencing. Patients presented with a range of phenotypes suggesting potential genetic causes. Approximately 80% were children with neurologic phenotypes. Insurance coverage was similar to that for established genetic tests. We identified 86 mutated alleles that were highly likely to be causative in 62 of the 250 patients, achieving a 25% molecular diagnostic rate (95% confidence interval, 20 to 31). Among the 62 patients, 33 had autosomal dominant disease, 16 had autosomal recessive disease, and 9 had X-linked disease. A total of 4 probands received two nonoverlapping molecular diagnoses, which potentially challenged the clinical diagnosis that had been made on the basis of history and physical examination. A total of 83% of the autosomal dominant mutant alleles and 40% of the X-linked mutant alleles occurred de novo. Recurrent clinical phenotypes occurred in patients with mutations that were highly likely to be causative in the same genes and in different genes responsible for genetically heterogeneous disorders. CONCLUSIONS: Whole-exome sequencing identified the underlying genetic defect in 25% of consecutive patients referred for evaluation of a possible genetic condition. (Funded by the National Human Genome Research Institute.).

A major X-linked locus affects kidney function in mice.

Chronic kidney disease is a common disease with increasing prevalence in the western population. One common reason for chronic kidney failure is diabetic nephropathy. Diabetic nephropathy and hyperglycemia are characteristics of the mouse inbred strain KK/HlJ, which is predominantly used as a model for metabolic syndrome due to its inherited glucose intolerance and insulin resistance. We used KK/HlJ, an albuminuria-sensitive strain, and C57BL/6J, an albuminuria-resistant strain, to perform a quantitative trait locus (QTL) cross to identify the genetic basis for chronic kidney failure. Albumin-creatinine ratio (ACR) was measured in 130 F2 male offspring. One significant QTL was identified on chromosome (Chr) X and four suggestive QTL were found on Chrs 6, 7, 12, and 13. Narrowing of the QTL region was focused on the X-linked QTL and performed by incorporating genotype and expression analyses for genes located in the region. From the 485 genes identified in the X-linked QTL region, a few candidate genes were identified using a combination of bioinformatic evidence based on genomic comparison of the parental strains and known function in urine homeostasis. Finally, this study demonstrates the significance of the X chromosome in the genetic determination of albuminuria.

Using bioinformatics and systems genetics to dissect HDL-cholesterol genetics in an MRL/MpJ x SM/J intercross.

A higher incidence of coronary artery disease is associated with a lower level of HDL-cholesterol. We searched for genetic loci influencing HDL-cholesterol in F2 mice from a cross between MRL/MpJ and SM/J mice. Quantitative trait loci (QTL) mapping revealed one significant HDL QTL (Apoa2 locus), four suggestive QTL on chromosomes 10, 11, 13, and 18 and four additional QTL on chromosomes 1 proximal, 3, 4, and 7 after adjusting HDL for the strong Apoa2 locus. A novel nonsynonymous polymorphism supports Lipg as the QTL gene for the chromosome 18 QTL, and a difference in Abca1 expression in liver tissue supports it as the QTL gene for the chromosome 4 QTL. Using weighted gene co-expression network analysis, we identified a module that after adjustment for Apoa2, correlated with HDL, was genetically determined by a QTL on chromosome 11, and overlapped with the HDL QTL. A combination of bioinformatics tools and systems genetics helped identify several candidate genes for both the chromosome 11 HDL and module QTL based on differential expression between the parental strains, cis regulation of expression, and causality modeling. We conclude that integrating systems genetics to a more-traditional genetics approach improves the power of complex trait gene identification.

Integration of QTL and bioinformatic tools to identify candidate genes for triglycerides in mice.

To identify genetic loci influencing lipid levels, we performed quantitative trait loci (QTL) analysis between inbred mouse strains MRL/MpJ and SM/J, measuring triglyceride levels at 8 weeks of age in F2 mice fed a chow diet. We identified one significant QTL on chromosome (Chr) 15 and three suggestive QTL on Chrs 2, 7, and 17. We also carried out microarray analysis on the livers of parental strains of 282 F2 mice and used these data to find cis-regulated expression QTL. We then narrowed the list of candidate genes under significant QTL using a "toolbox" of bioinformatic resources, including haplotype analysis; parental strain comparison for gene expression differences and nonsynonymous coding single nucleotide polymorphisms (SNP); cis-regulated eQTL in livers of F2 mice; correlation between gene expression and phenotype; and conditioning of expression on the phenotype. We suggest Slc25a7 as a candidate gene for the Chr 7 QTL and, based on expression differences, five genes (Polr3 h, Cyp2d22, Cyp2d26, Tspo, and Ttll12) as candidate genes for Chr 15 QTL. This study shows how bioinformatics can be used effectively to reduce candidate gene lists for QTL related to complex traits.

The mouse QTL map helps interpret human genome-wide association studies for HDL cholesterol.

Genome-wide association (GWA) studies represent a powerful strategy for identifying susceptibility genes for complex diseases in human populations but results must be confirmed and replicated. Because of the close homology between mouse and human genomes, the mouse can be used to add evidence to genes suggested by human studies. We used the mouse quantitative trait loci (QTL) map to interpret results from a GWA study for genes associated with plasma HDL cholesterol levels. We first positioned single nucleotide polymorphisms (SNPs) from a human GWA study on the genomic map for mouse HDL QTL. We then used mouse bioinformatics, sequencing, and expression studies to add evidence for one well-known HDL gene (Abca1) and three newly identified genes (Galnt2, Wwox, and Cdh13), thus supporting the results of the human study. For GWA peaks that occur in human haplotype blocks with multiple genes, we examined the homologous regions in the mouse to prioritize the genes using expression, sequencing, and bioinformatics from the mouse model, showing that some genes were unlikely candidates and adding evidence for candidate genes Mvk and Mmab in one haplotype block and Fads1 and Fads2 in the second haplotype block. Our study highlights the value of mouse genetics for evaluating genes found in human GWA studies.

A Bayesian framework for inference of the genotype-phenotype map for segregating populations.

Complex genetic interactions lie at the foundation of many diseases. Understanding the nature of these interactions is critical to developing rational intervention strategies. In mammalian systems hypothesis testing in vivo is expensive, time consuming, and often restricted to a few physiological endpoints. Thus, computational methods that generate causal hypotheses can help to prioritize targets for experimental intervention. We propose a Bayesian statistical method to infer networks of causal relationships among genotypes and phenotypes using expression quantitative trait loci (eQTL) data from genetically randomized populations. Causal relationships between network variables are described with hierarchical regression models. Prior distributions on the network structure enforce graph sparsity and have the potential to encode prior biological knowledge about the network. An efficient Monte Carlo method is used to search across the model space and sample highly probable networks. The result is an ensemble of networks that provide a measure of confidence in the estimated network topology. These networks can be used to make predictions of system-wide response to perturbations. We applied our method to kidney gene expression data from an MRL/MpJ × SM/J intercross population and predicted a previously uncharacterized feedback loop in the local renin-angiotensin system.