STAG2 mutations regulate 3D genome organization, chromatin loops, and Polycomb signaling in glioblastoma multiforme.

Inactivating mutations of genes encoding the cohesin complex are common in a wide range of human cancers. STAG2 is the most commonly mutated subunit. Here we report the impact of stable correction of endogenous, naturally occurring STAG2 mutations on gene expression, 3D genome organization, chromatin loops, and Polycomb signaling in glioblastoma multiforme (GBM). In two GBM cell lines, correction of their STAG2 mutations significantly altered the expression of ∼10% of all expressed genes. Virtually all the most highly regulated genes were negatively regulated by STAG2 (i.e., expressed higher in STAG2-mutant cells), and one of them-HEPH-was regulated by STAG2 in uncultured GBM tumors as well. While STAG2 correction had little effect on large-scale features of 3D genome organization (A/B compartments, TADs), STAG2 correction did alter thousands of individual chromatin loops, some of which controlled the expression of adjacent genes. Loops specific to STAG2-mutant cells, which were regulated by STAG1-containing cohesin complexes, were very large, supporting prior findings that STAG1-containing cohesin complexes have greater loop extrusion processivity than STAG2-containing cohesin complexes and suggesting that long loops may be a general feature of STAG2-mutant cancers. Finally, STAG2 mutation activated Polycomb activity leading to increased H3K27me3 marks, identifying Polycomb signaling as a potential target for therapeutic intervention in STAG2-mutant GBM tumors. Together, these findings illuminate the landscape of STAG2-regulated genes, A/B compartments, chromatin loops, and pathways in GBM, providing important clues into the largely still unknown mechanism of STAG2 tumor suppression.

The evolution of synaptic and cognitive capacity: Insights from the nervous system transcriptome of Aplysia.

The gastropod mollusk Aplysia is an important model for cellular and molecular neurobiological studies, particularly for investigations of molecular mechanisms of learning and memory. We developed an optimized assembly pipeline to generate an improved Aplysia nervous system transcriptome. This improved transcriptome enabled us to explore the evolution of cognitive capacity at the molecular level. Were there evolutionary expansions of neuronal genes between this relatively simple gastropod Aplysia (20,000 neurons) and Octopus (500 million neurons), the invertebrate with the most elaborate neuronal circuitry and greatest behavioral complexity? Are the tremendous advances in cognitive power in vertebrates explained by expansion of the synaptic proteome that resulted from multiple rounds of whole genome duplication in this clade? Overall, the complement of genes linked to neuronal function is similar between Octopus and Aplysia. As expected, a number of synaptic scaffold proteins have more isoforms in humans than in Aplysia or Octopus. However, several scaffold families present in mollusks and other protostomes are absent in vertebrates, including the Fifes, Lev10s, SOLs, and a NETO family. Thus, whereas vertebrates have more scaffold isoforms from select families, invertebrates have additional scaffold protein families not found in vertebrates. This analysis provides insights into the evolution of the synaptic proteome. Both synaptic proteins and synaptic plasticity evolved gradually, yet the last deuterostome-protostome common ancestor already possessed an elaborate suite of genes associated with synaptic function, and critical for synaptic plasticity.

The Hypervariable Tpr Multigene Family of Theileria Parasites, Defined by a Conserved, Membrane-Associated, C-Terminal Domain, Includes Several Copies with Defined Orthology Between Species.

Multigene families often play an important role in host-parasite interactions. One of the largest multigene families in Theileria parva, the causative agent of East Coast fever, is the T. parva repeat (Tpr) gene family. The function of the putative Tpr proteins remains unknown. The initial publication of the T. parva reference genome identified 39 Tpr family open reading frames (ORFs) sharing a conserved C-terminal domain. Twenty-eight of these are clustered in a central region of chromosome 3, termed the "Tpr locus", while others are dispersed throughout all four nuclear chromosomes. The Tpr locus contains three of the four assembly gaps remaining in the genome, suggesting the presence of additional, as yet uncharacterized, Tpr gene copies. Here, we describe the use of long-read sequencing to attempt to close the gaps in the reference assembly of T. parva (located among multigene families clusters), characterize the full complement of Tpr family ORFs in the T. parva reference genome, and evaluate their evolutionary relationship with Tpr homologs in other Theileria species. We identify three new Tpr family genes in the T. parva reference genome and show that sequence similarity among paralogs in the Tpr locus is significantly higher than between genes outside the Tpr locus. We also identify sequences homologous to the conserved C-terminal domain in five additional Theileria species. Using these sequences, we show that the evolution of this gene family involves conservation of a few orthologs across species, combined with gene gains/losses, and species-specific expansions.

Comparative Analysis of Genome of Ehrlichia sp. HF, a Model Bacterium to Study Fatal Human Ehrlichiosis.

BACKGROUND: The genus Ehrlichia consists of tick-borne obligatory intracellular bacteria that can cause deadly diseases of medical and agricultural importance. Ehrlichia sp. HF, isolated from Ixodes ovatus ticks in Japan [also referred to as I. ovatus Ehrlichia (IOE) agent], causes acute fatal infection in laboratory mice that resembles acute fatal human monocytic ehrlichiosis caused by Ehrlichia chaffeensis. As there is no small laboratory animal model to study fatal human ehrlichiosis, Ehrlichia sp. HF provides a needed disease model. However, the inability to culture Ehrlichia sp. HF and the lack of genomic information have been a barrier to advance this animal model. In addition, Ehrlichia sp. HF has several designations in the literature as it lacks a taxonomically recognized name. RESULTS: We stably cultured Ehrlichia sp. HF in canine histiocytic leukemia DH82 cells from the HF strain-infected mice, and determined its complete genome sequence. Ehrlichia sp. HF has a single double-stranded circular chromosome of 1,148,904 bp, which encodes 866 proteins with a similar metabolic potential as E. chaffeensis. Ehrlichia sp. HF encodes homologs of all virulence factors identified in E. chaffeensis, including 23 paralogs of P28/OMP-1 family outer membrane proteins, type IV secretion system apparatus and effector proteins, two-component systems, ankyrin-repeat proteins, and tandem repeat proteins. Ehrlichia sp. HF is a novel species in the genus Ehrlichia, as demonstrated through whole genome comparisons with six representative Ehrlichia species, subspecies, and strains, using average nucleotide identity, digital DNA-DNA hybridization, and core genome alignment sequence identity. CONCLUSIONS: The genome of Ehrlichia sp. HF encodes all known virulence factors found in E. chaffeensis, substantiating it as a model Ehrlichia species to study fatal human ehrlichiosis. Comparisons between Ehrlichia sp. HF and E. chaffeensis will enable identification of in vivo virulence factors that are related to host specificity, disease severity, and host inflammatory responses. We propose to name Ehrlichia sp. HF as Ehrlichia japonica sp. nov. (type strain HF), to denote the geographic region where this bacterium was initially isolated.

Treg-inducing capacity of genomic DNA of Bifidobacterium longum subsp. infantis.

Background: Allergic and autoimmune diseases comprise a group of inflammatory disorders caused by aberrant immune responses in which CD25+ forkhead box P3-positive regulatory T cells (Treg) cells that normally suppress inflammatory events are often poorly functioning. This has stimulated an intensive investigative effort to find ways of increasing Tregs as a method of therapy for these conditions. Commensal microbiota known to have health benefits in humans include the lactic acid-producing, probiotic bacteria B. longum subsp. infantis and Lactobacillus rhamnosus. Mechanistically, several mechanisms have been proposed to explain how probiotics may favorably affect host immunity, including the induction of Tregs. Analysis of emerging data from several laboratories, including our own, suggest that DNA methylation may be an important determinant of immune reactivity responsible for Treg induction. Although methylated CpG moieties in normal mammalian DNA are both noninflammatory and lack immunogenicity, unmethylated CpGs, found largely in microbial DNA, are immunostimulatory and display proinflammatory properties. Objective: We hypothesize that microbiota with more DNA methylation may potentiate Treg induction to a greater degree than microbiota with a lower content of methylation. The purpose of the present study was to test this hypothesis by studying the methylation status of whole genomic DNA (gDNA) and the Treg-inducing capacity of purified gDNA in each of the probiotic bacteria B. longum subsp. infantis and L. rhamnosus, and a pathogenic Escherichia coli strain B. Results: We showed that gDNA from B. longum subsp. infantis is a potent Treg inducer that displays a dose-dependent response pattern at a dose threshold of 20 µg of gDNA. No similar Treg-inducing responses were observed with the gDNA from L. rhamnosus or E. coli. We identified a unique CpG methylated motif in the gDNA sequencing of B. longum subsp. infantis which was not found in L. rhamnosus or E. coli strain B. Conclusion: Although the literature indicates that both B. longum subsp. infantis and L. rhamnosus strains contribute to health, our data suggest that they do so by different mechanisms. Further, because of its small molecular size, low cost, ease of synthesis, and unique Treg-inducing feature, this methylated CpG oligodeoxynucleotide (ODN) from B. longum would offer many attractive features for an ideal novel therapeutic vaccine candidate for the treatment of immunologic diseases, such as the allergic and autoimmune disorders, in which Treg populations are diminished.

Streptococcus pneumoniae in the heart subvert the host response through biofilm-mediated resident macrophage killing.

For over 130 years, invasive pneumococcal disease has been associated with the presence of extracellular planktonic pneumococci, i.e. diplococci or short chains in affected tissues. Herein, we show that Streptococcus pneumoniae that invade the myocardium instead replicate within cellular vesicles and transition into non-purulent biofilms. Pneumococci within mature cardiac microlesions exhibited salient biofilm features including intrinsic resistance to antibiotic killing and the presence of an extracellular matrix. Dual RNA-seq and subsequent principal component analyses of heart- and blood-isolated pneumococci confirmed the biofilm phenotype in vivo and revealed stark anatomical site-specific differences in virulence gene expression; the latter having major implications on future vaccine antigen selection. Our RNA-seq approach also identified three genomic islands as exclusively expressed in vivo. Deletion of one such island, Region of Diversity 12, resulted in a biofilm-deficient and highly inflammogenic phenotype within the heart; indicating a possible link between the biofilm phenotype and a dampened host-response. We subsequently determined that biofilm pneumococci released greater amounts of the toxin pneumolysin than did planktonic or RD12 deficient pneumococci. This allowed heart-invaded wildtype pneumococci to kill resident cardiac macrophages and subsequently subvert cytokine/chemokine production and neutrophil infiltration into the myocardium. This is the first report for pneumococcal biofilm formation in an invasive disease setting. We show that biofilm pneumococci actively suppress the host response through pneumolysin-mediated immune cell killing. As such, our findings contradict the emerging notion that biofilm pneumococci are passively immunoquiescent.

The Streptococcus agalactiae Stringent Response Enhances Virulence and Persistence in Human Blood.

Streptococcus agalactiae (group B Streptococcus [GBS]) causes serious infections in neonates. We previously reported a transposon sequencing (Tn-seq) system for performing genomewide assessment of gene fitness in GBS. In order to identify molecular mechanisms required for GBS to transition from a mucosal commensal lifestyle to bloodstream invasion, we performed Tn-seq on GBS strain A909 with human whole blood. Our analysis identified 16 genes conditionally essential for GBS survival in blood, of which 75% were members of the capsular polysaccharide (cps) operon. Among the non-cps genes identified as conditionally essential was relA, which encodes an enzyme whose activity is central to the bacterial stringent response-a conserved adaptation to environmental stress. We used blood coincubation studies of targeted knockout strains to confirm the expected growth defects of GBS deficient in capsule or stringent response activation. Unexpectedly, we found that the relA knockout strains demonstrated decreased expression of ß-hemolysin/cytolysin, an important cytotoxin implicated in facilitating GBS invasion. Furthermore, chemical activation of the stringent response with serine hydroxamate increased ß-hemolysin/cytolysin expression. To establish a mechanism by which the stringent response leads to increased cytotoxicity, we performed transcriptome sequencing (RNA-seq) on two GBS strains grown under stringent response or control conditions. This revealed a conserved decrease in the expression of genes in the arginine deiminase pathway during stringent response activation. Through coincubation with supplemental arginine and the arginine antagonist canavanine, we show that arginine availability is a determinant of GBS cytotoxicity and that the pathway between stringent response activation and increased virulence is arginine dependent.

A comparative analysis of library prep approaches for sequencing low input translatome samples.

BACKGROUND: Cell type-specific ribosome-pulldown has become an increasingly popular method for analysis of gene expression. It allows for expression analysis from intact tissues and monitoring of protein synthesis in vivo. However, while its utility has been assessed, technical aspects related to sequencing of these samples, often starting with a smaller amount of RNA, have not been reported. In this study, we evaluated the performance of five library prep protocols for ribosome-associated mRNAs when only 250 pg-4 ng of total RNA are used. RESULTS: We obtained total and RiboTag-IP RNA, in three biological replicates. We compared 5 methods of library preparation for Illumina Next Generation sequencing: NuGEN Ovation RNA-Seq system V2 Kit, TaKaRa SMARTer Stranded Total RNA-Seq Kit, TaKaRa SMART-Seq v4 Ultra Low Input RNA Kit, Illumina TruSeq RNA Library Prep Kit v2 and NEBNext® Ultra™ Directional RNA Library Prep Kit using slightly modified protocols each with 4 ng of total RNA. An additional set of samples was processed using the TruSeq kit with 70 ng, as a 'gold standard' control and the SMART-Seq v4 with 250 pg of total RNA. TruSeq-processed samples had the best metrics overall, with similar results for the 4 ng and 70 ng samples. The results of the SMART-Seq v4 processed samples were similar to TruSeq (Spearman correlation > 0.8) despite using lower amount of input RNA. All RiboTag-IP samples had an increase in the intronic reads compared with the corresponding whole tissue, suggesting that the IP captures some immature mRNAs. The SMARTer-processed samples had a higher representation of ribosomal and non-coding RNAs leading to lower representation of protein coding mRNA. The enrichment or depletion of IP samples compared to corresponding input RNA was similar across all kits except for SMARTer kit. CONCLUSION: RiboTag-seq can be performed successfully with as little as 250 pg of total RNA when using the SMART-Seq v4 kit and 4 ng when using the modified protocols of other library preparation kits. The SMART-Seq v4 and TruSeq kits resulted in the highest quality libraries. RiboTag IP RNA contains some immature transcripts.

Neuraminidase A-Exposed Galactose Promotes Streptococcus pneumoniae Biofilm Formation during Colonization.

Streptococcus pneumoniae is an opportunistic pathogen that colonizes the nasopharynx. Herein we show that carbon availability is distinct between the nasopharynx and bloodstream of adult humans: glucose is absent from the nasopharynx, whereas galactose is abundant. We demonstrate that pneumococcal neuraminidase A (NanA), which cleaves terminal sialic acid residues from host glycoproteins, exposed galactose on the surface of septal epithelial cells, thereby increasing its availability during colonization. We observed that S. pneumoniae mutants deficient in NanA and ß-galactosidase A (BgaA) failed to form biofilms in vivo despite normal biofilm-forming abilities in vitro Subsequently, we observed that glucose, sucrose, and fructose were inhibitory for biofilm formation, whereas galactose, lactose, and low concentrations of sialic acid were permissive. Together these findings suggested that the genes involved in biofilm formation were under some form of carbon catabolite repression (CCR), a regulatory network in which genes involved in the uptake and metabolism of less-preferred sugars are silenced during growth with preferred sugars. Supporting this notion, we observed that a mutant deficient in pyruvate oxidase, which converts pyruvate to acetyl-phosphate under non-CCR-inducing growth conditions, was unable to form biofilms. Subsequent comparative transcriptome sequencing (RNA-seq) analyses of planktonic and biofilm-grown pneumococci showed that metabolic pathways involving the conversion of pyruvate to acetyl-phosphate and subsequently leading to fatty acid biosynthesis were consistently upregulated during diverse biofilm growth conditions. We conclude that carbon availability in the nasopharynx impacts pneumococcal biofilm formation in vivo Additionally, biofilm formation involves metabolic pathways not previously appreciated to play an important role.

The essential genome of Streptococcus agalactiae.

BACKGROUND: Next-generation sequencing of transposon-genome junctions from a saturated bacterial mutant library (Tn-seq) is a powerful tool that permits genome-wide determination of the contribution of genes to fitness of the organism under a wide range of experimental conditions. We report development, testing, and results from a Tn-seq system for use in Streptococcus agalactiae (group B Streptococcus; GBS), an important cause of neonatal sepsis. METHODS: Our method uses a Himar1 mini-transposon that inserts at genomic TA dinucleotide sites, delivered to GBS on a temperature-sensitive plasmid that is subsequently cured from the bacterial population. In order to establish the GBS essential genome, we performed Tn-seq on DNA collected from three independent mutant libraries-with at least 135,000 mutants per library-at serial 24 h time points after outgrowth in rich media. RESULTS: After statistical analysis of transposon insertion density and distribution, we identified 13.5 % of genes as essential and 1.2 % as critical, with high levels of reproducibility. Essential and critical genes are enriched for fundamental cellular housekeeping functions, such as acyl-tRNA biosynthesis, nucleotide metabolism, and glycolysis. We further validated our system by comparing fitness assignments of homologous genes in GBS and a close bacterial relative, Streptococcus pyogenes, which demonstrated 93 % concordance. Finally, we used our fitness assignments to identify signal transduction pathway components predicted to be essential or critical in GBS. CONCLUSIONS: We believe that our baseline fitness assignments will be a valuable tool for GBS researchers and that our system has the potential to reveal key pathogenesis gene networks and potential therapeutic/preventative targets.

Drosophila anti-nematode and antibacterial immune regulators revealed by RNA-Seq.

BACKGROUND: Drosophila melanogaster activates a variety of immune responses against microbial infections. However, information on the Drosophila immune response to entomopathogenic nematode infections is currently limited. The nematode Heterorhabditis bacteriophora is an insect parasite that forms a mutualistic relationship with the gram-negative bacteria Photorhabdus luminescens. Following infection, the nematodes release the bacteria that quickly multiply within the insect and produce several toxins that eventually kill the host. Although we currently know that the insect immune system interacts with Photorhabdus, information on interaction with the nematode vector is scarce. RESULTS: Here we have used next generation RNA-sequencing to analyze the transcriptional profile of wild-type adult flies infected by axenic Heterorhabditis nematodes (lacking Photorhabdus bacteria), symbiotic Heterorhabditis nematodes (carrying Photorhabdus bacteria), and Photorhabdus bacteria alone. We have obtained approximately 54 million reads from the different infection treatments. Bioinformatic analysis shows that infection with Photorhabdus alters the transcription of a large number of Drosophila genes involved in translational repression as well in response to stress. However, Heterorhabditis infection alters the transcription of several genes that participate in lipidhomeostasis and metabolism, stress responses, DNA/protein synthesis and neuronal functions. We have also identified genes in the fly with potential roles in nematode recognition, anti-nematode activity and nociception. CONCLUSIONS: These findings provide fundamental information on the molecular events that take place in Drosophila upon infection with the two pathogens, either separately or together. Such large-scale transcriptomic analyses set the stage for future functional studies aimed at identifying the exact role of key factors in the Drosophila immune response against nematode-bacteria complexes.

Genetic variation in vitro and in vivo of an attenuated Lassa vaccine candidate.

The attenuated Lassa vaccine candidate ML29 is a laboratory-produced reassortant between Lassa and Mopeia viruses, two Old World arenaviruses that differ by 40% in nucleic acid sequence. In our previous studies, ML29 elicited sterilizing immunity against Lassa virus challenge in guinea pigs and marmosets and virus-specific cell-mediated immunity in both simian immunodeficiency virus (SIV)-infected and uninfected rhesus macaques. Here, we show that ML29 is stable after 12 passages in vitro without losing its plaque morphology or its attenuated phenotype in suckling mice. Additionally, we used deep sequencing to characterize the viral population comprising the original stock of ML29, the stock of ML29 after 12 passages in Vero cells, and the ML29 isolates obtained from vaccinated animals. Twenty-seven isolates bore approximately 77 mutations that exceeded 20% of the single-nucleotide polymorphism (SNP) changes at any single locus. Of these 77 mutations, 5 appeared to be host specific, for example, appearing in mice but not in primates. None of these mutations were reversions of ML29 to the sequences of the parental Lassa and Mopeia viruses. The host-specific mutations indicate viral adaptations to virus-host interactions, and such interactions make reasonable targets for antiviral approaches. Variants capable of chronic infection did not emerge from any of the primate infections, even in immune-deficient animals, indicating that the ML29 reassortant is reasonably stable in vivo. In conclusion, the preclinical studies of ML29 as a Lassa virus vaccine candidate have been advanced, showing high levels of protection in nonhuman primates and acceptable stability both in vitro and in vivo.

Genomic diversity of 2010 Haitian cholera outbreak strains.

The millions of deaths from cholera during the past 200 y, coupled with the morbidity and mortality of cholera in Haiti since October 2010, are grim reminders that Vibrio cholerae, the etiologic agent of cholera, remains a scourge. We report the isolation of both V. cholerae O1 and non-O1/O139 early in the Haiti cholera epidemic from samples collected from victims in 18 towns across eight Arrondissements of Haiti. The results showed two distinct populations of V. cholerae coexisted in Haiti early in the epidemic. As non-O1/O139 V. cholerae was the sole pathogen isolated from 21% of the clinical specimens, its role in this epidemic, either alone or in concert with V. cholerae O1, cannot be dismissed. A genomic approach was used to examine similarities and differences among the Haitian V. cholerae O1 and V. cholerae non-O1/O139 strains. A total of 47 V. cholerae O1 and 29 V. cholerae non-O1/O139 isolates from patients and the environment were sequenced. Comparative genome analyses of the 76 genomes and eight reference strains of V. cholerae isolated in concurrent epidemics outside Haiti and 27 V. cholerae genomes available in the public database demonstrated substantial diversity of V. cholerae and ongoing flux within its genome.

Extensive duplication of the Wolbachia DNA in chromosome four of Drosophila ananassae.

BACKGROUND: Lateral gene transfer (LGT) from bacterial Wolbachia endosymbionts has been detected in ~20% of arthropod and nematode genome sequencing projects. Many of these transfers are large and contain a substantial part of the Wolbachia genome. RESULTS: Here, we re-sequenced three D. ananassae genomes from Asia and the Pacific that contain large LGTs from Wolbachia. We find that multiple copies of the Wolbachia genome are transferred to the Drosophila nuclear genome in all three lines. In the D. ananassae line from Indonesia, the copies of Wolbachia DNA in the nuclear genome are nearly identical in size and sequence yielding an even coverage of mapped reads over the Wolbachia genome. In contrast, the D. ananassae lines from Hawaii and India show an uneven coverage of mapped reads over the Wolbachia genome suggesting that different parts of these LGTs are present in different copy numbers. In the Hawaii line, we find that this LGT is underrepresented in third instar larvae indicative of being heterochromatic. Fluorescence in situ hybridization of mitotic chromosomes confirms that the LGT in the Hawaii line is heterochromatic and represents ~20% of the sequence on chromosome 4 (dot chromosome, Muller element F). CONCLUSIONS: This collection of related lines contain large lateral gene transfers composed of multiple Wolbachia genomes that constitute >2% of the D. ananassae genome (~5 Mbp) and partially explain the abnormally large size of chromosome 4 in D. ananassae.

Single molecule sequencing and genome assembly of a clinical specimen of Loa loa, the causative agent of loiasis.

BACKGROUND: More than 20% of the world's population is at risk for infection by filarial nematodes and >180 million people worldwide are already infected. Along with infection comes significant morbidity that has a socioeconomic impact. The eight filarial nematodes that infect humans are Wuchereria bancrofti, Brugia malayi, Brugia timori, Onchocerca volvulus, Loa loa, Mansonella perstans, Mansonella streptocerca, and Mansonella ozzardi, of which three have published draft genome sequences. Since all have humans as the definitive host, standard avenues of research that rely on culturing and genetics have often not been possible. Therefore, genome sequencing provides an important window into understanding the biology of these parasites. The need for large amounts of high quality genomic DNA from homozygous, inbred lines; the availability of only short sequence reads from next-generation sequencing platforms at a reasonable expense; and the lack of random large insert libraries has limited our ability to generate high quality genome sequences for these parasites. However, the Pacific Biosciences single molecule, real-time sequencing platform holds great promise in reducing input amounts and generating sufficiently long sequences that bypass the need for large insert paired libraries. RESULTS: Here, we report on efforts to generate a more complete genome assembly for L. loa using genetically heterogeneous DNA isolated from a single clinical sample and sequenced on the Pacific Biosciences platform. To obtain the best assembly, numerous assemblers and sequencing datasets were analyzed, combined, and compared. Quiver-informed trimming of an assembly of only Pacific Biosciences reads by HGAP2 was selected as the final assembly of 96.4 Mbp in 2,250 contigs. This results in ~9% more of the genome in ~85% fewer contigs from ~80% less starting material at a fraction of the cost of previous Roche 454-based sequencing efforts. CONCLUSIONS: The result is the most complete filarial nematode assembly produced thus far and demonstrates the utility of single molecule sequencing on the Pacific Biosciences platform for genetically heterogeneous metazoan genomes.

Rapid transcriptome sequencing of an invasive pest, the brown marmorated stink bug Halyomorpha halys.

BACKGROUND: Halyomorpha halys (Stål) (Insecta:Hemiptera;Pentatomidae), commonly known as the Brown Marmorated Stink Bug (BMSB), is an invasive pest of the mid-Atlantic region of the United States, causing economically important damage to a wide range of crops. Native to Asia, BMSB was first observed in Allentown, PA, USA, in 1996, and this pest is now well-established throughout the US mid-Atlantic region and beyond. In addition to the serious threat BMSB poses to agriculture, BMSB has become a nuisance to homeowners, invading home gardens and congregating in large numbers in human-made structures, including homes, to overwinter. Despite its significance as an agricultural pest with limited control options, only 100 bp of BMSB sequence data was available in public databases when this project began. RESULTS: Transcriptome sequencing was undertaken to provide a molecular resource to the research community to inform the development of pest control strategies and to provide molecular data for population genetics studies of BMSB. Using normalized, strand-specific libraries, we sequenced pools of all BMSB life stages on the Illumina HiSeq. Trinity was used to assemble 200,000 putative transcripts in >100,000 components. A novel bioinformatic method that analyzed the strand-specificity of the data reduced this to 53,071 putative transcripts from 18,573 components. By integrating multiple other data types, we narrowed this further to 13,211 representative transcripts. CONCLUSIONS: Bacterial endosymbiont genes were identified in this dataset, some of which have a copy number consistent with being lateral gene transfers between endosymbiont genomes and Hemiptera, including ankyrin-repeat related proteins, lysozyme, and mannanase. Such genes and endosymbionts may provide novel targets for BMSB-specific biocontrol. This study demonstrates the utility of strand-specific sequencing in generating shotgun transcriptomes and that rapid sequencing shotgun transcriptomes is possible without the need for extensive inbreeding to generate homozygous lines. Such sequencing can provide a rapid response to pest invasions similar to that already described for disease epidemiology.

GAGE-B: an evaluation of genome assemblers for bacterial organisms.

MOTIVATION: A large and rapidly growing number of bacterial organisms have been sequenced by the newest sequencing technologies. Cheaper and faster sequencing technologies make it easy to generate very high coverage of bacterial genomes, but these advances mean that DNA preparation costs can exceed the cost of sequencing for small genomes. The need to contain costs often results in the creation of only a single sequencing library, which in turn introduces new challenges for genome assembly methods. RESULTS: We evaluated the ability of multiple genome assembly programs to assemble bacterial genomes from a single, deep-coverage library. For our comparison, we chose bacterial species spanning a wide range of GC content and measured the contiguity and accuracy of the resulting assemblies. We compared the assemblies produced by this very high-coverage, one-library strategy to the best assemblies created by two-library sequencing, and we found that remarkably good bacterial assemblies are possible with just one library. We also measured the effect of read length and depth of coverage on assembly quality and determined the values that provide the best results with current algorithms. CONTACT: salzberg@jhu.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A Longitudinal Metagenomic Comparative Analysis of Oral Microbiome Shifts in Patients Receiving Proton Radiation versus Photon Radiation for Head and Neck Cancer.

Introduction: Due to the radiation-sparing effects on salivary gland acini, changes in the composition of the oral microbiome may be a driver for improved outcomes in patients receiving proton radiation, with potentially worse outcomes in patients exposed to photon radiation therapy. To date, a head-to-head comparison of oral microbiome changes at a metagenomic level with longitudinal sampling has yet to be performed in these patient cohorts. Methods and Materials: To comparatively analyze oral microbiome shifts during head and neck radiation therapy, a prospective pilot cohort study was performed at the Maryland Proton Treatment Center and the University of Maryland Marlene and Stewart Greenebaum Comprehensive Cancer Center. A longitudinal metagenomic comparative analysis of oral microbiome shifts was performed at three time points (pre-radiation, during radiation, and immediately post-radiation). Head and neck cancer patients receiving proton radiation (n = 4) were compared to photon radiation (n = 4). Additional control groups included healthy age- and sex-matched controls (n = 5), head and neck cancer patients who never received radiation therapy (n = 8), and patients with oral inflammatory disease (n = 3). Results: Photon therapy patients presented with lower microbial alpha diversity at all timepoints, and there was a trend towards reduced species richness as compared with proton therapy. Healthy controls and proton patients exhibited overall higher and similar diversity. A more dysbiotic state was observed in patients receiving photon therapy as compared to proton therapy, in which oral microbial homeostasis was maintained. Mucositis was observed in 3/4 photon patients and was not observed in any proton patients during radiation therapy. The bacterial de novo pyrimidine biosynthesis pathway and the nitrate reduction V pathway were comparatively higher following photon exposure. These functional changes in bacterial metabolism may suggest that photon exposure produces a more permissive environment for the proliferation of pathogenic bacteria. Conclusion: Oral microbiome dysbiosis in patients receiving photon radiation may be associated with increased mucositis occurrence. Proton radiation therapy for head and neck cancer demonstrates a safer side effect profile in terms of oral complications, oral microbiome dysbiosis, and functional metabolic status.

Extensively duplicated and transcriptionally active recent lateral gene transfer from a bacterial Wolbachia endosymbiont to its host filarial nematode Brugia malayi.

BACKGROUND: Lymphatic filariasis is a neglected tropical disease afflicting more than 120 million people, while another 1.3 billion people are at risk of infection. The nematode worm Brugia malayi is one of the causative agents of the disease and exists in a mutualistic symbiosis with Wolbachia bacteria. Since extensive lateral gene transfer occurs frequently between Wolbachia and its hosts, we sought to measure the extent of such LGT in B. malayi by whole genome sequencing of Wolbachia-depleted worms. RESULTS: A considerable fraction (at least 115.4-kbp, or 10.6%) of the 1.08-Mbp Wolbachia wBm genome has been transferred to its nematode host and retains high levels of similarity, including 227 wBm genes and gene fragments. Complete open reading frames were transferred for 32 of these genes, meaning they have the potential to produce functional proteins. Moreover, four transfers have evidence of life stage-specific regulation of transcription at levels similar to other nematode transcripts, strengthening the possibility that they are functional. CONCLUSIONS: There is extensive and ongoing transfer of Wolbachia DNA to the worm genome and some transfers are transcribed in a stage-specific manner at biologically relevant levels.

Genome sequencing and comparative genomics of honey bee microsporidia, Nosema apis reveal novel insights into host-parasite interactions.

BACKGROUND: The microsporidia parasite Nosema contributes to the steep global decline of honey bees that are critical pollinators of food crops. There are two species of Nosema that have been found to infect honey bees, Nosema apis and N. ceranae. Genome sequencing of N. apis and comparative genome analysis with N. ceranae, a fully sequenced microsporidia species, reveal novel insights into host-parasite interactions underlying the parasite infections. RESULTS: We applied the whole-genome shotgun sequencing approach to sequence and assemble the genome of N. apis which has an estimated size of 8.5 Mbp. We predicted 2,771 protein- coding genes and predicted the function of each putative protein using the Gene Ontology. The comparative genomic analysis led to identification of 1,356 orthologs that are conserved between the two Nosema species and genes that are unique characteristics of the individual species, thereby providing a list of virulence factors and new genetic tools for studying host-parasite interactions. We also identified a highly abundant motif in the upstream promoter regions of N. apis genes. This motif is also conserved in N. ceranae and other microsporidia species and likely plays a role in gene regulation across the microsporidia. CONCLUSIONS: The availability of the N. apis genome sequence is a significant addition to the rapidly expanding body of microsprodian genomic data which has been improving our understanding of eukaryotic genome diversity and evolution in a broad sense. The predicted virulent genes and transcriptional regulatory elements are potential targets for innovative therapeutics to break down the life cycle of the parasite.