Approaches to Fungal Genome Annotation

Chromosome-level genome assembly and annotation of a potential model organism Gossypium arboreum ZB-1

Recent advancements in plant regeneration and synthetic polyploid creation have been documented in Gossypium arboreum ZB-1. These developments make ZB-1 a potential model within the Gossypium genus for investigating gene function and polyploidy. This work generated the sequence and annotation of the ZB-1 genome. The contig-level genome was constructed using the PacBio high-fidelity reads, encompassing 81 contigs with an N50 length of 112.12 Mb. The Hi-C data assisted the construction of the chromosome-level genome, which consists of 13 pseudo-chromosomes and 39 un-anchored contigs, with a total length of about 1.67 Gb. Repetitive sequences accounted for about 69.7% of the genome in length. Based on ab initio and evidence-based prediction, we have identified 48,021 protein-coding genes in the ZB-1 genome. Comparative genomics analysis revealed conserved gene content and arrangement between ZB-1 and G. arboreum SXY1. The single nucleotide polymorphism occurrence rate between ZB-1 and SXY1 was about 0.54 per 1,000 nucleotides. This study enriched the genomic resources for further exploration into cotton regeneration and polyploidy mechanisms.

Genome Analyses Reveal the Secondary Metabolites that Potentially Influence the Geographical Distribution of Fusarium pseudograminearum Populations

Fusarium crown rot (FCR), caused by Fusarium pseudograminearum, significantly impacts wheat yield and quality in China's Huanghuai region. The rapid F. pseudograminearum epidemic and FCR outbreak within a decade remain unexplained. In this study, two high-quality, chromosome-level genomes of F. pseudograminearum strains producing 3-acetyl-deoxynivalenol (3AcDON) and 15-acetyl-deoxynivalenol (15AcDON) toxins were assembled. Additionally, 38 related strains were resequenced. Genomic differences such as single nucleotide polymorphisms (SNPs), insertions/deletions (indels), and structural variations (SVs) among F. pseudograminearum strains were analyzed. The whole-genome SNP locus-based population classification mirrored the toxin chemotype (3AcDON and 15AcDON)-based classification, indicating the presence of genes associated with the trichothecene toxin gene cluster. Further analysis of differential SNP, indel, and SV loci between the 3AcDON and 15AcDON populations revealed a predominant connection to secondary metabolite synthesis genes. Notably, the majority of the secondary metabolite biosynthesis gene cluster loci were located in SNP-dense genomic regions, suggesting high mutability and a possible contribution to F. pseudograminearum population structure and environmental adaptability. This study provides insightful perspectives on the distribution and evolution of F. pseudograminearum and for forecasting the spread of wheat FCR, thereby aiding in the development of preventive measures and control strategies.

Diploid genome assembly of the Malbec grapevine cultivar enables haplotype-aware analysis of transcriptomic differences underlying clonal phenotypic variation

To preserve their varietal attributes, established grapevine cultivars (Vitis vinifera L. ssp. vinifera) must be clonally propagated, due to their highly heterozygous genomes. Malbec is a France-originated cultivar appreciated for producing high-quality wines and is the offspring of cultivars Prunelard and Magdeleine Noire des Charentes. Here, we have built a diploid genome assembly of Malbec, after trio binning of PacBio long reads into the two haploid complements inherited from either parent. After haplotype-aware deduplication and corrections, complete assemblies for the two haplophases were obtained with a very low haplotype switch-error rate (<0.025). The haplophase alignment identified > 25% of polymorphic regions. Gene annotation including RNA-seq transcriptome assembly and ab initio prediction evidence resulted in similar gene model numbers for both haplophases. The annotated diploid assembly was exploited in the transcriptomic comparison of four clonal accessions of Malbec that exhibited variation in berry composition traits. Analysis of the ripening pericarp transcriptome using either haplophases as a reference yielded similar results, although some differences were observed. Particularly, among the differentially expressed genes identified only with the Magdeleine-inherited haplotype as reference, we observed an over-representation of hypothetically hemizygous genes. The higher berry anthocyanin content of clonal accession 595 was associated with increased abscisic acid responses, possibly leading to the observed overexpression of phenylpropanoid metabolism genes and deregulation of genes associated with abiotic stress response. Overall, the results highlight the importance of producing diploid assemblies to fully represent the genomic diversity of highly heterozygous woody crop cultivars and unveil the molecular bases of clonal phenotypic variation.

Annotated genome and transcriptome of the endangered Caribbean mountainous star coral (Orbicella faveolata) using PacBio long-read sequencing

Long-read sequencing is revolutionizing de-novo genome assemblies, with continued advancements making it more readily available for previously understudied, non-model organisms. Stony corals are one such example, with long-read de-novo genome assemblies now starting to be publicly available, opening the door for a wide array of 'omics-based research. Here we present a new de-novo genome assembly for the endangered Caribbean star coral, Orbicella faveolata, using PacBio circular consensus reads. Our genome assembly improved the contiguity (51 versus 1,933 contigs) and complete and single copy BUSCO orthologs (93.6% versus 85.3%, database metazoa_odb10), compared to the currently available reference genome generated using short-read methodologies. Our new de-novo assembled genome also showed comparable quality metrics to other coral long-read genomes. Telomeric repeat analysis identified putative chromosomes in our scaffolded assembly, with these repeats at either one, or both ends, of scaffolded contigs. We identified 32,172 protein coding genes in our assembly through use of long-read RNA sequencing (ISO-seq) of additional O. faveolata fragments exposed to a range of abiotic and biotic treatments, and publicly available short-read RNA-seq data. With anthropogenic influences heavily affecting O. faveolata, as well as its increasing incorporation into reef restoration activities, this updated genome resource can be used for population genomics and other 'omics analyses to aid in the conservation of this species.

De novo assembly and annotation of the singing mouse genome

BACKGROUND

Developing genomic resources for a diverse range of species is an important step towards understanding the mechanisms underlying complex traits. Specifically, organisms that exhibit unique and accessible phenotypes-of-interest allow researchers to address questions that may be ill-suited to traditional model organisms. We sequenced the genome and transcriptome of Alston's singing mouse (Scotinomys teguina), an emerging model for social cognition and vocal communication. In addition to producing advertisement songs used for mate attraction and male-male competition, these rodents are diurnal, live at high-altitudes, and are obligate insectivores, providing opportunities to explore diverse physiological, ecological, and evolutionary questions.

RESULTS

Using PromethION, Illumina, and PacBio sequencing, we produced an annotated genome and transcriptome, which were validated using gene expression and functional enrichment analyses. To assess the usefulness of our assemblies, we performed single nuclei sequencing on cells of the orofacial motor cortex, a brain region implicated in song coordination, identifying 12 cell types.

CONCLUSIONS

These resources will provide the opportunity to identify the molecular basis of complex traits in singing mice as well as to contribute data that can be used for large-scale comparative analyses.

Clonal reproduction of Moniliophthora roreri and the emergence of unique lineages with distinct genomes during range expansion

The basidiomycete Moniliophthora roreri causes frosty pod rot of cacao (Theobroma cacao) in the western hemisphere. Moniliophthora roreri is considered asexual and haploid throughout its hemibiotrophic life cycle. To understand the processes driving genome modification, using long-read sequencing technology, we sequenced and assembled 5 high-quality M. roreri genomes out of a collection of 99 isolates collected throughout the pathogen's range. We obtained chromosome-scale assemblies composed of 11 scaffolds. We used short-read technology to sequence the genomes of 22 similarly chosen isolates. Alignments among the 5 reference assemblies revealed inversions, translocations, and duplications between and within scaffolds. Isolates at the front of the pathogens' expanding range tend to share lineage-specific structural variants, as confirmed by short-read sequencing. We identified, for the first time, 3 new mating type A locus alleles (5 in total) and 1 new potential mating type B locus allele (3 in total). Currently, only 2 mating type combinations, A1B1 and A2B2, are known to exist outside of Colombia. A systematic survey of the M. roreri transcriptome across 2 isolates identified an expanded candidate effector pool and provided evidence that effector candidate genes unique to the Moniliophthoras are preferentially expressed during the biotrophic phase of disease. Notably, M. roreri isolates in Costa Rica carry a chromosome segment duplication that has doubled the associated gene complement and includes secreted proteins and candidate effectors. Clonal reproduction of the haploid M. roreri genome has allowed lineages with unique genome structures and compositions to dominate as it expands its range, displaying a significant founder effect.

Convergent and complementary selection shaped gains and losses of eusociality in sweat bees

Sweat bees have repeatedly gained and lost eusociality, a transition from individual to group reproduction. Here we generate chromosome-length genome assemblies for 17 species and identify genomic signatures of evolutionary trade-offs associated with transitions between social and solitary living. Both young genes and regulatory regions show enrichment for these molecular patterns. We also identify loci that show evidence of complementary signals of positive and relaxed selection linked specifically to the convergent gains and losses of eusociality in sweat bees. This includes two pleiotropic proteins that bind and transport juvenile hormone (JH)-a key regulator of insect development and reproduction. We find that one of these proteins is primarily expressed in subperineurial glial cells that form the insect blood-brain barrier and that brain levels of JH vary by sociality. Our findings are consistent with a role of JH in modulating social behaviour and suggest that eusocial evolution was facilitated by alteration of the proteins that bind and transport JH, revealing how an ancestral developmental hormone may have been co-opted during one of life's major transitions. More broadly, our results highlight how evolutionary trade-offs have structured the molecular basis of eusociality in these bees and demonstrate how both directional selection and release from constraint can shape trait evolution.

Whole genome sequencing of Wilsonomyces carpophilus, an incitant of shot hole disease in stone fruits: insights into secreted proteins of a necrotrophic fungal repository

BACKGROUND

Shot hole is one of the important fungal diseases in stone fruits viz., peach, plum, apricot and cherry caused by Wilsonomyces carpophilus and almond among nut crops. Fungicides significantly decrease the disease. Pathogenicity studies proved a wide host range of the pathogen infecting all stone fruits and almond among the nut crops, however, the mechanism underlying host-pathogen interaction is still unknown. Molecular detection of the pathogen using polymerase chain reaction (PCR) based simple sequence repeat (SSR) markers is also unknown due to the unavailability of the pathogen genome.

METHODS AND RESULTS

We examined the morphology, pathology and genomics of the Wilsonomyces carpophilus. Whole genome sequencing of the W. carpophilus was carried out by Illumina HiSeq and PacBio high throughput sequencing plate-forms through hybrid assembly. Constant selection pressure alters the molecular mechanism of the pathogen causing disease. The studies revealed that the necrotrophs are more lethal with a complex pathogenicity mechanism and little-understood effector repositories. The different isolates of necrotrophic fungus W. carpophilus causing shot hole in stone fruits namely peach, plum, apricot and cherry, and almonds among the nut crops showed a significant variation in their morphology, however, the probability value (p = 0.29) suggests in-significant difference in the pathogenicity. Here, we reported draft genome of W. carpophilus of size 29.9 Mb (Accession number: PRJNA791904). A total of 10,901 protein-coding genes were predicted, including heterokaryon incompatibility genes, cytochrome-p450 genes, kinases, sugar transporters among others. We found 2851 simple sequence repeats (SSRs), tRNAs, rRNAs and pseudogenes in the genome. The most prominent proteins showing necrotrophic lifestyle of the pathogen were hydrolases, polysaccharide-degrading enzymes, esterolytic, lipolytic, and proteolytic enzymes accounted for 225 released proteins. Among the 223 fungal species, top-hit species distribution revealed the majority of hits against the Pyrenochaeta species followed by Ascochyta rabiei and Alternaria alternata.

CONCLUSION

Draft genome of W. carpophilus is 29.9 Mb based on Illumina HiSeq and PacBio hybrid assembly. The necrotrophs are more lethal with a complex pathogenicity mechanism. A significant variation in morphology was observed in different pathogen isolates. A total of 10,901 protein-coding genes were predicted in the pathogen genome including heterokaryon incompatibility, cytochrome-p450 genes, kinases and sugar transporters. We found 2851 SSRs, tRNAs, rRNAs and pseudogenes, and prominent proteins showing necrotrophic lifestyle such as hydrolases, polysaccharide-degrading enzymes, esterolytic, lipolytic and proteolytic enzymes. The top-hit species distribution were against the Pyrenochaeta spp. followed by Ascochyta rabiei.

Infection Process and Genome Assembly Provide Insights into the Pathogenic Mechanism of Destructive Mycoparasite Calcarisporium cordycipiticola with Host Specificity

Calcarisporium cordycipiticola is the pathogen in the white mildew disease of Cordyceps militaris, one of the popular mushrooms. This disease frequently occurs and there is no effective method for disease prevention and control. In the present study, C. militaris is found to be the only host of C. cordycipiticola, indicating strict host specificity. The infection process was monitored by fluorescent labeling and scanning and transmission electron microscopes. C. cordycipiticola can invade into the gaps among hyphae of the fruiting bodies of the host and fill them gradually. It can degrade the hyphae of the host by both direct contact and noncontact. The parasitism is initially biotrophic, and then necrotrophic as mycoparasitic interaction progresses. The approximate chromosome-level genome assembly of C. cordycipiticola yielded an N50 length of 5.45 Mbp and a total size of 34.51 Mbp, encoding 10,443 proteins. Phylogenomic analysis revealed that C. cordycipiticola is phylogenetically close to its specific host, C. militaris. A comparative genomic analysis showed that the number of CAZymes of C. cordycipiticola was much less than in other mycoparasites, which might be attributed to its host specificity. Secondary metabolite cluster analysis disclosed the great biosynthetic capabilities and potential mycotoxin production capability. This study provides insights into the potential pathogenesis and interaction between mycoparasite and its host.

High-quality reannotation of the king scallop genome reveals no 'gene-rich' feature and evolution of toxin resistance

The king scallop, Pecten maximus is a well-known, commercially important scallop species and is featured with remarkable tolerance to potent phytotoxins such as domoic acid. A high-quality genome can shed light on its biology and innovative evolution of toxin resistance. A reference genome has recently been published for P. maximus, however, it is suspicious that over 67,700 genes are annotated in this genome, which is unexpectedly larger than its close relatives of pectinids. Herein, we provide an improved high-quality chromosome-level reference genome assembly and annotation for the king scallop P. maximus. A final set of 26,995 genes is annotated after carefully checking and curation of the predicted gene models, which significantly improves the accuracy of gene structure information. The large number of gene duplicates in the previous genome is mainly distorted by the fragmented annotation. Through integrated genomic, evolutionary and transcriptomic analyses, we reveal that the Phi subfamily of ionotropic glutamate receptors (iGluRs) are well preserved in molluscs, and P. maximus experienced the rapid expansion of the Phi class of iGluR (GluF) gene family. The GluF genes exhibit ubiquitously high expression and altered sequence characteristics for ligand selectivity, which may contribute to the remarkable tolerance to neurotoxins in P. maximus. Taken together, our study disapproves the previous claim of the 'gene-rich' genome of this species and provides a high-quality genome assembly for further understanding of its biology and evolution of toxin resistance.

Comparative Transcriptome and Endophytic Bacterial Community Analysis of Morchella conica SH

The precious rare edible fungus Morchella conica is popular worldwide for its rich nutrition, savory flavor, and varieties of bioactive components. Due to its high commercial, nutritional, and medicinal value, it has always been a hot spot. However, the molecular mechanism and endophytic bacterial communities in M. conica were poorly understood. In this study, we sequenced, assembled, and analyzed the genome of M. conica SH. Transcriptome analysis reveals significant differences between the mycelia and fruiting body. As shown in this study, 1,329 and 2,796 genes were specifically expressed in the mycelia and fruiting body, respectively. The Gene Ontology (GO) enrichment showed that RNA polymerase II transcription activity-related genes were enriched in the mycelium-specific gene cluster, and nucleotide binding-related genes were enriched in the fruiting body-specific gene cluster. Further analysis of differentially expressed genes in different development stages resulted in finding two groups with distinct expression patterns. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment displays that glycan degradation and ABC transporters were enriched in the group 1 with low expressed level in the mycelia, while taurine and hypotaurine metabolismand tyrosine metabolism-related genes were significantly enriched in the group 2 with high expressed level in mycelia. Moreover, a dynamic shift of bacterial communities in the developing fruiting body was detected by 16S rRNA sequencing, and co-expression analysis suggested that bacterial communities might play an important role in regulating gene expression. Taken together, our study provided a better understanding of the molecular biology of M. conica SH and direction for future research on artificial cultivation.

Physiological, Genomic and Transcriptomic Analyses Reveal the Adaptation Mechanisms of Acidiella bohemica to Extreme Acid Mine Drainage Environments

Fungi in acid mine drainage (AMD) environments are of great concern due to their potentials of decomposing organic carbon, absorbing heavy metals and reducing AMD acidity. Based on morphological analysis and ITS/18S high-throughput sequencing technology, previous studies have provided deep insights into the diversity and community composition of fungi in AMD environments. However, knowledge about physiology, metabolic potential and transcriptome profiles of fungi inhabiting AMD environments is still scarce. Here, we reported the physiological, genomic, and transcriptomic characterization of Acidiella bohemica SYSU C17045 to improve our understanding of the physiological, genomic, and transcriptomic mechanisms underlying fungal adaptation to AMD environments. A. bohemica was isolated from an AMD environment, which has been proved to be an acidophilic fungus in this study. The surface of A. bohemica cultured in AMD solutions was covered with a large number of minerals such as jarosite. We thus inferred that the A. bohemica might have the potential of biologically induced mineralization. Taking advantage of PacBio single-molecule real-time sequencing, we obtained the high-quality genome sequences of A. bohemica (50 Mbp). To our knowledge, this was the first attempt to employ a third-generation sequencing technology to explore the genomic traits of fungi isolated from AMD environments. Moreover, our transcriptomic analysis revealed that a series of genes in the A. bohemica genome were related to its metabolic pathways of C, N, S, and Fe as well as its adaptation mechanisms, including the response to acid stress and the resistance to heavy metals. Overall, our physiological, genomic, and transcriptomic data provide a foundation for understanding the metabolic potential and adaptation mechanisms of fungi in AMD environments.

Application of an optimized annotation pipeline to the Cryptococcus deuterogattii genome reveals dynamic primary metabolic gene clusters and genomic impact of RNAi loss

Evaluating the quality of a de novo annotation of a complex fungal genome based on RNA-seq data remains a challenge. In this study, we sequentially optimized a Cufflinks-CodingQuary-based bioinformatics pipeline fed with RNA-seq data using the manually annotated model pathogenic yeasts Cryptococcus neoformans and Cryptococcus deneoformans as test cases. Our results show that the quality of the annotation is sensitive to the quantity of RNA-seq data used and that the best quality is obtained with 5–10 million reads per RNA-seq replicate. We also showed that the number of introns predicted is an excellent a priori indicator of the quality of the final de novo annotation. We then used this pipeline to annotate the genome of the RNAi-deficient species Cryptococcus deuterogattii strain R265 using RNA-seq data. Dynamic transcriptome analysis revealed that intron retention is more prominent in C. deuterogattii than in the other RNAi-proficient species C. neoformans and C. deneoformans. In contrast, we observed that antisense transcription was not higher in C. deuterogattii than in the two other Cryptococcus species. Comparative gene content analysis identified 21 clusters enriched in transcription factors and transporters that have been lost. Interestingly, analysis of the subtelomeric regions in these three annotated species identified a similar gene enrichment, reminiscent of the structure of primary metabolic clusters. Our data suggest that there is active exchange between subtelomeric regions, and that other chromosomal regions might participate in adaptive diversification of Cryptococcus metabolite assimilation potential.

Transcriptome Analysis of Tetraploid and Octoploid Common Reed (Phragmites australis)

Polyploidization in plants is thought to have occurred as coping mechanism with environmental stresses. Polyploidization-driven adaptation is often achieved through interplay of gene networks involved in differentially expressed genes, which triggers the plant to evolve special phenotypic traits for survival. Phragmites australis is a cosmopolitan species with highly variable phenotypic traits and high adaptation capacity to various habitats. The species' ploidy level varies from 3x to 12x, thus it is an ideal organism to investigate the molecular evolution of polyploidy and gene regulation mediated by different numbers of chromosome copies. In this study, we used high-throughput RNAseq data as a tool, to analyze the gene expression profiles in tetraploid and octoploid P. australis. The estimated divergence time between tetraploid and octoploid P. australis was dated to the border between Pliocene and Pleistocene. This study identified 439 up- and 956 down-regulated transcripts in tetraploids compared to octoploids. Gene ontology and pathway analysis revealed that tetraploids tended to express genes responsible for reproduction and seed germination to complete the reproduction cycle early, and expressed genes related to defense against UV-B light and fungi, whereas octoploids expressed mainly genes related to thermotolerance. Most differentially expressed genes were enriched in chaperones, folding catalysts and protein processing in endoplasmic reticulum pathways. Multiple biased isoform usage of the same gene was detected in differentially expressed genes, and the ones upregulated in octoploids were related to reduced DNA methylation. Our study provides new insights into the role of polyploidization on environmental responses and potential stress tolerance in grass species.

Genome annotation of disease-causing microorganisms

Humans have coexisted with pathogenic microorganisms throughout its history of evolution. We have never halted the exploration of pathogenic microorganisms. With the improvement of genome-sequencing technology and the continuous reduction of sequencing costs, an increasing number of complete genome sequences of pathogenic microorganisms have become available. Genome annotation of this massive sequence information has become a daunting task in biological research. This paper summarizes the approaches to the genome annotation of pathogenic microorganisms and the available popular genome annotation tools for prokaryotes, eukaryotes and viruses. Furthermore, real-world comparisons of different annotation tools using 12 genomes from prokaryotes, eukaryotes and viruses were conducted. Current challenges and problems were also discussed.

A comprehensive rat transcriptome built from large scale RNA-seq-based annotation

The rat is an important model organism in biomedical research for studying human disease mechanisms and treatments, but its annotated transcriptome is far from complete. We constructed a Rat Transcriptome Re-annotation named RTR using RNA-seq data from 320 samples in 11 different organs generated by the SEQC consortium. Totally, there are 52 807 genes and 114 152 transcripts in RTR. Transcribed regions and exons in RTR account for ∼42% and ∼6.5% of the genome, respectively. Of all 73 074 newly annotated transcripts in RTR, 34 213 were annotated as high confident coding transcripts and 24 728 as high confident long noncoding transcripts. Different tissues rather than different stages have a significant influence on the expression patterns of transcripts. We also found that 11 715 genes and 15 852 transcripts were expressed in all 11 tissues and that 849 house-keeping genes expressed different isoforms among tissues. This comprehensive transcriptome is freely available at http://www.unimd.org/rtr/. Our new rat transcriptome provides essential reference for genetics and gene expression studies in rat disease and toxicity models.

Quantitative global studies reveal differential translational control by start codon context across the fungal kingdom

Eukaryotic protein synthesis generally initiates at a start codon defined by an AUG and its surrounding Kozak sequence context, but the quantitative importance of this context in different species is unclear. We tested this concept in two pathogenic Cryptococcus yeast species by genome-wide mapping of translation and of mRNA 5' and 3' ends. We observed thousands of AUG-initiated upstream open reading frames (uORFs) that are a major contributor to translation repression. uORF use depends on the Kozak sequence context of its start codon, and uORFs with strong contexts promote nonsense-mediated mRNA decay. Transcript leaders in Cryptococcus and other fungi are substantially longer and more AUG-dense than in Saccharomyces. Numerous Cryptococcus mRNAs encode predicted dual-localized proteins, including many aminoacyl-tRNA synthetases, in which a leaky AUG start codon is followed by a strong Kozak context in-frame AUG, separated by mitochondrial-targeting sequence. Analysis of other fungal species shows that such dual-localization is also predicted to be common in the ascomycete mould, Neurospora crassa. Kozak-controlled regulation is correlated with insertions in translational initiation factors in fidelity-determining regions that contact the initiator tRNA. Thus, start codon context is a signal that quantitatively programs both the expression and the structures of proteins in diverse fungi.

The draft nuclear genome sequence and predicted mitochondrial proteome of Andalucia godoyi, a protist with the most gene-rich and bacteria-like mitochondrial genome

Background

Comparative analyses have indicated that the mitochondrion of the last eukaryotic common ancestor likely possessed all the key core structures and functions that are widely conserved throughout the domain Eucarya. To date, such studies have largely focused on animals, fungi, and land plants (primarily multicellular eukaryotes); relatively few mitochondrial proteomes from protists (primarily unicellular eukaryotic microbes) have been examined. To gauge the full extent of mitochondrial structural and functional complexity and to identify potential evolutionary trends in mitochondrial proteomes, more comprehensive explorations of phylogenetically diverse mitochondrial proteomes are required. In this regard, a key group is the jakobids, a clade of protists belonging to the eukaryotic supergroup Discoba, distinguished by having the most gene-rich and most bacteria-like mitochondrial genomes discovered to date.

Results

In this study, we assembled the draft nuclear genome sequence for the jakobid Andalucia godoyi and used a comprehensive in silico approach to infer the nucleus-encoded portion of the mitochondrial proteome of this protist, identifying 864 candidate mitochondrial proteins. The A. godoyi mitochondrial proteome has a complexity that parallels that of other eukaryotes, while exhibiting an unusually large number of ancestral features that have been lost particularly in opisthokont (animal and fungal) mitochondria. Notably, we find no evidence that the A. godoyi nuclear genome has or had a gene encoding a single-subunit, T3/T7 bacteriophage-like RNA polymerase, which functions as the mitochondrial transcriptase in all eukaryotes except the jakobids.

Conclusions

As genome and mitochondrial proteome data have become more widely available, a strikingly punctuate phylogenetic distribution of different mitochondrial components has been revealed, emphasizing that the pathways of mitochondrial proteome evolution are likely complex and lineage-specific. Unraveling this complexity will require comprehensive comparative analyses of mitochondrial proteomes from a phylogenetically broad range of eukaryotes, especially protists. The systematic in silico approach described here offers a valuable adjunct to direct proteomic analysis (e.g., via mass spectrometry), particularly in cases where the latter approach is constrained by sample limitation or other practical considerations.

Mycoparasitism illuminated by genome and transcriptome sequencing of Coniothyrium minitans, an important biocontrol fungus of the plant pathogen Sclerotinia sclerotiorum

Coniothyrium minitans is a mycoparasite of the notorious plant pathogen Sclerotinia sclerotiorum. To further understand the parasitism of C. minitans, we assembled and analysed its genome and performed transcriptome analyses. The genome of C. minitans strain ZS-1 was assembled into 350 scaffolds and had a size of 39.8 Mb. A total of 11 437 predicted genes and proteins were annotated, and 30.8 % of the blast hits matched proteins encoded by another member of the Pleosporales, Paraphaeosphaeria sporulosa, a worldwide soilborne fungus with biocontrol ability. The transcriptome of strain ZS-1 during the early interaction with S. sclerotiorum at 0, 4 and 12 h was analysed. The detected expressed genes were involved in responses to host defenses, including cell-wall-degrading enzymes, transporters, secretory proteins and secondary metabolite productions. Seventeen differentially expressed genes (DEGs) of fungal cell-wall-degrading enzymes (FCWDs) were up-regulated during parasitism, with only one down-regulated. Most of the monocarboxylate transporter genes of the major facilitator superfamily and all the detected ABC transporters, especially the heavy metal transporters, were significantly up-regulated. Approximately 8 % of the 11 437 proteins in C. minitans were predicted to be secretory proteins with catalytic activity. In the molecular function category, hydrolase activity, peptidase activity and serine hydrolase activity were enriched. Most genes involved in serine hydrolase activity were significantly up-regulated. This genomic analysis and genome-wide expression study demonstrates that the mycoparasitism process of C. minitans is complex and a broad range of proteins are deployed by C. minitans to successfully invade its host. Our study provides insights into the mechanisms of the mycoparasitism between C. minitans and S. sclerotiorum and identifies potential secondary metabolites from C. minitans for application as a biocontrol agent.

Reviving the Transcriptome Studies: An Insight Into the Emergence of Single-Molecule Transcriptome Sequencing

Advances in transcriptomics have provided an exceptional opportunity to study functional implications of the genetic variability. Technologies such as RNA-Seq have emerged as state-of-the-art techniques for transcriptome analysis that take advantage of high-throughput next-generation sequencing. However, similar to their predecessors, these approaches continue to impose major challenges on full-length transcript structure identification, primarily due to inherent limitations of read length. With the development of single-molecule sequencing (SMS) from PacBio, a growing number of studies on the transcriptome of different organisms have been reported. SMS has emerged as advantageous for comprehensive genome annotation including identification of novel genes/isoforms, long non-coding RNAs and fusion transcripts. This approach can be used across a broad spectrum of species to better interpret the coding information of the genome, and facilitate the biological function study. We provide an overview of SMS platform and its diverse applications in various biological studies, and our perspective on the challenges associated with the transcriptome studies.

Pectinase secreted by psychrotolerant fungi: identification, molecular characterization and heterologous expression of a cold-active polygalacturonase from Tetracladium sp

Background

Pectinolytic enzymes, which are used in several industries, especially in the clarification process during wine and fruit juice production, represent approximately 10% of the global enzyme market. To prevent the proliferation of undesired microorganisms, to retain labile and volatile flavor compounds, and to save energy, the current trend is to perform this process at low temperatures. However, the commercially available pectinases are highly active at temperatures approximately 50 °C and poorly active at temperatures below 35 °C, which is the reason why there is a constant search for cold-active pectinases. In preliminary studies, pectinolytic activity was detected in cold-adapted yeasts and yeast-like microorganisms isolated from Antarctica. The aim of the present work was to characterize pectinases secreted by these microorganisms and to express the best candidate in Pichia pastoris.

Results

Degradation of pectin by extracellular protein extracellular extracts obtained from 12 yeast cultures were assayed in plates at 4 °C to 37 °C and pH from 5.4 to 7.0, obtaining positive results in samples obtained from Dioszegia sp., Phenoliferia glacialis and Tetracladium sp. An enzyme was purified from Tetracladium sp., analyzed by peptide mass fingerprinting and compared to genome and transcriptome data from the same microorganism. Thus, the encoding gene was identified corresponding to a polygalacturonase-encoding gene. The enzyme was expressed in Pichia pastoris, and the recombinant polygalacturonase displayed higher activity at 15 °C than a mesophilic counterpart.

Conclusions

Extracellular pectinase activity was found in three yeast and yeast-like microorganisms from which the highest activity was displayed by Tetracladium sp., and the enzyme was identified as a polygalacturonase. The recombinant polygalacturonase produced in P. pastoris showed high activity at 15 °C, representing an attractive candidate to be applied in clarification processes in the production of fermented beverages and fruit juices.

Electronic supplementary material

The online version of this article (10.1186/s12934-019-1092-2) contains supplementary material, which is available to authorized users.

Long-Read Annotation: Automated Eukaryotic Genome Annotation Based on Long-Read cDNA Sequencing

Single-molecule full-length complementary DNA (cDNA) sequencing can aid genome annotation by revealing transcript structure and alternative splice forms, yet current annotation pipelines do not incorporate such information. Here we present long-read annotation (LoReAn) software, an automated annotation pipeline utilizing short- and long-read cDNA sequencing, protein evidence, and ab initio prediction to generate accurate genome annotations. Based on annotations of two fungal genomes (Verticillium dahliae and Plicaturopsis crispa) and two plant genomes (Arabidopsis [Arabidopsis thaliana] and Oryza sativa), we show that LoReAn outperforms popular annotation pipelines by integrating single-molecule cDNA-sequencing data generated from either the Pacific Biosciences or MinION sequencing platforms, correctly predicting gene structure, and capturing genes missed by other annotation pipelines.

Pathways of Pathogenicity: Transcriptional Stages of Germination in the Fatal Fungal Pathogen Rhizopus delemar

Rhizopus delemar is an invasive fungal pathogen responsible for the frequently fatal disease mucormycosis. Germination, a crucial mechanism by which infectious spores of Rhizopus delemar cause disease, is a key developmental process that transforms the dormant spore state into a vegetative one. The molecular mechanisms that underpin this transformation may be key to controlling mucormycosis; however, the regulation of germination remains poorly understood. This study describes the phenotypic and transcriptional changes that take place over the course of germination. This process is characterized by four distinct stages: dormancy, isotropic swelling, germ tube emergence, and hyphal growth. Dormant spores are shown to be transcriptionally unique, expressing a subset of transcripts absent in later developmental stages. A large shift in the expression profile is prompted by the initiation of germination, with genes involved in respiration, chitin, cytoskeleton, and actin regulation appearing to be important for this transition. A period of transcriptional consistency can be seen throughout isotropic swelling, before the transcriptional landscape shifts again at the onset of hyphal growth. This study provides a greater understanding of the regulation of germination and highlights processes involved in transforming Rhizopus delemar from a single-cellular to multicellular organism.IMPORTANCE Germination is key to the growth of many organisms, including fungal spores. Mucormycete spores exist abundantly within the environment and germinate to form hyphae. These spores are capable of infecting immunocompromised individuals, causing the disease mucormycosis. Germination from spore to hyphae within patients leads to angioinvasion, tissue necrosis, and often fatal infections. This study advances our understanding of how spore germination occurs in the mucormycetes, identifying processes we may be able to inhibit to help prevent or treat mucormycosis.

Integrated proteomics, genomics, metabolomics approaches reveal oxalic acid as pathogenicity factor in Tilletia indica inciting Karnal bunt disease of wheat

Tilletia indica incites Karnal bunt (KB) disease in wheat. To date, no KB resistant wheat cultivar could be developed due to non-availability of potential biomarkers related to pathogenicity/virulence for screening of resistant wheat genotypes. The present study was carried out to compare the proteomes of T. indica highly (TiK) and low (TiP) virulent isolates. Twenty one protein spots consistently observed as up-regulated/differential in the TiK proteome were selected for identification by MALDI-TOF/TOF. Identified sequences showed homology with fungal proteins playing essential role in plant infection and pathogen survival, including stress response, adhesion, fungal penetration, invasion, colonization, degradation of host cell wall, signal transduction pathway. These results were integrated with T. indica genome sequence for identification of homologs of candidate pathogenicity/virulence related proteins. Protein identified in TiK isolate as malate dehydrogenase that converts malate to oxaloacetate which is precursor of oxalic acid. Oxalic acid is key pathogenicity factor in phytopathogenic fungi. These results were validated by GC-MS based metabolic profiling of T. indica isolates indicating that oxalic acid was exclusively identified in TiK isolate. Thus, integrated omics approaches leads to identification of pathogenicity/virulence factor(s) that would provide insights into pathogenic mechanisms of fungi and aid in devising effective disease management strategies.

Whole genome and transcriptome analysis reveal adaptive strategies and pathogenesis of Calonectria pseudoreteaudii to Eucalyptus

BACKGROUND

Leaf blight caused by Calonectria spp. is one of the most destructive diseases to affect Eucalyptus nurseries and plantations. These pathogens mainly attack Eucalyptus, a tree with a diversity of secondary metabolites employed as defense-related phytoalexins. To unravel the fungal adaptive mechanisms to various phytoalexins, we examined the genome of C. pseudoreteaudii, which is one of the most aggressive pathogens in southeast Asia.

RESULTS

A 63.7 Mb genome with 14,355 coding genes of C. pseudoreteaudii were assembled. Genomic comparisons identified 1785 species-specific gene families in C. pseudoreteaudii. Most of them were not annotated and those annotated genes were enriched in peptidase activity, pathogenesis, oxidoreductase activity, etc. RNA-seq showed that 4425 genes were differentially expressed on the eucalyptus(the resistant cultivar E. grandis×E.camaldulensis M1) tissue induced medium. The annotation of GO term and KEGG pathway indicated that some of the differential expression genes were involved in detoxification and transportation, such as genes encoding ABC transporters, degrading enzymes of aromatic compounds and so on.

CONCLUSIONS

Potential genomic determinants of phytoalexin detoxification were identified in C. pseudoreteaudii by comparison with 13 other fungi. This pathogen seems to employ membrane transporters and degradation enzymes to detoxify Eucalyptus phytoalexins. Remarkably, the Calonectria genome possesses a surprising number of secondary metabolism backbone enzyme genes involving toxin biosynthesis. It is also especially suited for cutin and lignin degradation. This indicates that toxin and cell wall degrading enzymes may act important roles in the establishment of Calonectria leaf blight. This study provides further understanding on the mechanism of pathogenesis in Calonectria.

Whole-Genome Analysis Illustrates Global Clonal Population Structure of the Ubiquitous Dermatophyte Pathogen Trichophyton rubrum

Dermatophytes include fungal species that infect humans, as well as those that also infect other animals or only grow in the environment. The dermatophyte species Trichophyton rubrum is a frequent cause of skin infection in immunocompetent individuals. While members of the T. rubrum species complex have been further categorized based on various morphologies, their population structure and ability to undergo sexual reproduction are not well understood. In this study, we analyze a large set of T. rubrum and T. interdigitale isolates to examine mating types, evidence of mating, and genetic variation. We find that nearly all isolates of T. rubrum are of a single mating type, and that incubation with T. rubrum "morphotype" megninii isolates of the other mating type failed to induce sexual development. While the region around the mating type locus is characterized by a higher frequency of SNPs compared to other genomic regions, we find that the population is remarkably clonal, with highly conserved gene content, low levels of variation, and little evidence of recombination. These results support a model of recent transition to asexual growth when this species specialized to growth on human hosts.

Intraspecific comparative genomics of isolates of the Norway spruce pathogen (Heterobasidion parviporum) and identification of its potential virulence factors

BACKGROUND

Heterobasidion parviporum is an economically most important fungal forest pathogen in northern Europe, causing root and butt rot disease of Norway spruce (Picea abies (L.) Karst.). The mechanisms underlying the pathogenesis and virulence of this species remain elusive. No reference genome to facilitate functional analysis is available for this species.

RESULTS

To better understand the virulence factor at both phenotypic and genomic level, we characterized 15 H. parviporum isolates originating from different locations across Finland for virulence, vegetative growth, sporulation and saprotrophic wood decay. Wood decay capability and latitude of fungal origins exerted interactive effects on their virulence and appeared important for H. parviporum virulence. We sequenced the most virulent isolate, the first full genome sequences of H. parviporum as a reference genome, and re-sequenced the remaining 14 H. parviporum isolates. Genome-wide alignments and intrinsic polymorphism analysis showed that these isolates exhibited overall high genomic similarity with an average of at least 96% nucleotide identity when compared to the reference, yet had remarkable intra-specific level of polymorphism with a bias for CpG to TpG mutations. Reads mapping coverage analysis enabled the classification of all predicted genes into five groups and uncovered two genomic regions exclusively present in the reference with putative contribution to its higher virulence. Genes enriched for copy number variations (deletions and duplications) and nucleotide polymorphism were involved in oxidation-reduction processes and encoding domains relevant to transcription factors. Some secreted protein coding genes based on the genome-wide selection pressure, or the presence of variants were proposed as potential virulence candidates.

CONCLUSION

Our study reported on the first reference genome sequence for this Norway spruce pathogen (H. parviporum). Comparative genomics analysis gave insight into the overall genomic variation among this fungal species and also facilitated the identification of several secreted protein coding genes as putative virulence factors for the further functional analysis. We also analyzed and identified phenotypic traits potentially linked to its virulence.

Genome analysis reveals evolutionary mechanisms of adaptation in systemic dimorphic fungi

Dimorphic fungal pathogens cause a significant human disease burden and unlike most fungal pathogens affect immunocompetent hosts. To examine the origin of virulence of these fungal pathogens, we compared genomes of classic systemic, opportunistic, and non-pathogenic species, including Emmonsia and two basal branching, non-pathogenic species in the Ajellomycetaceae, Helicocarpus griseus and Polytolypa hystricis. We found that gene families related to plant degradation, secondary metabolites synthesis, and amino acid and lipid metabolism are retained in H. griseus and P. hystricis. While genes involved in the virulence of dimorphic pathogenic fungi are conserved in saprophytes, changes in the copy number of proteases, kinases and transcription factors in systemic dimorphic relative to non-dimorphic species may have aided the evolution of specialized gene regulatory programs to rapidly adapt to higher temperatures and new nutritional environments. Notably, both of the basal branching, non-pathogenic species appear homothallic, with both mating type locus idiomorphs fused at a single locus, whereas all related pathogenic species are heterothallic. These differences revealed that independent changes in nutrient acquisition capacity have occurred in the Onygenaceae and Ajellomycetaceae, and underlie how the dimorphic pathogens have adapted to the human host and decreased their capacity for growth in environmental niches.

Genome-Enhanced Detection and Identification (GEDI) of plant pathogens

Plant diseases caused by fungi and Oomycetes represent worldwide threats to crops and forest ecosystems. Effective prevention and appropriate management of emerging diseases rely on rapid detection and identification of the causal pathogens. The increase in genomic resources makes it possible to generate novel genome-enhanced DNA detection assays that can exploit whole genomes to discover candidate genes for pathogen detection. A pipeline was developed to identify genome regions that discriminate taxa or groups of taxa and can be converted into PCR assays. The modular pipeline is comprised of four components: (1) selection and genome sequencing of phylogenetically related taxa, (2) identification of clusters of orthologous genes, (3) elimination of false positives by filtering, and (4) assay design. This pipeline was applied to some of the most important plant pathogens across three broad taxonomic groups: Phytophthoras (Stramenopiles, Oomycota), Dothideomycetes (Fungi, Ascomycota) and Pucciniales (Fungi, Basidiomycota). Comparison of 73 fungal and Oomycete genomes led the discovery of 5,939 gene clusters that were unique to the targeted taxa and an additional 535 that were common at higher taxonomic levels. Approximately 28% of the 299 tested were converted into qPCR assays that met our set of specificity criteria. This work demonstrates that a genome-wide approach can efficiently identify multiple taxon-specific genome regions that can be converted into highly specific PCR assays. The possibility to easily obtain multiple alternative regions to design highly specific qPCR assays should be of great help in tackling challenging cases for which higher taxon-resolution is needed.

A Near-Complete Haplotype-Phased Genome of the Dikaryotic Wheat Stripe Rust Fungus Puccinia striiformis f. sp. tritici Reveals High Interhaplotype Diversity

A long-standing biological question is how evolution has shaped the genomic architecture of dikaryotic fungi. To answer this, high-quality genomic resources that enable haplotype comparisons are essential. Short-read genome assemblies for dikaryotic fungi are highly fragmented and lack haplotype-specific information due to the high heterozygosity and repeat content of these genomes. Here, we present a diploid-aware assembly of the wheat stripe rust fungus Puccinia striiformis f. sp. tritici based on long reads using the FALCON-Unzip assembler. Transcriptome sequencing data sets were used to infer high-quality gene models and identify virulence genes involved in plant infection referred to as effectors. This represents the most complete Puccinia striiformis f. sp. tritici genome assembly to date (83 Mb, 156 contigs, N50 of 1.5 Mb) and provides phased haplotype information for over 92% of the genome. Comparisons of the phase blocks revealed high interhaplotype diversity of over 6%. More than 25% of all genes lack a clear allelic counterpart. When we investigated genome features that potentially promote the rapid evolution of virulence, we found that candidate effector genes are spatially associated with conserved genes commonly found in basidiomycetes. Yet, candidate effectors that lack an allelic counterpart are more distant from conserved genes than allelic candidate effectors and are less likely to be evolutionarily conserved within the P. striiformis species complex and Pucciniales In summary, this haplotype-phased assembly enabled us to discover novel genome features of a dikaryotic plant-pathogenic fungus previously hidden in collapsed and fragmented genome assemblies.IMPORTANCE Current representations of eukaryotic microbial genomes are haploid, hiding the genomic diversity intrinsic to diploid and polyploid life forms. This hidden diversity contributes to the organism's evolutionary potential and ability to adapt to stress conditions. Yet, it is challenging to provide haplotype-specific information at a whole-genome level. Here, we take advantage of long-read DNA sequencing technology and a tailored-assembly algorithm to disentangle the two haploid genomes of a dikaryotic pathogenic wheat rust fungus. The two genomes display high levels of nucleotide and structural variations, which lead to allelic variation and the presence of genes lacking allelic counterparts. Nonallelic candidate effector genes, which likely encode important pathogenicity factors, display distinct genome localization patterns and are less likely to be evolutionary conserved than those which are present as allelic pairs. This genomic diversity may promote rapid host adaptation and/or be related to the age of the sequenced isolate since last meiosis.

Vegetative compatibility groups partition variation in the virulence of Verticillium dahliae on strawberry

Verticillium dahliae infection of strawberry (Fragaria x ananassa) is a major cause of disease-induced wilting in soil-grown strawberries across the world. To understand what components of the pathogen are affecting disease expression, the presence of the known effector VdAve1 was screened in a sample of Verticillium dahliae isolates. Isolates from strawberry were found to contain VdAve1 and were divided into two major clades, based upon their vegetative compatibility groups (VCG); no UK strawberry isolates contained VdAve1. VC clade was strongly related to their virulence levels. VdAve1-containing isolates pathogenic on strawberry were found in both clades, in contrast to some recently published findings. On strawberry, VdAve1-containing isolates had significantly higher virulence during early infection, which diminished in significance as the infection progressed. Transformation of a virulent non-VdAve1 containing isolate, with VdAve1 was found neither to increase nor decrease virulence when inoculated on a susceptible strawberry cultivar. There are therefore virulence factors that are epistatic to VdAve1 and potentially multiple independent routes to high virulence on strawberry in V. dahliae lineages. Genome sequencing a subset of isolates across the two VCGs revealed that isolates were differentiated at the whole genome level and contained multiple changes in putative effector content, indicating that different clonal VCGs may have evolved different strategies for infecting strawberry, leading to different virulence levels in pathogenicity tests. It is therefore important to consider both clonal lineage and effector complement as the adaptive potential of each lineage will differ, even if they contain the same race determining effector.

Genome sequencing and comparative genomics reveal the potential pathogenic mechanism of Cercospora sojina Hara on soybean

Frogeye leaf spot, caused by Cercospora sojina Hara, is a common disease of soybean in most soybean-growing countries of the world. In this study, we report a high-quality genome sequence of C. sojina by Single Molecule Real-Time sequencing method. The 40.8-Mb genome encodes 11,655 predicated genes, and 8,474 genes are revealed by RNA sequencing. Cercospora sojina genome contains large numbers of gene clusters that are involved in synthesis of secondary metabolites, including mycotoxins and pigments. However, much less carbohydrate-binding module protein encoding genes are identified in C. sojina genome, when compared with other phytopathogenic fungi. Bioinformatics analysis reveals that C. sojina harbours about 752 secreted proteins, and 233 of them are effectors. During early infection, the genes for metabolite biosynthesis and effectors are significantly enriched, suggesting that they may play essential roles in pathogenicity. We further identify 13 effectors that can inhibit BAX-induced cell death. Taken together, our results provide insights into the infection mechanisms of C. sojina on soybean.

The genome sequence of the soft-rot fungus Penicillium purpurogenum reveals a high gene dosage for lignocellulolytic enzymes

The high lignocellulolytic activity displayed by the soft-rot fungus Penicillium purpurogenum has made it a target for the study of novel lignocellulolytic enzymes. We have obtained a reference genome of 36.2 Mb of non-redundant sequence (11,057 protein-coding genes). The 49 largest scaffolds cover 90% of the assembly, and Core Eukaryotic Genes Mapping Approach (CEGMA) analysis reveals that our assembly captures almost all protein-coding genes. RNA-seq was performed and 93.1% of the reads aligned to the assembled genome. These data, plus the independent sequencing of a set of genes of lignocellulose-degrading enzymes, validate the quality of the genome sequence. P. purpurogenum shows a higher number of proteins with CAZy motifs, transcription factors and transporters as compared to other sequenced Penicillia. These results demonstrate the great potential for lignocellulolytic activity of this fungus and the possible use of its enzymes in related industrial applications.

Fungal Genome Annotation

The term "genome annotation" includes identification of protein-coding and noncoding sequences (e.g., repeats, rDNA, and ncRNA) in genome assemblies and attaching functional information (metadata) to these annotated features. Here, we describe the basic outline of fungal nuclear and mitochondrial genome annotation as performed at the US Department of Energy Joint Genome Institute (JGI).

Comparative ribosome profiling uncovers a dominant role for translational control in Toxoplasma gondii

Background

The lytic cycle of the protozoan parasite Toxoplasma gondii, which involves a brief sojourn in the extracellular space, is characterized by defined transcriptional profiles. For an obligate intracellular parasite that is shielded from the cytosolic host immune factors by a parasitophorous vacuole, the brief entry into the extracellular space is likely to exert enormous stress. Due to its role in cellular stress response, we hypothesize that translational control plays an important role in regulating gene expression in Toxoplasma during the lytic cycle. Unlike transcriptional profiles, insights into genome-wide translational profiles of Toxoplasma gondii are lacking.

Methods

We have performed genome-wide ribosome profiling, coupled with high throughput RNA sequencing, in intracellular and extracellular Toxoplasma gondii parasites to investigate translational control during the lytic cycle.

Results

Although differences in transcript abundance were mostly mirrored at the translational level, we observed significant differences in the abundance of ribosome footprints between the two parasite stages. Furthermore, our data suggest that mRNA translation in the parasite is potentially regulated by mRNA secondary structure and upstream open reading frames.

Conclusion

We show that most of the Toxoplasma genes that are dysregulated during the lytic cycle are translationally regulated.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4362-6) contains supplementary material, which is available to authorized users.

The mitochondrial genome of the plant-pathogenic fungus Stemphylium lycopersici uncovers a dynamic structure due to repetitive and mobile elements

Stemphylium lycopersici (Pleosporales) is a plant-pathogenic fungus that has been associated with a broad range of plant-hosts worldwide. It is one of the causative agents of gray leaf spot disease in tomato and pepper. The aim of this work was to characterize the mitochondrial genome of S. lycopersici CIDEFI-216, to use it to trace taxonomic relationships with other fungal taxa and to get insights into the evolutionary history of this phytopathogen. The complete mitochondrial genome was assembled into a circular double-stranded DNA molecule of 75,911 bp that harbors a set of 37 protein-coding genes, 2 rRNA genes (rns and rnl) and 28 tRNA genes, which are transcribed from both sense and antisense strands. Remarkably, its gene repertoire lacks both atp8 and atp9, contains a free-standing gene for the ribosomal protein S3 (rps3) and includes 13 genes with homing endonuclease domains that are mostly located within its 15 group I introns. Strikingly, subunits 1 and 2 of cytochrome oxidase are encoded by a single continuous open reading frame (ORF). A comparative mitogenomic analysis revealed the large extent of structural rearrangements among representatives of Pleosporales, showing the plasticity of their mitochondrial genomes. Finally, an exhaustive phylogenetic analysis of the subphylum Pezizomycotina based on mitochondrial data reconstructed their relationships in concordance with several studies based on nuclear data. This is the first report of a mitochondrial genome belonging to a representative of the family Pleosporaceae.

Rapid Expansion of Immune-Related Gene Families in the House Fly, Musca domestica

The house fly, Musca domestica, occupies an unusual diversity of potentially septic niches compared with other sequenced Dipteran insects and is a vector of numerous diseases of humans and livestock. In the present study, we apply whole-transcriptome sequencing to identify genes whose expression is regulated in adult flies upon bacterial infection. We then combine the transcriptomic data with analysis of rates of gene duplication and loss to provide insight into the evolutionary dynamics of immune-related genes. Genes up-regulated after bacterial infection are biased toward being evolutionarily recent innovations, suggesting the recruitment of novel immune components in the M. domestica or ancestral Dipteran lineages. In addition, using new models of gene family evolution, we show that several different classes of immune-related genes, particularly those involved in either pathogen recognition or pathogen killing, are duplicating at a significantly accelerated rate on the M. domestica lineage relative to other Dipterans. Taken together, these results suggest that the M. domestica immune response includes an elevated diversity of genes, perhaps as a consequence of its lifestyle in septic environments.

Harnessing Whole Genome Sequencing in Medical Mycology

Purpose of Review

Comparative genome sequencing studies of human fungal pathogens enable identification of genes and variants associated with virulence and drug resistance. This review describes current approaches, resources, and advances in applying whole genome sequencing to study clinically important fungal pathogens.

Recent Findings

Genomes for some important fungal pathogens were only recently assembled, revealing gene family expansions in many species and extreme gene loss in one obligate species. The scale and scope of species sequenced is rapidly expanding, leveraging technological advances to assemble and annotate genomes with higher precision. By using iteratively improved reference assemblies or those generated de novo for new species, recent studies have compared the sequence of isolates representing populations or clinical cohorts. Whole genome approaches provide the resolution necessary for comparison of closely related isolates, for example, in the analysis of outbreaks or sampled across time within a single host.

Summary

Genomic analysis of fungal pathogens has enabled both basic research and diagnostic studies. The increased scale of sequencing can be applied across populations, and new metagenomic methods allow direct analysis of complex samples.

Divergent and convergent modes of interaction between wheat and Puccinia graminis f. sp. tritici isolates revealed by the comparative gene co-expression network and genome analyses

BACKGROUND

Two opposing evolutionary constraints exert pressure on plant pathogens: one to diversify virulence factors in order to evade plant defenses, and the other to retain virulence factors critical for maintaining a compatible interaction with the plant host. To better understand how the diversified arsenals of fungal genes promote interaction with the same compatible wheat line, we performed a comparative genomic analysis of two North American isolates of Puccinia graminis f. sp. tritici (Pgt).

RESULTS

The patterns of inter-isolate divergence in the secreted candidate effector genes were compared with the levels of conservation and divergence of plant-pathogen gene co-expression networks (GCN) developed for each isolate. Comprative genomic analyses revealed substantial level of interisolate divergence in effector gene complement and sequence divergence. Gene Ontology (GO) analyses of the conserved and unique parts of the isolate-specific GCNs identified a number of conserved host pathways targeted by both isolates. Interestingly, the degree of inter-isolate sub-network conservation varied widely for the different host pathways and was positively associated with the proportion of conserved effector candidates associated with each sub-network. While different Pgt isolates tended to exploit similar wheat pathways for infection, the mode of plant-pathogen interaction varied for different pathways with some pathways being associated with the conserved set of effectors and others being linked with the diverged or isolate-specific effectors.

CONCLUSIONS

Our data suggest that at the intra-species level pathogen populations likely maintain divergent sets of effectors capable of targeting the same plant host pathways. This functional redundancy may play an important role in the dynamic of the "arms-race" between host and pathogen serving as the basis for diverse virulence strategies and creating conditions where mutations in certain effector groups will not have a major effect on the pathogen's ability to infect the host.

Comparative Analysis Highlights Variable Genome Content of Wheat Rusts and Divergence of the Mating Loci

Three members of the Puccinia genus, Pucciniatriticina (Pt), Pstriiformis f.sp. tritici (Pst), and Pgraminis f.sp. tritici (Pgt), cause the most common and often most significant foliar diseases of wheat. While similar in biology and life cycle, each species is uniquely adapted and specialized. The genomes of Pt and Pst were sequenced and compared to that of Pgt to identify common and distinguishing gene content, to determine gene variation among wheat rust pathogens, other rust fungi, and basidiomycetes, and to identify genes of significance for infection. Pt had the largest genome of the three, estimated at 135 Mb with expansion due to mobile elements and repeats encompassing 50.9% of contig bases; in comparison, repeats occupy 31.5% for Pst and 36.5% for Pgt We find all three genomes are highly heterozygous, with Pst [5.97 single nucleotide polymorphisms (SNPs)/kb] nearly twice the level detected in Pt (2.57 SNPs/kb) and that previously reported for Pgt Of 1358 predicted effectors in Pt, 784 were found expressed across diverse life cycle stages including the sexual stage. Comparison to related fungi highlighted the expansion of gene families involved in transcriptional regulation and nucleotide binding, protein modification, and carbohydrate degradation enzymes. Two allelic homeodomain pairs, HD1 and HD2, were identified in each dikaryotic Puccinia species along with three pheromone receptor (STE3) mating-type genes, two of which are likely representing allelic specificities. The HD proteins were active in a heterologous Ustilago maydis mating assay and host-induced gene silencing (HIGS) of the HD and STE3 alleles reduced wheat host infection.

Genome Sequence of Spizellomyces punctatus

Spizellomyces punctatus is a basally branching chytrid fungus that is found in the Chytridiomycota phylum. Spizellomyces species are common in soil and of importance in terrestrial ecosystems. Here, we report the genome sequence of S. punctatus, which will facilitate the study of this group of early diverging fungi.

Genome Sequence of the Edible Cultivated Mushroom Lentinula edodes (Shiitake) Reveals Insights into Lignocellulose Degradation

Lentinula edodes, one of the most popular, edible mushroom species with a high content of proteins and polysaccharides as well as unique aroma, is widely cultivated in many Asian countries, especially in China, Japan and Korea. As a white rot fungus with lignocellulose degradation ability, L. edodes has the potential for application in the utilization of agriculture straw resources. Here, we report its 41.8-Mb genome, encoding 14,889 predicted genes. Through a phylogenetic analysis with model species of fungi, the evolutionary divergence time of L. edodes and Gymnopus luxurians was estimated to be 39 MYA. The carbohydrate-active enzyme genes in L. edodes were compared with those of the other 25 fungal species, and 101 lignocellulolytic enzymes were identified in L. edodes, similar to other white rot fungi. Transcriptome analysis showed that the expression of genes encoding two cellulases and 16 transcription factor was up-regulated when mycelia were cultivated for 120 minutes in cellulose medium versus glucose medium. Our results will foster a better understanding of the molecular mechanism of lignocellulose degradation and provide the basis for partial replacement of wood sawdust with agricultural wastes in L. edodes cultivation.

Comparative genomics and transcriptomics of Pichia pastoris

Background

Pichia pastoris has emerged as an important alternative host for producing recombinant biopharmaceuticals, owing to its high cultivation density, low host cell protein burden, and the development of strains with humanized glycosylation. Despite its demonstrated utility, relatively little strain engineering has been performed to improve Pichia, due in part to the limited number and inconsistent frameworks of reported genomes and transcriptomes. Furthermore, the co-mingling of genomic, transcriptomic and fermentation data collected about Komagataella pastoris and Komagataella phaffii, the two strains co-branded as Pichia, has generated confusion about host performance for these genetically distinct species. Generation of comparative high-quality genomes and transcriptomes will enable meaningful comparisons between the organisms, and potentially inform distinct biotechnological utilies for each species.

Results

Here, we present a comprehensive and standardized comparative analysis of the genomic features of the three most commonly used strains comprising the tradename Pichia: K. pastoris wild-type, K. phaffii wild-type, and K. phaffii GS115. We used a combination of long-read (PacBio) and short-read (Illumina) sequencing technologies to achieve over 1000X coverage of each genome. Construction of individual genomes was then performed using as few as seven individual contigs to create gap-free assemblies. We found substantial syntenic rearrangements between the species and characterized a linear plasmid present in K. phaffii. Comparative analyses between K. phaffii genomes enabled the characterization of the mutational landscape of the GS115 strain. We identified and examined 35 non-synonomous coding mutations present in GS115, many of which are likely to impact strain performance. Additionally, we investigated transcriptomic profiles of gene expression for both species during cultivation on various carbon sources. We observed that the most highly transcribed genes in both organisms were consistently highly expressed in all three carbon sources examined. We also observed selective expression of certain genes in each carbon source, including many sequences not previously reported as promoters for expression of heterologous proteins in yeasts.

Conclusions

Our studies establish a foundation for understanding critical relationships between genome structure, cultivation conditions and gene expression. The resources we report here will inform and facilitate rational, organism-wide strain engineering for improved utility as a host for protein production.

Electronic supplementary material

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

Computational Identification of Novel Genes: Current and Future Perspectives

While it has long been thought that all genomic novelties are derived from the existing material, many genes lacking homology to known genes were found in recent genome projects. Some of these novel genes were proposed to have evolved de novo, ie, out of noncoding sequences, whereas some have been shown to follow a duplication and divergence process. Their discovery called for an extension of the historical hypotheses about gene origination. Besides the theoretical breakthrough, increasing evidence accumulated that novel genes play important roles in evolutionary processes, including adaptation and speciation events. Different techniques are available to identify genes and classify them as novel. Their classification as novel is usually based on their similarity to known genes, or lack thereof, detected by comparative genomics or against databases. Computational approaches are further prime methods that can be based on existing models or leveraging biological evidences from experiments. Identification of novel genes remains however a challenging task. With the constant software and technologies updates, no gold standard, and no available benchmark, evaluation and characterization of genomic novelty is a vibrant field. In this review, the classical and state-of-the-art tools for gene prediction are introduced. The current methods for novel gene detection are presented; the methodological strategies and their limits are discussed along with perspective approaches for further studies.

Daphnia magna transcriptome by RNA-Seq across 12 environmental stressors

The full exploration of gene-environment interactions requires model organisms with well-characterized ecological interactions in their natural environment, manipulability in the laboratory and genomic tools. The waterflea Daphnia magna is an established ecological and toxicological model species, central to the food webs of freshwater lentic habitats and sentinel for water quality. Its tractability and cyclic parthenogenetic life-cycle are ideal to investigate links between genes and the environment. Capitalizing on this unique model system, the STRESSFLEA consortium generated a comprehensive RNA-Seq data set by exposing two inbred genotypes of D. magna and a recombinant cross of these genotypes to a range of environmental perturbations. Gene models were constructed from the transcriptome data and mapped onto the draft genome of D. magna using EvidentialGene. The transcriptome data generated here, together with the available draft genome sequence of D. magna and a high-density genetic map will be a key asset for future investigations in environmental genomics.

DemaDb: an integrated dematiaceous fungal genomes database

Many species of dematiaceous fungi are associated with allergic reactions and potentially fatal diseases in human, especially in tropical climates. Over the past 10 years, we have isolated more than 400 dematiaceous fungi from various clinical samples. In this study, DemaDb, an integrated database was designed to support the integration and analysis of dematiaceous fungal genomes. A total of 92 072 putative genes and 6527 pathways that identified in eight dematiaceous fungi (Bipolaris papendorfii UM 226, Daldinia eschscholtzii UM 1400, D. eschscholtzii UM 1020, Pyrenochaeta unguis-hominis UM 256, Ochroconis mirabilis UM 578, Cladosporium sphaerospermum UM 843, Herpotrichiellaceae sp. UM 238 and Pleosporales sp. UM 1110) were deposited in DemaDb. DemaDb includes functional annotations for all predicted gene models in all genomes, such as Gene Ontology, EuKaryotic Orthologous Groups, Kyoto Encyclopedia of Genes and Genomes (KEGG), Pfam and InterProScan. All predicted protein models were further functionally annotated to Carbohydrate-Active enzymes, peptidases, secondary metabolites and virulence factors. DemaDb Genome Browser enables users to browse and visualize entire genomes with annotation data including gene prediction, structure, orientation and custom feature tracks. The Pathway Browser based on the KEGG pathway database allows users to look into molecular interaction and reaction networks for all KEGG annotated genes. The availability of downloadable files containing assembly, nucleic acid, as well as protein data allows the direct retrieval for further downstream works. DemaDb is a useful resource for fungal research community especially those involved in genome-scale analysis, functional genomics, genetics and disease studies of dematiaceous fungi. Database URL: http://fungaldb.um.edu.my.

Draft Genome Sequences of Fungus Aspergillus calidoustus

Here, we report the draft genome sequence of Aspergillus calidoustus (strain SF006504). The functional annotation of A. calidoustus predicts a relatively large number of secondary metabolite gene clusters. The presented genome sequence builds the basis for further genome mining.

Genome analysis of three Pneumocystis species reveals adaptation mechanisms to life exclusively in mammalian hosts

Pneumocystis jirovecii is a major cause of life-threatening pneumonia in immunosuppressed patients including transplant recipients and those with HIV/AIDS, yet surprisingly little is known about the biology of this fungal pathogen. Here we report near complete genome assemblies for three Pneumocystis species that infect humans, rats and mice. Pneumocystis genomes are highly compact relative to other fungi, with substantial reductions of ribosomal RNA genes, transporters, transcription factors and many metabolic pathways, but contain expansions of surface proteins, especially a unique and complex surface glycoprotein superfamily, as well as proteases and RNA processing proteins. Unexpectedly, the key fungal cell wall components chitin and outer chain N-mannans are absent, based on genome content and experimental validation. Our findings suggest that Pneumocystis has developed unique mechanisms of adaptation to life exclusively in mammalian hosts, including dependence on the lungs for gas and nutrients and highly efficient strategies to escape both host innate and acquired immune defenses.

Evolution of regulatory networks in Candida glabrata: learning to live with the human host

The opportunistic human fungal pathogen Candida glabrata is second only to C. albicans as the cause of Candida infections and yet is more closely related to Saccharomyces cerevisiae. Recent advances in functional genomics technologies and computational approaches to decipher regulatory networks, and the comparison of these networks among these and other Ascomycete species, have revealed both unique and shared strategies in adaptation to a human commensal/opportunistic pathogen lifestyle and antifungal drug resistance in C. glabrata. Recently, several C. glabrata sister species in the Nakeseomyces clade representing both human associated (commensal) and environmental isolates have had their genomes sequenced and analyzed. This has paved the way for comparative functional genomics studies to characterize the regulatory networks in these species to identify informative patterns of conservation and divergence linked to phenotypic evolution in the Nakaseomyces lineage.