An efficient genotyper and star-allele caller for pharmacogenomics.

High-throughput sequencing provides sufficient means for determining genotypes of clinically important pharmacogenes that can be used to tailor medical decisions to individual patients. However, pharmacogene genotyping, also known as star-allele calling, is a challenging problem that requires accurate copy number calling, structural variation identification, variant calling, and phasing within each pharmacogene copy present in the sample. Here we introduce Aldy 4, a fast and efficient tool for genotyping pharmacogenes that uses combinatorial optimization for accurate star-allele calling across different sequencing technologies. Aldy 4 adds support for long reads and uses a novel phasing model and improved copy number and variant calling models. We compare Aldy 4 against the current state-of-the-art star-allele callers on a large and diverse set of samples and genes sequenced by various sequencing technologies, such as whole-genome and targeted Illumina sequencing, barcoded 10x Genomics, and Pacific Biosciences (PacBio) HiFi. We show that Aldy 4 is the most accurate star-allele caller with near-perfect accuracy in all evaluated contexts, and hope that Aldy remains an invaluable tool in the clinical toolbox even with the advent of long-read sequencing technologies.

Targeted sequencing identifies a missense variant in the BEST3 gene associated with antihypertensive response to hydrochlorothiazide.

Chromosome 12q15 was identified in Genetic Epidemiology of Response Assessment (GERA) and replicated in Pharmacogenomic Evaluation of Antihypertensive Responses (PEAR) for its association with blood pressure (BP) response to hydrochlorothiazide (HCTZ). However, the functional variant is unknown and we aimed to identify the likely functional variants through targeted sequencing. The chromosome 12q15 region was sequenced in 397 best and worst responders to HCTZ in PEAR (N=199) and GERA (N=198) hypertensive study participants. Logistic regression was used for the association analysis adjusting for age, sex, race, and principal components 1 and 2. For validation, the significant single nucleotide polymorphism was tested for association with the change in systolic (ΔSBP) and diastolic BP (ΔDBP) post-treatment in the entire PEAR (N=370) and GERA (N=570) cohorts. A novel missense polymorphism (G>A, Pro383Leu) in BEST3, rs61747221, was significantly associated with better HCTZ response (P=0.0021, odds ratio=2.05). It was validated in the entire cohort of PEAR (ΔSBP: P=0.021, ß=-1.60, ΔDBP: P=0.023, ß=-1.08) and GERA (ΔSBP: P=0.028, ß=-1.95, ΔDBP: P=0.032, ß=-1.28). BEST3 encodes the calcium sensitive chloride channel in the vascular smooth muscle implicated in the regulation of BP, especially in response to vasoconstrictors like angiotensin II. These results suggest that BEST3 is involved in the chronic BP lowering mechanism of thiazides and highlight its importance as a genetic predictor of the BP response to thiazide diuretics.

PGRNseq: a targeted capture sequencing panel for pharmacogenetic research and implementation.

OBJECTIVES: Although the costs associated with whole-genome and whole-exome next-generation sequencing continue to decline, they remain prohibitively expensive for large-scale studies of genetic variation. As an alternative, custom-target sequencing has become a common methodology on the basis of its favorable balance between cost, throughput, and deep coverage. METHODS: We have developed PGRNseq, a custom-capture panel of 84 genes with associations to pharmacogenetic phenotypes, as a tool to explore the relationship between drug response and genetic variation, both common and rare. We utilized a set of 32 diverse HapMap trios and two clinical cohorts to assess platform performance, accuracy, and ability to discover novel variation. RESULTS: We found that PGRNseq generates ultra-deep coverage data (mean=496x) that are over 99.8% concordant with orthogonal datasets. In addition, in our testing sets, PGRNseq identified many novel, rare variants of interest, underscoring its value in both research and clinical settings. CONCLUSION: PGRNseq is an ideal platform for carrying out sequencing-based analyses of pharmacogenetic variation in large cohorts. In addition, the high accuracy associated with genotypes from PGRNseq highlight its utility as a clinical test.

Genome-wide copy number analysis uncovers a new HSCR gene: NRG3.

Hirschsprung disease (HSCR) is a congenital disorder characterized by aganglionosis of the distal intestine. To assess the contribution of copy number variants (CNVs) to HSCR, we analysed the data generated from our previous genome-wide association study on HSCR patients, whereby we identified NRG1 as a new HSCR susceptibility locus. Analysis of 129 Chinese patients and 331 ethnically matched controls showed that HSCR patients have a greater burden of rare CNVs (p = 1.50 × 10(-5)), particularly for those encompassing genes (p = 5.00 × 10(-6)). Our study identified 246 rare-genic CNVs exclusive to patients. Among those, we detected a NRG3 deletion (p = 1.64 × 10(-3)). Subsequent follow-up (96 additional patients and 220 controls) on NRG3 revealed 9 deletions (combined p = 3.36 × 10(-5)) and 2 de novo duplications among patients and two deletions among controls. Importantly, NRG3 is a paralog of NRG1. Stratification of patients by presence/absence of HSCR-associated syndromes showed that while syndromic-HSCR patients carried significantly longer CNVs than the non-syndromic or controls (p = 1.50 × 10(-5)), non-syndromic patients were enriched in CNV number when compared to controls (p = 4.00 × 10(-6)) or the syndromic counterpart. Our results suggest a role for NRG3 in HSCR etiology and provide insights into the relative contribution of structural variants in both syndromic and non-syndromic HSCR. This would be the first genome-wide catalog of copy number variants identified in HSCR.

Carcinoembryonic antigen-related cell adhesion molecule 16 interacts with alpha-tectorin and is mutated in autosomal dominant hearing loss (DFNA4).

We report on a secreted protein found in mammalian cochlear outer hair cells (OHC) that is a member of the carcinoembryonic antigen-related cell adhesion molecule (CEACAM) family of adhesion proteins. Ceacam16 mRNA is expressed in OHC, and its protein product localizes to the tips of the tallest stereocilia and the tectorial membrane (TM). This specific localization suggests a role in maintaining the integrity of the TM as well as in the connection between the OHC stereocilia and TM, a linkage essential for mechanical amplification. In agreement with this role, CEACAM16 colocalizes and coimmunoprecipitates with the TM protein α-tectorin. In addition, we show that mutation of CEACAM16 leads to autosomal dominant nonsyndromic deafness (ADNSHL) at the autosomal dominant hearing loss (DFNA4) locus. In aggregate, these data identify CEACAM16 as an α-tectorin-interacting protein that concentrates at the point of attachment of the TM to the stereocilia and, when mutated, results in ADNSHL at the DFNA4 locus.

Rare APOA5 promoter variants associated with paradoxical HDL cholesterol decrease in response to fenofibric acid therapy.

Individuals with mixed dyslipidemia, including high triglycerides (TGs) and low high density lipoprotein cholesterol (HDL-C), have increased risk for coronary events. We examined the effect of rare genetic variants in the APOA5 gene region on plasma HDL-C, apolipoprotein A-I (apoA-I), and TG response to fenofibric acid monotherapy and in combination with statins. The APOA5 gene region was sequenced in 1,612 individuals with mixed dyslipidemia in a randomized trial of fenofibric acid alone and in combination with statins. Student's t-test and rare variant burden tests were used to examine plasma HDL-C, apoA-I, and TG response. Rare APOA5 promoter region variants were associated with decreased HDL-C and apoA-I levels in response to fenofibric acid therapy; rare missense variants were associated with increased TG response to combination therapy. Further study is needed to examine the effect of these rare variants on coronary outcomes in this population in response to fenofibric acid monotherapy or combined with statins.

Comprehensive genetic testing for hereditary hearing loss using massively parallel sequencing.

The extreme genetic heterogeneity of nonsyndromic hearing loss (NSHL) makes genetic diagnosis expensive and time consuming using available methods. To assess the feasibility of target-enrichment and massively parallel sequencing technologies to interrogate all exons of all genes implicated in NSHL, we tested nine patients diagnosed with hearing loss. Solid-phase (NimbleGen) or solution-based (SureSelect) sequence capture, followed by 454 or Illumina sequencing, respectively, were compared. Sequencing reads were mapped using GSMAPPER, BFAST, and BOWTIE, and pathogenic variants were identified using a custom-variant calling and annotation pipeline (ASAP) that incorporates publicly available in silico pathogenicity prediction tools (SIFT, BLOSUM, Polyphen2, and Align-GVGD). Samples included one negative control, three positive controls (one biological replicate), and six unknowns (10 samples total), in which we genotyped 605 single nucleotide polymorphisms (SNPs) by Sanger sequencing to measure sensitivity and specificity for SureSelect-Illumina and NimbleGen-454 methods at saturating sequence coverage. Causative mutations were identified in the positive controls but not in the negative control. In five of six idiopathic hearing loss patients we identified the pathogenic mutation. Massively parallel sequencing technologies provide sensitivity, specificity, and reproducibility at levels sufficient to perform genetic diagnosis of hearing loss.

Novel Anaplastic Thyroid Cancer PDXs and Cell Lines: Expanding Preclinical Models of Genetic Diversity.

CONTEXT: Anaplastic thyroid cancer (ATC) is a rare, aggressive, and deadly disease. Robust preclinical thyroid cancer models are needed to adequately develop and study novel therapeutic agents. Patient-derived xenograft (PDX) models may resemble patient tumors by recapitulating key genetic alterations and gene expression patterns, making them excellent preclinical models for drug response evaluation. OBJECTIVE: We developed distinct ATC PDX models concurrently with cell lines and characterized them in vitro and in vivo. METHODS: Fresh thyroid tumor from patients with a preoperative diagnosis of ATC was surgically collected and divided for concurrent cell line and PDX model development. Cell lines were created by generating single cells through enzymatic digestion. PDX models were developed following direct subcutaneous implantation of fresh tumor on the flank of immune compromised/athymic mice. RESULTS: Six ATC PDX models and 4 cell lines were developed with distinct genetic profiles. Mutational characterization showed one BRAF/TP53/CDKN2A, one BRAF/CDKN2A, one BRAF/TP53, one TP53 only, one TERT-promoter/HRAS, and one TERT-promoter/KRAS/TP53/NF2/NFE2L2 mutated phenotype. Hematoxylin-eosin staining comparing the PDX models to the original patient surgical specimens show remarkable resemblance, while immunohistochemistry stains for important biomarkers were in full concordance (cytokeratin, TTF-1, PAX8, BRAF). Short tandem repeats DNA fingerprinting analysis of all PDX models and cell lines showed strong concordance with the original tumor. PDX successful establishment rate was 32%. CONCLUSION: We have developed and characterized 6 novel ATC PDX models with 4 matching cell lines. Each PDX model harbors a distinct genetic profile, making them excellent tools for preclinical therapeutic trials.

A High-throughput Approach to Identify Effective Systemic Agents for the Treatment of Anaplastic Thyroid Carcinoma.

BACKGROUND: Despite the use of aggressive multimodality treatment, most anaplastic thyroid carcinoma (ATC) patients die within a year of diagnosis. Although the combination of BRAF and MEK inhibitors has recently been approved for use in BRAF-mutated ATC, they remain effective in a minority of patients who are likely to develop drug resistance. There remains a critical clinical need for effective systemic agents for ATC with a reasonable toxicity profile to allow for rapid translational development. MATERIAL AND METHODS: Twelve human thyroid cancer cell lines with comprehensive genomic characterization were used in a high-throughput screening (HTS) of 257 compounds to select agents with maximal growth inhibition. Cell proliferation, colony formation, orthotopic thyroid models, and patient-derived xenograft (PDX) models were used to validate the selected agents. RESULTS: Seventeen compounds were effective, and docetaxel, LBH-589, and pralatrexate were selected for additional in vitro and in vivo analysis as they have been previously approved by the US Food and Drug Administration for other cancers. Significant tumor growth inhibition (TGI) was detected in all tested models treated with LBH-589; pralatrexate demonstrated significant TGI in the orthotopic papillary thyroid carcinoma model and 2 PDX models; and docetaxel demonstrated significant TGI only in the context of mutant TP53. CONCLUSIONS: HTS identified classes of systemic agents that demonstrate preferential effectiveness against aggressive thyroid cancers, particularly those with mutant TP53. Preclinical validation in both orthotopic and PDX models, which are accurate in vivo models mimicking tumor microenvironment, may support initiation of early-phase clinical trials in non-BRAF mutated or refractory to BRAF/MEK inhibition ATC.

Allelic decomposition and exact genotyping of highly polymorphic and structurally variant genes.

High-throughput sequencing provides the means to determine the allelic decomposition for any gene of interest-the number of copies and the exact sequence content of each copy of a gene. Although many clinically and functionally important genes are highly polymorphic and have undergone structural alterations, no high-throughput sequencing data analysis tool has yet been designed to effectively solve the full allelic decomposition problem. Here we introduce a combinatorial optimization framework that successfully resolves this challenging problem, including for genes with structural alterations. We provide an associated computational tool Aldy that performs allelic decomposition of highly polymorphic, multi-copy genes through using whole or targeted genome sequencing data. For a large diverse sequencing data set, Aldy identifies multiple rare and novel alleles for several important pharmacogenes, significantly improving upon the accuracy and utility of current genotyping assays. As more data sets become available, we expect Aldy to become an essential component of genotyping toolkits.

Frequent amplification of a chr19q13.41 microRNA polycistron in aggressive primitive neuroectodermal brain tumors.

We discovered a high-level amplicon involving the chr19q13.41 microRNA (miRNA) cluster (C19MC) in 11/45 ( approximately 25%) primary CNS-PNET, which results in striking overexpression of miR-517c and 520g. Constitutive expression of miR-517c or 520g promotes in vitro and in vivo oncogenicity, modulates cell survival, and robustly enhances growth of untransformed human neural stem cells (hNSCs) in part by upregulating WNT pathway signaling and restricting differentiation of hNSCs. Remarkably, the C19MC amplicon, which is very rare in other brain tumors (1/263), identifies an aggressive subgroup of CNS-PNET with distinct gene-expression profiles, characteristic histology, and dismal survival. Our data implicate miR-517c and 520g as oncogenes and promising biological markers for CNS-PNET and provide important insights into oncogenic properties of the C19MC locus.

Distribution of polymorphic and non-polymorphic microsatellite repeats in Xenopus tropicalis.

The results of our bioinformatics analysis have found over 91,000 di-, tri-, and tetranucleotide microsatellites in our survey of 25% of the X. tropicalis genome, suggesting there may be over 360,000 within the entire genome. Within the X. tropicalis genome, dinucleotide (78.7%) microsatellites vastly out numbered tri- and tetranucleotide microsatellites. Similarly, AT-rich repeats are overwhelmingly dominant. The four AT-only motifs (AT, AAT, AAAT, and AATT) account for 51,858 out of 91,304 microsatellites found. Individually, AT microsatellites were the most common repeat found, representing over half of all di-, tri-, and tetranucleotide microsatellites. This contrasts with data from other studies, which show that AC is the most frequent microsatellite in vertebrate genomes (Toth et al. 2000). In addition, we have determined the rate of polymorphism for 5,128 non-redundant microsatellites, embedded in unique sequences. Interestingly, this subgroup of microsatellites was determined to have significantly longer repeats than genomic microsatellites as a whole. In addition, microsatellite loci with tandem repeat lengths more than 30 bp exhibited a significantly higher degree of polymorphism than other loci. Pairwise comparisons show that tetranucleotide microsatellites have the highest polymorphic rates. In addition, AAT and ATC showed significant higher polymorphism than other trinucleotide microsatellites, while AGAT and AAAG were significantly more polymorphic than other tetranucleotide microsatellites.

Epigenetic regulation of HYAL-1 hyaluronidase expression. identification of HYAL-1 promoter.

HYAL-1 (hyaluronoglucosaminidase-1) belongs to the hyaluronidase family of enzymes that degrade hyaluronic acid. HYAL-1 is a marker for cancer diagnosis and a molecular determinant of tumor growth, invasion, and angiogenesis. The regulation of HYAL-1 expression is unknown. Real time reverse transcription-PCR using 11 bladder and prostate cancer cells and 69 bladder tissues showed that HYAL-1 mRNA levels are elevated 10-30-fold in cells/tissues that express high hyaluronidase activity. Although multiple transcription start sites (TSS) for HYAL-1 mRNA were detected in various tissues, the major TSS in many tissues, including bladder and prostate, was at nucleotide 27274 in the cosmid clone LUCA13 (AC002455). By analyzing the 1532 base sequence 5' to this TSS, using cloning and luciferase reporter assays, we identified a TACAAA sequence at position -31 and the minimal promoter region between nucleotides -93 and -38. Mutational analysis identified that nucleotides -73 to -50 (which include overlapping binding consensus sites for SP1, Egr-1, and AP-2), bases C(-71) and C(-59), and an NFkappaB-binding site (at position -15) are necessary for promoter activity. The chromatin immunoprecipitation assay identified that Egr-1, AP-2, and NFkappaB bind to the promoter in HYAL-1-expressing cells, whereas SP1 binds to the promoter in non-HYAL-1-expressing cells. 5-Aza-2'-deoxycytidine treatment, bisulfite DNA sequencing, and methylation-specific PCR revealed that HYAL-1 expression is regulated by methylation at C(-71) and C(-59); both Cs are part of the SP1/Egr-1-binding sites. Thus, HYAL-1 expression is epigenetically regulated by the binding of different transcription factors to the methylated and unmethylated HYAL-1 promoter.

MYH11 mutations result in a distinct vascular pathology driven by insulin-like growth factor 1 and angiotensin II.

Non-syndromic thoracic aortic aneurysms and dissections (TAADs) are inherited in an autosomal dominant manner in approximately 20% of cases. Familial TAAD is genetically heterogeneous and four loci have been mapped for this disease to date, including a locus at 16p for TAAD associated with patent ductus arteriosus (PDA). The defective gene at the 16p locus has recently been identified as the smooth muscle cell (SMC)-specific myosin heavy chain gene (MYH11). On sequencing MYH11 in 93 families with TAAD alone and three families with TAAD/PDA, we identified novel mutations in two families with TAAD/PDA, but none in families with TAAD alone. Histopathological analysis of aortic sections from two individuals with MYH11 mutations revealed SMC disarray and focal hyperplasia of SMCs in the aortic media. SMC hyperplasia leading to significant lumen narrowing in some of the vessels of the adventitia was also observed. Insulin-like growth factor-1 (IGF-1) was upregulated in mutant aortas as well as explanted SMCs, but no increase in transforming growth factor-beta expression or downstream targets was observed. Enhanced expression of angiotensin-converting enzyme and markers of Angiotensin II (Ang II) vascular inflammation (macrophage inflammatory protein-1alpha and beta) were also found. These data suggest that MYH11 mutations are likely to be specific to the phenotype of TAAD/PDA and result in a distinct aortic and occlusive vascular pathology potentially driven by IGF-1 and Ang II.

Mutations in smooth muscle alpha-actin (ACTA2) lead to thoracic aortic aneurysms and dissections.

The major function of vascular smooth muscle cells (SMCs) is contraction to regulate blood pressure and flow. SMC contractile force requires cyclic interactions between SMC alpha-actin (encoded by ACTA2) and the beta-myosin heavy chain (encoded by MYH11). Here we show that missense mutations in ACTA2 are responsible for 14% of inherited ascending thoracic aortic aneurysms and dissections (TAAD). Structural analyses and immunofluorescence of actin filaments in SMCs derived from individuals heterozygous for ACTA2 mutations illustrate that these mutations interfere with actin filament assembly and are predicted to decrease SMC contraction. Aortic tissues from affected individuals showed aortic medial degeneration, focal areas of medial SMC hyperplasia and disarray, and stenotic arteries in the vasa vasorum due to medial SMC proliferation. These data, along with the previously reported MYH11 mutations causing familial TAAD, indicate the importance of SMC contraction in maintaining the structural integrity of the ascending aorta.