A novel hybrid gene prediction method employing protein multiple sequence alignments

The pangenome of hexaploid bread wheat

There is an increasing understanding that variation in gene presence-absence plays an important role in the heritability of agronomic traits; however, there have been relatively few studies on variation in gene presence-absence in crop species. Hexaploid wheat is one of the most important food crops in the world and intensive breeding has reduced the genetic diversity of elite cultivars. Major efforts have produced draft genome assemblies for the cultivar Chinese Spring, but it is unknown how well this represents the genome diversity found in current modern elite cultivars. In this study we build an improved reference for Chinese Spring and explore gene diversity across 18 wheat cultivars. We predict a pangenome size of 140 500 ± 102 genes, a core genome of 81 070 ± 1631 genes and an average of 128 656 genes in each cultivar. Functional annotation of the variable gene set suggests that it is enriched for genes that may be associated with important agronomic traits. In addition to variation in gene presence, more than 36 million intervarietal single nucleotide polymorphisms were identified across the pangenome. This study of the wheat pangenome provides insight into genome diversity in elite wheat as a basis for genomics-based improvement of this important crop. A wheat pangenome, GBrowse, is available at http://appliedbioinformatics.com.au/cgi-bin/gb2/gbrowse/WheatPan/, and data are available to download from http://wheatgenome.info/wheat_genome_databases.php.

Comparative analysis of the predicted secretomes of Rosaceae scab pathogens Venturia inaequalis and V. pirina reveals expanded effector families and putative determinants of host range

BACKGROUND

Fungal plant pathogens belonging to the genus Venturia cause damaging scab diseases of members of the Rosaceae. In terms of economic impact, the most important of these are V. inaequalis, which infects apple, and V. pirina, which is a pathogen of European pear. Given that Venturia fungi colonise the sub-cuticular space without penetrating plant cells, it is assumed that effectors that contribute to virulence and determination of host range will be secreted into this plant-pathogen interface. Thus the predicted secretomes of a range of isolates of Venturia with distinct host-ranges were interrogated to reveal putative proteins involved in virulence and pathogenicity.

RESULTS

Genomes of Venturia pirina (one European pear scab isolate) and Venturia inaequalis (three apple scab, and one loquat scab, isolates) were sequenced and the predicted secretomes of each isolate identified. RNA-Seq was conducted on the apple-specific V. inaequalis isolate Vi1 (in vitro and infected apple leaves) to highlight virulence and pathogenicity components of the secretome. Genes encoding over 600 small secreted proteins (candidate effectors) were identified, most of which are novel to Venturia, with expansion of putative effector families a feature of the genus. Numerous genes with similarity to Leptosphaeria maculans AvrLm6 and the Verticillium spp. Ave1 were identified. Candidates for avirulence effectors with cognate resistance genes involved in race-cultivar specificity were identified, as were putative proteins involved in host-species determination. Candidate effectors were found, on average, to be in regions of relatively low gene-density and in closer proximity to repeats (e.g. transposable elements), compared with core eukaryotic genes.

CONCLUSIONS

Comparative secretomics has revealed candidate effectors from Venturia fungal plant pathogens that attack pome fruit. Effectors that are putative determinants of host range were identified; both those that may be involved in race-cultivar and host-species specificity. Since many of the effector candidates are in close proximity to repetitive sequences this may point to a possible mechanism for the effector gene family expansion observed and a route to diversification via transposition and repeat-induced point mutation.

QTL mapping for flowering-time and photoperiod insensitivity of cotton Gossypium darwinii Watt

Most wild and semi-wild species of the genus Gossypium are exhibit photoperiod-sensitive flowering. The wild germplasm cotton is a valuable source of genes for genetic improvement of modern cotton cultivars. A bi-parental cotton population segregating for photoperiodic flowering was developed by crossing a photoperiod insensitive irradiation mutant line with its pre-mutagenesis photoperiodic wild-type G. darwinii Watt genotype. Individuals from the F2 and F3 generations were grown with their parental lines and F1 hybrid progeny in the long day and short night summer condition (natural day-length) of Uzbekistan to evaluate photoperiod sensitivity, i.e., flowering-time during the seasons 2008-2009. Through genotyping the individuals of this bi-parental population segregating for flowering-time, linkage maps were constructed using 212 simple-sequence repeat (SSR) and three cleaved amplified polymorphic sequence (CAPS) markers. Six QTLs directly associated with flowering-time and photoperiodic flowering were discovered in the F2 population, whereas eight QTLs were identified in the F3 population. Two QTLs controlling photoperiodic flowering and duration of flowering were common in both populations. In silico annotations of the flanking DNA sequences of mapped SSRs from sequenced cotton (G. hirsutum L.) genome database has identified several potential 'candidate' genes that are known to be associated with regulation of flowering characteristics of plants. The outcome of this research will expand our understanding of the genetic and molecular mechanisms of photoperiodic flowering. Identified markers should be useful for marker-assisted selection in cotton breeding to improve early flowering characteristics.

Primer Sets Developed for Functional Genes Reveal Shifts in Functionality of Fungal Community in Soils

Phylogenetic diversity of soil microbes is a hot topic at the moment. However, the molecular tools for the assessment of functional diversity in the fungal community are less developed than tools based on genes encoding the ribosomal operon. Here 20 sets of primers targeting genes involved mainly in carbon cycling were designed and/or validated and the functioning of soil fungal communities along a chronosequence of land abandonment from agriculture was evaluated using them. We hypothesized that changes in fungal community structure during secondary succession would lead to difference in the types of genes present in soils and that these changes would be directional. We expected an increase in genes involved in degradation of recalcitrant organic matter in time since agriculture. Out of the investigated genes, the richness of the genes related to carbon cycling was significantly higher in fields abandoned for longer time. The composition of six of the genes analyzed revealed significant differences between fields abandoned for shorter and longer time. However, all genes revealed significant variance over the fields studied, and this could be related to other parameters than the time since agriculture such as pH, organic matter, and the amount of available nitrogen. Contrary to our initial hypothesis, the genes significantly different between fields were not related to the decomposition of more recalcitrant matter but rather involved in degradation of cellulose and hemicellulose.

Simultaneous gene finding in multiple genomes

MOTIVATION

As the tree of life is populated with sequenced genomes ever more densely, the new challenge is the accurate and consistent annotation of entire clades of genomes. We address this problem with a new approach to comparative gene finding that takes a multiple genome alignment of closely related species and simultaneously predicts the location and structure of protein-coding genes in all input genomes, thereby exploiting negative selection and sequence conservation. The model prefers potential gene structures in the different genomes that are in agreement with each other, or-if not-where the exon gains and losses are plausible given the species tree. We formulate the multi-species gene finding problem as a binary labeling problem on a graph. The resulting optimization problem is NP hard, but can be efficiently approximated using a subgradient-based dual decomposition approach.

RESULTS

The proposed method was tested on whole-genome alignments of 12 vertebrate and 12 Drosophila species. The accuracy was evaluated for human, mouse and Drosophila melanogaster and compared to competing methods. Results suggest that our method is well-suited for annotation of (a large number of) genomes of closely related species within a clade, in particular, when RNA-Seq data are available for many of the genomes. The transfer of existing annotations from one genome to another via the genome alignment is more accurate than previous approaches that are based on protein-spliced alignments, when the genomes are at close to medium distances.

AVAILABILITY AND IMPLEMENTATION

The method is implemented in C ++ as part of Augustus and available open source at http://bioinf.uni-greifswald.de/augustus/ CONTACT: stefaniekoenig@ymail.com or mario.stanke@uni-greifswald.deSupplementary information: Supplementary data are available at Bioinformatics online.

Draft genome of the leopard gecko, Eublepharis macularius

Background

Geckos are among the most species-rich reptile groups and the sister clade to all other lizards and snakes. Geckos possess a suite of distinctive characteristics, including adhesive digits, nocturnal activity, hard, calcareous eggshells, and a lack of eyelids. However, one gecko clade, the Eublepharidae, appears to be the exception to most of these ‘rules’ and lacks adhesive toe pads, has eyelids, and lays eggs with soft, leathery eggshells. These differences make eublepharids an important component of any investigation into the underlying genomic innovations contributing to the distinctive phenotypes in ‘typical’ geckos.

Findings

We report high-depth genome sequencing, assembly, and annotation for a male leopard gecko, Eublepharis macularius (Eublepharidae). Illumina sequence data were generated from seven insert libraries (ranging from 170 to 20 kb), representing a raw sequencing depth of 136X from 303 Gb of data, reduced to 84X and 187 Gb after filtering. The assembled genome of 2.02 Gb was close to the 2.23 Gb estimated by k-mer analysis. Scaffold and contig N50 sizes of 664 and 20 kb, respectively, were comparable to the previously published Gekko japonicus genome. Repetitive elements accounted for 42 % of the genome. Gene annotation yielded 24,755 protein-coding genes, of which 93 % were functionally annotated. CEGMA and BUSCO assessment showed that our assembly captured 91 % (225 of 248) of the core eukaryotic genes, and 76 % of vertebrate universal single-copy orthologs.

Conclusions

Assembly of the leopard gecko genome provides a valuable resource for future comparative genomic studies of geckos and other squamate reptiles.

Electronic supplementary material

The online version of this article (doi:10.1186/s13742-016-0151-4) contains supplementary material, which is available to authorized users.

Genome Sequencing of the Pyruvate-producing Strain Candida glabrata CCTCC M202019 and Genomic Comparison with Strain CBS138

Candida glabrata CCTCC M202019 as an industrial yeast strain that is widely used to produce α-oxocarboxylic acid. Strain M202019 has been proven to have a higher pyruvate-producing capacity than the reference strain CBS138. To characterize the genotype of the M202019 strain, we generated a draft sequence of its genome, which has a size of 12.1 Mbp and a GC content of 38.47%. Evidence accumulated during genome annotation suggests that strain M202019 has strong capacities for glucose transport and pyruvate biosynthesis, defects in pyruvate catabolism, as well as variations in genes involved in nutrient and dicarboxylic acid transport, oxidative phosphorylation, and other relevant aspects of carbon metabolism, which might promote pyruvate accumulation. In addition to differences in its central carbon metabolism, a genomic analysis revealed genetic differences in adhesion metabolism. Forty-nine adhesin-like proteins of strain M202019 were identified classified into seven subfamilies. Decreased amounts of adhesive proteins, and deletions or changes of low-complexity repeats and functional domains might lead to lower adhesion and reduced pathogenicity. Further virulence experiments validated the biological safety of strain M202019. Analysis of the C. glabrata CCTCC M202019 genome sequence provides useful insights into its genetic context, physical characteristics, and potential metabolic capacity.

A De Novo Genome Sequence Assembly of the Arabidopsis thaliana Accession Niederzenz-1 Displays Presence/Absence Variation and Strong Synteny

Arabidopsis thaliana is the most important model organism for fundamental plant biology. The genome diversity of different accessions of this species has been intensively studied, for example in the 1001 genome project which led to the identification of many small nucleotide polymorphisms (SNPs) and small insertions and deletions (InDels). In addition, presence/absence variation (PAV), copy number variation (CNV) and mobile genetic elements contribute to genomic differences between A. thaliana accessions. To address larger genome rearrangements between the A. thaliana reference accession Columbia-0 (Col-0) and another accession of about average distance to Col-0, we created a de novo next generation sequencing (NGS)-based assembly from the accession Niederzenz-1 (Nd-1). The result was evaluated with respect to assembly strategy and synteny to Col-0. We provide a high quality genome sequence of the A. thaliana accession (Nd-1, LXSY01000000). The assembly displays an N50 of 0.590 Mbp and covers 99% of the Col-0 reference sequence. Scaffolds from the de novo assembly were positioned on the basis of sequence similarity to the reference. Errors in this automatic scaffold anchoring were manually corrected based on analyzing reciprocal best BLAST hits (RBHs) of genes. Comparison of the final Nd-1 assembly to the reference revealed duplications and deletions (PAV). We identified 826 insertions and 746 deletions in Nd-1. Randomly selected candidates of PAV were experimentally validated. Our Nd-1 de novo assembly allowed reliable identification of larger genic and intergenic variants, which was difficult or error-prone by short read mapping approaches alone. While overall sequence similarity as well as synteny is very high, we detected short and larger (affecting more than 100 bp) differences between Col-0 and Nd-1 based on bi-directional comparisons. The de novo assembly provided here and additional assemblies that will certainly be published in the future will allow to describe the pan-genome of A. thaliana.

Identification of pathogenicity-related genes in Fusarium oxysporum f. sp. cepae

Pathogenic isolates of Fusarium oxysporum, distinguished as formae speciales (f. spp.) on the basis of their host specificity, cause crown rots, root rots and vascular wilts on many important crops worldwide. Fusarium oxysporum f. sp. cepae (FOC) is particularly problematic to onion growers worldwide and is increasing in prevalence in the UK. We characterized 31 F. oxysporum isolates collected from UK onions using pathogenicity tests, sequencing of housekeeping genes and identification of effectors. In onion seedling and bulb tests, 21 isolates were pathogenic and 10 were non-pathogenic. The molecular characterization of these isolates, and 21 additional isolates comprising other f. spp. and different Fusarium species, was carried out by sequencing three housekeeping genes. A concatenated tree separated the F. oxysporum isolates into six clades, but did not distinguish between pathogenic and non-pathogenic isolates. Ten putative effectors were identified within FOC, including seven Secreted In Xylem (SIX) genes first reported in F. oxysporum f. sp. lycopersici. Two highly homologous proteins with signal peptides and RxLR motifs (CRX1/CRX2) and a gene with no previously characterized domains (C5) were also identified. The presence/absence of nine of these genes was strongly related to pathogenicity against onion and all were shown to be expressed in planta. Different SIX gene complements were identified in other f. spp., but none were identified in three other Fusarium species from onion. Although the FOC SIX genes had a high level of homology with other f. spp., there were clear differences in sequences which were unique to FOC, whereas CRX1 and C5 genes appear to be largely FOC specific.

A Phylogenomic Framework to Study the Diversity and Evolution of Stramenopiles (=Heterokonts)

Stramenopiles or heterokonts constitute one of the most speciose and diverse clades of protists. It includes ecologically important algae (such as diatoms or large multicellular brown seaweeds), as well as heterotrophic (e.g., bicosoecids, MAST groups) and parasitic (e.g., Blastocystis, oomycetes) species. Despite their evolutionary and ecological relevance, deep phylogenetic relationships among stramenopile groups, inferred mostly from small-subunit rDNA phylogenies, remain unresolved, especially for the heterotrophic taxa. Taking advantage of recently released stramenopile transcriptome and genome sequences, as well as data from the genomic assembly of the MAST-3 species Incisomonas marina generated in our laboratory, we have carried out the first extensive phylogenomic analysis of stramenopiles, including representatives of most major lineages. Our analyses, based on a large data set of 339 widely distributed proteins, strongly support a root of stramenopiles lying between two clades, Bigyra and Gyrista (Pseudofungi plus Ochrophyta). Additionally, our analyses challenge the Phaeista-Khakista dichotomy of photosynthetic stramenopiles (ochrophytes) as two groups previously considered to be part of the Phaeista (Pelagophyceae and Dictyochophyceae), branch with strong support with the Khakista (Bolidophyceae and Diatomeae). We propose a new classification of ochrophytes within the two groups Chrysista and Diatomista to reflect the new phylogenomic results. Our stramenopile phylogeny provides a robust phylogenetic framework to investigate the evolution and diversification of this group of ecologically relevant protists.

Draft Genome Sequence of Capniomyces stellatus, the Obligate Gut Fungal Symbiont of Stonefly

Capniomyces stellatus is a host-specific endosymbiotic fungus, living in the hindgut of stoneflies (especially in Allocapnia). Here, we present the first draft genome sequence of the fungus, as well as the ab initio gene prediction and function analyses, which will facilitate the study and comparative analyses of insect-associated fungi.

The Genome of Undifilum oxytropis Provides Insights into Swainsonine Biosynthesis and Locoism

Undifilum oxytropis is a fungal endophyte of locoweeds. It produces swainsonine, which is the principal toxic ingredient of locoweeds. However, the genes, pathways and mechanisms of swainsonine biosynthesis are not known. In this study, the genome of U. oxytropis was firstly sequenced and assembled into a 70.05 megabases (Mb) draft genome, which encoded 11,057 protein-coding genes, and 54% of them were similar to current publicly available sequences. U. oxytropis genes were annotated and 164 putative genes were annotated into enzymes, such as Saccharopine dehydrogenase, Saccharopine oxidase, and Pyrroline-5-carboxylate reductase, hypothesized to be involved in the biosynthesis pathway of swainsonine. The genome sequence and gene annotation of U. oxytropis will provide new insights into functional analyses. The characterization of genes in swainsonine biosynthesis will greatly facilitate locoweed poisoning research and help direct locoism management.

De Novo Assembly and Characterization of the Transcriptome of Grasshopper Shirakiacris shirakii

Background: The grasshopper Shirakiacris shirakii is an important agricultural pest and feeds mainly on gramineous plants, thereby causing economic damage to a wide range of crops. However, genomic information on this species is extremely limited thus far, and transcriptome data relevant to insecticide resistance and pest control are also not available. Methods: The transcriptome of S. shirakii was sequenced using the Illumina HiSeq platform, and we de novo assembled the transcriptome. Results: Its sequencing produced a total of 105,408,878 clean reads, and the de novo assembly revealed 74,657 unigenes with an average length of 680 bp and N50 of 1057 bp. A total of 28,173 unigenes were annotated for the NCBI non-redundant protein sequences (Nr), NCBI non-redundant nucleotide sequences (Nt), a manually-annotated and reviewed protein sequence database (Swiss-Prot), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Based on the Nr annotation results, we manually identified 79 unigenes encoding cytochrome P450 monooxygenases (P450s), 36 unigenes encoding carboxylesterases (CarEs) and 36 unigenes encoding glutathione S-transferases (GSTs) in S. shirakii. Core RNAi components relevant to miroRNA, siRNA and piRNA pathways, including Pasha, Loquacious, Argonaute-1, Argonaute-2, Argonaute-3, Zucchini, Aubergine, enhanced RNAi-1 and Piwi, were expressed in S. shirakii. We also identified five unigenes that were homologous to the Sid-1 gene. In addition, the analysis of differential gene expressions revealed that a total of 19,764 unigenes were up-regulated and 4185 unigenes were down-regulated in larvae. In total, we predicted 7504 simple sequence repeats (SSRs) from 74,657 unigenes. Conclusions: The comprehensive de novo transcriptomic data of S. shirakii will offer a series of valuable molecular resources for better studying insecticide resistance, RNAi and molecular marker discovery in the transcriptome.

Genome and Transcriptome Sequences Reveal the Specific Parasitism of the Nematophagous Purpureocillium lilacinum 36-1

Purpureocillium lilacinum is a promising nematophagous ascomycete able to adapt diverse environments and it is also an opportunistic fungus that infects humans. A microbial inoculant of P. lilacinum has been registered to control plant parasitic nematodes. However, the molecular mechanism of the toxicological processes is still unclear because of the relatively few reports on the subject. In this study, using Illumina paired-end sequencing, the draft genome sequence and the transcriptome of P. lilacinum strain 36-1 infecting nematode-eggs were determined. Whole genome alignment indicated that P. lilacinum 36-1 possessed a more dynamic genome in comparison with P. lilacinum India strain. Moreover, a phylogenetic analysis showed that the P. lilacinum 36-1 had a closer relation to entomophagous fungi. The protein-coding genes in P. lilacinum 36-1 occurred much more frequently than they did in other fungi, which was a result of the depletion of repeat-induced point mutations (RIP). Comparative genome and transcriptome analyses revealed the genes that were involved in pathogenicity, particularly in the recognition, adhesion of nematode-eggs, downstream signal transduction pathways and hydrolase genes. By contrast, certain numbers of cellulose and xylan degradation genes and a lack of polysaccharide lyase genes showed the potential of P. lilacinum 36-1 as an endophyte. Notably, the expression of appressorium-formation and antioxidants-related genes exhibited similar infection patterns in P. lilacinum strain 36-1 to those of the model entomophagous fungi Metarhizium spp. These results uncovered the specific parasitism of P. lilacinum and presented the genes responsible for the infection of nematode-eggs.

Genomic characterization of plant cell wall degrading enzymes and in silico analysis of xylanases and polygalacturonases of Fusarium virguliforme

BACKGROUND

Plant cell wall degrading enzymes (PCWDEs) are a subset of carbohydrate-active enzymes (CAZy) produced by plant pathogens to degrade plant cell walls. To counteract PCWDEs, plants release PCWDEs inhibitor proteins (PIPs) to reduce their impact. Several transgenic plants expressing exogenous PIPs that interact with fungal glycoside hydrolase (GH)11-type xylanases or GH28-type polygalacturonase (PG) have been shown to enhance disease resistance. However, many plant pathogenic Fusarium species were reported to escape PIPs inhibition. Fusarium virguliforme is a soilborne pathogen that causes soybean sudden death syndrome (SDS). Although the genome of F. virguliforme was sequenced, there were limited studies focused on the PCWDEs of F. virguliforme. Our goal was to understand the genomic CAZy structure of F. viguliforme, and determine if exogenous PIPs could be theoretically used in soybean to enhance resistance against F. virguliforme.

RESULTS

F. virguliforme produces diverse CAZy to degrade cellulose and pectin, similar to other necrotorphic and hemibiotrophic plant pathogenic fungi. However, some common CAZy of plant pathogenic fungi that catalyze hemicellulose, such as GH29, GH30, GH44, GH54, GH62, and GH67, were deficient in F. virguliforme. While the absence of these CAZy families might be complemented by other hemicellulases, F. virguliforme contained unique families including GH131, polysaccharide lyase (PL) 9, PL20, and PL22 that were not reported in other plant pathogenic fungi or oomycetes. Sequence analysis revealed two GH11 xylanases of F. virguliforme, FvXyn11A and FvXyn11B, have conserved residues that allow xylanase inhibitor protein I (XIP-I) binding. Structural modeling suggested that FvXyn11A and FvXyn11B could be blocked by XIP-I that serves as good candidate for developing transgenic soybeans. In contrast, one GH28 PG, FvPG2, contains an amino acid substitution that is potentially incompatible with the bean polygalacturonase-inhibitor protein II (PvPGIP2).

CONCLUSIONS

Identification and annotation of CAZy provided advanced understanding of genomic composition of PCWDEs in F. virguliforme. Sequence and structural analyses of FvXyn11A and FvXyn11B suggested both xylanases were conserved in residues that allow XIP-I inhibition, and expression of both xylanases were detected during soybean roots infection. We postulate that a transgenic soybean expressing wheat XIP-I may be useful for developing root rot resistance to F. virguliforme.

Genome-Wide Survey of Gut Fungi (Harpellales) Reveals the First Horizontally Transferred Ubiquitin Gene from a Mosquito Host

Harpellales, an early-diverging fungal lineage, is associated with the digestive tracts of aquatic arthropod hosts. Concurrent with the production and annotation of the first four Harpellales genomes, we discovered that Zancudomyces culisetae, one of the most widely distributed Harpellales species, encodes an insect-like polyubiquitin chain. Ubiquitin and ubiquitin-like proteins are universally involved in protein degradation and regulation of immune response in eukaryotic organisms. Phylogenetic analyses inferred that this polyubiquitin variant has a mosquito origin. In addition, its amino acid composition, animal-like secondary structure, as well as the fungal nature of flanking genes all further support this as a horizontal gene transfer event. The single-copy polyubiquitin gene from Z. culisetae has lower GC ratio compared with homologs of insect taxa, which implies homogenization of the gene since its putatively ancient transfer. The acquired polyubiquitin gene may have served to improve important functions within Z. culisetae, by perhaps exploiting the insect hosts' ubiquitin-proteasome systems in the gut environment. Preliminary comparisons among the four Harpellales genomes highlight the reduced genome size of Z. culisetae, which corroborates its distinguishable symbiotic lifestyle. This is the first record of a horizontally transferred ubiquitin gene from disease-bearing insects to the gut-dwelling fungal endobiont and should invite further exploration in an evolutionary context.

Quality Assessment of Domesticated Animal Genome Assemblies

The era of high-throughput sequencing has made it relatively simple to sequence genomes and transcriptomes of individuals from many species. In order to analyze the resulting sequencing data, high-quality reference genome assemblies are required. However, this is still a major challenge, and many domesticated animal genomes still need to be sequenced deeper in order to produce high-quality assemblies. In the meanwhile, ironically, the extent to which RNAseq and other next-generation data is produced frequently far exceeds that of the genomic sequence. Furthermore, basic comparative analysis is often affected by the lack of genomic sequence. Herein, we quantify the quality of the genome assemblies of 20 domesticated animals and related species by assessing a range of measurable parameters, and we show that there is a positive correlation between the fraction of mappable reads from RNAseq data and genome assembly quality. We rank the genomes by their assembly quality and discuss the implications for genotype analyses.

IMA Genome-F 6: Draft genome sequences of Armillaria fuscipes, Ceratocystiopsis minuta, Ceratocystis adiposa, Endoconidiophora laricicola, E. polonica and Penicillium freii DAOMC 242723

The genomes of Armillaria fuscipes, Ceratocystiopsis minuta, Ceratocystis adiposa, Endoconidiophora laricicola, E. polonica, and Penicillium freii DAOMC 242723 are presented in this genome announcement. These six genomes are from plant pathogens and otherwise economically important fungal species. The genome sizes range from 21 Mb in the case of Ceratocystiopsis minuta to 58 Mb for the basidiomycete Armillaria fuscipes. These genomes include the first reports of genomes for the genus Endoconidiophora. The availability of these genome data will provide opportunities to resolve longstanding questions regarding the taxonomy of species in these genera. In addition these genome sequences through comparative studies with closely related organisms will increase our understanding of how these pathogens cause disease.

GeneValidator: identify problems with protein-coding gene predictions

: Genomes of emerging model organisms are now being sequenced at very low cost. However, obtaining accurate gene predictions remains challenging: even the best gene prediction algorithms make substantial errors and can jeopardize subsequent analyses. Therefore, many predicted genes must be time-consumingly visually inspected and manually curated. We developed GeneValidator (GV) to automatically identify problematic gene predictions and to aid manual curation. For each gene, GV performs multiple analyses based on comparisons to gene sequences from large databases. The resulting report identifies problematic gene predictions and includes extensive statistics and graphs for each prediction to guide manual curation efforts. GV thus accelerates and enhances the work of biocurators and researchers who need accurate gene predictions from newly sequenced genomes.

AVAILABILITY AND IMPLEMENTATION

GV can be used through a web interface or in the command-line. GV is open-source (AGPL), available at https://wurmlab.github.io/tools/genevalidator

CONTACT

: y.wurm@qmul.ac.uk

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Crystal structure and molecular mechanism of an aspartate/glutamate racemase from Escherichia coli O157

EcL-DER, the aspartate/glutamate racemase from the pathogen Escherichia coli O157, exhibits racemase activity for l-aspartate and l-glutamate. This study reports the crystal structures of apo-EcL-DER, the EcL-DER-l-aspartate and the EcL-DER-d-aspartate complexes. The EcL-DER structure contains two domains, forming pseudo-mirror symmetry in the active site. A unique catalytic pair consisting of Thr(83) and Cys(197) exists in the active site. The characteristic conformations of l-Asp and d-Asp in the active site provide a straight structural evidence for the racemization mechanism of EcL-DER. In addition, the diversity of catalytic pairs implies that PLP-independent amino acid racemases adopt various catalytic mechanisms and are classified into different subgroups.

Functional elucidation of the non-coding RNAs of Kluyveromyces marxianus in the exponential growth phase

BACKGROUND

Non-coding RNAs (ncRNAs), which perform diverse regulatory roles, have been found in organisms from all superkingdoms of life. However, there have been limited numbers of studies on the functions of ncRNAs, especially in nonmodel organisms such as Kluyveromyces marxianus that is widely used in the field of industrial biotechnology.

RESULTS

In this study, we measured changes in transcriptome at three time points during the exponential growth phase of K. marxianus by using strand-specific RNA-seq. We found that approximately 60% of the transcriptome consists of ncRNAs transcribed from antisense and intergenic regions of the genome that were transcribed at lower levels than mRNA. In the transcriptome, a substantial number of long antisense ncRNAs (lancRNAs) are differentially expressed and enriched in carbohydrate and energy metabolism pathways. Furthermore, this enrichment is evolutionarily conserved, at least in yeast. Particularly, the mode of regulation of mRNA/lancRNA pairs is associated with mRNA transcription levels; the correlation between the pairs is positive at high mRNA transcriptional levels and negative at low levels. In addition, significant induction of mRNA and coverage of more than half of the mRNA sequence by a lancRNA strengthens the positive correlation between mRNA/lancRNA pairs.

CONCLUSIONS

Transcriptome sequencing of K. marxianus in the exponential growth phase reveals pervasive transcription of ncRNAs with evolutionarily conserved functions. Studies of the mode of regulation of mRNA/lancRNA pairs suggest that induction of lancRNA may be associated with switch-like behavior of mRNA/lancRNA pairs and efficient regulation of the carbohydrate and energy metabolism pathways in the exponential growth phase of K. marxianus being used in industrial applications.

Identification of Multiple Phytotoxins Produced by Fusarium virguliforme Including a Phytotoxic Effector (FvNIS1) Associated With Sudden Death Syndrome Foliar Symptoms

Sudden death syndrome (SDS) of soybean is caused by a soilborne pathogen, Fusarium virguliforme. Phytotoxins produced by F. virguliforme are translocated from infected roots to leaves, in which they cause SDS foliar symptoms. In this study, additional putative phytotoxins of F. virguliforme were identified, including three secondary metabolites and 11 effectors. While citrinin, fusaric acid, and radicicol induced foliar chlorosis and wilting, Soybean mosaic virus (SMV)-mediated overexpression of F. virguliforme necrosis-inducing secreted protein 1 (FvNIS1) induced SDS foliar symptoms that mimicked the development of foliar symptoms in the field. The expression level of fvnis1 remained steady over time, although foliar symptoms were delayed compared with the expression levels. SMV::FvNIS1 also displayed genotype-specific toxicity to which 75 of 80 soybean cultivars were susceptible. Genome-wide association mapping further identified three single nucleotide polymorphisms at two loci, where three leucine-rich repeat receptor-like protein kinase (LRR-RLK) genes were found. Culture filtrates of fvnis1 knockout mutants displayed a mild reduction in phytotoxicity, indicating that FvNIS1 is one of the phytotoxins responsible for SDS foliar symptoms and may contribute to the quantitative susceptibility of soybean by interacting with the LRR-RLK genes.

Optimizing and benchmarking de novo transcriptome sequencing: from library preparation to assembly evaluation

Background

RNA-seq enables gene expression profiling in selected spatiotemporal windows and yields massive sequence information with relatively low cost and time investment, even for non-model species. However, there remains a large room for optimizing its workflow, in order to take full advantage of continuously developing sequencing capacity.

Method

Transcriptome sequencing for three embryonic stages of Madagascar ground gecko (Paroedura picta) was performed with the Illumina platform. The output reads were assembled de novo for reconstructing transcript sequences. In order to evaluate the completeness of transcriptome assemblies, we prepared a reference gene set consisting of vertebrate one-to-one orthologs.

Result

To take advantage of increased read length of >150 nt, we demonstrated shortened RNA fragmentation time, which resulted in a dramatic shift of insert size distribution. To evaluate products of multiple de novo assembly runs incorporating reads with different RNA sources, read lengths, and insert sizes, we introduce a new reference gene set, core vertebrate genes (CVG), consisting of 233 genes that are shared as one-to-one orthologs by all vertebrate genomes examined (29 species)., The completeness assessment performed by the computational pipelines CEGMA and BUSCO referring to CVG, demonstrated higher accuracy and resolution than with the gene set previously established for this purpose. As a result of the assessment with CVG, we have derived the most comprehensive transcript sequence set of the Madagascar ground gecko by means of assembling individual libraries followed by clustering the assembled sequences based on their overall similarities.

Conclusion

Our results provide several insights into optimizing de novo RNA-seq workflow, including the coordination between library insert size and read length, which manifested in improved connectivity of assemblies. The approach and assembly assessment with CVG demonstrated here would be applicable to transcriptome analysis of other species as well as whole genome analyses.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2007-1) contains supplementary material, which is available to authorized users.

Draft genome sequence of the silver pomfret fish, Pampus argenteus

Silver pomfret, Pampus argenteus, is a fish species from coastal waters. Despite its high commercial value, this edible fish has not been sequenced. Hence, its genetic and genomic studies have been limited. We report the first draft genome sequence of the silver pomfret obtained using a Next Generation Sequencing (NGS) technology. We assembled 38.7 Gb of nucleotides into scaffolds of 350 Mb with N50 of about 1.5 kb, using high quality paired end reads. These scaffolds represent 63.7% of the estimated silver pomfret genome length. The newly sequenced and assembled genome has 11.06% repetitive DNA regions, and this percentage is comparable to that of the tilapia genome. The genome analysis predicted 16 322 genes. About 91% of these genes showed homology with known proteins. Many gene clusters were annotated to protein and fatty-acid metabolism pathways that may be important in the context of the meat texture and immune system developmental processes. The reference genome can pave the way for the identification of many other genomic features that could improve breeding and population-management strategies, and it can also help characterize the genetic diversity of P. argenteus.

Genome sequence of a laccase producing fungus Trametes sp. AH28-2

Trametes sp. AH28-2 (CCTCC AF 2015027) is a white rot fungus isolated from rotting wood in China. Primary study indicated that this strain can be induced by kraft lignin to secrete high levels of extracellular laccase, and differentially express laccase genes upon addition of different phenolic compounds. Here we report the complete genome sequence of Trametes sp. AH28-2 and its genetic basis for lignin degradation and phenolic xenobiotics metabolism.

NCBI BLAST+ integrated into Galaxy

Background

The NCBI BLAST suite has become ubiquitous in modern molecular biology and is used for small tasks such as checking capillary sequencing results of single PCR products, genome annotation or even larger scale pan-genome analyses. For early adopters of the Galaxy web-based biomedical data analysis platform, integrating BLAST into Galaxy was a natural step for sequence comparison workflows.

Findings

The command line NCBI BLAST+ tool suite was wrapped for use within Galaxy. Appropriate datatypes were defined as needed. The integration of the BLAST+ tool suite into Galaxy has the goal of making common BLAST tasks easy and advanced tasks possible.

Conclusions

This project is an informal international collaborative effort, and is deployed and used on Galaxy servers worldwide. Several examples of applications are described here.

Beyond the Whole-Genome Duplication: Phylogenetic Evidence for an Ancient Interspecies Hybridization in the Baker's Yeast Lineage

Whole-genome duplications have shaped the genomes of several vertebrate, plant, and fungal lineages. Earlier studies have focused on establishing when these events occurred and on elucidating their functional and evolutionary consequences, but we still lack sufficient understanding of how genome duplications first originated. We used phylogenomics to study the ancient genome duplication occurred in the yeast Saccharomyces cerevisiae lineage and found compelling evidence for the existence of a contemporaneous interspecies hybridization. We propose that the genome doubling was a direct consequence of this hybridization and that it served to provide stability to the recently formed allopolyploid. This scenario provides a mechanism for the origin of this ancient duplication and the lineage that originated from it and brings a new perspective to the interpretation of the origin and consequences of whole-genome duplications.

BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs

MOTIVATION

Genomics has revolutionized biological research, but quality assessment of the resulting assembled sequences is complicated and remains mostly limited to technical measures like N50.

RESULTS

We propose a measure for quantitative assessment of genome assembly and annotation completeness based on evolutionarily informed expectations of gene content. We implemented the assessment procedure in open-source software, with sets of Benchmarking Universal Single-Copy Orthologs, named BUSCO.

AVAILABILITY AND IMPLEMENTATION

Software implemented in Python and datasets available for download from http://busco.ezlab.org.

CONTACT

evgeny.zdobnov@unige.ch

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Transmission and Demographic Dynamics of Coxsackievirus B1

The infectious activity of coxsackievirus B1 (CV-B1) in Taiwan was high from 2008 to 2010, following an alarming increase in severe neonate disease in the United States (US). To examine the relationship between CV-B1 strains isolated in Taiwan and those from other parts of the world, we performed a phylodynamic study using VP1 and partial 3Dpol (414 nt) sequences from 22 strains of CV-B1 isolated in Taiwan (1989-2010) and compared them to sequences from strains isolated worldwide. Phylogenetic trees were constructed by neighbor-joining, maximum likelihood, and Bayesian Monte Carlo Markov Chain methods. Four genotypes (GI-IV) in the VP1 region of CV-B1 and three genotypes (GA-C) in the 3Dpol region of enterovirus B were identified and had high support values. The phylogenetic analysis indicates that the GI and GIII strains in VP1 were geographically distributed in Taiwan (1993-1994) and in India (2007-2009). On the other hand, the GII and GIV strains appear to have a wider spatiotemporal distribution and ladder-like topology A stair-like phylogeny was observed in the VP1 region indicating that the phylogeny of the virus may be affected by different selection pressures in the specified regions. Further, most of the GI and GII strains in the VP1 tree were clustered together in GA in the 3D tree, while the GIV strains diverged into GB and GC. Taken together, these data provide important insights into the population dynamics of CV-B1 and indicate that incongruencies in specific gene regions may contribute to spatiotemporal patterns of epidemicity for this virus.

Combining GWAS and RNA-Seq Approaches for Detection of the Causal Mutation for Hereditary Junctional Epidermolysis Bullosa in Sheep

In this study, we demonstrate the use of a genome-wide association mapping together with RNA-seq in a reduced number of samples, as an efficient approach to detect the causal mutation for a Mendelian disease. Junctional epidermolysis bullosa is a recessive genodermatosis that manifests with neonatal mechanical fragility of the skin, blistering confined to the lamina lucida of the basement membrane and severe alteration of the hemidesmosomal junctions. In Spanish Churra sheep, junctional epidermolysis bullosa (JEB) has been detected in two commercial flocks. The JEB locus was mapped to Ovis aries chromosome 11 by GWAS and subsequently fine-mapped to an 868-kb homozygous segment using the identical-by-descent method. The ITGB4, which is located within this region, was identified as the best positional and functional candidate gene. The RNA-seq variant analysis enabled us to discover a 4-bp deletion within exon 33 of the ITGB4 gene (c.4412_4415del). The c.4412_4415del mutation causes a frameshift resulting in a premature stop codon at position 1472 of the integrin β4 protein. A functional analysis of this deletion revealed decreased levels of mRNA in JEB skin samples and the absence of integrin β4 labeling in immunohistochemical assays. Genotyping of c.4412_4415del showed perfect concordance with the recessive mode of the disease phenotype. Selection against this causal mutation will now be used to solve the problem of JEB in flocks of Churra sheep. Furthermore, the identification of the ITGB4 mutation means that affected sheep can be used as a large mammal animal model for the human form of epidermolysis bullosa with aplasia cutis. Our approach evidences that RNA-seq offers cost-effective alternative to identify variants in the species in which high resolution exome-sequencing is not straightforward.

MicroRNA-like small RNAs prediction in the development of Antrodia cinnamomea

Antrodia cinnamomea, a precious, host-specific brown-rot fungus that has been used as a folk medicine in Taiwan for centuries is known to have diverse bioactive compounds with potent pharmaceutical activity. In this study, different fermentation states of A. cinnamomea (wild-type fruiting bodies and liquid cultured mycelium) were sequenced using the next-generation sequencing (NGS) technique. A 45.58 Mb genome encoding 6,522 predicted genes was obtained. High quality reads were assembled into a total of 13,109 unigenes. Using a previously constructed pipeline to search for microRNAs (miRNAs), we then identified 4 predicted conserved miRNA and 63 novel predicted miRNA-like small RNA (milRNA) candidates. Target prediction revealed several interesting proteins involved in tri-terpenoid synthesis, mating type recognition, chemical or physical sensory protein and transporters predicted to be regulated by the miRNAs and milRNAs.

Current methods for automated annotation of protein-coding genes

We review software tools for gene prediction - the identification of protein-coding genes and their structure in genome sequences. The discussed approaches include methods based on RNA-Seq and current methods based on homology - comparative gene prediction and protein spliced alignments. Many methods require that their parameters are adjusted to the target species or its broader clade. These include ab initio gene finders, integrated approaches with ab initio components and some aligners. We also review current automatic methods for training for the common case that a bona fide training set of gene structures is not available before annotation.

Decatransin, a new natural product inhibiting protein translocation at the Sec61/SecYEG translocon

A new cyclic decadepsipeptide was isolated from Chaetosphaeria tulasneorum with potent bioactivity on mammalian and yeast cells. Chemogenomic profiling in S. cerevisiae indicated that the Sec61 translocon complex, the machinery for protein translocation and membrane insertion at the endoplasmic reticulum, is the target. The profiles were similar to those of cyclic heptadepsipeptides of a distinct chemotype (including HUN-7293 and cotransin) that had previously been shown to inhibit cotranslational translocation at the mammalian Sec61 translocon. Unbiased, genome-wide mutagenesis followed by full-genome sequencing in both fungal and mammalian cells identified dominant mutations in Sec61p (yeast) or Sec61α1 (mammals) that conferred resistance. Most, but not all, of these mutations affected inhibition by both chemotypes, despite an absence of structural similarity. Biochemical analysis confirmed inhibition of protein translocation into the endoplasmic reticulum of both co- and post-translationally translocated substrates by both chemotypes, demonstrating a mechanism independent of a translating ribosome. Most interestingly, both chemotypes were found to also inhibit SecYEG, the bacterial Sec61 translocon homolog. We suggest 'decatransin' as the name for this new decadepsipeptide translocation inhibitor.

Development and bin mapping of gene-associated interspecific SNPs for cotton (Gossypium hirsutum L.) introgression breeding efforts

BACKGROUND

Cotton (Gossypium spp.) is the largest producer of natural fibers for textile and is an important crop worldwide. Crop production is comprised primarily of G. hirsutum L., an allotetraploid. However, elite cultivars express very small amounts of variation due to the species monophyletic origin, domestication and further bottlenecks due to selection. Conversely, wild cotton species harbor extensive genetic diversity of prospective utility to improve many beneficial agronomic traits, fiber characteristics, and resistance to disease and drought. Introgression of traits from wild species can provide a natural way to incorporate advantageous traits through breeding to generate higher-producing cotton cultivars and more sustainable production systems. Interspecific introgression efforts by conventional methods are very time-consuming and costly, but can be expedited using marker-assisted selection.

RESULTS

Using transcriptome sequencing we have developed the first gene-associated single nucleotide polymorphism (SNP) markers for wild cotton species G. tomentosum, G. mustelinum, G. armourianum and G. longicalyx. Markers were also developed for a secondary cultivated species G. barbadense cv. 3-79. A total of 62,832 non-redundant SNP markers were developed from the five wild species which can be utilized for interspecific germplasm introgression into cultivated G. hirsutum and are directly associated with genes. Over 500 of the G. barbadense markers have been validated by whole-genome radiation hybrid mapping. Overall 1,060 SNPs from the five different species have been screened and shown to produce acceptable genotyping assays.

CONCLUSIONS

This large set of 62,832 SNPs relative to cultivated G. hirsutum will allow for the first high-density mapping of genes from five wild species that affect traits of interest, including beneficial agronomic and fiber characteristics. Upon mapping, the markers can be utilized for marker-assisted introgression of new germplasm into cultivated cotton and in subsequent breeding of agronomically adapted types, including cultivar development.

Development and bin mapping of gene-associated interspecific SNPs for cotton (Gossypium hirsutum L.) introgression breeding efforts

Background

Cotton (Gossypium spp.) is the largest producer of natural fibers for textile and is an important crop worldwide. Crop production is comprised primarily of G. hirsutum L., an allotetraploid. However, elite cultivars express very small amounts of variation due to the species monophyletic origin, domestication and further bottlenecks due to selection. Conversely, wild cotton species harbor extensive genetic diversity of prospective utility to improve many beneficial agronomic traits, fiber characteristics, and resistance to disease and drought. Introgression of traits from wild species can provide a natural way to incorporate advantageous traits through breeding to generate higher-producing cotton cultivars and more sustainable production systems. Interspecific introgression efforts by conventional methods are very time-consuming and costly, but can be expedited using marker-assisted selection.

Results

Using transcriptome sequencing we have developed the first gene-associated single nucleotide polymorphism (SNP) markers for wild cotton species G. tomentosum, G. mustelinum, G. armourianum and G. longicalyx. Markers were also developed for a secondary cultivated species G. barbadense cv. 3–79. A total of 62,832 non-redundant SNP markers were developed from the five wild species which can be utilized for interspecific germplasm introgression into cultivated G. hirsutum and are directly associated with genes. Over 500 of the G. barbadense markers have been validated by whole-genome radiation hybrid mapping. Overall 1,060 SNPs from the five different species have been screened and shown to produce acceptable genotyping assays.

Conclusions

This large set of 62,832 SNPs relative to cultivated G. hirsutum will allow for the first high-density mapping of genes from five wild species that affect traits of interest, including beneficial agronomic and fiber characteristics. Upon mapping, the markers can be utilized for marker-assisted introgression of new germplasm into cultivated cotton and in subsequent breeding of agronomically adapted types, including cultivar development.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-945) contains supplementary material, which is available to authorized users.

Flexible and accessible workflows for improved proteogenomic analysis using the Galaxy framework

Proteogenomics combines large-scale genomic and transcriptomic data with mass-spectrometry-based proteomic data to discover novel protein sequence variants and improve genome annotation. In contrast with conventional proteomic applications, proteogenomic analysis requires a number of additional data processing steps. Ideally, these required steps would be integrated and automated via a single software platform offering accessibility for wet-bench researchers as well as flexibility for user-specific customization and integration of new software tools as they emerge. Toward this end, we have extended the Galaxy bioinformatics framework to facilitate proteogenomic analysis. Using analysis of whole human saliva as an example, we demonstrate Galaxy’s flexibility through the creation of a modular workflow incorporating both established and customized software tools that improve depth and quality of proteogenomic results. Our customized Galaxy-based software includes automated, batch-mode BLASTP searching and a Peptide Sequence Match Evaluator tool, both useful for evaluating the veracity of putative novel peptide identifications. Our complex workflow (approximately 140 steps) can be easily shared using built-in Galaxy functions, enabling their use and customization by others. Our results provide a blueprint for the establishment of the Galaxy framework as an ideal solution for the emerging field of proteogenomics.

Widespread genome reorganization of an obligate virus mutualist

The family Polydnaviridae is of interest because it provides the best example of viruses that have evolved a mutualistic association with their animal hosts. Polydnaviruses in the genus Bracovirus are strictly associated with parasitoid wasps in the family Braconidae, and evolved ∼100 million years ago from a nudivirus. Each wasp species relies on its associated bracovirus to parasitize hosts, while each bracovirus relies on its wasp for vertical transmission. Prior studies establish that bracovirus genomes consist of proviral segments and nudivirus-like replication genes, but how these components are organized in the genomes of wasps is unknown. Here, we sequenced the genome of the wasp Microplitis demolitor to characterize the proviral genome of M. demolitor bracovirus (MdBV). Unlike nudiviruses, bracoviruses produce virions that package multiple circular, double-stranded DNAs. DNA segments packaged into MdBV virions resided in eight dispersed loci in the M. demolitor genome. Each proviral segment was bounded by homologous motifs that guide processing to form mature viral DNAs. Rapid evolution of proviral segments obscured homology between other bracovirus-carrying wasps and MdBV. However, some domains flanking MdBV proviral loci were shared with other species. All MdBV genes previously identified to encode proteins required for replication were identified. Some of these genes resided in a multigene cluster but others, including subunits of the RNA polymerase that transcribes structural genes and integrases that process proviral segments, were widely dispersed in the M. demolitor genome. Overall, our results indicate that genome dispersal is a key feature in the evolution of bracoviruses into mutualists.

A draft genome for the African crocodilian trypanosome Trypanosoma grayi

The availability of genome sequence data has greatly enhanced our understanding of the adaptations of trypanosomatid parasites to their respective host environments. However, these studies remain somewhat restricted by modest taxon sampling, generally due to focus on the most important pathogens of humans. To address this problem, at least in part, we are releasing a draft genome sequence for the African crocodilian trypanosome, Trypanosoma grayi ANR4. This dataset comprises genomic DNA sequences assembled de novo into contigs, encompassing over 10,000 annotated putative open reading frames and predicted protein products. Using phylogenomic approaches we demonstrate that T. grayi is more closely related to Trypanosoma cruzi than it is to the African trypanosomes T. brucei, T. congolense and T. vivax, despite the fact T. grayi and the African trypanosomes are each transmitted by tsetse flies. The data are deposited in publicly accessible repositories where we hope they will prove useful to the community in evolutionary studies of the trypanosomatids.

Draft sequencing and analysis of the genome of pufferfish Takifugu flavidus

The pufferfish Takifugu flavidus is an important economic species due to its outstanding flavour and high market value. It has been regarded as an excellent model of genetic study for decades as well. In the present study, three mate-pair libraries of T. flavidus genome were sequenced by the SOLiD 4 next-generation sequencing platform, and the draft genome was constructed with the short reads using an assisted assembly strategy. The draft consists of 50,947 scaffolds with an N50 value of 305.7 kb, and the average GC content was 45.2%. The combined length of repetitive sequences was 26.5 Mb, which accounted for 6.87% of the genome, indicating that the compactness of T. flavidus genome was approximative with that of T. rubripes genome. A total of 1,253 non-coding RNA genes and 30,285 protein-encoding genes were assigned to the genome. There were 132,775 and 394 presumptive genes playing roles in the colour pattern variation, the relatively slow growth and the lipid metabolism, respectively. Among them, genes involved in the microtubule-dependent transport system, angiogenesis, decapentaplegic pathway and lipid mobilization were significantly expanded in the T. flavidus genome. This draft genome provides a valuable resource for understanding and improving both fundamental and applied research with pufferfish in the future.

SnowyOwl: accurate prediction of fungal genes by using RNA-Seq and homology information to select among ab initio models

BACKGROUND

Locating the protein-coding genes in novel genomes is essential to understanding and exploiting the genomic information but it is still difficult to accurately predict all the genes. The recent availability of detailed information about transcript structure from high-throughput sequencing of messenger RNA (RNA-Seq) delineates many expressed genes and promises increased accuracy in gene prediction. Computational gene predictors have been intensively developed for and tested in well-studied animal genomes. Hundreds of fungal genomes are now or will soon be sequenced. The differences of fungal genomes from animal genomes and the phylogenetic sparsity of well-studied fungi call for gene-prediction tools tailored to them.

RESULTS

SnowyOwl is a new gene prediction pipeline that uses RNA-Seq data to train and provide hints for the generation of Hidden Markov Model (HMM)-based gene predictions and to evaluate the resulting models. The pipeline has been developed and streamlined by comparing its predictions to manually curated gene models in three fungal genomes and validated against the high-quality gene annotation of Neurospora crassa; SnowyOwl predicted N. crassa genes with 83% sensitivity and 65% specificity. SnowyOwl gains sensitivity by repeatedly running the HMM gene predictor Augustus with varied input parameters and selectivity by choosing the models with best homology to known proteins and best agreement with the RNA-Seq data.

CONCLUSIONS

SnowyOwl efficiently uses RNA-Seq data to produce accurate gene models in both well-studied and novel fungal genomes. The source code for the SnowyOwl pipeline (in Python) and a web interface (in PHP) is freely available from http://sourceforge.net/projects/snowyowl/.

Building the sugarcane genome for biotechnology and identifying evolutionary trends

Background

Sugarcane is the source of sugar in all tropical and subtropical countries and is becoming increasingly important for bio-based fuels. However, its large (10 Gb), polyploid, complex genome has hindered genome based breeding efforts. Here we release the largest and most diverse set of sugarcane genome sequences to date, as part of an on-going initiative to provide a sugarcane genomic information resource, with the ultimate goal of producing a gold standard genome.

Results

Three hundred and seventeen chiefly euchromatic BACs were sequenced. A reference set of one thousand four hundred manually-annotated protein-coding genes was generated. A small RNA collection and a RNA-seq library were used to explore expression patterns and the sRNA landscape. In the sucrose and starch metabolism pathway, 16 non-redundant enzyme-encoding genes were identified. One of the sucrose pathway genes, sucrose-6-phosphate phosphohydrolase, is duplicated in sugarcane and sorghum, but not in rice and maize. A diversity analysis of the s6pp duplication region revealed haplotype-structured sequence composition. Examination of hom(e)ologous loci indicate both sequence structural and sRNA landscape variation. A synteny analysis shows that the sugarcane genome has expanded relative to the sorghum genome, largely due to the presence of transposable elements and uncharacterized intergenic and intronic sequences.

Conclusion

This release of sugarcane genomic sequences will advance our understanding of sugarcane genetics and contribute to the development of molecular tools for breeding purposes and gene discovery.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-540) contains supplementary material, which is available to authorized users.

Predicting the fungal CUG codon translation with Bagheera

Background

Many eukaryotes have been shown to use alternative schemes to the universal genetic code. While most Saccharomycetes, including Saccharomyces cerevisiae, use the standard genetic code translating the CUG codon as leucine, some yeasts, including many but not all of the “Candida”, translate the same codon as serine. It has been proposed that the change in codon identity was accomplished by an almost complete loss of the original CUG codons, making the CUG positions within the extant species highly discriminative for the one or other translation scheme.

Results

In order to improve the prediction of genes in yeast species by providing the correct CUG decoding scheme we implemented a web server, called Bagheera, that allows determining the most probable CUG codon translation for a given transcriptome or genome assembly based on extensive reference data. As reference data we use 2071 manually assembled and annotated sequences from 38 cytoskeletal and motor proteins belonging to 79 yeast species. The web service includes a pipeline, which starts with predicting and aligning homologous genes to the reference data. CUG codon positions within the predicted genes are analysed with respect to amino acid similarity and CUG codon conservation in related species. In addition, the tRNA_CAG gene is predicted in genomic data and compared to known leu-tRNA_CAG and ser-tRNA_CAG genes. Bagheera can also be used to evaluate any mRNA and protein sequence data with the codon usage of the respective species. The usage of the system has been demonstrated by analysing six genomes not included in the reference data.

Conclusions

Gene prediction and consecutive comparison with reference data from other Saccharomycetes are sufficient to predict the most probable decoding scheme for CUG codons. This approach has been implemented into Bagheera (http://www.motorprotein.de/bagheera).

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-411) contains supplementary material, which is available to authorized users.

Anaerobic gut fungi: Advances in isolation, culture, and cellulolytic enzyme discovery for biofuel production

Anaerobic gut fungi are an early branching family of fungi that are commonly found in the digestive tract of ruminants and monogastric herbivores. It is becoming increasingly clear that they are the primary colonizers of ingested plant biomass, and that they significantly contribute to the decomposition of plant biomass into fermentable sugars. As such, anaerobic fungi harbor a rich reservoir of undiscovered cellulolytic enzymes and enzyme complexes that can potentially transform the conversion of lignocellulose into bioenergy products. Despite their unique evolutionary history and cellulolytic activity, few species have been isolated and studied in great detail. As a result, their life cycle, cellular physiology, genetics, and cellulolytic metabolism remain poorly understood compared to aerobic fungi. To help address this limitation, this review briefly summarizes the current body of knowledge pertaining to anaerobic fungal biology, and describes progress made in the isolation, cultivation, molecular characterization, and long-term preservation of these microbes. We also discuss recent cellulase- and cellulosome-discovery efforts from gut fungi, and how these interesting, non-model microbes could be further adapted for biotechnology applications.