Two long read-based genome assembly and annotation of polyploidy woody plants, Hibiscus syriacus L. using PacBio and Nanopore platforms.

Improvements in long read DNA sequencing and related techniques facilitated the generation of complex eukaryotic genomes. Despite these advances, the quality of constructed plant reference genomes remains relatively poor due to the large size of genomes, high content of repetitive sequences, and wide variety of ploidy. Here, we developed the de novo sequencing and assembly of high polyploid plant genome, Hibiscus syriacus, a flowering plant species of the Malvaceae family, using the Oxford Nanopore Technologies and Pacific Biosciences Sequel sequencing platforms. We investigated an efficient combination of high-quality and high-molecular-weight DNA isolation procedure and suitable assembler to achieve optimal results using long read sequencing data. We found that abundant ultra-long reads allow for large and complex polyploid plant genome assemblies with great recovery of repetitive sequences and error correction even at relatively low depth Nanopore sequencing data and polishing compared to previous studies. Collectively, our combination provides cost effective methods to improve genome continuity and quality compared to the previously reported reference genome by accessing highly repetitive regions. The application of this combination may enable genetic research and breeding of polyploid crops, thus leading to improvements in crop production.

A comparative genomic analysis of lichen-forming fungi reveals new insights into fungal lifestyles.

Lichen-forming fungi are mutualistic symbionts of green algae or cyanobacteria. We report the comparative analysis of six genomes of lichen-forming fungi in classes Eurotiomycetes and Lecanoromycetes to identify genomic information related to their symbiotic lifestyle. The lichen-forming fungi exhibited genome reduction via the loss of dispensable genes encoding plant-cell-wall-degrading enzymes, sugar transporters, and transcription factors. The loss of these genes reflects the symbiotic biology of lichens, such as the absence of pectin in the algal cell wall and obtaining specific sugars from photosynthetic partners. The lichens also gained many lineage- and species-specific genes, including those encoding small secreted proteins. These genes are primarily induced during the early stage of lichen symbiosis, indicating their significant roles in the establishment of lichen symbiosis.Our findings provide comprehensive genomic information for six lichen-forming fungi and novel insights into lichen biology and the evolution of symbiosis.

Alternative splicing diversifies the transcriptome and proteome of the rice blast fungus during host infection.

Alternative splicing (AS) contributes to diversifying and regulating cellular responses to environmental conditions and developmental cues by differentially producing multiple mRNA and protein isoforms from a single gene. Previous studies on AS in pathogenic fungi focused on profiling AS isoforms under a limited number of conditions. We analysed AS profiles in the rice blast fungus Magnaporthe oryzae, a global threat to rice production, using high-quality transcriptome data representing its vegetative growth (mycelia) and multiple host infection stages. We identified 4,270 AS isoforms derived from 2,413 genes, including 499 genes presumably regulated by infection-specific AS. AS appears to increase during infection, with 32.7% of the AS isoforms being produced during infection but absent in mycelia. Analysis of the isoforms observed at each infection stage showed that 636 AS isoforms were more abundant than corresponding annotated mRNAs, especially after initial hyphal penetration into host cell. Many such dominant isoforms were predicted to encode regulatory proteins such as transcription factors and phospho-transferases. We also identified the genes encoding distinct proteins via AS and confirmed the translation of some isoforms via a proteomic analysis, suggesting potential AS-mediated neo-functionalization of some genes during infection. Comprehensive profiling of the pattern of genome-wide AS during multiple stages of rice-M. oryzae interaction established a foundational resource that will help investigate the role and regulation of AS during rice infection.

Virus Isolation from the First Patient with SARS-CoV-2 in Korea.

Novel coronavirus (SARS-CoV-2) is found to cause a large outbreak started from Wuhan since December 2019 in China and SARS-CoV-2 infections have been reported with epidemiological linkage to China in 25 countries until now. We isolated SARS-CoV-2 from the oropharyngeal sample obtained from the patient with the first laboratory-confirmed SARS-CoV-2 infection in Korea. Cytopathic effects of SARS-CoV-2 in the Vero cell cultures were confluent 3 days after the first blind passage of the sample. Coronavirus was confirmed with spherical particle having a fringe reminiscent of crown on transmission electron microscopy. Phylogenetic analyses of whole genome sequences showed that it clustered with other SARS-CoV-2 reported from Wuhan.

Transcriptome Profiling of the Rice Blast Fungus Magnaporthe oryzae and Its Host Oryza sativa During Infection.

The rice blast (fungal pathogen: Magnaporthe oryzae and host: Oryza sativa) is one of the most important model pathosystems for understanding plant-microbe interactions. Although both genome sequences were published as the first cases of pathogen and host, only a few in planta transcriptome data during infection are available. Due to technical difficulties, previously reported fungal transcriptome data are not highly qualified to comprehensively profile the expression of fungal genes during infection. Here, we report the high-quality transcriptomes of M. oryzae and rice during infection using a sheath infection-based RNA sequencing approach. This comprehensive expression profiling of the fungal pathogen and its host will provide a better platform for understanding the plant-microbe interactions at the genomic level and serve as a valuable resource for the research community.

Comparative analysis of genome sequences of the conifer tree pathogen, Heterobasidion annosum s.s.

The causal agent of root and butt rot of conifer trees, Heterobasidion annosum, is widespread in boreal forests and economically responsible for annual loss of approximately 50 million euros to forest industries in Finland alone and much more at European level. In order to further understand the pathobiology of this fungus at the genome level, a Finnish isolate of H. annosum sensu stricto (isolate 03012) was sequenced and analyzed with the genome sequences of 23 white-rot and 13 brown-rot fungi. The draft genome assembly of H. annosum has a size of 31.01 Mb, containing 11,453 predicted genes. Whole genome alignment showed that 84.38% of H. annosum genome sequences were aligned with those of previously sequenced H. irregulare TC 32-1 counterparts. The result is further supported by the protein sequence clustering analysis which revealed that the two genomes share 6719 out of 8647 clusters. When sequencing reads of H. annosum were aligned against the genome sequences of H. irregulare, six single nucleotide polymorphisms were found in every 1 kb, on average. In addition, 98.68% of SNPs were found to be homo-variants, suggesting that the two species have long evolved from different niches. Gene family analysis revealed that most of the white-rot fungi investigated had more gene families involved in lignin degradation or modification, including laccases and peroxidase. Comparative analysis of the two Heterobasidion spp. as well as white-/brown-rot fungi would provide new insights for understanding the pathobiology of the conifer tree pathogen.

Draft Genome Sequence of the Fungus Associated with Oak Wilt Mortality in South Korea, Raffaelea quercus-mongolicae KACC44405.

The fungus Raffaelea quercus-mongolicae is the causal agent of Korean oak wilt, a disease associated with mass mortality of oak trees (e.g., Quercus spp.). The fungus is vectored and dispersed by the ambrosia beetle, Platypus koryoensis Here, we present the 27.0-Mb draft genome sequence of R. quercus-mongolicae strain KACC44405.

Whole-genome de novo sequencing, combined with RNA-Seq analysis, reveals unique genome and physiological features of the amylolytic yeast Saccharomycopsis fibuligera and its interspecies hybrid.

BACKGROUND: Genomic studies on fungal species with hydrolytic activity have gained increased attention due to their great biotechnological potential for biomass-based biofuel production. The amylolytic yeast Saccharomycopsis fibuligera has served as a good source of enzymes and genes involved in saccharification. Despite its long history of use in food fermentation and bioethanol production, very little is known about the basic physiology and genomic features of S. fibuligera. RESULTS: We performed whole-genome (WG) de novo sequencing and complete assembly of S. fibuligera KJJ81 and KPH12, two isolates from wheat-based Nuruk in Korea. Intriguingly, the KJJ81 genome (~38 Mb) was revealed as a hybrid between the KPH12 genome (~18 Mb) and another unidentified genome sharing 88.1% nucleotide identity with the KPH12 genome. The seven chromosome pairs of KJJ81 subgenomes exhibit highly conserved synteny, indicating a very recent hybridization event. The phylogeny inferred from WG comparisons showed an early divergence of S. fibuligera before the separation of the CTG and Saccharomycetaceae clades in the subphylum Saccharomycotina. Reconstructed carbon and sulfur metabolic pathways, coupled with RNA-Seq analysis, suggested a marginal Crabtree effect under high glucose and activation of sulfur metabolism toward methionine biosynthesis under sulfur limitation in this yeast. Notably, the lack of sulfate assimilation genes in the S. fibuligera genome reflects a unique phenotype for Saccharomycopsis clades as natural sulfur auxotrophs. Extended gene families, including novel genes involved in saccharification and proteolysis, were identified. Moreover, comparative genome analysis of S. fibuligera ATCC 36309, an isolate from chalky rye bread in Germany, revealed that an interchromosomal translocation occurred in the KPH12 genome before the generation of the KJJ81 hybrid genome. CONCLUSIONS: The completely sequenced S. fibuligera genome with high-quality annotation and RNA-Seq analysis establishes an important foundation for functional inference of S. fibuligera in the degradation of fermentation mash. The gene inventory facilitates the discovery of new genes applicable to the production of novel valuable enzymes and chemicals. Moreover, as the first gapless genome assembly in the genus Saccharomycopsis including members with desirable traits for bioconversion, the unique genomic features of S. fibuligera and its hybrid will provide in-depth insights into fungal genome dynamics as evolutionary adaptation.

Whole genome de novo sequencing and genome annotation of the world popular cultivated edible mushroom, Lentinula edodes.

Lentinula edodes, the popular shiitake mushroom, is one of the most important cultivated edible mushrooms. It is used as a food and for medicinal purposes. Here, we present the 46.1 Mb draft genome of L. edodes, comprising 13,028 predicted gene models. The genome assembly consists of 31 scaffolds. Gene annotation provides key information about various signaling pathways and secondary metabolites. This genomic information should help establish the molecular genetic markers for MAS/MAB and increase our understanding of the genome structure and function.

Genome-Wide Analysis of Hypoxia-Responsive Genes in the Rice Blast Fungus, Magnaporthe oryzae.

Rice blast fungus, Magnaporthe oryzae, is the most destructive pathogen in the rice-growing area. This fungus has a biotrophic phase early in infection and later switches to a necrotrophic lifestyle. During the biotrophic phase, the fungus competes with its host for nutrients and oxygen. Continuous uptake of oxygen is essential for successful establishment of blast disease of this pathogen. Here, we report transcriptional responses of the fungus to oxygen limitation. Transcriptome analysis using RNA-Seq identified that 1,047 genes were up-regulated in response to hypoxia. Those genes are involved in mycelial development, sterol biosynthesis, and metal ion transport based on hierarchical GO terms, and are well-conserved among three fungal species. In addition, null mutants of two hypoxia-responsive genes were generated and their roles in fungal development and pathogenicity tested. The mutant for the sterol regulatory element-binding protein gene, MoSRE1, exhibited increased sensitivity to a hypoxia-mimicking agent, increased conidiation, and delayed invasive growth within host cells, which is suggestive of important roles in fungal development. However, such defects did not cause any significant decrease in disease severity. The other null mutant, for the alcohol dehydrogenase gene MoADH1, showed no defect in the hypoxia-mimicking condition (using cobalt chloride) and fungal development. Taken together, this comprehensive transcriptional profiling in response to a hypoxic condition with experimental validations would provide new insights into fungal development and pathogenicity in plant pathogenic fungi.

Genome-wide profiling of DNA methylation provides insights into epigenetic regulation of fungal development in a plant pathogenic fungus, Magnaporthe oryzae.

DNA methylation is an important epigenetic modification that regulates development of plants and mammals. To investigate the roles of DNA methylation in fungal development, we profiled genome-wide methylation patterns at single-nucleotide resolution during vegetative growth, asexual reproduction, and infection-related morphogenesis in a model plant pathogenic fungus, Magnaporthe oryzae. We found that DNA methylation occurs in and around genes as well as transposable elements and undergoes global reprogramming during fungal development. Such reprogramming of DNA methylation suggests that it may have acquired new roles other than controlling the proliferation of TEs. Genetic analysis of DNA methyltransferase deletion mutants also indicated that proper reprogramming in methylomes is required for asexual reproduction in the fungus. Furthermore, RNA-seq analysis showed that DNA methylation is associated with transcriptional silencing of transposable elements and transcript abundance of genes in context-dependent manner, reinforcing the role of DNA methylation as a genome defense mechanism. This comprehensive approach suggests that DNA methylation in fungi can be a dynamic epigenetic entity contributing to fungal development and genome defense. Furthermore, our DNA methylomes provide a foundation for future studies exploring this key epigenetic modification in fungal development and pathogenesis.

dbCRY: a Web-based comparative and evolutionary genomics platform for blue-light receptors.

Cryptochromes are flavoproteins that play a central role in the circadian oscillations of all living organisms except archaea. Cryptochromes are clustered into three subfamilies: plant-type cryptochromes, animal-type cryptochromes and cryptochrome-DASH proteins. These subfamilies are composed of photolyase/cryptochrome superfamily with 6-4 photolyase and cyclobutane pyrimidine dimer photolyase. Cryptochromes have conserved domain architectures with two distinct domains, an N-terminal photolyase-related domain and a C-terminal domain. Although the molecular function and domain architecture of cryptochromes are conserved, their molecular mechanisms differ between plants and animals. Thus, cryptochromes are one of the best candidates for comparative and evolutionary studies. Here, we have developed a Web-based platform for comparative and evolutionary studies of cryptochromes, dbCRY (http://www.dbcryptochrome.org/). A pipeline built upon the consensus domain profile was applied to 1438 genomes and identified 1309 genes. To support comparative and evolutionary genomics studies, the Web interface provides diverse functions such as (i) browsing by species, (ii) protein domain analysis, (iii) multiple sequence alignment, (iv) homology search and (v) extended analysis opportunities through the implementation of 'Favorite Browser' powered by the Comparative Fungal Genomics Platform 2.0 (CFGP 2.0; http://cfgp.snu.ac.kr/). dbCRY would serve as a standardized and systematic solution for cryptochrome genomics studies. Database URL: http://www.dbcryptochrome.org/

Draft Genome Sequence of Endocarpon pusillum Strain KoLRILF000583.

Lichen Endocarpon pusillum is a highly desiccation-tolerant and dominant species in biological soil crusts in arid and semi-arid regions. We report the draft genome sequence of a lichen-forming fungus, E. pusillum strain KoLRILF000583. The draft genome assembly has a size of 37,173,200 bp with a GC content of 49.71%, consisting of 40 scaffolds.

Draft Genome Sequence of Umbilicaria muehlenbergii KoLRILF000956, a Lichen-Forming Fungus Amenable to Genetic Manipulation.

Umbilicaria muehlenbergii strain KoLRILF000956 is amenable to Agrobacterium tumefaciens-mediated transformation (ATMT), making it the only known genetically tractable lichen-forming fungus to date. We report another advancement in lichen genetics, a draft genome assembly for U. muehlenbergii with a size of 34,812,353 bp and a GC content of 47.12%, consisting of seven scaffolds.

Draft Genome Sequence of Lichen-Forming Fungus Cladonia metacorallifera Strain KoLRI002260.

The lichen-forming fungus Cladonia metacorallifera strain KoLRI002260 is capable of producing a number of secondary metabolites, including usnic, didymic, and squamatic acids, which have antitumor, antioxidant, and antibiotic activities. The draft genome assembly has a size of 36,682,060 bp, with a G+C content of 44.91%, and consists of 30 scaffolds.

Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species.

Hot pepper (Capsicum annuum), one of the oldest domesticated crops in the Americas, is the most widely grown spice crop in the world. We report whole-genome sequencing and assembly of the hot pepper (Mexican landrace of Capsicum annuum cv. CM334) at 186.6× coverage. We also report resequencing of two cultivated peppers and de novo sequencing of the wild species Capsicum chinense. The genome size of the hot pepper was approximately fourfold larger than that of its close relative tomato, and the genome showed an accumulation of Gypsy and Caulimoviridae family elements. Integrative genomic and transcriptomic analyses suggested that change in gene expression and neofunctionalization of capsaicin synthase have shaped capsaicinoid biosynthesis. We found differential molecular patterns of ripening regulators and ethylene synthesis in hot pepper and tomato. The reference genome will serve as a platform for improving the nutritional and medicinal values of Capsicum species.

Experimental evolution reveals genome-wide spectrum and dynamics of mutations in the rice blast fungus, Magnaporthe oryzae.

Knowledge on mutation processes is central to interpreting genetic analysis data as well as understanding the underlying nature of almost all evolutionary phenomena. However, studies on genome-wide mutational spectrum and dynamics in fungal pathogens are scarce, hindering our understanding of their evolution and biology. Here, we explored changes in the phenotypes and genome sequences of the rice blast fungus Magnaporthe oryzae during the forced in vitro evolution by weekly transfer of cultures on artificial media. Through combination of experimental evolution with high throughput sequencing technology, we found that mutations accumulate rapidly prior to visible phenotypic changes and that both genetic drift and selection seem to contribute to shaping mutational landscape, suggesting the buffering capacity of fungal genome against mutations. Inference of mutational effects on phenotypes through the use of T-DNA insertion mutants suggested that at least some of the DNA sequence mutations are likely associated with the observed phenotypic changes. Furthermore, our data suggest oxidative damages and UV as major sources of mutation during subcultures. Taken together, our work revealed important properties of original source of variation in the genome of the rice blast fungus. We believe that these results provide not only insights into stability of pathogenicity and genome evolution in plant pathogenic fungi but also a model in which evolution of fungal pathogens in natura can be comparatively investigated.

Global expression profiling of transcription factor genes provides new insights into pathogenicity and stress responses in the rice blast fungus.

Because most efforts to understand the molecular mechanisms underpinning fungal pathogenicity have focused on studying the function and role of individual genes, relatively little is known about how transcriptional machineries globally regulate and coordinate the expression of a large group of genes involved in pathogenesis. Using quantitative real-time PCR, we analyzed the expression patterns of 206 transcription factor (TF) genes in the rice blast fungus Magnaporthe oryzae under 32 conditions, including multiple infection-related developmental stages and various abiotic stresses. The resulting data, which are publicly available via an online platform, provided new insights into how these TFs are regulated and potentially work together to control cellular responses to a diverse array of stimuli. High degrees of differential TF expression were observed under the conditions tested. More than 50% of the 206 TF genes were up-regulated during conidiation and/or in conidia. Mutations in ten conidiation-specific TF genes caused defects in conidiation. Expression patterns in planta were similar to those under oxidative stress conditions. Mutants of in planta inducible genes not only exhibited sensitive to oxidative stress but also failed to infect rice. These experimental validations clearly demonstrated the value of TF expression patterns in predicting the function of individual TF genes. The regulatory network of TF genes revealed by this study provides a solid foundation for elucidating how M. oryzae regulates its pathogenesis, development, and stress responses.

CFGP 2.0: a versatile web-based platform for supporting comparative and evolutionary genomics of fungi and Oomycetes.

In 2007, Comparative Fungal Genomics Platform (CFGP; http://cfgp.snu.ac.kr/) was publicly open with 65 genomes corresponding to 58 fungal and Oomycete species. The CFGP provided six bioinformatics tools, including a novel tool entitled BLASTMatrix that enables search homologous genes to queries in multiple species simultaneously. CFGP also introduced Favorite, a personalized virtual space for data storage and analysis with these six tools. Since 2007, CFGP has grown to archive 283 genomes corresponding to 152 fungal and Oomycete species as well as 201 genomes that correspond to seven bacteria, 39 plants and 105 animals. In addition, the number of tools in Favorite increased to 27. The Taxonomy Browser of CFGP 2.0 allows users to interactively navigate through a large number of genomes according to their taxonomic positions. The user interface of BLASTMatrix was also improved to facilitate subsequent analyses of retrieved data. A newly developed genome browser, Seoul National University Genome Browser (SNUGB), was integrated into CFGP 2.0 to support graphical presentation of diverse genomic contexts. Based on the standardized genome warehouse of CFGP 2.0, several systematic platforms designed to support studies on selected gene families have been developed. Most of them are connected through Favorite to allow of sharing data across the platforms.

Remodeling of global transcription patterns of Cryptococcus neoformans genes mediated by the stress-activated HOG signaling pathways.

The ability to sense and adapt to a hostile host environment is a crucial element for virulence of pathogenic fungi, including Cryptococcus neoformans. These cellular responses are evoked by diverse signaling cascades, including the stress-activated HOG pathway. Despite previous analysis of central components of the HOG pathway, its downstream signaling network is poorly characterized in C. neoformans. Here we performed comparative transcriptome analysis with HOG signaling mutants to explore stress-regulated genes and their correlation with the HOG pathway in C. neoformans. In this study, we not only provide important insights into remodeling patterns of global gene expression for counteracting external stresses but also elucidate novel characteristics of the HOG pathway in C. neoformans. First, inhibition of the HOG pathway increases expression of ergosterol biosynthesis genes and cellular ergosterol content, conferring a striking synergistic antifungal activity with amphotericin B and providing an excellent opportunity to develop a novel therapeutic method for treatment of cryptococcosis. Second, a number of cadmium-sensitive genes are differentially regulated by the HOG pathway, and their mutation causes resistance to cadmium. Finally, we have discovered novel stress defense and HOG-dependent genes, which encode a sodium/potassium efflux pump, protein kinase, multidrug transporter system, and elements of the ubiquitin-dependent system.