Pharmacogenomics of 17-alpha hydroxyprogesterone caproate for recurrent preterm birth: a case-control study.

OBJECTIVE: To compare maternal genotypes between women with and without significant prolongation of pregnancy in the setting of 17-alpha hydroxyprogesterone caproate (17-P) administration for the prevention of recurrent preterm birth (PTB). DESIGN: Case-control. SETTING: Three tertiary-care centres across the USA. POPULATION: Women (n = 99) with ≥ 1 prior singleton spontaneous PTB, receiving 17-P. METHODS: Women were classified as having successful prolongation of pregnancy during the 17-P treated pregnancy, in two ways: (1) Definition A: success/non-success based on difference in gestational age at delivery between 17-P-treated and untreated pregnancies (success: delivered ≥ 3 weeks later with 17-P) and (2) Definition B: success/non-success based on reaching term (success: delivered at term with 17-P). MAIN OUTCOME MEASURES: To assess genetic variation, all women underwent whole exome sequencing. Between-group sequence variation was analysed with the Variant Annotation, Analysis, and Search Tool (VAAST). Genes scored by VAAST with P < 0.05 were then analysed with two online tools: (1) Protein ANalysis THrough Evolutionary Relationships (PANTHER) and (2) Database for Annotation, Visualization, and Integrated Discovery (DAVID). RESULTS: Using Definition A, there were 70 women with successful prolongation and 29 without; 1375 genes scored by VAAST had P < 0.05. Using Definition B, 47 women had successful prolongation and 52 did not; 1039 genes scored by VAAST had P < 0.05. PANTHER revealed key differences in gene ontology pathways. Many genes from definition A were classified as prematurity genes (P = 0.026), and those from definition B as pharmacogenetic genes (P = 0.0018); (P, non-significant after Bonferroni correction). CONCLUSION: A novel analytic approach revealed several genetic differences among women delivering early vs later with 17-P. TWEETABLE ABSTRACT: Several key genetic differences are present in women with recurrent preterm birth despite 17-P treatment.

A general approach to single-nucleotide polymorphism discovery.

Single-nucleotide polymorphisms (SNPs) are the most abundant form of human genetic variation and a resource for mapping complex genetic traits. The large volume of data produced by high-throughput sequencing projects is a rich and largely untapped source of SNPs (refs 2, 3, 4, 5). We present here a unified approach to the discovery of variations in genetic sequence data of arbitrary DNA sources. We propose to use the rapidly emerging genomic sequence as a template on which to layer often unmapped, fragmentary sequence data and to use base quality values to discern true allelic variations from sequencing errors. By taking advantage of the genomic sequence we are able to use simpler yet more accurate methods for sequence organization: fragment clustering, paralogue identification and multiple alignment. We analyse these sequences with a novel, Bayesian inference engine, POLYBAYES, to calculate the probability that a given site is polymorphic. Rigorous treatment of base quality permits completely automated evaluation of the full length of all sequences, without limitations on alignment depth. We demonstrate this approach by accurate SNP predictions in human ESTs aligned to finished and working-draft quality genomic sequences, a data set representative of the typical challenges of sequence-based SNP discovery.

The sequence of the human genome.

A 2.91-billion base pair (bp) consensus sequence of the euchromatic portion of the human genome was generated by the whole-genome shotgun sequencing method. The 14.8-billion bp DNA sequence was generated over 9 months from 27,271,853 high-quality sequence reads (5.11-fold coverage of the genome) from both ends of plasmid clones made from the DNA of five individuals. Two assembly strategies-a whole-genome assembly and a regional chromosome assembly-were used, each combining sequence data from Celera and the publicly funded genome effort. The public data were shredded into 550-bp segments to create a 2.9-fold coverage of those genome regions that had been sequenced, without including biases inherent in the cloning and assembly procedure used by the publicly funded group. This brought the effective coverage in the assemblies to eightfold, reducing the number and size of gaps in the final assembly over what would be obtained with 5.11-fold coverage. The two assembly strategies yielded very similar results that largely agree with independent mapping data. The assemblies effectively cover the euchromatic regions of the human chromosomes. More than 90% of the genome is in scaffold assemblies of 100,000 bp or more, and 25% of the genome is in scaffolds of 10 million bp or larger. Analysis of the genome sequence revealed 26,588 protein-encoding transcripts for which there was strong corroborating evidence and an additional approximately 12,000 computationally derived genes with mouse matches or other weak supporting evidence. Although gene-dense clusters are obvious, almost half the genes are dispersed in low G+C sequence separated by large tracts of apparently noncoding sequence. Only 1.1% of the genome is spanned by exons, whereas 24% is in introns, with 75% of the genome being intergenic DNA. Duplications of segmental blocks, ranging in size up to chromosomal lengths, are abundant throughout the genome and reveal a complex evolutionary history. Comparative genomic analysis indicates vertebrate expansions of genes associated with neuronal function, with tissue-specific developmental regulation, and with the hemostasis and immune systems. DNA sequence comparisons between the consensus sequence and publicly funded genome data provided locations of 2.1 million single-nucleotide polymorphisms (SNPs). A random pair of human haploid genomes differed at a rate of 1 bp per 1250 on average, but there was marked heterogeneity in the level of polymorphism across the genome. Less than 1% of all SNPs resulted in variation in proteins, but the task of determining which SNPs have functional consequences remains an open challenge.

The genome sequence of Drosophila melanogaster.

The fly Drosophila melanogaster is one of the most intensively studied organisms in biology and serves as a model system for the investigation of many developmental and cellular processes common to higher eukaryotes, including humans. We have determined the nucleotide sequence of nearly all of the approximately 120-megabase euchromatic portion of the Drosophila genome using a whole-genome shotgun sequencing strategy supported by extensive clone-based sequence and a high-quality bacterial artificial chromosome physical map. Efforts are under way to close the remaining gaps; however, the sequence is of sufficient accuracy and contiguity to be declared substantially complete and to support an initial analysis of genome structure and preliminary gene annotation and interpretation. The genome encodes approximately 13,600 genes, somewhat fewer than the smaller Caenorhabditis elegans genome, but with comparable functional diversity.

Comparative genomics of the eukaryotes.

A comparative analysis of the genomes of Drosophila melanogaster, Caenorhabditis elegans, and Saccharomyces cerevisiae-and the proteins they are predicted to encode-was undertaken in the context of cellular, developmental, and evolutionary processes. The nonredundant protein sets of flies and worms are similar in size and are only twice that of yeast, but different gene families are expanded in each genome, and the multidomain proteins and signaling pathways of the fly and worm are far more complex than those of yeast. The fly has orthologs to 177 of the 289 human disease genes examined and provides the foundation for rapid analysis of some of the basic processes involved in human disease.

RNA sequencing confirms similarities between PPI-responsive oesophageal eosinophilia and eosinophilic oesophagitis.

BACKGROUND: Although current American guidelines distinguish proton pump inhibitor-responsive oesophageal eosinophilia (PPI-REE) from eosinophilic oesophagitis (EoE), these entities are broadly similar. While two microarray studies showed that they have similar transcriptomes, more extensive RNA sequencing studies have not been done previously. AIM: To determine whether RNA sequencing identifies genetic markers distinguishing PPI-REE from EoE. METHODS: We retrospectively examined 13 PPI-REE and 14 EoE biopsies, matched for tissue eosinophil content, and 14 normal controls. Patients and controls were not PPI-treated at the time of biopsy. We did RNA sequencing on formalin-fixed, paraffin-embedded tissue, with differential expression confirmation by quantitative polymerase chain reaction (PCR). We validated the use of formalin-fixed, paraffin-embedded vs RNAlater-preserved tissue, and compared our formalin-fixed, paraffin-embedded EoE results to a prior EoE study. RESULTS: By RNA sequencing, no genes were differentially expressed between the EoE and PPI-REE groups at the false discovery rate (FDR) ≤0.01 level. Compared to normal controls, 1996 genes were differentially expressed in the PPI-REE group and 1306 genes in the EoE group. By less stringent criteria, only MAPK8IP2 was differentially expressed between PPI-REE and EoE (FDR = 0.029, 2.2-fold less in EoE than in PPI-REE), with similar results by PCR. KCNJ2, which was differentially expressed in a prior study, was similar in the EoE and PPI-REE groups by both RNA sequencing and real-time PCR. CONCLUSION: Eosinophilic oesophagitis and PPI-REE have comparable transcriptomes, confirming that they are part of the same disease continuum.

Age-related mutations and chronic myelomonocytic leukemia.

Chronic myelomonocytic leukemia (CMML) is a hematologic malignancy nearly confined to the elderly. Previous studies to determine incidence and prognostic significance of somatic mutations in CMML have relied on candidate gene sequencing, although an unbiased mutational search has not been conducted. As many of the genes commonly mutated in CMML were recently associated with age-related clonal hematopoiesis (ARCH) and aged hematopoiesis is characterized by a myelomonocytic differentiation bias, we hypothesized that CMML and aged hematopoiesis may be closely related. We initially established the somatic mutation landscape of CMML by whole exome sequencing followed by gene-targeted validation. Genes mutated in ⩾10% of patients were SRSF2, TET2, ASXL1, RUNX1, SETBP1, KRAS, EZH2, CBL and NRAS, as well as the novel CMML genes FAT4, ARIH1, DNAH2 and CSMD1. Most CMML patients (71%) had mutations in ⩾2 ARCH genes and 52% had ⩾7 mutations overall. Higher mutation burden was associated with shorter survival. Age-adjusted population incidence and reported ARCH mutation rates are consistent with a model in which clinical CMML ensues when a sufficient number of stochastically acquired age-related mutations has accumulated, suggesting that CMML represents the leukemic conversion of the myelomonocytic-lineage-biased aged hematopoietic system.

Identifying rare variants for genetic risk through a combined pedigree and phenotype approach: application to suicide and asthma.

Suicidal behavior is a complex disorder, with evidence for genetic risk independent of other genetic risk factors including psychiatric disorders. Since 1996, over 3000 DNA samples from Utah suicide decedents have been collected and banked for research use through the Utah Medical Examiner. In addition, over 12,000 Utah suicides were identified through examination of death certificates back to 1904. By linking this data with the Utah Population Database, we have identified multiple extended pedigrees with increased risk for suicide completion. A number of medical conditions co-occur with suicide, including asthma, and this study was undertaken to identify genetic risk common to asthma and suicide. This study tests the hypothesis that a particular comorbid condition may identify a more homogeneous genetic subgroup, facilitating the identification of specific genetic risk factors in that group. From pedigrees at increased risk for suicide, we identified three pedigrees also at significantly increased familial risk for asthma. Five suicide decedents from each of these pedigrees, plus an additional three decedents not from these pedigrees with diagnosed asthma, and 10 decedents with close relatives with asthma were genotyped. Results were compared with 183 publicly available unaffected control exomes from 1000 Genomes and CEPH (Centre d'etude du polymorphisme humain) samples genotyped on the same platform. A further 432 suicide decedents were also genotyped as non-asthma suicide controls. Genotyping was done using the Infinium HumanExome BeadChip. For analysis, we used the pedigree extension of Variant Annotation, Analysis and Search Tool (pVAAST) to calculate the disease burden of each gene. The Phenotype Driven Variant Ontological Re-ranking tool (Phevor) then re-ranked our pVAAST results in context of the phenotype. Using asthma as a seed phenotype, Phevor traversed biomedical ontologies and identified genes with similar biological properties to those known to result in asthma. Our top associated genes included those related to neurodevelopment or neural signaling (brain-derived neurotrophic factor (BDNF), neutral sphingomyelinase 2 (SMPD2), homeobox b2 (HOXB2), neural cell adhesion molecule (NCAM2), heterogeneous nuclear ribonucleoprotein A0 (HNRNPA0)), inflammation (free fatty acid receptor 2 (FFAR2)) and inflammation with additional evidence of neuronal involvement (oxidized low density lipoprotein receptor 1 (OLR1), toll-like receptor 3 (TLR3)). Of particular interest, BDNF has been previously implicated in both psychiatric disorders and asthma. Our results demonstrate the utility of combining pedigree and co-occurring phenotypes to identify rare variants associated with suicide risk in conjunction with specific co-occurring conditions.

Identification of key concepts in biomedical literature using a modified Markov heuristic.

MOTIVATION: The recent explosion of interest in mining the biomedical literature for associations between defined entities such as genes, diseases and drugs has made apparent the need for robust methods of identifying occurrences of these entities in biomedical text. Such concept-based indexing is strongly dependent on the availability of a comprehensive ontology or lexicon of biomedical terms. However, such ontologies are very difficult and expensive to construct, and often require extensive manual curation to render them suitable for use by automatic indexing programs. Furthermore, the use of statistically salient noun phrases as surrogates for curated terminology is not without difficulties, due to the lack of high-quality part-of-speech taggers specific to medical nomenclature. RESULTS: We describe a method of improving the quality of automatically extracted noun phrases by employing prior knowledge during the HMM training procedure for the tagger. This enhancement, when combined with appropriate training data, can greatly improve the quality and relevance of the extracted phrases, thereby enabling greater accuracy in downstream literature mining tasks.

Trimethylpsoralen induces small deletion mutations in Caenorhabditis elegans.

To examine the mutagenic spectrum of 4,5',8-trimethylpsoralen (TMP) in Caenorhabditis elegans, we isolated mutations in the unc-22 and pal-1 genes following TMP mutagenesis and analyzed them for restriction fragment length polymorphisms by Southern blot. Eleven of 21 unc-22 mutations exhibited restriction fragment length polymorphisms, 8 of which were deletions of between 0.10 and 15 kb in length. Both of two pal-1 mutations were also small deletions within this size range. Comparison of our results with previous studies on mutagenesis by gamma-rays and x-rays suggests that the mutagenic spectrum of TMP may be similar. TMP should be useful in generating mutations that cause complete loss of function of single genes and that are likely to result in allele-specific DNA polymorphisms.

A BMP homolog acts as a dose-dependent regulator of body size and male tail patterning in Caenorhabditis elegans.

We cloned the dbl-1 gene, a C. elegans homolog of Drosophila decapentaplegic and vertebrate BMP genes. Loss-of-function mutations in dbl-1 cause markedly reduced body size and defective male copulatory structures. Conversely, dbl-1 overexpression causes markedly increased body size and partly complementary male tail phenotypes, indicating that DBL-1 acts as a dose-dependent regulator of these processes. Evidence from genetic interactions indicates that these effects are mediated by a Smad signaling pathway, for which DBL-1 is a previously unidentified ligand. Our study of the dbl-1 expression pattern suggests a role for neuronal cells in global size regulation as well as male tail patterning.

Identification and characterization of a p53 homologue in Drosophila melanogaster.

The tumor suppressor gene p53 in mammalian cells plays a critical role in safeguarding the integrity of genome. It functions as a sequence-specific transcription factor. Upon activation by a variety of cellular stresses, p53 transactivates downstream target genes, through which it regulates cell cycle and apoptosis. However, little is known about p53 in invertebrates. Here we report the identification and characterization of a Drosophila p53 homologue gene, dp53. dp53 encodes a 385-amino acid protein with significant homology to human p53 (hp53) in the region of the DNA-binding domain, and to a lesser extent the tetramerization domain. Purified dp53 DNA-binding domain protein was shown to bind to the consensus hp53-binding site by gel mobility analysis. In transient transfection assays, expression of dp53 in Schneider cells transcriptionally activated promoters that contained consensus hp53-responsive elements. Moreover, a mutant dp53 (Arg-155 to His-155), like its hp53 counterpart mutant, exerted a dominant-negative effect on transactivation. Ectopic expression of dp53 in Drosophila eye disk caused cell death and led to a rough eye phenotype. dp53 is expressed throughout the development of Drosophila with highest expression levels in early embryogenesis, which has a maternal component. Consistent with this, dp53 RNA levels were high in the nurse cells of the ovary. It appears that p53 is structurally and functionally conserved from flies to mammals. Drosophila will provide a useful genetic system to the further study of the p53 network.