Monomorphic subtelomeric DNA in the filamentous fungus, Metarhizium anisopliae,contains a RecQ helicase-like gene

A deep dive into the RecQ interactome: something old and something new

RecQ family helicases are found in all domains of life and play roles in multiple processes that underpin genomic integrity. As such, they are often referred to as guardians or caretakers of the genome. Despite their importance, however, there is still much we do not know about their basic functions in vivo, nor do we fully understand how they interact in organisms that encode more than one RecQ family member. We recently took a multi-omics approach to better understand the Saccharomyces cerevisiae Hrq1 helicase and its interaction with Sgs1, with these enzymes being the functional homologs of the disease-linked RECQL4 and BLM helicases, respectively. Using synthetic genetic array analyses, immuno-precipitation coupled to mass spectrometry, and RNA-seq, we found that Hrq1 and Sgs1 likely participate in many pathways outside of the canonical DNA recombination and repair functions for which they are already known. For instance, connections to transcription, ribosome biogenesis, and chromatin/chromosome organization were uncovered. These recent results are briefly detailed with respect to current knowledge in the field, and possible follow-up experiments are suggested. In this way, we hope to gain a wholistic understanding of these RecQ helicases and how their mutation leads to genomic instability.

Telomere length de novo assembly of all 7 chromosomes and mitogenome sequencing of the model entomopathogenic fungus, Metarhizium brunneum, by means of a novel assembly pipeline

Background

More accurate and complete reference genomes have improved understanding of gene function, biology, and evolutionary mechanisms. Hybrid genome assembly approaches leverage benefits of both long, relatively error-prone reads from third-generation sequencing technologies and short, accurate reads from second-generation sequencing technologies, to produce more accurate and contiguous de novo genome assemblies in comparison to using either technology independently. In this study, we present a novel hybrid assembly pipeline that allowed for both mitogenome de novo assembly and telomere length de novo assembly of all 7 chromosomes of the model entomopathogenic fungus, Metarhizium brunneum.

Results

The improved assembly allowed for better ab initio gene prediction and a more BUSCO complete proteome set has been generated in comparison to the eight current NCBI reference Metarhizium spp. genomes. Remarkably, we note that including the mitogenome in ab initio gene prediction training improved overall gene prediction. The assembly was further validated by comparing contig assembly agreement across various assemblers, assessing the assembly performance of each tool. Genomic synteny and orthologous protein clusters were compared between Metarhizium brunneum and three other Hypocreales species with complete genomes, identifying core proteins, and listing orthologous protein clusters shared uniquely between the two entomopathogenic fungal species, so as to further facilitate the understanding of molecular mechanisms underpinning fungal-insect pathogenesis.

Conclusions

The novel assembly pipeline may be used for other haploid fungal species, facilitating the need to produce high-quality reference fungal genomes, leading to better understanding of fungal genomic evolution, chromosome structuring and gene regulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07390-y.

Telomere-to-telomere assembled and centromere annotated genomes of the two main subspecies of the button mushroom Agaricus bisporus reveal especially polymorphic chromosome ends

Agaricus bisporus, the most cultivated edible mushroom worldwide, is represented mainly by the subspecies var. bisporus and var. burnettii. var. bisporus has a secondarily homothallic life cycle with recombination restricted to chromosome ends, while var. burnettii is heterothallic with recombination seemingly equally distributed over the chromosomes. To better understand the relationship between genomic make-up and different lifestyles, we have de novo sequenced a burnettii homokaryon and synchronised gene annotations with updated versions of the published genomes of var. bisporus. The genomes were assembled into telomere-to-telomere chromosomes and a consistent set of gene predictions was generated. The genomes of both subspecies were largely co-linear, and especially the chromosome ends differed in gene model content between the two subspecies. A single large cluster of repeats was found on each chromosome at the same respective position in all strains, harbouring nearly 50% of all repeats and likely representing centromeres. Repeats were all heavily methylated. Finally, a mapping population of var. burnettii confirmed an even distribution of crossovers in meiosis, contrasting the recombination landscape of var. bisporus. The new findings using the exceptionally complete and well annotated genomes of this basidiomycete demonstrate the importance for unravelling genetic components underlying the different life cycles.

A complete annotation of the chromosomes of the cellulase producer Trichoderma reesei provides insights in gene clusters, their expression and reveals genes required for fitness

BACKGROUND

Investigations on a few eukaryotic model organisms showed that many genes are non-randomly distributed on chromosomes. In addition, chromosome ends frequently possess genes that are important for the fitness of the organisms. Trichoderma reesei is an industrial producer of enzymes for food, feed and biorefinery production. Its seven chromosomes have recently been assembled, thus making an investigation of its chromosome architecture possible.

RESULTS

We manually annotated and mapped 9194 ORFs on their respective chromosomes and investigated the clustering of the major gene categories and of genes encoding carbohydrate-active enzymes (CAZymes), and the relationship between clustering and expression. Genes responsible for RNA processing and modification, amino acid metabolism, transcription, translation and ribosomal structure and biogenesis indeed showed loose clustering, but this had no impact on their expression. A third of the genes encoding CAZymes also occurred in loose clusters that also contained a high number of genes encoding small secreted cysteine-rich proteins. Five CAZyme clusters were located less than 50 kb apart from the chromosome ends. These genes exhibited the lowest basal (but not induced) expression level, which correlated with an enrichment of H3K9 methylation in the terminal 50 kb areas indicating gene silencing. No differences were found in the expression of CAZyme genes present in other parts of the chromosomes. The putative subtelomeric areas were also enriched in genes encoding secreted proteases, amino acid permeases, enzyme clusters for polyketide synthases (PKS)-non-ribosomal peptide synthase (NRPS) fusion proteins (PKS-NRPS) and proteins involved in iron scavenging. They were strongly upregulated during conidiation and interaction with other fungi.

CONCLUSIONS

Our findings suggest that gene clustering on the T. reesei chromosomes occurs but generally has no impact on their expression. CAZyme genes, located in subtelomers, however, exhibited a much lower basal expression level. The gene inventory of the subtelomers suggests a major role of competition for nitrogen and iron supported by antibiosis for the fitness of T. reesei. The availability of fully annotated chromosomes will facilitate the use of genetic crossings in identifying still unknown genes responsible for specific traits of T. reesei.

The choanoflagellate Monosiga brevicollis karyotype revealed by the genome sequence: telomere-linked helicase genes resemble those of some fungi

The approximately 42 Mbp assembled genome sequence for the choanoflagellate Monosiga brevicollis reveals that most of the large scaffolds of 300-2,600 kb represent entire chromosomes or chromosome arms. Telomeres are partially assembled at the termini of 37 scaffolds, while another 43 scaffolds end in telomere-associated regions containing distinctive gene sets. Potential centromeric regions were identified on 39 scaffolds. Together, these observations suggest a karyotype of approximately 40 metacentric and submetacentric chromosomes averaging 1 Mbp in size. Genes encoding RecQ family DNA helicases, along with ankyrin-domain proteins and serine/threonine kinases, are associated with most telomeres, a feature shared with some fungi. This telomere-linked helicase gene arrangement might be ancestral to both fungi and choanoflagellates in the super-kingdom Opisthokonta; however, the great lability of telomere architecture suggests that it could also be a convergent feature.

The telomere-linked helicase (TLH) gene family in Magnaporthe oryzae: revised gene structure reveals a novel TLH-specific protein motif

Telomere-linked RecQ helicase (TLH) genes have been identified in several fungi, where they occur as small gene families with each member copy residing within ~10 kb of a telomere. Here we describe the characterization of all 11 TLH gene copies in the reference strain of the fungus Magnaporthe oryzae. A consensus gene prediction revealed that the previously reported TLH1 gene is actually a mutated copy, and the full-length gene is almost two times longer. Only four full-length TLH genes were present in the strain that was analyzed, with the remaining copies containing premature stops caused by point mutations, indels, transposon insertions, and a telomere truncation. Interestingly, all of the TLH gene copies possessed numerous mutations indicative of the action of the repeat-induced point mutation process. However, there was evidence of purifying selection indicating maintenance of gene function. Alignment of full-length proteins from M. oryzae, Schizosaccharomyces pombe and M. anisopliae revealed the presence of a novel, highly conserved protein motif which suggests that the telomere-linked helicases have different functions and/or substrates to the RecQ helicases encoded by "internal" genes.

Role of Blm and collaborating factors in recombination and survival following replication stress in Ustilago maydis

Inactivation of the structural gene for the RecQ family member, BLM in human, Sgs1 in budding yeast, or Rqh1 in fission yeast leads to inappropriate recombination, chromosome abnormalities, and disturbed replication fork progression. Studies with yeasts have demonstrated that auxiliary gene functions can contribute in overlapping ways with Sgs1 or Rqh1 to circumvent or overcome lesions in DNA caused by certain genotoxic agents. In the combined absence of these functions, recombination-mediated processes lead to severe loss of fitness. Here we performed a genetic study to determine the role of the Ustilago maydis Blm homolog in DNA repair and in alleviating replication stress. We characterized the single mutant as well as double mutants additionally deleted of genes encoding Srs2, Fbh1, Mus81, or Exo1. Unlike yeasts, neither the blm srs2, blm exo1, nor blm mus81 double mutant exhibited extreme loss of fitness. Inactivation of Brh2, the BRCA2 homolog, suppressed toxicity to hydroxyurea caused by loss of Blm function. However, differential suppression by Brh2 derivatives lacking the canonical DNA-binding region suggests that the particular domain structure comprising this DNA-binding region may be instrumental in promoting the observed hydroxyurea toxicity.

The Tolypocladium inflatum CPA element encodes a RecQ helicase-like gene

Previously, a repetitive CPA element was discovered in the genome of the filamentous fungus Tolypocladium inflatum; however, no further characterization was technically possible at that time. In this study, PCR amplification was used to detect a 4 kb conserved portion of the CPA element that appeared to be present in most, if not all, genomic CPA elements. The amplicons included a large open reading frame that was most similar to a RecQ helicase-like gene from Metarhizium anisopliae. The repetitive nature of the CPA element suggests that it is related to the eukaryotic Helitron class of transposable elements.

Telomere organization in the ligninolytic basidiomycete Pleurotus ostreatus

Telomeres are structural and functional chromosome regions that are essential for the cell cycle to proceed normally. They are, however, difficult to map genetically and to identify in genome-wide sequence programs because of their structure and repetitive nature. We studied the telomeric and subtelomeric organization in the basidiomycete Pleurotus ostreatus using a combination of molecular and bioinformatics tools that permitted us to determine 19 out of the 22 telomeres expected in this fungus. The telomeric repeating unit in P. ostreatus is TTAGGG, and the numbers of repetitions of this unit range between 25 and 150. The mapping of the telomere restriction fragments to linkage groups 6 and 7 revealed polymorphisms compatible with those observed by pulsed field gel electrophoresis separation of the corresponding chromosomes. The subtelomeric regions in Pleurotus contain genes similar to those described in other eukaryotic systems. The presence of a cluster of laccase genes in chromosome 6 and a bipartite structure containing a Het-related protein and an alcohol dehydrogenase are especially relevant; this bipartite structure is characteristic of the Pezizomycotina fungi Neurospora crassa and Aspergillus terreus. As far as we know, this is the first report describing the presence of such structures in basidiomycetes and the location of a laccase gene cluster in the subtelomeric region, where, among others, species-specific genes allowing the organism to adapt rapidly to the environment usually map.

Characterization of chromosome ends in the filamentous fungus Neurospora crassa

Telomeres and subtelomere regions have vital roles in cellular homeostasis and can facilitate niche adaptation. However, information on telomere/subtelomere structure is still limited to a small number of organisms. Prior to initiation of this project, the Neurospora crassa genome assembly contained only 3 of the 14 telomeres. The missing telomeres were identified through bioinformatic mining of raw sequence data from the genome project and from clones in new cosmid and plasmid libraries. Their chromosomal locations were assigned on the basis of paired-end read information and/or by RFLP mapping. One telomere is attached to the ribosomal repeat array. The remaining chromosome ends have atypical structures in that they lack distinct subtelomere domains or other sequence features that are associated with telomeres in other organisms. Many of the chromosome ends terminate in highly AT-rich sequences that appear to be products of repeat-induced point mutation, although most are not currently repeated sequences. Several chromosome termini in the standard Oak Ridge wild-type strain were compared to their counterparts in an exotic wild type, Mauriceville. This revealed that the sequences immediately adjacent to the telomeres are usually genome specific. Finally, despite the absence of many features typically found in the telomere regions of other organisms, the Neurospora chromosome termini still retain the dynamic nature that is characteristic of chromosome ends.

Telomeres in the rice blast fungus Magnaporthe oryzae: the world of the end as we know it

The subtelomeres of many microbial eukaryotes are highly enriched in genes with roles in niche adaptation. Host and cultivar specificity genes in the rice blast fungus Magnaporthe oryzae also tend to be located near telomeres. In addition, the M. oryzae telomeres are highly variable chromosome regions. These observations suggested that plant pathogenic fungi might also use subtelomere regions for the amplification of genes with adaptive significance. Targeted sequencing of the M. oryzae telomeres provided an opportunity to test this hypothesis, and has yielded valuable insights into the organization and dynamics of these important chromosome regions.

Organization of chromosome ends in the rice blast fungus, Magnaporthe oryzae

Eukaryotic pathogens of humans often evade the immune system by switching the expression of surface proteins encoded by subtelomeric gene families. To determine if plant pathogenic fungi use a similar mechanism to avoid host defenses, we sequenced the 14 chromosome ends of the rice blast pathogen, Magnaporthe oryzae. One telomere is directly joined to ribosomal RNA-encoding genes, at the end of the ∼2 Mb rDNA array. Two are attached to chromosome-unique sequences, and the remainder adjoin a distinct subtelomere region, consisting of a telomere-linked RecQ-helicase (TLH) gene flanked by several blocks of tandem repeats. Unlike other microbes, M.oryzae exhibits very little gene amplification in the subtelomere regions—out of 261 predicted genes found within 100 kb of the telomeres, only four were present at more than one chromosome end. Therefore, it seems unlikely that M.oryzae uses switching mechanisms to evade host defenses. Instead, the M.oryzae telomeres have undergone frequent terminal truncation, and there is evidence of extensive ectopic recombination among transposons in these regions. We propose that the M.oryzae chromosome termini play more subtle roles in host adaptation by promoting the loss of terminally-positioned genes that tend to trigger host defenses.

The RecQ gene family in plants

RecQ helicases are conserved throughout all kingdoms of life regarding their overall structure and function. They are 3'-5' DNA helicases resolving different recombinogenic DNA structures. The RecQ helicases are key factors in a number of DNA repair and recombination pathways involved in the maintenance of genome integrity. In eukaryotes the number of RecQ genes and the structure of RecQ proteins vary strongly between organisms. Therefore, they have been named RecQ-like genes. Knockouts of several RecQ-like genes cause severe diseases in animals or harmful cellular phenotypes in yeast. Until now the largest number of RecQ-like genes per organism has been found in plants. Arabidopsis and rice possess seven different RecQ-like genes each. In the almost completely sequenced genome of the moss Physcomitrella patens at least five RecQ-like genes are present. One of the major present and future research aims is to define putative plant-specific functions and to assign their roles in DNA repair and recombination pathways in relation to RecQ genes from other eukaryotes. Regarding their intron positions, the structures of six RecQ-like genes of dicots and monocots are virtually identical indicating a conservation over a time scale of 150 million years. In contrast to other eukaryotes one gene (RecQsim) exists exclusively in plants. It possesses an interrupted helicase domain but nevertheless seems to have maintained the RecQ function. Owing to a recent gene duplication besides the AtRecQl4A gene an additional RecQ-like gene (AtRecQl4B) exists in the Brassicaceae only. Genetic studies indicate that a AtRecQl4A knockout results in sensitivity to mutagens as well as an hyper-recombination phenotype. Since AtRecQl4B was still present, both genes must have non-redundant roles. Analysis of plant RecQ-like genes will not only increase the knowledge on DNA repair and recombination, but also on the evolution and radiation of protein families.

Evolutionary-conserved telomere-linked helicase genes of fission yeast are repressed by silencing factors, RNAi components and the telomere-binding protein Taz1

In Schizosaccharomyces pombe the RNAi machinery and proteins mediating heterochromatin formation regulate the transcription of non-coding centromeric repeats. These repeats share a high sequence similarity with telomere-linked helicase (tlh) genes, implying an ancestral relationship between the two types of elements and suggesting that transcription of the tlh genes might be regulated by the same factors as centromeric repeats. Indeed, we found that mutants lacking the histone methyltransferase Clr4, the Pcu4 cullin, Clr7 or Clr8, accumulate high levels of tlh forward and reverse transcripts. Mutations and conditions perturbing histone acetylation had similar effects further demonstrating that the tlh genes are normally repressed by heterochromatin. In contrast, mutations in the RNAi factors Dcr1, Ago1 or Rdp1 led only to a modest derepression of the tlh genes indicating an alternate pathway recruits heterochromatin components to telomeres. The telomere-binding protein Taz1 might be part of such a redundant pathway, tlh transcripts being present at low levels in Deltataz1 mutants and at higher levels in Deltataz1 Deltadcr1 double mutants. Surprisingly, the chromodomain protein Chp1, a component of the Ago1-containing RITS complex, contributes more to tlh repression than Ago1, indicating the repressive effects of Chp1 are partially independent of RITS. The tlh genes are found in the subtelomeric regions of several other fungi raising the intriguing possibility of conserved regulation and function.

The dawn of fungal pathogen genomics

Recent advances in sequencing technologies have led to a remarkable increase in the number of sequenced fungal genomes. Several important plant pathogenic fungi are among those that have been sequenced or are being sequenced. Additional fungal pathogens are likely to be sequenced in the near future. Analysis of the available genomes has provided useful information about genes that may be important for plant infection and colonization. Genome features, such as repetitive sequences, telomeres, conserved syntenic blocks, and expansion of pathogenicity-related genes, are discussed in detail with Magnaporthe oryzae (M. grisea) and Fusarium graminearum as examples. Functional and comparative genomic studies in plant pathogenic fungi, although still in the early stages and limited to a few pathogens, have enormous potential to improve our understanding of the molecular mechanisms involved in host-pathogen interactions. Development of advanced genomics tools and infrastructure is critical for efficient utilization of the vast wealth of available genome sequence information and will form a solid foundation for systems biology studies of plant pathogenic fungi.