Liang Junmin, Li Yuanjie, Dodds Peter N, Figueroa Melania, Sperschneider Jana, Han Shiling, Tsui Clement K M, Zhang Keyu, Li Leifu, Ma Zhanhong, Cai Lei
State Key Laboratory of Mycology, Institute of Microbiology, Chinese Academy of Sciences, Beijing, China.
College of Life Sciences, University of Chinese Academy of Sciences, Beijing, China.
Mol Ecol Resour. 2023 Apr;23(3):601-620. doi: 10.1111/1755-0998.13739. Epub 2022 Dec 8.
Rust fungi are characterized by large genomes with high repeat content and have two haploid nuclei in most life stages, which makes achieving high-quality genome assemblies challenging. Here, we described a pipeline using HiFi reads and Hi-C data to assemble a gigabase-sized fungal pathogen, Puccinia polysora f.sp. zeae, to haplotype-phased and chromosome-scale. The final assembled genome is 1.71 Gbp, with ~850 Mbp and 18 chromosomes in each haplotype, being currently one of the two giga-scale fungi assembled to chromosome level. Transcript-based annotation identified 47,512 genes for the dikaryotic genome with a similar number for each haplotype. A high level of interhaplotype variation was found with 10% haplotype-specific BUSCO genes, 5.8 SNPs/kbp, and structural variation accounting for 3% of the genome size. The P. polysora genome displayed over 85% repeat contents, with genome-size expansion and copy number increasing of species-specific orthogroups. Interestingly, these features did not affect overall synteny with other Puccinia species having smaller genomes. Fine-time-point transcriptomics revealed seven clusters of coexpressed secreted proteins that are conserved between two haplotypes. The fact that candidate effectors interspersed with all genes indicated the absence of a "two-speed genome" evolution in P. polysora. Genome resequencing of 79 additional isolates revealed a clonal population structure of P. polysora in China with low geographic differentiation. Nevertheless, a minor population differentiated from the major population by having mutations on secreted proteins including AvrRppC, indicating the ongoing virulence to evade recognition by RppC, a major resistance gene in Chinese corn cultivars. The high-quality assembly provides valuable genomic resources for future studies on disease management and the evolution of P. polysora.
锈菌的特点是基因组大、重复序列含量高,并且在大多数生命阶段有两个单倍体细胞核,这使得获得高质量的基因组组装具有挑战性。在这里,我们描述了一种使用高保真 reads 和 Hi-C 数据的流程,用于将千兆碱基大小的真菌病原体玉米多堆柄锈菌组装到单倍型定相和染色体水平。最终组装的基因组大小为 1.71 Gbp,每个单倍型约为 850 Mbp 且有 18 条染色体,是目前已组装到染色体水平的两个千兆规模真菌之一。基于转录本的注释确定双核基因组有 47,512 个基因,每个单倍型的基因数量相似。发现单倍型间变异水平较高,有 10% 的单倍型特异性 BUSCO 基因、每千碱基 5.8 个 SNP,结构变异占基因组大小的 3%。玉米多堆柄锈菌基因组的重复序列含量超过 85%,基因组大小扩展且物种特异性直系同源组的拷贝数增加。有趣的是,这些特征并不影响与其他基因组较小的柄锈菌物种的整体共线性。精细时间点转录组学揭示了两个单倍型之间保守的七组共表达分泌蛋白。候选效应子散布在所有基因中的事实表明玉米多堆柄锈菌不存在“双速基因组”进化。对另外 79 个分离株的基因组重测序揭示了中国玉米多堆柄锈菌的克隆群体结构,地理分化较低。然而,一个小群体与主要群体有所分化,其分泌蛋白包括 AvrRppC 发生了突变,这表明它正在进化出逃避中国玉米品种中主要抗性基因 RppC 识别的毒力。高质量的组装为未来玉米多堆柄锈菌的病害管理和进化研究提供了有价值的基因组资源。