Inferring rare disease risk variants based on exact probabilities of sharing by multiple affected relatives

Statistics to prioritize rare variants in family-based sequencing studies with disease subtypes

Family-based sequencing studies are increasingly used to find rare genetic variants of high risk for disease traits with familial clustering. In some studies, families with multiple disease subtypes are collected and the exomes of affected relatives are sequenced for shared rare variants (RVs). Since different families can harbor different causal variants and each family harbors many RVs, tests to detect causal variants can have low power in this study design. Our goal is rather to prioritize shared variants for further investigation by, for example, pathway analyses or functional studies. The transmission-disequilibrium test prioritizes variants based on departures from Mendelian transmission in parent-child trios. Extending this idea to families, we propose methods to prioritize RVs shared in affected relatives with two disease subtypes, with one subtype more heritable than the other. Global approaches condition on a variant being observed in the study and assume a known probability of carrying a causal variant. In contrast, local approaches condition on a variant being observed in specific families to eliminate the carrier probability. Our simulation results indicate that global approaches are robust to misspecification of the carrier probability and prioritize more effectively than local approaches even when the carrier probability is misspecified.

Analysis of high-risk pedigrees identifies 11 candidate variants for Alzheimer's disease

Introduction

Analysis of sequence data in high‐risk pedigrees is a powerful approach to detect rare predisposition variants.

Methods

Rare, shared candidate predisposition variants were identified from exome sequencing 19 Alzheimer's disease (AD)‐affected cousin pairs selected from high‐risk pedigrees. Variants were further prioritized by risk association in various external datasets. Candidate variants emerging from these analyses were tested for co‐segregation to additional affected relatives of the original sequenced pedigree members.

Results

AD‐affected high‐risk cousin pairs contained 564 shared rare variants. Eleven variants spanning 10 genes were prioritized in external datasets: rs201665195 (ABCA7), and rs28933981 (TTR) were previously implicated in AD pathology; rs141402160 (NOTCH3) and rs140914494 (NOTCH3) were previously reported; rs200290640 (PIDD1) and rs199752248 (PIDD1) were present in more than one cousin pair; rs61729902 (SNAP91), rs140129800 (COX6A2, AC026471), and rs191804178 (MUC16) were not present in a longevity cohort; and rs148294193 (PELI3) and rs147599881 (FCHO1) approached significance from analysis of AD‐related phenotypes. Three variants were validated via evidence of co‐segregation to additional relatives (PELI3, ABCA7, and SNAP91).

Discussion

These analyses support ABCA7 and TTR as AD risk genes, expand on previously reported NOTCH3 variant identification, and prioritize seven additional candidate variants.

An intronic variant in the CELF4 gene is associated with risk for colorectal cancer

BACKGROUND

Germline predisposition variants associated with colorectal cancer (CRC) have been identified but all are not yet identified. We sought to identify the responsible predisposition germline variant in an extended high-risk CRC pedigree that exhibited evidence of linkage to the 18q12.2 region (TLOD = +2.81).

METHODS

DNA from two distantly related carriers of the hypothesized predisposition haplotype on 18q12.2 was sequenced to identify candidate variants. The candidate rare variants shared by the related sequenced subjects were screened in 3,094 CRC cases and 5x population-matched controls from UKBiobank to test for association. Further segregation of the variant was tested via Taqman assay in other sampled individuals in the pedigree.

RESULTS

Analysis of whole genome sequence data for the two related hypothesized predisposition haplotype carriers, restricted to the shared haplotype boundaries, identified multiple (n = 6) rare candidate non-coding variants that were tested for association with CRC risk in UKBiobank. A rare intronic variant ofCELF4 gene, rs568643870, was significantly associated with CRC (p = 0.004, OR = 5.0), and segregated with CRC in other members of the linked pedigree.

CONCLUSION

Evidence of segregation in a high-risk pedigree, case-control association in an external dataset, and identification of additional CRC-affected carriers in the linked pedigree support a role for a rareCELF4 intronic variant in CRC risk.

A role for the MEGF6 gene in predisposition to osteoporosis

Osteoporosis is a common skeletal disorder characterized by deterioration of bone tissue. The set of genetic factors contributing to osteoporosis is not completely specified. High-risk osteoporosis pedigrees were analyzed to identify genes that may confer susceptibility to disease. Candidate predisposition variants were identified initially by whole exome sequencing of affected-relative pairs, approximately cousins, from 10 pedigrees. Variants were filtered on the basis of population frequency, concordance between pairs of cousins, affecting a gene associated with osteoporosis, and likelihood to have functionally damaging, pathogenic consequences. Subsequently, variants were tested for segregation in 68 additional relatives of the index carriers. A rare variant in MEGF6 (rs755467862) showed strong evidence of segregation with the disease phenotype. Predicted protein folding indicated the variant (Cys200Tyr) may disrupt structure of an EGF-like calcium-binding domain of MEGF6. Functional analyses demonstrated that complete loss of the paralogous genes megf6a and megf6b in zebrafish resulted in significant delay of cartilage and bone formation. Segregation analyses, in silico protein structure modeling, and functional assays support a role for MEGF6 in predisposition to osteoporosis.

FANCM c5791C>T stopgain mutation (rs144567652) is a familial colorectal cancer risk factor

Purpose

While familial aggregation of colorectal cancer (CRC) is recognized, the majority of the germline predisposition factors remain unidentified, and many high‐risk CRC pedigrees remain unexplained by known risk variants. Fanconi Anemia genes have been recognized to be associated with cancer risk. Notably, FANCM (OMIM 609644) variants have been reported to confer risk for CRC and breast cancer.

Methods

Exome sequencing of CRC‐affected cousins in a set of 47 independent extended high‐risk CRC pedigrees identified a candidate set of rare, shared variants. Variants were tested for association with risk in 744 Utah CRC cases and 1525 controls, and for segregation with CRC in affected relatives.

Results

A FANCM stopgain variant was observed in two CRC‐affected cousin pairs, each from an independent Utah high‐risk pedigree, and yielded a nonsignificant, but elevated OR = 2.05 in a set of Utah cases and controls. Segregation of the variant to other related CRC‐affected cases was observed in the two extended pedigrees.

Conclusion

A rare stopgain variant in FANCM (rs144567652) that is recognized as a breast cancer predisposition variant, and that has previously been proposed, but not confirmed, as a CRC predisposition variant, is validated here as a risk factor for familial CRC.

Identification and genomic analysis of pedigrees with exceptional longevity identifies candidate rare variants

BACKGROUND

Longevity as a phenotype entails living longer than average and typically includes living without chronic age-related diseases. Recently, several common genetic components to longevity have been identified. This study aims to identify additional genetic variants associated with longevity using unique and powerful analyses of pedigrees with a statistical excess of healthy elderly individuals identified in the Utah Population Database (UPDB).

METHODS

From an existing biorepository of Utah pedigrees, six independent cousin pairs were selected from four extended pedigrees that exhibited an excess of healthy elderly individuals; whole exome sequencing (WES) was performed on two elderly individuals from each pedigree who were either first cousins or first cousins once removed. Rare (<.01 population frequency) variants shared by at least one elderly cousin pair in a region likely to be identical by descent were identified as candidates. Ingenuity Variant Analysis was used to prioritize putative causal variants based on quality control, frequency, and gain or loss of function. The variant frequency was compared in healthy cohorts and in an Alzheimer's disease cohort. Remaining variants were filtered based on their presence in genes reported to have an effect on the aging process, aging of cells, or the longevity process. Validation of these candidate variants included tests of segregation on other elderly relatives.

RESULTS

Fifteen rare candidate genetic variants spanning 17 genes shared within cousins were identified as having passed prioritization criteria. Of those variants, six were present in genes that are known or predicted to affect the aging process: rs78408340 (PAM), rs112892337 (ZFAT), rs61737629 (ESPL1), rs141903485 (CEBPE), rs144369314 (UTP4), and rs61753103 (NUP88 and RABEP1). ESPL1 rs61737629 and CEBPE rs141903485 show additional evidence of segregation with longevity in expanded pedigree analyses (p-values = .001 and .0001, respectively).

DISCUSSION

This unique pedigree analysis efficiently identified several novel rare candidate variants that may affect the aging process and added support to seven genes that likely contribute to longevity. Further analyses showed evidence for segregation for two rare variants, ESPL1 rs61737629 and CEBPE rs141903485, in the original longevity pedigrees in which they were initially observed. These candidate genes and variants warrant further investigation.

Genome-wide Enrichment of De Novo Coding Mutations in Orofacial Cleft Trios

Although de novo mutations (DNMs) are known to increase an individual's risk of congenital defects, DNMs have not been fully explored regarding orofacial clefts (OFCs), one of the most common human birth defects. Therefore, whole-genome sequencing of 756 child-parent trios of European, Colombian, and Taiwanese ancestry was performed to determine the contributions of coding DNMs to an individual's OFC risk. Overall, we identified a significant excess of loss-of-function DNMs in genes highly expressed in craniofacial tissues, as well as genes associated with known autosomal dominant OFC syndromes. This analysis also revealed roles for zinc-finger homeobox domain and SOX2-interacting genes in OFC etiology.

Detection of rare disease variants in extended pedigrees using RVS

Summary

Family-based sequencing studies enable researchers to identify highly penetrant genetic variants too rare to be tested in conventional case-control studies, by studying co-segregation of variant and disease phenotypes. When multiple affected subjects in a family are sequenced, the probability that a variant or a set of variants is shared identical-by-descent by some or all affected relatives provides evidence against the null hypothesis of complete absence of linkage and association. The Rare Variant Sharing software package RVS implements a suite of tools to assess association and linkage between rare genetic variants and a dichotomous disease indicator in family pedigrees.

Availability and Implementation

RVS is available as open source software from the Bioconductor webpage at https://bioconductor.org/packages/release/bioc/html/RVS.html.

Supplementary information

Supplementary data are available at Bioinformatics online.

Contribution of common and rare variants to bipolar disorder susceptibility in extended pedigrees from population isolates

Current evidence from case/control studies indicates that genetic risk for psychiatric disorders derives primarily from numerous common variants, each with a small phenotypic impact. The literature describing apparent segregation of bipolar disorder (BP) in numerous multigenerational pedigrees suggests that, in such families, large-effect inherited variants might play a greater role. To identify roles of rare and common variants on BP, we conducted genetic analyses in 26 Colombia and Costa Rica pedigrees ascertained for bipolar disorder 1 (BP1), the most severe and heritable form of BP. In these pedigrees, we performed microarray SNP genotyping of 838 individuals and high-coverage whole-genome sequencing of 449 individuals. We compared polygenic risk scores (PRS), estimated using the latest BP1 genome-wide association study (GWAS) summary statistics, between BP1 individuals and related controls. We also evaluated whether BP1 individuals had a higher burden of rare deleterious single-nucleotide variants (SNVs) and rare copy number variants (CNVs) in a set of genes related to BP1. We found that compared with unaffected relatives, BP1 individuals had higher PRS estimated from BP1 GWAS statistics (P = 0.001 ~ 0.007) and displayed modest increase in burdens of rare deleterious SNVs (P = 0.047) and rare CNVs (P = 0.002 ~ 0.033) in genes related to BP1. We did not observe rare variants segregating in the pedigrees. These results suggest that small-to-moderate effect rare and common variants are more likely to contribute to BP1 risk in these extended pedigrees than a few large-effect rare variants.

Identification of a Rare PSEN1 Mutation (Thr119Ile) in Late-Onset Alzheimer's Disease With Early Presentation of Behavioral Disturbance

Alzheimer's disease (AD) is the most common form of neurodegenerative dementia. In this study, whole genome sequencing identifies one rare and likely pathogenic mutation in the presenilin 1 (PSEN1) gene (c.356C > T, p.T119I) associated with a frontal variant of AD. Affected individuals in the kindred developed late-onset cognitive decline accompanied with early presentation of psychiatric symptoms. Positive amyloid PiB PET tracing suggested presence of pathophysiological biomarker for AD. Whole genome sequencing analysis evaluated rare coding mutations in susceptible genes for various types of dementia and supported the role of PSEN1 as a causal gene. Identification of this T119I variant in PSEN1 might broaden the spectrum of genetic basis and clinical diversity of familial AD.

A Frameshift Variant in the CHST9 Gene Identified by Family-Based Whole Genome Sequencing Is Associated with Schizophrenia in Chinese Population

Recent studies imply that rare variants contribute to the risk of schizophrenia, however, the exact variants or genes responsible for this condition are largely unknown. In this study, we conducted whole genome sequencing (WGS) of 20 Chinese families. Each family consisted of at least two affected siblings diagnosed with schizophrenia and at least one unaffected sibling. We examined functional variants that were found in affected sibling(s) but not in unaffected sibling(s) within a family. Matching this criterion, a frameshift heterozygous deletion of CA (–/CA) at chromosome 18:24722722, also referred to as rs752084147, in the Carbohydrate Sulfotransferase 9 (CHST9) gene, was detected in two families. This deletion was confirmed by PCR-based Sanger sequencing. With the observed frequency of 0.00076 in Han Chinese population, we performed both case-control and family-based analyses to evaluate its association with schizophrenia. In the case-control analyses, Chi-square test P-value was 6.80e-12 and the P-value was 0.0008 after one million simulations. In family-based segregation analyses, segregation P-value was 7.72e-7 and simulated P-value was 5.70e-6. For both the case-control and family-based analyses, the CA deletion was significantly associated with schizophrenia in the Chinese population. Further investigation of this gene  is warranted in the development of schizophrenia by utilizing larger and more ethnically diverse samples.

Inferring disease risk genes from sequencing data in multiplex pedigrees through sharing of rare variants

We previously demonstrated how sharing of rare variants (RVs) in distant affected relatives can be used to identify variants causing a complex and heterogeneous disease. This approach tested whether single RVs were shared by all sequenced affected family members. However, as with other study designs, joint analysis of several RVs (e.g., within genes) is sometimes required to obtain sufficient statistical power. Further, phenocopies can lead to false negatives for some causal RVs if complete sharing among affected is required. Here, we extend our methodology (Rare Variant Sharing, RVS) to address these issues. Specifically, we introduce gene-based analyses, a partial sharing test based on RV sharing probabilities for subsets of affected relatives and a haplotype-based RV definition. RVS also has the desirable feature of not requiring external estimates of variant frequency or control samples, provides functionality to assess and address violations of key assumptions, and is available as open source software for genome-wide analysis. Simulations including phenocopies, based on the families of an oral cleft study, revealed the partial and complete sharing versions of RVS achieved similar statistical power compared with alternative methods (RareIBD and the Gene-Based Segregation Test), and had superior power compared with the pedigree Variant Annotation, Analysis, and Search Tool (pVAAST) linkage statistic. In studies of multiplex cleft families, analysis of rare single nucleotide variants in the exome of 151 affected relatives from 54 families revealed no significant excess sharing in any one gene, but highlighted different patterns of sharing revealed by the complete and partial sharing tests.

Simulating pedigrees ascertained for multiple disease-affected relatives

Background

Studies that ascertain families containing multiple relatives affected by disease can be useful for identification of causal, rare variants from next-generation sequencing data.

Results

We present the R package SimRVPedigree, which allows researchers to simulate pedigrees ascertained on the basis of multiple, affected relatives. By incorporating the ascertainment process in the simulation, SimRVPedigree allows researchers to better understand the within-family patterns of relationship amongst affected individuals and ages of disease onset.

Conclusions

Through simulation, we show that affected members of a family segregating a rare disease variant tend to be more numerous and cluster in relationships more closely than those for sporadic disease. We also show that the family ascertainment process can lead to apparent anticipation in the age of onset. Finally, we use simulation to gain insight into the limit on the proportion of ascertained families segregating a causal variant. SimRVPedigree should be useful to investigators seeking insight into the family-based study design through simulation.

Electronic supplementary material

The online version of this article (10.1186/s13029-018-0069-6) contains supplementary material, which is available to authorized users.

Exome sequencing of a large family identifies potential candidate genes contributing risk to bipolar disorder

Bipolar disorder is a mental illness with lifetime prevalence of about 1%. Previous genetic studies have identified multiple chromosomal linkage regions and candidate genes that might be associated with bipolar disorder. The present study aimed to identify potential susceptibility variants for bipolar disorder using 6 related case samples from a four-generation family. A combination of exome sequencing and linkage analysis was performed to identify potential susceptibility variants for bipolar disorder. Our study identified a list of five potential candidate genes for bipolar disorder. Among these five genes, GRID1(Glutamate Receptor Delta-1 Subunit), which was previously reported to be associated with several psychiatric disorders and brain related traits, is particularly interesting. Variants with functional significance in this gene were identified from two cousins in our bipolar disorder pedigree. Our findings suggest a potential role for these genes and the related rare variants in the onset and development of bipolar disorder in this one family. Additional research is needed to replicate these findings and evaluate their patho-biological significance.

Gene-based segregation method for identifying rare variants in family-based sequencing studies

Whole-exome sequencing using family data has identified rare coding variants in Mendelian diseases or complex diseases with Mendelian subtypes, using filters based on variant novelty, functionality, and segregation with the phenotype within families. However, formal statistical approaches are limited. We propose a gene-based segregation test (GESE) that quantifies the uncertainty of the filtering approach. It is constructed using the probability of segregation events under the null hypothesis of Mendelian transmission. This test takes into account different degrees of relatedness in families, the number of functional rare variants in the gene, and their minor allele frequencies in the corresponding population. In addition, a weighted version of this test allows incorporating additional subject phenotypes to improve statistical power. We show via simulations that the GESE and weighted GESE tests maintain appropriate type I error rate, and have greater power than several commonly used region-based methods. We apply our method to whole-exome sequencing data from 49 extended pedigrees with severe, early-onset chronic obstructive pulmonary disease (COPD) in the Boston Early-Onset COPD study (BEOCOPD) and identify several promising candidate genes. Our proposed methods show great potential for identifying rare coding variants of large effect and high penetrance for family-based sequencing data. The proposed tests are implemented in an R package that is available on CRAN (https://cran.r-project.org/web/packages/GESE/).

Whole exome association of rare deletions in multiplex oral cleft families

By sequencing the exomes of distantly related individuals in multiplex families, rare mutational and structural changes to coding DNA can be characterized and their relationship to disease risk can be assessed. Recently, several rare single nucleotide variants (SNVs) were associated with an increased risk of nonsyndromic oral cleft, highlighting the importance of rare sequence variants in oral clefts and illustrating the strength of family-based study designs. However, the extent to which rare deletions in coding regions of the genome occur and contribute to risk of nonsyndromic clefts is not well understood. To identify putative structural variants underlying risk, we developed a pipeline for rare hemizygous deletions in families from whole exome sequencing and statistical inference based on rare variant sharing. Among 56 multiplex families with 115 individuals, we identified 53 regions with one or more rare hemizygous deletions. We found 45 of the 53 regions contained rare deletions occurring in only one family member. Members of the same family shared a rare deletion in only eight regions. We also devised a scalable global test for enrichment of shared rare deletions.

Genetic factors influencing risk to orofacial clefts: today's challenges and tomorrow's opportunities

Orofacial clefts include cleft lip (CL), cleft palate (CP), and cleft lip and palate (CLP), which combined represent the largest group of craniofacial malformations in humans with an overall prevalence of one per 1,000 live births. Each of these birth defects shows strong familial aggregation, suggesting a major genetic component to their etiology. Genetic studies of orofacial clefts extend back centuries, but it has proven difficult to define any single etiologic mechanism because many genes appear to influence risk. Both linkage and association studies have identified several genes influencing risk, but these differ across families and across populations. Genome-wide association studies have identified almost two dozen different genes achieving genome-wide significance, and there are broad classes of 'causal genes' for orofacial clefts: a few genes strongly associated with risk and possibly directly responsible for Mendelian syndromes which include orofacial clefts as a key phenotypic feature of the syndrome, and multiple genes with modest individual effects on risk but capable of disrupting normal craniofacial development under the right circumstances (which may include exposure to environmental risk factors). Genomic sequencing studies are now underway which will no doubt reveal additional genes/regions where variants (sequence and structural) can play a role in controlling risk to orofacial clefts. The real challenge to medicine and public health is twofold: to identify specific genes and other etiologic factors in families with affected members and then to devise effective interventions for these different biological mechanisms controlling risk to complex and heterogeneous birth defects such as orofacial clefts.

Genome-wide association of familial prostate cancer cases identifies evidence for a rare segregating haplotype at 8q24.21

Previous genome-wide association studies (GWAS) of prostate cancer risk focused on cases unselected for family history and have reported over 100 significant associations. The International Consortium for Prostate Cancer Genetics (ICPCG) has now performed a GWAS of 2511 (unrelated) familial prostate cancer cases and 1382 unaffected controls from 12 member sites. All samples were genotyped on the Illumina 5M+exome single nucleotide polymorphism (SNP) platform. The GWAS identified a significant evidence for association for SNPs in six regions previously associated with prostate cancer in population-based cohorts, including 3q26.2, 6q25.3, 8q24.21, 10q11.23, 11q13.3, and 17q12. Of note, SNP rs138042437 (p = 1.7e(-8)) at 8q24.21 achieved a large estimated effect size in this cohort (odds ratio = 13.3). 116 previously sampled affected relatives of 62 risk-allele carriers from the GWAS cohort were genotyped for this SNP, identifying 78 additional affected carriers in 62 pedigrees. A test for an excess number of affected carriers among relatives exhibited strong evidence for co-segregation of the variant with disease (p = 8.5e(-11)). The majority (92 %) of risk-allele carriers at rs138042437 had a consistent estimated haplotype spanning approximately 100 kb of 8q24.21 that contained the minor alleles of three rare SNPs (dosage minor allele frequencies <1.7 %), rs183373024 (PRNCR1), previously associated SNP rs188140481, and rs138042437 (CASC19). Strong evidence for co-segregation of a SNP on the haplotype further characterizes the haplotype as a prostate cancer predisposition locus.

FamPipe: An Automatic Analysis Pipeline for Analyzing Sequencing Data in Families for Disease Studies

In disease studies, family-based designs have become an attractive approach to analyzing next-generation sequencing (NGS) data for the identification of rare mutations enriched in families. Substantial research effort has been devoted to developing pipelines for automating sequence alignment, variant calling, and annotation. However, fewer pipelines have been designed specifically for disease studies. Most of the current analysis pipelines for family-based disease studies using NGS data focus on a specific function, such as identifying variants with Mendelian inheritance or identifying shared chromosomal regions among affected family members. Consequently, some other useful family-based analysis tools, such as imputation, linkage, and association tools, have yet to be integrated and automated. We developed FamPipe, a comprehensive analysis pipeline, which includes several family-specific analysis modules, including the identification of shared chromosomal regions among affected family members, prioritizing variants assuming a disease model, imputation of untyped variants, and linkage and association tests. We used simulation studies to compare properties of some modules implemented in FamPipe, and based on the results, we provided suggestions for the selection of modules to achieve an optimal analysis strategy. The pipeline is under the GNU GPL License and can be downloaded for free at http://fampipe.sourceforge.net.

Novel variants in ZNF34 and other brain-expressed transcription factors are shared among early-onset MDD relatives

There are no known genetic variants with large effects on susceptibility to major depressive disorder (MDD). Although one proposed study approach is to increase sensitivity by increasing sample sizes, another is to focus on families with multiple affected individuals to identify genes with rare or novel variants with strong effects. Choosing the family-based approach, we performed whole-exome analysis on affected individuals (n = 12) across five MDD families, each with at least five affected individuals, early onset, and prepubertal diagnoses. We identified 67 genes where novel deleterious variants were shared among affected relatives. Gene ontology analysis shows that of these 67 genes, 18 encode transcriptional regulators, eight of which are expressed in the human brain, including four KRAB-A box-containing Zn(2+) finger repressors. One of these, ZNF34, has been reported as being associated with bipolar disorder and as differentially expressed in bipolar disorder patients compared to healthy controls. We found a novel variant-encoding a non-conservative P17R substitution in the conserved repressor domain of ZNF34 protein-segregating completely with MDD in all available individuals in the family in which it was discovered. Further analysis showed a common ZNF34 coding indel segregating with MDD in a separate family, possibly indicating the presence of an unobserved, linked, rare variant in that particular family. Our results indicate that genes encoding transcription factors expressed in the brain might be an important group of MDD candidate genes and that rare variants in ZNF34 might contribute to susceptibility to MDD and perhaps other affective disorders.

Whole exome sequencing of distant relatives in multiplex families implicates rare variants in candidate genes for oral clefts

A dozen genes/regions have been confirmed as genetic risk factors for oral clefts in human association and linkage studies, and animal models argue even more genes may be involved. Genomic sequencing studies should identify specific causal variants and may reveal additional genes as influencing risk to oral clefts, which have a complex and heterogeneous etiology. We conducted a whole exome sequencing (WES) study to search for potentially causal variants using affected relatives drawn from multiplex cleft families. Two or three affected second, third, and higher degree relatives from 55 multiplex families were sequenced. We examined rare single nucleotide variants (SNVs) shared by affected relatives in 348 recognized candidate genes. Exact probabilities that affected relatives would share these rare variants were calculated, given pedigree structures, and corrected for the number of variants tested. Five novel and potentially damaging SNVs shared by affected distant relatives were found and confirmed by Sanger sequencing. One damaging SNV in CDH1, shared by three affected second cousins from a single family, attained statistical significance (P = 0.02 after correcting for multiple tests). Family-based designs such as the one used in this WES study offer important advantages for identifying genes likely to be causing complex and heterogeneous disorders.