Comparison of two statistical methods for rare variant association analysis

Studies have shown that rare genetic variants have stronger effects in predisposing common diseases, and several statistical methods have been developed for association studies involving rare variants. In order to better understand how these statistical methods perform, we seek to compare two recently developed rare variant statistical methods (VT and C-alpha) on 10,000 simulated re-sequencing data sets with disease status and the corresponding 10,000 simulated null data sets. The SLC1A1 gene has been suggested to be associated with diastolic blood pressure (DBP) in previous studies. In the current study, we applied VT and C-alpha methods to the empirical re-sequencing data for the SLC1A1 gene from 300 whites and 200 blacks. We found that VT method obtains higher power and performs better than C-alpha method with the simulated data we used. The type I errors were well-controlled for both methods. In addition, both VT and C-alpha methods suggested no statistical evidence for the association between the SLC1A1 gene and DBP. Overall, our findings provided an important comparison of the two statistical methods for future reference and provided preliminary and pioneer findings on the association between the SLC1A1 gene and blood pressure.^

FZD6, MATN2 AND SLC25A32, POSSIBLE CANDIDATE GENES IN NONSYNDROMIC CLEFT LIP AND PALATE

Nonsyndromic cleft lip with or without cleft palate (NSCLP) is a common birth defect with a multifactorial etiology. Despite decades of research, the genetic underpinnings of NSCLP still remain largely unexplained. A genome wide association study (GWAS) of a large NSCLP African American family with seven affected individuals across three generations found evidence for linkage at 8q21.3-24.12 (LOD = 2.98). This region contained three biologically relevant candidate genes: Frizzled-6 (FZD6) (LOD = 2.8), Matrilin-2 (MATN2) (LOD = 2.3), and Solute Carrier Family 25, Member 32 (SLC26A32) (LOD = 1.6). Sequencing of the coding regions and the 5’ and 3’ UTRs of these genes in two affected family members identified a rare intronic variant, rs138557689 (c.-153+432A>C), in FZD6. The rs138557689/C allele segregated with the NSCLP phenotype; in silico analysis predicted and EMSA analysis showed that the 138557689/C allele creates new DNA binding sites. FZD6 is part of the WNT pathway, which is involved in craniofacial development, including midface development and upper lip fusion. Our novel findings suggest that an alteration in FZD6 gene regulation may perturb this tightly controlled biological pathway and in turn contribute to the development of NSCLP in this family. Studies are underway to further define how the rs138557689/C variant affects expression of FZD6.

Altered tumorigenic phenotype in mice with a hypomorphic p53 mutant allele

Mutations in the p53 tumor suppressor gene are found in over 50% of human tumors and in the germline of Li-Fraumeni syndrome families. About 80% of these mutations are missense in nature. In order to study how p53 missense mutations affect tumorigenesis in vivo, we focused on the murine p53 arg-to-his mutation at amino acid 172, which corresponds to the human hot spot mutation at amino acid 175. The double replacement procedure was employed to introduce the p53 R172H mutation into the p53 locus of ES cells and mice were generated. An additional 1bp deletion in the intron 2 splice acceptor site was detected in the same allele in mice. We named this allele p53R172HΔg. This allele makes a small amount of full length p53 mutant protein. ^ Spontaneous tumor formation and survival were studied in these mice. Mice heterozygous for the p53R172HΔg allele showed 50% survival at 17 months of age, similar to the p53+/− mice. Moreover, the p53R172HΔg/+ mice showed a distinct tumor spectrum: 55% sarcomas, including osteosarcoms, fibrosarcomas and angiosarcomas; 27% carcinomas, including lung adenocarcinomas, squamous cell carcinomas, hepatocellular carcinomas and islet cell carcinomas; and 18% lymphomas. Compared to the p53+/− mice, there was a clear increase in the frequency of carcinoma development and a decrease in lymphoma incidence. Among the sarcomas that developed, fibrosarcomas in the skin were also more frequently observed. More importantly, osteosarcomas and carinomas that developed in the p53R172HΔg/+ mice metastasized at very high frequency (64% and 67%, respectively) compared with less than 10% in the p53+/− mice. The metastatic lesions were usually found in lung and liver, and less frequently in other tissues. The altered tumor spectrum in the mice and increased metastatic potential of the tumors suggested that the p53R172H mutation represents a gain-of-function. ^ Mouse embryonic fibroblasts (MEFs) from the mice homozygous and heterozygous for the p53R172HΔg allele were studied for growth characteristics, immortalization potential and genomic instability. All of the p53R172HΔg /+ MEF lines are immortalized under a 3T3 protocol while under the same protocol p53+/− MEFs are not immortalized. Karyotype analysis showed a persistent appearance of chromosome end-to-end fusion in the MEFs both homozygous and heterozygous for the p53R172HΔg allele. These observations suggest that increased genomic instability in the cells may cause the altered tumor phenotypes. ^