Search

Abstract Background Numerous methods exist for enriching bacterial or mammalian mRNA prior to transcriptome experiments. Yet there persists a need for methods to enrich for mRNA in non-mammalian animal systems. For example, insects contain many important and interesting obligate intracellular bacteria, including endosymbionts and vector-borne pathogens. Such obligate intracellular bacteria are difficult to study by traditional methods. Therefore, genomics has greatly increased our understanding of these bacteria. Efficient subtraction methods are needed for removing both bacteria and insect rRNA in these systems to enable transcriptome-based studies. Findings A method is described that efficiently removes >95% of insect rRNA from total RNA samples, as determined by microfluidics and transcriptome sequencing. This subtraction yielded a 6.2-fold increase in mRNA abundance. Such a host rRNA-depletion strategy, in combination with bacterial rRNA depletion, is necessary to analyze transcription of obligate intracellular bacteria. Here, transcripts were identified that arise from a lateral gene transfer of an entire Wolbachia bacterial genome into a Drosophila ananassae chromosome. In this case, an rRNA depletion strategy is preferred over polyA-based enrichment since transcripts arising from bacteria-to-animal lateral gene transfer may not be poly-adenylated. Conclusions This enrichment method yields a significant increase in mRNA abundance when poly-A selection is not suitable. It can be used in combination with bacterial rRNA subtraction to enable experiments to simultaneously measure bacteria and insect mRNA in vector and endosymbiont biology experiments.

Abstract Background Tetrahymena thermophila, a widely studied model for cellular and molecular biology, is a binucleated single-celled organism with a germline micronucleus (MIC) and somatic macronucleus (MAC). The recent draft MAC genome assembly revealed low sequence repetitiveness, a result of the epigenetic removal of invasive DNA elements found only in the MIC genome. Such low repetitiveness makes complete closure of the MAC genome a feasible goal, which to achieve would require standard closure methods as well as removal of minor MIC contamination of the MAC genome assembly. Highly accurate preliminary annotation of Tetrahymena's coding potential was hindered by the lack of both comparative genomic sequence information from close relatives and significant amounts of cDNA evidence, thus limiting the value of the genomic information and also leaving unanswered certain questions, such as the frequency of alternative splicing. Results We addressed the problem of MIC contamination using comparative genomic hybridization with purified MIC and MAC DNA probes against a whole genome oligonucleotide microarray, allowing the identification of 763 genome scaffolds likely to contain MIC-limited DNA sequences. We also employed standard genome closure methods to essentially finish over 60% of the MAC genome. For the improvement of annotation, we have sequenced and analyzed over 60,000 verified EST reads from a variety of cellular growth and development conditions. Using this EST evidence, a combination of automated and manual reannotation efforts led to updates that affect 16% of the current protein-coding gene models. By comparing EST abundance, many genes showing apparent differential expression between these conditions were identified. Rare instances of alternative splicing and uses of the non-standard amino acid selenocysteine were also identified. Conclusion We report here significant progress in genome closure and reannotation of Tetrahymena thermophila. Our experience to date suggests that complete closure of the MAC genome is attainable. Using the new EST evidence, automated and manual curation has resulted in substantial improvements to the over 24,000 gene models, which will be valuable to researchers studying this model organism as well as for comparative genomics purposes.

Abstract Background Bracoviruses (BVs), a group of double-stranded DNA viruses with segmented genomes, are mutualistic endosymbionts of parasitoid wasps. Virus particles are replication deficient and are produced only by female wasps from proviral sequences integrated into the wasp genome. Virus particles are injected along with eggs into caterpillar hosts, where viral gene expression facilitates parasitoid survival and therefore perpetuation of proviral DNA. Here we describe a 223 kbp region of Glyptapanteles indiensis genomic DNA which contains a part of the G. indiensis bracovirus (GiBV) proviral genome. Results Eighteen of ~24 GiBV viral segment sequences are encoded by 7 non-overlapping sets of BAC clones, revealing that some proviral segment sequences are separated by long stretches of intervening DNA. Two overlapping BACs, which contain a locus of 8 tandemly arrayed proviral segments flanked on either side by ~35 kbp of non-packaged DNA, were sequenced and annotated. Structural and compositional analyses of this cluster revealed it exhibits a G+C and nucleotide composition distinct from the flanking DNA. By analyzing sequence polymorphisms in the 8 GiBV viral segment sequences, we found evidence for widespread selection acting on both protein-coding and non-coding DNA. Comparative analysis of viral and proviral segment sequences revealed a sequence motif involved in the excision of proviral genome segments which is highly conserved in two other bracoviruses. Conclusion Contrary to current concepts of bracovirus proviral genome organization our results demonstrate that some but not all GiBV proviral segment sequences exist in a tandem array. Unexpectedly, non-coding DNA in the 8 proviral genome segments which typically occupies ~70% of BV viral genomes is under selection pressure suggesting it serves some function(s). We hypothesize that selection acting on GiBV proviral sequences maintains the genetic island-like nature of the cluster of proviral genome segments described herein. In contrast to large differences in the predicted gene composition of BV genomes, sequences that appear to mediate processes of viral segment formation, such as proviral segment excision and circularization, appear to be highly conserved, supporting the hypothesis of a single origin for BVs.

Abstract Background Carrot (Daucus carota) is a major food crop in the US and worldwide. Its capacity for storage and its lifecycle as a biennial make it an attractive species for the introduction of foreign genes, especially for oral delivery of vaccines and other therapeutic proteins. Until recently efforts to express recombinant proteins in carrot have had limited success in terms of protein accumulation in the edible tap roots. Plastid genetic engineering offers the potential to overcome this limitation, as demonstrated by the accumulation of BADH in chromoplasts of carrot taproots to confer exceedingly high levels of salt resistance. The complete plastid genome of carrot provides essential information required for genetic engineering. Additionally, the sequence data add to the rapidly growing database of plastid genomes for assessing phylogenetic relationships among angiosperms. Results The complete carrot plastid genome is 155,911 bp in length, with 115 unique genes and 21 duplicated genes within the IR. There are four ribosomal RNAs, 30 distinct tRNA genes and 18 intron-containing genes. Repeat analysis reveals 12 direct and 2 inverted repeats ≥ 30 bp with a sequence identity ≥ 90%. Phylogenetic analysis of nucleotide sequences for 61 protein-coding genes using both maximum parsimony (MP) and maximum likelihood (ML) were performed for 29 angiosperms. Phylogenies from both methods provide strong support for the monophyly of several major angiosperm clades, including monocots, eudicots, rosids, asterids, eurosids II, euasterids I, and euasterids II. Conclusion The carrot plastid genome contains a number of dispersed direct and inverted repeats scattered throughout coding and non-coding regions. This is the first sequenced plastid genome of the family Apiaceae and only the second published genome sequence of the species-rich euasterid II clade. Both MP and ML trees provide very strong support (100% bootstrap) for the sister relationship of Daucus with Panax in the euasterid II clade. These results provide the best taxon sampling of complete chloroplast genomes and the strongest support yet for the sister relationship of Caryophyllales to the asterids. The availability of the complete plastid genome sequence should facilitate improved transformation efficiency and foreign gene expression in carrot through utilization of endogenous flanking sequences and regulatory elements.

Abstract Background Cotton (Gossypium hirsutum) is the most important fiber crop grown in 90 countries. In 2004–2005, US farmers planted 79% of the 5.7-million hectares of nuclear transgenic cotton. Unfortunately, genetically modified cotton has the potential to hybridize with other cultivated and wild relatives, resulting in geographical restrictions to cultivation. However, chloroplast genetic engineering offers the possibility of containment because of maternal inheritance of transgenes. The complete chloroplast genome of cotton provides essential information required for genetic engineering. In addition, the sequence data were used to assess phylogenetic relationships among the major clades of rosids using cotton and 25 other completely sequenced angiosperm chloroplast genomes. Results The complete cotton chloroplast genome is 160,301 bp in length, with 112 unique genes and 19 duplicated genes within the IR, containing a total of 131 genes. There are four ribosomal RNAs, 30 distinct tRNA genes and 17 intron-containing genes. The gene order in cotton is identical to that of tobacco but lacks rpl22 and infA. There are 30 direct and 24 inverted repeats 30 bp or longer with a sequence identity ≥ 90%. Most of the direct repeats are within intergenic spacer regions, introns and a 72 bp-long direct repeat is within the psaA and psaB genes. Comparison of protein coding sequences with expressed sequence tags (ESTs) revealed nucleotide substitutions resulting in amino acid changes in ndhC, rpl23, rpl20, rps3 and clpP. Phylogenetic analysis of a data set including 61 protein-coding genes using both maximum likelihood and maximum parsimony were performed for 28 taxa, including cotton and five other angiosperm chloroplast genomes that were not included in any previous phylogenies. Conclusion Cotton chloroplast genome lacks rpl22 and infA and contains a number of dispersed direct and inverted repeats. RNA editing resulted in amino acid changes with significant impact on their hydropathy. Phylogenetic analysis provides strong support for the position of cotton in the Malvales in the eurosids II clade sister to Arabidopsis in the Brassicales. Furthermore, there is strong support for the placement of the Myrtales sister to the eurosid I clade, although expanded taxon sampling is needed to further test this relationship.

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., Competing Interests: Conflict of interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

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., Competing Interests: The author(s) declare(s) that there is no conflict of interest., (Copyright: © 2024 Meiller T, et al.)

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., (© 2023. The Author(s).)

Theileria parva is an economically important, intracellular, tick-transmitted parasite of cattle. A live vaccine against the parasite is effective against challenge from cattle-transmissible T. parva but not against genotypes originating from the African Cape buffalo, a major wildlife reservoir, prompting the need to characterize genome-wide variation within and between cattle- and buffalo-associated T. parva populations. Here, we describe a capture-based target enrichment approach that enables, for the first time, de novo assembly of nearly complete T. parva genomes derived from infected host cell lines. This approach has exceptionally high specificity and sensitivity and is successful for both cattle- and buffalo-derived T. parva parasites. De novo genome assemblies generated for cattle genotypes differ from the reference by ~54K single nucleotide polymorphisms (SNPs) throughout the 8.31 Mb genome, an average of 6.5 SNPs/kb. We report the first buffalo-derived T. parva genome, which is ~20 kb larger than the genome from the reference, cattle-derived, Muguga strain, and contains 25 new potential genes. The average non-synonymous nucleotide diversity (πN) per gene, between buffalo-derived T. parva and the Muguga strain, was 1.3%. This remarkably high level of genetic divergence is supported by an average Wright's fixation index (FST), genome-wide, of 0.44, reflecting a degree of genetic differentiation between cattle- and buffalo-derived T. parva parasites more commonly seen between, rather than within, species. These findings present clear implications for vaccine development, further demonstrated by the ability to assemble nearly all known antigens in the buffalo-derived strain, which will be critical in design of next generation vaccines. The DNA capture approach used provides a clear advantage in specificity over alternative T. parva DNA enrichment methods used previously, such as those that utilize schizont purification, is less labor intensive, and enables in-depth comparative genomics in this apicomplexan parasite., Competing Interests: The authors have declared that no competing interests exist.

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.

Surface/ground waters could be polluted when rain-water and/or snow-melt water infiltrate through waste rock piles at mine sites and dissolve secondary minerals (salts) from rock surfaces. It is important to reduce solute loading by the optimal configuration of waste rock piles. This requires the proper definition and determination of the characteristic properties of waste rock piles in terms of metal leaching and, in particular, rate control mechanisms and scaling laws, and their dependence upon configuration variables. For revealing these characteristic properties this paper proposes a pile-scale C-Q relation: C = C s (1 - e -P/Q ), (P ≡ kλβψ), where C and C s are respectively solute concentration and particle's saturation concentration, Q is the flow rate of the water through a waste rock pile, k represents the effective or average dissolution coefficient of a mineral specie from rock surfaces, β represents rock pile depth, λ represents the ratio of the sum of the surface areas of rocks to the volume that the rocks occupy, and ψ is the sum of the cross-sections of water-flow channels in a waste rock pile. The two characteristic properties revealed by the C-Q relation are: (1) P, the product of k, λ, β, and ψ (P ≡ kλβψ), which is the characteristic property of a waste rock pile in terms of metal leaching, named here the solute production potential; and (2) the ratio of P to Q, P/Q, a non-dimensional number, designated as α (α ≡ P/Q), named here the rate control quotient, which is the scaling law and the rate control mechanism indicator. The value of α quantitatively indicates what controls the rate of mineral dissolution, and it also relates smaller-scale metal-leaching testing results to their corresponding full scales. When α becomes small, say α < 0.5, the rate of solute production potential P becomes in control, and the solute loading is nearly independent of Q; when α becomes larger, say α > 2.5, solute concentration would become close to its saturation concentration C s , and Q determines solute loading (that is, the solute loading is proportional to Q). When 0.5 < α < 2.5, both Q and P are in control, a mixed control mechanism. The 20 years of measurements of mine drainage chemistry from the main waste rock piles at the Equity Silver mine, BC, Canada, are used to illustrate how to determine the two characteristic properties P and α, and how well they are able to describe the waste rock piles in terms of metal leaching., (Crown Copyright © 2019. Published by Elsevier B.V. All rights reserved.)

FDA proactively invests in tools to support innovation of emerging technologies, such as infectious disease next generation sequencing (ID-NGS). Here, we introduce FDA-ARGOS quality-controlled reference genomes as a public database for diagnostic purposes and demonstrate its utility on the example of two use cases. We provide quality control metrics for the FDA-ARGOS genomic database resource and outline the need for genome quality gap filling in the public domain. In the first use case, we show more accurate microbial identification of Enterococcus avium from metagenomic samples with FDA-ARGOS reference genomes compared to non-curated GenBank genomes. In the second use case, we demonstrate the utility of FDA-ARGOS reference genomes for Ebola virus target sequence comparison as part of a composite validation strategy for ID-NGS diagnostic tests. The use of FDA-ARGOS as an in silico target sequence comparator tool combined with representative clinical testing could reduce the burden for completing ID-NGS clinical trials.

Introduction: Available literature has repeatedly ascertained the inverse relationship of physical activity to mental health. Engaging in regular physical activity decreases symptoms of anxiety, depression, and stress in a dose-responsive manner., Methods: Cross-sectional, national data from the 2015 Behavioral Risk Factor Surveillance System (BRFSS) were analyzed to establish the relationship between physical activity and days of poor mental health. The association between physical activity and days of poor mental health was assessed by multiple linear regression using SAS., Results: Findings support the assumption of an inverse dose-response relationship between physical activity levels and self-reported days of poor mental health. Individuals who engage in higher levels of physical activity reported significantly fewer days of poor mental health when compared to those who reported low or no physical activity (e.g. -2.23, 95%CI (-2.31, -2.15)) for highly active vs. inactive individuals (unadjusted model)). Furthermore, individuals who only met the guidelines for aerobic exercise reported significantly fewer days of poor mental health than those who only met the muscle strengthening guidelines (-0.85 (95% CI: -0.98, -0.73))., Limitations: Due to the cross-sectional study design, the direction of the relationship is unknown and reverse causation is possible., Conclusion: Higher physical activity levels were associated with fewer self-reported days of poor mental health, which reiterates the important association between physical activity and mental health with national data. The findings of this study could prove to be valuable for developing interventions aimed at improving mental health. Nonetheless, several research gaps persist that warrant further empirical inquiry., (Copyright © 2019. Published by Elsevier B.V.)

Babesia microti, a tick-transmitted, intraerythrocytic protozoan parasite circulating mainly among small mammals, is the primary cause of human babesiosis. While most cases are transmitted by Ixodes ticks, the disease may also be transmitted through blood transfusion and perinatally. A comprehensive analysis of genome composition, genetic diversity, and gene expression profiling of seven B. microti isolates revealed that genetic variation in isolates from the Northeast United States is almost exclusively associated with genes encoding the surface proteome and secretome of the parasite. Furthermore, we found that polymorphism is restricted to a small number of genes, which are highly expressed during infection. In order to identify pathogen-encoded factors involved in host-parasite interactions, we screened a proteome array comprised of 174 B. microti proteins, including several predicted members of the parasite secretome. Using this immuno-proteomic approach we identified several novel antigens that trigger strong host immune responses during the onset of infection. The genomic and immunological data presented herein provide the first insights into the determinants of B. microti interaction with its mammalian hosts and their relevance for understanding the selective pressures acting on parasite evolution.

It has been suggested that urinary PCA3 and TMPRSS2:ERG fusion tests and serum PHI correlate to cancer aggressiveness-related pathological criteria at prostatectomy. To evaluate and compare their ability in predicting prostate cancer aggressiveness, PHI and urinary PCA3 and TMPRSS2:ERG (T2) scores were assessed in 154 patients who underwent radical prostatectomy for biopsy-proven prostate cancer. Univariate and multivariate analyses using logistic regression and decision curve analyses were performed. All three markers were predictors of a tumor volume≥0.5 mL. Only PHI predicted Gleason score≥7. T2 score and PHI were both independent predictors of extracapsular extension(≥pT3), while multifocality was only predicted by PCA3 score. Moreover, when compared to a base model (age, digital rectal examination, serum PSA, and Gleason sum at biopsy), the addition of both PCA3 score and PHI to the base model induced a significant increase (+12%) when predicting tumor volume>0.5 mL. PHI and urinary PCA3 and T2 scores can be considered as complementary predictors of cancer aggressiveness at prostatectomy.

High throughput sequencing has accelerated the determination of genome sequences for thousands of human infectious disease pathogens and dozens of their vectors. The scale and scope of these data are enabling genotype-phenotype association studies to identify genetic determinants of pathogen virulence and drug/insecticide resistance, and phylogenetic studies to track the origin and spread of disease outbreaks. To maximize the utility of genomic sequences for these purposes, it is essential that metadata about the pathogen/vector isolate characteristics be collected and made available in organized, clear, and consistent formats. Here we report the development of the GSCID/BRC Project and Sample Application Standard, developed by representatives of the Genome Sequencing Centers for Infectious Diseases (GSCIDs), the Bioinformatics Resource Centers (BRCs) for Infectious Diseases, and the U.S. National Institute of Allergy and Infectious Diseases (NIAID), part of the National Institutes of Health (NIH), informed by interactions with numerous collaborating scientists. It includes mapping to terms from other data standards initiatives, including the Genomic Standards Consortium's minimal information (MIxS) and NCBI's BioSample/BioProjects checklists and the Ontology for Biomedical Investigations (OBI). The standard includes data fields about characteristics of the organism or environmental source of the specimen, spatial-temporal information about the specimen isolation event, phenotypic characteristics of the pathogen/vector isolated, and project leadership and support. By modeling metadata fields into an ontology-based semantic framework and reusing existing ontologies and minimum information checklists, the application standard can be extended to support additional project-specific data fields and integrated with other data represented with comparable standards. The use of this metadata standard by all ongoing and future GSCID sequencing projects will provide a consistent representation of these data in the BRC resources and other repositories that leverage these data, allowing investigators to identify relevant genomic sequences and perform comparative genomics analyses that are both statistically meaningful and biologically relevant.

First identified in 1982, Escherichia coli O157:H7 is the dominant enterohemorrhagic serotype underlying food-borne human infections in North America. Here, we report the genomes of twenty-six strains derived from patients and the bovine reservoir. These resources enable detailed whole-genome comparisons and permit investigations of genotypic and phenotypic plasticity.

The freshwater snail Biomphalaria glabrata is closely associated with the transmission of human schistosomiasis. An ecologically sound method has been proposed to control schistosomiasis using genetically modified snails to displace endemic, susceptible ones. To assess the viability of this form of biological control, studies towards understanding the molecular makeup of the snail relative to the presence of endogenous mobile genetic elements are being undertaken since they can be exploited for genetic transformation studies. We previously cloned a 1.95kb BamHI fragment in B. glabrata (BGR2) with sequence similarity to the human long interspersed nuclear element (LINE or L1). A contiguous, full-length sequence corresponding to BGR2, hereafter-named nimbus (BgI), has been identified from a B. glabrata bacterial artificial chromosome (BAC) library. Sequence analysis of the 65,764bp BAC insert contained one full-length, complete nimbus (BgI) element (element I), two full-length elements (elements II and III) containing deletions and flanked by target site duplications and 10 truncated copies. The intact nimbus (BgI) contained two open-reading frames (ORFs 1 and 2) encoding the characteristic hallmark domains found in non-long terminal repeat retrotransposons belonging to the I-clade; a nucleic acid binding protein in ORF1 and an apurinic/apyrimidinic endonuclease, reverse transcriptase and RNase H in ORF2. Phylogenetic analysis revealed that nimbus (BgI) is closely related to Drosophila (I factor), mosquito Aedes aegypti (MosquI) and chordate ascidian Ciona intestinalis (CiI) retrotransposons. Nimbus (BgI) represents the first complete mobile element characterised from a mollusk that appears to be transcriptionally active and is widely distributed in snails of the neotropics and the Old World.

Parasitic nematodes that cause elephantiasis and river blindness threaten hundreds of millions of people in the developing world. We have sequenced the approximately 90 megabase (Mb) genome of the human filarial parasite Brugia malayi and predict approximately 11,500 protein coding genes in 71 Mb of robustly assembled sequence. Comparative analysis with the free-living, model nematode Caenorhabditis elegans revealed that, despite these genes having maintained little conservation of local synteny during approximately 350 million years of evolution, they largely remain in linkage on chromosomal units. More than 100 conserved operons were identified. Analysis of the predicted proteome provides evidence for adaptations of B. malayi to niches in its human and vector hosts and insights into the molecular basis of a mutualistic relationship with its Wolbachia endosymbiont. These findings offer a foundation for rational drug design.

The complete genome of Aeromonas hydrophila ATCC 7966(T) was sequenced. Aeromonas, a ubiquitous waterborne bacterium, has been placed by the Environmental Protection Agency on the Contaminant Candidate List because of its potential to cause human disease. The 4.7-Mb genome of this emerging pathogen shows a physiologically adroit organism with broad metabolic capabilities and considerable virulence potential. A large array of virulence genes, including some identified in clinical isolates of Aeromonas spp. or Vibrio spp., may confer upon this organism the ability to infect a wide range of hosts. However, two recognized virulence markers, a type III secretion system and a lateral flagellum, that are reported in other A. hydrophila strains are not identified in the sequenced isolate, ATCC 7966(T). Given the ubiquity and free-living lifestyle of this organism, there is relatively little evidence of fluidity in terms of mobile elements in the genome of this particular strain. Notable aspects of the metabolic repertoire of A. hydrophila include dissimilatory sulfate reduction and resistance mechanisms (such as thiopurine reductase, arsenate reductase, and phosphonate degradation enzymes) against toxic compounds encountered in polluted waters. These enzymes may have bioremediative as well as industrial potential. Thus, the A. hydrophila genome sequence provides valuable insights into its ability to flourish in both aquatic and host environments.

Since its introduction a decade ago, whole-genome shotgun sequencing (WGS) has been the main approach for producing cost-effective and high-quality genome sequence data. Until now, the Sanger sequencing technology that has served as a platform for WGS has not been truly challenged by emerging technologies. The recent introduction of the pyrosequencing-based 454 sequencing platform (454 Life Sciences, Branford, CT) offers a very promising sequencing technology alternative for incorporation in WGS. In this study, we evaluated the utility and cost-effectiveness of a hybrid sequencing approach using 3730xl Sanger data and 454 data to generate higher-quality lower-cost assemblies of microbial genomes compared to current Sanger sequencing strategies alone.

The genetic relatedness of the Bacillus anthracis typing phages Gamma and Cherry was determined by nucleotide sequencing and comparative analysis. The genomes of these two phages were identical except at three variable loci, which showed heterogeneity within individual lysates and among Cherry, Wbeta, Fah, and four Gamma bacteriophage sequences.

Rice (Oryza sativa L.) chromosome 3 is evolutionarily conserved across the cultivated cereals and shares large blocks of synteny with maize and sorghum, which diverged from rice more than 50 million years ago. To begin to completely understand this chromosome, we sequenced, finished, and annotated 36.1 Mb ( approximately 97%) from O. sativa subsp. japonica cv Nipponbare. Annotation features of the chromosome include 5915 genes, of which 913 are related to transposable elements. A putative function could be assigned to 3064 genes, with another 757 genes annotated as expressed, leaving 2094 that encode hypothetical proteins. Similarity searches against the proteome of Arabidopsis thaliana revealed putative homologs for 67% of the chromosome 3 proteins. Further searches of a nonredundant amino acid database, the Pfam domain database, plant Expressed Sequence Tags, and genomic assemblies from sorghum and maize revealed only 853 nontransposable element related proteins from chromosome 3 that lacked similarity to other known sequences. Interestingly, 426 of these have a paralog within the rice genome. A comparative physical map of the wild progenitor species, Oryza nivara, with japonica chromosome 3 revealed a high degree of sequence identity and synteny between these two species, which diverged approximately 10,000 years ago. Although no major rearrangements were detected, the deduced size of the O. nivara chromosome 3 was 21% smaller than that of japonica. Synteny between rice and other cereals using an integrated maize physical map and wheat genetic map was strikingly high, further supporting the use of rice and, in particular, chromosome 3, as a model for comparative studies among the cereals.

African trypanosomes cause human sleeping sickness and livestock trypanosomiasis in sub-Saharan Africa. We present the sequence and analysis of the 11 megabase-sized chromosomes of Trypanosoma brucei. The 26-megabase genome contains 9068 predicted genes, including approximately 900 pseudogenes and approximately 1700 T. brucei-specific genes. Large subtelomeric arrays contain an archive of 806 variant surface glycoprotein (VSG) genes used by the parasite to evade the mammalian immune system. Most VSG genes are pseudogenes, which may be used to generate expressed mosaic genes by ectopic recombination. Comparisons of the cytoskeleton and endocytic trafficking systems with those of humans and other eukaryotic organisms reveal major differences. A comparison of metabolic pathways encoded by the genomes of T. brucei, T. cruzi, and Leishmania major reveals the least overall metabolic capability in T. brucei and the greatest in L. major. Horizontal transfer of genes of bacterial origin has contributed to some of the metabolic differences in these parasites, and a number of novel potential drug targets have been identified.

The complete 2,343,479-bp genome sequence of the gram-negative, pathogenic oral bacterium Porphyromonas gingivalis strain W83, a major contributor to periodontal disease, was determined. Whole-genome comparative analysis with other available complete genome sequences confirms the close relationship between the Cytophaga-Flavobacteria-Bacteroides (CFB) phylum and the green-sulfur bacteria. Within the CFB phyla, the genomes most similar to that of P. gingivalis are those of Bacteroides thetaiotaomicron and B. fragilis. Outside of the CFB phyla the most similar genome to P. gingivalis is that of Chlorobium tepidum, supporting the previous phylogenetic studies that indicated that the Chlorobia and CFB phyla are related, albeit distantly. Genome analysis of strain W83 reveals a range of pathways and virulence determinants that relate to the novel biology of this oral pathogen. Among these determinants are at least six putative hemagglutinin-like genes and 36 previously unidentified peptidases. Genome analysis also reveals that P. gingivalis can metabolize a range of amino acids and generate a number of metabolic end products that are toxic to the human host or human gingival tissue and contribute to the development of periodontal disease.

The mosquito-borne malaria parasite Plasmodium falciparum kills an estimated 0.7-2.7 million people every year, primarily children in sub-Saharan Africa. Without effective interventions, a variety of factors-including the spread of parasites resistant to antimalarial drugs and the increasing insecticide resistance of mosquitoes-may cause the number of malaria cases to double over the next two decades. To stimulate basic research and facilitate the development of new drugs and vaccines, the genome of Plasmodium falciparum clone 3D7 has been sequenced using a chromosome-by-chromosome shotgun strategy. We report here the nucleotide sequences of chromosomes 10, 11 and 14, and a re-analysis of the chromosome 2 sequence. These chromosomes represent about 35% of the 23-megabase P. falciparum genome.

The genome of the flowering plant Arabidopsis thaliana has five chromosomes. Here we report the sequence of the largest, chromosome 1, in two contigs of around 14.2 and 14.6 megabases. The contigs extend from the telomeres to the centromeric borders, regions rich in transposons, retrotransposons and repetitive elements such as the 180-base-pair repeat. The chromosome represents 25% of the genome and contains about 6,850 open reading frames, 236 transfer RNAs (tRNAs) and 12 small nuclear RNAs. There are two clusters of tRNA genes at different places on the chromosome. One consists of 27 tRNA(Pro) genes and the other contains 27 tandem repeats of tRNA(Tyr)-tRNA(Tyr)-tRNA(Ser) genes. Chromosome 1 contains about 300 gene families with clustered duplications. There are also many repeat elements, representing 8% of the sequence.

Arabidopsis thaliana is an important model system for plant biologists. In 1996 an international collaboration (the Arabidopsis Genome Initiative) was formed to sequence the whole genome of Arabidopsis and in 1999 the sequence of the first two chromosomes was reported. The sequence of the last three chromosomes and an analysis of the whole genome are reported in this issue. Here we present the sequence of chromosome 3, organized into four sequence segments (contigs). The two largest (13.5 and 9.2 Mb) correspond to the top (long) and the bottom (short) arms of chromosome 3, and the two small contigs are located in the genetically defined centromere. This chromosome encodes 5,220 of the roughly 25,500 predicted protein-coding genes in the genome. About 20% of the predicted proteins have significant homology to proteins in eukaryotic genomes for which the complete sequence is available, pointing to important conserved cellular functions among eukaryotes.

Arabidopsis thaliana (Arabidopsis) is unique among plant model organisms in having a small genome (130-140 Mb), excellent physical and genetic maps, and little repetitive DNA. Here we report the sequence of chromosome 2 from the Columbia ecotype in two gap-free assemblies (contigs) of 3.6 and 16 megabases (Mb). The latter represents the longest published stretch of uninterrupted DNA sequence assembled from any organism to date. Chromosome 2 represents 15% of the genome and encodes 4,037 genes, 49% of which have no predicted function. Roughly 250 tandem gene duplications were found in addition to large-scale duplications of about 0.5 and 4.5 Mb between chromosomes 2 and 1 and between chromosomes 2 and 4, respectively. Sequencing of nearly 2 Mb within the genetically defined centromere revealed a low density of recognizable genes, and a high density and diverse range of vestigial and presumably inactive mobile elements. More unexpected is what appears to be a recent insertion of a continuous stretch of 75% of the mitochondrial genome into chromosome 2.