Molecular and characterization of cassava zinc finger-homeodomain (ZF-HD) transcription factors reveals their role in disease resistance

Abstract

The production of cassava (Manihot esculenta Crantz) is constantly threatened by cassava bacterial blight (CBB), caused by Xanthomonas phaseoli pv. manihotis (Xpm). Zinc finger homeodomain (ZF-HD) belongs to a family of homozygous heterotypic cassette genes widely implicated in various developmental and physiological processes in plants. Despite their importance, a comprehensive analysis of ZF-HD genes, particularly those involved in disease resistance, has not been performed for cassava. In the present study, we utilized bioinformatics methods to identify 21 ZF-HD genes distributed across 11 chromosomes of cassava genome, with the majority exhibiting gene structure without introns. Phylogenetic analysis categorized these genes into two major groups (MIF and ZHD) with five subgroups. We observed fourteen pairs of duplicated genes, suggesting that segmental duplication has likely facilitated the expansion of the cassava ZF-HD gene family. Comparative orthologous analyses between cassava and other plant species shed light on the evolutionary trajectory of this gene family. Promoter analyses revealed multiple hormone- and stress-related elements, indicative of a functional role in stress responses. Expression profiling through RNA-seq and qRT-PCR demonstrated that certain cassava ZF-HD genes are up-regulated in response to Xpm11 infection, suggesting their involvement in defense mechanisms. Notably, MeZHD7 gene was identified via virus induced gene silencing (VIGS) as potentially crucial in conferring resistance against CBB. Results from subcellular localization experiments indicated that MeZHD7 was localized in the nucleus. The Luciferase reporter assay demonstrated an interaction between MeZHD7 and MeMIF5. These findings may lay the foundation for further cloning and functional analyses of cassava ZF-HD genes, particularly those associated with pathogen resistance.