Transposable elements in filamentous fungi

Update on the Basic Understanding of Fusarium graminearum Virulence Factors in Common Wheat Research

Wheat is one of the most important food crops, both in China and worldwide. Wheat production is facing extreme stresses posed by different diseases, including Fusarium head blight (FHB), which has recently become an increasingly serious concerns. FHB is one of the most significant and destructive diseases affecting wheat crops all over the world. Recent advancements in genomic tools provide a new avenue for the study of virulence factors in relation to the host plants. The current review focuses on recent progress in the study of different strains of Fusarium infection. The presence of genome-wide repeat-induced point (RIP) mutations causes genomic mutations, eventually leading to host plant susceptibility against Fusarium invasion. Furthermore, effector proteins disrupt the host plant resistance mechanism. In this study, we proposed systematic modification of the host genome using modern biological tools to facilitate plant resistance against foreign invasion. We also suggested a number of scientific strategies, such as gene cloning, developing more powerful functional markers, and using haplotype marker-assisted selection, to further improve FHB resistance and associated breeding methods.

An Antarctic lichen isolate (Cladonia borealis) genome reveals potential adaptation to extreme environments

Cladonia borealis is a lichen that inhabits Antarctica's harsh environment. We sequenced the whole genome of a C. borealis culture isolated from a specimen collected in Antarctica using long-read sequencing technology to identify specific genetic elements related to its potential environmental adaptation. The final genome assembly produced 48 scaffolds, the longest being 2.2 Mbp, a 1.6 Mbp N50 contig length, and a 36 Mbp total length. A total of 10,749 protein-coding genes were annotated, containing 33 biosynthetic gene clusters and 102 carbohydrate-active enzymes. A comparative genomics analysis was conducted on six Cladonia species, and the genome of C. borealis exhibited 45 expanded and 50 contracted gene families. We identified that C. borealis has more Copia transposable elements and expanded transporters (ABC transporters and magnesium transporters) compared to other Cladonia species. Our results suggest that these differences contribute to C. borealis' remarkable adaptability in the Antarctic environment. This study also provides a useful resource for the genomic analysis of lichens and genetic insights into the survival of species isolated from Antarctica.

Mutations in the Second Alternative Oxidase Gene: A New Approach to Group Aspergillus niger Strains

Alternative oxidase is a terminal oxidase in the branched mitochondrial electron transport chain of most fungi including Aspergillus niger (subgenus Circumdati, section Nigri). A second, paralogous aox gene (aoxB) is extant in some A. niger isolates but also present in two divergent species of the subgenus Nidulantes-A. calidoustus and A. implicatus-as well as in Penicillium swiecickii. Black aspergilli are cosmopolitan opportunistic fungi that can cause diverse mycoses and acute aspergillosis in immunocompromised individuals. Amongst the approximately 75 genome-sequenced A. niger strains, aoxB features considerable sequence variation. Five mutations were identified that rationally affect transcription or function or terminally modify the gene product. One mutant allele that occurs in CBS 513.88 and A. niger neotype strain CBS 554.65 involves a chromosomal deletion that removes exon 1 and intron 1 from aoxB. Another aoxB allele results from retrotransposon integration. Three other alleles result from point mutations: a missense mutation of the start codon, a frameshift, and a nonsense mutation. A. niger strain ATCC 1015 has a full-length aoxB gene. The A. niger sensu stricto complex can thus be subdivided into six taxa according to extant aoxB allele, which may facilitate rapid and accurate identification of individual species.

Recent transposable element bursts are associated with the proximity to genes in a fungal plant pathogen

The activity of transposable elements (TEs) contributes significantly to pathogen genome evolution. TEs often destabilize genome integrity but may also confer adaptive variation in pathogenicity or resistance traits. De-repression of epigenetically silenced TEs often initiates bursts of transposition activity that may be counteracted by purifying selection and genome defenses. However, how these forces interact to determine the expansion routes of TEs within a pathogen species remains largely unknown. Here, we analyzed a set of 19 telomere-to-telomere genomes of the fungal wheat pathogen Zymoseptoria tritici. Phylogenetic reconstruction and ancestral state estimates of individual TE families revealed that TEs have undergone distinct activation and repression periods resulting in highly uneven copy numbers between genomes of the same species. Most TEs are clustered in gene poor niches, indicating strong purifying selection against insertions near coding sequences, or as a consequence of insertion site preferences. TE families with high copy numbers have low sequence divergence and strong signatures of defense mechanisms (i.e., RIP). In contrast, small non-autonomous TEs (i.e., MITEs) are less impacted by defense mechanisms and are often located in close proximity to genes. Individual TE families have experienced multiple distinct burst events that generated many nearly identical copies. We found that a Copia element burst was initiated from recent copies inserted substantially closer to genes compared to older copies. Overall, TE bursts tended to initiate from copies in GC-rich niches that escaped inactivation by genomic defenses. Our work shows how specific genomic environments features provide triggers for TE proliferation in pathogen genomes.

Investigating the origin of subtelomeric and centromeric AT-rich elements in Aspergillus flavus

An in silico study of Aspergillus flavus genome stability uncovered significant variations in both coding and non-coding regions. The non-coding insertions uniformly consisted of AT-rich sequences that are evolutionarily maintained, albeit distributed at widely different sites in an array of A. flavus strains. A survey of ≥ 2kb AT-rich elements (AT ≥ 70%; ATEs) in non-centromeric regions uncovered two major categories of ATEs. The first category is composed of homologous insertions at ectopic, non-allelic sites that contain homology to transposable elements (TEs; Classes B, C, D, and E). Strains differed significantly in frequency, position, and TE type, but displayed a common enrichment in subtelomeric regions. The TEs were heavily mutated, with patterns consistent with the ancestral activity of repeat-induced point mutations (RIP). The second category consists of a conserved set of novel subtelomeric ATE repeats (Classes A, G, G, H, I and J) which lack discernible TEs and, unlike TEs, display a constant polarity relative to the telomere. Members of one of these classes are derivatives of a progenitor ATE that is predicted to have undergone extensive homologous recombination during evolution. A third category of ATEs consists of ~100 kb regions at each centromere. Centromeric ATEs and TE clusters within these centromeres display a high level of sequence identity between strains. These studies suggest that transposition and RIP are forces in the evolution of subtelomeric and centromeric structure and function.

The first discovery of Tc1 transposons in yeast

Background

Identification of transposons without close homologs is still a difficult task. IS630/Tc1/mariner transposons, classified into a superfamily, are probably the most widespread DNA transposons in nature. Tc1/mariner transposons have been discovered in animals, plants, and filamentous fungi, however, not in yeast.

Results

In the present study, we report the discovery of two intact Tc1 transposons in yeast and filamentous fungi, respectively. The first one, named Tc1-OP1 (DD40E), represents Tc1 transposons in Ogataea parapolymorpha. The second one, named Tc1-MP1 (DD34E), represents Tc1 transposons in the Rhizopodaceae and Mucoraceae families. As a homolog of Tc1-OP1 and Tc1-MP1, IS630-AB1 (DD34E) was discovered as an IS630 transposon in Acinetobacter spp.

Conclusion

Tc1-OP1 is not only the first reported Tc1 transposon in yeast, but also the first reported nonclassical Tc1 transposon. Tc1-OP1 is the largest of IS630/Tc1/mariner transposons reported to date and significantly different from others. Notably, Tc1-OP1 encodes a serine-rich domain and a transposase, extending the current knowledge of Tc1 transposons. The phylogenetic relationships of Tc1-OP1, Tc1-MP1 and IS630-AB1 indicated that these transposons had evolved from a common ancestor. Tc1-OP1, Tc1-MP1 and IS630-AB1 can be used as reference sequences to facilitate the identification of IS630/Tc1/mariner transposons. More Tc1/mariner transposons will be identified in yeast, following our discovery.

De Novo Long-Read Whole-Genome Assemblies and the Comparative Pan-Genome Analysis of Ascochyta Blight Pathogens Affecting Field Pea

Ascochyta Blight (AB) is a major disease of many cool-season legumes globally. In field pea, three fungal pathogens have been identified to be responsible for this disease in Australia, namely Peyronellaea pinodes, Peyronellaea pinodella and Phoma koolunga. Limited genomic resources for these pathogens have been generated, which has hampered the implementation of effective management strategies and breeding for resistant cultivars. Using Oxford Nanopore long-read sequencing, we report the first high-quality, fully annotated, near-chromosome-level nuclear and mitochondrial genome assemblies for 18 isolates from the Australian AB complex. Comparative genome analysis was performed to elucidate the differences and similarities between species and isolates using phylogenetic relationships and functional diversity. Our data indicated that P. pinodella and P. koolunga are heterothallic, while P. pinodes is homothallic. More homology and orthologous gene clusters are shared between P. pinodes and P. pinodella compared to P. koolunga. The analysis of the repetitive DNA content showed differences in the transposable repeat composition in the genomes and their expression in the transcriptomes. Significant repeat expansion in P. koolunga’s genome was seen, with strong repeat-induced point mutation (RIP) activity being evident. Phylogenetic analysis revealed that genetic diversity can be exploited for species marker development. This study provided the much-needed genetic resources and characterization of the AB species to further drive research in key areas such as disease epidemiology and host–pathogen interactions.

Chromosomal polymorphism of the Ceratocystis fimbriata species complex in Brazil

Ceratocystis fimbriata is an important pathogen that causes wilt in several plant species. Despite the importance of this pathogen, knowledge about its karyotypic polymorphism and genomic architecture is limited. The main objective of this study was to investigate the karyotype of isolates of the C. fimbriata species complex from different host plants and geographical origins in Brazil. First, the identity of the isolates was confirmed conducting multilocus sequence analysis (MLSA) phylogeny using β-tubulin (TUBB), translation elongation factor 1α (TEF-1α) and mating-type (MAT1 and MAT2) gene sequences. To investigate the chromosomal polymorphism, two conditions of pulsed-field gel electrophoresis (PFGE) were used and the karyotypes of the isolates obtained. The retrotransposon-microsatellite amplified polymorphism (REMAP) molecular marker was utilized to assess the genetic variability among isolates. In the MLSA utilizing the concatenated gene sequences, Ceratocystis cacaofunesta and C. fimbriata formed separate clades, but considerable variation among C. fimbriata isolates was observed. Polymorphism in chromosome number and size was found, indicating the existence of genomic differences among isolates and occurrence of chromosomal rearrangements in the species complex. The number of chromosomes varied from seven to nine and the estimated minimum chromosome sizes were estimated to be between 2.7 to 6.0 Mbp. Small polymorphic chromosomes ranging from 1.2 to 1.8 Mbp were observed in all isolates, raising the hypothesis that they could be supernumerary chromosomes. REMAP analysis revealed a high genetic variability and that isolates from the same host tend to group together in a same cluster. Our results bring new insights into the chromosomal diversity and genome organization of the C. fimbriata complex.

Comparative Genome Analyses of Plant Rust Pathogen Genomes Reveal a Confluence of Pathogenicity Factors to Quell Host Plant Defense Responses

Switchgrass rust caused by Puccinia novopanici (P. novopanici) has the ability to significantly affect the biomass yield of switchgrass, an important biofuel crop in the United States. A comparative genome analysis of P. novopanici with rust pathogen genomes infecting monocot cereal crops wheat, barley, oats, maize and sorghum revealed the presence of larger structural variations contributing to their genome sizes. A comparative alignment of the rust pathogen genomes resulted in the identification of collinear and syntenic relationships between P. novopanici and P. sorghi; P. graminis tritici 21–0 (Pgt 21) and P. graminis tritici Ug99 (Pgt Ug99) and between Pgt 21 and P. triticina (Pt). Repeat element analysis indicated a strong presence of retro elements among different Puccinia genomes, contributing to the genome size variation between ~1 and 3%. A comparative look at the enriched protein families of Puccinia spp. revealed a predominant role of restriction of telomere capping proteins (RTC), disulfide isomerases, polysaccharide deacetylases, glycoside hydrolases, superoxide dismutases and multi-copper oxidases (MCOs). All the proteomes of Puccinia spp. share in common a repertoire of 75 secretory and 24 effector proteins, including glycoside hydrolases cellobiohydrolases, peptidyl-propyl isomerases, polysaccharide deacetylases and protein disulfide-isomerases, that remain central to their pathogenicity. Comparison of the predicted effector proteins from Puccinia spp. genomes to the validated proteins from the Pathogen–Host Interactions database (PHI-base) resulted in the identification of validated effector proteins PgtSR1 (PGTG_09586) from P. graminis and Mlp124478 from Melampsora laricis across all the rust pathogen genomes.

Exploring a novel genomic safe-haven site in the human pathogenic mould Aspergillus fumigatus

Aspergillus fumigatus is the most important airborne fungal pathogen and allergen of humans causing high morbidity and mortality worldwide. The factors that govern pathogenicity of this organism are multi-factorial and are poorly understood. Molecular tools to dissect the mechanisms of pathogenicity in A. fumigatus have improved significantly over the last 20 years however many procedures have not been standardised for A. fumigatus. Here, we present a new genomic safe-haven locus at the site of an inactivated transposon, named SH-aft4, which can be used to insert DNA sequences in the genome of this fungus without impacting its phenotype. We show that we are able to effectively express a transgene construct from the SH-aft4 and that natural regulation of promoter function is conserved at this site. Furthermore, the SH-aft4 locus is highly conserved in the genome of a wide range of clinical and environmental isolates including the isolates commonly used by many laboratories CEA10, Af293 and ATCC46645, allowing a wide range of isolates to be manipulated. Our results show that the aft4 locus can serve as a site for integration of a wide range of genetic constructs to aid functional genomics studies of this important human fungal pathogen.

Reciprocal recombination genomic signatures in the symbiotic arbuscular mycorrhizal fungi Rhizophagus irregularis

Arbuscular mycorrhizal fungi (AMF) are part of the most widespread fungal-plant symbiosis. They colonize at least 80% of plant species, promote plant growth and plant diversity. These fungi are multinucleated and contain either one or two haploid nuclear genotypes (monokaryon and dikaryon) identified by the alleles at a putative mating-type locus. This taxon has been considered as an ancient asexual scandal because of the lack of observable sexual structures. Despite identification of a putative mating-type locus and functional activation of genes related to mating when two isolates co-exist, it remains unknown if the AMF life cycle involves a sexual or parasexual stage. We used publicly available genome sequences to test if Rhizophagus irregularis dikaryon genomes display signatures of sexual reproduction in the form of reciprocal recombination patterns, or if they display exclusively signatures of parasexual reproduction involving gene conversion. We used short-read and long-read sequence data to identify nucleus-specific alleles within dikaryons and then compared them to orthologous gene sequences from related monokaryon isolates displaying the same putative MAT-types as the dikaryon. We observed that the two nucleus-specific alleles of the dikaryon A5 are more related to the homolog sequences of monokaryon isolates displaying the same putative MAT-type than between each other. We also observed that these nucleus-specific alleles displayed reciprocal recombination signatures. These results confirm that dikaryon and monokaryon isolates displaying the same putative MAT-type are related in their life-cycle. These results suggest that a genetic exchange mechanism, involving reciprocal recombination in dikaryon genomes, allows AMF to generate genetic diversity.

Intra-Species Genomic Variation in the Pine Pathogen Fusarium circinatum

Fusarium circinatum is an important global pathogen of pine trees. Genome plasticity has been observed in different isolates of the fungus, but no genome comparisons are available. To address this gap, we sequenced and assembled to chromosome level five isolates of F. circinatum. These genomes were analysed together with previously published genomes of F. circinatum isolates, FSP34 and KS17. Multi-sample variant calling identified a total of 461,683 micro variants (SNPs and small indels) and a total of 1828 macro structural variants of which 1717 were copy number variants and 111 were inversions. The variant density was higher on the sub-telomeric regions of chromosomes. Variant annotation revealed that genes involved in transcription, transport, metabolism and transmembrane proteins were overrepresented in gene sets that were affected by high impact variants. A core genome representing genomic elements that were conserved in all the isolates and a non-redundant pangenome representing all genomic elements is presented. Whole genome alignments showed that an average of 93% of the genomic elements were present in all isolates. The results of this study reveal that some genomic elements are not conserved within the isolates and some variants are high impact. The described genome-scale variations will help to inform novel disease management strategies against the pathogen.

Transposable Elements in the Genome of the Lichen-Forming Fungus Umbilicaria pustulata and Their Distribution in Different Climate Zones along Elevation

Transposable elements (TEs) are an important source of genome plasticity across the tree of life. Drift and natural selection are important forces shaping TE distribution and accumulation. Fungi, with their multifaceted phenotypic diversity and relatively small genome size, are ideal models to study the role of TEs in genome evolution and their impact on the host's ecological and life history traits. Here we present an account of all TEs found in a high-quality reference genome of the lichen-forming fungus Umbilicaria pustulata, a macrolichen species comprising two climatic ecotypes: Mediterranean and cold temperate. We trace the occurrence of the newly identified TEs in populations along three elevation gradients using a Pool-Seq approach to identify TE insertions of potential adaptive significance. We found that TEs cover 21.26% of the 32.9 Mbp genome, with LTR Gypsy and Copia clades being the most common TEs. We identified 28 insertions displaying consistent insertion frequency differences between the two host ecotypes across the elevation gradients. Most of the highly differentiated insertions were located near genes, indicating a putative function. This pioneering study of the content and climate niche-specific distribution of TEs in a lichen-forming fungus contributes to understanding the roles of TEs in fungal evolution.

Cryptic genetic structure and copy-number variation in the ubiquitous forest symbiotic fungus Cenococcum geophilum

Ectomycorrhizal (ECM) fungi associated with plants constitute one of the most successful symbiotic interactions in forest ecosystems. ECM support trophic exchanges with host plants and are important factors for the survival and stress resilience of trees. However, ECM clades often harbour morpho-species and cryptic lineages, with weak morphological differentiation. How this relates to intraspecific genome variability and ecological functioning is poorly known. Here, we analysed 16 European isolates of the ascomycete Cenococcum geophilum, an extremely ubiquitous forest symbiotic fungus with no known sexual or asexual spore-forming structures but with a massively enlarged genome. We carried out whole-genome sequencing to identify single-nucleotide polymorphisms. We found no geographic structure at the European scale but divergent lineages within sampling sites. Evidence for recombination was restricted to specific cryptic lineages. Lineage differentiation was supported by extensive copy-number variation. Finally, we confirmed heterothallism with a single MAT1 idiomorph per genome. Synteny analyses of the MAT1 locus revealed substantial rearrangements and a pseudogene of the opposite MAT1 idiomorph. Our study provides the first evidence for substantial genome-wide structural variation, lineage-specific recombination and low continent-wide genetic differentiation in C. geophilum. Our study provides a foundation for targeted analyses of intra-specific functional variation in this major symbiosis.

Internally Symmetrical Stwintrons and Related Canonical Introns in Hypoxylaceae Species

Spliceosomal introns are pervasive in eukaryotes. Intron gains and losses have occurred throughout evolution, but the origin of new introns is unclear. Stwintrons are complex intervening sequences where one of the sequence elements (5′-donor, lariat branch point element or 3′-acceptor) necessary for excision of a U2 intron (external intron) is itself interrupted by a second (internal) U2 intron. In Hypoxylaceae, a family of endophytic fungi, we uncovered scores of donor-disrupted stwintrons with striking sequence similarity among themselves and also with canonical introns. Intron–exon structure comparisons suggest that these stwintrons have proliferated within diverging taxa but also give rise to proliferating canonical introns in some genomes. The proliferated (stw)introns have integrated seamlessly at novel gene positions. The recently proliferated (stw)introns appear to originate from a conserved ancestral stwintron characterised by terminal inverted repeats (45–55 nucleotides), a highly symmetrical structure that may allow the formation of a double-stranded intron RNA molecule. No short tandem duplications flank the putatively inserted intervening sequences, which excludes a DNA transposition-based mechanism of proliferation. It is tempting to suggest that this highly symmetrical structure may have a role in intron proliferation by (an)other mechanism(s).

Historical genomics reveals the evolutionary mechanisms behind multiple outbreaks of the host-specific coffee wilt pathogen Fusarium xylarioides

BACKGROUND

Nearly 50% of crop yields are lost to pests and disease, with plants and pathogens locked in an amplified co-evolutionary process of disease outbreaks. Coffee wilt disease, caused by Fusarium xylarioides, decimated coffee production in west and central Africa following its initial outbreak in the 1920s. After successful management, it later re-emerged and by the 2000s comprised two separate epidemics on arabica coffee in Ethiopia and robusta coffee in east and central Africa.

RESULTS

Here, we use genome sequencing of six historical culture collection strains spanning 52 years to identify the evolutionary processes behind these repeated outbreaks. Phylogenomic reconstruction using 13,782 single copy orthologs shows that the robusta population arose from the initial outbreak, whilst the arabica population is a divergent sister clade to the other strains. A screen for putative effector genes involved in pathogenesis shows that the populations have diverged in gene content and sequence mainly by vertical processes within lineages. However, 15 putative effector genes show evidence of horizontal acquisition, with close homology to genes from F. oxysporum. Most occupy small regions of homology within wider scaffolds, whereas a cluster of four genes occupy a 20Kb scaffold with strong homology to a region on a mobile pathogenicity chromosome in F. oxysporum that houses known effector genes. Lacking a match to the whole mobile chromosome, we nonetheless found close associations with DNA transposons, especially the miniature impala type previously proposed to facilitate horizontal transfer of pathogenicity genes in F. oxysporum. These findings support a working hypothesis that the arabica and robusta populations partly acquired distinct effector genes via transposition-mediated horizontal transfer from F. oxysporum, which shares coffee as a host and lives on other plants intercropped with coffee.

CONCLUSION

Our results show how historical genomics can help reveal mechanisms that allow fungal pathogens to keep pace with our efforts to resist them. Our list of putative effector genes identifies possible future targets for fungal control. In turn, knowledge of horizontal transfer mechanisms and putative donor taxa might help to design future intercropping strategies that minimize the risk of transfer of effector genes between closely-related Fusarium taxa.

Phylogenetic and phenotypic characterization of Fusarium oxysporum f. sp. niveum isolates from Florida-grown watermelon

Fusarium wilt of watermelon (Citrullus lanatus) caused by Fusarium oxysporum f. sp. niveum (Fon), has become an increasing concern of farmers in the southeastern USA, especially in Florida. Management of this disease, most often through the use of resistant cultivars and crop rotation, requires an accurate understanding of an area’s pathogen population structure and phenotypic characteristics. This study improved the understanding of the state’s pathogen population by completing multilocus sequence analysis (MLSA) of two housekeeping genes (BT and TEF) and two loci (ITS and IGS), aggressiveness and race-determining bioassays on 72 isolates collected between 2011 and 2015 from major watermelon production areas in North, Central, and South Florida. Multilocus sequence analysis (MLSA) failed to group race 3 isolates into a single large clade; moreover, clade membership was not apparently correlated with aggressiveness (which varied both within and between clades), and only slightly with sampling location. The failure of multilocus sequence analysis using four highly conserved housekeeping genes and loci to clearly group and delineate known Fon races provides justification for future whole genome sequencing efforts whose more robust genomic comparisons will provide higher resolution of intra-species genetic distinctions. Consequently, these results suggest that identification of Fon isolates by race determination alone may fail to detect economically important phenotypic characteristics such as aggressiveness leading to inaccurate risk assessment.

Advances in Genetic Engineering Technology and Its Application in the Industrial Fungus Aspergillus oryzae

The filamentous fungus Aspergillus oryzae is an important strain in the traditional fermentation and food processing industries and is often used in the production of soy sauce, soybean paste, and liquor-making. In addition, A. oryzae has a strong capacity to secrete large amounts of hydrolytic enzymes; therefore, it has also been used in the enzyme industry as a cell factory for the production of numerous native and heterologous enzymes. However, the production and secretion of foreign proteins by A. oryzae are often limited by numerous bottlenecks that occur during transcription, translation, protein folding, translocation, degradation, transport, secretion, etc. The existence of these problems makes it difficult to achieve the desired target in the production of foreign proteins by A. oryzae. In recent years, with the decipherment of the whole genome sequence, basic research and genetic engineering technologies related to the production and utilization of A. oryzae have been well developed, such as the improvement of homologous recombination efficiency, application of selectable marker genes, development of large chromosome deletion technology, utilization of hyphal fusion techniques, and application of CRISPR/Cas9 genome editing systems. The development and establishment of these genetic engineering technologies provided a great deal of technical support for the industrial production and application of A. oryzae. This paper reviews the advances in basic research and genetic engineering technologies of the fermentation strain A. oryzae mentioned above to open up more effective ways and research space for the breeding of A. oryzae production strains in the future.

Comparative Genomics: Insights on the Pathogenicity and Lifestyle of Rhizoctonia solani

Proper management of agricultural disease is important to ensure sustainable food security. Staple food crops like rice, wheat, cereals, and other cash crops hold great export value for countries. Ensuring proper supply is critical; hence any biotic or abiotic factors contributing to the shortfall in yield of these crops should be alleviated. Rhizoctonia solani is a major biotic factor that results in yield losses in many agriculturally important crops. This paper focuses on genome informatics of our Malaysian Draft R. solani AG1-IA, and the comparative genomics (inter- and intra- AG) with four AGs including China AG1-IA (AG1-IA_KB317705.1), AG1-IB, AG3, and AG8. The genomic content of repeat elements, transposable elements (TEs), syntenic genomic blocks, functions of protein-coding genes as well as core orthologous genic information that underlies R. solani’s pathogenicity strategy were investigated. Our analyses show that all studied AGs have low content and varying profiles of TEs. All AGs were dominant for Class I TE, much like other basidiomycete pathogens. All AGs demonstrate dominance in Glycoside Hydrolase protein-coding gene assignments suggesting its importance in infiltration and infection of host. Our profiling also provides a basis for further investigation on lack of correlation observed between number of pathogenicity and enzyme-related genes with host range. Despite being grouped within the same AG with China AG1-IA, our Draft AG1-IA exhibits differences in terms of protein-coding gene proportions and classifications. This implies that strains from similar AG do not necessarily have to retain similar proportions and classification of TE but must have the necessary arsenal to enable successful infiltration and colonization of host. In a larger perspective, all the studied AGs essentially share core genes that are generally involved in adhesion, penetration, and host colonization. However, the different infiltration strategies will depend on the level of host resilience where this is clearly exhibited by the gene sets encoded for the process of infiltration, infection, and protection from host.

Lineage-Specific Genes and Cryptic Sex: Parallels and Differences between Arbuscular Mycorrhizal Fungi and Fungal Pathogens

Arbuscular mycorrhizal fungi (AMF) live as obligate root symbionts on almost all land plants. They have long been regarded as ancient asexuals that have propagated clonally for millions of years. However, genomic studies in Rhizophagus irregularis and other AMF revealed many features indicative of sex. Surprisingly, comparative genomics of conspecific isolates of R. irregularis revealed an unexpected interstrain diversity, suggesting that AMF carry a high number of lineage-specific (LS) genes. Intriguingly, cryptic sex and LS genomic regions have previously been reported in a number of fungal pathogens of plants and humans. Here, we discuss these genomic similarities and highlight their potential relevance for AMF adaptation to the environment and for symbiotic functioning.

Degenerated Virulence and Irregular Development of Fusarium oxysporum f. sp. niveum Induced by Successive Subculture

Successive cultivation of fungi on artificial media has been reported to cause the sectorization, which leads to degeneration of developmental phenotype, and virulence. Fusarium oxysporum f. sp. niveum (Fon), the causal agent of watermelon Fusarium wilt, forms degenerated sectors after successive cultivation. In the present research, we demonstrated that subculture with aged mycelia increased the incidence of degenerations. To further investigate the differences between the Fon wild type (sporodochial type, ST) and variants (MT: mycelial type and PT: pionnotal type), developmental phenotypes and pathogenicity to watermelon were examined. Results in variants (PT2, PT3, PT11, and MT6) were different from ST with mycelia growth, conidia production and chlamydospore formation. Virulence of degenerated variants on susceptible watermelon Grand Baby (GB) cultivar was determined after inoculation with Fon variants and Fon ST. In root dipping methods, Fon variants showed no significant differences in disease progress compared with ST. Fon variants showed a significant decrease in disease progression compared with ST through infested soil inoculation. The contrasting results of two inoculation methods suggest that the degenerative changes due to repeated successive cultivation may lead to the loss of pathogen virulence-related factors of the early stage of Fon infection process. Therefore, cell wall-degrading enzymes (CWDEs; cellulase, pectinase, and xylanase) activities of different variants were analyzed. All Fon degenerated variants demonstrated significant decreases of CWDEs activities compared with ST. Additionally, transcript levels of 9 virulence-related genes (fmk1, fgb1, pacC, xlnR, pl1, rho1, gas1, wc1, and fow1) were assessed in normal state. The degenerated variants demonstrated a significantly low level of tested virulence-related gene transcripts except for fmk1, xlnR, and fow1. In summary, the degeneration of Fon is triggered with successive subculture through aged mycelia. The degeneration showed significant impacts on virulence to watermelon, which was correlated with the reduction of CWDEs activities and declining expression of a set of virulence-related genes.

Advances in Molecular Tools and In Vivo Models for the Study of Human Fungal Pathogenesis

Pathogenic fungi represent an increasing infectious disease threat to humans, especially with an increasing challenge of antifungal drug resistance. Over the decades, numerous tools have been developed to expedite the study of pathogenicity, initiation of disease, drug resistance and host-pathogen interactions. In this review, we highlight advances that have been made in the use of molecular tools using CRISPR technologies, RNA interference and transposon targeted mutagenesis. We also discuss the use of animal models in modelling disease of human fungal pathogens, focusing on zebrafish, the silkworm, Galleria mellonella and the murine model.

The utility of iPBS retrotransposons markers to analyze genetic variation in yeast

Yeasts are one of the main organisms in the food industry and effective components of many ecosystems. The method for identifying and detecting certain yeast species or strains is a crucial step for the food industry and should be simple, reliable, fast, and inexpensive. In our study, inter-priming binding sites (iPBS) retrotransposon marker system was employed to elucidate the genetic variability at intraspecies and interspecies levels among 112 yeast strains belonging to eight species previously obtained from fermented foods. The molecular identification of yeast strains was firstly confirmed by sequencing the D1/D2 domain of the 26S rRNA. The eight selected retrotransposon-based primers produced 278 bands, all of which were polymorphic with an average of 34.75 polymorphic fragments per primer. The averages of polymorphism information contents and the resolving power values for the iPBS marker system were 0.23 and 10.11, respectively. The genetic parameters within each yeast species obtained from iPBS markers were observed as; the percentage of polymorphic loci for each species ranging from 19.23% to 71.21%, Nei's gene diversity from 0.085 to 0.228, while Shannon's information index values ranging from 0.125 to 0.349. The value of gene flow (0.09) and genetic variation among the populations (0.85) showed higher genetic variation among the species. UPGMA analyses demonstrated considerable genetic variability in the yeast strains, clustered them according to their species, and revealed the intraspecific variation. Each of the selected iPBS primer provided enough species-discrimination. Present evaluations suggest the utility of iPBS marker system to estimate the genetic variation of yeast strains. This study is a preliminary point for further studies on the identification methodology, and population genetics of yeast species having importance in the food industry with iPBS markers.

Transposon mobilization in the human fungal pathogen Cryptococcus is mutagenic during infection and promotes drug resistance in vitro

When transitioning from the environment, pathogenic microorganisms must adapt rapidly to survive in hostile host conditions. This is especially true for environmental fungi that cause opportunistic infections in immunocompromised patients since these microbes are not well adapted human pathogens. Cryptococcus species are yeastlike fungi that cause lethal infections, especially in HIV-infected patients. Using Cryptococcus deneoformans in a murine model of infection, we examined contributors to drug resistance and demonstrated that transposon mutagenesis drives the development of 5-fluoroorotic acid (5FOA) resistance. Inactivation of target genes URA3 or URA5 primarily reflected the insertion of two transposable elements (TEs): the T1 DNA transposon and the TCN12 retrotransposon. Consistent with in vivo results, increased rates of mutagenesis and resistance to 5FOA and the antifungal drugs rapamycin/FK506 (rap/FK506) and 5-fluorocytosine (5FC) were found when Cryptococcus was incubated at 37° compared to 30° in vitro, a condition that mimics the temperature shift that occurs during the environment-to-host transition. Inactivation of the RNA interference (RNAi) pathway, which suppresses TE movement in many organisms, was not sufficient to elevate TE movement at 30° to the level observed at 37°. We propose that temperature-dependent TE mobilization in Cryptococcus is an important mechanism that enhances microbial adaptation and promotes pathogenesis and drug resistance in the human host.

Repeat-Induced Point Mutations Drive Divergence between Fusarium circinatum and Its Close Relatives

The Repeat-Induced Point (RIP) mutation pathway is a fungal-specific genome defense mechanism that counteracts the deleterious effects of transposable elements. This pathway permanently mutates its target sequences by introducing cytosine to thymine transitions. We investigated the genome-wide occurrence of RIP in the pitch canker pathogen, Fusarium circinatum, and its close relatives in the Fusarium fujikuroi species complex (FFSC). Our results showed that the examined fungi all exhibited hallmarks of RIP, but that they differed in terms of the extent to which their genomes were affected by this pathway. RIP mutations constituted a large proportion of all the FFSC genomes, including both core and dispensable chromosomes, although the latter were generally more extensively affected by RIP. Large RIP-affected genomic regions were also much more gene sparse than the rest of the genome. Our data further showed that RIP-directed sequence diversification increased the variability between homologous regions of related species, and that RIP-affected regions can interfere with homologous recombination during meiosis, thereby contributing to post-mating segregation distortion. Taken together, these findings suggest that RIP can drive the independent divergence of chromosomes, alter chromosome architecture, and contribute to the divergence among F. circinatum and other members of this economically important group of fungi.

Transposon-Mediated Horizontal Transfer of the Host-Specific Virulence Protein ToxA between Three Fungal Wheat Pathogens

This work dissects the tripartite horizontal transfer of ToxA, a gene that has a direct negative impact on global wheat yields. Defining the extent of horizontally transferred DNA is important because it can provide clues to the mechanisms that facilitate HGT. Our analysis of ToxA and its surrounding 14 kb suggests that this gene was horizontally transferred in two independent events, with one event likely facilitated by a type II DNA transposon. These horizontal transfer events are now in various processes of decay in each species due to the repeated insertion of new transposons and subsequent rounds of targeted mutation by a fungal genome defense mechanism known as repeat induced point mutation. This work highlights the role that HGT plays in the evolution of host adaptation in eukaryotic pathogens. It also increases the growing body of evidence indicating that transposons facilitate adaptive HGT events between fungi present in similar environments and hosts.

Most known examples of horizontal gene transfer (HGT) between eukaryotes are ancient. These events are identified primarily using phylogenetic methods on coding regions alone. Only rarely are there examples of HGT where noncoding DNA is also reported. The gene encoding the wheat virulence protein ToxA and the surrounding 14 kb is one of these rare examples. ToxA has been horizontally transferred between three fungal wheat pathogens (Parastagonospora nodorum, Pyrenophora tritici-repentis, and Bipolaris sorokiniana) as part of a conserved ∼14 kb element which contains coding and noncoding regions. Here we used long-read sequencing to define the extent of HGT between these three fungal species. Construction of near-chromosomal-level assemblies enabled identification of terminal inverted repeats on either end of the 14 kb region, typical of a type II DNA transposon. This is the first description of ToxA with complete transposon features, which we call ToxhAT. In all three species, ToxhAT resides in a large (140-to-250 kb) transposon-rich genomic island which is absent in isolates that do not carry the gene (annotated here as toxa⁻). We demonstrate that the horizontal transfer of ToxhAT between P. tritici-repentis and P. nodorum occurred as part of a large (∼80 kb) HGT which is now undergoing extensive decay. In B. sorokiniana, in contrast, ToxhAT and its resident genomic island are mobile within the genome. Together, these data provide insight into the noncoding regions that facilitate HGT between eukaryotes and into the genomic processes which mask the extent of HGT between these species.

Genomic analysis of the aggressive tree pathogen Ceratocystis albifundus

The overall goal of this study was to determine whether the genome of an important plant pathogen in Africa, Ceratocystis albifundus, is structured into subgenomic compartments, and if so, to establish how these compartments are distributed across the genome. For this purpose, the publicly available genome of C. albifundus was complemented with the genome sequences for four additional isolates using the Illumina HiSeq platform. In addition, a reference genome for one of the individuals was assembled using both PacBio and Illumina HiSeq technologies. Our results showed a high degree of synteny between the five genomes, although several regions lacked detectable long-range synteny. These regions were associated with the presence of accessory genes, lower genetic similarity, variation in read-map depth, as well as transposable elements and genes associated with host-pathogen interactions (e.g. effectors and CAZymes). Such patterns are regarded as hallmarks of accelerated evolution, particularly of accessory subgenomic compartments in fungal pathogens. Our findings thus showed that the genome of C. albifundus is made-up of core and accessory subgenomic compartments, which is an important step towards characterizing its pangenome. This study also highlights the value of comparative genomics for understanding mechanisms that may underly and influence the biology and evolution of pathogens.

Phenotypic and Genetic Characterization of Botrytis cinerea Population from Kiwifruit in Sichuan Province, China

Botrytis cinerea (anamorph of Botryotinia fuckeliana) causes gray mold on numerous plants, including kiwifruit. The primary aim of this study was to investigate the phenotypic and genetic characteristics of the Botrytis cinerea population from kiwifruit in Sichuan Province, China. In all, 176 isolates were collected from kiwifruit orchards from eight geographic regions in Sichuan. All isolates were identified as B. cinerea sensu stricto based on the combined datasets, including morphological criteria, determination of the Bc-hch allele, and phylogenetic analysis of the genes RPB2, G3PDH, and HSP60. Three colony types (i.e., sclerotial, mycelial, and conidial) were observed on potato dextrose agar after 2 weeks, with sclerotial isolates, the predominant category, accounting for 40.91%. No obvious differences in microscopic characteristics were observed among the three types. Three genotypes of transposable elements were identified in the B. cinerea population: boty, flipper, and transposa types. The most prevalent genotype from different geographic populations of B. cinerea was transposa; in contrast, the flipper genotype accounted for only 3.98% of the total population, whereas the vacuma genotype was absent. According to MAT locus amplification, 87 and 89 isolates are MAT1-1 and MAT1-2 type, respectively, and the two mating types were found to be balanced overall in the population. Forty-eight representative isolates were all able to cause gray mold to some extent, and disease severities were significantly different between the cultivars Hongyang and Hort16A (P < 0.01). Disease severity was significantly greater on young leaves than on mature leaves (P < 0.01). No significant relationship was found between pathogenicity and geographical region, colony type, or transposon distribution. The results obtained in the present study suggest a relatively uniform species diversity of Botrytis but rich phenotypic and genetic differentiation within the B. cinerea population on kiwifruit in China. Utilizing resistant cultivars and rain-shelter cultivation instead of fungicides may be an effective approach to delaying pathogen variability.

Current Status of Fusarium oxysporum Formae Speciales and Races

The Fusarium oxysporum species complex includes both plant pathogenic and nonpathogenic strains, which are commonly found in soils. F. oxysporum has received considerable attention from plant pathologists for more than a century owing to its broad host range and the economic losses it causes. The narrow host specificity of pathogenic strains has led to the concept of formae speciales, each forma specialis grouping strains with the same host range. Initially restricted to one plant species, this host range was later found to be broader for many formae speciales. In addition, races were identified in some formae speciales, generally with cultivar-level specialization. In 1981, Armstrong and Armstrong listed 79 F. oxysporum formae speciales and mentioned races in 16 of them. Since then, the known host range of F. oxysporum has considerably increased, and many new formae speciales and races have been identified. We carried out a comprehensive search of the literature to propose this review of F. oxysporum formae speciales and races. We recorded 106 well-characterized formae speciales, together with 37 insufficiently documented ones, and updated knowledge on races and host ranges. We also recorded 58 plant species/genera susceptible to F. oxysporum but for which a forma specialis has not been characterized yet. This review raises issues regarding the nomenclature and the description of F. oxysporum formae speciales and races.

Phenotypic instability in fungi

Fungi are prone to phenotypic instability, that is, the vegetative phase of these organisms, be they yeasts or molds, undergoes frequent switching between two or more behaviors, often with different morphologies, but also sometime having different physiologies without any obvious morphological outcome. In the context of industrial utilization of fungi, this can have a negative impact on the maintenance of strains and/or on their productivity. Instabilities have been shown to result from various mechanisms, either genetic or epigenetic. This chapter will review different types of instabilities and discuss some lesser-known ones, mostly in filamentous fungi, while it will direct readers to additional literature in the case of well-known phenomena such as the amyloid prions or fungal senescence. It will present in depth the "white/opaque" switch of Candida albicans and the "crippled growth" degeneration of the model fungus Podospora anserina. These are two of the most thoroughly studied epigenetic phenotypic switches. I will also discuss the "sectors" presented by many filamentous ascomycetes, for which a prion-based model exists but is not demonstrated. Finally, I will also describe intriguing examples of phenotypic instability for which an explanation has yet to be provided.

Epigenetically facilitated mutational assimilation: epigenetics as a hub within the inclusive evolutionary synthesis

After decades of debate about the existence of non‐genetic inheritance, the focus is now slowly shifting towards dissecting its underlying mechanisms. Here, we propose a new mechanism that, by integrating non‐genetic and genetic inheritance, may help build the long‐sought inclusive vision of evolution. After briefly reviewing the wealth of evidence documenting the existence and ubiquity of non‐genetic inheritance in a table, we review the categories of mechanisms of parent–offspring resemblance that underlie inheritance. We then review several lines of argument for the existence of interactions between non‐genetic and genetic components of inheritance, leading to a discussion of the contrasting timescales of action of non‐genetic and genetic inheritance. This raises the question of how the fidelity of the inheritance system can match the rate of environmental variation. This question is central to understanding the role of different inheritance systems in evolution. We then review and interpret evidence indicating the existence of shifts from inheritance systems with low to higher transmission fidelity. Based on results from different research fields we propose a conceptual hypothesis linking genetic and non‐genetic inheritance systems. According to this hypothesis, over the course of generations, shifts among information systems allow gradual matching between the rate of environmental change and the inheritance fidelity of the corresponding response. A striking conclusion from our review is that documented shifts between types of inherited non‐genetic information converge towards epigenetics (i.e. inclusively heritable molecular variation in gene expression without change in DNA sequence). We then interpret the well‐documented mutagenicity of epigenetic marks as potentially generating a final shift from epigenetic to genetic encoding. This sequence of shifts suggests the existence of a relay in inheritance systems from relatively labile ones to gradually more persistent modes of inheritance, a relay that could constitute a new mechanistic basis for the long‐proposed, but still poorly documented, hypothesis of genetic assimilation. A profound difference between the genocentric and the inclusive vision of heredity revealed by the genetic assimilation relay proposed here lies in the fact that a given form of inheritance can affect the rate of change of other inheritance systems. To explore the consequences of such inter‐connection among inheritance systems, we briefly review published theoretical models to build a model of genetic assimilation focusing on the shift in the engraving of environmentally induced phenotypic variation into the DNA sequence. According to this hypothesis, when environmental change remains stable over a sufficient number of generations, the relay among inheritance systems has the potential to generate a form of genetic assimilation. In this hypothesis, epigenetics appears as a hub by which non‐genetically inherited environmentally induced variation in traits can become genetically encoded over generations, in a form of epigenetically facilitated mutational assimilation. Finally, we illustrate some of the major implications of our hypothetical framework, concerning mutation randomness, the central dogma of molecular biology, concepts of inheritance and the curing of inherited disorders, as well as for the emergence of the inclusive evolutionary synthesis.

New insights into Phakopsora pachyrhizi infection based on transcriptome analysis in planta

Asian soybean rust (ASR) is one of the most destructive diseases affecting soybeans. The causative agent of ASR, the fungus Phakopsora pachyrhizi, presents characteristics that make it difficult to study in vitro, limiting our knowledge of plant-pathogen dynamics. Therefore, this work used leaf lesion laser microdissection associated with deep sequencing to determine the pathogen transcriptome during compatible and incompatible interactions with soybean. The 36,350 generated unisequences provided an overview of the main genes and biological pathways that were active in the fungus during the infection cycle. We also identified the most expressed transcripts, including sequences similar to other fungal virulence and signaling proteins. Enriched P. pachyrhizi transcripts in the resistant (PI561356) soybean genotype were related to extracellular matrix organization and metabolic signaling pathways and, among infection structures, in amino acid metabolism and intracellular transport. Unisequences were further grouped into gene families along predicted sequences from 15 other fungi and oomycetes, including rust fungi, allowing the identification of conserved multigenic families, as well as being specific to P. pachyrhizi. The results revealed important biological processes observed in P. pachyrhizi, contributing with information related to fungal biology and, consequently, a better understanding of ASR.

Opposite Polarity Monospore Genome De Novo Sequencing and Comparative Analysis Reveal the Possible Heterothallic Life Cycle of Morchella importuna

Morchella is a popular edible fungus worldwide due to its rich nutrition and unique flavor. Many research efforts were made on the domestication and cultivation of Morchella all over the world. In recent years, the cultivation of Morchella was successfully commercialized in China. However, the biology is not well understood, which restricts the further development of the morel fungus cultivation industry. In this paper, we performed de novo sequencing and assembly of the genomes of two monospores with a different mating type (M04M24 and M04M26) isolated from the commercially cultivated strain M04. Gene annotation and comparative genome analysis were performed to study differences in CAZyme (Carbohydrate-active enzyme) enzyme content, transcription factors, duplicated sequences, structure of mating type sites, and differences at the gene and functional levels between the two monospore strains of M. importuna. Results showed that the de novo assembled haploid M04M24 and M04M26 genomes were 48.98 and 51.07 Mb, respectively. A complete fine physical map of M. importuna was obtained from genome coverage and gene completeness evaluation. A total of 10,852 and 10,902 common genes and 667 and 868 endemic genes were identified from the two monospore strains, respectively. The Gene Ontology (GO) and KAAS (KEGG Automatic Annotation Serve) enrichment analyses showed that the endemic genes performed different functions. The two monospore strains had 99.22% collinearity with each other, accompanied with certain position and rearrangement events. Analysis of complete mating-type loci revealed that the two monospore M. importuna strains contained an independent mating-type structure and remained conserved in sequence and location. The phylogenetic and divergence time of M. importuna was analyzed at the whole-genome level for the first time. The bifurcation time of morel and tuber was estimated to be 201.14 million years ago (Mya); the two monospore strains with a different mating type represented the evolution of different nuclei, and the single copy homologous genes between them were also different due to a genetic differentiation distance about 0.65 Mya. Compared with truffles, M. importuna had an extension of 28 clusters of orthologous genes (COGs) and a contraction of two COGs. The two different polar nuclei with different degrees of contraction and expansion suggested that they might have undergone different evolutionary processes. The different mating-type structures, together with the functional clustering and enrichment analysis results of the endemic genes of the two different polar nuclei, imply that M. importuna might be a heterothallic fungus and the interaction between the endemic genes may be necessary for its complete life history. Studies on the genome of M. importuna facilitate a better understanding of morel biology and evolution.

Multiple independent origins for a subtelomeric locus associated with growth rate in Fusarium circinatum

Fusarium is a diverse assemblage that includes a large number of species of considerable medical and agricultural importance. Not surprisingly, whole genome sequences for many Fusarium species have been published or are in the process of being determined, the availability of which is invaluable for deciphering the genetic basis of key phenotypic traits. Here we investigated the distribution, genic composition, and evolutionary history of a locus potentially determining growth rate in the pitch canker pathogen F. circinatum. We found that the genomic region underlying this locus is highly conserved amongst F. circinatum and its close relatives, except for the presence of a 12 000 base pair insertion in all of the examined isolates of F. circinatum. This insertion encodes for five genes and our phylogenetic analyses revealed that each was most likely acquired through horizontal gene transfer from polyphyletic origins. Our data further showed that this region is located in a region low in G+C content and enriched for repetitive sequences and transposable elements, which is situated near the telomere of Chromosome 3 of F. circinatum. As have been shown for other fungi, these findings thus suggest that the emergence of the unique 12 000 bp region in F. circinatum is linked to the dynamic evolutionary processes associated with subtelomeres that, in turn, have been implicated in the ecological adaptation of fungal pathogens.

Massive Expansion of Gypsy-Like Retrotransposons in Microbotryum Fungi

Transposable elements (TEs) are selfish, autonomously replicating DNA sequences that constitute a major component of eukaryotic genomes and contribute to genome evolution through their movement and amplification. Many fungal genomes, including the anther-smut fungi in the basidiomycete genus Microbotryum, have genome defense mechanisms, such as repeat-induced point mutation (RIP), which hypermutate repetitive DNA and limit TE activity. Little is known about how hypermutation affects the tempo of TE activity and their sequence evolution. Here we report the identification of a massive burst-like expansion of Gypsy-like retrotransposons in a strain of Microbotryum. This TE expansion evidently occurred in the face of RIP-like hypermutation activity. By examining the fitness of individual TE insertion variants, we found that RIP-like mutations impair TE fitness and limit proliferation. Our results provide evidence for a punctuated pattern of TE expansion in a fungal genome, similar to that observed in animals and plants. While targeted hypermutation is often thought of as an effective protection against mobile element activity, our findings suggest that active TEs can persist and undergo selection while they proliferate in genomes that have RIP-like defenses.

Analysis of Transposable Elements in Coccidioides Species

Coccidioides immitis and C. posadasii are primary pathogenic fungi that cause disease in immunologically-normal animals and people. The organism is found exclusively in arid regions of the Southwestern United States, Mexico, and South America, but not in other parts of the world. This study is a detailed analysis of the transposable elements (TE) in Coccidioides spp. As is common in most fungi, Class I and Class II transposons were identified and the LTR Gypsy superfamily is the most common. The minority of Coccidioides Gypsy transposons contained regions highly homologous to polyprotein domains. Phylogenetic analysis of the integrase and reverse transcriptase sequences revealed that many, but not all, of the Gypsy reverse transcriptase and integrase domains clustered by species suggesting extensive transposition after speciation of the two Coccidiodies spp. The TEs were clustered and the distribution is enriched for the ends on contigs. Analysis of gene expression data from C. immitis found that protein-coding genes within 1 kB of hAT or Gypsy TEs were poorly expressed. The expression of C. posadasii genes within 1 kB of Gypsy TEs was also significantly lower compared to all genes but the difference in expression was smaller than C. immitis. C. posadasii orthologs of C. immitis Gyspsy-associated genes were also likely to be TE-associated. In both C. immitis and C. posadasii the TEs were preferentially associated with genes annotated with protein kinase gene ontology terms. These observations suggest that TE may play a role in influencing gene expression in Coccidioides spp. Our hope is that these bioinformatic studies of the potential TE influence on expression and evolution of Coccidioides will prompt the development of testable hypotheses to better understand the role of TEs in the biology and gene regulation of Coccidioides spp.

Comparative genomics reveals cotton-specific virulence factors in flexible genomic regions in Verticillium dahliae and evidence of horizontal gene transfer from Fusarium

Verticillium dahliae isolates are most virulent on the host from which they were originally isolated. Mechanisms underlying these dominant host adaptations are currently unknown. We sequenced the genome of V. dahliae Vd991, which is highly virulent on its original host, cotton, and performed comparisons with the reference genomes of JR2 (from tomato) and VdLs.17 (from lettuce). Pathogenicity-related factor prediction, orthology and multigene family classification, transcriptome analyses, phylogenetic analyses, and pathogenicity experiments were performed. The Vd991 genome harbored several exclusive, lineage-specific (LS) genes within LS regions (LSRs). Deletion mutants of the seven genes within one LSR (G-LSR2) in Vd991 were less virulent only on cotton. Integration of G-LSR2 genes individually into JR2 and VdLs.17 resulted in significantly enhanced virulence on cotton but did not affect virulence on tomato or lettuce. Transcription levels of the seven LS genes in Vd991 were higher during the early stages of cotton infection, as compared with other hosts. Phylogenetic analyses suggested that G-LSR2 was acquired from Fusarium oxysporum f. sp. vasinfectum through horizontal gene transfer. Our results provide evidence that horizontal gene transfer from Fusarium to Vd991 contributed significantly to its adaptation to cotton and may represent a significant mechanism in the evolution of an asexual plant pathogen.

Cut-and-Paste Transposons in Fungi with Diverse Lifestyles

Transposable elements (TEs) shape genomes via recombination and transposition, lead to chromosomal rearrangements, create new gene neighborhoods, and alter gene expression. They play key roles in adaptation either to symbiosis in Amanita genus or to pathogenicity in Pyrenophora tritici-repentis. Despite growing evidence of their importance, the abundance and distribution of mobile elements replicating in a "cut-and-paste" fashion is barely described so far. In order to improve our knowledge on this old and ubiquitous class of transposable elements, 1,730 fungal genomes were scanned using both de novo and homology-based approaches. DNA TEs have been identified across the whole data set and display uneven distribution from both DNA TE classification and fungal taxonomy perspectives. DNA TE content correlates with genome size, which confirms that many transposon families proliferate simultaneously. In contrast, it is independent from intron density, average gene distance and GC content. TE count is associated with species' lifestyle and tends to be elevated in plant symbionts and decreased in animal parasites. Lastly, we found that fungi with both RIP and RNAi systems have more total DNA TE sequences but less elements retaining a functional transposase, what reflects stringent control over transposition.

Deciphering Pathogenicity of Fusarium oxysporum From a Phylogenomics Perspective

Fusarium oxysporum is a large species complex of both plant and human pathogens that attack a diverse array of species in a host-specific manner. Comparative genomic studies have revealed that the host-specific pathogenicity of the F. oxysporum species complex (FOSC) was determined by distinct sets of supernumerary (SP) chromosomes. In contrast to common vertical transfer, where genetic materials are transmitted via cell division, SP chromosomes can be transmitted horizontally between phylogenetic lineages, explaining the polyphyletic nature of the host-specific pathogenicity of the FOSC. The existence of a diverse array of SP chromosomes determines the broad host range of this species complex, while the conserved core genome maintains essential house-keeping functions. Recognition of these SP chromosomes enables the functional and structural compartmentalization of F. oxysporum genomes. In this review, we examine the impact of this group of cross-kingdom pathogens on agricultural productivity and human health. Focusing on the pathogenicity of F. oxysporum in the phylogenomic framework of the genus Fusarium, we elucidate the evolution of pathogenicity within the FOSC. We conclude that a population genomics approach within a clearly defined phylogenomic framework is essential not only for understanding the evolution of the pathogenicity mechanism but also for identifying informative candidates associated with pathogenicity that can be developed as targets in disease management programs.

Comparative genomics of Cryptococcus neoformans var. grubii associated with meningitis in HIV infected and uninfected patients in Vietnam

The vast burden of cryptococcal meningitis occurs in immunosuppressed patients, driven by HIV, and is caused by Cryptococcus neoformansvar. grubii. We previously reported cryptococcal meningitis in Vietnam arising atypically in HIV uninfected, apparently immunocompetent patients, caused by a single amplified fragment length polymorphism (AFLP) cluster of C. neoformansvar. grubii (VNIγ). This variant was less common in HIV infected individuals; it remains unclear why this lineage is associated with apparently immunocompetent patients. To study this host tropism we aimed to further our understanding of clinical phenotype and genomic variation within Vietnamese C. neoformansvar. grubii. After performing MLST on C. neoformans clinical isolates we identified 14 sequence types (STs); ST5 correlated with the VNIγ cluster. We next compared clinical phenotype by lineage and found HIV infected patients with cryptococcal meningitis caused by ST5 organisms were significantly more likely to have lymphadenopathy (11% vs. 4%, p = 0.05 Fisher’s exact test) and higher blood lymphocyte count (median 0.76 versus 0.55 X10⁹ cells/L, p = 0.001, Kruskal-Wallis test). Furthermore, survivors of ST5 infections had evidence of worse disability outcomes at 70 days (72.7% (40/55) in ST5 infections versus 57.1% (52/91) non-ST5 infections (OR 2.11, 95%CI 1.01 to 4.41), p = 0.046). To further investigate the relationship between strain and disease phenotype we performed genome sequencing on eight Vietnamese C. neoformansvar. grubii. Eight genome assemblies exhibited >99% nucleotide sequence identity and we identified 165 kbp of lineage specific to Vietnamese isolates. ST5 genomes harbored several strain specific regions, incorporating 19 annotated coding sequences and eight hypothetical proteins. These regions included a phenolic acid decarboxylase, a DEAD-box ATP-dependent RNA helicase 26, oxoprolinases, a taurine catabolism dioxygenase, a zinc finger protein, membrane transport proteins and various drug transporters. Our work outlines the complexity of genomic pathogenicity in cryptococcal infections and identifies a number of gene candidates that may aid the disaggregation of the pathways associated with the pathogenesis of Cryptococcus neoformansvar. grubii.

Cryptococcal meningitis is a brain infection caused by a yeast, Cryptococcus neoformans, and results in an estimated 600 000 deaths each year. Disease usually only occurs in patients who have some problem with their immune systems—most commonly Human Immunodeficiency Virus (HIV) infection. However, it is increasingly recognized that disease can occur, particularly in southeast and east Asia, in patients with apparently normal immune systems (‘immunocompetent’). We previously showed that almost all infections in immunocompetent patients in Vietnam are due to just one small ‘family’ (or lineage) of Cryptococcus neoformans var. grubii, which we called VNIγ. This is in contrast to disease in HIV infected patients, which can be caused by a number of different families. This suggests that VNIγ strains have an increased ability to cause disease. Here, we define the pattern of disease caused by VNIγ infections compared with other strains in HIV infected patients, and use whole genome sequencing—comparing the entire genetic codes from different strains—to try and understand which genes give the VNIγ family this special ability to cause disease in immunocompetent patients.

Filamentous ascomycete genomes provide insights into Copia retrotransposon diversity in fungi

Background

The relative scarcity of Copia retrotransposons has been recently characterized in metazoans in comparison with the other superfamilies of LTR elements. Furthermore, Copia retrotransposons have often a particular dynamics that results in a highly predominant single clade of elements within a large host taxon, such as the GalEa-like retrotransposons in crustaceans. Taking advantage of the skyrocketing amount of genomic data available for fungi, we carried out the first large-scale comparative genomic analysis of the Copia clades in filamentous ascomycetes.

Results

Screening 30 completely sequenced genomes allowed us to identify more than 2500 Copia copies with conserved LTR, which are distributed in 138 families. Their characterization revealed that fungal Copia diversity is much broader than previously thought with at least 27 clades, 23 of which likely correspond to new ones. While the Copia copy number is low in most species, the two clades GalEa and FunCo1 are widely distributed and highly dominate Copia content as they both account for 80% of the detected sequences.

Conclusions

In Fungi, GalEa retrotransposons are restricted to Pezizomycotina in which they can make up an outstandingly high proportion of the genome (up to 10% in Cenococcum geophilum). At last, we revealed that fungal GalEa elements structurally differ from all other Copia elements with an absence of Primer Binding Site. These elements however harbor a Conserved Hairpin Site which is probably essential for their transposition.

Electronic supplementary material

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

Mapping Transposon Insertions in Bacterial Genomes by Arbitrarily Primed PCR

Transposons can be used to easily generate and label the location of mutations throughout bacterial and other genomes. Transposon insertion mutants may be screened for a phenotype as individual isolates, or by selection applied to a pool of thousands of mutants. Identifying the location of a transposon insertion is critical for connecting phenotype to the genetic lesion. In this unit, we present an easy and detailed approach for mapping transposon insertion sites using arbitrarily-primed PCR (AP-PCR). Two rounds of PCR are used to (1) amplify DNA spanning the transposon insertion junction, and (2) increase the specific yield of transposon insertion junction fragments for sequence analysis. The resulting sequence is mapped to a bacterial genome to identify the site of transposon insertion. In this protocol, AP-PCR as it is routinely used to map sites of transposon insertion within Staphylococcus aureus, is used to illustrate the principle. Guidelines are provided for adapting this protocol for mapping insertions in other bacterial genomes. Mapping transposon insertions using this method is typically achieved in 2 to 3 days if starting from a culture of the transposon insertion mutant. © 2017 by John Wiley & Sons, Inc.

Biology, dynamics, and applications of transposable elements in basidiomycete fungi

The phylum Basidiomycota includes filamentous fungi and yeast species with different ecological and genomic characteristics. Transposable elements (TEs) are abundant components of most eukaryotic genomes, and their transition from being genomic parasites to key drivers of genomic architecture, functionality, and evolution is a subject receiving much attention. In light of the abundant genomic information released during the last decade, the aims of this mini-review are to discuss the dynamics and impact of TEs in basidiomycete fungi. To do this, we surveyed and explored data from 75 genomes, which encompass the phylogenetic diversity of the phylum Basidiomycota. We describe annotation approaches and analyze TE distribution in the context of species phylogeny and genome size. Further, we review the most relevant literature about the role of TEs in species lifestyle, their impact on genome architecture and functionality, and the defense mechanisms evolved to control their proliferation. Finally, we discuss potential applications of TEs that can drive future innovations in fungal research.

Specific PCR Detection of Fusarium oxysporum f. sp. vasinfectum California Race 4 Based on a Unique Tfo1 Insertion Event in the PHO Gene

A highly virulent race 4 (Cal race 4) of Fusarium oxysporum f. sp. vasinfectum was identified in California cotton fields in 2001, and has since been found in increasing numbers of fields. Cal race 4 isolates contain a unique Tfo1 transposon insertion in the PHO gene that was not found in other F. oxysporum f. sp. vasinfectum genotypes. Based on this insertion, a multiplex polymerase chain reaction method was developed to detect the Cal race 4 pathogen. A panel of F. oxysporum f. sp. vasinfectum isolates representing different vegetative compatibility groups (VCG) and DNA sequence types was assembled to test the specificity of the detection method. In all, 16 of 17 Cal race 4 isolates produced a 583-bp amplicon; the other isolate produced a 396-bp amplicon reflecting the absence of the Tfo1 insertion. This isolate was a moderately virulent pathogen among Cal race 4 isolates. In total, 80 other F. oxysporum isolates associated with cotton and 11 other formae speciales of F. oxysporum produced only the 396-bp amplicon. The method also distinguished Cal race 4 isolates from India race 4 isolates and China race 7 isolates, which did not possess the unique Tfo1 insertion but otherwise had identical DNA sequences, and all belong to VCG0114. The method is capable of detecting the pathogen directly from infected stem tissues even before external symptom appears and, thus, provides an effective tool for timely identification of infested fields and seed lots, and should help reduce dissemination of Cal race 4 in the U.S. Cotton Belt.