Phytophthora theobromicola is an emerging cacao pathogen recently identified in Brazil as an aggressive agent of black pod rot. We generated genome assemblies for two P. theobromicola isolates using long-read sequencing and five additional isolates using short reads. Comparative analysis revealed a genome size and predicted gene content comparable to P. citrophthora, a closely related species with a broad host range that includes both citrus and cacao. An intraspecies sequence-graph analysis revealed a highly dynamic genome structure with high proportion of variable effectors. Syntenic orthology analysis across 13 Phytophthora species identified orthologous gene groups conserved only in cacao pathogens and others specific to P. theobromicola. RxLR effectors and CAZymes were particularly enriched among lineage-specific syntenic groups, with RxLRs preferentially located near transposable elements and within gene-sparse, repeat-rich regions. Transcriptome analysis of infected cacao tissues showed that 88% of predicted effectors were expressed, with pods exhibiting the highest number of upregulated genes. Notably, several RxLRs classified as P. theobromicola-specific syntenic orthologs were highly expressed in infected tissues, suggesting that these lineage-specific effectors may play key roles in host-pathogen interactions unique to cacao. Together, our findings highlight the dynamic architecture and functional plasticity of the P. theobromicola genome, providing foundational insights into its virulence strategies and supporting future studies on host adaptation and effector evolution in emerging cacao pathogens.

Marasmius crinis-equi, known as the horse hair fungus, is a litter-decomposing fungus that is opportunistically pathogenic on various plants. We sequenced the genome of M. crinis-equi strain GH-76 isolated from thread blight disease-infected cacao leaves in Ghana. The genome is 72.8 Mbp in size and 91.2% complete via benchmarking universal single-copy orthologs (BUSCO).

The oomycete Phytophthora palmivora infects the fruit of cacao trees (Theobroma cacao) causing black pod rot and reducing yields. Cacao genotypes vary in their resistance levels to P. palmivora, yet our understanding of how cacao fruit respond to the pathogen at the molecular level during disease establishment is limited. To address this issue, disease development and RNA-Seq studies were conducted on pods of seven cacao genotypes (ICS1, WFT, Gu133, Spa9, CCN51, Sca6 and Pound7) to better understand their reactions to the post-penetration stage of P. palmivora infection. The pod tissue-P. palmivora pathogen assay resulted in the genotypes being classified as susceptible (ICS1, WFT, Gu133 and Spa9) or resistant (CCN51, Sca6 and Pound7). The number of differentially expressed genes (DEGs) ranged from 1625 to 6957 depending on genotype. A custom gene correlation approach identified 34 correlation groups. De novo motif analysis was conducted on upstream promoter sequences of differentially expressed genes, identifying 76 novel motifs, 31 of which were over-represented in the upstream sequences of correlation groups and associated with gene ontology terms related to oxidative stress response, defense against fungal pathogens, general metabolism and cell function. Genes in one correlation group (Group 6) were strongly induced in all genotypes and enriched in genes annotated with defense-responsive terms. Expression pattern profiling revealed that genes in Group 6 were induced to higher levels in the resistant genotypes. An additional analysis allowed the identification of 17 candidate cis-regulatory modules likely to be involved in cacao defense against P. palmivora. This study is a comprehensive exploration of the cacao pod transcriptional response to P. palmivora spread after infection. We identified cacao genes, promoter motifs, and promoter motif combinations associated with post-penetration resistance to P. palmivora in cacao pods and provide this information as a resource to support future and ongoing efforts to breed P. palmivora-resistant cacao.

Thread blight disease (TBD) is an emerging foliar disease on cacao (Theobroma cacao), the source of chocolate. Multiple species of pathogenic fungi in the Marasmiinae suborder of Agaricales are responsible for causing TBD in cacao. A recent phylogenetic study showed that Marasmius palmivorus, the second most frequently isolated TBD-causing fungus in Ghana, was phylogenetically distinct from the Marasmiaceae genera Marasmius and Moniliophthora, leading to the naming of a new genus, Paramarasmius. The genus currently contains two species, the type species P. palmivorus and P. mesosporus. This report describes our release of a draft nuclear genome assembly and annotation of P. palmivorus strain GH-12 derived from long- and short-read sequencing. This high-quality draft is the first released genome for P. palmivorus and provides insight into the pathogenicity and evolution of an emerging pathogen.

Members of the Marasmiaceae family of basidiomycetous fungi cause thread blight disease (TBD) in cacao, Theobroma cacao, and are a threat to chocolate production worldwide. The most cited causal agent of TBD is Marasmiellus scandens; however, among recent TBD-causing Marasmiaceae isolations in Ghana, only 8% were M. scandens. Nonetheless, owing to the emerging threat of TBD to cacao, generating a high-quality genome assembly of the pathogen M. scandens is of considerable importance. This announcement describes our assembly, annotation, and release of a draft nuclear genome for M. scandens strain GH-19 using a hybrid assembly approach combining long and short read sequencing data.

Dieback disease is one of the constraints to cocoa production in Indonesia. A newly identified dieback symptom has been observed in cocoa farms in East Luwu since 2015 and later observed in Luwu, Pinrang, Enrekang and Soppeng Regency, South Sulawesi. The potential causal agent identity was determined using morphological characteristics, DNA sequences, and phylogenetic analyses. A total of four fungal isolates out of seven isolates that were recovered from diseased stems morphologically resembled Lasiodiplodia species. ITS and EF1 alpha sequencing confirmed that three of the Lasiodiplodia isolates were Lasiodiplodia theobromae, and one isolate was Lasiodiplodia pseudotheobromae. The aggresiveness of all fungal isolates and control was further determined by a leaf disk assay. Experiments showed that L. theobromae isolate CAS0321 was the most aggressive isolate in a leaf disk assay, while other Lasiodiplodia isolates, Fusarium isolates and Diaporthe isolate, were less aggressive. L. theobromae isolate CAS0321 was further used for a pathogenicity test on cocoa seedlings. Two months after inoculation, L. theobromae isolate CAS0321 caused elongated streaks on seedling stems with dark brown to black streaking of vascular tissues similar to that observed in natural infections. This is the first report of a virulent isolate of L. theobromae associated with cocoa dieback in Sulawesi.

Theobroma cacao, the chocolate tree, is indigenous to the Amazon basin, the greatest biodiversity hotspot on earth. Recent advancement in plant genomics highlights the importance of de novo sequencing of multiple reference genomes to capture the genome diversity present in different cacao populations. In this study, three high-quality chromosome-level genomes of wild cacao were constructed, de novo assembled with HiFi long reads sequencing, and scaffolded using a reference-free strategy. These genomes represent the three most important genetic clusters of cacao trees from the Upper Amazon region. The three wild cacao genomes were compared with two reference genomes of domesticated cacao. The five cacao genetic clusters were inferred to have diverged in the early and middle Pleistocene period, approximately 1.83-0.69 million years ago. The results shown here serve as an example of understanding how the Amazonian biodiversity was developed. The three wild cacao genomes provide valuable resources for studying genetic diversity and advancing genetic improvement of this species.

Here, we report a hybrid genome of Marasmius tenuissimus strain MS-2, a cacao thread blight disease causing isolate that was collected from cacao leaves in Tafo, Eastern region, Ghana. The final assembly consists of 2,083 contigs spanning 69,843,039 bp, with 49.21% GC content, 92.6% BUSCO completeness, and scaffold N50 186,871.

The basidiomycete Moniliophthora roreri causes frosty pod rot of cacao (Theobroma cacao) in the western hemisphere. Moniliophthora roreri is considered asexual and haploid throughout its hemibiotrophic life cycle. To understand the processes driving genome modification, using long-read sequencing technology, we sequenced and assembled 5 high-quality M. roreri genomes out of a collection of 99 isolates collected throughout the pathogen's range. We obtained chromosome-scale assemblies composed of 11 scaffolds. We used short-read technology to sequence the genomes of 22 similarly chosen isolates. Alignments among the 5 reference assemblies revealed inversions, translocations, and duplications between and within scaffolds. Isolates at the front of the pathogens' expanding range tend to share lineage-specific structural variants, as confirmed by short-read sequencing. We identified, for the first time, 3 new mating type A locus alleles (5 in total) and 1 new potential mating type B locus allele (3 in total). Currently, only 2 mating type combinations, A1B1 and A2B2, are known to exist outside of Colombia. A systematic survey of the M. roreri transcriptome across 2 isolates identified an expanded candidate effector pool and provided evidence that effector candidate genes unique to the Moniliophthoras are preferentially expressed during the biotrophic phase of disease. Notably, M. roreri isolates in Costa Rica carry a chromosome segment duplication that has doubled the associated gene complement and includes secreted proteins and candidate effectors. Clonal reproduction of the haploid M. roreri genome has allowed lineages with unique genome structures and compositions to dominate as it expands its range, displaying a significant founder effect.

Frosty pod rot (FPR), caused by Moniliophthora roreri (MR), is one of the most important diseases of cacao, causing losses of up to 90% in areas where it is present. A limited number of resistant clones have been identified, but additional sources of disease resistance are being sought by breeding programs. RNA sequencing was used to identify differences in constitutive (mock-inoculated) and pathogen-induced gene expression in cacao pods from trees resistant and susceptible to MR. The greatest number of differentially expressed genes (DEGs) ( n = 2,039) was identified in susceptible material following inoculation with MR (7 days postinoculation), with the most common pathway categories being cell growth/death ( n = 102), signaling ( n = 308), and carbohydrate metabolism ( n = 141). Resistant material showed little overall change in gene expression (five DEGs) following pathogen inoculation, suggesting that infection had not become established. For this reason, differences in constitutive gene expression are likely to provide the most relevant information on FPR resistance in this study. In the absence of the pathogen, 144 genes were expressed at a higher level in the resistant material and 265 at a higher level in the susceptible material. The genes most highly upregulated in resistant material encode for disease resistance proteins in the leucine-rich repeat family and cysteine-rich receptor-like protein kinases. Ten DEGs were located within quantitative trait loci on chromosomes 2, 7, and 9 previously associated with FPR resistance in this population. This study provides the first analysis of constitutively expressed genes contributing to FPR resistance and identifies significant genes within known disease resistance quantitative trait loci. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY 4.0 International license .

Phytophthora species are oomycete plant pathogens that cause great economic and ecological impacts. The Phytophthora genus includes over 180 known species, infecting a wide range of plant hosts, including crops, trees, and ornamentals. We sequenced the genomes of 31 individual Phytophthora species and 24 individual transcriptomes to study genetic relationships across the genus. De novo genome assemblies revealed variation in genome sizes, numbers of predicted genes, and in repetitive element content across the Phytophthora genus. A genus-wide comparison evaluated orthologous groups of genes. Predicted effector gene counts varied across Phytophthora species by effector family, genome size, and plant host range. Predicted numbers of apoplastic effectors increased as the host range of Phytophthora species increased. Predicted numbers of cytoplasmic effectors also increased with host range but leveled off or decreased in Phytophthora species that have enormous host ranges. With extensive sequencing across the Phytophthora genus, we now have the genomic resources to evaluate horizontal gene transfer events across the oomycetes. Using a machine-learning approach to identify horizontally transferred genes with bacterial or fungal origin, we identified 44 candidates over 36 Phytophthora species genomes. Phylogenetic reconstruction indicates that the transfers of most of these 44 candidates happened in parallel to major advances in the evolution of the oomycetes and Phytophthora spp. We conclude that the 31 genomes presented here are essential for investigating genus-wide genomic associations in genus Phytophthora. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Thread blight disease (TBD) is a serious emerging threat to global cacao production, causing leaf blight symptoms and necrosis in the canopy of infected trees. Though previously thought to be caused primarily by the fungus Marasmiellus scandens, recent work has revealed that a variety of pathogenic Marasmiaceae are responsible for TBD, with Marasmius tenuissimus identified as the dominant causal agent of TBD in Ghanaian cacao plantations. M. tenuissimus has also been reported in remote organic cacao plantations in the Amazonas Department of Peru, infecting 90 to 100% of the trees surveyed. Here, we have assembled and annotated the nuclear genome of M. tenuissimus isolate GH-37. This resource is a first attempt at a high-quality draft genome utilizing both second- and third-generation sequencing for M. tenuissimus and extends the current understanding of the evolution, ecology, and virulence of this pathogen. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

The flowering plant genus Erythroxylum contains approximately 300 species, including the economically and socially consequential crops called coca. We present the genome sequences of Erythroxylum coca and E. novogranatense, two cultigens produced for medicinal and quotidian use in the Andes and Amazon regions of South America, as well as the international cocaine industry. Sequencing was performed on an Illumina X-Ten platform, and reads were assembled by a de novo method followed by finishing via comparison with several species from the same genus. The BioProject, raw and assembled data can be accessed in GenBank for E. coca (PRJNA676123; JAJMLV000000000) and E. novogranatense (PRJNA675212; JAJKBF000000000)

We present the whole genome sequences of 56 wild Erythroxylum species from Africa, China, and the American tropics. Deep Illumina sequencing was performed on a single leaf of each voucher. We de novo assembled sequence reads and then identified and used conserved regions across all preassemblies join contigs in a finishing step. The raw and assembled data is publicly available via Genbank.

Black pod rot, caused by Phytophthora palmivora, is a devastating disease of Theobroma cacao L. (cacao) leading to huge losses for farmers and limiting chocolate industry supplies. To understand resistance responses of cacao leaves to P. palmivora, Stage 2 leaves of genotypes Imperial College Selection 1 (ICS1), Colección Castro Naranjal 51 (CCN51), and Pound7 were inoculated with zoospores and monitored for symptoms up to 48 h. Pound7 consistently showed less necrosis than ICS1 and CCN51 48 h after inoculation. RNA-Seq was carried out on samples 24 h post inoculation. A total of 24,672 expressed cacao genes were identified, and 2,521 transcripts showed induction in at least one P. palmivora-treated genotype compared to controls. There were 115 genes induced in the P. palmivora-treated samples in all three genotypes. Many of the differentially expressed genes were components of KEGG pathways important in plant defense signal perception (the plant MAPK signaling pathway, plant hormone signal transduction, and plant pathogen interactions), and plant defense metabolite biosynthesis (phenylpropanoid biosynthesis, α-linolenic acid metabolism, ethylene biosynthesis, and terpenoid backbone biosynthesis). A search of putative cacao resistance genes within the cacao transcriptome identified 89 genes with prominent leucine-rich repeat (LRR) domains, 170 protein kinases encoding genes, 210 genes with prominent NB-ARC domains, 305 lectin-related genes, and 97 cysteine-rich RK genes. We further analyzed the cacao leaf transcriptome in detail focusing on gene families-encoding proteins important in signal transduction (MAP kinases and transcription factors) and direct plant defense (Germin-like, ubiquitin-associated, lectin-related, pathogenesis-related, glutathione-S-transferases, and proteases). There was a massive reprogramming of defense gene processes in susceptible cacao leaf tissue after infection, which was restricted in the resistant genotype Pound7. Most genes induced in Pound7 were induced in ICS1/CCN51. The level of induction was not always proportional to the infection level, raising the possibility that genes are responding to infection more strongly in Pound7. There were also defense-associated genes constitutively differentially expressed at higher levels in specific genotypes, possibly providing a prepositioned defense. Many of the defense genes occur in blocks where members are constitutively expressed at different levels, and some members are induced by Ppal infection. With further study, the identified candidate genes and gene blocks may be useful as markers for breeding disease-resistant cacao genotypes against P. palmivora.

Thread blight disease has recently been described as an emerging disease on cacao (Theobroma cacao) in Ghana. In Ghana, thread blight disease is caused by multiple species of the Marasmiaceae family: Marasmius tenuissimus, M. crinis-equi, M. palmivorus, and Marasmiellus scandens. Interestingly, two additional members of the Marasmiaceae; Moniliophthora roreri (frosty pod rot) and Moniliophthora perniciosa (witches’ broom disease), are major pathogens of cacao in the Western hemisphere. It is important to accurately characterize the genetic relationships among these economically important species in support of their disease management. We used data from Illumina NGS-based genome sequencing efforts to study the mitochondrial genomes (mitogenomes) of the four cacao thread blight associated pathogens from Ghana and compared them with published mitogenomes of Mon. roreri and Mon. perniciosa. There is a remarkable interspecies variation in mitogenome size within the six cacao-associated Marasmiaceae species, ranging from 43,121 to 109,103 bp. The differences in genome lengths are primarily due to the number and lengths of introns, differences in intergenic space, and differences in the size and numbers of unidentified ORFs (uORF). Among seven M. tenuissimus mitogenomes sequenced, there is variation in size and sequence pointing to divergent evolution patterns within the species. The intronic regions show a high degree of sequence variation compared to the conserved sequences of the 14 core genes. The intronic ORFs identified, regardless of species, encode GIY-YIG or LAGLIDADG domain-containing homing endonuclease genes. Phylogenetic relationships using the 14 core proteins largely mimic the phylogenetic relationships observed in gene order patterns, grouping M. tenuissimus with M. crinis-equi, and M. palmivorus with Mon. roreri and Mon. perniciosa, leaving Mar. scandens as an outlier. The results from this study provide evidence of independent expansion/contraction events and sequence diversification in each species and establish a foundation for further exploration of the evolutionary trajectory of the fungi in Marasmiaceae family.

Economically, cacao ( Theobroma cacao L.) is a major tropical commodity for the Americas; however, severe losses due to Moniliophthora roreri (Cif. and Par.), which causes frosty pod rot (FPR), and Phytophthora spp., which causes black pod rot (BPR), have reduced cacao production in the Americas. The objectives of this study are to (i) re-confirm the QTL using different marker set; (ii) discover new QTL associated with FPR and BPR resistance using SNP markers; and (iii) find genes in the candidate QTL regions. At CATIE in Turrialba, Costa Rica, an F 1 mapping population of cacao was obtained by crossing “POUND 7,” a clone moderately susceptible to FPR and resistant to BPR, with “UF 273,” resistant to FPR and highly susceptible to BPR. A total of 179 F 1 progeny were fingerprinted with 5149 SNP markers and a dense linkage map composed of 10 linkage groups was developed using 2910 polymorphic SNP markers. Also segregating F 1 trees were screened for resistance to FPR and BPR diseases. Seven QTL previously reported on chromosomes 2, 7, and 8 for FPR resistance and on chromosomes 4, 8, and 10 for BPR resistance were confirmed. Additionally, eight QTL were identified for FPR resistance (chromosomes 4, 9, and 10) and BPR resistance (chromosome 2). The expression of genes commonly associated with plant defense and disease resistance that are located within the identified QTL was confirmed.

Frosty pod rot (FPR) (Moniliophthora roreri), along with black pod rot (Phytophthora species) and witches' broom disease (Moniliophthora perniciosa) constitute the main phytosanitary problems limiting cacao (Theobroma cacao) production causing severe yield losses. One of the main sought after methods of pod rot management is the selection of tolerant cacao genotypes. Typically, the selection is carried out through the quantification of the percentage of diseased pods (PDP). However, PDP does not consider the relative productivity, or production potential (PT) of the genotype. Production potential can vary among cacao genotypes. Consequently, genotypes with similar PT can have similar or vastly different disease tolerance levels as measured by PDP. The disease and production index (DPI) was developed to integrate a genotype's tolerance to M. roreri and other diseases as measured by PDP, with its PT. Here, we evaluated the number of healthy pods, number of diseased pods, and weight of fresh seed for 29 clones grown in replicated five-tree plots over 4 years. The data obtained was used to calculate PDP and DPI for each clone for three different disease combinations: frosty pod rot alone, pod rots other than frosty pod rot, and the combination of all pod rots. Multivariate analysis verified that DPI discriminated between clones based on productivity and disease tolerance. Surprisingly, there was a close ranking of clones between resistance to FPR and resistance to all other pod rots. The DPI can be used in breeding programmes focused on the selection of high yielding disease-tolerant cacao genotypes.

Theobroma cacao (chocolate tree) is currently under serious threat from thread blight disease (TBD), which has been attributed to the causal agent Marasmiellus scandens in other regions of the world. TBD in Ghana has similar symptomology but variable signs. This study sought to determine whether TBD in Ghana was caused by a single agent and whether Marasmiellus scandens was a significant agent of TBD. Forty-eight isolates were collected from eight geographical locations in Ghana for morphological and molecular characterization. Disease signs occurred as vegetative rhizomorphs or hyphal aggregates, which were classified into five morphotypes: A, abundant thin, black, “horse hair”-type rhizomorphs; B, scattered brown rhizomorphs; C, whitish to brownish-white; D, faint cream or dull white; and E, aggregates of shiny or silky white hyphae. Sequencing and analyses of three loci—the internal transcribed spacer region of the nuclear ribosomal repeat, nuclear large subunit, and mitochondrial small subunit—detected four species, all members of the Marasmiaceae, causing TBD-like disease. These were identified as Marasmius crinis-equi (morphotype A), Marasmius tenuissimus (morphotypes B and C), Marasmiellus palmivorus (morphotype E), and Marasmiellus scandens (morphotype D). Marasmius tenuissimus, the most frequently isolated TBD fungus in this study, is primarily an Asian fungus and not previously associated with diseases of cacao. Marasmiellus palmivorus, the second most frequently isolated fungus, is a pan-tropical pathogen with a broad host range; this is the first report of the fungus causing TBD on cacao. Marasmius crinis-equi also has a broad pan-tropical distribution and host range and causes thread blight on several tropical tree crops. Surprisingly, Marasmiellus scandens, the most frequently cited agent of TBD in cacao, made up only 8% of the isolates.