Genome sequencing and analysis of Kloeckera apiculata strain 34-9, a biocontrol agent against postharvest pathogens in citrus

PCR-based gene targeting in Hanseniaspora uvarum

Lack of gene-function analyses tools limits studying the biology of Hanseniaspora uvarum, one of the most abundant yeasts on grapes and in must. We investigated a rapid PCR-based gene targeting approach for one-step gene replacement in this diploid yeast. To this end, we generated and validated two synthetic antibiotic resistance genes, pFA-hygXL and pFA-clnXL, providing resistance against hygromycin and nourseothricin, respectively, for use with H. uvarum. Addition of short flanking-homology regions of 56-80 bp to these selection markers via PCR was sufficient to promote gene targeting. We report here the deletion of the H. uvarum LEU2 and LYS2 genes with these marker genes via two rounds of consecutive transformations, each resulting in the generation of auxotrophic strains (leu2/leu2; lys2/lys2). The hereby constructed leucine auxotrophic leu2/leu2 strain was subsequently complemented in a targeted manner, thereby further validating this approach. PCR-based gene targeting in H. uvarum was less efficient than in Saccharomyces cerevisiae. However, this approach, combined with the availability of two marker genes, provides essential tools for directed gene manipulations in H. uvarum.

Chromosome Genome Sequencing and Comparative Transcriptome-Based Analyses of Kloeckera apiculata 34-9 Unveil the Potential Biocontrol Mechanisms Against Citrus Green Mold

Biological control is an environmentally friendly, safe, and replaceable strategy for disease management. Genome sequences of a certain biocontrol agent could lay a solid foundation for the research of molecular biology, and the more refined the reference genome, the more information it provides. In the present study, a higher resolution genome of Kloeckera apiculata 34-9 was assembled using high-throughput chromosome conformation capture (Hi-C) technology. A total of 8.07 M sequences of K. apiculata 34-9 genome was anchored onto 7 pesudochromosomes, which accounting for about 99.51% of the whole assembled sequences, and 4,014 protein-coding genes were annotated. Meanwhile, the detailed gene expression changes of K. apiculata 34-9 were obtained under low temperature and co-incubation with Penicillium digitatum treatments, respectively. Totally 254 differentially expressed genes (DEGs) were detected with low temperature treatment, of which 184 and 70 genes were upregulated and downregulated, respectively. Some candidate genes were significantly enriched in ribosome biosynthesis in eukaryotes and ABC transporters. The expression of gene Kap003732 and Kap001595 remained upregulated and downregulated through the entire time-points, respectively, indicating that they might be core genes for positive and negative response to low temperature stress. When co-incubation with P. digitatum, a total of 2,364 DEGs were found, and there were 1,247 upregulated and 1,117 downregulated genes, respectively. Biosynthesis of lysine and arginine, and phenylalanine metabolism were the highest enrichment of the cluster and KEGG analyses of the co-DEGs, the results showed that they might be involved in the positive regulation of K. apiculata 34-9 response to P. digitatum. The completeness of K. apiculata 34-9 genome and the transcriptome data presented here are essential for providing a high-quality genomic resource and it might serve as valuable molecular properties for further studies on yeast genome, expression pattern of biocontrol system, and postharvest citrus storage and preservation.

Differentiation of Penicillium roqueforti from Closely Related Species Contaminating Cheeses and Dairy Environment

Currently, Penicillium roqueforti and the closely related P. carneum and P. paneum are identified based on their macromorphology, micromorphology, and molecular properties, the determination of which involves time-consuming procedures. Culture collections focused on dairy isolates of P. roqueforti require quick and efficient tools for routine applications to identify the (a) taxonomy affiliation and (b) morphological properties of strains that influence the sensory properties of blue-veined cheeses. Here, we assessed the morphological variability of P. roqueforti, P. carneum, P. paneum, and P.crustosum on artificial, Edam-like, and Roquefort-like media. Molecular tools were used to test P. roqueforti strains and clones effectively. A novel primer, PrsF, was tested for specificity within strains and isolates of P. roqueforti compared to P. carneum, P. paneum, and P. crustosum. The results reveal that PrsF was specific to the P. roqueforti samples and did not amplify the other tested Penicillium species. Identification based simultaneously on the specificity of the PrsF primer pair and cultivation of P. roqueforti strains on Roquefort-like medium represents an effective method for expanding the collections and practical use of P. roqueforti in the dairy industry.

Development of Genetic Modification Tools for Hanseniasporauvarum

Apiculate yeasts belonging to the genus Hanseniaspora are commonly isolated from viticultural settings and often dominate the initial stages of grape must fermentations. Although considered spoilage yeasts, they are now increasingly becoming the focus of research, with several whole-genome sequencing studies published in recent years. However, tools for their molecular genetic manipulation are still lacking. Here, we report the development of a tool for the genetic modification of Hanseniaspora uvarum. This was employed for the disruption of the HuATF1 gene, which encodes a putative alcohol acetyltransferase involved in acetate ester formation. We generated a synthetic marker gene consisting of the HuTEF1 promoter controlling a hygromycin resistance open reading frame (ORF). This new marker gene was used in disruption cassettes containing long-flanking (1000 bp) homology regions to the target locus. By increasing the antibiotic concentration, transformants were obtained in which both alleles of the putative HuATF1 gene were deleted in a diploid H. uvarum strain. Phenotypic characterisation including fermentation in Müller-Thurgau must showed that the null mutant produced significantly less acetate ester, particularly ethyl acetate. This study marks the first steps in the development of gene modification tools and paves the road for functional gene analyses of this yeast.

Genome-wide identification, characterization and expression analysis of lineage-specific genes within Hanseniaspora yeasts

Lineage-specific genes (LSGs) are defined as genes with sequences that are not significantly similar to those in any other lineage. LSGs have been proposed, and sometimes shown, to have significant effects in the evolution of biological function. In this study, two sets of Hanseniaspora spp. LSGs were identified by comparing the sequences of the Kloeckera apiculata genome and of 80 other yeast genomes. This study identified 344 Hanseniaspora-specific genes (HSGs) and 109 genes ('orphan genes') specific to K. apiculata. Three thousand three hundred thirty-one K. apiculata genes that showed significant similarity to at least one sequence outside the Hanseniaspora were classified into evolutionarily conserved genes. We analyzed their sequence features, functional categories, gene origin, gene structure and gene expression. We also investigated the predicted cellular roles and Gene Ontology categories of the LSGs using functional inference. The patterns of the functions of LSGs do not deviate significantly from genome-wide average. The results showed that a few LSGs were formed by gene duplication, followed by rapid sequence divergence. Many of the HSGs and orphan genes exhibited altered expression in response to abiotic stress. Studying these LSGs might be helpful for understanding the molecular mechanism of yeast adaption.

The Shared and Specific Genes and a Comparative Genomics Analysis within Three Hanseniaspora Strains

Kloeckera apiculata plays an important role in the inhibition of citrus postharvest blue and green mould diseases. This study was based on the previous genome sequencing of K. apiculata strain 34-9. After homologous comparison, scaffold 27 was defined as the mitochondrial (mt) sequence of K. apiculata 34-9. The comparison showed a high level of sequence identity between scaffold 27 and the known mtDNA of Hanseniaspora uvarum. The genome sequence of H. vineae T02/19AF showed several short and discontinuous fragments homologous to the mtDNA of H. uvarum. The shared and specific genes of K. apiculata, H. uvarum, and H. vineae were analysed by family using the TreeFam methodology. GO analysis was used to classify the shared and specific genes. Most of the gene families were classified into the functional categories of cellular component and metabolic processes. The whole-genome phylogram and genome synteny analysis showed that K. apiculata was more closely related to H. uvarum than to H. vineae. The genomic comparisons clearly displayed the locations of the homologous regions in each genome. This analysis could contribute to discovering the genomic similarities and differences within the genus Hanseniaspora. In addition, some regions were not collinearity-matched in the genome of K. apiculata compared with that of H. uvarum or H. vineae, and these sequences might have resulted from evolutionary variations.