Development and Validation of an RNA Sequencing Assay for Gene Fusion Detection in Formalin-Fixed, Paraffin-Embedded Tumors.

RNA sequencing (RNA-seq) is a well-validated tool for detecting gene fusions in fresh-frozen tumors; however, RNA-seq is much more challenging to use with formalin-fixed, paraffin-embedded (FFPE) tumor samples. We evaluated the performance of RNA-seq to detect gene fusions in clinical FFPE tumor samples. Our assay identified all 15 spiked-in NTRK fusions from RNA reference material and six known fusions from five cancer cell lines. Limit of detection for the assay was assessed with a series of dilutions of RNA from the cell line H2228. These fusions can be detected when the dilution is down to 10%. Good intra-assay and interassay reproducibility was observed in three specimens. For clinical validation, the assay detected 10 of 12 fusions initially identified by a DNA panel (covering 23 gene fusions) in clinical specimens (83.3% sensitivity), whereas one fusion (MET fusion) was identified in another 34 fusion-negative tumor specimens as determined by the DNA panel (negative prediction value of 94.3%). This MET fusion was confirmed by RT-PCR and Sanger sequencing, which found that this is a false-negative result for the DNA panel. The assay also identified 26 extra fusions not covered by the DNA panel, 20 (76.9%) of which were validated by RT-PCR and Sanger sequencing. Therefore, this RNA assay has reasonable performance and could complement DNA-based next-generation sequencing assays.

Bacterial Diversity of Intestinal Microbiota in Patients with Substance Use Disorders Revealed by 16S rRNA Gene Deep Sequencing.

Substance abuse and addiction are worldwide concerns. In China, populated with over 1.3 billion people, emerging studies show a steady increase in substance abuse and substance-related problems. Some of the major challenges include a lack of an effective evaluation platform to determine the health status of substance-addicted subjects. It is known that the intestinal microbiota is associated to the occurrence and development of human diseases. However, the changes of bacterial diversity of intestinal microbiota in substance-addicted subjects have not been clearly characterized. Herein, we examined the composition and diversity of intestinal microbiota in 45 patients with substance use disorders (SUDs) and in 48 healthy controls (HCs). The results show that the observed species diversity index and the abundance of Thauera, Paracoccus, and Prevotella are significantly higher in SUDs compared to HCs. The functional diversity of the putative metagenomes analysis reveals that pathways including translation, DNA replication and repair, and cell growth and death are over-represented while cellular processes and signaling, and metabolism are under-represented in SUDs. Overall, the analyses show that there seem to be changes in the microbiota that are associated with substance use across an array of SUDs, providing fundamental knowledge for future research in substance-addiction assessment tests.

Genome-wide open chromatin regions and their effects on the regulation of silk protein genes in Bombyx mori.

Nucleosome-depleted open chromatin regions (OCRs) often harbor transcription factor (TF) binding sites that are associated with active DNA regulatory elements. To investigate the regulation of silk-protein genes, DNA molecules isolated from the silk glands of third-day fifth-instar silkworm larvae and embryo-derived (BmE) cells were subjected to formal dehyde-assisted isolation of regulatory elements (FAIRE) and high-throughput sequencing. In total, 68,000 OCRs were identified, and a number of TF-binding motifs were predicted. In particular, OCRs located near silk-protein genes contained potential binding sites for functional TFs. Moreover, many TFs were found to bind to clusters of OCRs upstream of silk-protein genes, and to regulate the expression of these genes. The expression of silk protein genes may be related not only to regulating TFs (such as fkh, Bmdimm, and Bmsage), but also to developmental and hormone-induced TFs (such as zen, eve, Br, and eip74ef). Elucidation of genome-wide OCRs and their regulatory motifs in silk protein genes will provide valuable data and clues for characterizing the mechanisms of transcriptional control of silk protein genes.

Expression map of a complete set of gustatory receptor genes in chemosensory organs of Bombyx mori.

Most lepidopteran species are herbivores, and interaction with host plants affects their gene expression and behavior as well as their genome evolution. Gustatory receptors (Grs) are expected to mediate host plant selection, feeding, oviposition and courtship behavior. However, due to their high diversity, sequence divergence and extremely low level of expression it has been difficult to identify precisely a complete set of Grs in Lepidoptera. By manual annotation and BAC sequencing, we improved annotation of 43 gene sequences compared with previously reported Grs in the most studied lepidopteran model, the silkworm, Bombyx mori, and identified 7 new tandem copies of BmGr30 on chromosome 7, bringing the total number of BmGrs to 76. Among these, we mapped 68 genes to chromosomes in a newly constructed chromosome distribution map and 8 genes to scaffolds; we also found new evidence for large clusters of BmGrs, especially from the bitter receptor family. RNA-seq analysis of diverse BmGr expression patterns in chemosensory organs of larvae and adults enabled us to draw a precise organ specific map of BmGr expression. Interestingly, most of the clustered genes were expressed in the same tissues and more than half of the genes were expressed in larval maxillae, larval thoracic legs and adult legs. For example, BmGr63 showed high expression levels in all organs in both larval and adult stages. By contrast, some genes showed expression limited to specific developmental stages or organs and tissues. BmGr19 was highly expressed in larval chemosensory organs (especially antennae and thoracic legs), the single exon genes BmGr53 and BmGr67 were expressed exclusively in larval tissues, the BmGr27-BmGr31 gene cluster on chr7 displayed a high expression level limited to adult legs and the candidate CO2 receptor BmGr2 was highly expressed in adult antennae, where few other Grs were expressed. Transcriptional analysis of the Grs in B. mori provides a valuable new reference for finding genes involved in plant-insect interactions in Lepidoptera and establishing correlations between these genes and vital insect behaviors like host plant selection and courtship for mating.

Genomic adaptation to polyphagy and insecticides in a major East Asian noctuid pest.

The tobacco cutworm, Spodoptera litura, is among the most widespread and destructive agricultural pests, feeding on over 100 crops throughout tropical and subtropical Asia. By genome sequencing, physical mapping and transcriptome analysis, we found that the gene families encoding receptors for bitter or toxic substances and detoxification enzymes, such as cytochrome P450, carboxylesterase and glutathione-S-transferase, were massively expanded in this polyphagous species, enabling its extraordinary ability to detect and detoxify many plant secondary compounds. Larval exposure to insecticidal toxins induced expression of detoxification genes, and knockdown of representative genes using short interfering RNA (siRNA) reduced larval survival, consistent with their contribution to the insect's natural pesticide tolerance. A population genetics study indicated that this species expanded throughout southeast Asia by migrating along a South India-South China-Japan axis, adapting to wide-ranging ecological conditions with diverse host plants and insecticides, surviving and adapting with the aid of its expanded detoxification systems. The findings of this study will enable the development of new pest management strategies for the control of major agricultural pests such as S. litura.

Transcriptome analysis of interactions between silkworm and cytoplasmic polyhedrosis virus.

Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) specifically infects silkworm midgut (MG) and multiplication occurs mainly in posterior midgut (PM). In this study, MG and fat body (FB) were extracted at 0, 3, 24, and 72 h after BmCPV infection. The total sequence reads of each sample were more than 1510000, and the mapping ratio exceeded 95.3%. Upregulated transcripts increased in MG during the infection process. Gene ontology (GO) categories showed that antioxidants were all upregulated in FB but not in MG. BGI001299, BGI014434, BGI012068, and BGI009201 were MG-specific genes with transmembrane transport function, the expression of which were induced by BmCPV. BGI001299, BGI014434, and BGI012068 expressed in entire MG and may be involved in BmCPV invasion. BGI009201 expressed only in PM and may be necessary for BmCPV proliferation. BmPGRP-S2 and BGI012452 (a putative serine protease) were induced by BmCPV and may be involved in immune defense against BmCPV. The expression level of BmCPV S1, S2, S3, S6, and S7 was high and there was no expression of S9 in MG 72 h, implying that the expression time of structural protein coding genes is earlier. These results provide insights into the mechanism of BmCPV infection and host defense.

Identification and characterization of toxins in the venom gland of the Chinese bird spider, Haplopelma hainanum, by transcriptomic analysis.

Tarantula venoms provide a model system for studying toxin selectivity, structure-activity relationships and molecular evolution of peptide toxins. Previous studies have identified a large number of peptide toxins in the venom of the Chinese bird spider Haplopelma hainanum, generally regarded as a highly venomous spider. However, the lack of available RNA-seq transcriptomic and genomic data is an obstacle to understanding its venom at the molecular level. In this study, we investigated the venom gland transcriptome of H. hainanum by RNA-seq, in the absence of an available genomic sequence. We identified 201 potential toxins among 57 181 de novo assembled transcripts, including knottins, Kunitz-type toxins, enzymes and other proteins. We systematically identified most of the knottins and Kunitz-type toxins, some of which showed strongly biased expression in the venom gland, including members of the huwentoxin-1, huwentoxin-2 and magi-1 families. We also discovered several novel potential toxins. These data demonstrate the high molecular and structural diversity in the venom toxins of H. hainanum. This study offers a useful strategy for exploring the complex components of spider venoms.

A comprehensive analysis of the chorion locus in silkmoth.

Despite more than 40 years of intense study, essential features of the silkmoth chorion (eggshell) are still not fully understood. To determine the precise structure of the chorion locus, we performed extensive EST analysis, constructed a bacterial artificial chromosome (BAC) contig, and obtained a continuous genomic sequence of 871,711 base pairs. We annotated 127 chorion genes in two segments interrupted by a 164 kb region with 5 non-chorion genes, orthologs of which were on chorion bearing scaffolds in 4 ditrysian families. Detailed transcriptome analysis revealed expression throughout choriogenesis of most chorion genes originally categorized as "middle", and evidence for diverse regulatory mechanisms including cis-elements, alternative splicing and promoter utilization, and antisense RNA. Phylogenetic analysis revealed multigene family associations and faster evolution of early chorion genes and transcriptionally active pseudogenes. Proteomics analysis identified 99 chorion proteins in the eggshell and micropyle localization of 1 early and 6 Hc chorion proteins.

Construction, complete sequence, and annotation of a BAC contig covering the silkworm chorion locus.

The silkmoth chorion was studied extensively by F.C. Kafatos' group for almost 40 years. However, the complete structure of the chorion locus was not obtained in the genome sequence of Bombyx mori published in 2008 due to repetitive sequences, resulting in gaps and an incomplete view of the locus. To obtain the complete sequence of the chorion locus, expressed sequence tags (ESTs) derived from follicular epithelium cells were used as probes to screen a bacterial artificial chromosome (BAC) library. Seven BACs were selected to construct a contig which covered the whole chorion locus. By Sanger sequencing, we successfully obtained complete sequences of the chorion locus spanning 871,711 base pairs on chromosome 2, where we annotated 127 chorion genes. The dataset reported here will recruit more researchers to revisit one of the oldest model systems which has been used to study developmentally regulated gene expression. It also provides insights into egg development and fertilization mechanisms and is relevant to applications related to improvements in breeding procedures and transgenesis.

Complete Genome Sequence of Bacillus bombysepticus, a Pathogen Leading to Bombyx mori Black Chest Septicemia.

Bacillus bombysepticus is a Gram-positive spore-forming bacterium. Here, we announce the first complete genome sequence of this organism isolated from the cadavers of silkworm larvae that had been sick. The genome contains a single circular chromosome and a circular plasmid. Analyses of the B. bombysepticus genome will provide insights into its pathomechanisms and biology.