Hypomorphic alleles pose challenges in rare disease genomic variant interpretation.

Exon skipping associated with an ATP7B intronic variant in a patient with Wilson's disease. (A) Sashimi plot visualization of aligned RNA sequencing data from proband liver tissue at ATP7B exons 14-13-12. The red track shows traditional RNA-seq data; the blue track shows RNA-seq enriched with exon capture (cDNA-cap) which achieves higher depth of protein-coding transcripts. The histogram indicates overall sequencing depth while arcs tabulate the number of junction-spanning reads supporting exon pairs. (B) The domain structure (top) and exon structure (bottom) of ATP7B. Loss of exon 13 (dashed box) would remove a transmembrane domain and disrupt the first phosphorylation domain.

Gene by stress genome-wide interaction analysis and path analysis identify EBF1 as a cardiovascular and metabolic risk gene.

We performed gene-environment interaction genome-wide association analysis (G × E GWAS) to identify SNPs whose effects on metabolic traits are modified by chronic psychosocial stress in the Multi-Ethnic Study of Atherosclerosis (MESA). In Whites, the G × E GWAS for hip circumference identified five SNPs within the Early B-cell Factor 1 (EBF1) gene, all of which were in strong linkage disequilibrium. The gene-by-stress interaction (SNP × STRESS) term P-values were genome-wide significant (Ps = 7.14E-09 to 2.33E-08, uncorrected; Ps = 1.99E-07 to 5.18E-07, corrected for genomic control). The SNP-only (without interaction) model P-values (Ps = 0.011-0.022) were not significant at the conventional genome-wide significance level. Further analysis of related phenotypes identified gene-by-stress interaction effects for waist circumference, body mass index (BMI), fasting glucose, type II diabetes status, and common carotid intimal-medial thickness (CCIMT), supporting a proposed model of gene-by-stress interaction that connects cardiovascular disease (CVD) risk factor endophenotypes such as central obesity and increased blood glucose or diabetes to CVD itself. Structural equation path analysis suggested that the path from chronic psychosocial stress to CCIMT via hip circumference and fasting glucose was larger (estimate = 0.26, P = 0.033, 95% CI = 0.02-0.49) in the EBF1 rs4704963 CT/CC genotypes group than the same path in the TT group (estimate = 0.004, P = 0.34, 95% CI = -0.004-0.012). We replicated the association of the EBF1 SNPs and hip circumference in the Framingham Offspring Cohort (gene-by-stress term P-values = 0.007-0.012) as well as identified similar path relationships. This observed and replicated interaction between psychosocial stress and variation in the EBF1 gene may provide a biological hypothesis for the complex relationship between psychosocial stress, central obesity, diabetes, and cardiovascular disease.

Age-dependent association between pulmonary tuberculosis and common TOX variants in the 8q12-13 linkage region.

Only a small fraction of individuals infected with Mycobacterium tuberculosis develop clinical tuberculosis (TB) in their lifetime. Genetic epidemiological evidence suggests a genetic determinism of pulmonary TB (PTB), but the molecular basis of genetic predisposition to PTB remains largely unknown. We used a positional-cloning approach to carry out ultrafine linkage-disequilibrium mapping of a previously identified susceptibility locus in chromosomal region 8q12-13 by genotyping 3,216 SNPs in a family-based Moroccan sample including 286 offspring with PTB. We observed 44 PTB-associated SNPs (p < 0.01), which were genotyped in an independent set of 317 cases and 650 controls from Morocco. A single signal, consisting of two correlated SNPs close to TOX, rs1568952 and rs2726600 (combined p = 1.1 × 10(-5) and 9.2 × 10(-5), respectively), was replicated. Stronger evidence of association was found in individuals who developed PTB before the age of 25 years (combined p for rs1568952 = 4.4 × 10(-8); odds ratio of PTB for AA versus AG/GG = 3.09 [1.99-4.78]). The association with rs2726600 (p = 0.04) was subsequently replicated in PTB-affected subjects under 25 years in a study of 243 nuclear families from Madagascar. Stronger evidence of replication in Madagascar was obtained for additional SNPs in strong linkage disequilibrium with the two initial SNPs (p = 0.003 for rs2726597), further confirming the signal. We thus identified around rs1568952 and rs2726600 a cluster of SNPs strongly associated with early-onset PTB in Morocco and Madagascar. SNP rs2726600 is located in a transcription-factor binding site in the 3' region of TOX, and further functional explorations will focus on CD4 T lymphocytes.

Fine mapping of a linkage peak with integration of lipid traits identifies novel coronary artery disease genes on chromosome 5.

BACKGROUND: Coronary artery disease (CAD), and one of its intermediate risk factors, dyslipidemia, possess a demonstrable genetic component, although the genetic architecture is incompletely defined. We previously reported a linkage peak on chromosome 5q31-33 for early-onset CAD where the strength of evidence for linkage was increased in families with higher mean low density lipoprotein-cholesterol (LDL-C). Therefore, we sought to fine-map the peak using association mapping of LDL-C as an intermediate disease-related trait to further define the etiology of this linkage peak. The study populations consisted of 1908 individuals from the CATHGEN biorepository of patients undergoing cardiac catheterization; 254 families (N = 827 individuals) from the GENECARD familial study of early-onset CAD; and 162 aorta samples harvested from deceased donors. Linkage disequilibrium-tagged SNPs were selected with an average of one SNP per 20 kb for 126.6-160.2 MB (region of highest linkage) and less dense spacing (one SNP per 50 kb) for the flanking regions (117.7-126.6 and 160.2-167.5 MB) and genotyped on all samples using a custom Illumina array. Association analysis of each SNP with LDL-C was performed using multivariable linear regression (CATHGEN) and the quantitative trait transmission disequilibrium test (QTDT; GENECARD). SNPs associated with the intermediate quantitative trait, LDL-C, were then assessed for association with CAD (i.e., a qualitative phenotype) using linkage and association in the presence of linkage (APL; GENECARD) and logistic regression (CATHGEN and aortas). RESULTS: We identified four genes with SNPs that showed the strongest and most consistent associations with LDL-C and CAD: EBF1, PPP2R2B, SPOCK1, and PRELID2. The most significant results for association of SNPs with LDL-C were: EBF1, rs6865969, p = 0.01; PPP2R2B, rs2125443, p = 0.005; SPOCK1, rs17600115, p = 0.003; and PRELID2, rs10074645, p = 0.0002). The most significant results for CAD were EBF1, rs6865969, p = 0.007; PPP2R2B, rs7736604, p = 0.0003; SPOCK1, rs17170899, p = 0.004; and PRELID2, rs7713855, p = 0.003. CONCLUSION: Using an intermediate disease-related quantitative trait of LDL-C we have identified four novel CAD genes, EBF1, PRELID2, SPOCK1, and PPP2R2B. These four genes should be further examined in future functional studies as candidate susceptibility loci for cardiovascular disease mediated through LDL-cholesterol pathways.

Accounting for genetic heterogeneity in homozygosity mapping: application to Mendelian susceptibility to mycobacterial disease.

INTRODUCTION: Genome-wide homozygosity mapping is a powerful method for locating rare recessive Mendelian mutations. However, statistical power decreases dramatically in the presence of genetic heterogeneity. METHODS: The authors applied an empirical approach to test for linkage accounting for genetic heterogeneity by calculating the sum of positive per-family multipoint LOD scores (S) across all positions, and obtaining corresponding empirical p values (EmpP) through permutations. RESULTS: The statistical power of the approach was found to be consistently higher than the classical heterogeneity LOD by simulations. Among 21 first-cousin matings with a single affected child, for five families linked to a locus of interest and 16 families to other loci, S/EmpP achieved a power of 40% versus 28% for heterogeneity LOD at an α level of 0.001. The mean size of peak linkage regions was markedly higher for true loci than false positive regions. The S/EmpP approach was applied to a sample of 17 consanguineous families with Mendelian susceptibility to mycobacterial disease, leading to the identification of two mutations in IL12RB1 and TYK2 from the largest of six linkage regions at p<10(-3). CONCLUSIONS: The S/EmpP approach is a flexible and powerful approach that can be applied to linkage analysis of families with suspected Mendelian disorders.

Correlation of intergenerational family sizes suggests a genetic component of reproductive fitness.

Reproductive fitness is a complex phenotype that is a direct measure of Darwinian selection. Estimation of the genetic contribution to this phenotype in human populations is confounded by within-family correlations of sociocultural, economic, and other nongenetic factors that influence family sizes. Here, we report an intergenerational correlation in reproductive success in the Hutterites, a human population that is relatively homogeneous with respect to sociocultural factors that influence fertility. We introduce an estimator of this correlation that takes into account the presence of multiple parent-offspring pairs from the same nuclear family. Statistical significance of the estimated correlation is assessed by a permutation test that maintains the overall structure of the pedigree. Further, temporal trends in fertility within this population are accounted for. Applying these methods to the S-Leut Hutterites yields a correlation in effective family size of 0.29 between couples and their sons and 0.18 between couples and their daughters, with empirical P<1x10-6 and P=.0041, respectively. Similar results were obtained for completed families (0.31 between couples and their sons and 0.23 between couples and their daughters; empirical P<1x10-6 and P=.00059, respectively). We interpret these results as indicating a significant genetic component to reproductive fitness in the Hutterites.