Genetic variation among 82 pharmacogenes: The PGRNseq data from the eMERGE network

Genetic variation can affect drug response in multiple ways, although it remains unclear how rare genetic variants affect drug response. The electronic Medical Records and Genomics (eMERGE) Network, collaborating with the Pharmacogenomics Research Network, began eMERGE‐PGx, a targeted sequencing study to assess genetic variation in 82 pharmacogenes critical for implementation of “precision medicine.” The February 2015 eMERGE‐PGx data release includes sequence‐derived data from ∼5,000 clinical subjects. We present the variant frequency spectrum categorized by variant type, ancestry, and predicted function. We found 95.12% of genes have variants with a scaled Combined Annotation‐Dependent Depletion score above 20, and 96.19% of all samples had one or more Clinical Pharmacogenetics Implementation Consortium Level A actionable variants. These data highlight the distribution and scope of genetic variation in relevant pharmacogenes, identifying challenges associated with implementing clinical sequencing for drug treatment at a broader level, underscoring the importance for multifaceted research in the execution of precision medicine.

Genetic variation in the HLA region is associated with susceptibility to herpes zoster

Herpes zoster, commonly referred to as shingles, is caused by the varicella zoster virus (VZV). VZV initially manifests as chicken pox, most commonly in childhood, can remain asymptomatically latent in nerve tissues for many years and often re-emerges as shingles. Although reactivation may be related to immune suppression, aging and female sex, most inter-individual variability in re-emergence risk has not been explained to date. We performed a genome-wide association analyses in 22 981 participants (2280 shingles cases) from the electronic Medical Records and Genomics Network. Using Cox survival and logistic regression, we identified a genomic region in the combined and European ancestry groups that has an age of onset effect reaching genome-wide significance (P>1.0 × 10⁻⁸). This region tags the non-coding gene HCP5 (HLA Complex P5) in the major histocompatibility complex. This gene is an endogenous retrovirus and likely influences viral activity through regulatory functions. Variants in this genetic region are known to be associated with delay in development of AIDS in people infected by HIV. Our study provides further suggestion that this region may have a critical role in viral suppression and could potentially harbor a clinically actionable variant for the shingles vaccine.

Design and anticipated outcomes of the eMERGE-PGx project: a multicenter pilot for preemptive pharmacogenomics in electronic health record systems

We describe here the design and initial implementation of the eMERGE-PGx project. eMERGE-PGx, a partnership of the eMERGE and PGRN consortia, has three objectives : 1) Deploy PGRNseq, a next-generation sequencing platform assessing sequence variation in 84 proposed pharmacogenes, in nearly 9,000 patients likely to be prescribed drugs of interest in a 1–3 year timeframe across several clinical sites; 2) Integrate well-established clinically-validated pharmacogenetic genotypes into the electronic health record with associated clinical decision support and assess process and clinical outcomes of implementation; and 3) Develop a repository of pharmacogenetic variants of unknown significance linked to a repository of EHR-based clinical phenotype data for ongoing pharmacogenomics discovery. We describe site-specific project implementation and anticipated products, including genetic variant and phenotype data repositories, novel variant association studies, clinical decision support modules, clinical and process outcomes, approaches to manage incidental findings, and patient and clinician education methods.

Genome-wide association scan for childhood caries implicates novel genes

Dental caries is the most common chronic disease in children and a major public health concern due to its increasing incidence, serious health and social co-morbidities, and socio-demographic disparities in disease burden. We performed the first genome-wide association scan for dental caries to identify associated genetic loci and nominate candidate genes affecting tooth decay in 1305 US children ages 3-12 yrs. Affection status was defined as 1 or more primary teeth with evidence of decay based on intra-oral examination. No associations met strict criteria for genome-wide significance (p < 10E-7); however, several loci (ACTN2, MTR, and EDARADD, MPPED2, and LPO) with plausible biological roles in dental caries exhibited suggestive evidence for association. Analyses stratified by home fluoride level yielded additional suggestive loci, including TFIP11 in the low-fluoride group, and EPHA7 and ZMPSTE24 in the sufficient-fluoride group. Suggestive loci were tested but not significantly replicated in an independent sample (N = 1695, ages 2-7 yrs) after adjustment for multiple comparisons. This study reinforces the complexity of dental caries, suggesting that numerous loci, mostly having small effects, are involved in cariogenesis. Verification/replication of suggestive loci may highlight biological mechanisms and/or pathways leading to a fuller understanding of the genetic risks for dental caries.

A study of linkage and association of body mass index in the Old Order Amish

Obesity is thought to have a genetic component with the estimates of heritability ranging from 0.25-0.40. As part of an ongoing study of obesity in the Old Order Amish, seven two- and three-generation families (157 individuals) were assessed for 21 traits related to obesity, including body mass index (BMI) and BMI-percentile (a standardized distribution of BMI adjusted for age and sex). Genotyping was performed using a panel of 384 short-tandem repeat markers. In this sample, the estimates of heritability ranged from 0.16-0.31 for BMI and from 0.40-0.52 for BMI-percentile. Model-independent linkage analysis identified candidate regions on chromosomes 1, 5, 7, 8, and 11. Given that several markers on 7q were significant for both BMI and BMI-percentile (P < or = 0.001) and that the structural locus for leptin was located on 7q, this region was considered to be the primary candidate region. Subsequent typing of additional flanking markers on 7q corroborated the original findings. Tests of intrafamilial association for alleles at markers in this candidate region were significant at similar levels. Although there is some evidence for linkage and association in the region containing leptin, there appears to be stronger evidence for linkage (P < or = 0.001) and association (P < or = 0.00001) with BMI in a region 10-15 cM further downstream of leptin, flanked by markers D7S1804 and D7S3070 with peak values from D7S495-D7S1798. Evidence from linkage and association studies suggests that this region (D7S1804-D7S3070) may be responsible, at least in part, for variation in BMI and BMI-percentile in the Old Order Amish.

Autosomal dominant macrothrombocytopenia with leukocyte inclusions (May-Hegglin anomaly) is linked to chromosome 22q12-13

Macrothrombocytopenia with leukocyte inclusions (May-Hegglin anomaly) is a rare autosomal dominant disorder characterized by thrombocytopenia, giant platelets, and Döhle body-like inclusions in leukocytes. To determine the genetic basis of this disorder, we performed a genome-wide screen for linkage in three families with May-Hegglin anomaly. For the pooled analysis of the three families, three markers on chromosome 22 had two-point logarithm-of-difference (lod) scores greater than 3, with a maximum lod score of 3.91 at a recombination fraction (theta) of 0.076 for marker D22S683. Within the largest family (MHA-1), the maximum lod score was 5.36 at theta=0 at marker D22S445. Fine mapping of recombination events using eight adjacent markers indicated that the minimal disease region of family MHA-1 alone is in the approximately 26 cM region from D22S683 to the telomere. The maximum lod score for the three families combined was 5.84 at theta=0 for marker IL2RB. With the assumption of locus homogeneity, haplotype analysis of family MHA-4 indicated the disease region is centromeric to marker D22S1045. These data best support a minimal disease region from D22S683 to D22S1045, a span of about 1 Mb of DNA that contains 17 known genes and 4 predicted genes. Further analysis of this region will identify the genetic basis of May-Hegglin anomaly, facilitating subsequent characterization of the biochemical role of the disease gene in platelet formation.

Cryptic terminal rearrangement of chromosome 22q13.32 detected by FISH in two unrelated patients

Two unrelated patients with cryptic subtelomeric deletions of 22q13.3 were identified using FISH with the commercially available Oncor probe, D22S39. Proband 1 was found to have a derivative chromosome 22 resulting from the unbalanced segregation of a t(1;22)(q44;q13.32) in her mother. Additional FISH analysis of proband 1 and her mother placed the breakpoint on chromosome 22 in this family proximal to D22S55 and D22S39 and distal to D22S45. We have mapped D22S39 to within 170 kb of D22S21 using pulsed field gel electrophoresis. D22S21 is genetically mapped between D22S55 and D22S45. These data indicate that the deletion in proband 1 is smaller than in eight of nine reported del(22)(q13.3) patients. Probands 1 and 2 share features of hypotonia, developmental delay, and expressive language delay, also seen in previously reported del(22)(q13.3) patients, although proband 1 appears to be more mildly affected. Proband 1 is also trisomic for the region 1q44-->qter. This very small duplication has been previously reported only once and the patient had idiopathic mental retardation. This is the first report where 22q13.3 terminal deletion patients have been identified through the use of FISH, and the first report of a deletion of this region occurring because of missegregation of a parental balanced cryptic translocation. We feel that investigation of the frequency of del(22)(q13.3) in the idiopathic mentally retarded population is warranted and may be aided by the ability to use a commercially available probe (D22S39), which is already currently in use in a large number of cytogenetic laboratories.

Segregation of a familial balanced (12;10) insertion resulting in Dup(10)(q21.2q22.1) and Del(10)(q21.2q22.1) in first cousins

An interchromosomal insertion in 3 generations of a family was ascertained through two developmentally delayed first cousins. Cytogenetic analysis using G-banding and chromosome painting showed an apparently balanced direct insertion of chromosome 10 material into chromosome 12, ins(12;10)(q15;q21.2q22.1), in the mothers and grandfather of these children. The proposita inherited only the derivative 10 chromosome, resulting in deletion of 10q21.2 --> 22.1 while her cousin inherited only the derivative 12, resulting in duplication of 10q21.2 --> 22.1. A comparison of the proposita with published deletion cases suggests a pattern of anomalies attributable to deletion of the 10q21 --> q22 region: developmental delay, hypotonia, a heart murmur, telecanthus, broad nasal root and ear abnormalities. This is the first report of a nontandem duplication of the 10q21 --> q22 region. The phenotype of the cousin with the duplication does not overlap greatly with published tandem 10q duplications. Finally, this report reaffirms the importance of obtaining family studies of patients with interstitial chromosomal abnormalities.

Cognitive, behavioral, and neuroanatomical assessment of two unrelated male children expressing FRAXE

Standardized cognitive, behavioral, and neuroanatomical data are presented on 2 unrelated boys with the FRAXE (FMR2) GCC expansion mutation. In the context of normal IQ, both boys had a history of developmental delay, including significant problems with communication, attention, and overactivity. Additionally, one child was diagnosed with autistic disorder. Data from these 2 cases are compared to analogous information from previous reports about individuals with the FRAXE or FRAXA (FMR1) mutation. These comparisons support the idea that FRAXE is associated with nonspecific developmental delay and possibly high-functioning autism.

Two genes required for the binding of an essential Saccharomyces cerevisiae kinetochore complex to DNA

Kinetochores are DNA-protein structures that assemble on centromeric DNA and attach chromosomes to spindle microtubules. Because of their simplicity, the 125-bp centromeres of Saccharomyces cerevisiae are particularly amenable to molecular analysis. Budding yeast centromeres contain three sequence elements of which centromere DNA sequence element III (CDEIII) appears to be particularly important. cis-acting mutations in CDEIII and trans-acting mutations in genes encoding subunits of the CDEIII-binding complex (CBF3) prevent correct chromosome transmission. Using temperature-sensitive mutations in CBF3 subunits, we show a strong correlation between DNA-binding activity measured in vitro and kinetochore activity in vivo. We extend previous findings by Goh and Kilmartin [Goh, P.-Y. & Kilmartin, J.V. (1993) J. Cell Biol. 121, 503-512] to argue that DNA-bound CBF3 may be involved in the operation of a mitotic checkpoint but that functional CBF3 is not required for the assembly of a bipolar spindle.

Identification of essential components of the S. cerevisiae kinetochore

We have designed and utilized two in vivo assays of kinetochore integrity in S. cerevisiae. One assay detects relaxation of a transcription block formed at centromeres; the other detects an increase in the mitotic stability of a dicentric test chromosome. ctf13-30 and ctf14-42 were identified as putative kinetochore mutants by both assays. CTF14 is identical to NDC10/CBF2, a recently identified essential gene that encodes a 110 kd kinetochore component. CTF13 is an essential gene that encodes a predicted 478 amino acid protein with no homology to known proteins. ctf13 mutants missegregate chromosomes at permissive temperature and transiently arrest at nonpermissive temperature as large-budded cells with a G2 DNA content and a short spindle. Antibodies recognizing epitope-tagged CTF13 protein decrease the electrophoretic mobility of a CEN DNA-protein complex formed in vitro. Together, the genetic and biochemical data indicate that CTF13 is an essential kinetochore protein.