Search

Variation in the expression level and activity of genes involved in drug disposition and action ('pharmacogenes') can affect drug response and toxicity, especially when in tissues of pharmacological importance. Previous studies have relied primarily on microarrays to understand gene expression differences, or have focused on a single tissue or small number of samples. The goal of this study was to use RNA-sequencing (RNA-seq) to determine the expression levels and alternative splicing of 389 Pharmacogenomics Research Network pharmacogenes across four tissues (liver, kidney, heart and adipose) and lymphoblastoid cell lines, which are used widely in pharmacogenomics studies. Analysis of RNA-seq data from 139 different individuals across the 5 tissues (20-45 individuals per tissue type) revealed substantial variation in both expression levels and splicing across samples and tissue types. Comparison with GTEx data yielded a consistent picture. This in-depth exploration also revealed 183 splicing events in pharmacogenes that were previously not annotated. Overall, this study serves as a rich resource for the research community to inform biomarker and drug discovery and use.

Objective: Educational offerings to fill the bioinformatics knowledge gap are a key component to enhancing access and use of health data from the All of Us Research Program. We developed a Train the Trainer-based, innovative training series including project-based learning, modular on-demand demonstrations, and unstructured tutorial time as a model for educational engagement in the All of Us community., Materials and Methods: We highlight our training modules and content, with training survey data informing cycles of development in the creation of a 6-module training series with modular demonstrations., Results: We have conducted 2 public iterations of the Train the Trainer (Tx3) Series based on survey feedback while training over 300 registered researchers to access and analyze data on the All of Us Researcher Workbench., Discussion and Conclusion: Future directions of the Tx3 Series include enhanced focus on project-based learning and learner requests for modularity and asynchronous materials access., (© The Author(s) 2024. Published by Oxford University Press on behalf of the American Medical Informatics Association.)

Despite the large public health toll of smoking, genetic studies of smoking cessation have been limited with few discoveries of risk or protective loci. We investigated common and rare variant associations with success in quitting smoking using a cohort from 8 randomized controlled trials involving 2231 participants and a total of 10,020 common and 24,147 rare variants. We identified 14 novel markers including 6 mapping to genes previously related to psychiatric and substance use disorders, 4 of which were protective (CYP2B6 (rs1175607105), HTR3B (rs1413172952; rs1204720503), rs80210037 on chr15), and 2 of which were associated with reduced cessation (PARP15 (rs2173763), SCL18A2 (rs363222)). The others mapped to areas associated with cancer including FOXP1 (rs1288980) and ZEB1 (rs7349). Network analysis identified significant canonical pathways for the serotonin receptor signaling pathway, nicotine and bupropion metabolism, and several related to tumor suppression. Two novel markers (rs6749438; rs6718083) on chr2 are flanked by genes associated with regulation of bodyweight. The identification of novel loci in this study can provide new targets of pharmacotherapy and inform efforts to develop personalized treatments based on genetic profiles., (© 2024. The Author(s).)

Background: Mosaicism for chromosomal structural abnormalities, other than marker or ring chromosomes, is rarely inherited., Methods: We performed cytogenetics studies and breakpoint analyses on a family with transmission of mosaicism for a derivative chromosome 8 (der(8)), resulting from an unbalanced translocation between the long arms of chromosomes 8 and 21 over three generations., Results: The proband and his maternal half-sister had mosaicism for a der(8) cell line leading to trisomy of the distal 21q, and both had Down syndrome phenotypic features. Mosaicism for a cell line with the der(8) and a normal cell line was also detected in a maternal half-cousin. The der(8) was inherited from the maternal grandmother who had four abnormal cell lines containing the der(8), in addition to a normal cell line. One maternal half-aunt had the der(8) and an isodicentric chromosome 21 (idic(21)). Sequencing studies revealed microhomologies at the junctures of the der(8) and idic(21) in the half-aunt, suggesting a replicative mechanism in the rearrangement formation. Furthermore, interstitial telomeric sequences (ITS) were identified in the juncture between chromosomes 8 and 21 in the der(8)., Conclusion: Mosaicism in the proband, his half-sister and half-cousin resulting from loss of chromosome 21 material from the der(8) appears to be a postzygotic event due to the genomic instability of ITS and associated with selective growth advantage of normal cells. The reversion of the inherited der(8) to a normal chromosome 8 in this family resembles revertant mosaicism of point mutations. We propose that ITS could mediate recurring revertant mosaicism for some constitutional chromosomal structural abnormalities., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2023. No commercial re-use. See rights and permissions. Published by BMJ.)

Phlebotomine sand flies are of global significance as important vectors of human disease, transmitting bacterial, viral, and protozoan pathogens, including the kinetoplastid parasites of the genus Leishmania, the causative agents of devastating diseases collectively termed leishmaniasis. More than 40 pathogenic Leishmania species are transmitted to humans by approximately 35 sand fly species in 98 countries with hundreds of millions of people at risk around the world. No approved efficacious vaccine exists for leishmaniasis and available therapeutic drugs are either toxic and/or expensive, or the parasites are becoming resistant to the more recently developed drugs. Therefore, sand fly and/or reservoir control are currently the most effective strategies to break transmission. To better understand the biology of sand flies, including the mechanisms involved in their vectorial capacity, insecticide resistance, and population structures we sequenced the genomes of two geographically widespread and important sand fly vector species: Phlebotomus papatasi, a vector of Leishmania parasites that cause cutaneous leishmaniasis, (distributed in Europe, the Middle East and North Africa) and Lutzomyia longipalpis, a vector of Leishmania parasites that cause visceral leishmaniasis (distributed across Central and South America). We categorized and curated genes involved in processes important to their roles as disease vectors, including chemosensation, blood feeding, circadian rhythm, immunity, and detoxification, as well as mobile genetic elements. We also defined gene orthology and observed micro-synteny among the genomes. Finally, we present the genetic diversity and population structure of these species in their respective geographical areas. These genomes will be a foundation on which to base future efforts to prevent vector-borne transmission of Leishmania parasites., Competing Interests: The authors have declared that no competing interests exist., (Copyright: This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.)

Background: Children with persistent, isolated microscopic hematuria typically undergo a limited diagnostic workup and are monitored for signs of kidney disease in long-term longitudinal follow-up, which can delay diagnosis and allow disease progression in some cases., Methods: To determine the clinical utility of genetic screening in this population, we performed targeted genetic testing using a custom, 32-gene next-generation sequencing panel for progressive kidney disease on children referred to the Texas Children's Hospital Pediatric Nephrology clinic for persistent, microscopic hematuria (n = 30; cohort 1). Patients with microscopic hematuria identified by urinalysis on at least two separate occasions were eligible for enrollment, but those with other evidence of kidney disease were excluded. Results were analyzed for sequence variants using the American College of Medical Genetics and Genomics (ACMG) guideline for data interpretation and were validated using a secondary analysis of a dataset of children with hematuria and normal kidney function who had undergone genetic testing as part of an industry-sponsored program (cohort 2; n = 67)., Results: In cohort 1 33% of subjects (10/30) had pathogenic or likely pathogenic (P/LP) variants in the type IV collagen genes (COL4A3/A4/A5), and 10% (3/30) had variants of uncertain significance in these genes. The high diagnostic rate in type IV collagen genes was confirmed in cohort 2, where 27% (18/67) of subjects had P/LP variants in COL4A3/A4/A5 genes., Conclusions: Children with persistent, isolated microscopic hematuria have a high likelihood of having pathogenic variants in type IV collagen genes and genetic screening should be considered. A higher resolution version of the Graphical abstract is available as Supplementary information., (© 2022. The Author(s), under exclusive licence to International Pediatric Nephrology Association.)

Objectives: We evaluated multiple genotyping/sequencing approaches in a homologous region of chromosome 19, and investigated associations of two common 3'-UTR CYP2A6 variants with activity in vivo., Methods: Individuals (n = 1704) of European and African ancestry were phenotyped for the nicotine metabolite ratio (NMR), an index of CYP2A6 activity, and genotyped/sequenced using deep amplicon exon sequencing, SNP array, genotype imputation and targeted capture sequencing. Amplicon exon sequencing was the gold standard to which other methods were compared within-individual for CYP2A6, CYP2A7, CYP2A13, and CYP2B6 exons to identify highly discordant positions. Linear regression models evaluated the association of CYP2A6*1B and rs8192733 genotypes (coded additively) with logNMR., Results: All approaches were ≤2.6% discordant with the gold standard; discordant calls were concentrated at few positions. Fifteen positions were discordant in >10% of individuals, with 12 appearing in regions of high identity between homologous genes (e.g. CYP2A6 and CYP2A7). For six, allele frequencies in our study and online databases were discrepant, suggesting errors in online sources. In the European-ancestry group (n = 935), CYP2A6*1B and rs8192733 were associated with logNMR (P < 0.001). A combined model found main effects of both variants on increasing logNMR. Similar trends were found in those of African ancestry (n = 506)., Conclusion: Multiple genotyping/sequencing approaches used in this chromosome 19 region contain genotyping/sequencing errors, as do online databases. Gene-specific primers and SNP array probes must consider gene homology; short-read sequencing of related genes in a single reaction should be avoided. Using improved sequencing approaches, we characterized two gain-of-function 3'-UTR variants, including the relatively understudied rs8192733., (Copyright © 2022 Wolters Kluwer Health, Inc. All rights reserved.)

Purpose: The Mayo-Baylor RIGHT 10K Study enabled preemptive, sequence-based pharmacogenomics (PGx)-driven drug prescribing practices in routine clinical care within a large cohort. We also generated the tools and resources necessary for clinical PGx implementation and identified challenges that need to be overcome. Furthermore, we measured the frequency of both common genetic variation for which clinical guidelines already exist and rare variation that could be detected by DNA sequencing, rather than genotyping., Methods: Targeted oligonucleotide-capture sequencing of 77 pharmacogenes was performed using DNA from 10,077 consented Mayo Clinic Biobank volunteers. The resulting predicted drug response-related phenotypes for 13 genes, including CYP2D6 and HLA, affecting 21 drug-gene pairs, were deposited preemptively in the Mayo electronic health record., Results: For the 13 pharmacogenes of interest, the genomes of 79% of participants carried clinically actionable variants in 3 or more genes, and DNA sequencing identified an average of 3.3 additional conservatively predicted deleterious variants that would not have been evident using genotyping., Conclusion: Implementation of preemptive rather than reactive and sequence-based rather than genotype-based PGx prescribing revealed nearly universal patient applicability and required integrated institution-wide resources to fully realize individualized drug therapy and to show more efficient use of health care resources., Competing Interests: Conflict of Interest Liewei Wang, John Logan Black III, and Richard M. Weinshilboum are cofounders of and stockholders in OneOme, LLC, which was used only to return results to the study participants. Additionally, John Logan Black III and Mayo Clinic Ventures have applied for a patent on the CNVAR software cited in this study as well as the methodology upon which the software is based. All other authors declare no conflicts of interest., (Copyright © 2022 American College of Medical Genetics and Genomics. Published by Elsevier Inc. All rights reserved.)

Pharmacogenetic tests typically target selected sequence variants to identify haplotypes that are often defined by star (∗) allele nomenclature. Due to their design, these targeted genotyping assays are unable to detect novel variants that may change the function of the gene product and thereby affect phenotype prediction and patient care. In the current study, 137 DNA samples that were previously characterized by the Genetic Testing Reference Material (GeT-RM) program using a variety of targeted genotyping methods were recharacterized using targeted and whole genome sequencing analysis. Sequence data were analyzed using three genotype calling tools to identify star allele diplotypes for CYP2C8, CYP2C9, and CYP2C19. The genotype calls from next-generation sequencing (NGS) correlated well to those previously reported, except when novel alleles were present in a sample. Six novel alleles and 38 novel suballeles were identified in the three genes due to identification of variants not covered by targeted genotyping assays. In addition, several ambiguous genotype calls from a previous study were resolved using the NGS and/or long-read NGS data. Diplotype calls were mostly consistent between the calling algorithms, although several discrepancies were noted. This study highlights the utility of NGS for pharmacogenetic testing and demonstrates that there are many novel alleles that are yet to be discovered, even in highly characterized genes such as CYP2C9 and CYP2C19., (Copyright © 2022 Association for Molecular Pathology and American Society for Investigative Pathology. All rights reserved.)

Clinical pharmacogenomic testing typically uses targeted genotyping, which only detects variants included in the test design and may vary among laboratories. To evaluate the potential patient impact of genotyping compared with sequencing, which can detect common and rare variants, an in silico targeted genotyping panel was developed based on the variants most commonly included in clinical tests and applied to a cohort of 10,030 participants who underwent sequencing for CYP1A2, CYP2C19, CYP2C9, CYP2D6, CYP3A4, CYP3A5, DPYD, SLCO1B1, TPMT, UGT1A1, and VKORC1. The results of in silico targeted genotyping were compared with the clinically reported sequencing results. Of the 10,030 participants, 2780 (28%) had at least one potentially clinically relevant variant/allele identified by sequencing that would not have been detected in a standard targeted genotyping panel. The genes with the largest number of participants with variants only detected by sequencing were SLCO1B1, DPYD, and CYP2D6, which affected 13%, 6.3%, and 3.5% of participants, respectively. DPYD (112 variants) and CYP2D6 (103 variants) had the largest number of unique variants detected only by sequencing. Although targeted genotyping detects most clinically significant pharmacogenomic variants, sequencing-based approaches are necessary to detect rare variants that collectively affect many patients. However, efforts to establish pharmacogenomic variant classification systems and nomenclature to accommodate rare variants will be required to adopt sequencing-based pharmacogenomics., (Copyright © 2022 Association for Molecular Pathology and American Society for Investigative Pathology. Published by Elsevier Inc. All rights reserved.)

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Endocrine Society. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Endocrine Society. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Introduction: Alaska Native and American Indian (AN/AI) populations have higher tobacco use prevalence than other ethnic/racial groups. Pharmacogenetic testing to tailor tobacco cessation treatment may improve cessation rates. This study characterized polymorphic variations among AN/AI people in genes associated with metabolism of nicotine and drugs used for tobacco cessation., Methods: Recruitment of AN/AI individuals represented six subgroups, five geographic subgroups throughout Alaska and a subgroup comprised of AIs from the lower 48 states living in Alaska. We sequenced the CYP2A6 and CYP2B6 genes to identify known and novel gain, reduced, and loss-of-function alleles, including structural variation (eg, gene deletions, duplications, and hybridizations)., Results: Variant allele frequencies differed substantially between AN/AI subgroups. The gene deletion CYP2A6*4 and reduced function CYP2A6*9 alleles were found at high frequency in Northern/Western subgroups and in Lower 48/Interior subgroups, respectively. The reduced function CYP2B6*6 allele was observed in all subgroups and a novel, predicted reduced function CYP2B6 variant was found at relatively high frequency in the Southeastern subgroup., Conclusions: Diverse CYP2A6 and CYP2B6 variation among the subgroups highlight the need for comprehensive pharmacogenetic testing to guide tobacco cessation therapy for AN/AI populations., Implications: Nicotine metabolism is largely determined by CYP2A6 genotype, and variation in CYP2A6 activity has altered the treatment success in other populations. These findings suggest pharmacogenetic-guided smoking cessation drug treatment could provide benefit to this unique population seeking tobacco cessation therapy., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Research on Nicotine and Tobacco. All rights reserved.For permissions, please e-mail: journals.permissions@oup.com.)

CYP2C9 and CYP2C19 are highly polymorphic pharmacogenes; however, clinically actionable genetic variability in drug metabolism due to these genes has been limited to a few common alleles. The identification and functional characterization of less-common open reading frame sequence variation might help to individualize therapy with drugs that are substrates for the enzymes encoded by these genes. The present study identified seven uncharacterized variants each in CYP2C9 and CYP2C19 using next-generation sequence data for 1013 subjects, and functionally characterized the encoded proteins. Constructs were created and transiently expressed in COS-1 cells for the assay of protein concentration and enzyme activities using fluorometric substrates and liquid chromatography- tandem mass spectrometry with tolbutamide ( CYP2C9 ) and (S)-mephenytoin ( CYP2C19 ) as prototypic substrates. The results were compared with the SIFT, Polyphen, and Provean functional prediction software programs. Cytochrome P450 oxidoreductase (CPR) activities were also determined. Positive correlations were observed between protein content and fluorometric enzyme activity for variants of CYP2C9 ( P < 0.05) and CYP2C19 ( P < 0.0005). However, CYP2C9 709G>C and CYP2C19 65A>G activities were much lower than predicted based on protein content. Substrate intrinsic clearance values for CYP2C9 218C>T, 343A>C, and CYP2C19 337G>A, 518C>T, 556C>T, and 557G>A were less than 25% of wild-type allozymes. CPR activity levels were similar for all variants. In summary, sequencing of CYP2C9 and CYP2C19 in 1013 subjects identified low-frequency variants that had not previously been functionally characterized. In silico predictions were not always consistent with functional assay results. These observations emphasize the need for high-throughput methods for pharmacogene variant mutagenesis and functional characterization., (Copyright © 2019 by The American Society for Pharmacology and Experimental Therapeutics.)

We used whole-organ mapping to study the locoregional molecular changes in a human bladder containing multifocal cancer. Widespread DNA methylation changes were identified in the entire mucosa, representing the initial field effect. The field effect was associated with subclonal low-allele frequency mutations and a small number of DNA copy alterations. A founder mutation in the RNA splicing gene, ACIN1, was identified in normal mucosa and expanded clonally with an additional 21 mutations in progression to carcinoma. The patterns of mutations and copy number changes in carcinoma in situ and foci of carcinoma were almost identical, confirming their clonal origins. The pathways affected by the DNA copy alterations and mutations, including the Kras pathway, were preceded by the field changes in DNA methylation, suggesting that they reinforced mechanisms that had already been initiated by methylation. The results demonstrate that DNA methylation can serve as the initiator of bladder carcinogenesis., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

Background: Recently, 34 genes had been associated with differential expression relative to blood pressure (BP)/ hypertension (HTN). We hypothesize that some of the genes associated with BP/HTN are also associated with BP response to antihypertensive treatment with thiazide diuretics., Methods: We assessed these 34 genes for association with differential expression to BP response to thiazide diuretics with RNA sequencing in whole blood samples from 150 hypertensive participants from the Pharmacogenomic Evaluation of Antihypertensive Responses (PEAR) and PEAR-2 studies. PEAR white and PEAR-2 white and black participants (n = 50 for each group) were selected based on the upper and lower quartile of BP response to hydrochlorothiazide (HCTZ) and to chlorthalidone., Results: FOS, DUSP1 and PPP1R15A were differentially expressed across all cohorts (meta-analysis p-value < 2.0 × 10 - 6 ), and responders to HCTZ or chlorthalidone presented up-regulated transcripts. Rs11065987 in chromosome 12, a trans-eQTL for expression of FOS, PPP1R15A and other genes, is also associated with BP response to HCTZ in PEAR whites (SBP: β = - 2.1; p = 1.7 × 10 - 3 ; DBP: β = - 1.4; p = 2.9 × 10 - 3 )., Conclusions: These findings suggest FOS, DUSP1 and PPP1R15A as potential molecular determinants of antihypertensive response to thiazide diuretics., Trial Registration: NCT00246519 , NCT01203852 www.clinicaltrials.gov.

Background: Plasma renin is an important regulator of blood pressure (BP). Plasma renin activity (PRA) has been shown to correlate with variability in BP response to antihypertensive agents. We conducted a genome-wide association study to identify single-nucleotide polymorphisms (SNPs) associated with baseline PRA using data from the PEAR study (Pharmacogenomic Evaluation of Antihypertensive Responses)., Methods: Multiple linear regression analysis was performed in 461 whites and 297 blacks using an additive model, adjusting for age, sex, and ancestry-specific principal components. Top SNPs were prioritized by testing the expected direction of association for BP response to atenolol and hydrochlorothiazide. Top regions from the BP response prioritization were tested for functional evidence through differences in gene expression by genotype using RNA sequencing data. Regions with functional evidence were assessed for replication with baseline PRA in an independent study (PEAR-2)., Results: Our top SNP rs3784921 was in the SNN-TXNDC11 gene region. The G allele of rs3784921 was associated with higher baseline PRA (β=0.47; P =2.09×10 -6 ) and smaller systolic BP reduction in response to hydrochlorothiazide (β=2.97; 1-sided P =0.006). In addition, TXNDC11 expression differed by rs3784921 genotype ( P =0.007), and rs1802409, a proxy SNP for rs3784921 ( r 2 =0.98-1.00), replicated in PEAR-2 (β=0.15; 1-sided P =0.038). Additional SNPs associated with baseline PRA that passed BP response prioritization were in/near the genes CHD9, XIRP2, and GHR. CONCLUSIONS: We identified multiple regions associated with baseline PRA that were prioritized through BP response signals to 2 mechanistically different antihypertensive drugs., Clinical Trial Registration: URL: https://www.clinicaltrials.gov. Unique identifier: NCT00246519., (© 2018 American Heart Association, Inc.)

Objectives: Smoking patterns and cessation rates vary widely across smokers and can be influenced by variation in rates of nicotine metabolism [i.e. cytochrome P450 2A6 (CYP2A6), enzyme activity]. There is high heritability of CYP2A6-mediated nicotine metabolism (60-80%) owing to known and unidentified genetic variation in the CYP2A6 gene. We aimed to identify and characterize additional genetic variants at the CYP2A6 gene locus., Methods: A new CYP2A6-specific sequencing method was used to investigate genetic variation in CYP2A6. Novel variants were characterized in a White human liver bank that has been extensively phenotyped for CYP2A6. Linkage and haplotype structure for the novel single nucleotide polymorphisms (SNPs) were assessed. The association between novel five-SNP diplotypes and nicotine metabolism rate was investigated., Results: Seven high-frequency (minor allele frequencies ≥6%) noncoding SNPs were identified as important contributors to CYP2A6 phenotypes in a White human liver bank (rs57837628, rs7260629, rs7259706, rs150298687 (also denoted rs4803381), rs56113850, rs28399453, and rs8192733), accounting for two times more variation in in-vitro CYP2A6 activity relative to the four established functional CYP2A6 variants that are frequently tested in Whites (CYP2A6*2, *4, *9, and *12). Two pairs of novel SNPs were in high linkage disequilibrium, allowing us to establish five-SNP diplotypes that were associated with CYP2A6 enzyme activity (rate of nicotine metabolism) in-vitro in the liver bank and in-vivo among smokers., Conclusion: The novel five-SNP diplotype may be useful to incorporate into CYP2A6 genotype models for personalized prediction of nicotine metabolism rate, cessation success, and response to pharmacotherapies.

Thiazide diuretics (TD) are commonly prescribed anti-hypertensives worldwide. However, <40% of patients treated with thiazide monotherapy achieve BP control. This study uses whole transcriptome sequencing to identify novel molecular markers associated with BP response to TD. We assessed global RNA expression levels in whole blood samples from 150 participants, representing patients in the upper and lower quartile of BP response to TD from the Pharmacogenomic Evaluation of Antihypertensive Responses (PEAR) (50 whites) and from PEAR-2 (50 whites and 50 blacks). In each study cohort, we performed poly-A RNA-sequencing in baseline samples from 25 responders and 25 non-responders to hydrochlorothiazide (HCTZ) or chlorthalidone. At FDR adjusted p-value < 0.05, 29 genes were differentially expressed in relation to HCTZ or chlorthalidone BP response in whites. For each differentially expressed gene, replication was attempted in the alternate white group and PEAR-2 blacks. CEBPD (meta-analysis p = 1.8 × 10 -11 ) and TSC22D3 (p = 1.9 × 10 -9 ) were differentially expressed in all 3 cohorts, and explain, in aggregate, 21.9% of response variability to TD. This is the first report of the use of transcriptome-wide sequencing data to identify molecular markers of antihypertensive drug response. These findings support CEBPD and TSC22D3 as potential biomarkers of BP response to TD.

Our previous GWAS using samples from the NSABP P-1 and P-2 selective estrogen receptor modulator (SERM) breast cancer prevention trials identified SNPs in ZNF423 and near CTSO that were associated with breast cancer risk during SERM chemoprevention. We have now performed Next Generation DNA sequencing to identify additional SNPs that might contribute to breast cancer risk and to extend our observation that SNPs located hundreds of bp from estrogen response elements (EREs) can alter estrogen receptor alpha (ERα) binding in a SERM-dependent fashion. Our study utilized a nested case-control cohort selected from patients enrolled in the original GWAS, with 199 cases who developed breast cancer during SERM therapy and 201 matched controls who did not. We resequenced approximately 500 kb across both ZNF423 and CTSO , followed by functional genomic studies. We identified 4079 SNPs across ZNF423 and 3876 across CTSO , with 9 SNPs in ZNF423 and 12 in CTSO with p  < 1E-02 that were within 500 bp of an ERE motif. The rs746157 ( p  = 8.44E-04) and rs12918288 SNPs ( p  = 3.43E-03) in intron 5 of ZNF423 , were in linkage equilibrium and were associated with alterations in ER-binding to an ERE motif distant from these SNPs. We also studied all nonsynonymous SNPs in both genes and observed that one nsSNP in ZNF423 displayed decreased protein expression. In conclusion, we identified additional functional SNPs in ZNF423 that were associated with SNP and SERM-dependent alternations in ER binding and transcriptional regulation for an ERE at a distance from the SNPs, thus providing novel insight into mechanisms of SERM effect.

Background: Helicoverpa armigera and Helicoverpa zea are major caterpillar pests of Old and New World agriculture, respectively. Both, particularly H. armigera, are extremely polyphagous, and H. armigera has developed resistance to many insecticides. Here we use comparative genomics, transcriptomics and resequencing to elucidate the genetic basis for their properties as pests., Results: We find that, prior to their divergence about 1.5 Mya, the H. armigera/H. zea lineage had accumulated up to more than 100 more members of specific detoxification and digestion gene families and more than 100 extra gustatory receptor genes, compared to other lepidopterans with narrower host ranges. The two genomes remain very similar in gene content and order, but H. armigera is more polymorphic overall, and H. zea has lost several detoxification genes, as well as about 50 gustatory receptor genes. It also lacks certain genes and alleles conferring insecticide resistance found in H. armigera. Non-synonymous sites in the expanded gene families above are rapidly diverging, both between paralogues and between orthologues in the two species. Whole genome transcriptomic analyses of H. armigera larvae show widely divergent responses to different host plants, including responses among many of the duplicated detoxification and digestion genes., Conclusions: The extreme polyphagy of the two heliothines is associated with extensive amplification and neofunctionalisation of genes involved in host finding and use, coupled with versatile transcriptional responses on different hosts. H. armigera's invasion of the Americas in recent years means that hybridisation could generate populations that are both locally adapted and insecticide resistant.

Background: Thiazide diuretics are among the most commonly prescribed antihypertensives. However, <50% of thiazide-treated patients achieve blood pressure (BP) control. Herein, we used different omics (genomics and transcriptomics) to identify novel biomarkers of thiazide diuretics BP response., Methods and Results: Genome-wide analysis included 228 white hypertensives with BP determined at baseline and after 9 weeks of hydrochlorothiazide. Single-nucleotide polymorphisms with P <5×10 - 5 were prioritized according to their biological function, using RegulomeDB, haploreg, and Genome-Wide Annotation of Variants. The results from the prioritization approach revealed rs10995 as the most likely functional single-nucleotide polymorphism, among single-nucleotide polymorphisms tested, that has been associated with hydrochlorothiazide BP response. The rs10995 G-allele was associated with better BP response to hydrochlorothiazide versus noncarriers (Δ systolic BP/Δ diastolic BP: -12.3/-8.2 versus -6.8/-3.5 mm Hg, respectively, Δ systolic BP P=3×10 -4 , Δ diastolic BP P=5×10 -5 ). This association was replicated in independent participants treated with chlorthalidone. In addition, rs10995 G-allele was associated with increased mRNA expression of VASP (vasodilator-stimulated phosphoprotein). Moreover, baseline expression of the VASP mRNA was significantly higher in 25 good responders to hydrochlorothiazide compared with 25 poor responders (P=0.01). This finding was replicated in independent participants treated with chlorthalidone (P=0.04). Last, allelic-specific expression analysis revealed a significant but modest imbalance with rs10995 and rs10156, a single-nucleotide polymorphism in high linkage disequilibrium (r 2 =0.7) with rs10995, which both could contribute to the observed genetic effects by affecting VASP mRNA expression., Conclusions: This study highlights the strength of using different omics to identify novel biomarkers of drug response and suggests VASP as a potential determinant of thiazide diuretics BP response., Clinical Trial Registration: URL: http://www.clinicaltrials.gov. Unique identifiers: NCT00246519 and NCT01203852., (© 2017 American Heart Association, Inc.)

Manduca sexta, known as the tobacco hornworm or Carolina sphinx moth, is a lepidopteran insect that is used extensively as a model system for research in insect biochemistry, physiology, neurobiology, development, and immunity. One important benefit of this species as an experimental model is its extremely large size, reaching more than 10 g in the larval stage. M. sexta larvae feed on solanaceous plants and thus must tolerate a substantial challenge from plant allelochemicals, including nicotine. We report the sequence and annotation of the M. sexta genome, and a survey of gene expression in various tissues and developmental stages. The Msex&#95;1.0 genome assembly resulted in a total genome size of 419.4 Mbp. Repetitive sequences accounted for 25.8% of the assembled genome. The official gene set is comprised of 15,451 protein-coding genes, of which 2498 were manually curated. Extensive RNA-seq data from many tissues and developmental stages were used to improve gene models and for insights into gene expression patterns. Genome wide synteny analysis indicated a high level of macrosynteny in the Lepidoptera. Annotation and analyses were carried out for gene families involved in a wide spectrum of biological processes, including apoptosis, vacuole sorting, growth and development, structures of exoskeleton, egg shells, and muscle, vision, chemosensation, ion channels, signal transduction, neuropeptide signaling, neurotransmitter synthesis and transport, nicotine tolerance, lipid metabolism, and immunity. This genome sequence, annotation, and analysis provide an important new resource from a well-studied model insect species and will facilitate further biochemical and mechanistic experimental studies of many biological systems in insects., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

Each year in the United States, thousands of cases of sudden and unexpected deaths of infants, children, and young adults are assigned an undetermined cause of death after postmortem investigation and autopsy. Heritable genetic variants have been suggested as the cause of up to a third of sudden death (SD) cases. Elucidation of the genetic variants involved in SD cases is important to not only help establish cause and manner of death of these individuals, but to also aid in determining whether familial genetic testing should be considered. Previously, these types of postmortem screenings have not been a feasible option for most county medical examiners' and coroners' offices. We sequenced full exons of 64 genes associated with SD in the largest known cohort (351) of infant and young SD decedents using massively parallel sequencing at <$600 per sample. Genetic variants were assessed through literature review and clinical evaluation by a multidisciplinary consortium of experts. Thirteen individuals (3.7%), eight infants (2.8% of those <1 yr of age) and five children/young adults (7.0% of those >1 yr of age), were found to have a reportable genetic variant contributing to SD. These percentages represent an estimate lower than those previously reported. Overall yields and results likely vary between studies due to differences in evaluation techniques and reporting. Additionally, we recommend ongoing assessment of data, including nonreported novel variants, as technology and literature continually advance. This study demonstrates a strategy to implement molecular autopsies in medicolegal investigations of young SD decedents., (© 2016 Methner et al.; Published by Cold Spring Harbor Laboratory Press.)

We performed RNA sequencing (RNAseq) on peripheral blood mononuclear cells (PBMCs) to identify differentially expressed gene transcripts (DEGs) after kidney transplantation and after the start of immunosuppressive drugs. RNAseq is superior to microarray to determine DEGs because it's not limited to available probes, has increased sensitivity, and detects alternative and previously unknown transcripts. DEGs were determined in 32 adult kidney recipients, without clinical acute rejection (AR), treated with antibody induction, calcineurin inhibitor, mycophenolate, with and without steroids. Blood was obtained pre-transplant (baseline), week 1, months 3 and 6 post-transplant. PBMCs were isolated, RNA extracted and gene expression measured using RNAseq. Principal components (PCs) were computed using a surrogate variable approach. DEGs post-transplant were identified by controlling false discovery rate (FDR) at < 0.01 with at least a 2 fold change in expression from pre-transplant. The top 5 DEGs with higher levels of transcripts in blood at week 1 were TOMM40L, TMEM205, OLFM4, MMP8, and OSBPL9 compared to baseline. The top 5 DEGs with lower levels at week 1 post-transplant were IL7R, KLRC3, CD3E, CD3D, and KLRC2 (Striking Image) compared to baseline. The top pathways from genes with lower levels at 1 week post-transplant compared to baseline, were T cell receptor signaling and iCOS-iCOSL signaling while the top pathways from genes with higher levels than baseline were axonal guidance signaling and LXR/RXR activation. Gene expression signatures at month 3 were similar to week 1. DEGs at 6 months post-transplant create a different gene signature than week 1 or month 3 post-transplant. RNAseq analysis identified more DEGs with lower than higher levels in blood compared to baseline at week 1 and month 3. The number of DEGs decreased with time post-transplant. Further investigations to determine the specific lymphocyte(s) responsible for differential gene expression may be important in selecting and personalizing immune suppressant drugs and may lead to targeted therapies.

Background: The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats., Results: We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits., Conclusions: These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

Gall-forming arthropods are highly specialized herbivores that, in combination with their hosts, produce extended phenotypes with unique morphologies [1]. Many are economically important, and others have improved our understanding of ecology and adaptive radiation [2]. However, the mechanisms that these arthropods use to induce plant galls are poorly understood. We sequenced the genome of the Hessian fly (Mayetiola destructor; Diptera: Cecidomyiidae), a plant parasitic gall midge and a pest of wheat (Triticum spp.), with the aim of identifying genic modifications that contribute to its plant-parasitic lifestyle. Among several adaptive modifications, we discovered an expansive reservoir of potential effector proteins. Nearly 5% of the 20,163 predicted gene models matched putative effector gene transcripts present in the M. destructor larval salivary gland. Another 466 putative effectors were discovered among the genes that have no sequence similarities in other organisms. The largest known arthropod gene family (family SSGP-71) was also discovered within the effector reservoir. SSGP-71 proteins lack sequence homologies to other proteins, but their structures resemble both ubiquitin E3 ligases in plants and E3-ligase-mimicking effectors in plant pathogenic bacteria. SSGP-71 proteins and wheat Skp proteins interact in vivo. Mutations in different SSGP-71 genes avoid the effector-triggered immunity that is directed by the wheat resistance genes H6 and H9. Results point to effectors as the agents responsible for arthropod-induced plant gall formation., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history., Competing Interests: The authors have declared that no competing interests exist.

We describe here the design and initial implementation of the eMERGE-PGx project. eMERGE-PGx, a partnership of the Electronic Medical Records and Genomics Network and the Pharmacogenomics Research Network, has three objectives: (i) to 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- to 3-year time frame across several clinical sites; (ii) to integrate well-established clinically validated pharmacogenetic genotypes into the electronic health record with associated clinical decision support and to assess process and clinical outcomes of implementation; and (iii) to develop a repository of pharmacogenetic variants of unknown significance linked to a repository of electronic health record-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 managing incidental findings, and patient and clinician education methods.

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community., (© 2014 Chen et al.; Published by Cold Spring Harbor Laboratory Press.)

Background & Aims: Increasing grade of pancreatic intraepithelial neoplasia (PanIN) has been associated with progression to pancreatic ductal adenocarcinoma (PDAC). However, the mechanisms that control progression from PanINs to PDAC are not well understood. We investigated the genetic alterations involved in this process., Methods: Genomic DNA samples from laser-capture microdissected PDACs and adjacent PanIN2 and PanIN3 lesions from 10 patients with pancreatic cancer were analyzed by exome sequencing., Results: Similar numbers of somatic mutations were identified in PanINs and tumors, but the mutational load varied greatly among cases. Ten of the 15 isolated PanINs shared more than 50% of somatic mutations with associated tumors. Mutations common to tumors and clonally related PanIN2 and PanIN3 lesions were identified as genes that could promote carcinogenesis. KRAS and TP53 frequently were altered in PanINs and tumors, but few other recurrently modified genes were detected. Mutations in DNA damage response genes were prevalent in all samples. Genes that encode proteins involved in gap junctions, the actin cytoskeleton, the mitogen-activated protein kinase signaling pathway, axon guidance, and cell-cycle regulation were among the earliest targets of mutagenesis in PanINs that progressed to PDAC., Conclusions: Early stage PanIN2 lesions appear to contain many of the somatic gene alterations required for PDAC development., (Copyright © 2013 AGA Institute. Published by Elsevier Inc. All rights reserved.)

Background: The considerable genetic predisposition to deep vein thrombosis (DVT) is only partially accounted for by known genetic risk variants. Rare single-nucleotide variants (SNVs) of the coding areas of hemostatic genes may explain part of this missing heritability. The ADAMTS13 and VWF genes encode two interconnected proteins with fundamental hemostatic functions, the disruption of which may result in thrombosis., Objectives: To study the distribution and burden of rare coding SNVs of ADAMTS13 and VWF found by sequencing in cases and controls of DVT., Patients/methods: The protein-coding areas of 186 hemostatic/proinflammatory genes were sequenced by next-generation technology in 94 thrombophilia-negative patients with DVT and 98 controls. Gene-specific information on ADAMTS13 and VWF was used to study the association between DVT and rare coding SNVs of the two genes., Results: More than 70 billion base pairs of raw sequence data were produced to sequence the 700-kb target area with a median redundancy of × 45 in 192 individuals. Most of the 4366 SNVs identified were rare and non-synonymous, indicating pathogenetic potential. Rare (frequency of < 1%) and low-frequency (< 5%) coding SNVs of ADAMTS13 were associated with DVT (prevalence 17% vs. 4%; odds ratio [OR] 4.8 and 95% confidence interval [CI] 1.6-15.0 for rare coding; prevalence 36% vs. 23%, OR 1.9 and 95% CI 1.0-3.5 for low-frequency coding). Patients with rare coding SNVs of ADAMTS13 had lower plasma levels of ADAMTS-13 activity than patients without them. SNVs of VWF were not associated with DVT., Conclusions: We found an excess of rare coding SNVs of the ADAMTS13 gene in patients with DVT., (© 2013 International Society on Thrombosis and Haemostasis.)

Endometrial cancer is the most common gynecological malignancy, with more than 280,000 cases occurring annually worldwide. Although previous studies have identified important common somatic mutations in endometrial cancer, they have primarily focused on a small set of known cancer genes and have thus provided a limited view of the molecular basis underlying this disease. Here we have developed an integrated systems-biology approach to identifying novel cancer genes contributing to endometrial tumorigenesis. We first performed whole-exome sequencing on 13 endometrial cancers and matched normal samples, systematically identifying somatic alterations with high precision and sensitivity. We then combined bioinformatics prioritization with high-throughput screening (including both shRNA-mediated knockdown and expression of wild-type and mutant constructs) in a highly sensitive cell viability assay. Our results revealed 12 potential driver cancer genes including 10 tumor-suppressor candidates (ARID1A, INHBA, KMO, TTLL5, GRM8, IGFBP3, AKTIP, PHKA2, TRPS1, and WNT11) and two oncogene candidates (ERBB3 and RPS6KC1). The results in the "sensor" cell line were recapitulated by siRNA-mediated knockdown in endometrial cancer cell lines. Focusing on ARID1A, we integrated mutation profiles with functional proteomics in 222 endometrial cancer samples, demonstrating that ARID1A mutations frequently co-occur with mutations in the phosphatidylinositol 3-kinase (PI3K) pathway and are associated with PI3K pathway activation. siRNA knockdown in endometrial cancer cell lines increased AKT phosphorylation supporting ARID1A as a novel regulator of PI3K pathway activity. Our study presents the first unbiased view of somatic coding mutations in endometrial cancer and provides functional evidence for diverse driver genes and mutations in this disease.

Background: Up to 7% of patients with severe-to-profound deafness do not benefit from cochlear implantation. Given the high surgical implantation and clinical management cost of cochlear implantation (>$1 million lifetime cost), prospective identification of the worst performers would reduce unnecessary procedures and healthcare costs. Because cochlear implants bypass the membranous labyrinth but rely on the spiral ganglion for functionality, we hypothesize that cochlear implant (CI) performance is dictated in part by the anatomic location of the cochlear pathology that underlies the hearing loss. As a corollary, we hypothesize that because genetic testing can identify sites of cochlear pathology, it may be useful in predicting CI performance., Methods: 29 adult CI recipients with idiopathic adult-onset severe-to-profound hearing loss were studied. DNA samples were subjected to solution-based sequence capture and massively parallel sequencing using the OtoSCOPE(®) platform. The cohort was divided into three CI performance groups (good, intermediate, poor) and genetic causes of deafness were correlated with audiometric data to determine whether there was a gene-specific impact on CI performance., Results: The genetic cause of deafness was determined in 3/29 (10%) individuals. The two poor performers segregated mutations in TMPRSS3, a gene expressed in the spiral ganglion, while the good performer segregated mutations in LOXHD1, a gene expressed in the membranous labyrinth. Comprehensive literature review identified other good performers with mutations in membranous labyrinth-expressed genes; poor performance was associated with spiral ganglion-expressed genes., Conclusions: Our data support the underlying hypothesis that mutations in genes preferentially expressed in the spiral ganglion portend poor CI performance while mutations in genes expressed in the membranous labyrinth portend good CI performance. Although the low mutation rate in known deafness genes in this cohort likely relates to the ascertainment characteristics (postlingual hearing loss in adult CI recipients), these data suggest that genetic testing should be implemented as part of the CI evaluation to test this association prospectively., (Copyright © 2012 Elsevier B.V. All rights reserved.)

Background: Next-generation DNA sequencing is opening new avenues for genetic association studies in common diseases that, like deep vein thrombosis (DVT), have a strong genetic predisposition still largely unexplained by currently identified risk variants. In order to develop sequencing and analytical pipelines for the application of next-generation sequencing to complex diseases, we conducted a pilot study sequencing the coding area of 186 hemostatic/proinflammatory genes in 10 Italian cases of idiopathic DVT and 12 healthy controls., Results: A molecular-barcoding strategy was used to multiplex DNA target capture and sequencing, while retaining individual sequence information. Genomic libraries with barcode sequence-tags were pooled (in pools of 8 or 16 samples) and enriched for target DNA sequences. Sequencing was performed on ABI SOLiD-4 platforms. We produced > 12 gigabases of raw sequence data to sequence at high coverage (average: 42X) the 700-kilobase target area in 22 individuals. A total of 1876 high-quality genetic variants were identified (1778 single nucleotide substitutions and 98 insertions/deletions). Annotation on databases of genetic variation and human disease mutations revealed several novel, potentially deleterious mutations. We tested 576 common variants in a case-control association analysis, carrying the top-5 associations over to replication in up to 719 DVT cases and 719 controls. We also conducted an analysis of the burden of nonsynonymous variants in coagulation factor and anticoagulant genes. We found an excess of rare missense mutations in anticoagulant genes in DVT cases compared to controls and an association for a missense polymorphism of FGA (rs6050; p = 1.9 × 10(-5), OR 1.45; 95% CI, 1.22-1.72; after replication in > 1400 individuals)., Conclusions: We implemented a barcode-based strategy to efficiently multiplex sequencing of hundreds of candidate genes in several individuals. In the relatively small dataset of our pilot study we were able to identify bona fide associations with DVT. Our study illustrates the potential of next-generation sequencing for the discovery of genetic variation predisposing to complex diseases.

We demonstrate that phosphatidylinositol 3-kinase (PI3K) pathway aberrations occur in >80% of endometrioid endometrial cancers, with coordinate mutations of multiple PI3K pathway members being more common than predicted by chance. PIK3R1 (p85α) mutations occur at a higher rate in endometrial cancer than in any other tumor lineage, and PIK3R2 (p85β), not previously demonstrated to be a cancer gene, is also frequently mutated. The dominant activation event in the PI3K pathway appears to be PTEN protein loss. However, in tumors with retained PTEN protein, PI3K pathway mutations phenocopy PTEN loss, resulting in pathway activation. KRAS mutations are common in endometrioid tumors activating independent events from PI3K pathway aberrations. Multiple PIK3R1 and PIK3R2 mutations demonstrate gain of function, including disruption of a novel mechanism of pathway regulation wherein p85α dimers bind and stabilize PTEN. Taken together, the PI3K pathway represents a critical driver of endometrial cancer pathogenesis and a novel therapeutic target.

We present a genetic map for Xenopus tropicalis, consisting of 2886 Simple Sequence Length Polymorphism (SSLP) markers. Using a bioinformatics-based strategy, we identified unique SSLPs within the X. tropicalis genome. Scaffolds from X. tropicalis genome assembly 2.0 (JGI) were scanned for Simple Sequence Repeats (SSRs); unique SSRs were then tested for amplification and polymorphisms using DNA from inbred Nigerian and Ivory Coast individuals. Thus identified, the SSLPs were genotyped against a mapping cross panel of DNA samples from 190 F2 individuals. Nearly 4000 SSLPs were genotyped, yielding a 2886-marker genetic map consisting of 10 major linkage groups between 73 and 132cM in length, and 4 smaller linkage groups between 7 and 40cM. The total effective size of the map is 1658cM, and the average intermarker distance for each linkage group ranged from 0.27 to 0.75cM. Fluorescence In Situ Hybridization (FISH) was carried out using probes for genes located on mapped scaffolds to assign linkage groups to chromosomes. Comparisons of this map with the X. tropicalis genome Assembly 4.1 (JGI) indicate that the map provides representation of a minimum of 66% of the X. tropicalis genome, incorporating 758 of the approximately 1300 scaffolds over 100,000bp. The genetic map and SSLP marker database constitute an essential resource for genetic and genomic analyses in X. tropicalis., (Copyright © 2011 Elsevier Inc. All rights reserved.)