Search

Mutualistic intracellular symbiosis between bacteria and insects is a widespread phenomenon that has contributed to the global success of insects. The symbionts, by provisioning nutrients lacking from diets, allow various insects to occupy or dominate ecological niches that might otherwise be unavailable. One such insect is the glassy-winged sharpshooter (Homalodisca coagulata), which feeds on xylem fluid, a diet exceptionally poor in organic nutrients. Phylogenetic studies based on rRNA have shown two types of bacterial symbionts to be coevolving with sharpshooters: the gamma-proteobacterium Baumannia cicadellinicola and the Bacteroidetes species Sulcia muelleri. We report here the sequencing and analysis of the 686,192-base pair genome of B. cicadellinicola and approximately 150 kilobase pairs of the small genome of S. muelleri, both isolated from H. coagulata. Our study, which to our knowledge is the first genomic analysis of an obligate symbiosis involving multiple partners, suggests striking complementarity in the biosynthetic capabilities of the two symbionts: B. cicadellinicola devotes a substantial portion of its genome to the biosynthesis of vitamins and cofactors required by animals and lacks most amino acid biosynthetic pathways, whereas S. muelleri apparently produces most or all of the essential amino acids needed by its host. This finding, along with other results of our genome analysis, suggests the existence of metabolic codependency among the two unrelated endosymbionts and their insect host. This dual symbiosis provides a model case for studying correlated genome evolution and genome reduction involving multiple organisms in an intimate, obligate mutualistic relationship. In addition, our analysis provides insight for the first time into the differences in symbionts between insects (e.g., aphids) that feed on phloem versus those like H. coagulata that feed on xylem. Finally, the genomes of these two symbionts provide potential targets for controlling plant pathogens such as Xylella fastidiosa, a major agroeconomic problem, for which H. coagulata and other sharpshooters serve as vectors of transmission. © 2006 Wu et al.

We report here the sequencing and analysis of the genome of the thermophilic bacterium Carboxydothermus hydrogenoformans Z-2901. This species is a model for studies of hydrogenogens, which are diverse bacteria and archaea that grow anaerobically utilizing carbon monoxide (CO) as their sole carbon source and water as an electron acceptor, producing carbon dioxide and hydrogen as waste products. Organisms that make use of CO do so through carbon monoxide dehydrogenase complexes. Remarkably, analysis of the genome of C. hydrogenoformans reveals th e presence of at least five highly differentiated anaerobic carbon monoxide dehydrogenase complexes, which may in part explain how this species is able to grow so much more rapidly on CO than many other species. Analysis of the genome also has provided many general insights into the metabolism of this organism which should make it easier to use it as a source of biologically produced hydrogen gas. One surprising finding is the presence of many genes previously found only in sporulating species in the Firmicutes Phylum. Although this species is also a Firmicutes, it was not known to sporulate previously. Here we show that it does sporulate and because it is missing many of the genes involved in sporulation in other species, this organism may serve as a ''minimal'' model for sporulation studies. In addition, using phylogenetic profile analysis, we have identified many uncharacterized gene families found in all known sporulating Firmicutes, but not in any non-sporulating bacteria, including a sigma factor not known to be involved in sporulation previously. Copyright: ©2005 Wu et al.

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.

"Leaky gut," or high intestinal barrier permeability, is common in preterm newborns. The role of the microbiota in this process remains largely uncharacterized. We employed both short- and long-read sequencing of the 16S rRNA gene and metagenomes to characterize the intestinal microbiome of a longitudinal cohort of 113 preterm infants born between 24 0/7 and 32 6/7  weeks of gestation. Enabled by enhanced taxonomic resolution, we found that a significantly increased abundance of Bifidobacterium breve and a diet rich in mother's breastmilk were associated with intestinal barrier maturation during the first week of life. We combined these factors using genome-resolved metagenomics and identified a highly specialized genetic capability of the Bifidobacterium strains to assimilate human milk oligosaccharides and host-derived glycoproteins. Our study proposes mechanistic roles of breastmilk feeding and intestinal microbial colonization in postnatal intestinal barrier maturation; these observations are critical toward advancing therapeutics to prevent and treat hyperpermeable gut-associated conditions, including necrotizing enterocolitis (NEC). IMPORTANCE Despite improvements in neonatal intensive care, necrotizing enterocolitis (NEC) remains a leading cause of morbidity and mortality. "Leaky gut," or intestinal barrier immaturity with elevated intestinal permeability, is the proximate cause of susceptibility to NEC. Early detection and intervention to prevent leaky gut in "at-risk" preterm neonates are critical for decreasing the risk of potentially life-threatening complications like NEC. However, the complex interactions between the developing gut microbial community, nutrition, and intestinal barrier function remain largely uncharacterized. In this study, we reveal the critical role of a sufficient breastmilk feeding volume and the specialized carbohydrate metabolism capability of Bifidobacterium in the coordinated postnatal improvement of the intestinal barrier. Determining the clinical and microbial biomarkers that drive the intestinal developmental disparity will inform early detection and novel therapeutic strategies to promote appropriate intestinal barrier maturation and prevent NEC and other adverse health conditions in preterm infants.

Eukaryotic genomes can acquire bacterial DNA via lateral gene transfer (LGT). 1 A prominent source of LGT is Wolbachia, 2 a widespread endosymbiont of arthropods and nematodes that is transmitted maternally through female germline cells. 3 , 4 The DNA transfer from the Wolbachia endosymbiont wAna to Drosophila ananassae is extensive 5-7 and has been localized to chromosome 4, contributing to chromosome expansion in this lineage. 6 As has happened frequently with claims of bacteria-to-eukaryote LGT, the contribution of wAna transfers to the expanded size of D. ananassae chromosome 4 has been specifically contested 8 owing to an assembly where Wolbachia sequences were classified as contaminants and removed. 9 Here, long-read sequencing with DNA from a Wolbachia-cured line enabled assembly of 4.9 Mbp of nuclear Wolbachia transfers (nuwts) in D. ananassae and a 24-kbp nuclear mitochondrial transfer. The nuwts are <8,000 years old in at least two locations in chromosome 4 with at least one whole-genome integration followed by rapid extensive duplication of most of the genome with regions that have up to 10 copies. The genes in nuwts are accumulating small indels and mobile element insertions. Among the highly duplicated genes are cifA and cifB, two genes associated with Wolbachia-mediated Drosophila cytoplasmic incompatibility. The wAna strain that was the source of nuwts was subsequently replaced by a different wAna endosymbiont. Direct RNA Nanopore sequencing of Wolbachia-cured lines identified nuwt transcripts, including spliced transcripts, but functionality, if any, remains elusive., Competing Interests: Declaration of interests E.S.T. is currently affiliated with NCBI at the National Institutes of Health, Bethesda, MD, USA. M.G. is currently affiliated with Applied Physics Laboratory, Johns Hopkins University, Laurel, MD, USA. M.C. is currently affiliated with the National Institute for Allergy and Infectious Disease at the National Institutes of Health, Bethesda, MD, USA. B.C.S. is currently affiliated with AstraZeneca, Rockville, MD, USA., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

var genes encode Plasmodium falciparum erythrocyte membrane protein-1 (PfEMP1) antigens. These highly diverse antigens are displayed on the surface of infected erythrocytes and play a critical role in immune evasion and sequestration of infected erythrocytes. Studies of var expression using non-leukocyte-depleted blood are challenging because of the predominance of host genetic material and lack of conserved var segments. Our goal was to enrich for parasite RNA, allowing de novo assembly of var genes and detection of expressed novel variants. We used two overall approaches: (i) enriching for total mRNA in the sequencing library preparations and (ii) enriching for parasite RNA with a custom capture array based on Roche's SeqCap EZ enrichment system. The capture array was designed with probes based on the whole 3D7 reference genome and an additional >4,000 full-length var gene sequences from other P. falciparum strains. We tested each method on RNA samples from Malian children with severe or uncomplicated malaria infections. All reads mapping to the human genome were removed, the remaining reads were assembled de novo into transcripts, and from these, var -like transcripts were identified and annotated. The capture array produced the longest maximum length and largest numbers of var gene transcripts in each sample, particularly in samples with low parasitemia. Identifying the most-expressed var gene sequences in whole-blood clinical samples without the need for extensive processing or generating sample-specific reference genome data is critical for understanding the role of PfEMP1s in malaria pathogenesis. IMPORTANCE Malaria parasites display antigens on the surface of infected red blood cells in the human host that facilitate attachment to blood vessels, contributing to the severity of infection. These antigens are highly variable, allowing the parasite to evade the immune system. Identifying these expressed antigens is critical to understanding the development of severe malarial disease. However, clinical samples contain limited amounts of parasite genetic material, a challenge for sequencing efforts further compounded by the extreme diversity of the parasite surface antigens. We present a method that enriches for these antigen sequences in clinical samples using a custom capture array, requiring minimal processing in the field. While our results are focused on the malaria parasite Plasmodium falciparum, this approach has broad applicability to other highly diverse antigens from other parasites and pathogens such as those that cause giardiasis and leishmaniasis.

The sequence diversity of natural and laboratory populations of Brugia pahangi and Brugia malayi was assessed with Illumina resequencing followed by mapping in order to identify single nucleotide variants and insertions/deletions. In natural and laboratory Brugia populations, there is a lack of sequence diversity on chromosome X relative to the autosomes (πX/πA = 0.2), which is lower than the expected (πX/πA = 0.75). A reduction in diversity is also observed in other filarial nematodes with neo-X chromosome fusions in the genera Onchocerca and Wuchereria, but not those without neo-X chromosome fusions in the genera Loa and Dirofilaria. In the species with neo-X chromosome fusions, chromosome X is abnormally large, containing a third of the genetic material such that a sizable portion of the genome is lacking sequence diversity. Such profound differences in genetic diversity can be consequential, having been associated with drug resistance and adaptability, with the potential to affect filarial eradication., Competing Interests: The authors have declared that no competing interests exist.

Library preparation for high-throughput sequencing applications is a critical step in producing representative, unbiased sequencing data. The iGenomX Riptide High Throughput Rapid Library Prep Kit purports to provide high-quality sequencing data with lower costs compared to other Illumina library kits. To test these claims, we compared sequence data quality of Riptide libraries to libraries constructed with KAPA Hyper and NEBNext Ultra. Across several single-source genome samples, mapping performance and de novo assembly of Riptide libraries were similar to conventional libraries prepared with the same DNA. Poor performance of some libraries resulted in low sequencing depth. In particular, degraded DNA samples may be challenging to sequence with Riptide. There was little cross-well plate contamination with the overwhelming majority of reads belong to the proper source genomes. The sequencing of metagenome samples using different Riptide primer sets resulted in variable taxonomic assignment of reads. Increased adoption of the Riptide kit will decrease library preparation costs. However, this method might not be suitable for degraded DNA., (© 2021. The Author(s).)

The newest generation of DNA sequencing technology is highlighted by the ability to generate sequence reads hundreds of kilobases in length. Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) have pioneered competitive long read platforms, with more recent work focused on improving sequencing throughput and per-base accuracy. We used whole-genome sequencing data produced by three PacBio protocols (Sequel II CLR, Sequel II HiFi, RS II) and two ONT protocols (Rapid Sequencing and Ligation Sequencing) to compare assemblies of the bacteria Escherichia coli and the fruit fly Drosophila ananassae. In both organisms tested, Sequel II assemblies had the highest consensus accuracy, even after accounting for differences in sequencing throughput. ONT and PacBio CLR had the longest reads sequenced compared to PacBio RS II and HiFi, and genome contiguity was highest when assembling these datasets. ONT Rapid Sequencing libraries had the fewest chimeric reads in addition to superior quantification of E. coli plasmids versus ligation-based libraries. The quality of assemblies can be enhanced by adopting hybrid approaches using Illumina libraries for bacterial genome assembly or polishing eukaryotic genome assemblies, and an ONT-Illumina hybrid approach would be more cost-effective for many users. Genome-wide DNA methylation could be detected using both technologies, however ONT libraries enabled the identification of a broader range of known E. coli methyltransferase recognition motifs in addition to undocumented D. ananassae motifs. The ideal choice of long read technology may depend on several factors including the question or hypothesis under examination. No single technology outperformed others in all metrics examined., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

Long-read and strand-specific sequencing technologies together facilitate the de novo assembly of high-quality haplotype-resolved human genomes without parent-child trio data. We present 64 assembled haplotypes from 32 diverse human genomes. These highly contiguous haplotype assemblies (average minimum contig length needed to cover 50% of the genome: 26 million base pairs) integrate all forms of genetic variation, even across complex loci. We identified 107,590 structural variants (SVs), of which 68% were not discovered with short-read sequencing, and 278 SV hotspots (spanning megabases of gene-rich sequence). We characterized 130 of the most active mobile element source elements and found that 63% of all SVs arise through homology-mediated mechanisms. This resource enables reliable graph-based genotyping from short reads of up to 50,340 SVs, resulting in the identification of 1526 expression quantitative trait loci as well as SV candidates for adaptive selection within the human population., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Semaphorin 4D (Sema4D) is a glycoprotein that is expressed by several tumors and immune cells. It can function as a membrane bound protein or as a cleaved soluble protein (sSema4D). We sought to investigate the translational potential of plasma sSema4D as an immune marker in plasma of patients with head and neck squamous cell carcinoma (HNSCC). Paired peripheral blood and tumor tissue samples of 104 patients with HNSCC were collected at the same time point to allow for real time analysis. Scoring of the histological inflammatory subtype (HIS) was carried out using Sema4D immunohistochemistry on the tumor tissue. sSema4D was detected in plasma using direct ELISA assay. Defining elevated sSema4D as values above the 95 th percentile in healthy controls, our data showed that sSema4D levels in plasma were elevated in 25.0% (95% CI, 16.7-34.9%) of the patients with HNSCC and showed significant association with HIS immune excluded (HIS-IE) (p = 0.007), Sema4D +ve tumor cells (TCs) (p = 0.018) and PD-L1 +ve immune cells (ICs) (p = 0.038). A multi-variable logistic regression analysis showed that HIS was significantly (P = 0.004) associated with elevated sSema4D, an association not explained by available patient-level factors. Using the IO-360 nanoString platform, differential gene expression (DGE) analysis of 10 HNSCC tumor tissues showed that patients with high sSema4D in plasma (HsS4D) clustered as IFN- γ negative tumor immune signature and were mostly HIS-IE. The IC type in the HsS4D paired tumor tissue was predominantly myeloid, while the lymphoid compartment was higher in the low sSema4D (LsS4D). The Wnt signaling pathway was upregulated in the HsS4D group. Further analysis using the IO-360, 770 gene set, showed significant non-inflamed profile of the HsS4D tumors compared to the LsS4D. In conclusion, our data reveals an association between sSema4D and the histological inflammatory subtype., Competing Interests: TW is CEO of WebbCures, LLC, co-founder of Screen Therapeutics and serves as an advisor for Immunaccel Labs, LLC. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Younis, Ghita, Elnaggar, Chaisuparat, Theofilou, Dyalram, Ord, Davila, Tallon, Papadimitriou, Webb, Bentzen and Lubek.)

Quantification tools for RNA sequencing (RNA-Seq) analyses are often designed and tested using human transcriptomics data sets, in which full-length transcript sequences are well annotated. For prokaryotic transcriptomics experiments, full-length transcript sequences are seldom known, and coding sequences must instead be used for quantification steps in RNA-Seq analyses. However, operons confound accurate quantification of coding sequences since a single transcript does not necessarily equate to a single gene. Here, we introduce FADU (Feature Aggregate Depth Utility), a quantification tool designed specifically for prokaryotic RNA-Seq analyses. FADU assigns partial count values proportional to the length of the fragment overlapping the target feature. To assess the ability of FADU to quantify genes in prokaryotic transcriptomics analyses, we compared its performance to those of eXpress, featureCounts, HTSeq, kallisto, and Salmon across three paired-end read data sets of (i) Ehrlichia chaffeensis , (ii) Escherichia coli , and (iii) the Wolbachia endosymbiont w Bm. Across each of the three data sets, we find that FADU can more accurately quantify operonic genes by deriving proportional counts for multigene fragments within operons. FADU is available at https://github.com/IGS/FADU IMPORTANCE Most currently available quantification tools for transcriptomics analyses have been designed for human data sets, in which full-length transcript sequences, including the untranslated regions, are well annotated. In most prokaryotic systems, full-length transcript sequences have yet to be characterized, leading to prokaryotic transcriptomics analyses being performed based on only the coding sequences. In contrast to eukaryotes, prokaryotes contain polycistronic transcripts, and when genes are quantified based on coding sequences instead of transcript sequences, this leads to an increased abundance of improperly assigned ambiguous multigene fragments, specifically those mapping to multiple genes in operons. Here, we describe FADU, a quantification tool for prokaryotic RNA-Seq analyses designed to assign proportional counts with the purpose of better quantifying operonic genes while minimizing the pitfalls associated with improperly assigning fragment counts from ambiguous transcripts., (Copyright © 2021 Chung et al.)

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.

The 13,647-bp complete mitochondrial genome of Mansonella perstans was sequenced and is syntenic to the mitochondrial genome of Mansonella ozzardi Phylogenetic analysis of the mitochondrial genome is consistent with the known phylogeny of ONC5 group filarial nematodes., (Copyright © 2020 Chung et al.)

Brugia pahangi is a zoonotic parasite that is closely related to human-infecting filarial nematodes. Here, we report the nearly complete genome of Brugia pahangi , including assemblies of four autosomes and an X chromosome, with only seven gaps. The Y chromosome is still not completely assembled., (Copyright © 2020 Mattick et al.)

Lymphatic filariasis is a devastating disease caused by filarial nematode roundworms, which contain obligate Wolbachia endosymbionts. Here, we assembled the genome of w Bp, the Wolbachia endosymbiont of the filarial nematode Brugia pahangi , from Illumina, Pacific Biosciences, and Oxford Nanopore data. The complete, circular genome is 1,072,967 bp., (Copyright © 2020 Lebov et al.)

Bacterial vaginosis-associated bacterium 1 (BVAB1) is an as-yet uncultured bacterial species found in the human vagina that belongs to the family Lachnospiraceae within the order Clostridiales . As its name suggests, this bacterium is often associated with bacterial vaginosis (BV), a common vaginal disorder that has been shown to increase a woman's risk for HIV, Chlamydia trachomatis , and Neisseria gonorrhoeae infections as well as preterm birth. BVAB1 has been further associated with the persistence of BV following metronidazole treatment, increased vaginal inflammation, and adverse obstetrics outcomes. There is no available complete genome sequence of BVAB1, which has made it difficult to mechanistically understand its role in disease. We present here a circularized metagenome-assembled genome (cMAG) of BVAB1 as well as a comparative analysis including an additional six metagenome-assembled genomes (MAGs) of this species. These sequences were derived from cervicovaginal samples of seven separate women. The cMAG was obtained from a metagenome sequenced with long-read technology on a PacBio Sequel II instrument while the others were derived from metagenomes sequenced on the Illumina HiSeq platform. The cMAG is 1.649 Mb in size and encodes 1,578 genes. We propose to rename BVAB1 to " Candidatus Lachnocurva vaginae" based on phylogenetic analyses, and provide genomic and metabolomic evidence that this candidate species may metabolize D-lactate, produce trimethylamine (one of the chemicals responsible for BV-associated odor), and be motile. The cMAG and the six MAGs are valuable resources that will further contribute to our understanding of the heterogeneous etiology of bacterial vaginosis., (Copyright © 2020 Holm, France, Ma, McComb, Robinson, Mehta, Tallon, Brotman and Ravel.)

Background: Plasmodium falciparum (Pf) whole-organism sporozoite vaccines have been shown to provide significant protection against controlled human malaria infection (CHMI) in clinical trials. Initial CHMI studies showed significantly higher durable protection against homologous than heterologous strains, suggesting the presence of strain-specific vaccine-induced protection. However, interpretation of these results and understanding of their relevance to vaccine efficacy have been hampered by the lack of knowledge on genetic differences between vaccine and CHMI strains, and how these strains are related to parasites in malaria endemic regions., Methods: Whole genome sequencing using long-read (Pacific Biosciences) and short-read (Illumina) sequencing platforms was conducted to generate de novo genome assemblies for the vaccine strain, NF54, and for strains used in heterologous CHMI (7G8 from Brazil, NF166.C8 from Guinea, and NF135.C10 from Cambodia). The assemblies were used to characterize sequences in each strain relative to the reference 3D7 (a clone of NF54) genome. Strains were compared to each other and to a collection of clinical isolates (sequenced as part of this study or from public repositories) from South America, sub-Saharan Africa, and Southeast Asia., Results: While few variants were detected between 3D7 and NF54, we identified tens of thousands of variants between NF54 and the three heterologous strains. These variants include SNPs, indels, and small structural variants that fall in regulatory and immunologically important regions, including transcription factors (such as PfAP2-L and PfAP2-G) and pre-erythrocytic antigens that may be key for sporozoite vaccine-induced protection. Additionally, these variants directly contributed to diversity in immunologically important regions of the genomes as detected through in silico CD8 + T cell epitope predictions. Of all heterologous strains, NF135.C10 had the highest number of unique predicted epitope sequences when compared to NF54. Comparison to global clinical isolates revealed that these four strains are representative of their geographic origin despite long-term culture adaptation; of note, NF135.C10 is from an admixed population, and not part of recently formed subpopulations resistant to artemisinin-based therapies present in the Greater Mekong Sub-region., Conclusions: These results will assist in the interpretation of vaccine efficacy of whole-organism vaccines against homologous and heterologous CHMI.

To better understand the transcriptomic interplay of organisms associated with lymphatic filariasis, we conducted multispecies transcriptome sequencing (RNA-Seq) on the filarial nematode Brugia malayi , its Wolbachia endosymbiont w Bm, and its laboratory vector Aedes aegypti across the entire B. malayi life cycle. In w Bm, transcription of the noncoding 6S RNA suggests that it may be a regulator of bacterial cell growth, as its transcript levels correlate with bacterial replication rates. For A. aegypti , the transcriptional response reflects the stress that B. malayi infection exerts on the mosquito with indicators of increased energy demand. In B. malayi , expression modules associated with adult female samples consistently contained an overrepresentation of genes involved in chromatin remodeling, such as the bromodomain-containing proteins. All bromodomain-containing proteins encoded by B. malayi were observed to be upregulated in the adult female, embryo, and microfilaria life stages, including 2 members of the bromodomain and extraterminal (BET) protein family. The BET inhibitor JQ1(+), originally developed as a cancer therapeutic, caused lethality of adult worms in vitro , suggesting that it may be a potential therapeutic that can be repurposed for treating lymphatic filariasis. IMPORTANCE The current treatment regimen for lymphatic filariasis is mostly microfilaricidal. In an effort to identify new drug candidates for lymphatic filariasis, we conducted a three-way transcriptomics/systems biology study of one of the causative agents of lymphatic filariasis, Brugia malayi , its Wolbachia endosymbiont w Bm, and its vector host Aedes aegypti at 16 distinct B. malayi life stages. B. malayi upregulates the expression of bromodomain-containing proteins in the adult female, embryo, and microfilaria stages. In vitro , we find that the existing cancer therapeutic JQ1(+), which is a bromodomain and extraterminal protein inhibitor, has adulticidal activity in B. malayi ., (Copyright © 2019 Chung et al.)

Background: We sought to evaluate intratumor heterogeneity in squamous cell carcinoma of the oral cavity (OCC) and specifically determine the effect of physical separation and histologic differentiation within the same tumor., Methods: We performed whole exome sequencing on five biopsy sites-two from well-differentiated, two from poorly differentiated regions, and one from normal parenchyma-from five primary OCC specimens., Results: We found high levels of intratumor heterogeneity and, in four primary tumors, identified only 0 to 2 identical mutations in all subsites. We found that the heterogeneity inversely correlated with physical separation and that pairs of well-differentiated samples were more similar to each other than analogous poorly differentiated specimens. Only TP53 mutations, but not other purported "driver mutations" in head and neck squamous cell carcinoma, were found in multiple biopsy sites., Conclusion: These data highlight the challenges to characterization of the mutational landscape of OCC with single site biopsy and have implications for personalized medicine., (© 2019 Wiley Periodicals, Inc.)

Here, we present a comprehensive transcriptomics data set of Brugia malayi, its Wolbachia endosymbiont w Bm, and its vector host. This study samples from 16 stages across the entire B. malayi life cycle, including stage 1 through 4 larvae, adult males and females, embryos, immature microfilariae, and mature microfilariae.

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.

Enrichment methodologies enable the analysis of minor members in multi-species transcriptomic data. We compared the standard enrichment of bacterial and eukaryotic mRNA to a targeted enrichment using an Agilent SureSelect (AgSS) capture for Brugia malayi, Aspergillus fumigatus, and the Wolbachia endosymbiont of B. malayi (wBm). Without introducing significant systematic bias, the AgSS quantitatively enriched samples, resulting in more reads mapping to the target organism. The AgSS-enriched libraries consistently had a positive linear correlation with their unenriched counterparts (r 2  = 0.559-0.867). Up to a 2,242-fold enrichment of RNA from the target organism was obtained following a power law (r 2  = 0.90), with the greatest fold enrichment achieved in samples with the largest ratio difference between the major and minor members. While using a single total library for prokaryote and eukaryote enrichment from a single RNA sample could be beneficial for samples where RNA is limiting, we observed a decrease in reads mapping to protein coding genes and an increase in multi-mapping reads to rRNAs in AgSS enrichments from eukaryotic total RNA libraries compared to eukaryotic poly(A)-enriched libraries. Our results support a recommendation of using AgSS targeted enrichment on poly(A)-enriched libraries for eukaryotic captures, and total RNA libraries for prokaryotic captures, to increase the robustness of multi-species transcriptomic studies.

Staphylococcus aureus is a major human pathogen with well-characterized bacteriophage contributions to its virulence potential. Recently, we identified plasmidial and episomal prophages in S. aureus strains using an extra-chromosomal DNA (exDNA) isolation and sequencing approach, uncovering the plasmidial phage ϕBU01, which was found to encode important virulence determinants. Here, we expanded our extra-chromosomal sequencing of S. aureus , selecting 15 diverse clinical isolates with known chromosomal sequences for exDNA isolation and next-generation sequencing. We uncovered the presence of additional episomal prophages in 5 of 15 samples, but did not identify any plasmidial prophages. exDNA isolation was found to enrich for circular prophage elements, and qPCR characterization of the strains revealed that such prophage enrichment is detectable only in exDNA samples and would likely be missed in whole-genome DNA preparations (e.g., detection of episomal prophages did not correlate with higher prophage excision rates nor higher excised prophage copy numbers in qPCR experiments using whole-genome DNA). In S. aureus MSSA476, we found that enrichment and excision of the prophage ϕSa4ms into the cytoplasm was temporal and that episomal prophage localization did not appear to be a precursor to lytic cycle replication, suggesting ϕSa4ms excision into the cytoplasm may be part of a novel lysogenic switch. For example, we show that ϕSa4ms excision alters the promoter and transcription of htrA 2 , encoding a stress-response serine protease, and that alternative promotion of htrA 2 confers increased heat-stress survival in S. aureus COL. Overall, exDNA isolation and focused sequencing may offer a more complete genomic picture for bacterial pathogens, offering insights into important chromosomal dynamics likely missed with whole-genome DNA-based approaches.

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., (Copyright © 2017 Hooven et al.)

Bordetella holmesii causes respiratory and invasive diseases in humans, but its pathogenesis remains poorly understood. We report here the genome sequences of seven bacteremia isolates of B. holmesii , including the type strain. Comparative analysis of these sequences may aid studies of B. holmesii biology and assist in the development of species-specific diagnostic strategies., (Copyright © 2017 Tettelin et al.)

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.

As sequencing technologies have evolved, the tools to analyze these sequences have made similar advances. However, for multi-species samples, we observed important and adverse differences in alignment specificity and computation time for bwa- mem (Burrows-Wheeler aligner-maximum exact matches) relative to bwa-aln. Therefore, we sought to optimize bwa-mem for alignment of data from multi-species samples in order to reduce alignment time and increase the specificity of alignments. In the multi-species cases examined, there was one majority member (i.e. Plasmodium falciparum or Brugia malayi ) and one minority member (i.e. human or the Wolbachia endosymbiont w Bm) of the sequence data. Increasing bwa-mem seed length from the default value reduced the number of read pairs from the majority sequence member that incorrectly aligned to the reference genome of the minority sequence member. Combining both source genomes into a single reference genome increased the specificity of mapping, while also reducing the central processing unit (CPU) time. In Plasmodium , at a seed length of 18 nt, 24.1 % of reads mapped to the human genome using 1.7±0.1 CPU hours, while 83.6 % of reads mapped to the Plasmodium genome using 0.2±0.0 CPU hours (total: 107.7 % reads mapping; in 1.9±0.1 CPU hours). In contrast, 97.1 % of the reads mapped to a combined Plasmodium- human reference in only 0.7±0.0 CPU hours. Overall, the results suggest that combining all references into a single reference database and using a 23 nt seed length reduces the computational time, while maximizing specificity. Similar results were found for simulated sequence reads from a mock metagenomic data set. We found similar improvements to computation time in a publicly available human-only data set.

Neorickettsia helminthoeca, a type species of the genus Neorickettsia, is an endosymbiont of digenetic trematodes of veterinary importance. Upon ingestion of salmonid fish parasitized with infected trematodes, canids develop salmon poisoning disease (SPD), an acute febrile illness that is particularly severe and often fatal in dogs without adequate treatment. We determined and analysed the complete genome sequence of N. helminthoeca: a single small circular chromosome of 884 232 bp encoding 774 potential proteins. N. helminthoeca is unable to synthesize lipopolysaccharides and most amino acids, but is capable of synthesizing vitamins, cofactors, nucleotides and bacterioferritin. N. helminthoeca is, however, distinct from majority of the family Anaplasmataceae to which it belongs, as it encodes nearly all enzymes required for peptidoglycan biosynthesis, suggesting its structural hardiness and inflammatory potential. Using sera from dogs that were experimentally infected by feeding with parasitized fish or naturally infected in southern California, Western blot analysis revealed that among five predicted N. helminthoeca outer membrane proteins, P51 and strain-variable surface antigen were uniformly recognized. Our finding will help understanding pathogenesis, prevalence of N. helminthoeca infection among trematodes, canids and potentially other animals in nature to develop effective SPD diagnostic and preventive measures. Recent progresses in large-scale genome sequencing have been uncovering broad distribution of Neorickettsia spp., the comparative genomics will facilitate understanding of biology and the natural history of these elusive environmental bacteria., (© 2017 The Authors. Microbial Biotechnology published by John Wiley & Sons Ltd and Society for Applied Microbiology.)

We report here paired isogenic Burkholderia pseudomallei genomes obtained from three patients receiving intravenous meropenem for melioidosis treatment, with post-meropenem isolates developing decreased susceptibility. Two genomes were finished, and four were drafted to improved high-quality standard. These genomes will be used to identify meropenem resistance mechanisms in B. pseudomallei ., (Copyright © 2017 Price et al.)

Despite numerous advances in genomics and bioinformatics, technological hurdles remain to examine host-microbe transcriptomics. Sometimes the transcriptome of either or both can be ascertained merely by generating more sequencing reads. However, many cases exist where bacterial mRNA needs to be enriched further to enable cost-effective sequencing of the pathogen or endosymbiont. While a suitable method is commercially available for mammalian samples of this type, development of such methods has languished for invertebrate samples. Furthermore, a common method across multiple taxa would facilitate comparisons between bacteria in invertebrate vectors and their vertebrate hosts. Here, a method is described to concurrently remove polyadenylated transcripts, prokaryotic rRNA, and eukaryotic rRNA, including those with low amounts of starting material (e.g. 100 ng). In a Wolbachia-Drosophila system, this bacterial mRNA enrichment yielded a 3-fold increase in Wolbachia mRNA abundance and a concomitant 3.3-fold increase in the percentage of transcripts detected. More specifically, 70% of the genome could be recovered by transcriptome sequencing compared to 21% in the total RNA. Sequencing of similar bacterial mRNA-enriched samples generated from Ehrlichia-infected canine cells covers 93% of the Ehrlichia genome, suggesting ubiquitous transcription across the entire Ehrlichia chaffeensis genome. This technique can potentially be used to enrich bacterial mRNA in many studies of host-microbe interactions.

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., (Copyright © 2016, American Society for Microbiology. All Rights Reserved.)

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.

Peripheral T cell lymphoma (PTCL) is a heterogeneous malignancy with poor response to current therapeutic strategies and incompletely characterized genetics. We conducted whole exome sequencing of matched PTCL and non-malignant samples from 12 patients, spanning 8 subtypes, to identify potential oncogenic mutations in PTCL. Analysis of the mutations identified using computational algorithms, CHASM, PolyPhen2, PROVEAN, and MutationAssessor to predict the impact of these mutations on protein function and PTCL tumorigenesis, revealed 104 somatic mutations that were selected as high impact by all four algorithms. Our analysis identified recurrent somatic missense or nonsense mutations in 70 genes, 9 of which contained mutations predicted significant by all 4 algorithms: ATM, RUNX1T1, WDR17, NTRK3, TP53, TRMT12, CACNA2D1, INTS8, and KCNH8. We observed somatic mutations in ATM (ataxia telangiectasia-mutated) in 5 out of the 12 samples and mutations in the common gamma chain (γc) signaling pathway (JAK3, IL2RG, STAT5B) in 3 samples, all of which also harbored mutations in ATM. Our findings contribute insights into the genetics of PTCL and suggest a relationship between γc signaling and ATM in T cell malignancy.

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.

Macrolide resistance, emerging in Streptococcus pneumoniae and other Gram-positive bacteria, is increasingly due to efflux pumps encoded by mef/mel(msr) operons found on discrete mobile genetic elements. The regulation of mef/mel(msr) in these elements is not well understood. We identified the mef(E)/mel transcriptional start, localized the mef(E)/mel promoter, and demonstrated attenuation of transcription as a mechanism of regulation of macrolide-inducible mef-mediated macrolide resistance in S. pneumoniae. The mef(E)/mel transcriptional start site was a guanine 327 bp upstream of mef(E). Consensus pneumococcal promoter -10 (5'-TATACT-3') and -35 (5'-TTGAAC-3') boxes separated by 17 bp were identified 7 bp upstream of the start site. Analysis of the predicted secondary structure of the 327 5' region identified four pairs of inverted repeats R1-R8 predicted to fold into stem-loops, a small leader peptide [MTASMRLR, (Mef(E)L)] required for macrolide induction and a Rho-independent transcription terminator. RNA-seq analyses provided confirmation of transcriptional attenuation. In addition, expression of mef(E)L was also influenced by mef(E)L-dependent mRNA stability. The regulatory region 5' of mef(E) was highly conserved in other mef/mel(msr)-containing elements including Tn1207.1 and the 5612IQ complex in pneumococci and Tn1207.3 in Group A streptococci, indicating a regulatory mechanism common to a wide variety of Gram-positive bacteria containing mef/mel(msr) elements.

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.

Cassava is a major tropical food crop in the Euphorbiaceae family that has high carbohydrate production potential and adaptability to diverse environments. Here we present the draft genome sequences of a wild ancestor and a domesticated variety of cassava and comparative analyses with a partial inbred line. We identify 1,584 and 1,678 gene models specific to the wild and domesticated varieties, respectively, and discover high heterozygosity and millions of single-nucleotide variations. Our analyses reveal that genes involved in photosynthesis, starch accumulation and abiotic stresses have been positively selected, whereas those involved in cell wall biosynthesis and secondary metabolism, including cyanogenic glucoside formation, have been negatively selected in the cultivated varieties, reflecting the result of natural selection and domestication. Differences in microRNA genes and retrotransposon regulation could partly explain an increased carbon flux towards starch accumulation and reduced cyanogenic glucoside accumulation in domesticated cassava. These results may contribute to genetic improvement of cassava through better understanding of its biology.

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.

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.

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.

First identified in 1982 as a human pathogen, enterohemorrhagic Escherichia coli of the O157:H7 serotype is a major cause of food-borne acquired human infections. Here, we report the genome sequence of the first known strain of this serotype isolated in the United States.