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.

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.

Genome-wide site-specific differential methylation in the blood of individuals with Klinefelter syndrome

Klinefelter syndrome (KS) (47 XXY) is a common sex-chromosome aneuploidy with an estimated prevalence of one in every 660 male births. Investigations into the associations between DNA methylation and the highly variable clinical manifestations of KS have largely focused on the supernumerary X chromosome; systematic investigations of the epigenome have been limited. We obtained genome-wide DNA methylation data from peripheral blood using the Illumina HumanMethylation450K platform in 5 KS (47 XXY) versus 102 male (46 XY) and 113 female (46 XX) control subjects participating in the COPDGene Study. Empirical Bayes-mediated models were used to test for differential methylation by KS status. CpG sites with a false-discovery rate < 0.05 in the discovery cohort which were available on the first-generation HumanMethylation 27 K platform were further examined in an independent replication cohort of 2 KS subjects, 590 male, and 495 female controls drawn from the International COPD Genetics Network (ICGN). Differential methylation at sites throughout the genome were identified, including 86 CpG sites that were differentially methylated in KS subjects relative to both male and female controls. CpG sites annotated to the HEN1 methyltransferase homolog 1 (HENMT1), calcyclin-binding protein (CACYBP), and GTPase-activating protein (SH3 domain)-binding protein 1 (G3BP1) genes were among the "KS-specific" loci that were replicated in ICGN. We conclude that site-specific differential methylation exists throughout the genome in KS. The functional impact and clinical relevance of these differentially methylated loci should be explored in future studies.

Association between genes on chromosome 4p16 and non-syndromic oral clefts in four populations

Isolated cleft lip with or without cleft palate and cleft palate are among the most common human birth defects. Several candidate gene studies on MSX1 have shown significant association between markers in MSX1 and risk of oral clefts, and re-sequencing studies have identified multiple mutations in MSX1 in a small minority of cases, which may account for 1-2% of all isolated oral clefts cases. We explored the 2-Mb region around MSX1, using a marker map of 393 single nucleotide polymorphisms (SNPs) in 297 cleft lip, with or without cleft palate, case-parent trios and 84 cleft palate trios from Maryland, Taiwan, Singapore, and Korea. Both individual markers and haplotypes of two to five SNPs showed several regions yielding statistical evidence for linkage and disequilibrium. Two genes (STK32B and EVC) yielded consistent evidence from cleft lip, with or without cleft palate, trios in all four populations. These two genes plus EVC2 also yielded suggestive evidence for linkage and disequilibrium among cleft palate trios. This analysis suggests that several genes, not just MSX1, in this region may influence risk of oral clefts.

Application of transmission disequilibrium tests to nonsyndromic oral clefts: including candidate genes and environmental exposures in the models

Extensive epidemiological and genetic studies of the cause of oral clefts have demonstrated strong familial aggregation but have failed to yield definitive evidence of any single genetic mechanism. We used the transmission/disequilibrium test (TDT) to investigate the relationship between oral clefts and markers associated with five candidate genes by utilizing 160 parent-offspring trios. Conditional logistic regression models extended the TDT to include covariates as effect modifiers, thus permitting tests for gene-environment interactions. For four of these candidates [transforming growth factor alpha (TGFA), transforming growth factor beta 3 (TGFB3), retinoic acid receptor (RARA), and the proto-oncogene BCL3], we detected modestly elevated odds ratios for the transmission of one marker allele to cleft probands when all the trios were analyzed together. These odds ratios increased when information on type of cleft, race, family history, or maternal smoking were incorporated as effect modifiers. We detected significant interaction between maternal smoking and the transmission of alleles for markers near TGFA and TGFB3; excess transmission of allele 3 at BCL3 was most significant among cleft lip probands; and the odds ratios for transmission of alleles at D19S178 and THRA1 were significant when ethnic group was included in the model. We suggest that utilizing an analytical strategy that allows for stratification of data and incorporating environmental effects into a single analysis may be more effective for detecting genes of small effect.

Resin-process methods for improved isolation of organisms from blood and other body fluids

Two resin-based processes and conventional culturing methods were compared in isolating organisms from blood specimens of patients given antimicrobial agents. Of 20 patients whose specimens were found to contain clinically relevant organisms, one or both of the resin-process methods isolated organisms in eight after conventional culturing was unsuccessful. Conventional cultures were successful in five cases in which the resin-process methods did not work, and both resin-process and conventional culturing were successful in seven cases. In eight cases in which other body fluids were cultured, only a resin-process method was successful in four, although they are intended for use only with blood specimens. Data indicate that resin-process methods should be used along with conventional culturing methods when patients have already been given an antimicrobial. A review of the literature shows that the devices improve the recovery of organisms but does not establish that detection time is decreased.

Enhanced detection of polymicrobic bacteremia by repeat subculture of previously positive blood cultures

Blood subcultures repeated 3 days after the cultures were first identified as positives increased our detection of polymicrobic bacteremia in 9.1 to 27% of clinically significant patient episodes. Reincubation and repeated subculture of previously positive blood cultures had a direct impact on the therapeutic management of patients with polymicrobic bacteremia.