Pavel Neumann,Ludmila Oliveira,Tae-Soo Jang,Petr Novák,Andrea Koblížková,Veit Schubert,Andreas Houben,Jiří Macas
Disruption of the standard kinetochore in holocentric Cuscuta species
kinetochore,centromere,holocentric,monocentric,Cuscuta
Significance Segregation of chromosomes depends on the centromere and the kinetochore, which mediates its connection to spindle microtubules. In contrast to monocentric chromosomes with a single centromere domain, holocentric chromosomes, which evolved from monocentric ones independently several times in both animals and plants, usually have hundreds of centromere units distributed along the entire chromosome length. Here, we report that holocentric chromosomes in the parasitic plants of the genus Cuscuta divide properly during mitosis despite dramatic changes in their kinetochore composition, suggesting a unique mechanism of chromosome segregation. Our results, together with recent discoveries in other species, suggest that centromere determination and kinetochore composition are remarkably plastic in holocentric organisms and thus may provide important insights into the principles of chromosome segregation.
Disruption of the standard kinetochore in holocentric Cuscuta species

Segregation of chromosomes depends on the centromere and the kinetochore, which mediates its connection to spindle microtubules. In contrast to monocentric chromosomes with a single centromere domain, holocentric chromosomes, which evolved from monocentric ones independently several times in both animals and plants, usually have hundreds of centromere units distributed along the entire chromosome length. Here, we report that holocentric chromosomes in the parasitic plants of the genus Cuscuta divide properly during mitosis despite dramatic changes in their kinetochore composition, suggesting a unique mechanism of chromosome segregation. Our results, together with recent discoveries in other species, suggest that centromere determination and kinetochore composition are remarkably plastic in holocentric organisms and thus may provide important insights into the principles of chromosome segregation.

Faithful segregation of chromosomes during mitosis and meiosis depends on the centromere, a chromosomal domain that facilitates the attachment of chromosomes to spindle microtubules. In monocentric chromosomes, the centromere is localized at a single site per chromosome, which is morphologically discernible as a primary constriction. Holocentric chromosomes, on the contrary, lack this primary constriction and instead have the centromere domains distributed along almost the entire chromosome length. Holocentricity evolved from monocentric organization independently several times during the evolution of both plants and animals (1); however, the causes of the transitions are still enigmatic. This is primarily because only a few holocentric species have been studied so far and because most groups of holocentric species evolved from the monocentric ancestors a long time ago, making the factors involved in the transition elusive. In most species, the centromere is epigenetically determined by the presence of CENH3, a centromere-specific variant of histone H3 that replaces the canonical H3 histones in centromeric nucleosomes (2). At the same time, CENH3 serves as the basis for the kinetochore, a complex multiprotein structure that mediates the connection between centromeric chromatin and the microtubules of the mitotic spindle in most species. The backbone of the kinetochore consists of the constitutive centromere-associated network (CCAN), which connects the kinetochore with centromeric chromatin, and the KMN network, which constitutes an interface toward spindle microtubules (3, 4). The function of the kinetochore is regulated by additional proteins, the most studied of which belong to the spindle assembly checkpoint (SAC) (5, 6) and the chromosome passenger complex (CPC) (7–9). The role of CENH3 in centromere determination predicts that the transition from monocentric to holocentric centromere organization requires the formation of CENH3-containing domains along entire chromosomes. Indeed, in the few holocentric species studied to date, CENH3 is typically localized along the entire poleward surface of each chromatid where microtubules attach (10, 11). An exception are holocentric insects that lack CENH3 and use an alternative pathway of kinetochore assembly that depends on CENP-T protein (12–14). Recently, we identified the first exception in plants, in Cuscuta europaea, which belongs to the holocentric subgenus Cuscuta of the parasitic plant genus Cuscuta (Convolvulaceae) (15). In this species, the chromosomes restrict CENH3 to only one to three heterochromatin bands, despite being attached to the mitotic spindle along their entire length. This suggests that CENH3 has either lost its centromere function in this species or acts in parallel with an additional CENH3-independent mechanism of kinetochore assembly. Since monocentric relatives of C. europaea from the sister subgenus, Grammica, and the more distant subgenus, Monogynella, have CENH3 localized specifically in primary constrictions (16), it is plausible that the peculiar CENH3 localization in C. europaea resulted from changes in kinetochore assembly that were linked to the transition to holocentricity in the subgenus Cuscuta. However, how kinetochore assembly has changed and whether these changes are related to the transition to holocentricity remain unknown. In this study, we addressed these questions by comparing the repertoire of major structural and regulatory kinetochore proteins and their chromosomal localization between two Cuscuta species from the holocentric subgenus Cuscuta (C. europaea and C. epithymum), two monocentric Cuscuta species from the sister subgenus Grammica (C. australis and C. campestris), and Ipomoea nil, which was included as an outgroup Convolvulaceae species. To obtain high-quality data for gene identification in the two holocentric Cuscuta species, we sequenced both their genomes and transcriptomes. The chromosomal localization of kinetochore proteins was determined using antibodies developed against key proteins representing different subcomponents of the kinetochore. Comparison of the results between monocentric and holocentric species allowed us to uncover an unprecedented level of changes that occurred specifically in the holocentric species and thus likely played an important role in the transition to holocentricity in Cuscuta.

Sequencing of the holocentric species C. europaea and C. epithymum resulted in genome assemblies of 975.8 Mb (N50 = 17.9 Mb) and 997 Mb (N50 = 3.3 Mb), respectively (SI Appendix, Supporting text and Table S1). The completeness of gene space and quality of gene prediction were assessed using BUSCO and were comparable to genome assemblies previously published for the monocentric Cuscuta relatives C. australis and C. campestris (SI Appendix, Fig. S1). The quality of gene prediction in the genome assembly was also verified by the independent assembly of the transcriptomes, which showed similar results following BUSCO analysis (SI Appendix, Table S2). To identify kinetochore protein sequences in the species selected for this study, we created a sequence database of 29 structural and regulatory kinetochore proteins known in plants. First, we used the database as a query for blastp searches to identify homologous protein sequences in the monocentric species C. australis, C. campestris, and Ipomoea nil. The identified sequences were manually verified and corrected when needed and added to the database to improve its sensitivity for homologous protein recognition. The improved database was then used for blastp searches in the two holocentric Cuscuta species. Comparison of the identified kinetochore protein genes revealed that all the 29 tested genes are present and mostly intact in the monocentric species, whereas in the holocentric species, some of the genes are either absent, significantly truncated, or duplicated accompanied by a higher rate of sequence divergence (Fig. 1A, SI Appendix, Table S3, and Dataset S1). Transcriptome data analysis revealed that the intact and truncated genes identified in holocentric Cuscuta species are transcribed. The lost genes included both eudicotyledonous plant homologs of KNL2, referred to as αKNL2 and βKNL2 (25), and four of the eight spindle assembly checkpoint (SAC) genes, namely, BMF1, BMF2, BMF3, and MAD2 (Fig. 1A). Their absence was in all cases confirmed by comparison of genomic loci possessing these genes in C. australis with the orthologous loci in C. europaea and C. epithymum (SI Appendix, Figs. S2–S4), as well as by their absence in genome-independent transcriptome assemblies. The only exception was BMF1 whose transcriptionally inactive fragment remains in C. epithymum (SI Appendix, Fig. S3). Large gene truncations took place in three structural kinetochore protein genes, including CENP-C, KNL1, and ZWINT1, and the SAC gene MAD1 (Figs. 2 and 3 and SI Appendix, Figs. S5 and S6). Finally, the CENH3 gene in holocentric species was found to have duplicated once in the common ancestor of C. europaea and C. epithymum, and once independently in each of the two species. The diversification of the duplicated CENH3 genes in holocentric species resulted in considerably higher protein sequence variability for CENH3 compared with monocentric Grammica species, suggesting that they evolved more rapidly (SI Appendix, Figs. S7–S10). Given the function of proteins that are either missing or truncated, the changes could have affected kinetochore assembly and function at multiple levels, from CENH3 loading (absence of KNL2) and kinetochore assembly (truncation of CENP-C, KNL1, and ZWINT1), to regulation of its function (absence of several key proteins of SAC) (Fig. 1 B and C).

Since KNL2 is essential for the proper loading of CENH3 to centromeres (25, 33–35), the loss of both αKNL2 and βKNL2 in holocentric Cuscuta species is likely to have a serious impact on CENH3 localization. On holocentric chromosomes, CENH3 is expected to specifically localize along the poleward side of each chromatid. In contrast to this expectation, we have previously shown that CENH3 occurs in all but one prominent transversal heterochromatin band in C. europaea and that CENH3 distribution does not correlate with the distribution of mitotic spindle attachment sites detected with antibodies against α-tubulin (15 and Fig. 4 A and B). To determine the localization of CENH3 in C. epithymum, we developed two antibodies against different N-terminal sequence variants of the proteins. Although the antibodies were made to recognize all CENH3 protein sequence variants present in the tested plant, neither of them produced a signal on chromosomes and nuclei that could be distinguished from the background (SI Appendix, Figs. S11, S12 A and B, and S13). On the contrary, both antibodies developed for C. epithymum detected CENH3 in the heterochromatin domains in C. europaea (SI Appendix, Fig. S11 C and E), demonstrating that they were functional for in situ detection. These results suggest that CENH3 is either not present in the chromatin of C. epithymum or that its levels are considerably lower than those in C. europaea, and thus below the limits of detection for the applied in situ immunodetection technique. This was also confirmed using western blot detection of CENH3 in the two holocentric species (SI Appendix, Fig. S13B). Despite the absence of CENH3 signals, α-tubulin immunostaining revealed attachment of mitotic spindle microtubules to chromosomes along their poleward sides, confirming the holocentric nature of chromosomes in C. epithymum (Fig. 4C). This was in contrast to monocentric Cuscuta spp., which had microtubules attached only to CENH3-containing domains (Fig. 4D). These results suggest that CENH3 does not function as a foundational kinetochore protein in holocentric Cuscuta species.

The chromosomal distribution of CENH3 together with the truncation of three structural kinetochore proteins suggested that kinetochore assembly may be impaired in holocentric Cuscuta species. To test whether the kinetochore assembles along the poleward chromosome surface, as expected for holocentric chromosomes, we examined the localization of CENP-C, which is a linker between CENH3 and the KMN network, and of MIS12, KNL1, and NDC80, which represent the three complexes of the KMN network (Fig. 1B). Antibodies were developed against peptides designed from domains that were conserved in the holocentric species. However, owing to the high sequence similarity between species, it was likely that the antibodies against KNL1, NDC80, and MIS12 would also recognize homologous proteins from monocentric Cuscuta species. Indeed, when these antibodies were used for in situ detection, monocentromeres in C. australis as well as in C. reflexa from the more distant subgenus Monogynella were labeled, demonstrating the functionality of the antibodies (Fig. 4 E–G and SI Appendix, Fig. S14). The antibodies against KNL1 and NDC80 proved to be particularly versatile, functioning even in Rhynchospora pubera, an evolutionarily very distant plant species with holocentric chromosomes, where they detected holocentromere-characteristic signals for both proteins (Fig. 4 H and I). In agreement with the lack of CENH3 signal in C. epithymum, CENP-C, KNL1, and NDC80 were not detected on either mitotic chromosomes or in interphase nuclei in this species (SI Appendix, Fig. S12). In C. europaea, these three proteins were detected in small subdomains embedded within CENH3-containing heterochromatin during interphase but not on mitotic chromosomes (Fig. 4 J–L, SI Appendix, Fig. S15, and Movie S1). Simultaneous in situ detection of KNL1 with either CENP-C or NDC80 revealed that these proteins fully colocalized (Fig. 4 M and N and Movies S2 and S3). These results suggest that the assembly of the kinetochore during interphase in C. europaea still depends, at least in part, on the presence of CENH3, but that kinetochore organization is disrupted before cells enter mitosis. Strikingly, 2 to 16 (n = 100) discrete spots were labeled with the MIS12 antibody in the interphase nuclei of both holocentric species (Fig. 4 O and P). In C. europaea, these spots were always located away from the CENH3-containing heterochromatin (Fig. 4O and Movie S4), suggesting the absence or depletion of MIS12 in domains that contain CENP-C, KNL1, and NDC80. No signals were detected with the MIS12 antibody on mitotic chromosomes in either holocentric Cuscuta species (SI Appendix, Fig. S16).

To test whether the regulatory kinetochore complexes form on chromosomes in holocentric Cuscuta species despite the absence of the tested kinetochore proteins and the massive loss of the SAC genes observed, we raised antibodies against BUB3;1/2 and Borealin, which are components of the SAC and CPC, respectively. While the BUB3;1/2 antibodies produced monocentric-like signals on chromosomes in C. reflexa, and holocentromere-like signals in R. pubera, BUB3;1/2 was not detectable on chromosomes in holocentric Cuscuta species (Fig. 5 A and B). On the contrary, the antibodies against Borealin labeled the chromosomes in the region around areas of sister chromatid cohesion at centromeres in monocentric C. reflexa and along the entire chromosome length in both holocentric Cuscuta species (Fig. 5 C–E). These results indicate that the conventional SAC is abolished, while the CPC maintains at least some of its functions in holocentric Cuscuta species.

The peculiar CENH3 localization in C. europaea described in our previous study (15) suggested that the transition to holocentricity in the genus Cuscuta may have been associated with the formation of a CENH3-independent kinetochore assembly. In this study, we have demonstrated that the transition to holocentricity in Cuscuta species was associated with extensive changes in structural and regulatory kinetochore protein genes, and disruption of both standard kinetochore assembly and SAC regulation of mitotic chromosome segregation. This distinguishes holocentric Cuscuta species from both the holocentric nematode Caenorhabditis elegans, which uses the CENH3–CENP-C pathway of kinetochore assembly (36), and holocentric insects, in which the CENH3–CENP-C pathway of kinetochore assembly was lost and replaced by the CENP-T pathway (12–14) (Fig. 6). We hypothesize that one of the most important changes in the evolution of holocentric Cuscuta species was the loss of KNL2. In C. elegans, RNAi depletion of KNL2 leads to a reduction in the presence of CENH3 to levels undetectable by immunodetection, resulting in chromosome segregation defects and embryonic lethality (33, 34). Similar phenotypes have been observed in KNL2 mutants in other species, including A. thaliana, demonstrating the general importance of KNL2 for CENH3 loading (25, 35). Therefore, the depletion/absence of CENH3 in C. epithymum chromatin could be due to the absence of both KNL2 variants. On the contrary, it is puzzling that CENH3 accumulates in heterochromatin domains in C. europaea despite the loss of KNL2. Given that all heterochromatin domains that contain CENH3 possess the same repetitive sequences, whereas the heterochromatin domain that lacks these repeats also lacks CENH3 (15, 37), the incorporation of CENH3 into these domains could be DNA sequence dependent. In light of the importance of KNL2 and CENH3 for centromere determination and kinetochore assembly, it is surprising that the loss of KNL2 in both holocentric Cuscuta species, the depletion/absence of CENH3 on chromosomes in C. epithymum, and the peculiar CENH3 distribution on chromosomes in C. europaea are neither lethal nor cause chromosome segregation defects. The simplest explanation is that CENH3 is no longer necessary for correct chromosome segregation in holocentric Cuscuta species (SI Appendix, Fig. S17). The absence of detectable levels of structural kinetochore proteins on mitotic chromosomes in holocentric Cuscuta species is in contrast not only to monocentric Cuscuta species but also to the holocentric-like distribution of NDC80 and KNL1 in R. pubera (Cyperaceae), which was used as a holocentric control plant in this study (Fig. 4). This suggests that the formation of the standard kinetochore is disrupted in holocentric Cuscuta species. In C. epithymum, this could be primarily a direct consequence of the depletion/absence of CENH3 on the chromosomes. In C. europaea, the causes of kinetochore disruption must be different because CENH3-containing heterochromatin is present throughout the cell cycle and partially colocalizes with CENP-C, KNL1, and NDC80 proteins during interphase. The reasons why the putative complex of kinetochore proteins formed during interphase disappears at the onset of mitosis are not clear. Considering that three structural kinetochore proteins are truncated (Fig. 1A), one possibility is that the complex falls apart because of disrupted interactions between kinetochore components (Fig. 1C). We hypothesize that the N-terminal truncation of CENP-C may play an important role because CENP-C is the only protein known to link centromeric chromatin to the outer kinetochore in plants (Figs. 1 B and C and 2A). Although the N terminus of CENP-C is divergent in sequence between eukaryotes, it has been shown to bind MIS12c in both humans and yeast, indicating a conserved function (18, 38, 39). Given that this function is also conserved in plants, the N-terminal truncation of CENP-C in C. europaea should interfere with MIS12c binding. Consistent with this notion, the MIS12 antibody did not detect signals in CENH3-containing domains, suggesting that MIS12 is absent or depleted in domains containing CENP-C, KNL1, and NDC80 (Fig. 4O and Movie S4). Because we cannot rule out the possibility that the MIS12 antibody cross-reacts with other proteins, further evidence is needed to confirm that the signals in discrete CENH3-lacking domains are indeed due to MIS12 accumulation. While the colocalization of CENP-C, KNL1, and NDC80 suggests that the kinetochore assembles during interphase, despite the absence of MIS12, the complex may not be sufficiently stable to survive mitosis. The N-terminal truncation of CENP-C is, however, unlikely to cause the disappearance of the protein itself because the N terminus is not required for the binding of centromeric nucleosomes (Fig. 2A and 40). Although the internal portion of CENP-C contains a domain that binds centromeric nucleosomes in humans and yeast (Fig. 2A), the high sequence divergence of CENP-C prevented us from determining by a sequence similarity–based approach whether it overlaps with the region lost in C. europaea. On the contrary, the large size disparity between the domains containing CENH3 and CENP-C (Fig. 4 J and M and Movies S1 and S2) suggests that there is an imbalance between the levels of the two proteins that may reflect inefficient binding of CENP-C to CENH3. The results discussed above support a model in which holocentric Cuscuta species either use substantially reduced kinetochores lacking CENH3, CENP-C, KNL1, MIS12, and NDC80 or, more likely, have evolved a completely unique mechanism of chromosome attachment to the mitotic spindle. This conclusion is also supported by the degeneracy of SAC genes that would have been required had the kinetochore been present and functioning in a conventional manner. The absence of detectable amounts of KNL1, NDC80, and MIS12 on metaphase chromosomes distinguishes holocentric Cuscuta species from holocentric insects in which kinetochores lack CENH3 and CENP-C but have a complete KMN network connected to centromeric chromatin via CENP-T (12, 14). Although the existence of CENP-T in plants cannot be ruled out, it is unlikely that the potential presence of CENP-T in holocentric Cuscuta species could compensate for the absence/depletion of CENP-C and KMN proteins on mitotic chromosomes, unless the hypothetical CENP-T protein evolved to connect centromeric chromatin and microtubules either directly or via interaction with other proteins. Alternative kinetochores have already been described in Kinetoplastida, most of which have lost CENH3 and all CCAN and KMN genes. They consist of proteins that probably evolved from meiotic components of chromosome synapsis and homologous recombination machinery (41, 42). Moreover, kinetochore-independent chromosomal movement along the spindle, facilitated by kinesin motor proteins, has been described for acentric chromosomes in Drosophila neuroblasts (43, 44) and for chromatin knobs in maize (45, 46). Given that all holocentric Cuscuta species are limited to the subgenus Cuscuta, while monocentric species belong to either the sister subgenus Grammica or the evolutionarily more distant subgenus Monogynella, holocentricity in the genus Cuscuta has probably a monophyletic origin. Although C. europaea and C. epithymum diverged only 1.2 to 3.2 Mya, the subgenus Cuscuta is much older, having evolved 16.6 to 17.5 Mya (SI Appendix, Fig. S10 and 16). Therefore, studying only two closely related holocentric Cuscuta species can only provide limited insight into the evolution of holocentricity. In this regard, it will be particularly intriguing to identify changes shared by all holocentric Cuscuta species, as these may have played a crucial role in the transition to holocentricity. Overall, we have shown that the transition to holocentricity in Cuscuta species was unique among all species studied to date. It was accompanied, and perhaps even triggered, by the degeneration of standard kinetochore structure and regulation and the formation of a unique mechanism for chromosome attachment to microtubules. The insights gained in this study provide the basis for future studies aimed at uncovering the plasticity of kinetochore assembly and discovering yet unknown principles of chromosome segregation.

Seeds of C. europaea (serial number: 0101147) were obtained from the Royal Botanic Gardens (Ardingly, UK). C. epithymum plants were collected from a natural population at “U Cáby” (Kroclov, Czech Republic). Seeds of C. australis and C. campestris were provided by Prof. Jianqiang Wu (Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China) and Dr. Chnar Fathoulla (University of Salahaddin, Kurdistan Region, Iraq), respectively. C. reflexa Roxb. plant was obtained from the Botanic Gardens of the Rhenish Friedrich-Wilhelm University (Bonn, Germany). Cuscuta plants were cultivated on the following host plant species: Urtica dioica (C. europaea), Betonica officinalis and Coleus blumei (C. epithymum), Ocimum basilicum (C. australis and C. campestris), or Pelargonium zonale (C. reflexa). Plants of R. pubera were obtained from Dr. André Marques (Max Planck Institute for Plant Breeding Research, Cologne, Germany).

DNA for Illumina and Pac-Bio sequencing was isolated using the CTAB method from nuclei extracted from young shoots of C. europaea and C. epithymum as described previously (47). Shotgun Illumina paired-end sequencing of DNA was performed by the Brigham Young University (Provo, UT, USA) and Admera Health (South Plainfield, NJ, USA). High-molecular-weight nuclear DNA used for Oxford nanopore sequencing was isolated using a modified CTAB protocol as described previously (48). Nanopore sequencing was performed as described (37). Detailed information about all genome sequence datasets produced in this study is provided in SI Appendix, Table S4. Illumina paired-end reads and Oxford nanopore reads were assembled using MaSuRCA (49). PacBio HiFi reads were assembled using Hifiasm assembler [v0.15.5-r350; (50)] with default parameters for PacBio HiFi sequence reads. Since the quality of the HiFi-based assemblies was considerably better than those generated by MaSuRCA (SI Appendix, Table S1), they were selected for submission to European Nucleotide Archive (https://www.ebi.ac.uk/ena/browser/home; accession numbers: GCA_945859875 (C. europaea) and GCA_945859915 (C. epithymum)). Completeness and contiguity of assemblies were evaluated using BUSCO [v5.2.2; (51)] and QUAST [v5.0.2; (52)]. Genome characteristics were evaluated using kmer analysis and the jellyfish program (53) with kmer lengths 21 and 51 for Illumina and PacBio HIFI sequence reads, respectively. Heterozygosity was estimated using GenomeScope program (54).

Total RNA was isolated using the Trizol method. Preliminary sequencing for de-novo transcriptome assemblies of C. epithymum, C. europaea, and C. campestris was performed at GATC Biotech (Konstanz, Germany) using Illumina technology producing 50 bp paired-end reads. In each species, RNA was isolated from shoots and inflorescences, mixed in a 1:1 ratio, treated with DNase I (Ambion, Austin, TX, USA), and then enriched for poly-A fraction using the Dynabeads mRNA purification kit (Thermo Fisher Scientific, Waltham, MA, USA). Deep transcriptome sequencing of C. epithymum, C. europaea, and C. australis was done using RNA isolated from shoot tips, shoot internodia, or inflorescences at various stages of development. For each species and tissue, the RNA samples were produced in three biological replicates (samples from different plants collected at different times). Subtraction of poly-A RNA using NEBNext Ultra II with a Poly-A Selection kit (New England Biolabs, Ipswich, MA, USA) and poly-A RNA sequencing were performed at Admera Health (South Plainfield, NJ, USA). The sequencing generated more than 500 million 151 nt-long paired-end reads for each RNA sample, giving a total yield of about 5 billion reads per species (SI Appendix, Table S5). Transcriptomes were de-novo assembled using the Trinity program (55) with default options from paired-end reads. Sequences from individual replicates and tissue samples of each species were concatenated before their assembly. The presence of single-copy orthologs in the transcriptomes was evaluated using the BUSCO (v.5.2.2) program (51). To create gene models, paired-end RNA-Seq Illumina reads were aligned to genome assembly using the STAR program [v2.7.7a; (56)] with parameters --outSAMstrandField intronMotif --outSAMtype BAM SortedByCoordinate --alignIntronMax 20000. Each sample was aligned independently. Resulting alignments were merged into a single BAM file using samtools (57). Whole-length transcripts and genes were then reconstructed using the StringTie program [v2.1.7; (58)] with parameters -c 2 -f 0.05. Candidate coding regions within transcript sequences were identified using TransDecoder program (https://github.com/TransDecoder/TransDecoder) with default settings. Predicted protein sequences from C. europaea and C. epithymum were compared with published proteomes of C. campestris, C. australis, and Ipomoea nil using program OrthoFinder [v2.5.2; (59)] to identify orthologs and orthogroups. Genome assemblies and associated files containing detailed information about predicted gene models, protein, and CDS sequences were downloaded from http://plabipd.de/portal/cuscuta-campestris (C. campestris) or GenBank (https://www.ncbi.nlm.nih.gov/genbank/; C. australis: GCA_003260385.1; I. nil: GCF_001879475.1). RNA-Seq data for these species were downloaded from the Sequence Read Archive (SRA; https://www.ncbi.nlm.nih.gov/sra) from the following accession numbers: SRR6664647–SRR6664654 (C. australis), ERR1916345–ERR1916364 (C. campestris), and DRR024544–DRR024549 (I. nil). The RNA-Seq data produced in this study or downloaded from other studies were used to verify and correct automatically predicted gene models if needed. Manual verification and editing of gene models were performed using Apollo Genome Annotation Editor (60).

Structural and regulatory kinetochore protein sequences identified in A. thaliana were downloaded from uniprot database and published studies (27, 29, 61). These sequences were used for blastp searches to identify their homologs in genome assemblies of C. australis and C. campestris (62, 63), representing monocentric Cuscuta species, and in I. nil (64), selected as a monocentric nonparasitic genus of the family Convolvulaceae. All sequences with significant similarity hits were manually inspected to remove false positives, correct erroneous protein sequences, or add additional variants due to alternative splicing. Protein sequences from A. thaliana and the three Convolvulaceae species were combined into a reference dataset that was used for blastp and tblastn searches to find homologous kinetochore protein genes in holocentric C. epithymum and C. europaea. The searches were primarily performed in gene and protein sequences predicted using StringTie in the assembly produced from Pac-Bio reads, but the results were verified using the data from the parallel genome assemblies that were made from Illumina and nanopore reads as well as the transcriptome assemblies produced using Trinity. CENH3 sequences from additional Cuscuta species or other plants of the same Cuscuta species were obtained from our previous study [C. campestris, C. japonica (15)], identified in transcriptome shotgun assemblies (C. reflexa, C. campestris) or other available genome assembly (C. epithymum), amplified from RNA using RT-PCR or RACE methods (C. epithymum), or reconstructed from available next-generation genome sequence data using GRABb and GeneWise programs [C. americana, C. californica, C. pentagona; (65, 66)]. More detailed information about the sources of the CENH3 sequences is provided in SI Appendix, Table S6. Sequence alignments were performed using MUSCLE (67). Time trees were inferred using ITS and rbcL sequences and methods described in our previous study (16). ITS and rbcL sequences from C. australis were reconstructed from Illumina paired-end reads (SRA run accession number: SRR5851367) using RepeatExplorer (68). A search for conserved sequence motives was performed using MEME (69). Sequence logos were generated using WebLogo (70). The sources of CENP-C and ZWINT1 sequences used for MEME and WebLogo analyses are provided in SI Appendix, Table S7.

Antibodies to all kinetochore proteins used in this study were custom-produced by GenScript (Piscataway, NJ, USA) or Biomatik (Cambridge, ON, Canada) against peptides designed from regions that were most conserved among Cuscuta species and I. nil (SI Appendix, Fig. S18 and Table S8). Detailed information about the development and validation of the antibodies is provided in SI Appendix, Supporting text. The mouse monoclonal antibody to α-tubulin was purchased from Sigma-Aldrich (St. Louis, MO, USA; catalog number: T6199). Western blot and in situ immunodetection techniques are described in SI Appendix, Supporting text.

For conventional wide-field fluorescence microscopy, a Zeiss AxioImager.Z2 microscope equipped with an Axiocam 506 mono camera was used along with an Apotome2.0 device for better resolution in the z-axis, which was needed when the images were composed of multiple optical sections. Images were generated using the ZEN 3.2 software (Carl Zeiss GmbH). To capture signals at the super-resolution level (~120 nm using a 488-nm laser), spatial structured illumination microscopy (3D-SIM) was performed using a 63×/1.4 Oil Plan-Apochromat objective on an Elyra PS.1 microscope system, controlled by the ZENBlack software (Carl Zeiss GmbH). Images were captured using the 405-, 488-, and 561-nm laser lines for excitation and the appropriate emission filters (71). Three-dimensional movies were produced from 3D-SIM image stacks using the Imaris 9.7 (Bitplane) software.

Click here for additional data file. Click here for additional data file.