Search

We present a draft sequence of the genome of Aedes aegypti , the primary vector for yellow fever and dengue fever, which at ∼1376 million base pairs is about 5 times the size of the genome of the malaria vector Anopheles gambiae . Nearly 50% of the Ae. aegypti genome consists of transposable elements. These contribute to a factor of ∼4 to 6 increase in average gene length and in sizes of intergenic regions relative to An. gambiae and Drosophila melanogaster . Nonetheless, chromosomal synteny is generally maintained among all three insects, although conservation of orthologous gene order is higher (by a factor of ∼2) between the mosquito species than between either of them and the fruit fly. An increase in genes encoding odorant binding, cytochrome P450, and cuticle domains relative to An. gambiae suggests that members of these protein families underpin some of the biological differences between the two mosquito species.

How genomic diversity within bacterial populations originates and is maintained in the presence of frequent recombination is a central problem in understanding bacterial evolution. Natural populations of Borrelia burgdorferi, the bacterial agent of Lyme disease, consist of diverse genomic groups co-infecting single individual vertebrate hosts and tick vectors. To understand mechanisms of sympatric genome differentiation in B. burgdorferi, we sequenced and compared 23 genomes representing major genomic groups in North America and Europe. Linkage analysis of >13,500 single-nucleotide polymorphisms revealed pervasive horizontal DNA exchanges. Although three times more frequent than point mutation, recombination is localized and weakly affects genome-wide linkage disequilibrium. We show by computer simulations that, while enhancing population fitness, recombination constrains neutral and adaptive divergence among sympatric genomes through periodic selective sweeps. In contrast, simulations of frequency-dependent selection with recombination produced the observed pattern of a large number of sympatric genomic groups associated with major sequence variations at the selected locus. We conclude that negative frequency-dependent selection targeting a small number of surface-antigen loci (ospC in particular) sufficiently explains the maintenance of sympatric genome diversity in B. burgdorferi without adaptive divergence. We suggest that pervasive recombination makes it less likely for local B. burgdorferi genomic groups to achieve host specialization. B. burgdorferi genomic groups in the northeastern United States are thus best viewed as constituting a single bacterial species, whose generalist nature is a key to its rapid spread and human virulence.

Castor bean (Ricinus communis) is an oilseed crop that belongs to the spurge (Euphorbiaceae) family, which comprises approximately 6,300 species that include cassava (Manihot esculenta), rubber tree (Hevea brasiliensis) and physic nut (Jatropha curcas). It is primarily of economic interest as a source of castor oil, used for the production of high-quality lubricants because of its high proportion of the unusual fatty acid ricinoleic acid. However, castor bean genomics is also relevant to biosecurity as the seeds contain high levels of ricin, a highly toxic, ribosome-inactivating protein. Here we report the draft genome sequence of castor bean (4.6-fold coverage), the first for a member of the Euphorbiaceae. Whereas most of the key genes involved in oil synthesis and turnover are single copy, the number of members of the ricin gene family is larger than previously thought. Comparative genomics analysis suggests the presence of an ancient hexaploidization event that is conserved across the dicotyledonous lineage.

The human malaria parasite Plasmodium vivax is responsible for 25-40% of the approximately 515 million annual cases of malaria worldwide. Although seldom fatal, the parasite elicits severe and incapacitating clinical symptoms and often causes relapses months after a primary infection has cleared. Despite its importance as a major human pathogen, P. vivax is little studied because it cannot be propagated continuously in the laboratory except in non-human primates. We sequenced the genome of P. vivax to shed light on its distinctive biological features, and as a means to drive development of new drugs and vaccines. Here we describe the synteny and isochore structure of P. vivax chromosomes, and show that the parasite resembles other malaria parasites in gene content and metabolic potential, but possesses novel gene families and potential alternative invasion pathways not recognized previously. Completion of the P. vivax genome provides the scientific community with a valuable resource that can be used to advance investigation into this neglected species.

Parasitic nematodes that cause elephantiasis and river blindness threaten hundreds of millions of people in the developing world. We have sequenced the ∼90 megabase (Mb) genome of the human filarial parasite Brugia malayi and predict ∼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 ∼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.

Borrelia burgdorferi is a causative agent of Lyme disease in North America and Eurasia. The first complete genome sequence of B. burgdorferi strain 31, available for more than a decade, has assisted research on the pathogenesis of Lyme disease. Because a single genome sequence is not sufficient to understand the relationship between genotypic and geographic variation and disease phenotype, we determined the whole-genome sequences of 13 additional B. burgdorferi isolates that span the range of natural variation. These sequences should allow improved understanding of pathogenesis and provide a foundation for novel detection, diagnosis, and prevention strategies.

Background Endodontic infections are a leading cause of oro-facial pain and tooth loss in western countries, and may lead to severe life-threatening infections. These infections are polymicrobial with high bacterial diversity. Understanding the spatial transition of microbiota from normal oral cavities through the infected root canal to the acute periapical abscess can improve our knowledge of the pathogenesis of endodontic infections and lead to more effective treatment. We obtained samples from the oral cavity, infected root canal and periapical abscess of 8 patients (5 with localized and 3 with systemic infections). Microbial populations in these samples were analyzed using next-generation sequencing of 16S rRNA amplicons. Bioinformatics tools and statistical tests with rigorous criteria were used to elucidate the spatial transition of the microbiota from normal to diseased sites. Results On average, 10,000 partial 16S rRNA gene sequences were obtained from each sample. All sequences fell into 11 different bacterial phyla. The microbial diversity in root canal and abscess samples was significantly lower than in the oral samples. Streptococcus was the most abundant genus in oral cavities while Prevotella and Fusobacterium were most abundant in diseased samples. The microbiota community structures of root canal and abscess samples were, however, more similar to each other than to the oral cavity microbiota. Using rigorous criteria and novel bioinformatics tools, we found that Granulicatella adiacens, Eubacterium yurii, Prevotella melaninogenica, Prevotella salivae, Streptococcus mitis, and Atopobium rimae were over-represented in diseased samples. Conclusions We used a novel approach and high-throughput methodologies to characterize the microbiota associated normal and diseased oral sites in the same individuals.

Inflammatory bowel diseases (IBD), such as Crohn’s disease, are chronic, immunologically mediated disorders that have severe medical consequences. The current hypothesis is that these diseases are due to an overly aggressive immune response to a subset of commensal enteric bacteria. Studies to date on IBD have suggested that the disorder may be caused by a combination of bacteria and host susceptibility; however the etiologies of these diseases remain an enigma. In this application, we propose to develop and demonstrate the ability to profile Crohn’s disease at an unprecedented molecular level by elucidation of specific biomarkers (bacterial strains, genes, or proteins) that correlate to disease symptoms. To achieve this goal, we will employ a multidisciplinary approach based on metagenomic and metaproteomic molecular tools to elucidate the composition of the commensal microbiota in monozygotic twins that are either healthy or exhibit Crohn’s disease (for concordant, both are diseased; for discordant, one is healthy and one is diseased). The central hypotheses of this proposal are (1) that specific members and/or functional activities of the gastrointestinal (GI) microbiota differ in patients with Crohn’s disease as compared to healthy individuals, and (2) that it will be possible to elucidate microbial signatures which correlate with the occurrence and progression of this disease by integration of data obtained from 16S rRNA based molecular fingerprinting, metagenomics, and metaproteomics approaches. To address these hypotheses, three specific aims are proposed: 1) Obtain data on community gene content (metagenome) in a subset of healthy twins and twins with Crohn’s Disease to assess potential differences in the metabolic capabilities of the gut microbiota associated with CD, 2) Obtain data on community protein content (metaproteome) in a subset of healthy twins and twins with Crohn’s Disease to assess the state of expressed proteins associated with CD, 3) Apply various statistical clustering and classification methods to correlate/associate microbial community composition, gene and protein content with patient metadata, including metabolite profiles and clinical phenotype. The ultimate goal of these efforts is to identify novel biomarkers for non-invasive diagnostics of CD and to eventually identify drug targets (i.e. bacterial strains) for cure or suppression of disease symptoms.

Emerging evidence that the gut microbiota may contribute in important ways to human health and disease has led us and others to hypothesize that both symbiotic and pathological relationships between gut microbes and their host may by key contributors to obesity and the metabolic complications of obesity. Our “Thrifty Microbiome Hypothesis” posits that gut microbiota play a key role in human energy homeostasis. Specifically, constituents of the gut microbial community may introduce a survival advantage to its host in times of nutrient scarcity, promoting positive energy balance by increasing efficiency of nutrient absorption and improving metabolic efficiency and energy storage. However, in the presence of excess nutrients, fat accretion and obesity may result, and in genetically predisposed individuals, increased fat mass may result in preferential abdominal obesity, ectopic fat deposition (liver, muscle), and metabolic complications of obesity (insulin resistance, hypertension, hyperlipidemia). Furthermore, in the presence of excess nutrients, a pathological transition of the gut microbial community may occur, causing leakage of bacterial products into the intestinal lymphatics and portal circulation, thereby inducing an inflammatory state, further aggravating metabolic syndrome traits and accelerating atherosclerosis. This pathological transition and the extent to which antimicrobial leakage occurs and causes inflammatory and other maladaptive sequelae of obesity may also be influenced by host factors, including genetics. In the proposed study, we will directly test the Thrifty Mirobiome Hypothesis by performing detailed genomic and functional assessment of gut microbial communities in intensively phenotyped and genotyped human subjects before and after intentional manipulation of the gut microbiome. To address this hypothesis, we will enroll three age- and sex-matched groups from the Old Order Amish: (i) 50 obese subjects (BMI > 30 kg/m2) with one or more dysmetabolic manifestations (malignant obesity), (ii) 50 obese subjects (BMI > 30 kg/m2) without any dysmetabolic manifestations (benign obesity), and (iii) 50 nonobese subjects (BMI < 25 kg/m2) without any features of the metabolic syndrome and characterize the architecture of the gut microbiota from the subjects enrolled in this study by high-throughput sequencing of 16S rRNA genes. These studies will provide a deeper understanding of the role of gut microbes in terms of ‘who’s there?’, ‘what are they doing?’, and ‘how are they influencing host energy homeostasis, obesity and its metabolic complications?’

Bacterial diversity in endodontic infections has not been sufficiently studied. The use of modern pyrosequencing technology should allow for more comprehensive analysis than traditional Sanger sequencing. This study investigated bacterial diversity in endodontic infections through taxonomic classification based on 16S rRNA gene sequences generated by 454 GS-FLX pyrosequencing and conventional Sanger capillary sequencing technologies. Sequencings were performed on 7 specimens from endodontic infections. On average, 47 vs. 28,590 sequences were obtained per sample for Sanger sequencing vs. pyrosequencing, representing a 600-fold difference in “depth-of-coverage”. Based on Ribosomal Database Project (RDP II) Classifier analysis, pyrosequencing identified 179 bacterial genera in 13 phyla, which was significantly more than Sanger sequencing. The phylum Bacteroidetes was the most prevalent bacterial phylum. These results indicate that bacterial communities in endodontic infections are more diverse than previously demonstrated. In addition, deep-coverage pyrosequencing of the 16S rRNA gene revealed low-abundance micro-organisms with potential clinical implications.

Microbial pathogens have evolved sophisticated mechanisms for evasion of host innate and adaptive immunities. PFam54 is the largest paralogous gene family in the genomes of Borrelia burgdorferi, the Lyme disease bacterium. One member of PFam54, the complement-regulator acquiring surface proteins 1 (BbCrasp-1), is able to abort the alternative pathway of complement activation via binding human complement-regulator factor H (FH). The gene coding for BbCRASP-1 exists in a tandem array of PFam54 genes in the B. burgdorferi genome, a result apparently of repeated gene duplications. To help elucidate the functions of the large number of PFam54 genes, we performed phylogenomic and structural analyses of the PFam54 gene array from ten B. burgdorferi genomes. Analyses based on gene tree, genome synteny, and structural models revealed rapid adaptive evolution of this array through gene duplication, gene loss, and functional diversification. Individual PFam54 genes, however, do not show high intra-population sequence polymorphisms as genes providing evasion from adaptive immunity generally do. PFam54 members able to bind human FH are not monophyletic, suggesting that human FH affinity, however strong, is an incidental rather than main function of these PFam54 proteins. The large number of PFam54 genes existing in any single B. burgdorferi genome may target different innate-immunity proteins of a single host species or the same immune protein of a variety of host species. Genetic variability of the PFam54 gene array suggests that universally present PFam54 lineages such as BBA64, BBA65, BBA66, and BBA73 may be better candidates for the development of broad-spectrum vaccines or drugs than strain-restricted lineages such as BbCRASP-1.

Schistosoma mansoni is responsible for the neglected tropical disease schistosomiasis that affects 210 million people in 76 countries. We report here analysis of the 363 megabase nuclear genome of the blood fluke. It encodes at least 11,809 genes, with an unusual intron size distribution, and novel families of micro-exon genes that undergo frequent alternate splicing. As the first sequenced flatworm, and a representative of the lophotrochozoa, it offers insights into early events in the evolution of the animals, including the development of a body pattern with bilateral symmetry, and the development of tissues into organs. Our analysis has been informed by the need to find new drug targets. The deficits in lipid metabolism that make schistosomes dependent on the host are revealed, while the identification of membrane receptors, ion channels and more than 300 proteases, provide new insights into the biology of the life cycle and novel targets. Bioinformatics approaches have identified metabolic chokepoints while a chemogenomic screen has pinpointed schistosome proteins for which existing drugs may be active. The information generated provides an invaluable resource for the research community to develop much needed new control tools for the treatment and eradication of this important and neglected disease.

Background: Burkholderia mallei is an understudied biothreat agent responsible for glanders which can be lethal in humans and animals. Research with this pathogen has been hampered in part by constraints of Select Agent regulations for safety reasons. Whole genomic sequencing (WGS) is an apt approach to characterize newly discovered or poorly understood microbial pathogens. Methodology/Principal Findings: We performed WGS on a strain of B. mallei, SAVP1, previously pathogenic, that was experimentally infected in 6 equids (4 ponies, 1 mule, 1 donkey), natural hosts, for purposes of producing antibodies. Multiple high inocula were used in some cases. Unexpectedly SAVP1 appeared to be avirulent in the ponies and mule, and attenuated in the donkey, but induced antibodies. We determined the genome sequence of SAVP1 and compared it to a strain that was virulent in horses and a human. In comparison, this phenotypic avirulent SAVP1 strain was missing multiple genes including all the animal type III secretory system (T3SS) complex of genes demonstrated to be essential for virulence in mice and hamster models. The loss of these genes in the SAVP1 strain appears to be the consequence of a multiple gene deletion across insertion sequence (IS) elements in the B. mallei genome. Therefore, the strain by itself is unlikely to revert naturally to its virulent phenotype. There were other genes present in one strain and not the other and vice-versa. Conclusion/Significance: The discovery that this strain of B. mallei was both avirulent in the natural host ponies, and did not possess T3SS associated genes may be fortuitous to advance biodefense research. The deleted virulence-essential T3SS is not likely to be re-acquired naturally. These findings may provide a basis for exclusion of SAVP1 from the Select Agent regulation or at least discussion of what else would be required for exclusion. This exclusion could accelerate research by investigators not possessing BSL-3 facilities and facilitate the production of reagents such as antibodies without the restraints of Select Agent regulation.

We present the genome sequences of a new clinical isolate of the important human pathogen, Aspergillus fumigatus, A1163, and two closely related but rarely pathogenic species, Neosartorya fischeri NRRL181 and Aspergillus clavatus NRRL1. Comparative genomic analysis of A1163 with the recently sequenced A. fumigatus isolate Af293 has identified core, variable and up to 2% unique genes in each genome. While the core genes are 99.8% identical at the nucleotide level, identity for variable genes can be as low 40%. The most divergent loci appear to contain heterokaryon incompatibility (het) genes associated with fungal programmed cell death such as developmental regulator rosA. Cross-species comparison has revealed that 8.5%, 13.5% and 12.6%, respectively, of A. fumigatus, N. fischeri and A. clavatus genes are species-specific. These genes are significantly smaller in size than core genes, contain fewer exons and exhibit a subtelomeric bias. Most of them cluster together in 13 chromosomal islands, which are enriched for pseudogenes, transposons and other repetitive elements. At least 20% of A. fumigatus-specific genes appear to be functional and involved in carbohydrate and chitin catabolism, transport, detoxification, secondary metabolism and other functions that may facilitate the adaptation to heterogeneous environments such as soil or a mammalian host. Contrary to what was suggested previously, their origin cannot be attributed to horizontal gene transfer (HGT), but instead is likely to involve duplication, diversification and differential gene loss (DDL). The role of duplication in the origin of lineage-specific genes is further underlined by the discovery of genomic islands that seem to function as designated “gene dumps” and, perhaps, simultaneously, as “gene factories”., Author Summary Aspergillus is an extremely diverse genus of filamentous ascomycetous fungi (molds) found ubiquitously in soil and decomposing vegetation. Being supreme opportunists, aspergilli have adapted to overcome various chemical, physical, and biological stresses found in heterogeneous environments. While most species in the genus are saprophytes, a surprising number are able to infect wounded plants and animals. Remarkably, the allergic human host also responds abnormally to the aspergilli with lung and sinus disease. The advent of immunosuppressive agents and other medical advances have created a large worldwide pool of human hosts susceptible to some Aspergillus species, including the world's most harmful mold and the causative agent of invasive aspergillosis, Aspergillus fumigatus. In this study, we have used the power of comparative genomics to gain insight into genetic mechanisms that may contribute to the metabolic versatility and pathogenicity of this important human pathogen. Comparison of the genomes of two A. fumigatus clinical isolates and two closely related, but rarely pathogenic species showed that their genomes contain several large isolate- and species-specific chromosomal islands. The metabolic capabilities encoded by these highly labile regions are likely to contribute to their rapid adaptation to heterogeneous environments such as soil or a living host.

Anthrax, caused by the bacterium Bacillus anthracis, is a disease of historical and current importance that is found throughout the world. The basis of its historical transmission is anecdotal and its true global population structure has remained largely cryptic. Seven diverse B. anthracis strains were whole-genome sequenced to identify rare single nucleotide polymorphisms (SNPs), followed by phylogenetic reconstruction of these characters onto an evolutionary model. This analysis identified SNPs that define the major clonal lineages within the species. These SNPs, in concert with 15 variable number tandem repeat (VNTR) markers, were used to subtype a collection of 1,033 B. anthracis isolates from 42 countries to create an extensive genotype data set. These analyses subdivided the isolates into three previously recognized major lineages (A, B, and C), with further subdivision into 12 clonal sub-lineages or sub-groups and, finally, 221 unique MLVA15 genotypes. This rare genomic variation was used to document the evolutionary progression of B. anthracis and to establish global patterns of diversity. Isolates in the A lineage are widely dispersed globally, whereas the B and C lineages occur on more restricted spatial scales. Molecular clock models based upon genome-wide synonymous substitutions indicate there was a massive radiation of the A lineage that occurred in the mid-Holocene (3,064-6,127 ybp). On more recent temporal scales, the global population structure of B. anthracis reflects colonial-era importation of specific genotypes from the Old World into the New World, as well as the repeated industrial importation of diverse genotypes into developed countries via spore-contaminated animal products. These findings indicate humans have played an important role in the evolution of anthrax by increasing the proliferation and dispersal of this now global disease. Finally, the value of global genotypic analysis for investigating bioterrorist-mediated outbreaks of anthrax is demonstrated.

It has been known for decades that human Lyme disease is caused by the three spirochete species Borrelia burgdorferi , Borrelia afzelii , and Borrelia garinii . Recently, Borrelia valaisiana , Borrelia spielmanii , and Borrelia bissettii have been associated with Lyme disease. We report the complete genome sequences of B. valaisiana VS116, B. spielmanii A14S, and B. bissettii DN127.

Human Lyme disease is commonly caused by several species of spirochetes in the Borrelia genus. In Eurasia these species are largely Borrelia afzelii , B. garinii , B. burgdorferi , and B. bavariensis sp. nov. Whole-genome sequencing is an excellent tool for investigating and understanding the influence of bacterial diversity on the pathogenesis and etiology of Lyme disease. We report here the whole-genome sequences of four isolates from two of the Borrelia species that cause human Lyme disease, B. afzelii isolates ACA-1 and PKo and B. garinii isolates PBr and Far04.

Human Lyme disease is caused by a number of related Borrelia burgdorferi sensu lato species. We report here the complete genome sequence of Borrelia sp. isolate SV1 from Finland. This isolate is to date the closest known relative of B. burgdorferi sensu stricto, but it is sufficiently genetically distinct from that species that it and its close relatives warrant its candidacy for new-species status. We suggest that this isolate should be named “ Borrelia finlandensis .”

The pathogenic bacterium Bacillus anthracis has become the subject of intense study as a result of its use in a bioterrorism attack in the United States in September and October 2001. Previous studies suggested that B. anthracis Ames Ancestor, the original Ames fully virulent plasmid-containing isolate, was the ideal reference. This study describes the complete genome sequence of that original isolate, derived from a sample kept in cold storage since 1981.