
==== Front
Parasit Vectors
Parasit Vectors
Parasites & Vectors
1756-3305
BioMed Central London

6477
10.1186/s13071-024-06477-1
Research
Identification and function characterization of NcAP2XII-4 in Neospora caninum
Nan Huizhu 123
Lu Xin 12
Zhang Chao 12
Yang Xin 12
Ying Zhu 4
Ma Lei lmahappy@hebtu.edu.cn

123
1 https://ror.org/004rbbw49 grid.256884.5 0000 0004 0605 1239 Ministry of Education Key Laboratory of Molecular and Cellular Biology, College of Life Sciences, Hebei Normal University, Shijiazhuang, 050024 Hebei China
2 Hebei Collaborative Innovation Center for Eco-Environment, Shijiazhuang, Hebei China
3 https://ror.org/004rbbw49 grid.256884.5 0000 0004 0605 1239 Hebei Key Laboratory of Animal Physiology, Biochemistry and Molecular Biology, College of Life Sciences, Hebei Normal University, Shijiazhuang, Hebei China
4 https://ror.org/04v3ywz14 grid.22935.3f 0000 0004 0530 8290 National Animal Protozoa Laboratory, College of Veterinary Medicine, China Agricultural University, Beijing, 100193 China
14 9 2024
14 9 2024
2024
17 39223 5 2024
3 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
Background

Neospora caninum is a protozoan parasite in the Apicomplexa controlled by complex signaling pathways. Transcriptional control, an important way to regulate gene expression, has been almost absent in the N. caninum life process. However, to date, research on the transcriptional regulation of the AP2 family factors in N. caninum has been extremely limited. A prior study demonstrated that removing rhoptry protein 5 (ROP5), a significant virulence factor, resulted in abnormal expression levels of predicted NcAP2XII-4 in N. caninum, suggesting that the factor may regulate the function of ROP5. This study aimed to identify NcAP2XII-4 and its function in transcriptional regulation.

Methods

The NcAP2XII-4 gene was identified by analyzing the N. caninum genome. A polyclonal antibody against the protein was prepared and purified, and its expression and localization in the parasite were detected using western blot (WB) and immunofluorescence assay (IFA). The ΔNcAP2XII-4 strain was constructed from the Nc1 strain using CRISPR/Cas9 to study its effect on the growth and development of N. caninum, and DAP-Seq and electrophoretic mobility shift assay (EMSA) were used to verify the transcriptional regulatory functions of the gene.

Results

Bioinformatic analysis showed that NcAP2XII-4 consists of 11,976 bp and encodes 3991 amino acids, with a predicted molecular mass of 410 kDa. The protein has two AP2 domains, 1207aa-1251aa and 3453aa-3500aa, and is predicted to be located in the nucleus. The results of PCR, WB, and IFA were in accordance with the bioinformatics analysis. ΔNcAP2XII-4 was successfully constructed, but the strain could not be released and ultimately succumbed within parasitophorous vacuoles (PVs). Plaque assays demonstrated that parasites lacking this gene could not form plaques. One motif was successfully identified using DAP-Seq technique. Two prokaryotic expression vectors containing the AP2 domain of NcAP2XII-4 were successfully constructed, and two prokaryotic expression proteins, AP2-D1 and AP2-D2, and ROP5 biotinylated probes were prepared. Using EMSA, NcAP2XII-4 was shown to regulate ROP5 transcription by binding to its promoter.

Conclusions

NcAP2XII-4 is an essential gene in N. caninum. This study provides a foundation for further research on transcriptional regulation in N. caninum and identifies a new candidate factor for the development of vaccines against N. caninum.

Graphical Abstract

Keywords

Neospora caninum
NcAP2XII-4
DAP-Seq
Transcription factor
http://dx.doi.org/10.13039/501100001809 National Natural Science Foundation of China 31902292 Ma Lei http://dx.doi.org/10.13039/501100015286 Hebei Provincial Key Research Projects 22326703D Ma Lei issue-copyright-statement© BioMed Central Ltd., part of Springer Nature 2024
==== Body
pmcBackground

Neosporosis, caused by the protozoan parasite Neospora caninum, poses a significant threat to animal husbandry. This apicomplexan intracellular parasite primarily infects canids as its definitive host and a variety of other animals as intermediate hosts [1]. The disease manifests through symptoms, such as abortion, stillbirth, mummification, and neurological dysfunction [2]. Despite the severity of neosporosis, effective vaccines for N. caninum remain elusive. Inactivated vaccines, recombinant protein vaccines, and live vaccines [3–5] have not been used in the field. To combat this disease, further research is needed.

AP2 transcription family factors are important for plant growth, reproduction and interactions with the environment. Balaji et al. [6] found a domain in Apicomplexan protozoa that was similar to the AP2 superfamily protein in plants. The corresponding genes were subsequently named Apicomplexan AP2 (ApiAP2). Tracey et al. [7] comprehensively surveyed the DNA-binding specificities of 27 ApiAP2 protein families in Plasmodium falciparum and demonstrated that ApiAP2 proteins modulate key regulatory decisions throughout parasite development. For example, AP2-O is involved in biological processes including morphogenesis, locomotion, midgut penetration, protection against mosquito immunity, and preparation for subsequent oocyst development [8]. AP2-I is a key regulator of red blood cell invasion [9], while AP2-FG regulates female-specific gene expression [10]. AP2-G5 is essential to gametocyte maturation [11], and AP2-G triggers the initiation of differentiation and early development of gametocytes in the late trophozoite stage [12]. Moreover, AP2-O, AP2-SP, AP2-HC and AP2-G2 were studied in Plasmodium [13, 14]. Timothy et al. [15] described a novel set of antimalarial compounds that interact with the AP2 DNA-binding domains. ApiAP2 factors have been found in other apicomplexans and have been verified to bind to DNA containing specific motifs, acting as transcription factors. In Toxoplasma, 67 ApiAP2 proteins have been identified [16, 17]. Many studies have been conducted on the function of many ApiAP2 transcription factors in Toxoplasma gondii [18–20]. Many factors are related to parasite virulence; cell division cycle; development of tachyzoites, bradyzoites, and cysts; and sexual stage regulation of the parasite [17, 19, 21–25]. Fifty-three ApiAP2 transcription factors have been predicted in Eimeria tenella [26], but their functions have not been studied in depth. Members of the ApiAP2 transcription factor family play important roles in regulating life cycle transitions and other aspects of parasitic biology. Understanding the intricate roles of ApiAP2 factors in apicomplexan parasites holds great promise for advancing our knowledge of parasitic biology and developing effective interventions.

Deletion of the key virulence factor, rhoptry protein 5 (ROP5), leads to abnormal expression levels of predicted NcAP2XII-4, NcAP2XI-5, and NcAP2VII-2 in N. caninum, suggesting that these factors are related to the regulation of the parasite life cycle [27]. The function of these factors, especially in NcAP2XII-4, has been rarely studied in N. caninum, even across all parasites. Therefore, this study attempted to identify and explore the function of N. caninum NcAP2XII-4. This will lay the foundation for further research on the transcriptional regulation of N. caninum by ApiAP2 proteins and provide new directions for the development of anti-N. caninum drugs.

Methods

Parasite strains and cell culture

The N. caninum Nc1 strain used in this study was a gift from Professor Qun Liu (College of Veterinary Medicine, China Agricultural University) and stored at the College of Life Sciences, Hebei Normal University. All parasites were maintained in confluent monolayers of human foreskin fibroblasts (HFFs). HFF cells were cultured in DMEM supplemented with 20% fetal bovine serum (FBS). Neospora caninum tachyzoites were inoculated into cells and cultured with DMEM supplemented with 2% FBS at 37 ºC with 5% CO2.

Bioinformatic analysis of NcAP2XII-4

The protein sequence of NCLIV_063920 (NcAP2XII-4) was obtained from the ToxoDB (http://toxodb.org/toxo/) website. The physical and chemical properties of N. caninum NcAP2XII-4 were obtained using the online tool ExPASy-Protparam (https://web.expasy.org/protparam/), subcellular localization by Cello v2.5, transmembrane structure by TMHMM 2.0, signal peptide by SignalP-4.1 Server, the secondary structure by SOPMA, the tertiary structure by SWISS-MODEL, and phosphorylation site by Netphos 3.1 Server.

Amplification of the NcAP2XII-4 gene

Total RNA was extracted from 1 × 108 tachyzoites of N. caninum Nc1 using TRIzol reagent and converted to cDNA using an EasyScript First-Strand cDNA Synthesis SuperMix kit (TransGen, Beijing, China), in accordance with the manufacturer’s instructions. The open reading frame was amplified by PCR using the primers NcAP2XII-4-F/R (F: 5′-ATGGAGCACTGGGACTTGGGCGAGGAC-3′/R: 5′-GTGGAGGTCCAGGGTTGCAGATGCG-3′). The product was cloned into the pEASY®-Blunt Cloning Vector (TransGen Biotech, Beijing, China), and the recombinant plasmid was identified by PCR and sequencing.

Preparation of NcAP2XII-4 polypeptide antibody

One antigenic sequence containing 12 amino acids (C-PETSRAGSSGVA) was designed for NcAP2XII-4 polypeptides. The polypeptide was synthesized using the typical solid-phase synthesis and qualified by HPLC–MS to ensure 95% purity. The polypeptide was conjugated to keyhole limpet hemocyanin (KLH) at the N-terminus. New Zealand White rabbits were immunized with NcAP2XII-4-KLH conjugates to prepare anti-NcAP2XII-4 polypeptide antibodies. Two rabbits were each vaccinated with a dose of 500 µg polypeptide per session by a single subcutaneous injection at 2-week intervals. Serum was obtained after five immunizations and stored at −80 ºC after purification.

Western blot (WB) and immunofluorescence assay (IFA)

For WB, 1 × 108 tachyzoites were collected, and the nuclear protein was extracted using the ProteinExt® Mammalian Nuclear and Cytoplasmic Protein Extraction Kit (TransGen Biotech, Beijing, China), according to the manufacturer’s instructions. The proteins were separated using an 8% SDS-PAGE and then transferred onto a polyvinylidene fluoride (PVDF) membrane. The specific primary antibody, anti-NcAP2XII-4, was added at a dilution of 1:100 and incubated overnight at 4 ºC, followed by incubation with HRP-labeled goat anti-rabbit IgG (TransGen Biotech, Beijing, China). The bands were visualized using a highly sensitive ECL luminescence reagent (Sangon Biotech, Shanghai, China). IFA was used to detect the NcAP2XII-4 subcellular localization in N. caninum, as previously reported [27]. The difference was the use of primary and secondary antibodies. The primary antibodies used were rabbit anti-NcAP2XII-4 (1:100) and mouse anti-NcACP (1:100). The secondary antibodies used were Cy3-labeled donkey anti-rabbit IgG (1:100, TransGen Biotech, Beijing, China) and FITC-labeled goat anti-mouse IgG (1:100, TransGen Biotech, Beijing, China).

Phylogenetic analysis

Nucleotide sequences were subjected to correction before BLAST alignment (http://toxodb.org/toxo/ and https://blast.ncbi.nlm.nih.gov/Blast.cgi). All sequences downloaded from GenBank were aligned using MEGA 6.06 (https://www.megasoftware.net) and edited manually (available on request). Bayesian analysis was performed using MrBayes 3.0 (http://nbisweden.github.io/MrBayes/) with four MCMC strands, 100,000 generations, and an initial burn-in of 1000. Four categories of among-site rate variations were used, and trees were sampled every ten generations.

Construction of the NcAP2XII-4 knockout N. caninum strain

The CRISPR/Cas9 system was used to generate NcAP2XII-4 deficiency (∆NcAP2XII-4 parasites) as previously described [28]. The EuPaGDT Library in ToxoDB (https://toxodb.org) was used to design the gRNA targeting sites of the plasmids pCRISPR-CAS9-AP2XII-4. The Cas9 upstream and downstream fragments containing gRNA sequences were amplified and ligated using a multi-fragment connection kit (TransGen Biotech, Beijing, China). To disrupt the NcAP2XII-4 locus, the chloramphenicol resistance gene (CAT) was inserted into the 3′ and 5′ flanks of the NcAP2XII-4 region and ligated into the plasmid backbone carrying ampicillin resistance. To construct homologous recombinant plasmids of NcAP2XII-4 (pAP2XII-4-CAT), the 3′ flanking and 5′ flanking sequences of the NcAP2XII-4 gene were amplified from the genomic DNA of Nc1 parasites. The pCRISPR-CAS9-AP2XII-4 and pAP2XII-4-CAT plasmids were co-transfected into Nc1 parasites. Thereafter, the parasites were cultured and screened using chloramphenicol. The resulting strain was named ΔNcAP2XII-4.

Plaque assay

Experiments were performed as previously reported [27]. HFFs had been previously cultured in six-well plates. Subsequently, 300 Nc1 tachyzoites per well were inoculated into the cells. Concurrently, all ΔNcAP2XII-4 tachyzoites, which had been transfected with a knockout plasmid, were infected into the cells. The plate was incubated for 7 days in a 37 ºC incubator with 5% CO2, fixed with 4% paraformaldehyde for 20 min, and stained with 2% crystal violet staining (Sangon Biotech, Shanghai, China). The stained wells were visualized under a fluorescence microscope (Leica Microsystems, Germany). The plaque area measurement was performed as previously described [29].

Screening of NcAP2XII-4 binding sites

DAP-seq was performed as previously described [30]. Tachyzoites (1 × 109) were collected and purified, and then genomic DNA was extracted using a genomic DNA kit (TransGen Biotech, Beijing, China). A DNA library was constructed by fragmentation, recovery, end repair, and joining in turns. The pDAP-Halo-AP2XII-4 expression vector was constructed and subsequently expressed in a wheat germ expression system provided by Promega. Upon binding the NcAP2XII-4 protein to the prepared DNA library, high-throughput sequencing was performed using an Illumina HiSeq (Zoonbio, Jiangsu, China). The data obtained were aligned with the N. caninum genome using BWA-MEM. All reliable peaks were obtained in the two independent biological replicates by combining both the MACS call peak and IDR software. The peaks were annotated using HOMER software. The conserved motifs were analyzed using MEME software. The frequency distribution of reads near the TSS was analyzed using depTools2 software. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed on genes associated with the peaks using the Gene Ontology Resource (https://geneontology.org/) and KEGG (http://www.genome.jp/).

Electrophoretic mobility shift assay (EMSA)

To validate the binding of the NcAP2XII-4 to the ROP5 gene promoter, two sequences of the AP2 domain of NcAP2XII-4, AP2-D1, and AP2-D2 were amplified from N. caninum cDNA. The sequenced AP2-D1 and AP2-D2 sequences were subcloned separately into the pET30a expression vector containing His and Halo-tag to create the pET-Halo-AP2-D1 and pET-Halotag-AP2-D2 plasmids. The correct pET-Halo-AP2-D1 and pET-Halotag-AP2-D2 plasmids were transformed into chemically competent Transetta (DE3) cells (TransGen Biotech, Beijing, China) for protein production. Protein expression was then induced with 0.8 mM IPTG for 4 h at 37 °C. Following the protocol of the Ni–NTA resin (TransGen Biotech, Beijing, China), proteins were extracted and purified after the bacteria had been harvested. Subsequently, they were detected by SDS-PAGE. Biotin-labeled DNA probes measuring 28 bp in length (probe F:5′-GGCACTTCTTCAACTGACAACGTAATCA-3′, probe R:5′-TGATTACGTTGTCAGTTGAAGAAGTGCC-3′, probe F (cold): 5′- GGCACTTCTTCAACTGACAACGTAATCA-3′) were designed and synthesized for use in EMSA. EMSA reactions were performed according to the manufacturer’s instructions for the Lightshift Chemiluminescent EMSA Kit (Thermo Scientific) [31]. Finally, the signals were detected using a fluorescence/chemiluminescence imaging system (ChemiScope 6000Pro; CLiNX, Shanghai, China).

Results

Bioinformatics analysis of NcAP2XII-4

Bioinformatic analysis showed that the full-length sequence of N. caninum NcAP2XII-4 was 12,259 bp, including one intron and two exons. The length of coding gene was 11,976 bp, which encoded 3991 amino acids. The molecular weight of the protein was approximately 410 kDa, with an instability coefficient of 55.69 and hydrophilicity of – 1.063. These characteristics indicate that it is a hydrophilic protein with an isoelectric point of 7.06. This protein did not contain any signal peptides (Fig. 1a) or transmembrane regions (Fig. 1c). The phosphorylation sites of a protein play important roles in predicting its potential function. The predicted analysis indicated that there were 171 phosphorylation sites, consisting of 106 serine sites, 64 threonine sites, and one tyrosine site. Subcellular localization analysis suggested that the protein was located in the nucleus with a reliable index of 3.731 and high credibility. Sequence alignment and modeling analysis showed that N. caninum NcAP2XII-4 contains two AP2 domains, 1207aa-1251aa and 3453aa-3500aa, comprising one alpha-helix and three beta-folded structures in each domain (Fig. 1b). These structures exhibit similarities to the classical model PF140633 AP2 domain. The two AP2 domains are linked through an alpha helix to bind with DNA and exert transcriptional regulatory function, suggesting that N. caninum NcAP2XII-4 could play a certain role in the transcriptional regulation in the life cycle of the parasite.Fig. 1 Bioinformatic analysis of Neospora caninum transcription factor NcAP2XII-4. a Prediction of the signal peptide structure of NcAP2XII-4. b Prediction of the domain structure of NcAP2XII-4. c Prediction of the transmembrane region of the NcAP2XII-4

Identification of N. caninum NcAP2XII-4

Three fragments were successfully amplified using N. caninum cDNA as a template. Additionally, a fragment of approximately 12,000 bp was amplified, and the fragment matched the predicted gene length of NcAP2XII-4 (Fig. 2a). DNA sequence assays showed that the full-length N. caninum NcAP2XII-4 was successfully obtained, and the sequence exhibited 100% identity with the known sequence. According to the results of the phylogenetic analysis, N. caninum NcAP2XII-4 exhibits a greater amino acid sequence identity with Besnoitia besnoitai, Cystoisospora suis, and Hammondia hammondi (Fig. 2b). Western blot showed that the polypeptide antibody could specifically recognize the natural NcAP2XII-4 in the parasite, with a size of approximately 400 KDa, which was consistent with the expected size, indicating the good specificity of anti-NcAP2XII-4 antibody (Fig. 2c). The IFA test showed that NcAP2XII-4 was distributed in the nucleus, consistent with the predicted result (Fig. 2d). These results suggested that N. caninum NcAP2XII-4 was successfully identified.Fig. 2 Identification of NcAP2XII-4. a Amplification of NcAP2XII-4 gene. b Phylogenetic analysis of NcAP2XII-4 gene evolution in the apicomplexan classes. c Western blot results of the NcAP2XII-4 antibody. d Immunofluorescence localization of NcAP2XII-4 protein in Neospora caninum. Scale bar, 5 µm

NcAP2XII-4 is an essential gene in N. caninum

To investigate the function of N. caninum NcAP2XII-4, parasites deficient in NcAP2XII-4 were generated using the following method as illustrated in Fig. 3a. Using CRISPR/Cas9 technology, a pCRISPR/Cas9-AP2XII-4 plasmid with CAT was constructed and co-transfected into Nc1 with the donor DNA. Based on the preliminary evidence of specific recombination bands of 1940 bp and 2350 bp by PCR, the NcAP2XII-4 gene was disrupted in N. caninum (Fig. 3b). The results of IFA indicated that Nc1 strain had an evident NcAP2XII-4 signal, while the signal was undetectable in the NcAP2XII-4 KO strain (Fig. 3c). It is suggested that the parasite loss of NcAP2XII-4, which is also known as ΔNcAP2XII-4, was successful. It was discovered during the culture process of parasites. The PVs of ΔNcAP2XII-4 were formed in the host cells; however, these PVs were significantly smaller than those of Nc1, and the shapes of the parasites were abnormal (Fig. 3c). After 5 days of continuous culturing with CAT medium, it was confirmed that all parasites remained unreleased and eventually died. The plaque assay showed that Nc1 formed normal plaques, whereas ΔNcAP2XII-4 failed to produce plaques (Fig. 3d, e). Hence, we infer that NcAP2XII-4 is an essential gene for maintaining normal life activities in N. caninum.Fig. 3 Construction and identification of NcAP2XII-4 KO parasites. a Schematic diagram of NcAP2XII-4 KO parasite construction. b Gene identification of NcAP2XII-4 KO parasites. c Analysis of IFA results. The vacuole formation of Nc1 parasites with the signal of AP2XII-4 (red) was normal, whereas NcAP2XII-4 KO parasites without the signal of AP2XII-4 were abnormal. d Plaque assay comparing growth of Nc1 and NcAP2XII-4 KO parasites. e The plaque areas were counted by randomly choosing 30 plaques and using the pixel point; the data were compiled from three independent experiments. “Pi” indicates pixels of each plaque

Screening for the binding site of NcAP2XII-4

Genomic DNA of N. caninum was successfully extracted (Fig. 4a) and fragmented (Fig. 4b). The DNA fragment containing the AP2 domain of the NcAP2XII-4 gene was amplified to a size of approximately 2400 bp (Fig. 4c). The target gene was successfully connected to the cell-free expression vector pDAP, resulting in the construction of the expression vector pDAP-AP2XII-4, which was confirmed by sequencing (Fig. 4d). Western blot showed that a 129-kDa HaloTag-labeled target protein was successfully obtained (Fig. 4e).Fig. 4 Gene library and pDAP-AP2XII–4 expression vector. a Genomic DNA extraction. b Fragmentation of genomic DNA. c Amplification of NcAP2XII-4 gene. d Construction map of the expression vector. The green region represents the APXII-4 insertion gene. e Western blot analysis of the protein expression vector. M: Marker; 1: target protein

The entire gene functional region is divided into promoter, exon, intron, transcription start site (TSS), transcription termination site (TTS), and intergenic, but transcription factor-binding sites are generally enriched near the TSS region. Accordingly, our attention was directed towards the peaks that had been annotated as promoters and TSS. Two independent biological replicates, replicate 1 and replicate 2, identified 83 and 82 peaks with fold enrichment > 2.0, respectively, with 42 overlapping peaks (Fig. 5a). A binding motif sequence, HCRAAASYGAWK, was identified using MME software analysis (Fig. 5b). The effective peaks were distributed both upstream and downstream of TSS (Fig. 5c) and were located on chromosomes with different numbers (Fig. 5d, e). Following the annotation of the peaks, they were mainly distributed in the intergenic regions (57%) and promoters (30%).Fig. 5 DAP-seq data analysis. a Venn diagram showing the number of peaks identified in Replicate 1 and Replicate 2. b Sequences enriched around the summits of AP2XII-4 peaks. c Distribution of peaks near the transcription start sites. d Annotation statistics of peaks in the genomic region. Pie chart showing the functional categories of target genes (123 genes in total). Hypothetical protein-encoding genes are not shown. The number of members in each group is shown in the chart. e Distribution of peaks on chromosomes, with the height indicating the relative mass of the peaks

Identification of the potential biological process and signaling pathways of peak-related genes

To date, the function of NcAP2XII-4 in N. caninum has remained unclear. According to the GO and KEGG annotation data, all N. caninum genes were chosen as the background, and all genes (P < 0.05) related to all peaks were obtained to estimate the differential GO and KEGG pathways. GO can be divided into three parts, molecular function, biological process, and cellular components. Among these, the target genes, regulated by NcAP2XII-4, involved in biological processes were the most common, followed by those involved in cellular components and molecular function in N. caninum (Fig. 6a). The bubble chart analysis revealed that the genes that differ the most are those that are involved in the regulation of proteasomal protein catabolic process, sperm midpiece, and ubiquitin protein ligase binding, corresponding to the biological process, cellular component, and molecular function (Fig. 6b). Sperm midpieces may be associated with proliferation. KEGG is a database containing information on biochemical reactions, signaling pathways, metabolic pathways, and biological processes. KEGG pathway enrichment analysis was performed on the differentially expressed genes, and the results showed that the significantly differentially expressed genes were mainly enriched in the pathways of AGE-RAGE signaling pathway in diabetic complications, proteasome, aminoacyl-tRNA biosynthesis, purine metabolism, and oxidative phosphorylation (Fig. 6c, d).Fig. 6 GO and KEGG analyses of peak-associated genes. a Bar chart analysis of the GO of peak-associated genes. b Bubble chart analysis of the GO enrichment of the top 20 peak-associated genes corresponding to the promoter. c Bar chart of the KEGG analysis of peak-associated genes. d Bubble chart analysis of the KEGG enrichment of the top 20 peak-associated genes corresponding to the promoter

NcAP2XII-4 binding ROP5 promoter in N. caninum

Previous studies showed that the expression of NcAP2XII-4 is abnormal after deletion of ROP5 [27], suggesting a potential interaction between NcAP2XII-4 and ROP5. Analysis of DAP-seq data revealed that regulation of the ROP5 promoter was found among the NcAP2XII-4 regulatory factors, indicating that NcAP2XII-4 might bind to a motif in the ROP5 promoter. To verify this, we prepared a biotin-labeled ROP5 probe and two NcAP2XII-4 proteins containing the AP2 domain, AP2XII-4-D1 (Fig. 7a), and AP2XII-4-D2 (Fig. 7c), which was successfully verified by WB, respectively (Fig. 7b, d), to perform EMSA experiments. The results showed that the target bands appeared as expected in the lanes of probe + AP2XII-4-D1 and probe + AP2XII-4-D2 but were negative in the lanes of probe (cold) + AP2XII-4-D1 and probe (cold) + AP2XII-4-D2, suggesting the specificity of the probe (Fig. 7e). In other words, the probe specifically combined AP2XII-4-D1 and AP2XII-4-D2 proteins, and there was competition between the probe and the cold probe, implying that NcAP2XII-4 could bind to the ROP5 promoter and regulate the function of ROP5 in N. caninum.Fig. 7 NcAP2XII-4 binding to the ROP5 promoter using EMSA assay. a Analysis of the AP2-D1 protein by SDS-PAGE. b Analysis of AP2-D1 protein by WB. c Analysis of the AP2-D2 protein by SDS-PAGE. d Analysis of AP2-D2 protein by WB. M: Marker. e EMSA assay. 1: Positive probe; 2: positive probe + positive protein; 3: probe; 4: probe + protein AP2XII-4-D1; 5: probe + protein AP2XII-4-D1 + probe (cold); 6: probe 1 + protein AP2XII-4-D2; 7: probe + protein AP2XII-4-D2 + probe (cold); 8: protein AP2XII-4-D1; 9: protein AP2XII-4-D2; 10: probe + HaloTag-labeled protein. *Binding protein

Discussion

Neospora caninum is a protozoan that parasitizes animal cells and poses a significant threat to dairy cattle. Despite its importance, there are few reports on the transcriptional regulation of N. caninum. One notable report indicated that the deletion of ROP5, a key virulence factor of N. caninum, increases the transcription level of NcAP2XII-4 [27]. Building on this, our study successfully identified N. caninum NcAP2XII-4 at the gene, protein, and cell levels. We confirmed its transcriptional regulation function and described it as an essential gene of the parasite.

The AP2/ERF family of transcription factors is related to biological development and stress resistance. The N. caninum NcAP2XII-4 gene (NCLIV_063920) was discovered in the ToxoDB database during this study. Structural predictive analysis revealed two domains similar to PF14-0633. However, the predicted protein size of approximately 410 KDa presents a challenge to expressing the full-length protein and preparing the antibody. After numerous experiments, a polypeptide antibody against NcAP2XII-4 was successfully prepared and employed to the detect NcAP2XII-4 in N. caninum. Historically, in most cases, the full-length protein or fragment was selected as an antigen for as long as feasible. The fusion of numerous peptides could be utilized to prepare antibodies according to the results of this study. This method provides convenience in preparing antibodies against difficult antigens, such as those with complex structures, large molecular weights, and low codon usage frequencies.

In recent years, the interaction between proteins and DNA has been widely used to study transcriptional regulation in model organisms. Previously, the combination of ChIP-seq and RNA-seq was utilized to explore the interaction. However, this technology necessitates high-quality antibodies, which presents a challenge when studying proteins with large molecular weights and obtaining high-quality antibodies. DAP-seq not only has the advantages of ChIP-seq and RNA-seq but also does not require the preparation of high-quality antibodies and even has higher sensitivity than ChIP-seq and RNA-seq. DAP-seq is particularly well suited for studying the interactions between DNA and proteins that are difficult to obtain. Currently, DAP-seq is widely used in the study of protein and DNA interactions in model plants, such as maize and arabidopsis [32]. Ning Lihua et al. [33] applied DAP-seq technology to detect 317 maize transcription factors that combined with target genes and played a role in transcription regulation. Among them, 97 target genes were detected and reported by ChIP-seq, and 220 target genes were detected by DAP-seq. The results are highly reliable. He et al. [34] investigated the transcriptional regulation function of abscisic acid and synthetic auxin 1-naphthyl acetic acid (NAA) on grape norisoprenoid. They combined weighted gene co-expression network analysis (WGCNA) and DAP-seq and identified regulatory factors of norisoprenoid production. Zhu et al. [35] identified 52 key direct targets of OsbZIP09 using RNA-seq and DAP-seq analyses, including OsLOX2 and late embryogenesis-rich (LEA) family genes involved in controlling seed germination in rice. D'Inca et al. [36] explored the transcriptional regulatory mechanism of grape VviNAC33 in terminating photosynthetic activity and organ growth using the DAP-seq method.

DAP-seq has been extensively applied to investigate transcription factors in maize, Arabidopsis thaliana, poplar, strain, Chinese pine, and human, with maize and rice being the most prominent examples. Furthermore, DNA-seq technology has been employed to investigate the function of bacterial proteins. Trouillon et al. [37] identified 39 genes that are mainly involved in the metabolism of Pseudomonas aeruginosa and are regulated by eight transcription factors using DAP-seq technology. The success rate of the DAP-seq experiment is also related to the type of transcription factor; for example, the success rate of bZIP and NAC transcription factor families is obviously higher than that of bHLH and MADs-box families [35, 38]. ChIP-seq and RNA-seq have been widely used in parasite research in the past. Although DAP-seq technology has many advantages, it has not been applied to research on the protein function of parasites. In this study, we investigated the transcriptional regulation of NcAP2XII-4 by using DAP-Seq. Notably, the motif of NcAP2XII-4 has been successfully identified with good conservation and high reliability, thus confirming the feasibility of the DAP-seq technique in the study of transcription factors of N. caninum. This has provided a reference method for the study of transcriptional regulation of N. caninum and other apicomplexans. The motif obtained from the above analysis was analyzed, and its binding protein contained rhoptry protein family factors. Previous studies have shown that the deletion of ROP5, which is a key virulence factor of N. caninum, causes an increase in the transcription level of NcAP2XII-4. Thus, we have prepared a protein containing the AP2 domain in N. caninum NcAP2XII-4 and conducted the EMSA test with a ROP5 biotin-labeled probe. The results indicate that NcAP2XII-4 is capable of binding to the promoter of ROP5 and exerting its transcriptional activity. This not only demonstrates the reliability of the results of the DAP-seq technique but also further demonstrates the feasibility of the technique in the field of parasites. This study has filled a significant gap in the research field of transcription regulation factors of the ApiAP2 family of N. caninum. It has laid a foundation for the research on the transcription regulation function of N. caninum and provided new candidate factors for the research and development of N. caninum vaccine.

Conclusion

The transcription factors of N. caninum were identified at the gene, protein, and cell levels of NcAP2XII-4, and it was determined that the gene is an essential gene for N. caninum to complete its life cycle. The gene primarily functions in transcription regulation in N. caninum and regulates the activity of ROP5.

Acknowledgements

We are grateful to Prof. Qun Liu (National Animal Protozoa Laboratory, College of Veterinary Medicine, China Agricultural University, Beijing, China) for gifting the parasites, CRISPR/Cas9 plasmid, and mouse anti-NcSRS2 antibody and Jing Liu (National Animal Protozoa Laboratory, College of Veterinary Medicine, China Agricultural University, Beijing, China) for providing technical guidance.

Author contributions

LM, HN and XL contributed to the study design. HN, XL, XY, CZ and ZY contributed to collect and analyze data. HN and LM contributed to the manuscript development. All authors provided critical revisions and took responsibility for the work. All the authors have read and approved the final manuscript.

Funding

This work was supported by National Natural Science Foundation of China (No. 31902292) and the Hebei Provincial Key Research Projects, China (No. 22326703D). The funding sources had no role in the study design, data collection, analysis, interpretation, or writing of this report.

Availability of data and materials

No datasets were generated or analyzed during the current study.

Declarations

Ethics approval and consent to participate

This study complied with the guidelines of the Regulations for the Administration of Affairs Concerning Experimental Animals and was approved by the Biomedical Ethical Committee of Hebei Normal University, Hebei, China (ID: 2019SC14). No animals were harmed during the sampling process.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Huizhu Nan and Xin Lu contributed equally to this work.
==== Refs
References

1. Turlewicz-Podbielska H Ruszkowski JJ Wojciechowski J Pomorska-Mól M Seroprevalence of Neospora caninum in pet cats, dogs and rabbits from urban areas of Poland BMC Vet Res 2024 20 37 10.1186/s12917-024-03894-3 38297328
Turlewicz-Podbielska H, Ruszkowski JJ, Wojciechowski J, Pomorska-Mól M. Seroprevalence of Neospora caninum in pet cats, dogs and rabbits from urban areas of Poland. BMC Vet Res. 2024;20:37.38297328 10.1186/s12917-024-03894-3
2. Campero LM Basso W Moré G Fiorani F Hecker YP Echaide I Neosporosis in Argentina: Past, present and future perspectives Veterin Parasitol Reg Stud Rep 2023 41 100882
Campero LM, Basso W, Moré G, Fiorani F, Hecker YP, Echaide I, et al. Neosporosis in Argentina: Past, present and future perspectives. Veterin Parasitol Reg Stud Rep. 2023;41:100882.
3. Imhof D Pownall W Hänggeli KPA Monney C Román LR Ortega-Mora LM Immunization with a multivalent listeria monocytogenes vaccine leads to a strong reduction in vertical transmission and cerebral parasite burden in pregnant and non-pregnant mice infected with Neospora caninum Vaccines 2023 11 156 10.3390/vaccines11010156 36680001
Imhof D, Pownall W, Hänggeli KPA, Monney C, Román LR, Ortega-Mora LM, et al. Immunization with a multivalent listeria monocytogenes vaccine leads to a strong reduction in vertical transmission and cerebral parasite burden in pregnant and non-pregnant mice infected with Neospora caninum. Vaccines. 2023;11:156.36680001 10.3390/vaccines11010156
4. Marin MS Hecker YP Quintana S Pérez SE Leunda MR Cantón GJ Immunization with inactivated antigens of Neospora caninum induces toll-like receptors 3, 7, 8 and 9 in maternal-fetal interface of infected pregnant heifers Vet Parasitol 2017 243 12 17 10.1016/j.vetpar.2017.06.005 28807280
Marin MS, Hecker YP, Quintana S, Pérez SE, Leunda MR, Cantón GJ, et al. Immunization with inactivated antigens of Neospora caninum induces toll-like receptors 3, 7, 8 and 9 in maternal-fetal interface of infected pregnant heifers. Vet Parasitol. 2017;243:12–7.28807280 10.1016/j.vetpar.2017.06.005
5. Horcajo P Regidor-Cerrillo J Aguado-Martínez A Hemphill A Ortega-Mora LM Vaccines for bovine neosporosis: current status and key aspects for development Parasite Immunol 2016 38 709 723 10.1111/pim.12342 27314818
Horcajo P, Regidor-Cerrillo J, Aguado-Martínez A, Hemphill A, Ortega-Mora LM. Vaccines for bovine neosporosis: current status and key aspects for development. Parasite Immunol. 2016;38:709–23.27314818 10.1111/pim.12342
6. Balaji S Babu MM Iyer LM Aravind L Discovery of the principal specific transcription factors of Apicomplexa and their implication for the evolution of the AP2-integrase DNA binding domains Nucleic Acids Res 2005 33 3994 4006 10.1093/nar/gki709 16040597
Balaji S, Babu MM, Iyer LM, Aravind L. Discovery of the principal specific transcription factors of Apicomplexa and their implication for the evolution of the AP2-integrase DNA binding domains. Nucleic Acids Res. 2005;33:3994–4006.16040597 10.1093/nar/gki709
7. Campbell TL De Silva EK Olszewski KL Elemento O Llinás M Identification and genome-wide prediction of DNA binding specificities for the ApiAP2 family of regulators from the malaria parasite PLoS Pathog 2010 6 e1001165 10.1371/journal.ppat.1001165 21060817
Campbell TL, De Silva EK, Olszewski KL, Elemento O, Llinás M. Identification and genome-wide prediction of DNA binding specificities for the ApiAP2 family of regulators from the malaria parasite. PLoS Pathog. 2010;6:e1001165.21060817 10.1371/journal.ppat.1001165
8. Kaneko I Iwanaga S Kato T Kobayashi I Yuda M Genome-wide identification of the target genes of AP2-O, a Plasmodium AP2-family transcription factor PLoS Pathog 2015 11 e1004905 10.1371/journal.ppat.1004905 26018192
Kaneko I, Iwanaga S, Kato T, Kobayashi I, Yuda M. Genome-wide identification of the target genes of AP2-O, a Plasmodium AP2-family transcription factor. PLoS Pathog. 2015;11:e1004905.26018192 10.1371/journal.ppat.1004905
9. Santos JM Josling G Ross P Joshi P Orchard L Campbell T Red blood cell invasion by the malaria parasite is coordinated by the PfAP2-I transcription factor Cell Host Microbe 2017 21 731 741 10.1016/j.chom.2017.05.006 28618269
Santos JM, Josling G, Ross P, Joshi P, Orchard L, Campbell T, et al. Red blood cell invasion by the malaria parasite is coordinated by the PfAP2-I transcription factor. Cell Host Microbe. 2017;21:731–41.28618269 10.1016/j.chom.2017.05.006
10. Yuda M Kaneko I Iwanaga S Murata Y Kato T Female-specific gene regulation in malaria parasites by an AP2-family transcription factor Mol Microbiol 2020 113 40 51 10.1111/mmi.14334 31231888
Yuda M, Kaneko I, Iwanaga S, Murata Y, Kato T. Female-specific gene regulation in malaria parasites by an AP2-family transcription factor. Mol Microbiol. 2020;113:40–51.31231888 10.1111/mmi.14334
11. Shang X Shen S Tang J He X Zhao Y Wang C A cascade of transcriptional repression determines sexual commitment and development in Plasmodium falciparum Nucleic Acids Res 2021 49 9264 9279 10.1093/nar/gkab683 34365503
Shang X, Shen S, Tang J, He X, Zhao Y, Wang C, et al. A cascade of transcriptional repression determines sexual commitment and development in Plasmodium falciparum. Nucleic Acids Res. 2021;49:9264–79.34365503 10.1093/nar/gkab683
12. Yuda M Kaneko I Murata Y Iwanaga S Nishi T Mechanisms of triggering malaria gametocytogenesis by AP2-G Parasitol Int 2021 84 102403 10.1016/j.parint.2021.102403 34119684
Yuda M, Kaneko I, Murata Y, Iwanaga S, Nishi T. Mechanisms of triggering malaria gametocytogenesis by AP2-G. Parasitol Int. 2021;84:102403.34119684 10.1016/j.parint.2021.102403
13. Modrzynska K Pfander C Chappell L Yu L Suarez C Dundas K A Knockout screen of ApiAP2 genes reveals networks of interacting transcriptional regulators controlling the Plasmodium life cycle Cell Host Microbe 2017 21 11 22 10.1016/j.chom.2016.12.003 28081440
Modrzynska K, Pfander C, Chappell L, Yu L, Suarez C, Dundas K, et al. A Knockout screen of ApiAP2 genes reveals networks of interacting transcriptional regulators controlling the Plasmodium life cycle. Cell Host Microbe. 2017;21:11–22.28081440 10.1016/j.chom.2016.12.003
14. Shang X Wang C Fan Y Guo G Wang F Zhao Y Genome-wide landscape of ApiAP2 transcription factors reveals a heterochromatin-associated regulatory network during Plasmodium falciparum blood-stage development Nucleic Acids Res 2022 50 3413 3431 10.1093/nar/gkac176 35288749
Shang X, Wang C, Fan Y, Guo G, Wang F, Zhao Y, et al. Genome-wide landscape of ApiAP2 transcription factors reveals a heterochromatin-associated regulatory network during Plasmodium falciparum blood-stage development. Nucleic Acids Res. 2022;50:3413–31.35288749 10.1093/nar/gkac176
15. Russell TJ De Silva EK Crowley VM Shaw-Saliba K Dube N Josling G Inhibitors of ApiAP2 protein DNA binding exhibit multistage activity against Plasmodium parasites PLoS Pathog 2022 18 e1010887 10.1371/journal.ppat.1010887 36223427
Russell TJ, De Silva EK, Crowley VM, Shaw-Saliba K, Dube N, Josling G, et al. Inhibitors of ApiAP2 protein DNA binding exhibit multistage activity against Plasmodium parasites. PLoS Pathog. 2022;18:e1010887.36223427 10.1371/journal.ppat.1010887
16. Iyer LM Anantharaman V Wolf MY Aravind L Comparative genomics of transcription factors and chromatin proteins in parasitic protists and other eukaryotes Int J Parasitol 2008 38 1 31 10.1016/j.ijpara.2007.07.018 17949725
Iyer LM, Anantharaman V, Wolf MY, Aravind L, et al. Comparative genomics of transcription factors and chromatin proteins in parasitic protists and other eukaryotes. Int J Parasitol. 2008;38:1–31.17949725 10.1016/j.ijpara.2007.07.018
17. Srivastava S White MW Sullivan WJ Toxoplasma gondii AP2XII-2 contributes to proper progression through s-phase of the cell cycle mSphere 2020 5 e00520 e00542 10.1128/mSphere.00542-20 32611703
Srivastava S, White MW, Sullivan WJ. Toxoplasma gondii AP2XII-2 contributes to proper progression through s-phase of the cell cycle. mSphere. 2020;5:e00520–42.32611703 10.1128/mSphere.00542-20
18. Hong DP Radke JB White MW Opposing transcriptional mechanisms regulate Toxoplasma development mSphere. 2017 2 1659 1672 10.1128/mSphere.00347-16
Hong DP, Radke JB, White MW. Opposing transcriptional mechanisms regulate Toxoplasma development. mSphere. 2017;2:1659–72.10.1128/mSphere.00347-16
19. Radke JB Worth D Hong D Huang S Sullivan WJ Jr Wilson EH Transcriptional repression by ApiAP2 factors is central to chronic toxoplasmosis PLoS Pathog 2018 14 e1007035 10.1371/journal.ppat.1007035 29718996
Radke JB, Worth D, Hong D, Huang S, Sullivan WJ Jr, Wilson EH, et al. Transcriptional repression by ApiAP2 factors is central to chronic toxoplasmosis. PLoS Pathog. 2018;14:e1007035.29718996 10.1371/journal.ppat.1007035
20. Reid AJ Vermont SJ Cotton JA Harris D Hill-Cawthorne GA Könen-Waisman S Comparative genomics of the apicomplexan parasites Toxoplasma gondii and Neospora caninum: coccidia differing in host range and transmission strategy PLoS Pathog 2012 8 e1002567 10.1371/journal.ppat.1002567 22457617
Reid AJ, Vermont SJ, Cotton JA, Harris D, Hill-Cawthorne GA, Könen-Waisman S, et al. Comparative genomics of the apicomplexan parasites Toxoplasma gondii and Neospora caninum: coccidia differing in host range and transmission strategy. PLoS Pathog. 2012;8:e1002567.22457617 10.1371/journal.ppat.1002567
21. Lesage KM Huot L Mouveaux T Courjol F Saliou JM Gissot M Cooperative binding of ApiAP2 transcription factors is crucial for the expression of virulence genes in Toxoplasma gondii Nucleic Acids Res 2018 46 6057 6068 10.1093/nar/gky373 29788176
Lesage KM, Huot L, Mouveaux T, Courjol F, Saliou JM, Gissot M. Cooperative binding of ApiAP2 transcription factors is crucial for the expression of virulence genes in Toxoplasma gondii. Nucleic Acids Res. 2018;46:6057–68.29788176 10.1093/nar/gky373
22. Khelifa AS Guillen Sanchez C Lesage KM Huot L Mouveaux T Pericard P TgAP2IX-5 is a key transcriptional regulator of the asexual cell cycle division in Toxoplasma gondii Nat Commun 2021 12 116 10.1038/s41467-020-20216-x 33414462
Khelifa AS, Guillen Sanchez C, Lesage KM, Huot L, Mouveaux T, Pericard P, et al. TgAP2IX-5 is a key transcriptional regulator of the asexual cell cycle division in Toxoplasma gondii. Nat Commun. 2021;12:116.33414462 10.1038/s41467-020-20216-x
23. Zarringhalam K Ye S Lou J Rezvani Y Gubbels MJ Cell cycle-regulated ApiAP2s and parasite development: the Toxoplasma paradigm Curr Opin Microbiol 2023 76 102383 10.1016/j.mib.2023.102383 37898053
Zarringhalam K, Ye S, Lou J, Rezvani Y, Gubbels MJ. Cell cycle-regulated ApiAP2s and parasite development: the Toxoplasma paradigm. Curr Opin Microbiol. 2023;76:102383.37898053 10.1016/j.mib.2023.102383
24. Wang C Hu D Tang X Song X Wang S Zhang S Internal daughter formation of Toxoplasma gondii tachyzoites is coordinated by transcription factor TgAP2IX-5[J] Cell Microbiol 2021 23 e13291 10.1111/cmi.13291 33217129
Wang C, Hu D, Tang X, Song X, Wang S, Zhang S, et al. Internal daughter formation of Toxoplasma gondii tachyzoites is coordinated by transcription factor TgAP2IX-5[J]. Cell Microbiol. 2021;23:e13291.33217129 10.1111/cmi.13291
25. Srivastava S Holmes MJ White MW Sullivan WJ Jr Toxoplasma gondii AP2XII-2 contributes to transcriptional repression for sexual commitment mSphere 2023 8 e0060622 10.1128/msphere.00606-22 36786611
Srivastava S, Holmes MJ, White MW, Sullivan WJ Jr. Toxoplasma gondii AP2XII-2 contributes to transcriptional repression for sexual commitment. mSphere. 2023;8:e0060622.36786611 10.1128/msphere.00606-22
26. Chen L Tang X Sun P Hu D Zhang Y Wang C Comparative transcriptome profiling of Eimeria tenella in various developmental stages and functional analysis of an ApiAP2 transcription factor exclusively expressed during sporogony Parasit Vectors 2023 16 241 10.1186/s13071-023-05828-8 37468981
Chen L, Tang X, Sun P, Hu D, Zhang Y, Wang C, et al. Comparative transcriptome profiling of Eimeria tenella in various developmental stages and functional analysis of an ApiAP2 transcription factor exclusively expressed during sporogony. Parasit Vectors. 2023;16:241.37468981 10.1186/s13071-023-05828-8
27. Ma L Liu J Li M Fu Y Zhang X Liu Q Rhoptry protein 5 (ROP5) is a key virulence factor in Neospora caninum Front Microbiol 2017 8 370 10.3389/fmicb.2017.00370 28326073
Ma L, Liu J, Li M, Fu Y, Zhang X, Liu Q, et al. Rhoptry protein 5 (ROP5) is a key virulence factor in Neospora caninum. Front Microbiol. 2017;8:370.28326073 10.3389/fmicb.2017.00370
28. Wang F Wang X Song X Ma L Yang J Liu Q Function of Neospora caninum dense granule protein 7 in innate immunity in mice Parasitol Res 2021 120 197 207 10.1007/s00436-020-06961-4 33164154
Wang F, Wang X, Song X, Ma L, Yang J, Liu Q, et al. Function of Neospora caninum dense granule protein 7 in innate immunity in mice. Parasitol Res. 2021;120:197–207.33164154 10.1007/s00436-020-06961-4
29. Li M Wang H Liu J Hao P Ma L Liu Q The apoptotic role of metacaspase in Toxoplasma gondii Front Microbiol 2015 6 1560 26834715
Li M, Wang H, Liu J, Hao P, Ma L, Liu Q. The apoptotic role of metacaspase in Toxoplasma gondii. Front Microbiol. 2015;6:1560.26834715
30. Bartlett A O'Malley RC Huang SC Galli M Nery JR Gallavotti A Mapping genome-wide transcription factor binding sites using DAP-seq Nat Protoc 2017 12 1659 1672 10.1038/nprot.2017.055 28726847
Bartlett A, O’Malley RC, Huang SC, Galli M, Nery JR, Gallavotti A, et al. Mapping genome-wide transcription factor binding sites using DAP-seq. Nat Protoc. 2017;12:1659–72.28726847 10.1038/nprot.2017.055
31. Liu D Shi S Hao Z Xiong W Luo M OsbZIP81, a homologue of arabidopsis VIP1, may positively regulate JA levels by directly targetting the genes in JA signaling and metabolism pathway in rice Int J Mol Sci 2019 20 2360 10.3390/ijms20092360 31086007
Liu D, Shi S, Hao Z, Xiong W, Luo M. OsbZIP81, a homologue of arabidopsis VIP1, may positively regulate JA levels by directly targetting the genes in JA signaling and metabolism pathway in rice. Int J Mol Sci. 2019;20:2360.31086007 10.3390/ijms20092360
32. Hutin S Blanc-Mathieu R Rieu P Parcy F Lai X Zubieta C Identification of plant transcription factor DNA-binding sites using seq-DAP-seq Methods Mol Biol 2023 2698 119 145 10.1007/978-1-0716-3354-0_9 37682473
Hutin S, Blanc-Mathieu R, Rieu P, Parcy F, Lai X, Zubieta C. Identification of plant transcription factor DNA-binding sites using seq-DAP-seq. Methods Mol Biol. 2023;2698:119–45.37682473 10.1007/978-1-0716-3354-0_9
33. Ning L Wang J Shi P Wu C Dong Z Zhao H Mapping genome-wide binding sites of endosperm specific expression transcription factor O2 using DAP-Seq Chin Sci Bull 2019 64 2537 2548 10.1360/N972019-00334
Ning L, Wang J, Shi P, Wu C, Dong Z, Zhao H. Mapping genome-wide binding sites of endosperm specific expression transcription factor O2 using DAP-Seq. Chin Sci Bull. 2019;64:2537–48.10.1360/N972019-00334
34. He L Meng N Castellarin SD Wang Y Sun Q Li X Combined metabolite and transcriptome profiling reveals the norisoprenoid responses in grape berries to abscisic acid and synthetic auxin Int J Mol Sci 2021 22 1420 10.3390/ijms22031420 33572582
He L, Meng N, Castellarin SD, Wang Y, Sun Q, Li X, et al. Combined metabolite and transcriptome profiling reveals the norisoprenoid responses in grape berries to abscisic acid and synthetic auxin. Int J Mol Sci. 2021;22:1420.33572582 10.3390/ijms22031420
35. Zhu C Wang C Lu C Wang J Zhou Y Xiong M Genome-wide identification and expression analysis of OsbZIP09 target genes in rice reveal its mechanism of controlling seed germination Int J Mol Sci 2021 22 1661 10.3390/ijms22041661 33562219
Zhu C, Wang C, Lu C, Wang J, Zhou Y, Xiong M, et al. Genome-wide identification and expression analysis of OsbZIP09 target genes in rice reveal its mechanism of controlling seed germination. Int J Mol Sci. 2021;22:1661.33562219 10.3390/ijms22041661
36. D'Inca E Cazzaniga S Foresti C Vitulo N Bertini E Galli M VviNAC33 promotes organ de-greening and represses vegetative growth during the vegetative-to-mature phase transition in grapevine New Phytol 2021 231 726 746 10.1111/nph.17263 33567124
D’Inca E, Cazzaniga S, Foresti C, Vitulo N, Bertini E, Galli M, et al. VviNAC33 promotes organ de-greening and represses vegetative growth during the vegetative-to-mature phase transition in grapevine. New Phytol. 2021;231:726–46.33567124 10.1111/nph.17263
37. Trouillon J Ragno M Simon V Attrée I Elsen S Transcription inhibitors with XRE DNA-binding and cupin signal-sensing domains drive metabolic diversification in pseudomonas mSystems. 2021 6 e00753 e820 10.1128/mSystems.00753-20 33436508
Trouillon J, Ragno M, Simon V, Attrée I, Elsen S. Transcription inhibitors with XRE DNA-binding and cupin signal-sensing domains drive metabolic diversification in pseudomonas. mSystems. 2021;6:e00753-e820.33436508 10.1128/mSystems.00753-20
38. Cao Y Zeng H Ku L Ren Z Han Y Su H ZmIBH1-1 regulates plant architecture in maize J Exp Bot 2020 71 2943 2955 10.1093/jxb/eraa052 31990030
Cao Y, Zeng H, Ku L, Ren Z, Han Y, Su H, et al. ZmIBH1-1 regulates plant architecture in maize. J Exp Bot. 2020;71:2943–55.31990030 10.1093/jxb/eraa052
