Exome sequencing of gastric adenocarcinoma identifies recurrent somatic mutations in cell adhesion and chromatin remodelling genes

Gastric cancer is a major cause of global cancer mortality. We surveyed the spectrum of somatic alterations in gastric cancer by sequencing the exomes of 15 gastric adenocarcinomas and their matched normal DNAs. Frequently mutated genes in the adenocarcinomas included TP53 (11/15 tumors), PIK3CA (3/15) and ARID1A (3/15). Cell adhesion was the most enriched biological pathway among the frequently mutated genes. A prevalence screening confirmed mutations in FAT4, a cadherin family gene, in 5% of gastric cancers (6/110) and FAT4 genomic deletions in 4% (3/83) of gastric tumors. Frequent mutations in chromatin remodeling genes (ARID1A, MLL3 and MLL) also occurred in 47% of the gastric cancers. We detected ARID1A mutations in 8% of tumors (9/110), which were associated with concurrent PIK3CA mutations and microsatellite instability. In functional assays, we observed both FAT4 and ARID1A to exert tumor-suppressor activity. Somatic inactivation of FAT4 and ARID1A may thus be key tumorigenic events in a subset of gastric cancers.

Selective extracellular vesicle-mediated export of an overlapping set of microRNAs from multiple cell types

Background: MicroRNAs (miRNAs) are a class of small RNA molecules that regulate expression of specific mRNA targets. They can be released from cells, often encapsulated within extracellular vesicles (EVs), and therefore have the potential to mediate intercellular communication. It has been suggested that certain miRNAs may be selectively exported, although the mechanism has yet to be identified. Manipulation of the miRNA content of EVs will be important for future therapeutic applications. We therefore wished to assess which endogenous miRNAs are enriched in EVs and how effectively an overexpressed miRNA would be exported.

Results: Small RNA libraries from HEK293T cells and vesicles before or after transfection with a vector for miR-146a overexpression were analysed by deep sequencing. A subset of miRNAs was found to be enriched in EVs; pathway analysis of their predicted target genes suggests a potential role in regulation of endocytosis. RT-qPCR in additional cell types and analysis of publicly available data revealed that many of these miRNAs tend to be widely preferentially exported. Whilst overexpressed miR-146a was highly enriched both in transfected cells and their EVs, the cellular:EV ratios of endogenous miRNAs were not grossly altered. MiR-451 was consistently the most highly exported miRNA in many different cell types. Intriguingly, Argonaute2 (Ago2) is required for miR-451 maturation and knock out of Ago2 has been shown to decrease expression of other preferentially exported miRNAs (eg miR-150 and miR-142-3p).

Conclusion: The global expression data provided by deep sequencing confirms that specific miRNAs are enriched in EVs released by HEK293T cells. Observation of similar patterns in a range of cell types suggests that a common mechanism for selective miRNA export may exist.

Molecular cloning, expression, and regulation in Escherichia coli K-12 of a chromosome-mediated aerobactin iron transport system from a human invasive isolate of E. coli K1

We have cloned chromosomal genes determining the aerobactin iron transport system from the Escherichia coli K1 strain VW187. Mapping and hybridization experiments showed that the VW187 aerobactin region was identical to that of the plasmid ColV-K30. However, in the E. coli K-12 background, the biosynthesis of both siderophore and ferric aerobactin receptor encoded by the VW187-derived recombinant plasmids was not repressed by iron to the same extent found when a recombinant plasmid derived from pColV-K30 was used. RNA-DNA dot-blot hybridization experiments demonstrated that the aerobactin-specific mRNA synthesized by the VW187-derived clones was not iron regulated in E. coli K-12. In contrast, the synthesis of aerobactin and its receptor in strain VW187 was completely repressed by iron regardless of whether the recombinant plasmids originated from VW187 or pColV-K30. Similar results were obtained with gene fusions in which a promoterless lac operon was placed under the control of aerobactin promoter regions of either chromosome- or plasmid-mediated aerobactin systems. DNA sequencing of the chromosomal aerobactin promoter region showed changes in bases located immediately upstream to the -35 region compared with the corresponding region in pColV-K30, which is known to be part of the binding site for the Fur repressor protein.

Fast RNA conjugations on solid phase by strain-promoted cycloadditions

Strain promoted cycloaddition is presented as a tool for RNA conjugation on the solid phase; RNA-cyclooctyne conjugates are prepared by cycloaddition to both azide (strain-promoted azide-alkyne cycloaddition, SPAAC) and nitrile oxide dipoles (strain-promoted nitrile oxide-alkyne cycloaddition, SPNOAC). The conjugation is compatible with 2'-OMe blocks and with 2'-O-TBDMS protection on the ribose moieties of the sugar. Nitrile oxide dipoles are found to be more reactive click partners than azides. The conjugation proceeds within 10 min in aqueous solvents, at room temperature without any metal catalyst and tolerates dipoles of varying steric bulk and electronic demands, including pyrenyl, coumarin and dabcyl derivatives. © 2012 The Royal Society of Chemistry.

Detecting sequence dependent transcriptional pauses from RNA and protein number time series

Background: Evidence suggests that in prokaryotes sequence-dependent transcriptional pauses a?ect the dynamics of transcription and translation, as well as of small genetic circuits. So far, a few pause-prone sequences have been identi?ed from in vitro measurements of transcription elongation kinetics.

Results: Using a stochastic model of gene expression at the nucleotide and codon levels with realistic parameter values, we investigate three di?erent but related questions and present statistical methods for their analysis. First, we show that information from in vivo RNA and protein temporal numbers is su?cient to discriminate between models with and without a pause site in their coding sequence. Second, we demonstrate that it is possible to separate a large variety of models from each other with pauses of various durations and locations in the template by means of a hierarchical clustering and a random forest classi?er. Third, we introduce an approximate likelihood function that allows to estimate the location of a pause site.

Conclusions: This method can aid in detecting unknown pause-prone sequences from temporal measurements of RNA and protein numbers at a genome-wide scale and thus elucidate possible roles that these sequences play in the dynamics of genetic networks and phenotype.

A whole-genome massively parallel sequencing analysis of BRCA1 mutant oestrogen receptor-negative and -positive breast cancers

BRCA1 encodes a tumour suppressor protein that plays pivotal roles in homologous recombination (HR) DNA repair, cell-cycle checkpoints, and transcriptional regulation. BRCA1 germline mutations confer a high risk of early-onset breast and ovarian cancer. In more than 80% of cases, tumours arising in BRCA1 germline mutation carriers are oestrogen receptor (ER)-negative; however, up to 15% are ER-positive. It has been suggested that BRCA1 ER-positive breast cancers constitute sporadic cancers arising in the context of a BRCA1 germline mutation rather than being causally related to BRCA1 loss-of-function. Whole-genome massively parallel sequencing of ER-positive and ER-negative BRCA1 breast cancers, and their respective germline DNAs, was used to characterize the genetic landscape of BRCA1 cancers at base-pair resolution. Only BRCA1 germline mutations, somatic loss of the wild-type allele, and TP53 somatic mutations were recurrently found in the index cases. BRCA1 breast cancers displayed a mutational signature consistent with that caused by lack of HR DNA repair in both ER-positive and ER-negative cases. Sequencing analysis of independent cohorts of hereditary BRCA1 and sporadic non-BRCA1 breast cancers for the presence of recurrent pathogenic mutations and/or homozygous deletions found in the index cases revealed that DAPK3, TMEM135, KIAA1797, PDE4D, and GATA4 are potential additional drivers of breast cancers. This study demonstrates that BRCA1 pathogenic germline mutations coupled with somatic loss of the wild-type allele are not sufficient for hereditary breast cancers to display an ER-negative phenotype, and has led to the identification of three potential novel breast cancer genes (ie DAPK3, TMEM135, and GATA4).

RNA interference targeting cathepsin B of the carcinogenic liver fluke, Opisthorchis viverrini

Functional genomics have not been reported for Opisthorchis viverrini or the related fish-borne fluke, Clonorchis sinensis. Here we describe the introduction by square wave electroporation of Cy3-labeled small RNA into adult O. viverrini worms. Adult flukes were subjected to square wave electroporation employing a single pulse for 20 ms of 125V in the presence of 50 µg/ml of Cy3-siRNA. The parasites tolerated this manipulation and, at 24 and 48 h after electroporation, fluorescence from the Cy3-siRNA was evident throughout the parenchyma of the worms, with strong fluorescence evident in the guts and reproductive organs of the adult worms. Second, other worms were treated using the same electroporation settings with double stranded RNA targeting an endogenous papain-like cysteine protease, cathepsin B. This manipulation resulted in a significant reduction in specific mRNA levels encoding cathepsin B, and a significant reduction in cathepsin B activity against the diagnostic peptide, Z-Arg-Arg-AMC. This appears to be the first report of introduction of reporter genes into O. viverrini and the first report of experimental RNA interference (RNAi) in this fluke. The findings indicated the presence of an intact RNAi pathway in these parasites which, in turn, provides an opportunity to probe gene functions in this neglected tropical disease pathogen.

Long-range Transcriptome Sequencing Reveals Cancer Cell Growth Regulatory Chimeric mRNA

mRNA chimeras from chromosomal translocations often play a role as transforming oncogenes. However, cancer transcriptomes also contain mRNA chimeras that may play a role in tumor development, which arise as transcriptional or post-transcriptional events. To identify such chimeras, we developed a deterministic screening strategy for long-range sequence analysis. High-throughput, long-read sequencing was then performed on cDNA libraries from major tumor histotypes and corresponding normal tissues. These analyses led to the identification of 378 chimeras, with an unexpectedly high frequency of expression (˜2 x 10(-5) of all mRNA). Functional assays in breast and ovarian cancer cell lines showed that a large fraction of mRNA chimeras regulates cell replication. Strikingly, chimeras were shown to include both positive and negative regulators of cell growth, which functioned as such in a cell-type-specific manner. Replication-controlling chimeras were found to be expressed by most cancers from breast, ovary, colon, uterus, kidney, lung, and stomach, suggesting a widespread role in tumor development.

Old drug, New target:Ellipticines selectively inhibit RNA Polymerase I transcription

Transcription byRNApolymerase I (Pol-I) is the main driving force behind ribosome biogenesis, a fundamental cellular process that requires the coordinated transcription of all three nuclear polymerases. Increased Pol-I transcription and the concurrent increase in ribosome biogenesis has been linked to the high rates of proliferation in cancers. The ellipticine family contains a number of potent anticancer therapeutic agents, some having progressed to stage I and II clinical trials; however, the mechanism by which many of the compounds work remains unclear. It has long been thought that inhibition of Top2 is the main reason behind the drugs antiproliferative effects. Here we report that a number of the ellipticines, including 9-hydroxyellipticine, are potent and specific inhibitors of Pol-I transcription, with IC50 in vitro and in cells in the nanomolar range. Essentially, the drugs did not affect Pol-II and Pol-III transcription, demonstrating a high selectivity.Wehave shown that Pol-I inhibition occurs by a p53-, ATM/ATR-, and Top2-independent mechanism. We discovered that the drug influences the assembly and stability of preinitiation complexes by targeting the interaction between promoter recognition factor SL1 and the rRNA promoter. Our findings will have an impact on the design and development of novel therapeutic agents specifically targeting ribosome biogenesis.

Sequencing and its consequences:Path dependence and the relationships between genetics and medicalization

Both advocacy for and critiques of the Human Genome Project assume a self-sustaining relationship between genetics and. medicalization. However, this assumption ignores the ways in which the meanings of genetic research are conditional on its position in sequences of events. Based, on analyses of three conditions for which at least one putative gene or genetic marker has been identified, this article argues that critical junctures in the institutional stabilization of phenotypes and the mechanisms that sustain such classifications over time configure the practices and meanings of genetic research. Path dependence is critical to understanding the lack of consistent fit between genetics and medlcalization.

Exploring the Transcriptome of Atlantic Salmon (Salmo salar) Skin, a Major Defense Organ

The skin of fish is the first line of defense against pathogens and parasites. The skin transcriptome of the Atlantic salmon is poorly characterized, and currently only 2,089 expressed sequence tags (ESTs) out of a total of half a million sequences are generated from skin-derived cDNA libraries. The primary aim of this study was to enhance the transcriptomic knowledge of salmon skin by using next-generation sequencing (NGS) technology, namely the Roche-454 platform. An equimolar mixture of high-quality RNA from skin and epidermal samples of salmon reared in either freshwater or seawater was used for 454-sequencing. This technique yielded over 600,000 reads, which were assembled into 34,696 isotigs using Newbler. Of these isotigs, 12 % had not been sequenced in Atlantic salmon, hence representing previously unreported salmon mRNAs that can potentially be skin-specific. Many full-length genes have been acquired, representing numerous biological processes. Mucin proteins are the main structural component of mucus and we examined in greater detail the sequences we obtained for these genes. Several isotigs exhibited homology to mammalian mucins (MUC2, MUC5AC and MUC5B). Mucin mRNAs are generally > 10 kbp and contain large repetitive units, which pose a challenge towards full-length sequence discovery. To date, we have not unearthed any full-length salmon mucin genes with this dataset, but have both N- and C-terminal regions of a mucin type 5. This highlights the fact that, while NGS is indeed a formidable tool for sequence data mining of non-model species, it must be complemented with additional experimental and bioinformatic work to characterize some mRNA sequences with complex features.

Resolving candidate genes of mouse skeletal muscle QTL via RNA-Seq and expression network analyses

BACKGROUND:

We have recently identified a number of Quantitative Trait Loci (QTL) contributing to the 2-fold muscle weight difference between the LG/J and SM/J mouse strains and refined their confidence intervals. To facilitate nomination of the candidate genes responsible for these differences we examined the transcriptome of the tibialis anterior (TA) muscle of each strain by RNA-Seq.

13,726 genes were expressed in mouse skeletal muscle. Intersection of a set of 1061 differentially expressed transcripts with a mouse muscle Bayesian Network identified a coherent set of differentially expressed genes that we term the LG/J and SM/J Regulatory Network (LSRN). The integration of the QTL, transcriptome and the network analyses identified eight key drivers of the LSRN (Kdr, Plbd1, Mgp, Fah, Prss23, 2310014F06Rik, Grtp1, Stk10) residing within five QTL regions, which were either polymorphic or differentially expressed between the two strains and are strong candidates for quantitative trait genes (QTGs) underlying muscle mass. The insight gained from network analysis including the ability to make testable predictions is illustrated by annotating the LSRN with knowledge-based signatures and showing that the SM/J state of the network corresponds to a more oxidative state. We validated this prediction by NADH tetrazolium reductase staining in the TA muscle revealing higher oxidative potential of the SM/J compared to the LG/J strain (p<0.03).

Thus, integration of fine resolution QTL mapping, RNA-Seq transcriptome information and mouse muscle Bayesian Network analysis provides a novel and unbiased strategy for nomination of muscle QTGs.

The Interlaboratory Robustness Of Next-Generation Sequencing (IRON) Study Phase II: Deep-Sequencing Analyses Of Hematological Malignancies Performed In 8,867 Cases By An International Network Involving 27 Laboratories

Introduction: Amplicon deep-sequencing using second-generation sequencing technology is an innovative molecular diagnostic technique and enables a highly-sensitive detection of mutations. As an international consortium we had investigated previously the robustness, precision, and reproducibility of 454 amplicon next-generation sequencing (NGS) across 10 laboratories from 8 countries (Leukemia, 2011;25:1840-8).

Aims: In Phase II of the study, we established distinct working groups for various hematological malignancies, i.e. acute myeloid leukemia (AML), acute lymphoblastic leukemia (ALL), chronic lymphocytic leukemia (CLL), chronic myelogenous leukemia (CML), myelodysplastic syndromes (MDS), myeloproliferative neoplasms (MPN), and multiple myeloma. Currently, 27 laboratories from 13 countries are part of this research consortium. In total, 74 gene targets were selected by the working groups and amplicons were developed for a NGS deep-sequencing assay (454 Life Sciences, Branford, CT). A data analysis pipeline was developed to standardize mutation interpretation both for accessing raw data (Roche Amplicon Variant Analyzer, 454 Life Sciences) and variant interpretation (Sequence Pilot, JSI Medical Systems, Kippenheim, Germany).

Results: We will report on the design, standardization, quality control aspects, landscape of mutations, as well as the prognostic and predictive utility of this assay in a cohort of 8,867 cases. Overall, 1,146 primer sequences were designed and tested. In detail, for example in AML, 924 cases had been screened for CEBPA mutations. RUNX1 mutations were analyzed in 1,888 cases applying the deep-sequencing read counts to study the stability of such mutations at relapse and their utility as a biomarker to detect residual disease. Analyses of DNMT3A (n=1,041) were focused to perform landscape investigations and to address the prognostic relevance. Additionally, this working group is focusing on TET2, ASXL1, and TP53 analyses. A novel prognostic model is being developed allowing stratification of AML into prognostic subgroups based on molecular markers only. In ALL, 1,124 pediatric and adult cases have been screened, including 763 assays for TP53 mutations both at diagnosis and relapse of ALL. Pediatric and adult leukemia expert labs developed additional content to study the mutation incidence of other B and T lineage markers such as IKZF1, JAK2, IL7R, PAX5, EP300, LEF1, CRLF2, PHF6, WT1, JAK1, PTEN, AKT1, IL7R, NOTCH1, CREBBP, or FBXW7. Further, the molecular landscape of CLL is changing rapidly. As such, a separate working group focused on analyses including NOTCH1, SF3B1, MYD88, XPO1, FBXW7 and BIRC3. Currently, 922 cases were screened to investigate the range of mutational burden of NOTCH1 mutations for their prognostic relevance. In MDS, RUNX1 mutation analyses were performed in 977 cases. The prognostic relevance of TP53 mutations in MDS was assessed in additional 327 cases, including isolated deletions of chromosome 5q. Next, content was developed targeting genes of the cellular splicing component, e.g. SF3B1, SRSF2, U2AF1, and ZRSR2. In BCR-ABL1-negative MPN, nine genes of interest (JAK2, MPL, TET2, CBL, KRAS, EZH2, IDH1, IDH2, ASXL1) have been analyzed in a cohort of 155 primary myelofibrosis cases searching for novel somatic mutations and addressing their relevance for disease progression and leukemia transformation. Moreover, an assay was developed and applied to CMML cases allowing the simultaneous analysis of 25 leukemia-associated target genes in a single sequencing run using just 20 ng of starting DNA. Finally, nine laboratories are studying CML, applying ultra-deep sequencing of the BCR-ABL1 tyrosine kinase domain. Analyses were performed on 615 cases investigating the dynamics of expansion of mutated clones under various tyrosine kinase inhibitor therapies.

Conclusion: Molecular characterization of hematological malignancies today requires high diagnostic sensitivity and specificity. As part of the IRON-II study, a network of laboratories analyzed a variety of disease entities applying amplicon-based NGS assays. Importantly, the consortium not only standardized assay design for disease-specific panels, but also achieved consensus on a common data analysis pipeline for mutation interpretation. Distinct working groups have been forged to address scientific tasks and in total 8,867 cases had been analyzed thus far.