High-quality genome assembly and genetic transformation system of Lasiodiplodia theobromae strain LTTK16-3, a fungal pathogen of Chinese hickory

ABSTRACT Lasiodiplodia theobromae, as one of the causative agents associated with Chinese hickory trunk cankers, has caused huge economic losses to the Chinese hickory industry. Although the biological characteristics of this pathogen and the occurrence pattern of this disease have been well studied, few studies have addressed the related mechanisms due to the poor molecular and genetic study basis of this fungus. In this study, we sequenced and assembled L. theobromae strain LTTK16-3, isolated from a Chinese hickory tree (cultivar of Linan) in Linan, Zhejiang province, China. Phylogenetic analysis and comparative genomics analysis presented crucial cues in the prediction of LTTK16-3, which shared similar regulatory mechanisms of transcription, DNA replication, and DNA damage response with the other four Chinese hickory trunk canker-associated Botryosphaeria strains including, Botryosphaeria dothidea, Botryosphaeria fabicerciana, Botryosphaeria qingyuanensis, and Botryosphaeria corticis. Moreover, it contained 18 strain-specific protein clusters (not conserved in the other L. theobromae strains, AM2As and CITRA15), with potential roles in specific host-pathogen interactions during the Chinese hickory infection. Additionally, an efficient system for L. theobromae protoplast preparation and polyethylene glycol (PEG) -mediated genetic transformation was firstly established as the foundation for its future mechanisms study. Collectively, the high-quality genome data and the efficient transformation system of L. theobromae here set up the possibility of targeted molecular improvements for Chinese hickory canker control. IMPORTANCE Fungi with disparate genomic features are physiologically diverse, possessing species-specific survival strategies and environmental adaptation mechanisms. The high-quality genome data and related molecular genetic studies are the basis for revealing the mechanisms behind the physiological traits that are responsible for their environmental fitness. In this study, we sequenced and assembled the LTTK16-3 strain, the genome of Lasiodiplodia theobromae first obtained from a diseased Chinese hickory tree (cultivar of Linan) in Linan, Zhejiang province, China. Further phylogenetic analysis and comparative genomics analysis provide crucial cues in the prediction of the proteins with potential roles in specific host-pathogen interactions during the Chinese hickory infection. An efficient PEG-mediated genetic transformation system of L. theobromae was established as the foundation for the future mechanisms exploration. The above genetic information and tools set up valuable clues to study L. theobromae pathogenesis and assist in Chinese hickory canker control.


Genome sequencing, assembly, and annotation
To generate the basis for studying the origins and mechanisms behind the L. theobro mae survival strategies and pathogenicity mechanisms, L. theobromae strain LTTK16-3, isolated from a Chinese hickory tree (cultivar of Linan) in Linan, Zhejiang province of China (6), was used for genome sequencing.L. theobromae has an estimated genome size of 42.52 Mb based on 21 K-mer analysis, and the K-mer distributions followed a Poisson distribution (Fig. 1A).As shown in Table 1, approximately 6 Gb ONT reads were obtained (reads coverage 139×).The de novo genome assemblies found that the assembly size of LTTK16-3 strain was 42.82 Mb, the total contig numbers was 10, contig N 50 was 5.67 Mb, and the maximum contig length was 7.93 Mb (Table 1).The completeness of genome assemblies assessed by BUSCO v5.12 (https:// busco.ezlab.org/)found that the LTTK16-3 strain contains 98.71% complete orthologs at the Ascomycota level (n = 1,706) (Fig. 1B; Table 1).The telomere repeats determined at the start or end region of contigs (5′-TTAGGG-3′/5′-CCCTAA-3′) showed that the assembly of strain LTTK16-3 contained six contigs with (TTAGGG)n start, five contigs with (CCCTAA)n end, and three contigs reached telomere-to-telomere (T2T) chromoso mal level (Table 1).Repeats masked before gene prediction found that the repeat content of the LTTK16-3 strain was 3.37% (Table 1; Fig. 1B).The gene prediction showed that a total of 12,516 protein-coding genes were identified in the repeatmasked genome assembly of LTTK16-3 strain (Table 1).Additionally, as shown in Fig. 2, gene density ranged from one to eight genes per 100 kilobases (kb) across the chromosomes, and the GC contents of the total genomes were 54.57% (Table 1; Fig. 2).Intra-genomic syntenic analysis only detected nine syntenic blocks containing 75 pairs of homoeologous genes in the genome of LTTK16-3 strain, which is consistent with its relatively low genomic heterozygosity (Fig. 1A and 2).The gene functional annotation listed in Table 1 shows that the LTTK16-3 strain contained around 2,457 pathogenhost interaction genes, 237 carbohydrate-active enzymes, and 190 cytochrome P450 enzymes.Additionally, a total of 715 putative secreted proteins were identified using our previously defined pipeline (Table 1) (23), and 51 secondary metabolite biosyn thetic genes were identified using online antiSMASH (https://fungismash.secondaryme tabolites.org/#!/start) (Table 1) (24,25).

Phylogenetic analysis and comparative genomics analysis
The phylogenetic tree of the LTTK16-3 strain was made based on the result of orthogroups using the Species Tree inference from All Genes (STAG) and Species Tree Root Inference from gene Duplication Events (STRIDE) in OrthoFinder software (26).As shown in Fig. 3A, the phylogenetic tree revealed that LTTK16-3 clustered with reported L. theobromae strain CITRA15 isolated from citrus and L. theobromae strain AM2As

GWH accession a GWHBFSF00000000
ONT reads (Gb) cortices formed a single monophyletic clade individually, indicating that among these Chinese hickory trunk canker causative agents, LTTK16-3 strain was genetically far away from the other four Botryosphaeria strains.Consistently, orthologous protein cluster analysis conducted at the online web service OrthoVenn2 (https://orthovenn2.bioinfo toolkits.net/home)found that LTTK16-3 contained the most species-specific protein clusters, i.e., 131, while the other four Botryosphaeria strains only contained 0-21 species-specific protein clusters (Fig. 3B).Additionally, these five Chinese hickory trunk canker causative fungi shared 8,818 core orthologous protein clusters (includ ing 9,095 orthologous proteins) (Fig. 3B), and a Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis mapped above 8,818 core orthologous protein clusters to 137 KEGG pathways (Table S2) with seven significantly enriched KEGG pathways (P < 0.05), including the ribosome, basal transcription factor, RNA polymerase, nucleotide excision repair, RNA degradation, DNA replication, and mismatch repair (Fig. S1), which indicates the similar regulatory mechanisms of transcription, DNA replication, and DNA damage response among the five Chinese hickory trunk canker causative pathogens.Intriguingly, comparative analysis between LTTK16 and other previously reported L. These density metrics were calculated with 100 kb sliding windows.The syntenic genomic blocks are illustrated with innermost lines.
theobromae strains, including AM2As and CITRA15 (E-value cutoff is 1e-05 and inflation value cutoff is 2.0) showed that a total of 8,975 orthogroups were shared by LTTK16, AM2As, and CITRA15, while 18, 32, and 6 orthogroups were found to be specific to LTTK16, AM2As, and CITRA15, respectively (Fig. 3C).Additionally, bidirectional BLAST analysis found that most (72%) of the above 18 LTTK16-3 strain-specific protein clusters (including 39 LTTK16-3 strain-specific proteins) have orthologous proteins in the other four Chinese hickory trunk canker-related Botryosphaeria strains (Table S3).Further analysis of above the LTTK16-3 strain-specific proteins by GO (Table S4) and Pfam domain annotation (Table S5) found two specific kinases, LT_TK16_09351 and LT_TK16_09352, containing the Pkinase domain (PF00069) enriched in the GO term "intracellular signal transduction" (GO:0035556), which might contribute to the specific survival strategies of Chinese hickory pathogens.

Sensitivity of L. theobromae to geneticin (G-418) and protoplast Preparation
To establish a transformation system for L. theobromae, the sensitivity to G-418 was first conduced for subsequent selection of transformants carrying the functional G-418 resistant gene.The sensitivity of the wild-type L. theobromae strain LTTK16-40 was examined on potato dextrose agar (PDA) medium containing G-418 at concentra tions ranging from 0 to 400 µg/mL.As shown in Fig. 4, ≥100 µg/mL of G-418 was able to significantly inhibit the LTTK16-40 growth, and ≥200 µg/mL of G-418 com pletely inhibited the LTTK16-40 growth.Therefore, successful transformants carrying the functional G-418 resistant gene were screened using PDA with 200 µg/mL of G-418 (Fig. 4A and B).
To screen the best enzymatic digestion condition for protoplast preparation during the PEG-mediated protoplast transformation, the prepared fresh mycelia were placed in 10 mL of enzymatic digestion solution (0.3 g cellulase, 0.3 g lysozyme, 0.25 g lysing enzyme, and 0.08 g driselase) and shaken at 30°C, 100 rpm.The protoplasts began to form after 0.5 h of enzymatic treatment (Fig. S2).The amount of protoplasm production increased continuously during the first 3 h and then leveled off.After 4-h of enzymatic digestion, the quality of protoplasm production started to decrease, and the cell membranes of some protoplasm ruptured with irregular shape (Fig. 4C and D).Therefore, it was proved that the suitable enzymatic digestion time for L. theobromae protoplast preparation during transformation was 3-4 h (Fig. 4C and D).

Screening and detection of transformants
The two recombinant pYF11-neo plasmids containing LtActin-GFP fusion cassette and LtH1-GFP fusion cassette were transferred into the wild-type L. theobromae strain LTTK16-40, named LTTK16-40:: LtActin-GFP and LTTK16-40::LtH1-GFP, respectively.The corresponding transformants were selected using PDA selection medium containing 200 µg/mL G-418 and verified by PCR identification together with sequencing to ensure the accuracy of the in-frame fusion region (Fig. 5A through C).As shown in Fig. 5B and C, the confirmed transformants were able to grow on PDA selection medium, while the wild-type L. theobromae strain LTTK16-40 was unable to grow.Addition ally, further confocal microscopic examination showed that the above transformants exhibited intense fluorescence, suggesting that the designed PEG-mediated protoplast transformation system in this study works successfully for L. theobromae transforma tion (Fig. 5D).The recombinant pYF11-neo plasmids containing GFP-fused genes were successfully transferred into LTTK16-40 via this system, and the corresponding GFP-fused protein could stably express in L. theobromae (Fig. 5D).Additionally, growth and pathogenicity examination found that the growth rate and virulence of transform ants LTTK16-40:: LtActin-GFP and LTTK16-40::LtH1-GFP demonstrated no significant differences to that of the wild-type L. theobromae strain LTTK16-40 (Fig. 6), indicating this transformation system had no side effect on L. theobromae growth and virulence.

DISCUSSION
As one of the causative agents associated with Chinese hickory trunk cankers, L. theobromae has caused huge economic losses to the Chinese hickory industry due to its extremely strong virulence (6).Besides damaging Chinese hickory trees, this pathogen is becoming a significant threat to various woody plants with wide geographic distribution (8,9,12).Here, we sequenced and assembled the LTTK16-3 strain, the genome of L. theobromae first obtained from a diseased Chinese hickory tree (cultivar of Linan) in Linan, Zhejiang province of China.The phylogenetic tree revealed that among the five Chinese hickory trunk canker causative agents, the LTTK16-3 strain was genetically far away from the other four Botryosphaeria strains.Comparing the genomic information of L. theobromae strain LTTK16-3 with the published genomes of the other four Chinese hickory trunk canker-associated Botryosphaeria species (19,20) found that the genome assembly of strain LTTK16-3 had the minimum assembly size and contig number, as well as the largest contig N 50 and the maximum contig length (Table 1).Addition ally, the repeat contents of strain LTTK16-3 (3.37%) were the smallest one (3.37%)compared to that of the other four reported Chinese hickory canker causal agents, B. dothidea (BDLA16-7 strain, 7.96%), B. cortices (BCTK16-35 strain, 8.50%), B. fabicerciana (BFLG18-2 strain, 8.20%), and B. qingyuanensis (BQTK16-30 strain, 6.03%) (19,20).The orthologous protein cluster analysis found that LTTK16-3 contained 131 species-specific protein clusters that are not conserved among other Chinese hickory trunk canker causative Botryosphaeria species (Fig. 3B).The difference initially appeared in genomics is ultimately displayed at the cellular level as disparate strategies for survival and infection.Our previous study has reported that L. theobromae was the most virulent pathogen of Chinese hickory among the Botryosphaeriaceae species (6).Thus, we speculated that the above species-specific protein clusters identified here likely support that discrepancy in virulence.
In fungi, protein kinases (PKs), as well as the reversible protein phosphorylation, play critical roles in signal transduction, growth, sexual reproduction, and stress responses (27,28).For instance, Fusarium graminearum, a devastating fungal pathogen of cereal crops, has 116 PK genes, of which 20 of them appear to be essential, 42 of them are important for virulence, 32 of them are associated with vegetative growth, and 20 of them play important roles in sexual production, and the corresponding mutants fail to produce perithecia (28).Additionally, PKs, such as FgYak1, FgSsk2, FgPbs2, and FgHog1, are the key regulators in host-derived oxidative stress response during the infection (28).The kinase FgCak1 regulating the phosphorylation of FgNem1 at Thr72, Ser187, and Ser287 together with kinase FgTor and phosphatase FgPpg1 under the target of rapamycin pathway plays crucial roles in F. graminearum lipid droplet biogen esis, vegetative growth, sexual development, and virulence (29).In Botrytis cinerea, a fungal pathogen causing pre-and post-harvest diseases in more than 200 economically important crops, the mitogen-activated protein kinase BcMkk1 responsible for orches trating changes of fungal cell wall positively regulates pathogenicity via the cell wall integrity pathway, which controls vegetative growth, conidiation, and responses to cell wall and oxidative stresses (30).Thus, the PKs with specific functions contribute to the related infection and survival strategies of fungal pathogens (28)(29)(30).In the L. theobromae strain LTTK16-3, there are two predicted strain-specific PKs, LT_TK16_09351 and LT_TK16_09352, that are not conserved in the other L. theobromae strains (AM2As and CITRA15) and potentially involved in intracellular signal transduction.Therefore, investigating the potential functions of the above two predicted PKs could be a valuable way to further explore the mechanisms involved in the specific host-pathogen interac tions during the Chinese hickory infection.
The development of genetic transformation techniques is a precondition for gene function and related mechanism analysis (31,32).In filamentous fungi, the PEG-medi ated genetic transformation system has been widely demonstrated to be a simple and effective system for transformation (33,34).PEG is a crucial agent enhancing transforma tion efficiency (35).It can promote adhesion by forming the molecular bridge between cytomembrane and exogenous nucleic acids (18).Additionally, it can also alter the membrane permeability and facilitate the entry of DNA into cells by inducing disordered charges on the cytomembrane surface (18).The PEG-mediated genetic transformation mainly contains two steps, including protoplast preparation and PEG-mediated fusion of exogenous nucleic acids and protoplasts.Due to the high variety of the fungal cell wall components among different species (18), there is no universal transformation method that can be applied to different fungal strains.Here, we first established a PEG-medi ated genetic transformation system with an efficient protoplast preparation protocol for L. theobromae.Compared to the ATMT system for L. theobromae transformation, which employs pycnidiospores as transformation targets (21), the PEG-mediated genetic transformation system here directly targets fungal protoplasts.Yet, sporulation is a difficult task, and sporulation efficiency can be negatively compromised by temperature, light, and substrate regimes (21).Thus, the employment of protoplasts, rather than pycnidiospores, as the biological target in transformation is more appropriate, which is an improvement in advancing the transformation efficiency.
Overall, here we sequenced and assembled the LTTK16-3 strain, the genome of L. theobromae first obtained from a diseased Chinese hickory tree (cultivar of Linan) in Linan, Zhejiang province of China (Fig. 1 and 2).Phylogenetic analysis and comparative genomics analysis found that though L. theobromae LTTK16-3 strain was genetically far away from the other four Botryosphaeria strains, all five Chinese hickory trunk canker causative pathogens shared 8,818 core orthologous protein clusters (Fig. 3).KEGG pathway analysis of the above core orthologous protein clusters presented crucial cues in the prediction of the similar regulatory mechanisms of transcription, DNA replication, and DNA damage response among five pathogens (Fig. S1; Table S2).Orthologous protein cluster analysis of LTTK16 and other previously reported L. theobromae strains AM2As and CITRA15 showed that LTTK16-3 strain contained 18 LTTK16-3 strain-specific protein clusters (including 39 LTTK16-3 strain-specific proteins), and bidirectional BLAST analysis found that 72% of the above strain-specific proteins have orthologous proteins in the other four Chinese hickory trunk canker related-Botryosphaeria strains (Table S3), indicating their potential roles in specific host-pathogen interactions during the Chinese hickory infection.Thus, we developed an efficient PEG-mediated genetic transformation system establishing the foundation for future studies on the mechanisms of L. theobro mae and setting up the possibility of targeted molecular improvements for Chinese hickory canker control.

Fungal strain, culture conditions, and geneticin (G-418) sensitivity determina tion
The wild-type Lasiodiplodia theobromae strain LTTK16-3 used in genomic sequence and genetic transformation was obtained from a diseased Chinese hickory tree (cultivar of Linan) in Linan, Zhejiang province of China.The LTTK16-3 strain and the resulting transformants were stored on the potato dextrose agar (200 g potato, 20 g glucose, and 20 g agar per liter of pure water) slants at 4°C.To determine the sensitivity to G-418, 5-mm mycelial plugs of each strain taken from a 36-h-old colony edge were inoculated on PDA supplemented without/with each stress agent and then incubated at 25°C in the dark.The concentrations for G-418 were indicated in figure legends.Each experiment was repeated three times independently.

DNA extraction, genome sequencing, and assembly
For the preparation of genomic DNA used in genome sequencing and assembly, the L. theobromae strain LTTK16-3 was incubated on PDA plates at 25°C in the dark.After mycelia grew to cover nearly two-thirds of the PDA plate surfaces, the cultures were transferred to a mortar and ground with liquid nitrogen.The resultant powder was placed in a 2-mL centrifuge tube, and the genomic DNA was extracted using a Genomic DNA Kit (Sangon Biotech Co., Ltd., Shanghai) according to the manufac turer's instructions.After assessing the purity and concentration of the extracted genomic DNA with a NanoDrop One spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA), the genomic sequence was analyzed by the Oxford Nanopore Technologies (ONT) PromethION sequencing platform at Biomarker Technologies Co., Ltd (Beijing, China).Additionally, the RNA sequencing was also performed on the Illumina HiSeq4000 sequencing platform using the above mycelium materials collected from PDA media to provide transcription evidence.The related RNA-seq data were deposited in Genome Sequence Archive (GSA) under the accession number CRR332439.The lowquality readings were filtered and high-quality filtered sub-reads were assembled using NextDenovo v2.4.0 and NextPolish v1.3.1 (both available online at https://github.com/Nextomics).A 21-mer was selected for k-mer analysis, and the 21-mer depth fre quency distribution was calculated using jellyfish (version 1.1.12).The genome size and heterozygosity were visualized using GenomeScope (version 1.0).The completeness of genome assemblies was assessed using BUSCO v5.12 (https://busco.ezlab.org/).The telomere repeats were determined at the start or end region of contigs (5′-TTAGGG-3′/ 5′-CCCTAA-3′) to test if contig has reached T2T chromosomal level (36).Repeats were masked before gene prediction by RepeatMasker v4.1.2(http://www.repeatmasker.org/) using a de novo repeat library generated by RepeatModeler v2.01 (http://www.repeatmasker.org/RepeatModeler/).The raw sequence data are publicly available at GSA under the accessions number CRR332440, and the genome assembly is publicly available at the Genome Warehouse (GWH, https://ngdc.cncb.ac.cn/gwh) under the accession number GWHBFSF00000000 and under BioProject PRJCA005744 in National Genomics Data Center, China National Center for Bioinformation (CNCB-NGDC Members and Partners, 2021, https://ngdc.cncb.ac.cn).Also, the corresponding genome sequencing data have been uploaded to the National Center for Biotechnology Information (NCBI) website under the accession number PRJNA1056688.

Phylogenetic analysis
A homologous single-copy gene-based approach was applied to generate a phyloge netic tree of LTTK16-3 genome sequences.The STAG algorithm of OrthoFinder was used to reconstruct a phylogenetic tree based on homologous genes, and the STRIDE algorithm was used to root the species tree in OrthoFinder (25,26).The genomic data of Diplodia seriata (38), L. theobromae strain AM2As (12), and L. theobromae strain CITRA15 (8) used in phylogenetic analysis were obtained from the NCBI web portal.The genomic data of the other four Chinese hickory trunk cankers associated Botryosphaeria species (including B. dothidea, B. fabicerciana, B. qingyuanensis, and B. cortices) published in our previous study were obtained from the GWH (https://ngdc.cncb.ac.cn/gwh) (19,20).

Construction of recombinant plasmids
The corresponding native promoter and open reading frame of LtActin and LtH1 were amplified using the primer pairs listed in Table S1.For generating the native promoter-LtActin-GFP fusion construct and native promoter-LtH1-GFP fusion construct, the above native promoter-LtActin and promoter-LtH1 fragment were co-transformed with XhoI-digested pYF11 into yeast strain XK125, respectively.The pYF11 contains the neomycin resistance gene, encoding aminoglycoside phosphotransferase, which allows the selection of the resulting fungal transformants using geneticin (G-418), and its map is provided in Fig. S3.The yeast plasmids pYF11-LtActin-GFP and pYF11-LtH1-GFP were rescued from the resulting Trp + yeast transformants.Then, the recombinant plasmids were extracted and transformed into Escherichia coli DH5α for large-scale amplification.The resulting recombinant plasmids were verified by sequencing to ensure the accuracy of the in-frame fusion region.The verified recombinant plasmids were selected to construct the fluorescent protein-labeled strains LTTK16-40:: LtActin-GFP and LTTK16-40::LtH1-GFP by transforming these two recombinant plasmids into the protoplasts of the LTTK16-3 strain.

Fluorescent protein-labeled strain construction
To generate protoplasts, mycelial plugs cut from 36-h-old colony edge were placed into 250-mL flasks containing 100 mL yeast extract peptone dextrose (1% yeast extract, 2% peptone, and 2% glucose) culture.Then, the flasks were kept on a rotary shaker (175 rpm, 25°C) for 16 h, and mycelia were collected using a sterile nylon filter and washed with distilled water thrice to remove medium residue.The harvested mycelia were re-suspen ded in 10 mL enzymatic digestion solution containing 0.7 M NaCl, 0.3 g cellulase (Ryon Bio-Tech, Shanghai, China, RM1030), 0.3 g lysozyme (Ryon Bio-Tech, Shanghai, China, RM1027), 0.25 g lysing enzyme (Sigma, St. Louis, MO, USA, L1412), and 0.08 g driselase (Sigma-Aldrich, St. Louis, MO, USA, D9515) at 30°C and 85 rpm.The state of protoplasts waThe verified recombinant s monitored during the enzymatic digestion via microscopy.Protoplasts were separated from cell debris via filtration through three layers of lens-cleaning tissue (2105-802, GE Healthcare, Chicago, IL, USA) and gently rinsed twice with 0.7 M NaCl buffer.Protoplasts were collected via centrifuging at 5,000 rpm at 4°C for 8 min and re-suspended in 5 mL STC buffer [0.8 M sorbitol and 50 mM Tris-HCl (pH 8.0)] twice.For transformation, 750 µL of protoplasts, 200 µL of SPTC (STC with 40%, wt/vol PEG 6000) buffer, 5 µL of heparin (5 mg/mL), and 100 µL (>200 µg/mL) of recombinant plasmids were mixed and incubated on ice for 30 min; then, 400 µL of SPTC was added into the above suspension and incubated at 25°C for 20 min.Transformed protoplasts were mixed into 100 mL RM medium (1 g yeast extract, 1 g casein hydrolysate, 274 g sucrose, and 16 g agar powder per 1 L pure water) at 35°C, poured into petri plates, and incubated at 25°C in dark.After 12 h, RM plates were overlaid with 10 mL of SRM medium (1 g yeast extract, 1 g casein hydrolysate, 342 g sucrose, and 12 g agar per 1 L pure water) containing G-418 for transformant selection and incubated at 25°C.After 2-6 days, geneticin-resistant colonies appeared, and individual transformants were transferred to PDA plates containing G-418 for fluorescence examination.The transformants expressing the GFP-fused protein in right fluorescence localizations were transferred to PDA without G-418 to grow for PCR identification and sequencing.After the PCR identification and sequencing, the right transformants were saved in a PDA slant without G-418 to be used as an inoculation source in the subsequent experiments.

Microscopic examinations
Fluorescence signals were examined with a Zeiss LSM780 confocal microscope (Gottingen, Niedersachsen, Germany).The subcellular locations of LtActin-GFP and LtH1-GFP were examined using the following confocal microscopy settings: laser 488 nm at 50% power, pinhole 90 µm, and master gain 580.

FIG 3
FIG 3 Phylogenetic analysis and comparative genomics analysis of L. theobromae strain LTTK16-3.(A) The phylogenetic tree of the LTTK16-3 strain.The species tree was inferred as part of the OrthoFinder pipeline using the STAG algorithm and rooted by STRIDE.STRIDE probability values are shown at internal nodes.Scale bar indicates the number of substitutions per site.(B) Venn diagram depicting the number of shared and specific protein clusters among five Botryosphaeriaceae species.(C) Venn diagram depicting the number of shared and specific protein clusters among three L. theobromae strains.