Chromosome-Level Assembly Reveals a Fifteen-Chromosome Aneuploid Genome and Environmental Adaptation Strategy of Chinese Traditional Medical Fungus Wolfiporia hoelen

The sclerotia of Wolfiporia hoelen are one of the most important traditional Chinese medicines and foods commonly used in China, Japan, Korea, and other Asian countries. To provide a high-quality reference genome and deepen our understanding of the genome of W. hoelen to elucidate various biological phenomena. In this study, we assembled three genomes of W. hoelen using a combination of Nanopore and Illumina sequencing strategies. The fifteen-chromosome genome L7 of W. hoelen was assembled with two-sided telomere and rDNA sequences for the first time. The chromosome count was subsequently confirmed through collinearity analysis, correcting the previous belief that W. hoelen had only fourteen chromosomes. Moreover, the aneuploid genome was discovered in W. hoelen for the first time through sequencing depth analysis of different chromosomes, and only some strains of W. hoelen exhibit aneuploid genomes. According to the genome analysis of homokaryotic offspring and protoplast-isolated strains, a potential variation in chromosome allocation patterns was revealed. Moreover, the gene function enrichment analysis of genes on reduplicated chromosomes demonstrated that aneuploidy in the genome may be the result of environmental adaptation for W. hoelen. The discovery of an aneuploid genome also provides new ideas for genetic improvement of W. hoelen.


Introduction
Genomics is the core and foundation of life science research and the pioneering direction of future scientific and technological innovation.With the advancement of sequencing and assembly technologies, numerous genomes of edible and medicinal macrofungi have been assembled and published.Not only widely cultivated mushrooms like Lentinula edodes [1], Auricularia heimuer [2], Agaricus bisporus [3], and others [4,5], but also some newly domesticated species such as Oudemansiella raphanipes [6] and Naematelia aurantialba [7] have had their genomes published.Meanwhile, the genome assembly level has risen significantly.Some species have been assembled to the chromosome or near-chromosome level [8,9], and some have been assembled to the T2T level [10].Meanwhile, the chromosome number of more and more species has been revealed [11].Unlike plants and animals, the chromosome number in fungi is difficult to observe due to their small size [12].For the study of chromosome numbers, genetic linkage maps, macroscopic observations, pulsed-field gel electrophoresis, and high-throughput chromosome conformation capture can all be used [13][14][15][16][17][18].
Eukaryotes organize their genomes into chromosomes, with each organism having its own specific karyotype, originally defined by the number and appearance of the chromosomes in the nucleus [19].Euploidy refers to a chromosome number that is an integral multiple of the chromosome number for a haploid cell [20].Aneuploidy is a deviation from a balanced genome caused by either the gain or loss of chromosomes, leading to an unbalanced genomic state [21].The aneuploid phenomenon was discovered in plants, animals, and also fungi [22][23][24].The yeast fungi Saccharomyces cerevisiae, Candida albicans, and Cryptococcus neoformans have all been demonstrated to have widely aneuploid genomes [25][26][27], which can be artificially induced under adverse environmental conditions [28][29][30].Recently, the aneuploid phenomenon was discovered in the dimorphic edible mushroom Tremella fuciformis, which includes yeast and mycelial forms.This phenomenon was demonstrated to be related to the scarce formation of basidiospores [31].
Wofiporia hoelen (Fr.)Y.C. Dai and V. Papp [32,33] is a traditional Chinese medicine known as 'Fuling' in China.It grows underground on the roots of pines and is distributed in East Asia [34].The sclerotia of W. hoelen are edible and medicinal tissues widely used as traditional crude drugs in China, South Korea, and Japan.They can also be consumed as food [35,36].With the advancement of pharmacy, various pharmacological activities of the asexual sclerotia of W. hoelen have been demonstrated, including diuresis, tranquilization [37], antioxidant effects [38], antitumor properties [39], anti-inflammatory effects [40], and immunomodulatory effects [41].
As an important and unique species capable of forming sclerotia.To further reveal the evolution and formation mechanism of sclerotia, and mechanisms of function metabolite synthesis, its genome was published in 2020 for the first time [42].Subsequently, Cao et al. [43], Kim et al. [44], and Zhang et al. [45] successively published the genome of W. hoelen.Meanwhile, we also assembled and published the homokaryotic genome SS20 of W. hoelen based on the establishment of homokaryotic strain isolation and identification methods [46,47].Strain SS20 is a homokaryotic strain, and it was also the first homozygous genome assembled in our previous study [8].In addition, Cao et al. first proposed the 14 chromosomes of W. hoelen based on microscopic observations [43].And we further agreed with the results based on the telomere numbers and Hi-C-assisted assembly.This is the first chromosome-level genome of W. hoelen.Even though the basic information of the W. hoelen genome has been revealed, the relation between genome and biological phenomenon is still little known, such as the different strains have distinct homokaryotic offspring types [48], and hybrid strains often demonstrate poorer growth performance compared to the parent strains.
In the present study, we assembled three different genomes of W. hoelen using a combination of Nanopore long-read sequencing and Illumina short-read sequencing strategies to achieve chromosomal-level genomes.The complete chromosomal-level genome of W. hoelen with 15 chromosomes was assembled here.This study corrected previous findings that W. hoelen has 14 chromosomes.The aneuploid genomes of W. hoelen were discovered, and the aneuploid phenomenon of various W. hoelen strains and homokaryotic offspring was analyzed.This study revealed a potential variation in chromosome allocation patterns.Gene function analysis of genes on reduplicated chromosomes revealed the possible reasons for the formation of aneuploid genomes in W. hoelen.This study will deepen the understanding of genomics, genetics, and the evolution of macrofungi, and provide new insights for the genetic improvement of W. hoelen.

Chromosome-Level Genome Assembly of Wolfiporia hoelen Strain L7
To further assemble a high-quality genome of W. hoelen based on a previous study [8], the homokaryotic strain L7 was utilized for genome sequencing.The sequencing was performed using the Illumina NovaSeq 6000 platform and the PromethION platform.Homokaryotic strain L7 was isolated from strain L12, which was derived from a distinct strain with a previously published SS20 genome [8].A de novo genome assembly was performed by integrating Nanopore long-read sequencing and Illumina short-read sequencing (Figure 1).The predicted genome size is 57.13 Mb.There was no apparent heterozygous peak, and the heterozygosity was low at 0.106% based on the k-mer (k = 17) analysis (Figure 2(B1)).
Furthermore, the telomere and rDNA analysis showed 14 contigs of genome L7 containing double-sided telomeres, and 1 contig containing single-sided telomeres and rDNA, indicating that the genome L7 has been assembled to the chromosome level (Figure 1).In total, 96.98% of sequences were deemed as chromosomes (57.88 Mb), with the longest chromosome 5.46 Mb and the shortest 2.15 Mb.A genome of 59.68 Mb was assembled based on 294,029 Nanopore reads (~85×, 5.08 Gb data size), which consists of 28 contigs with a contig N50 ~4.11 Mb, with the longest contig being 5.46 Mb.The mapping rate of Illumina data is 99.83%, the genome coverage is 99.94%, and 97.9% of complete BUSCOs all indicate the high genome integrity of L7.Furthermore, the telomere and rDNA analysis showed 14 contigs of genome L7 containing double-sided telomeres, and 1 contig containing single-sided telomeres and rDNA, indicating that the genome L7 has been assembled to the chromosome level (Figure 1).In total, 96.98% of sequences were deemed as chromosomes (57.88 Mb), with the longest chromosome 5.46 Mb and the shortest 2.15 Mb.

Genome Comparison and Collinearity Analysis of Different Wolfioria hoelen Strains
Meanwhile, the genomes of Pr2C and PrT were also sequenced and assembled, which were all homokaryotic strains and have low genome heterozygosity (Figure 2  This study This study This study [8] When comparing the three genomes, it is important to consider factors such as genome size, contig N50, percentage of repeat sequences, genome coverage, and the presence of complete BUSCOs.These parameters provide insights into the accuracy of genome assembly.The mapping ratio of Illumina data for L7 is higher than that of Pr2C and PrT, indicating the high quality of the L7 genome.Furthermore, 13 contigs with doublesided telomeres, 1 contig with single-sided telomere and single-sided rDNA, and 2 contigs with single-sided telomeres were discovered in the Pr2C genome, while 11 contigs with double-sided telomeres and 5 contigs with single-sided telomeres were identified in the PrT genome (Table S1).In addition, the distribution of sequencing depth in genome L7 showed high consistency, but in the genomes Pr2C and PrT, the sudden fluctuations in sequencing depth were sporadic (Figure 2).This also indicates the high-quality assembly of genome L7.
Compared to the SS20 genome we published previously as the highest quality genome [8], the number of contigs has significantly decreased from 78 to 27/28 contigs.The contig N50 was significantly extended from 3760 kb to over 4100 kb (Table 1).The complete chromosomes contain telomeres on both ends.Genome L7 (15 chromosomes) has more chromosomes than genome Pr2C (13 chromosomes), which has more chromosomes than genome PrT (11 chromosomes), and more chromosomes than genome SS20 (9 chromosomes) (Table S1).This demonstrates the higher quality of the newly assembled genomes compared to genome SS20.Meanwhile, the longest chromosome of the SS20 genome, anchored according to Hi-C data, was 7.47 Mb [8], significantly larger than the largest chromosome found here.This suggests that chromosome 1 of SS20 may have been incorrectly assembled.
Moreover, the collinearity analysis was conducted for the three genomes based on protein identity (Figure 3).The majority of chromosomes in L7 exhibit high collinearity with the contigs of Pr2C and PrT.In total, 13 contigs of the Pr2C genome exhibit high collinearity with the L7 genome, and 10 contigs of the PrT genome show high collinearity with the L7 genome.This indicates that the corresponding contigs have been assembled at the chromosome level.Meanwhile, some contigs exhibit high collinearity with genome L7 but lack a single-sided telomere, such as contig3 of genome Pr2C, and contig9 and contig12 of genome PrT.The most challenging chromosome to assemble is Chr15, as the corresponding contigs of the Pr2C and PrT genomes were scattered.Contig15, contig17, and contig 20 of the PrT genome, as well as contig 16 and contig 15 of the Pr2C genome, exhibit high collinearity with Chr15 in the L7 genome.The collinearity analysis revealed no breaks in the chromosomes of genome L7.This suggests that W. hoelen has 15 chromosomes, which were exclusively assembled using the OLC (Overlap-Layout-Consensus) algorithm.While there are some assembly mistakes, contig1 of the PrT genome has three telomeres, and contig16 of Pr2C is significantly shorter than the corresponding Chr15, and has doublesided telomeres.The 15 chromosomes of W. hoelen revealed here correct the previous results of 14 chromosomes [8].This study This study This study [8] Moreover, the collinearity analysis was conducted for the three genomes based on protein identity (Figure 3).The majority of chromosomes in L7 exhibit high collinearity with the contigs of Pr2C and PrT.In total, 13 contigs of the Pr2C genome exhibit high collinearity with the L7 genome, and 10 contigs of the PrT genome show high collinearity with the L7 genome.This indicates that the corresponding contigs have been assembled at the chromosome level.Meanwhile, some contigs exhibit high collinearity with genome L7 but lack a single-sided telomere, such as contig3 of genome Pr2C, and contig9 and contig12 of genome PrT.The most challenging chromosome to assemble is Chr15, as the corresponding contigs of the Pr2C and PrT genomes were scattered.Contig15, contig17, and contig 20 of the PrT genome, as well as contig 16 and contig 15 of the Pr2C genome, exhibit high collinearity with Chr15 in the L7 genome.The collinearity analysis revealed no breaks in the chromosomes of genome L7.This suggests that W. hoelen has 15 chromosomes, which were exclusively assembled using the OLC (Overlap-Layout-Consensus) algorithm.While there are some assembly mistakes, contig1 of the PrT genome has three telomeres, and contig16 of Pr2C is significantly shorter than the corresponding Chr15, and has double-sided telomeres.The 15 chromosomes of W. hoelen revealed here correct the previous results of 14 chromosomes [8].

Aneuploidy Genome of Wolfiporia hoelen
Based on determining the number of chromosomes, the ploidy of the W. hoelen genome was further studied.Sequencing depth of reads can reflect the phenomenon of aneuploid genomes.The sequencing depth of genomes L7, Pr2C, and PrT was analyzed based on Nanopore sequencing reads (Figure 2).The obvious aneuploidy phenomenon was discovered in W. hoelen.Chromosome 15 of genome L7 has 2× (fold) sequencing depth compared to other chromosomes.Chromosomes 14 and 15 of genome Pr2C have 1.5× and 2× sequencing depth, respectively, compared to other chromosomes.Chromosome 1 of genome PrT has an additional 1.5× sequencing depth compared to other chromosomes based on genome Pr2C.This is the first discovery of genomic aneuploidy in W. hoelen.

Aneuploidy Genome of Wolfiporia hoelen
Based on determining the number of chromosomes, the ploidy of the W. hoelen genome was further studied.Sequencing depth of reads can reflect the phenomenon of aneuploid genomes.The sequencing depth of genomes L7, Pr2C, and PrT was analyzed based on Nanopore sequencing reads (Figure 2).The obvious aneuploidy phenomenon was discovered in W. hoelen.Chromosome 15 of genome L7 has 2× (fold) sequencing depth compared to other chromosomes.Chromosomes 14 and 15 of genome Pr2C have 1.5× and 2× sequencing depth, respectively, compared to other chromosomes.Chromosome 1 of genome PrT has an additional 1.5× sequencing depth compared to other chromosomes based on genome Pr2C.This is the first discovery of genomic aneuploidy in W. hoelen.
Further, the sequencing depth of genomes L7, Pr2C, and PrT was analyzed based on Illumina sequencing reads.When comparing the sequencing depth analyzed with Nanopore reads, the sequencing depth analyzed with Illumina reads exhibited less consistency, with more abrupt fluctuations between high and low sequencing depth regions.While the Illumina sequencing depth can also reflect the aneuploidy phenomenon and yield similar results (Figure S1).
Due to the strain L7 being isolated with a single spore, the aneuploid genome of L7 may have been caused by abnormal chromosome allocation during meiosis or by the parent strain L12 also having an aneuploid genome.To strains Pr2C and PrT, which were isolated with protoplasts, the aneuploidy phenomenon reflects that the original strain L14 also had an aneuploid genome.Therefore, some parental strains were used to conduct aneuploidy genome analysis.

Aneuploidy Characters for Different-Type Heterokaryotic Strains
To further reveal the aneuploidy characteristics in species W. hoelen, the heterokaryotic strains isolated from cultivated or wild sclerotia were sequenced.The sequencing depth of different chromosomes was analyzed, and genome L7 was used as the reference genome.
In the six sequenced strains, three strains (L2, L12, L14) exhibit an aneuploidy phenomenon, where the duplicated chromosomes are all from chromosome 15.The other three strains have euploid genomes (Figure 4).Chromosome 15 of genomes L2 and L12 has 1.5× sequencing depth compared to other chromosomes, while chromosome 15 of genome L14 has 3× sequencing depth compared to other chromosomes (Figure 4).In the species W. hoelen, some strains have euploid genomes, while others have aneuploid genomes.
While the Illumina sequencing depth can also reflect the aneuploidy phenomenon and yield similar results (Figure S1).
Due to the strain L7 being isolated with a single spore, the aneuploid genome of L7 may have been caused by abnormal chromosome allocation during meiosis or by the parent strain L12 also having an aneuploid genome.To strains Pr2C and PrT, which were isolated with protoplasts, the aneuploidy phenomenon reflects that the original strain L14 also had an aneuploid genome.Therefore, some parental strains were used to conduct aneuploidy genome analysis.

Aneuploidy Characters for Different-Type Heterokaryotic Strains
To further reveal the aneuploidy characteristics in species W. hoelen, the heterokaryotic strains isolated from cultivated or wild sclerotia were sequenced.The sequencing depth of different chromosomes was analyzed, and genome L7 was used as the reference genome.
In the six sequenced strains, three strains (L2, L12, L14) exhibit an aneuploidy phenomenon, where the duplicated chromosomes are all from chromosome 15.The other three strains have euploid genomes (Figure 4).Chromosome 15 of genomes L2 and L12 has 1.5× sequencing depth compared to other chromosomes, while chromosome 15 of genome L14 has 3× sequencing depth compared to other chromosomes (Figure 4).In the species W. hoelen, some strains have euploid genomes, while others have aneuploid genomes.In our previous results, the strains can be divided into two types, each with a different homokaryotic strain type [49].For Type I strains, the homokaryotic strain was significantly different from its parent strain, showing a slower growth rate and less aerial mycelia.The homokaryotic strain of Type II strains shows no significant difference from its parent strain.It grows rapidly and has thick aerial mycelia.Among the six strains used in this study, L2, L4, and L14 were type I strains, while L11, L12, and L43 were type II strains.Notably, L2, L12, and L14 exhibited aneuploid genomes, whereas L4, L11, and L43 exhibited euploid genomes (Figure 4).This indicates that the strain types are not directly related to euploid or aneuploid genomes.

Aneuploidy Characters for Homokaryotic Offspring of L12 and L14
To determine if chromosome allocation is normal during meiosis, the genomes of homokaryotic offspring were sequenced.The sequencing depth of different chromosomes was analyzed.For homokaryotic offspring of strain L12, the genome L7 was used as the reference genome, and the genome Pr2C was used as the reference genome for homokaryotic offspring of strain L14.
Eight homokaryotic offspring of strain L12 were further sequenced and analyzed.Among them, three of which have 2× sequencing depth for chromosome 15 compared to other chromosomes, while the other five strains exhibited a euploid genome with similar sequencing depth across all chromosomes (Figure 5).Meanwhile, twenty homokaryotic offspring of strain L14 were further sequenced and analyzed.Among them, fifteen strains exhibited aneuploid genomes, while five strains showed euploid genomes.Among the fifteen offspring strains with aneuploid genomes, some strains only exhibit one chromosome duplication.Specifically, ten strains have 2× sequencing depth for chromosome 15 compared to other chromosomes, four strains have 3× sequencing depth for chromosome 15 compared to other chromosomes, one strain has 4× sequencing depth for chromosome 15 compared to other chromosomes, and one strain has 2× sequencing depth for chromosome 13 compared to other chromosomes.Moreover, some strains have chromosome duplication for more than one chromosome.For strain S10, chromosome 14 has a sequencing depth of 1.5×, and chromosome 15 has a sequencing depth of 2×.For strain S13, chromosome 13 has a sequencing depth of 1.5×, and chromosome 15 has a sequencing depth of 2×.Lastly, for strain S15, chromosome 11 has a sequencing depth of 1.5×, and chromosome 15 has a sequencing depth of 2× (Figure 6).
According to the law of normal distribution of chromosomes.If the genome of L12 contains two types of nuclei, one containing a euploid genome (15 chromosomes), and the other containing an aneuploid genome (15 chromosomes + Chr15).The homokaryotic offspring would only include strains with a euploid genome (15 chromosomes) or with an aneuploid genome (15 chromosomes + Chr15).This was in accordance with the results of genome analysis for nine homokaryotic offspring of strain L12, among which five have euploid genomes, and four have aneuploid genomes (including strain L7).Whereas strain L14 exhibited reduplication only on chromosome 15, the homokaryotic offspring showed reduplication on chromosomes 11, 13, and 14 except chromosome 15.That presents the strain L14 may have a disordered chromosome distribution.That may be related to different types of W. hoelen strains [48].In general, if a chromosome were to be reduplicated, the corresponding nucleus would lack a chromosome, leading to death due to chromosome loss.That may be the reason why strains with absent chromosomes were not obtained.
However, for the strains Pr2C and PrT, isolated with protoplasts, chromosome 14 was reduplicated, and in PrT, chromosome 1 was also reduplicated, except for chromosome 15 (Figure 2(A2,A3)).The homokaryotic strains isolated with protoplasts may have more than one nucleus with the same mating type.The results showed that chromosomes 1 and 14 may also be reduplicated in some nuclei.Even though the L14 strain only exhibited reduplication of chromosome 15, the various nuclei in strain L15 may possess different genome types, resulting in diverse genomes with varying chromosome reduplication patterns.To reveal the actual chromosome distribution law and genome variations among different nuclei of the same strain, conducting additional genome analyses on strains isolated with protoplasts may provide the answers in follow-up studies.
than one nucleus with the same mating type.The results showed that chromosomes 1 and 14 may also be reduplicated in some nuclei.Even though the L14 strain only exhibited reduplication of chromosome 15, the various nuclei in strain L15 may possess differen genome types, resulting in diverse genomes with varying chromosome reduplication pat terns.To reveal the actual chromosome distribution law and genome variations among different nuclei of the same strain, conducting additional genome analyses on strains iso lated with protoplasts may provide the answers in follow-up studies.

Gene Function Enrichment of Genes on Chromosome 15
In the aneuploid genome, the gene dosage effects were prominent [21].In the isolated strains from wild or cultivated sclerotia, only chromosome 15 was duplicated, albeit with

Gene Function Enrichment of Genes on Chromosome 15
In the aneuploid genome, the gene dosage effects were prominent [21].In the isolated strains from wild or cultivated sclerotia, only chromosome 15 was duplicated, albeit with varying multiples.Gene function enrichment was conducted for genes located on chromosome 15.On the whole, 54 genes are enriched in the 'biological process', 120 genes are enriched in the 'molecular function', and 107 genes are enriched in the 'cellular components'.For biological processes, the majority of genes are enriched in metabolic process (16), cellular process (15), and single-organism process (10).For molecular function, the majority of genes are enriched in catalytic activity (62) and binding (52).For cellular components, the majority of genes are enriched in the cell (19), membrane (17), and organelle (13) (Figure S2).
According to the results, the enriched Gene Ontology (GO) terms of genes on chromosome 15 are mainly associated with environmental stress responses, such as salt stress, high or low-temperature stress, drought stress, heavy metal stress, and pathogen infection (Table 2).Moreover, chitin being one of the cell wall components, chitin biosynthesis was also enriched (Figure 7).As the cell wall is closely related to the outer environment, chitin biosynthesis must play an important role in resisting adverse environments.The study demonstrated that chromosome reduplication in W. hoelen primarily contributes to environmental adaptation.On the one hand, W. hoelen is widely distributed or grown from the northeast to the south of China, and it can adapt to various environments.On the other hand, the sclerotia grow underground, coming into contact and interacting with numerous other microorganisms.The chromosome reduplication and aneuploid genome formation of W. hoelen may indicate an adaptation to the environments.
In our previous results, chromosome 14 of the SS20 genome was found to play an important role in adapting to harsh environments.The genome features a high ratio of transposons, low gene density, low GC content, absence of single-copy homologous genes, and a collection of resistance genes [8].Because Chr15 was considered to have similar functions, the gene density, repetitive sequence ratio, and GC content were analyzed.It was found that chromosome 15 has a high ratio of repetitive sequences, low gene density, low GC content, and a concentration of resistance genes (Table S2).The reduplicated chromosome 15 was identified as the homologous chromosome of chromosome 14 in the SS20 genome, which was classified as an accessory chromosome in the species W. hoelen [8].This further demonstrated the relationship between chromosome 15 and environmental adaptation.Table 2. GO enrichment of genes on chromosome 15 and related response process.

GO Term and Related Response Biological Process
Reference GO: 0003855 3-dehydroquinate dehydratase activity

Homokaryotic Strain Isolation
The strains CGMCC 5.545 (L14), CGMCC 5.55 (L12), and NBRC 30268 (L11) were kindly provided by Hubei Provincial Hospital of Traditional Chinese Medicine.Strains CGMCC 5.137 (L2), L12, and L14 were preserved at the General Microbiological Culture Collection Center (CGMCC).Strain L4 was isolated from sclerotia cultivated in Yunnan Province, while strain L43 was isolated from wild sclerotia collected from Hubei Province.The homokaryotic strain was isolated using our previous methods [47,85] and confirmed through specific polymorphic loci of the rpb2 gene [47].Homokaryotic strains Pr2C and PrT were isolated using protoplasts from strain CGMCC 5.545, while other homokaryotic strains were obtained through single-spore isolation.All the strains were maintained on potato dextrose agar medium (PDA: 200 g/L potato, 20 g/L glucose, 18 g/L agar) at 4 • C for future experiments.

Genome Sequencing, Assembly, and Quality Assessment
The mycelia sample was prepared using our previous methods [8].The total genomic DNA was extracted using the conventional cetyltrimethylammonium bromide method, and its quality was assessed using the Nanopore One spectrophotometer (NanoDrop Technology, Wilmington, DE, USA) and the Qubit 3.0 Fluorometer (Life Technologies, Carlsbad, CA, USA).For the Illumina and Oxford Nanopore library preparation, the VAHTS ® Universal Plus DNA Library Prep Kit for Illumina v2 and the SQK-LSK110, EXP-NBD104/114 kits were used separately.The library preparation followed standard protocols for each kit.The Illumina NovaSeq 6000 platform and PromethION platforms (Oxford Nanopore Technologies, Oxford, UK) were used for genome sequencing at Wuhan Benagen Tech Solutions Company Limited (Benagen, Wuhan, China).
Raw sequencing data generated by the PromethION platform underwent quality assessment and were processed using Oxford Nanopore GUPPY (version 0.3.0) to eliminate unsuccessful reads.The remaining passed reads were used for further analysis.For genome assembly, the software NECAT (version v0.0.1) (https://github.com/xiaochuanle/NECAT)(accessed on 10 August 2024) was used for preliminary genome assembly, and Racon (version 1.4.11) was utilized for error correction based on ONT sequencing data for two rounds.Further, the software Pilon (version 1.23) was used for error correction based on Illumina sequencing data for two rounds [86].The Circos graph of the genome was generated using the 'Advanced Circos' module of TBtools v 2.086 [87].The heterozygosity of the genome was estimated based on analysis of k-mer depth distribution using Jellyfish5 [88] and GenomeScope6 [89] for Illumina data.The integrity of chromosome-level genome assembly was estimated by BUSCO [90].

Annotation of Gene Structure and Function
The combination of ab initio prediction, homology-based prediction, and transcriptomeassisted prediction was used to identify the genes encoding proteins.The software Exonerate (version 2.4.0) was used for homology-based prediction, Augustus (version 3.

Whole-Genome Collinearity Analysis
To explore the number of chromosomes, whole-genome collinearity was conducted among the three genomes of W. hoelen sequenced and assembled in this study: Pr2C, PrT, and L7.The 'One Step MCScanX' and 'Multiple Synteny Plot' commands in the 'Comparative Genomics module' of the TBtools software (v 2.086) were used for the analysis and visualization of collinearity results [87].

Genome Resequencing
The sample preparation methods are the same as those in our previous study [46].A total of 34 strains were subjected to genome resequencing, including 6 heterokaryotic strains isolated from sclerotia, 20 homokaryotic offspring isolated from strain L14, and 8 homokaryotic offspring isolated from strain L12.The genome was sequenced using the high-throughput sequencing instrument DNBSEQ-T7 (BGI, Shenzhen, China).

Sequencing Depth Analysis
The filtered clean data were mapped to the assembled genome using Minimap2 v2. 1 [93] for data generated by Nanopore sequencing, while BWA 0.7.17 [94] was utilized for Illumina and DNBSEQ-T7 data.Samtools 0.1.19was used to transfer BAM files, sort, and index them [95].Furthermore, the sequencing depth was calculated using Sambamba 0.6.6 for all 10 kb bins of the genome [96].The sequencing depth data were processed and graphed using WPS Excel version 17147 (WPS Office, Beijing, China) and GraphPad Prism 8.0.2 (Boston, MA, USA).

Conclusions
In conclusion, we have successfully assembled a complete aneuploid chromosomelevel reference genome of W. hoelen for the first time and revealed that the formation of aneuploid genomes may be related to environmental adaptation.Firstly, the first chromosomallevel genome, L7, was assembled here using only the Overlap-Layout-Consensus algorithm.This study revealed 15 chromosomes with double-sided telomere and rDNA sequences for the first time, correcting the previous findings that W. hoelen has 14 chromosomes.Secondly, the aneuploid genomes of W. hoelen were discovered for the first time, and the ploidy of homokaryotic offspring and different W. hoelen strains were studied.The study revealed aneuploidy phenomena in different W. hoelen strains and discovered potential variations in chromosome allocation during the meiosis process of W. hoelen.Thirdly, the gene enrichment analysis was conducted for the genes on reduplicated chromosomes, revealing that the formation of aneuploid genomes of W. hoelen may be associated with environmental adaptation.

Figure 1 .
Figure 1.Circos graph of genome characteristics of Wolfiporia hoelen.From the outer ring to the inner ring are 15 chromosomes, transposable elements, retrotransposons, gene density, and large fragment duplication.The red arrow points to telomere sequences, and the green arrow points to rDNA sequences.

Figure 1 .
Figure 1.Circos graph of genome characteristics of Wolfiporia hoelen.From the outer ring to the inner ring are 15 chromosomes, transposable elements, retrotransposons, gene density, and large fragment duplication.The red arrow points to telomere sequences, and the green arrow points to rDNA sequences.

Figure 4 .
Figure 4. Sequencing depth of different chromosomes for different heterokaryotic strains mapping to genome L7.The red chromosome code presents the reduplicated chromosome, and the red number indicates the reduplicated folds.

Figure 5 .
Figure 5.Sequencing depth of different chromosomes for homokaryotic offspring of L12.The red chromosome code presents the reduplicated chromosome, and the red number indicates the redu plicated folds.

Figure 5 .
Figure 5.Sequencing depth of different chromosomes for homokaryotic offspring of L12.The red chromosome code presents the reduplicated chromosome, and the red number indicates the reduplicated folds.

Figure 6 .
Figure 6.Sequencing depth of different chromosomes for homokaryotic offspring of L14.The red chromosome code presents the reduplicated chromosome, and the red number indicates the reduplicated folds.

Figure 6 .
Figure 6.Sequencing depth of different chromosomes for homokaryotic offspring of L14.The red chromosome code presents the reduplicated chromosome, and the red number indicates the reduplicated folds.

Figure 7 .
Figure 7. GO enrichment analysis of genes on chromosome 15.The rectangle represents enriche first ten GO items, and the ellipse represents other items.Red represents obvious enriched, followe by orange and yellow.

Figure 7 .
Figure 7. GO enrichment analysis of genes on chromosome 15.The rectangle represents enriched first ten GO items, and the ellipse represents other items.Red represents obvious enriched, followed by orange and yellow.

Table 2 .
GO enrichment of genes on chromosome 15 and related response process.