Multi-platform next-generation sequencing of the domestic turkey (Meleagris gallopavo): genome assembly and analysis.

A synergistic combination of two next-generation sequencing platforms with a detailed comparative BAC physical contig map provided a cost-effective assembly of the genome sequence of the domestic turkey (Meleagris gallopavo). Heterozygosity of the sequenced source genome allowed discovery of more than 600,000 high quality single nucleotide variants. Despite this heterozygosity, the current genome assembly (∼1.1 Gb) includes 917 Mb of sequence assigned to specific turkey chromosomes. Annotation identified nearly 16,000 genes, with 15,093 recognized as protein coding and 611 as non-coding RNA genes. Comparative analysis of the turkey, chicken, and zebra finch genomes, and comparing avian to mammalian species, supports the characteristic stability of avian genomes and identifies genes unique to the avian lineage. Clear differences are seen in number and variety of genes of the avian immune system where expansions and novel genes are less frequent than examples of gene loss. The turkey genome sequence provides resources to further understand the evolution of vertebrate genomes and genetic variation underlying economically important quantitative traits in poultry. This integrated approach may be a model for providing both gene and chromosome level assemblies of other species with agricultural, ecological, and evolutionary interest.

Genomic comparison between a virulent type A1 strain of Francisella tularensis and its attenuated O-antigen mutant.

We report the complete genome sequences of TI0902, a highly virulent type A1 strain, and TIGB03, a related, attenuated chemical mutant strain. Compared to the wild type, the mutant strain had 45 point mutations and a 75.9-kb duplicated region that had not been previously observed in Francisella species.

Identification of an Attenuated Substrain of Francisella tularensis SCHU S4 by Phenotypic and Genotypic Analyses.

Pneumonic tularemia is a highly debilitating and potentially fatal disease caused by inhalation of Francisella tularensis. Most of our current understanding of its pathogenesis is based on the highly virulent F. tularensis subsp. tularensis strain SCHU S4. However, multiple sources of SCHU S4 have been maintained and propagated independently over the years, potentially generating genetic variants with altered virulence. In this study, the virulence of four SCHU S4 stocks (NR-10492, NR-28534, NR-643 from BEI Resources and FTS-635 from Battelle Memorial Institute) along with another virulent subsp. tularensis strain, MA00-2987, were assessed in parallel. In the Fischer 344 rat model of pneumonic tularemia, NR-643 and FTS-635 were found to be highly attenuated compared to NR-10492, NR-28534, and MA00-2987. In the NZW rabbit model of pneumonic tularemia, NR-643 caused morbidity but not mortality even at a dose equivalent to 500x the LD50 for NR-10492. Genetic analyses revealed that NR-10492 and NR-28534 were identical to each other, and nearly identical to the reference SCHU S4 sequence. NR-643 and FTS-635 were identical to each other but were found to have nine regions of difference in the genomic sequence when compared to the published reference SCHU S4 sequence. Given the genetic differences and decreased virulence, NR-643/FTS-635 should be clearly designated as a separate SCHU S4 substrain and no longer utilized in efficacy studies to evaluate potential vaccines and therapeutics against tularemia.

Characterization of the hypothalamic transcriptome in response to food deprivation reveals global changes in long noncoding RNA, and cell cycle response genes.

The hypothalamus integrates energy balance information from the periphery using different neuronal subtypes within each of the hypothalamic areas. However, the effects of prandial state on global mRNA, microRNA and long noncoding (lnc) RNA expression within the whole hypothalamus are largely unknown. In this study, mice were given either a 24-h fast, or ad libitum access to food. RNA samples were analyzed by microarray, and then a subset was confirmed using quantitative real-time PCR (QPCR). A total of 540 mRNAs were either up- or down-regulated with food deprivation. Since gene ontology enrichment analyses identified several categories of mRNAs related to cell cycle processes, ten cell-cycle-related genes were further analyzed using QPCR with six confirmed to be significantly up-regulated and one down-regulated in response to 24-h fasting. While 22 independent microRNAs were differentially expressed by microarray, secondary analysis by QPCR failed to confirm significant changes with fasting. There were 622 lncRNAs identified as differentially expressed, and of three tested by QPCR, two were confirmed. Overall, this is the first time that expression of hypothalamic lncRNAs has been shown to be responsive to food deprivation. In addition, this study is the first to identify a list of lncRNAs with high expression in RNA extracted from hypothalamus. Individual contributions from specific miRNA, lncRNA and mRNAs to the food deprivation response can now be further studied at the physiological and biochemical levels.

OpaR controls a network of downstream transcription factors in Vibrio parahaemolyticus BB22OP.

Vibrio parahaemolyticus is an emerging world-wide human pathogen that is associated with food-borne gastroenteritis when raw or undercooked seafood is consumed. Expression of virulence factors in this organism is modulated by the phenomenon known as quorum sensing, which permits differential gene regulation at low versus high cell density. The master regulator of quorum sensing in V. parahaemolyticus is OpaR. OpaR not only controls virulence factor gene expression, but also the colony and cellular morphology associated with growth on a surface and biofilm formation. Whole transcriptome Next Generation sequencing (RNA-Seq) was utilized to determine the OpaR regulon by comparing strains BB22OP (opaR+, LM5312) and BB22TR (∆opaR1, LM5674). This work, using the published V. parahaemolyticus BB22OP genome sequence, confirms and expands upon a previous microarray analysis for these two strains that used an Affymetrix GeneChip designed from the closely related V. parahaemolyticus RIMD2210633 genome sequence. Overall there was excellent correlation between the microarray and RNA-Seq data. Eleven transcription factors under OpaR control were identified by both methods and further confirmed by quantitative reverse transcription PCR (qRT-PCR) analysis. Nine of these transcription factors were demonstrated to be direct OpaR targets via in vitro electrophoretic mobility shift assays with purified hexahistidine-tagged OpaR. Identification of the direct and indirect targets of OpaR, including small RNAs, will enable the construction of a network map of regulatory interactions important for the switch between the nonpathogenic and pathogenic states.

Complete genome sequence of Salmonella enterica serovar Typhimurium VNP20009, a strain engineered for tumor targeting.

A mutagenized and genetically modified derivative of Salmonella enterica serovar Typhimurium 14028S, VNP20009 (ATCC 202165), is attenuated in host virulence and accumulates preferentially in tumors. Here, we report the complete genome of this anticancer therapy agent consisting of one chromosome and one virulence plasmid. The major genetic features that distinguish VNP20009 from its parental strain are: an engineered msbB deletion, a 3'-extension in pykA, a Tn10-caused 16.6-kbp inversion leading to the disruption of purM, a 108-kbp Suwwan deletion resulting in the loss of 128 genes, and 50 non-synonymous SNPs.

Complete Genome Sequence of Prepandemic Vibrio parahaemolyticus BB22OP.

The number of inflammatory gastroenteritis outbreaks due to the food-borne pathogen Vibrio parahaemolyticus is rising sharply worldwide and in the United States in particular. Here we report the complete, annotated genome sequence of the prepandemic V. parahaemolyticus strain BB22OP and make some initial comparisons to the complete genome sequence for pandemic strain RIMD2210633.

Analysis of high-throughput sequencing data.

Next-generation sequencing has revolutionized biology by exponentially increasing sequencing output while dramatically lowering costs. High-throughput sequence data with shorter reads has opened up new applications such as whole genome resequencing, indel and SNP detection, transcriptome sequencing, etc. Several tools are available for the analysis of high-throughput sequencing data. In this chapter, we describe the use of an ultrafast alignment program, Bowtie, to align short-read sequence (SRS) data against the Arabidopsis reference genome. The alignment files generated from Bowtie will be used to identify SNPs and indels using Maq.