Bardet-Biedl syndrome (BBS) is a defining ciliopathy, notable for extensive allelic and genetic heterogeneity, almost all of which has been identified through sequencing. Recent data have suggested that copy-number variants (CNVs) also contribute to BBS. We used a custom oligonucleotide array comparative genomic hybridization (aCGH) covering 20 genes that encode intraflagellar transport (IFT) components and 74 ciliopathy loci to screen 92 unrelated individuals with BBS, irrespective of their known mutational burden. We identified 17 individuals with exon-disruptive CNVs (18.5%), including 13 different deletions in eight BBS genes (BBS1, BBS2, ARL6/BBS3, BBS4, BBS5, BBS7, BBS9, and NPHP1) and a deletion and a duplication in other ciliopathy-associated genes (ALMS1 and NPHP4, respectively). By contrast, we found a single heterozygous exon-disruptive event in a BBS-associated gene (BBS9) in 229 control subjects. Superimposing these data with resequencing revealed CNVs to (1) be sufficient to cause disease, (2) Mendelize heterozygous deleterious alleles, and (3) contribute oligogenic alleles by combining point mutations and exonic CNVs in multiple genes. Finally, we report a deletion and a splice site mutation in IFT74, inherited under a recessive paradigm, defining a candidate BBS locus. Our data suggest that CNVs contribute pathogenic alleles to a substantial fraction of BBS-affected individuals and highlight how either deletions or point mutations in discrete splice isoforms can induce hypomorphic mutations in genes otherwise intolerant to deleterious variation. Our data also suggest that CNV analyses and resequencing studies unbiased for previous mutational burden is necessary to delineate the complexity of disease architecture.

Patterns of amino acid conservation have served as a tool for understanding protein evolution1. The same principles have also found broad application in human genomics, driven by the need to interpret the pathogenic potential of variants in patients2. Here we performed a systematic comparative genomics analysis of human disease-causing missense variants. We found that an appreciable fraction of disease-causing alleles are fixed in the genomes of other species, suggesting a role for genomic context. We developed a model of genetic interactions that predicts most of these to be simple pairwise compensations. Functional testing of this model on two known human disease genes3,4 revealed discrete cis amino acid residues that, although benign on their own, could rescue the human mutations in vivo. This approach was also applied to ab initio gene discovery to support the identification of a de novo disease driver in BTG2 that is subject to protective cis-modification in more than 50 species. Finally, on the basis of our data and models, we developed a computational tool to predict candidate residues subject to compensation. Taken together, our data highlight the importance of cis-genomic context as a contributor to protein evolution; they provide an insight into the complexity of allele effect on phenotype; and they are likely to assist methods for predicting allele pathogenicity5,6.

Cardiac two-pore domain potassium channels (K2P) exist in organisms from Drosophila to humans; however, their role in cardiac function is not known. We identified a K2P gene, CG8713 (sandman), in a Drosophila genetic screen and show that sandman is critical to cardiac function. Mice lacking an ortholog of sandman, TWIK-related potassium channel (TREK-1, also known Kcnk2), exhibit exaggerated pressure overload-induced concentric hypertrophy and alterations in fetal gene expression, yet retain preserved systolic and diastolic cardiac function. While cardiomyocyte-specific deletion of TREK-1 in response to in vivo pressure overload resulted in cardiac dysfunction, TREK-1 deletion in fibroblasts prevented deterioration in cardiac function. The absence of pressure overload-induced dysfunction in TREK-1-KO mice was associated with diminished cardiac fibrosis and reduced activation of JNK in cardiomyocytes and fibroblasts. These findings indicate a central role for cardiac fibroblast TREK-1 in the pathogenesis of pressure overload-induced cardiac dysfunction and serve as a conceptual basis for its inhibition as a potential therapy.

Rapid advances and cost erosion in exome and genome analysis of patients with both rare and common genetic disorders has accelerated gene discovery and illuminated fundamental biological mechanisms. The thrill of discovery has been accompanied, however, with the sobering appreciation that human genomes are burdened with a large number of rare and ultra rare variants, thereby posing a significant challenge in dissecting both the effect of such alleles on protein function and the biological relevance of these events to patient pathology. In an effort to develop model systems that are able to generate surrogates of human pathologies, a powerful suite of tools have been developed in zebrafish, taking advantage of the relatively small (compared to invertebrate models) evolutionary distance of that genome to humans, the orthology of several organs and signaling processes, and the suitability of this organism for medium and high throughput phenotypic screening. Here we will review the use of this model organism in dissecting human genetic disorders; we will highlight how diverse strategies have informed disease causality and genetic architecture; and we will discuss relative strengths and limitations of these approaches in the context of medical genome sequencing.