Consensus coding sequence (CCDS) database: a standardized set of human and mouse protein-coding regions supported by expert curation

The Consensus Coding Sequence (CCDS) project provides a dataset of protein-coding regions that are identically annotated on the human and mouse reference genome assembly in genome annotations produced independently by NCBI and the Ensembl group at EMBL-EBI. This dataset is the product of an international collaboration that includes NCBI, Ensembl, HUGO Gene Nomenclature Committee, Mouse Genome Informatics and University of California, Santa Cruz. Identically annotated coding regions, which are generated using an automated pipeline and pass multiple quality assurance checks, are assigned a stable and tracked identifier (CCDS ID). Additionally, coordinated manual review by expert curators from the CCDS collaboration helps in maintaining the integrity and high quality of the dataset. The CCDS data are available through an interactive web page (https://www.ncbi.nlm.nih.gov/CCDS/CcdsBrowse.cgi) and an FTP site (ftp://ftp.ncbi.nlm.nih.gov/pub/CCDS/). In this paper, we outline the ongoing work, growth and stability of the CCDS dataset and provide updates on new collaboration members and new features added to the CCDS user interface. We also present expert curation scenarios, with specific examples highlighting the importance of an accurate reference genome assembly and the crucial role played by input from the research community.

Whole genome comparison of Aspergillus flavus L-morphotype strain NRRL 3357 (type) and S-morphotype strain AF70

Aspergillus flavus is a saprophytic fungus that infects corn, peanuts, tree nuts and other agriculturally important crops. Once the crop is infected the fungus has the potential to secrete one or more mycotoxins, the most carcinogenic of which is aflatoxin. Aflatoxin contaminated crops are deemed unfit for human or animal consumption, which results in both food and economic losses. Within A. flavus, two morphotypes exist: the S strains (small sclerotia) and L strains (large sclerotia). Significant morphological and physiological differences exist between the two morphotypes. For example, the S-morphotypes produces sclerotia that are smaller (< 400 μm), greater in quantity, and contain higher concentrations of aflatoxin than the L-morphotypes (>400 μm). The morphotypes also differ in pigmentation, pH homeostasis in culture and the number of spores produced. Here we report the first full genome sequence of an A. flavus S morphotype, strain AF70. We provide a comprehensive comparison of the A. flavus S-morphotype genome sequence with a previously sequenced genome of an L-morphotype strain (NRRL 3357), including an in-depth analysis of secondary metabolic clusters and the identification SNPs within their aflatoxin gene clusters.

Mobile elements and mitochondrial genome expansion in the soil fungus and potato pathogen Rhizoctonia solani AG-3

The soil fungus Rhizoctonia solani is an economically important pathogen of agricultural and forestry crops. Here, we present the complete sequence and analysis of the mitochondrial genome of R. solani, field isolate Rhs1AP. The genome (235 849 bp) is the largest mitochondrial genome of a filamentous fungus sequenced to date and exhibits a rich accumulation of introns, novel repeat sequences, homing endonuclease genes, and hypothetical genes. Stable secondary structures exhibited by repeat sequences suggest that they comprise functional, possibly catalytic RNA elements. RNA-Seq expression profiling confirmed that the majority of homing endonuclease genes and hypothetical genes are transcriptionally active. Comparative analysis suggests that the mitochondrial genome of R. solani is an example of a dynamic history of expansion in filamentous fungi.

Sequencing of mitochondrial genomes of nine Aspergillus and Penicillium species identifies mobile introns and accessory genes as main sources of genome size variability

BACKGROUND

The genera Aspergillus and Penicillium include some of the most beneficial as well as the most harmful fungal species such as the penicillin-producer Penicillium chrysogenum and the human pathogen Aspergillus fumigatus, respectively. Their mitochondrial genomic sequences may hold vital clues into the mechanisms of their evolution, population genetics, and biology, yet only a handful of these genomes have been fully sequenced and annotated.

RESULTS

Here we report the complete sequence and annotation of the mitochondrial genomes of six Aspergillus and three Penicillium species: A. fumigatus, A. clavatus, A. oryzae, A. flavus, Neosartorya fischeri (A. fischerianus), A. terreus, P. chrysogenum, P. marneffei, and Talaromyces stipitatus (P. stipitatum). The accompanying comparative analysis of these and related publicly available mitochondrial genomes reveals wide variation in size (25-36 Kb) among these closely related fungi. The sources of genome expansion include group I introns and accessory genes encoding putative homing endonucleases, DNA and RNA polymerases (presumed to be of plasmid origin) and hypothetical proteins. The two smallest sequenced genomes (A. terreus and P. chrysogenum) do not contain introns in protein-coding genes, whereas the largest genome (T. stipitatus), contains a total of eleven introns. All of the sequenced genomes have a group I intron in the large ribosomal subunit RNA gene, suggesting that this intron is fixed in these species. Subsequent analysis of several A. fumigatus strains showed low intraspecies variation. This study also includes a phylogenetic analysis based on 14 concatenated core mitochondrial proteins. The phylogenetic tree has a different topology from published multilocus trees, highlighting the challenges still facing the Aspergillus systematics.

CONCLUSIONS

The study expands the genomic resources available to fungal biologists by providing mitochondrial genomes with consistent annotations for future genetic, evolutionary and population studies. Despite the conservation of the core genes, the mitochondrial genomes of Aspergillus and Penicillium species examined here exhibit significant amount of interspecies variation. Most of this variation can be attributed to accessory genes and mobile introns, presumably acquired by horizontal gene transfer of mitochondrial plasmids and intron homing.

Genome sequencing and analysis of the filamentous fungus Penicillium chrysogenum

Industrial penicillin production with the filamentous fungus Penicillium chrysogenum is based on an unprecedented effort in microbial strain improvement. To gain more insight into penicillin synthesis, we sequenced the 32.19 Mb genome of P. chrysogenum Wisconsin54-1255 and identified numerous genes responsible for key steps in penicillin production. DNA microarrays were used to compare the transcriptomes of the sequenced strain and a penicillinG high-producing strain, grown in the presence and absence of the side-chain precursor phenylacetic acid. Transcription of genes involved in biosynthesis of valine, cysteine and alpha-aminoadipic acid-precursors for penicillin biosynthesis-as well as of genes encoding microbody proteins, was increased in the high-producing strain. Some gene products were shown to be directly controlling beta-lactam output. Many key cellular transport processes involving penicillins and intermediates remain to be characterized at the molecular level. Genes predicted to encode transporters were strongly overrepresented among the genes transcriptionally upregulated under conditions that stimulate penicillinG production, illustrating potential for future genomics-driven metabolic engineering.

What can comparative genomics tell us about species concepts in the genus Aspergillus?

Understanding the nature of species" boundaries is a fundamental question in evolutionary biology. The availability of genomes from several species of the genus Aspergillus allows us for the first time to examine the demarcation of fungal species at the whole-genome level. Here, we examine four case studies, two of which involve intraspecific comparisons, whereas the other two deal with interspecific genomic comparisons between closely related species. These four comparisons reveal significant variation in the nature of species boundaries across Aspergillus. For example, comparisons between A. fumigatus and Neosartorya fischeri (the teleomorph of A. fischerianus) and between A. oryzae and A. flavus suggest that measures of sequence similarity and species-specific genes are significantly higher for the A. fumigatus - N. fischeri pair. Importantly, the values obtained from the comparison between A. oryzae and A. flavus are remarkably similar to those obtained from an intra-specific comparison of A. fumigatus strains, giving support to the proposal that A. oryzae represents a distinct ecotype of A. flavus and not a distinct species. We argue that genomic data can aid Aspergillus taxonomy by serving as a source of novel and unprecedented amounts of comparative data, as a resource for the development of additional diagnostic tools, and finally as a knowledge database about the biological differences between strains and species.

Bioinformatics-enabled identification of the HrpL regulon and type III secretion system effector proteins of Pseudomonas syringae pv. phaseolicola 1448A

The ability of Pseudomonas syringae pv. phaseolicola to cause halo blight of bean is dependent on its ability to translocate effector proteins into host cells via the hypersensitive response and pathogenicity (Hrp) type III secretion system (T3SS). To identify genes encoding type III effectors and other potential virulence factors that are regulated by the HrpL alternative sigma factor, we used a hidden Markov model, weight matrix model, and type III targeting-associated patterns to search the genome of P. syringae pv. phaseolicola 1448A, which recently was sequenced to completion. We identified 44 high-probability putative Hrp promoters upstream of genes encoding the core T3SS machinery, 27 candidate effectors and related T3SS substrates, and 10 factors unrelated to the Hrp system. The expression of 13 of these candidate HrpL regulon genes was analyzed by real-time polymerase chain reaction, and all were found to be upregulated by HrpL. Six of the candidate type III effectors were assayed for T3SS-dependent translocation into plant cells using the Bordetella pertussis calmodulin-dependent adenylate cyclase (Cya) translocation reporter, and all were translocated. PSPPH1855 (ApbE-family protein) and PSPPH3759 (alcohol dehydrogenase) have no apparent T3SS-related function; however, they do have homologs in the model strain P. syringae pv. tomato DC3000 (PSPTO2105 and PSPTO0834, respectively) that are similarly upregulated by HrpL. Mutations were constructed in the DC3000 homologs and found to reduce bacterial growth in host Arabidopsis leaves. These results establish the utility of the bioinformatic or candidate gene approach to identifying effectors and other genes relevant to pathogenesis in P. syringae genomes.

The Arabidopsis thaliana JASMONATE INSENSITIVE 1 gene is required for suppression of salicylic acid-dependent defenses during infection by Pseudomonas syringae

Many plant pathogens suppress antimicrobial defenses using virulence factors that modulate endogenous host defenses. The Pseudomonas syringae phytotoxin coronatine (COR) is believed to promote virulence by acting as a jasmonate analog, because COR-insensitive 1 (coil) Arabidopsis thaliana and tomato mutants are impaired in jasmonate signaling and exhibit reduced susceptibility to P. syringae. To further investigate the role of jasmonate signaling in disease development, we analyzed several jasmonate-insensitive A. thaliana mutants for susceptibility to P. syringae pv. tomato strain DC3000 and sensitivity to COR. Jasmonate-insensitive 1 (jin1) mutants exhibit both reduced susceptibility to P. syringae pv. tomato DC3000 and reduced sensitivity to COR, whereas jasmonate-resistant 1 (jar1) plants exhibit wild-type responses to both COR and P. syringae pv. tomato DC3000. A jin1 jar1 double mutant does not exhibit enhanced jasmonate insensitivity, suggesting that JIN1 functions downstream of jasmonic acid-amino acid conjugates synthesized by JAR1. Reduced disease susceptibility in jin1 mutants is correlated with elevated expression of pathogenesis-related 1 (PR-1) and is dependent on accumulation of salicylic acid (SA). We also show that JIN1 is required for normal P. syringae pv. tomato DC3000 symptom development through an SA-independent mechanism. Thus, P. syringae pv. tomato DC3000 appears to utilize COR to manipulate JIN1-dependent jasmonate signaling both to suppress SA-mediated defenses and to promote symptom development.

UV Irradiation and Desiccation Modulate the Three-dimensional Extracellular Matrix of Nostoc commune (Cyanobacteria)

Cyanobacterium Nostoc commune can tolerate the simultaneous stresses of desiccation, UV irradiation, and oxidation. Acidic WspA, of approximately 33.6 kDa, is secreted to the three-dimensional extracellular matrix and accounts for greater than 70% of the total soluble protein. The wspA gene of N. commune strain DRH1 was cloned and found in a single genomic copy, in a monocistronic operon. Transcription of wspA and sodF (superoxide dismutase), and synthesis and secretion of WspA, were induced upon desiccation or UV-A/B irradiation of cells. Recombinant WspA binds the UV-A/B absorbing pigment scytonemin through non-covalent interactions. WspA peptide polymorphism, and heterogeneity of multiple wspA sequences within cells of a single colony, account for distinct WspA isoforms. WspA has no similarity to entries in the sequence databases and wspA, a possible xenolog, is restricted to a subset of strains in the "form species" N. commune characterized through group I intron phylogeny. We hypothesize that WspA plays a central role in the global stress response of N. commune through modulation of the structure and function of the three-dimensional extracellular matrix, particularly the transport, distribution, and/or macromolecular architecture of mycosporine and scytonemin UV-A/B absorbing pigment complexes.

Whole-genome sequence analysis of Pseudomonas syringae pv. phaseolicola 1448A reveals divergence among pathovars in genes involved in virulence and transposition

Pseudomonas syringae pv. phaseolicola, a gram-negative bacterial plant pathogen, is the causal agent of halo blight of bean. In this study, we report on the genome sequence of P. syringae pv. phaseolicola isolate 1448A, which encodes 5,353 open reading frames (ORFs) on one circular chromosome (5,928,787 bp) and two plasmids (131,950 bp and 51,711 bp). Comparative analyses with a phylogenetically divergent pathovar, P. syringae pv. tomato DC3000, revealed a strong degree of conservation at the gene and genome levels. In total, 4,133 ORFs were identified as putative orthologs in these two pathovars using a reciprocal best-hit method, with 3,941 ORFs present in conserved, syntenic blocks. Although these two pathovars are highly similar at the physiological level, they have distinct host ranges; 1448A causes disease in beans, and DC3000 is pathogenic on tomato and Arabidopsis. Examination of the complement of ORFs encoding virulence, fitness, and survival factors revealed a substantial, but not complete, overlap between these two pathovars. Another distinguishing feature between the two pathovars is their distinctive sets of transposable elements. With access to a fifth complete pseudomonad genome sequence, we were able to identify 3,567 ORFs that likely comprise the core Pseudomonas genome and 365 ORFs that are P. syringae specific.

Lineage-specific regions in Pseudomonas syringae pv tomato DC3000

SUMMARY Comparative analyses of the chromosome of Pseudomonas syringae pv tomato DC3000 with the finished, complete genomes of Pseudomonas aeruginosa PAO1, an animal pathogen, and the non-pathogenic soil inhabitant Pseudomonas putida KT2440 revealed a high degree of sequence conservation in genes involved in 'housekeeping functions'. However, divergence is present among these three fluorescent pseudomonads, yielding 'suites' of species-specific genes that may provide the genetic basis for adaptation to an ecological niche and lifestyle. For DC3000, 1053 genes located on the chromosome were specific to DC3000 and not present in PAO1 or KT2440. The majority of these DC3000-specific genes either lack a known function or are mobile genetic elements. However, these genes do share features among themselves such as association with regions of atypical trinucleotides, unusual G+C content and localization within large tracts of DC3000-specific sequence, suggestive of lateral gene transfer events. Indeed, a comparison of syntenic blocks among these three complete Pseudomonas genomes revealed that a substantial portion (533) of the DC3000-specific chromosomal genes (1053) were located in lineage-specific regions (defined as being larger than 2 kb and enriched in mobile genetic elements and/or genes specific to DC3000 in this three-way comparison). A large proportion of mobile genetic elements (199 of 318 genes; 63%), which are highly enriched in DC3000, were present within such regions. Similarly, most of the genes encoding type III secretion system virulence effectors were located in lineage-specific regions. Consistent with the plasticity of the DC3000 genome, a putative chromosomal inversion mediated by identical copies of ISPsy6 involving 2838 kb (44%) of the DC3000 genome was detected. These data suggest that a substantial portion of the differentiation of DC3000, a plant pathogen, from an animal pathogen and a soil inhabitant has involved transfer of a large number of novel genes coupled with amplification of mobile genetic elements.

The complete genome sequence of the Arabidopsis and tomato pathogen Pseudomonas syringae pv. tomato DC3000

We report the complete genome sequence of the model bacterial pathogen Pseudomonas syringae pathovar tomato DC3000 (DC3000), which is pathogenic on tomato and Arabidopsis thaliana. The DC3000 genome (6.5 megabases) contains a circular chromosome and two plasmids, which collectively encode 5,763 ORFs. We identified 298 established and putative virulence genes, including several clusters of genes encoding 31 confirmed and 19 predicted type III secretion system effector proteins. Many of the virulence genes were members of paralogous families and also were proximal to mobile elements, which collectively comprise 7% of the DC3000 genome. The bacterium possesses a large repertoire of transporters for the acquisition of nutrients, particularly sugars, as well as genes implicated in attachment to plant surfaces. Over 12% of the genes are dedicated to regulation, which may reflect the need for rapid adaptation to the diverse environments encountered during epiphytic growth and pathogenesis. Comparative analyses confirmed a high degree of similarity with two sequenced pseudomonads, Pseudomonas putida and Pseudomonas aeruginosa, yet revealed 1,159 genes unique to DC3000, of which 811 lack a known function.

Identification of novel hrp-regulated genes through functional genomic analysis of the Pseudomonas syringae pv. tomato DC3000 genome

Pseudomonas syringae pv. tomato (Pst) strain DC3000 infects the model plants Arabidopsis thaliana and tomato, causing disease symptoms characterized by necrotic lesions surrounded by chlorosis. One mechanism used by Pst DC3000 to infect host plants is the type III protein secretion system, which is thought to deliver multiple effector proteins to the plant cell. The exact number of type III effectors in Pst DC3000 or any other plant pathogenic bacterium is not known. All known type III effector genes of P. syringae are regulated by HrpS, an NtrC family protein, and the HrpL alternative sigma factor, which presumably binds to a conserved cis element (called the "hrp box") in the promoters of type III secretion-associated genes. In this study, we designed a search motif based on the promoter sequences conserved in 12 published hrp operons and putative effector genes in Pst DC3000. Seventy-three predicted genes were retrieved from the January 2001 release of the Pst DC3000 genome sequence, which had 95% genome coverage. The expression of the 73 genes was analysed by microarray and Northern blotting, revealing 24 genes/operons (including eight novel genes), the expression of which was consistently higher in hrp-inducing minimal medium than in nutrient-rich Luria-Bertani broth. Expression of all eight genes was dependent on the hrpS gene. Most were also dependent on the hrpL gene, but at least one was dependent on the hrpS gene, but not on the hrpL gene. An AvrRpt2-based type III translocation assay provides evidence that some of the hrpS-regulated novel genes encode putative effector proteins.

Identification of Pseudomonas syringae pv. tomato genes induced during infection of Arabidopsis thaliana

Phytopathogenic bacteria possess a large number of genes that allow them to grow and cause disease on plants. Many of these genes should be induced when the bacteria come in contact with plant tissue. We used a modified in vivo expression technology (IVET) approach to identify genes from the plant pathogen Pseudomonas syringae pv. tomato that are induced upon infection of Arabidopsis thaliana and isolated over 500 in planta-expressed (ipx) promoter fusions. Sequence analysis of 79 fusions revealed several known and potential virulence genes, including hrp/hrc, avr and coronatine biosynthetic genes. In addition, we identified metabolic genes presumably important for adaptation to growth in plant tissue, as well as several genes with unknown function that may encode novel virulence factors. Many ipx fusions, including several corresponding to novel genes, are dependent on HrpL, an alternative RNA polymerase sigma factor that regulates the expression of virulence genes. Expression analysis indicated that several ipx fusions are strongly induced upon inoculation into plant tissue. Disruption of one ipx gene, conserved effector locus (CEL) orf1, encoding a putative lytic murein transglycosylase, resulted in decreased virulence of P. syringae. Our results demonstrate that this screen can be used successfully to isolate genes that are induced in planta, including many novel genes potentially involved in pathogenesis.