Skip to main content

Constant ratio between the genomic components of bipartite begomoviruses during infection and transmission

Abstract

The genomic components of multipartite viruses are encapsidated in separate virus particles, and the frequencies of genomic components represent one of the key genetic features. Many begomoviruses of economic significance are bipartite, and the details of the association between their genomic components remain largely unexplored. We first analyzed the temporal dynamics of the quantities of DNA-A and DNA-B and the B/A ratio of the squash leaf curl China virus (SLCCNV) in plants and found that while the quantities of DNA-A and DNA-B varied significantly during infection, the B/A ratio remained constant. We then found that changes in the B/A ratio in agrobacteria inoculum may significantly alter the B/A ratio in plants at 6 days post inoculation, but the differences disappeared shortly thereafter. We next showed that while the quantities of DNA-A and DNA-B among plants infected by agrobacteria, sap transmission and whitefly-mediated transmission differed significantly, the B/A ratios were similar. Further analysis of gene expression revealed that the ratio of the expression of genes encoded by DNA-A and DNA-B varied significantly during infection. Finally, we monitored the temporal dynamics of the quantities of DNA-A and DNA-B and the B/A ratio of another bipartite begomovirus, and a constant B/A ratio was similarly observed. Our findings highlight the maintenance of a constant ratio between the two genomic components of bipartite begomoviruses during infection and transmission, and provide new insights into the biology of begomoviruses.

Introduction

Viruses pose considerable threats to human health and crop production worldwide [1, 2]. As a group of obligate parasites, most viruses are characterized by the encapsidation of viral genomic components by protein coats. Based on the architecture, organization and packaging of viral genomic components, viruses can be divided into three categories, namely monopartite, segmented and multipartite [3, 4]. Monopartite viruses have only one genomic component that is encapsidated in a protein coat. The genome of segmented viruses consists of two or more genomic components that are encapsidated together in a single virus particle. Multipartite viruses have several genomic components that are packaged into separate virus particles [3, 4]. Multipartitism is prevalent, especially in plant viruses, accounting for 30–40% of plant virus genera and families [3]. Studies in recent years have highlighted that the relative frequencies of the genomic components of multipartite viruses may vary significantly in virus life cycle [3, 5,6,7]. And more importantly, variation of the relative frequencies of viral genomic components may significantly impact the pathogenesis of viruses in their hosts [3, 4, 8,9,10].

Begomoviruses (family Geminiviridae) are a group of plant single-stranded DNA (ssDNA) viruses that are transmitted by whiteflies of the Bemisia tabaci complex [11, 12]. The genus Begomovirus is now the largest in the entire virosphere as recognized by the International Committee on Taxonomy of Viruses [11]. Moreover, in recent decades begomoviruses have exacted a heavy toll on the production of many solanaceous, cucurbitaceous, malvaceous and leguminous crops in warm and temperate regions around the globe [12]. Diseases caused by begomoviruses that are of particular significance include tomato yellow leaf curl, cotton leaf curl and cassava mosaic [12]. Begomoviruses can be either monopartite or bipartite depending on the number of genomic components [11]. The genome of bipartite begomoviruses consists of two circular ssDNA molecules of around 2.6 kb, referred to as DNA-A and DNA-B. Monopartite begomoviruses contain only one genomic component that resembles the DNA-A of bipartite begomoviruses [11]. Out of the 445 species proposed in the genus Begomovirus, 164 are bipartite, and the majority of them are distributed in the New World [13]. In recent years, an increase is witnessed in the incidence of crop diseases caused by bipartite begomoviruses in Asia such as squash leaf curl China virus (SLCCNV) and Sri Lankan cassava mosaic virus (SLCMV) [14,15,16,17]. Research efforts to disentangle the molecular biology of these viruses may help to combat them and thereby safeguard crop production.

Although the evolutionary origin of begomoviral DNA-A has not been unambiguously resolved, distinct origins of DNA-A and DNA-B have been proposed [18,19,20]. DNA-B was proposed to be originated as a satellite molecule that was later captured by monopartite begomoviruses or their progenitors [19]. Sequence analyses revealed that while always associated, DNA-A and DNA-B underwent distinct evolutionary histories, with DNA-B harboring more genetic variations resulted from mutation and recombination than DNA-A [19, 21]. DNA-A encodes proteins required for virus replication, gene expression regulation, encapsidation, whitefly transmission and host defense suppression, and DNA-B encodes two proteins that regulate intra- and intercellular virus movement in plants [11, 22]. Both DNA-A and DNA-B are required for virus infection in plants, with only a few exceptions [11, 23]. As for the movement of the two DNA components, it was found that DNA-A and DNA-B of abutilon mosaic virus, a bipartite begomovirus, spread within N. benthamiana plants independently during the early course of systemic infection, resulting in a stochastic distribution of DNA-A- and DNA-A/B-infected nuclei [24]. Nevertheless, while the importance of the association between DNA-A and DNA-B has been well established, the quantitative relationship between the two genomic components remain barely characterized. For example, how is the ratio between the two genomic components regulated during virus infection, and is this ratio affected by the process of virus transmission?

In this study, we examined the ratio between the two DNA components of bipartite begomoviruses during infection and transmission. First, we monitored the temporal dynamics of the quantities of SLCCNV DNA-A and DNA-B and their ratio in plants. Second, we examined the effects of B/A ratio in agrobacteria inoculum on the dynamics of the quantities of SLCCNV DNA-A and DNA-B and their ratio in plants. Third, we analyzed the impact of sap- and whitefly-mediated transmission on the quantities of SLCCNV DNA-A and DNA-B and their ratio in plants. Fourth, we analyzed the dynamics of the expression of genes encoded by SLCCNV DNA-A and DNA-B and their ratio in plants. Finally, we monitored the temporal dynamics of the quantities of DNA-A and DNA-B and the B/A ratio of another bipartite begomovirus SLCMV. Our results provide new insights into the biology of bipartite begomoviruses.

Materials and methods

Plants, viruses and whiteflies

Four species of plants were used, namely cotton (Gossypium hirsutum cv. Zhe-Mian 1793), zucchini (Cucurbita pepo cv. Faguodongkui), squash (C. moschata cv. Mibennangua), and tobacco (Nicotiana benthamiana). All plants were grown in insect-proof greenhouses under controlled temperature at 25 ± 3 ℃. Cotton plants were used for whitefly rearing, zucchini and squash plants were used for the experiments with SLCCNV, and tobacco plants were used for the experiments with SLCMV. For agrobacteria-mediated virus inoculation and virus transmission experiment, zucchini and squash plants at 1–2 true leaf stage were used, and tobacco plants at five true leaf stage were used. Two bipartite begomoviruses were used, namely SLCCNV (GenBank accession codes: MG525551 [DNA-A] and MG525552 [DNA-B]) and SLCMV (GenBank accession codes: OK571385 [DNA-A] and OK571386 [DNA-B]). Both SLCCNV isolate Guangxi2017 and SLCMV isolate Cambodia2015 were collected from field, and infectious clones were constructed as mentioned before [25, 26]. No viral satellite has been reported to be associated with the two bipartite begomoviruses in the field [14,15,16,17]. Agrobacteria containing infectious clones of DNA-A or DNA-B were first cultured separately until OD600 reached 1.5, and then pelleted and resuspended in resuspension buffer (10 mM MgCl2, 10 mM MES, 200 µM acetosyringone). Next, agrobacteria were mixed and introduced into all the true leaves and cotyledons of zucchini or squash plants, and the first true leaf at the bottom for tobacco plants. Agroinoculation was conducted with 1 ml syringe. For insects, a culture of Asia II 1 whiteflies of the B. tabaci complex (mtCOI GenBank accession codes: DQ309077) was used. Whiteflies were originally collected from field and thereafter reared on cotton plants in insect-proof cages in climate chambers (26 ± 2 ℃, 60–80% relative humidity and 14/10 h light/dark cycles).

Extraction of DNA and quantification of DNA-A and DNA-B in plants

DNA extraction was performed using Plant Genomic DNA Kit (Tiangen, China). Ensuring quantitative PCR(qPCR) analysis of viral DNA-A and DNA-B was performed using SYBR Green Premix Pro Taq HS qPCR Kit (Accurate Biology, China) and CFX96 Real-Time PCR Detection System (Bio-Rad, USA). The actin gene was used as an internal reference. Each gene was analyzed in duplicate technical repeats for each of the samples. The average threshold cycle (Ct) was calculated per sample. Primers are listed in Table 1.

Table 1 Primers used in this study

Extraction of RNA and analysis of gene expression level

Total RNAs were extracted with TRIzol following user manual. cDNA was synthesized using Evo M-MLV RT Kit with gDNA Clean for qPCR (Accurate Biology, China). qPCR was performed as mentioned above with primers listed in Table 1.

Analysis of temporal dynamics of the quantities of DNA-A and DNA-B and the expression of AV1 and BV1 in plants

To obtain virus-infected plants, equal quantity of agrobacteria containing infectious clones of DNA-A and DNA-B were mixed and then used for agroinoculation. Next, the first fully-expanded leaves of inoculated plants were collected at designated days post inoculation (dpi) and subjected to DNA extraction and the quantification of DNA-A and DNA-B. For SLCCNV-infected zucchini plants, at 12, 18, 24 and 30 dpi, RNA samples were similarly collected from the first fully-expanded leaves and subjected to RNA extraction and the quantification of AV1, BV1, and actin mRNA level. All plants were sampled only once.

Analysis of the effects of B/A ratio in agrobacteria inoculum on quantities of DNA-A and DNA-B and their ratio in plants

To prepare agrobacteria inoculum with different B/A ratio (4, 1 and 1/4), agrobacteria containing infectious clones of DNA-A were mixed with different quantity of agrobacteria containing infectious clones of DNA-B. In these mixed solutions, the final OD600 value of agrobacteria containing infectious clones of DNA-A was kept constant. These agrobacteria solutions were then used for inoculation in zucchini plants. At 6, 18 and 30 dpi, the first fully-expanded leaves of the inoculated plants were collected and subjected to DNA extraction and DNA-A and DNA-B quantification.

Analysis of the dynamics of quantities of DNA-A and DNA-B and their ratio in plants during sap- and whitefly-mediated virus transmission

SLCCNV-infected zucchini plants were obtained by agroinoculation using a mixture of equal quantity of agrobacteria containing infectious clones of DNA-A and DNA-B. SLCCNV-infected plants at 18 dpi were sampled and the remaining plants were used as the inoculum plants in sap- and whitefly-mediated virus transmission. For sap-mediated transmission, symptomatic leaves were collected, grounded in liquid nitrogen and dissolved in 0.01 M phosphate buffered saline (pH 7.4) (3:20, wt/vol). The resultant mixtures were filtered using medical gauze and then inoculated onto zucchini leaves by rubbing with carborundum powder (600 mesh). For whitefly-mediated transmission, whiteflies were collected and transferred to infected plants for a 4 days virus acquisition. Viruliferous whiteflies were then collected and released onto zucchini seedlings for a 4 days virus transmission using leaf-clip cages as mentioned before [26]. The number of whiteflies per test plants was ten. Post transmission, plants were all kept in climate chambers for 18 days, and then subjected to DNA extraction and DNA-A and DNA-B quantification.

Statistical analysis

QPCR data of DNA-A and DNA-B quantity, and AV1 and BV1 expression were normalized to that of plant actin using 2−∆Ct method. Next, B/A ratio or BV1/AV1 were calculated for each sample. One-way analysis of variance (ANOVA) along with Fisher’s least significant difference (LSD) was used for the analysis of statistical significance. All data were presented as the mean ± standard errors of mean (mean ± SEM). Statistical analyses were conducted using SPSS Statistics 20.0 and EXCEL.

Results

Temporal dynamics of the quantities of SLCCNV DNA-A and DNA-B and their ratio in plants

We analyzed the temporal dynamics of the quantities of SLCCNV DNA-A and DNA-B and their ratio in zucchini and squash, two susceptible cucurbitaceous host plants of SLCCNV [26]. When zucchini plants were inoculated with a mixture of equal quantities of agrobacteria containing infectious clones of DNA-A and DNA-B, the quantities of DNA-A and DNA-B increased from 6 to 18 dpi and then changed only marginally from 18 to 30 dpi (Fig. 1A and B). When squash plants were inoculated, the relative quantities of DNA-A and DNA-B increased from 6 to 12 dpi, decreased from 12 to 24 dpi, and did not show significant changes from 24 to 30 dpi (Fig. 2A and B). Interestingly, in both zucchini and squash, the B/A ratios showed no significant changes from 6 to 30 dpi (Figs. 1 and 2 C, Table 2).

Fig. 1
figure 1

Temporal dynamics of the quantities of SLCCNV DNA-A and DNA-B and their ratio in zucchini plants. Zucchini plants were inoculated with a mixture of equal quantities of agrobacteria containing infectious clones of SLCCNV DNA-A and DNA-B, and then sampled at designated dpi for the quantification of DNA-A (A), DNA-B (B) and B/A ratio (C). Values are means ± SEM. The number of plants tested was 10 at 6 dpi, 6 at 12 dpi, 9 at 18 dpi, 8 at 24 and 30 dpi. Different letters above the columns indicate significant differences (one-way ANOVA, P < 0.05)

Fig. 2
figure 2

Temporal dynamics of the quantities of SLCCNV DNA-A and DNA-B and their ratio in squash plants. Squash plants were inoculated with a mixture of equal quantities of agrobacteria containing infectious clones of SLCCNV DNA-A and DNA-B; relative quantities of DNA-A and DNA-B and the B/A ratio (A, B and C). Values are means ± SEM. The number of plants tested was 10 at 6 dpi, 12 at 12, 18, 24 and 30 dpi. Different letters above the columns indicate significant differences (one-way ANOVA, P < 0.05)

Table 2 Statistics of one-way ANOVA

Effects of B/A ratio in agrobacteria inoculum on quantities of SLCCNV DNA-A and DNA-B and their ratio in plants

We prepared three agrobacterium inocula containing different ratios of B/A (4, 1 and 1/4), and kept the level of DNA-A constant. For both DNA-A and DNA-B, B/A ratio in agrobacterium inocula did not significantly affected their quantities at 6, 18 and 30 dpi, respectively (Fig. 3A, B; Table 2). A significantly higher B/A ratio was found in plants inoculated with the inoculum with the highest ratio of B/A at 6 dpi. However, the difference disappeared at 18 and 30 dpi (Fig. 3C; Table 2).

Fig. 3
figure 3

Effects of the B/A ratio in agrobacteria inoculum on the temporal dynamics of quantities of SLCCNV DNA-A and DNA-B and their ratio in zucchini plants. Agrobacteria containing infectious clones of SLCCNV DNA-A were mixed with different quantity of agrobacteria containing infectious clones of DNA-B. The final OD600 value of agrobacteria containing infectious clones of DNA-A was kept constant in these solutions; relative quantities of DNA-A and DNA-B, and the B/A ratio (A, B and C). Values are means ± SEM. The number of plants tested for the three treatments was 8, 12 and 10 at 6 dpi, 12, 11 and 11 at 18 dpi, 12, 12 and 11 at 30 dpi. *** indicates significant differences (one-way ANOVA, P < 0.001)

Dynamics of the quantities of SLCCNV DNA-A and DNA-B and their ratio in plants during transmission

We compared the quantities of DNA-A and DNA-B among plants of three treatments of virus inoculation: agrobacteria-inoculation, sap-mediated transmission and whitefly-mediated transmission. The methods of virus inoculation significantly affected the quantities of DNA-A and DNA-B, but not the B/A ratio (Fig. 4; Table 2). Specifically, the highest quantity of DNA-A and DNA-B was found in agrobacteria-inoculated plants, followed by the whitefly-inoculated plants, and then the sap-inoculated plants (Fig. 4A and B).

Fig. 4
figure 4

Quantities of SLCCNV DNA-A and DNA-B and their ratio in zucchini plants that were inoculated by agrobacteria, sap and whitefly. Zucchini plants were inoculated with a mixture of equal quantities of agrobacteria containing infectious clones of SLCCNV DNA-A and DNA-B. At 18 dpi, some plants were sampled and the others were used as the source of inoculum for sap- and whitefly-mediated virus transmission. At 18 days post transmission, sap- and whitefly-inoculated plants were sampled; relative quantities of DNA-A and DNA-B and the B/A ratio (A, B and C). Values are means ± SEM. The number of plants tested was 15 for agrobacteria and sap treatments, 10 for whitefly treatment. Different letters above the columns indicate significant differences (one-way ANOVA, P < 0.05)

Temporal dynamics of the expression of SLCCNV AV1 and BV1 and their ratio in plants

At 6 dpi, very low level of AV1 and BV1 expression was detected and thus only the data at 12, 18, 24 and 30 dpi were analyzed and presented. The relative expression of AV1 was similar at 12 and 18 dpi, decreased significantly from 18 to 24 dpi, and did not show significant changes from 24 to 30 dpi (Fig. 5A). The relative expression of BV1 decreased from 12 to 24 dpi and then did not show significant change from 24 to 30 dpi (Fig. 5B). The ratio of BV1 expression to that of AV1 decreased significantly from 12 to 18 dpi, increased from 18 to 24 dpi and then did not show significant change from 24 to 30 dpi (Fig. 5C). These data indicate that dpi significantly affected the expression of AV1 and BV1 and their ratio (Fig. 5; Table 2).

Fig. 5
figure 5

Temporal dynamics of expression of SLCCNV AV1andBV1and their ratio in zucchini plants. Zucchini plants were inoculated with a mixture of equal quantities of agrobacteria containing infectious clones of SLCCNV DNA-A and DNA-B; relative expressions of AV1 and BV1 and BV1/AV1 expression ratio (A, B and C). Values are means ± SEM. The number of plants tested was 12 at 12, 18 and 24 dpi, and 11 at 30 dpi. Different letters above the columns indicate significant differences (one-way ANOVA, P < 0.05)

Temporal dynamics of the quantities of SLCMV DNA-A and DNA-B and their ratio in plants

When tobacco plants were inoculated with a mixture of equal quantities of agrobacteria containing infectious clones of SLCMV DNA-A and DNA-B, dpi significantly affected the quantities of DNA-A and DNA-B, but not their ratio (Fig. 6; Table 2). Specifically, DNA-A quantity increased from 6 to 12 dpi, and did not change significantly from 12 to 18 dpi, and then increased significantly from 18 to 24 dpi (Fig. 6A). DNA-B quantity increased from 6 to 12 dpi, and did not show significant change from 12 to 18 dpi, and then increased from 18 to 24 dpi (Fig. 6B).

Fig. 6
figure 6

Temporal dynamics of quantities of SLCMV DNA-A and DNA-B and their ratio in tobacco plants. Tobacco plants were inoculated with a mixture of equal quantities of agrobacteria containing infectious clones of SLCMV DNA-A and DNA-B, and then sampled at designated dpi for the quantification of DNA-A (A), DNA-B (B) and B/A ratio (C). Values are means ± SEM. The number of plants tested was 14 at 6, 12 and 18 dpi, and 19 at 24 dpi. Different letters above the columns indicate significant differences (one-way ANOVA, P < 0.05)

Discussion

In this study, we characterized the ratio between the two genomic components of bipartite begomoviruses. We first showed that when two plant species were inoculated with a mixture of equal quantities of agrobacteria inoculum containing infectious clones of SLCCNV DNA-A and DNA-B, the quantities of DNA-A and DNA-B varied significantly with time; However in each of the two plant species the B/A ratio remained constant during SLCCNV infection (Figs. 1 and 2; Table 2). We then inoculated zucchini plants with agrobacteria inocula of different B/A ratios of SLCCNV and observed the B/A ratio during infection. We demonstrated that while the variation of B/A ratio in agrobacteria inoculum may affect the B/A ratio in plants in the initial phase of infection, the effects subsequently disappeared (Fig. 3; Table 2). We further showed while the quantities of DNA-A and DNA-B in the plants infected by agrobacteria–inoculum, sap-mediated inoculation, and whitefly-mediated transmission differed significantly, the B/A ratios did not differ among plants infected with the three modes of virus inoculation (Fig. 4; Table 2). We next found that the expression of the gene BV1 encoded by DNA-B and the gene AV1 encoded by DNA-A and their ratio varied significantly during infection (Fig. 5; Table 2). Finally, we tested another bipartite begomovirus SLCMV, and similarly found that while the quantities of DNA-A and DNA-B varied significantly during infection, the B/A ratio remained constant (Fig. 6; Table 2).

For bipartite begomoviruses, DNA-A and DNA-B share no sequence similarity, except for a highly conserved common region of ∼200 nt that contains the viral origin of replication [18, 19]. During replication, Rep proteins encoded by DNA-A bind specifically to the replication origin of viral DNAs and initiate rolling circle replication [27, 28]. In this study, we found that B/A ratio of bipartite begomoviruses in plants is constant under various scenarios, indicating that the accumulation of DNA-A and DNA-B is tightly associated. The tightly-associated accumulation of the two genomic components may be attributable to the activity of Rep proteins, which are responsible for the replication of both DNA-A and DNA-B. Notably, B/A ratio of SLCCNV in two plant species is always higher than 1, suggesting that more DNA-Bs were constantly generated than DNA-As. For another bipartite begomovirus SLCMV, the B/A ratio in plant is ∼1. These results suggest that at least for some bipartite begomoviruses, Rep proteins do not replicate DNA-A and DNA-B randomly, as random replication of DNA-A and DNA-B is likely to result in similar quantities of the two genomic components. As the interaction between Rep proteins and common regions dictates viral DNA replication, we propose that higher accumulation of DNA-B than DNA-A may be due to the different affinity of the two common regions to Rep proteins. Additionally, the fact that DNA-A-encoded Rep proteins replicate DNA-B with efficiencies equal to or higher than that of DNA-A, indicates that the two genomic components of bipartite begomoviruses have adapted well to each other. Considering the disparate evolutionary origins of the two genomic components [18,19,20], it would be of interest to decipher how they have adapted to each other in the long-term association.

As the genomic components of viruses serve as the template for viral gene expression, changes in the copy number of certain viral genomic components may directly impact viral gene expression. In theory, multipartitism may promote the regulation of the expression of genes encoded by different genomic components [3, 4]. By changing the ratio of genomic components, namely the relative copy number of viral genes, the genes encoded by over-represented genomic components may be expressed preferentially to facilitate adaptation [3, 29]. For example, changes in the frequency of a given genomic component positively and linearly affected the expression of the corresponding genes in FBNSV, a multipartite nanovirus [29]. In our study, B/A ratio remained constant and between 3.42 and 4.26 in zucchini plants, indicating higher accumulation of DNA-B than that DNA-A. On the contrary, characterization of mRNA level shows higher expression of DNA-A-encoded AV1 than DNA-B-encoded BV1 at all the time points analyzed. While more investigations are required to explore the mechanisms underlying the preferred expression of AV1, we propose that the activity of promoters of the two viral genes may play a role as they directly dictate viral gene expression. More importantly, while the B/A ratio at DNA level remained constant, BV1/AV1 expression ratio changed significantly during infection. Together, it seems that unlike FBNSV, B/A ratio of bipartite begomoviruses in plants may play a minor role, if any, in modulating the expression of genes encoded by the two genomic components. The reason of divergence may be that while each of the eight genomic components of FBNSV encodes only one gene, each genomic component of bipartite begomoviruses encodes more than one gene [8, 11]. Under this scenario, an exquisite and finely-tuned gene expression regulation system, instead of the simple B/A ratio in plants, is required for the precise control of viral gene expression in bipartite begomoviruses. Additionally, these findings suggest that cautions should be taken when applying knowledge obtained from multipartite viruses to studies of bipartite begomoviruses.

Taken together, here we have found that the B/A ratio of bipartite begomoviruses in plants remains constant during infection. Our results show that bipartite begomoviruses exert control on the B/A ratio in plants as the variation of B/A ratio in agrobacteria inoculum, sap transmission and whitefly-mediated virus transmission can barely affect the B/A ratio in plants as the infection progresses. Our findings uncovered a new feature of the genomes of bipartite begomoviruses, namely a constant quantitative relationship between the two genomic components during infection and transmission. This feature may play a role in the biology of bipartite begomoviruses, and should be subjected to further investigation for its contribution to viral pathogenesis.

Data Availability

The datasets supporting the conclusions of this article are included within the article and its additional files.

References

  1. World Health Organization. https://www.who.int/emergencies/diseases/novel-coronavirus-2019. Accessed 5 April 2023.

  2. Jones RAC. Global plant virus disease pandemics and epidemics. Plants. 2021;10(2):233.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Sicard A, Michalakis Y, GutieÂrrez S, Blanc S. The strange lifestyle of multipartite viruses. PLoS Pathog. 2016;12(11):e1005819.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Zwart MP, Elena SF. Modeling multipartite virus evolution: the genome formula facilitates rapid adaptation to heterogeneous environments. Virus Evol. 2020;6(1):veaa022.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Hu ZY, Zhang XL, Liu W, Zhou Q, Zhang Q, Li GH, Yao Q. Genome segments accumulate with different frequencies in Bombyx mori bidensovirus. J Basic Microbiol. 2016;56(12):1338–43.

    Article  CAS  PubMed  Google Scholar 

  6. Wu BL, Zwart MP, Sánchez-Navarro JA, Elena SF. Within-host evolution of segments ratio for the tripartite genome of alfalfa mosaic virus. Sci Rep. 2017;7:5004.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Moreau Y, Gil P, Exbrayat A, Rakotoarivony I, Breard E, Sailleau C, Viarouge C, Zientara S, Savini G, Goffredo M, Mancini G, Loire E, Gutierrez S. The genome segments of bluetongue virus differ in copy number in a host-specific manner. J Virol. 2020;95(1):e01834–20.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Sicard A, Yvon M, Timchenko T, Gronenborn B, Michalakis Y, Gutierrez S, Blanc S. Gene copy number is differentially regulated in a multipartite virus. Nat Commun. 2013;4:2248.

    Article  PubMed  Google Scholar 

  9. Valdano E, Manrubia S, Gómez S, Arenas A. Endemicity and prevalence of multipartite viruses under heterogeneous between-host transmission. PLoS Comput Biol. 2019;15(3):e1006876.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Zhang YJ, Wu ZX, Holme P, Yang KC. Advantage of being multicomponent and spatial: multipartite viruses colonize structured populations with lower thresholds. Phys Rev Lett. 2019;123(13):138101.

    Article  CAS  PubMed  Google Scholar 

  11. Fiallo-Olivé E, Lett JM, Martin DP, Roumagnac P, Varsani A, Zerbini FM, Navas-Castillo J, ICTV Report Consortium. ICTV virus taxonomy profile: Geminiviridae 2021. J Gen Virol. 2021;102(12):001696.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Fiallo-Olivé E, Navas-Castillo J. 2023. Begomoviruses: what is the secret(s) of their success? Trends Plant Sci. 2023;8(6):715 – 27.

  13. The International Committee on Taxonomy of Viruses. https://talk.ictvonline.org/ictv-reports/ictv_online_report/ssdna-viruses/w/geminiviridae/392/genus-begomovirus. Accessed 15 November 2022.

  14. Wang D, Huang GX, Shi T, Wang GF, Fang FR, Zhang X, Ye J. Surveillance and distribution of the emergent sri lankan cassava mosaic virus in China. Phytopathol Res. 2020;2:18.

    Article  Google Scholar 

  15. Wu HJ, Li M, Ni H, Bin P, Gu QS. Molecular and biological characterization of melon-infecting squash leaf curl China virus in China. J Integr Agr. 2020;19:570–7.

    Article  CAS  Google Scholar 

  16. Saokham K, Hemniam N, Roekwan S, Hunsawattanakul S, Thawinampan J, Siriwan W. Survey and molecular detection of sri lankan cassava mosaic virus in Thailand. PLoS ONE. 2021;16:e0252846.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Qiu YH, Zhang HJ, Tian W, Fan LG, Du MX, Yuan GL, Wang DX, Wen CL, Xu XL. First report of squash leaf curl china virus infecting tomato in China. Plant Dis. 2022;106(9):2539.

    Article  Google Scholar 

  18. Nawaz-ul-Rehman MS, Fauquet CM. Evolution of geminiviruses and their satellites. FEBS Lett. 2009;583:1825–32.

    Article  CAS  PubMed  Google Scholar 

  19. Briddon RW, Patil BL, Bagewadi B, Nawaz-ul-Rehman MS, Fauquet CM. Distinct evolutionary histories of the DNA-A and DNA-B components of bipartite begomoviruses. BMC Evol Biol. 2010;10:97.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Saccardo F, Cettul E, Palmano S, Noris E, Firrao G. On the alleged origin of geminiviruses from extrachromosomal DNAs of phytoplasmas. BMC Evol Biol. 2011;11:185.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Xavier CAD, Godinho MT, Mar TB, Ferro CG, Sande OFL, Silva JC, Ramos-Sobrinho R, Nascimento RN, Assuncao I, Lima GSA, Lima ATM, Zerbini FM. Evolutionary dynamics of bipartite begomoviruses revealed by complete genome analysis. Mol Ecol. 2021;30(15):3747–67.

    Article  CAS  PubMed  Google Scholar 

  22. Noueiry AO, Lucas WJ, Gilbertson RL. Two proteins of a plant DNA virus coordinate nuclear and plasmodesmal transport. Cell. 1994;76(5):925–32.

    Article  CAS  PubMed  Google Scholar 

  23. Saunders K, Salim N, Mali VR, Malathi VG, Briddon R, Markham PG, Stanley J. Characterization of sri lankan cassava mosaic virus and indian cassava mosaic virus: evidence for acquisition of a DNA B component by a monopartite begomovirus. Virology. 2002;293(1):63–74.

    Article  CAS  PubMed  Google Scholar 

  24. Weigel K, Pohl JO, Wege C, Jeske H. A population genetics perspective on geminivirus infection. J Virol. 2015;89(23):11926–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Chi Y, Pan LL, Bouvaine S, Fan YY, Liu YQ, Liu SS, Seal S, Wang XW. Differential transmission of sri lankan cassava mosaic virus by three cryptic species of the whitefly Bemisia tabaci complex. Virology. 2019;540:141–49.

    Article  PubMed  Google Scholar 

  26. Pan LL, Chi Y, Liu C, Fan YY, Liu SS. Mutations in the coat protein of a begomovirus result in altered transmission by different species of whitefly vectors. Virus Evol. 2020;6(1):veaa014.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Fontes EPB, Eagle PA, Sipe PS, Luckow VA, Hanley-Bowdoin L. Interaction between a geminivirus replication protein and origin DNA is essential for viral replication. J Biol Chem. 1994;269(11):8459–65.

    Article  CAS  PubMed  Google Scholar 

  28. Laufs J, Traut W, Heyraud F, Matzeit V, Rogers SG, Schell J, Gronenborn B. In vitro cleavage and joining at the viral origin of replication by the replication initiator protein of tomato yellow leaf curl virus. Proc Natl Acad Sci U S A. 1995;92(9):3879–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Gallet R, Mattia JD, Rave S, Zeddam J, Vitalis R, Michalakis Y, Blanc S. Gene copy number variations at the within-host population level modulate gene expression in a multipartite virus. Virus Evol. 2022;8(2):1–11.

    Article  Google Scholar 

Download references

Acknowledgements

We thank Ms. Xiao-Xiao Feng from the Agricultural Experiment Station, Zhejiang University for greenhouse management.

Funding

Financial support for this study was provided by the National Key Research and Development Program (2021YFC2600100), the National Natural Science Foundation of China (31930092), the Student Research Training Program of Zhejiang University (Y202204167) and the earmarked fund for China Agriculture Research System (CARS-23-C05).

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization, YXX and LLP; Data curation, YXX and LLP; Formal analysis, YXX, YJW and LLP; Funding acquisition, SSL and LLP; Investigation, YXX, DL and YJW; Methodology, YXX, DL and YJW; Project administration, YXX, DL and YJW; Supervision, SSL and LLP; Validation, LLP; Writing—review & editing, SSL and LLP. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Li-Long Pan.

Ethics declarations

Competing interests

The authors declare no conflict of interest.

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors approved the manuscript and gave consent for publication.

Additional information

Publishers Note

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

Electronic supplementary material

Rights and permissions

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xiao, YX., Li, D., Wu, YJ. et al. Constant ratio between the genomic components of bipartite begomoviruses during infection and transmission. Virol J 20, 186 (2023). https://doi.org/10.1186/s12985-023-02148-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12985-023-02148-2

Keywords