Molecular Characterization of Isolates of the <em>Banana bunchy top virus</em> from the District of Chókwè, Mozambique

Banana bunchy top virus (BBTV) was recently detected in Mozambique and appears to be limited to the provinces of Gaza, Maputo and Zambezia, but it has great potential to spread to other provinces. Despite its importance, nothing is known about the BBTV isolates that occur in Mozambique. In this study, the sequences of the S and R genes of 40 isolates chosen as representatives of samples collected previously from eleven farms of the four administrative posts of the district of Chóckwè, Gaza province, were sequenced and analyzed. The S-DNA nucleotide sequences of the analyzed isolates were highly conserved, with identity ranging from 97% to 100%. The same was observed for the R-DNA sequences, with most identities ranging between 98% and 100% among the isolates from Chókwè and above 90% when compared to the isolates from GenBank. The phylogenetic analysis showed that the Mozambican BBTV isolates belong to the Pacific-Indian Oceans (PIO) group, showing greater proximity to the isolate JQ820453 from Malawi than to the isolates from sub-Saharan countries, which were grouped in a distinct subclade. This is the first study conducted to determine the molecular characteristics of BBTV isolates present in Mozambique.

Keywords:

nanovirus

;

Pacific-Indian Oceans (PIO)

;

South-east Asian (SEA)

;

Phylogeny

;

DNA S

;

DNA R

;

Sub-Saharan African (SSA)

Subject:

Biology and Life Sciences - Biology and Biotechnology

1. Introduction

Bananas and plantains are important food crops and daily sources of carbohydrates for at least 100 million people on the African continent [1,2,3]. The importance of bananas as a staple food is exemplified by Uganda, where per capita consumption reaches approximately 223 kg per person per year. Of the 20 million tons of bananas produced in East and Southern Africa, only 1% is exported, mainly by Zimbabwe, Uganda and South Africa. Approximately 4% of the 11 million tons produced in West and Central Africa are exported, mainly by Côte d’Ivoire and Cameroon [4,5].

In Africa, the production of bananas and plantains is often hampered by biological and environmental constraints. Many pests and diseases have been introduced to the continent through infected planting material. With the expansion of the banana industry in African countries, these diseases and pests spread to several new areas. In many areas, disease-causing pathogens occur in epidemic proportions. The enormous costs involved in controlling some of these diseases make it almost impossible to continue banana production in severely affected areas [6].

In Mozambique, bananas (Musa spp.) are produced mainly on small farms; however, some large farms are found around the provinces of Maputo Province, Gaza, Nampula and Manica. These farms mainly cultivate Cavendish varieties and supply large cities such as Maputo and Beira and the markets of South Africa, Botswana and the Kingdom of eSwatini. Currently, banana plantations (farms) occupy 62,000 ha in Mozambique, with an estimated production of 470,000 thousand tons per year [7].

Production on small farms in southern Mozambique (Maputo Province and Gaza Province) averages 45 tons per ha for a production area of 3216 ha. In addition to the various pests (nematodes and weevils) and diseases (Sigatoka, Fusarium disease or Mal de Panamá - Foc TR4 and currently bunchy top da Bananeira) that can affect banana crops, the lack of production knowledge, low performance, low soil fertility and long drought periods (particularly in the south of the country) are factors that hinder crop management [7]. Among the diseases cited, Fusariosis or Panama disease (FocTR4) and Banana bunchy top disease (BBTD) have been the most important.

BBTD, whose etiological agent is the Banana bunchy top virus (BBTV), was introduced and reported for the first time in Sub-Saharan Africa (SSA) in the 1960s in the Democratic Republic of Congo [8]. Since then, the spread of BBTV has been confirmed in 16 countries in the Central, Southern and West African regions, and it was first recorded in South Africa in 2015 [9] and in Mozambique in June 2016, in the Chókwè irrigated perimeter, more specifically in the Primeira Zona locality, in Gaza province [7].

BBTV belongs to the domain Monodnariria, kingdom Shotokuvirae, phylum Cresdnaviricota, class Arfiviricetes, order Mulpavirales, family Nanoviridae and genus Babuvirus [10]. BBTV affects only members of the genus Musa and is transmitted in a persistent circulatory manner by its only known vector, Pentalonia nigronervosa Coq. [11,12]. The BBTV genome has six circular single-stranded DNA (ssDNA) components (R-DNA, U3-DNA, S-DNA, M-DNA, C-DNA and N-DNA), each approximately 1 kilobase (kb) in size [13]. Each component potentially encodes a protein, except for R-DNA, which encodes two proteins linked to viral replication [14,15].

S-DNA encodes the capsid protein (CP), which forms isometric particles and aids in genome packaging [16]. C-, N- and M-DNA encode cell cycle binding proteins (Clink), nuclear transport proteins (NSPs) and movement proteins (MPs), respectively (17-19). The function of U3-DNA remains unknown for BBTV. Rep and CPs are evolutionarily more conserved than the other proteins and therefore serve as useful markers for the identification and classification of ssDNA viruses [20].

The R gene has been used to classify isolates from different countries located in distinct geographical regions [21]. There are two suggested groups, the Pacific-Indian Oceans (PIO) and South-east Asian (SEA) groups, characterized by an approximately 10% nucleotide difference, while the intragroup difference ranges from 1.9 to 10% [21].

In Mozambique, the genetic characteristics of BBTV isolates occurring in producing fields located in the south of the country have not yet been reported. These data could allow inferences about the possible geographical origin of the viral isolates present and support preventive measures related to the control of viral spread in Africa. In this study, the R-DNA and S-DNA of 40 virus isolates, representing the 175 isolates previously collected in 23 locations of 11 farms located in the district of Chókwè, were sequenced and analyzed. The results obtained are presented and discussed here.

2. Materials and Methods

2.1. BBTV Isolates Analyzed

The 40 BBTV isolates used in this study were chosen proportionally to the number of isolates from 23 fields located on 11 farms in the 4 administrative posts of the Chókwè district (initial site of the outbreak), Gaza Province (Figure 1b) in Mozambique (Figure 1a), with the help of researchers from the Institute of Agricultural Research of Mozambique - IIAM and commercial and family farmers involved in banana production. Located in southern Gaza Province, the district of Chókwè (Figure 1c) has the following geographical coordinates: 24°05′ and 24°48′ South latitude; 32°33′ and 33°35′ East Longitude [22].

2.2. Amplification and Analysis of the S and R Genes

Total plant DNA was extracted from leaf tissue in the laboratory of the Biotechnology Center of the Eduardo Mondlane University of Mozambique (UEM), following the protocol of Lodhi et al. [23]. For extraction, 0.35 g of fresh plant tissue was macerated in liquid nitrogen and then 1.5 mL of 2% CTAB buffer (100 mM Tris-HCl, pH 8.0; 20 mM EDTA; 1.4 M NaCl; 80 mM Na2SO3; 2% PVP-10 and 2% cetyl trimethyl ammonium bromide (CTAB) containing 0.2% β mercaptoethanol. After maceration, the tissue was transferred to a 1.5 mL sterile microcentrifuge tube. Then, 750 µl of extraction buffer with 2% PVP was added and incubated at 60°C for 30 minutes, and the tubes were mixed by inversion every ten minutes. A volume of 750 µl of chloroform:isoamyl alcohol solution (24:1) was added to the tubes and centrifuged at 12,000 rpm for 10 minutes at 23°C. The aqueous phase (supernatant) was transferred to new Eppendorf tubes, and the DNA was precipitated by adding 0.6 volumes of isopropanol and incubated at -20°C for one hour. After centrifugation at 12,000 rpm for 10 minutes, the supernatant was discarded, and the DNA was washed with 500 µl of 70% ethanol and resuspended in 100 µl of 1X TBE buffer. The tubes containing the extracted DNA were placed on dry ice and transported to Brazil to the Laboratory of Molecular Virology, Department of Plant Pathology, Federal University of Lavras (UFLA), for further studies.

In the DNA amplification reaction, primers specially designed to flank the S gene region and the R gene region were employed, as specified in Table 1.

DNA amplification was performed in a 25 μl reaction consisting of 2.5 μl of 10X PCR buffer containing 15 mM MgCl2, 2.0 μl of 10 mM dNTPs, 1.0 μl of forward primer and 1.0 μl of reverse primer at a concentration of 10 μM, 0.25 μl of Taq DNA polymerase (Sigma, USA), 1.0 μl of DNA and enough water for a total reaction volume of 25 μl. Amplification was performed in a Master Cycler Thermocycler (Eppendorf, Germany) with the following program: 94°C for 2 minutes followed by 35 cycles of 94°C for 30 seconds, 49°C (DNA S) or 52°C (DNA R) for 20 seconds and 72°C for 60 seconds, with a final extension of 72°C for 5 minutes. The amplified products were analyzed by electrophoresis in a 1% agarose gel and counterstained with Gel Red (Biotium). The bands were visualized and documented using an Alpha Imager (Alpha Innotech Corp., USA).

2.3. Sequencing and Analysis of S and R DNA Sequences

The PCR products for each of the R-DNA and S-DNA (CP) fragments were purified and sent for sequencing at the Brazilian company FIOCRUZ. The sequence data were analyzed using the BioEdit software program (version 7.0.90) and NCBI (National Center for Biotechnology Information). The multiple alignments of the nucleotide and amino acid sequences of the studied isolates with other viral isolates available in GenBank (Table 2), as well determination of the identity between them, were performed using the program CRUSTAL W2 (V.2.0). Isolates from the PIO group, including some isolates from sub-Saharan Africa, and from the SEA group were chosen for comparison. An isolate of pineapple bunchy top virus (ABTV). The phylogenetic trees were constructed by the neighbor-joining (NJ) method in the MEGA7.02 program [24] using 2,000 bootstrap replicates.

3. Results

3.1. Amplification of the S and R Genes

The 40 BBTV isolates from eleven farms of the four administrative posts of the district of Chóckwè, Gaza province, Mozambique, were sequenced and analyzed. The bands amplified from the S gene can be seen in Figure 2 and those from the R gene in Figure 3.

3.2. Analysis of the BBTV S-DNA Sequences

The nucleotide identities between the 40 isolates collected in Mozambique (Table 3) and the isolates available in GenBank (Table 2) showed that this gene is highly conserved, with more than 50% of the isolates having 100% identity, and the lowest identity was 98%. When compared to the GenBank isolates, the lowest identities observed (92-93%) occurred with isolates MT433376 and MT433375 from Indonesia, KM607468 from Taiwan, KM607469 from the Philippines, KM607536 and AF238876 from China and AB078023 from Japan. The highest identities occurred with isolates JQ820467 from Rwanda, KM607470 from Egypt, KM607505 from the Democratic Republic of Congo (DRC), and JF755980 and JQ820455 from Malawi, all from the African continent, although from different geographical regions. Only Malawi is found in the surrounding region, bordering western Mozambique. These isolates present in the African continent probably had a common geographical origin.

The amino acid identities observed among the Mozambican isolates were similar, with only seven isolates varying between 96 and 99% and the others showing 100% identity. When compared to isolates from GenBank, in addition to the higher identities already specified above, they also showed high identities with isolates from South and Southeast Asia and South Pacific: MK140621 from Pakistan, AB252642 from Miamar and KM607584 from Tonga. Notably, isolates BBTV_172_F11 from Mozambique and AB078023 from Japan showed 98% and 93% nucleotide identity, respectively, when compared to the other isolates. The amino acid identity of the Mozambican isolate was between 95 and 97%, indicating that the nucleotide substitutions were of the synonymous type, while the Japanese isolate presented amino acid identities between 95 and 98%, indicating that the substitutions were of the synonymous type.

Figure 4 and Figure 5 illustrate the phylogenetic trees constructed based on the nucleotide and amino acid sequences, respectively, of the S-DNA. There was a clear subdivision of isolates into two clades, one with isolates from China, Japan, the Philippines, Taiwan and Indonesia and one with Mozambican isolates and the remaining isolates from GenBank.

The Mozambican isolates were grouped into the same clade as the isolates from Myanmar, Malawi, Rwanda and DRC. Although isolates from Myanmar, Australia, Burundi, Egypt, Tonga, Pakistan and China were grouped in the same clade as the Mozambican isolates, they occupied different subclades.

In the tree constructed based on amino acid sequences, in addition to the isolates from Malawi and Rwanda, the isolates from Myanmar, Pakistan, Tonga and Egypt were also grouped with the Mozambican isolates. This must have occurred because the nucleotide substitutions are synonymous. In this tree, some isolates moved and occupied other subclades. Isolate BBTV_271_F1 occupied a subclade together with isolate AB252644 from Miami and isolate KM607505 from the DRC; isolate BBTV_292_F2 was placed in a subclade close to isolates from China, Japan, the Philippines, Indonesia and Taiwan, while isolates BBTV_286_F2 and BBTV_179_F11 were placed in different subclades and clades.

3.3. Analysis of the BBTV R DNA Sequences

Like the S-DNA, the R-DNA gene was also highly conserved. Except isolates BBTV_286_F2 and BBTV_273_F1, which presented identities between 93 and 96% when compared to each other, all the other isolates presented identities between 98% and 100%. When compared with the PIO isolates from GenBank, the identity of these two isolates ranged from 94% to 96%, while all other Mozambican isolates exhibited identities between 96% and 99%. On the other hand, when compared with the SEA isolates, the two isolates presented identities between 87 and 89%, while the others presented identities between 90 and 93%.

When comparing the identities between the PIO isolates from GenBank, it was observed that they have 97% to 99% identity among themselves and between 90 and 93% identities with the SEA isolates. These results allowed us to infer that the isolates that occur in Mozambique are of the PIO type, i.e., they belong to the Indian Pacific Ocean group.

The amino acid identities between the two isolates cited above (BBTV_273_F1 and BBTV_286_F2) and the other Mozambican isolates were lower than that of the nucleotides, between 91 and 93%, indicating that the nucleotide substitutions were nonsynonymous.

The same was observed with isolate 3, whose nucleotide identity ranged from 94% to 99% and amino acid identity ranged from 91% to 97%. Most isolates presented identities of 100% when compared to each other and from 97% to 100% when compared to the isolates from the PIO group and 91% to 95% when compared to the isolates from the SEA group.

The phylogenetic tree based on the nucleotide sequences is illustrated in Figure 6. The isolates were grouped into two distinct clades, one with the isolates of the PIO group and the other with the isolates of the SEA group. All Mozambican isolates were grouped into a subclade with the Malawi isolate JQ820453. The other isolates of the PIO group were grouped in the same clade but in a distinct subclade.

In the phylogenetic tree based on amino acid sequences (Figure 7), the distribution was slightly different. There was also a separation of the isolates into two clades, one with the isolates of the PIO group and the other with the isolates of the SEA group. However, in addition to isolate JQ820453 from Malawi, isolates KM607635 from DRC, JQ820465 from Rwanda and, a little further away, AF416467 from Tonga also grouped together with most of the Mozambican isolates. Another eight isolates from Mozambique were mixed with isolates from different parts of the world in different combinations. Both the identity and phylogenetic results support the classification of the studied isolates within the PIO group.

4. Discussion

In this study, the R-DNA and S-DNA genes of Mozambican BBTV isolates from the Chókwè district of Gaza were sequenced and analyzed for the first time, revealing new information about their classification and possible geographical origin.

The nucleotide identities of the R-DNA among the 40 analyzed isolates revealed that 38 of them presented a difference of at most 2%, showing the high conservation of this gene, which has already been observed by authors from other regions of the world. Adegbola [37] in Nigeria also found 100% identities among the local isolates. Another study conducted in Cameroon by Oben et al. [38] showed a similar result, that is, 100% homology between nucleotide sequences of the analyzed isolates.

In phylogenetic studies performed by Karan et al. [39] and reviewed by Yu et al. [20], the classification of BBTV isolates, based on nucleotide sequences, into two groups, South Pacific group and the Asian group, was proposed. The Mozambican isolates showed higher identities with isolates from the PIO group, between 97% and 99%, and the lowest identities with the isolates from the SEA group, between 90% and 92%, indicating that the Mozambican isolates belong to the first group. These results were corroborated by the distribution of the isolates in the phylogenetic tree based on nucleotide sequences, which showed a clear subdivision of the isolates into two distinct clades, with all the Mozambican isolates grouped together with the isolates of the PIO group. Other authors studying isolates from different African countries also reported that the isolates belonged to the PIO group [21,37,38,40,41,42,43].

Interestingly, the Mozambican isolates were very similar to isolate JQ820453 from Malawi, which borders Mozambique to the south and west. This demonstrates that there is a high probability that the Mozambican isolates came from this country. Kumar et al. [8] had already considered Malawi as an area for the dissemination of BBTV through the exchange of planting material of preferred cultivars among farmers in bordering countries. However, it is necessary to consider that South Africa also has the presence of BBTV and shares a border with southern Mozambique. As there is no information on the complete gene sequences of the R-DNA in GenBank, it could not be included in this study. However, South Africa cannot be ruled out as the origin of the inoculants that reached Malawi and Mozambique, as it is geographically closer to southern Mozambique than Malawi.

Research on the introduction of BBTV in sub-Saharan Africa (SSA) indicated that the arrival of BBTV in Africa may have occurred in two ways: in the DRC in 1950, it was probably through infected propagules brought from southern Asia or the South Pacific [44,45,46]; in Equatorial Guinea and Gabon, it may have been introduced by infected aphids (Pentalonia nigronervosa) brought by migrant workers from the Philippines. Then, it is assumed that it has spread to other African countries (Rwanda, Burundi and Central African Republic) [46,47].

This virus is known to be widely prevalent in Central African countries and in Malawi, Southern Africa [8]. To date, the occurrence of the disease has been reported in 16 countries of Central and Southern Africa (Angola, Benin, Egypt, Nigeria, Gabon, Burundi, Cameroon, Central African Republic, DRC, Congo, Equatorial Guinea, Malawi, Rwanda, Zambia, South Africa and Mozambique) [48].

In the nucleotide sequence analysis of the S gene of the isolates from Mozambique, it was observed that it is also highly conserved, with more than 50% of the isolates being 100% identical, and the lowest identity was 98%. Additionally, in the case of the S gene, the greatest nucleotide identities were observed with the isolates of the PIO group, especially the African isolates JQ820467 from Rwanda, KM607470 from Egypt, KM607505 from the DRC, and JF755980 and JQ820455 from Malawi. The phylogenetic tree constructed based on nucleotide sequences, similar to that observed for the R gene, also indicated similarity to the isolates of the PIO group, especially isolates RJQ820467 from Rwanda and JQ820455 from Malawi.

Considering the danger that this disease poses to the banana industry and to food safety, urgent measures are needed to combat it, such as the installation of indexing laboratories and the availability of positive controls to increase confidence in the diagnosis based on PCR. It is also necessary to train the local team involved in teaching and extension activities and the farmers in the recognition and control of the virus. In this context, there is an urgent need to make farmers and government officials aware of the disease and the need to implement control measures, including (i) large-scale production and supply of virus free planting materials of the varieties preferred by farmers to rehabilitate banana production in the affected regions, preventing the use of infected planting material; (ii) raising awareness among farmers of the need to destroy infected material; (iii) greater phytosanitary surveillance and implementation of measures including restrictions on the movement of planting materials from regions affected by the disease, especially those bordering areas affected by BBTD; (iv) implementation of awareness programs among farmers, extensionists and regulatory agencies; (v) training to improve disease monitoring and diagnosis capacity; (vi) application of quarantine standards and integrated farming practices (ICP) aimed at reducing the spread of the disease and the aphid vector Pentalonia niigronervosa; and (vii) promoting and developing robust research programs to identify resistant varieties, different strains of the virus and its origins as well as the genetic diversity existing in each banana production region in Mozambique. Such measures are also recommended by several researchers, such as James et al. [31], Oben et al. [38], Adegbola et al. [37] and Ximba et al. [49].

Considering that the presence of this virus in Mozambique and in neighboring countries constitutes a major threat to commercial crops and family banana producers, whose culture is considered subsistence for this segment of the population in Mozambique, these actions are essential and need to be implemented urgently and then supported by government agencies.

5. Conclusions

The results of this study support the classification of Mozambican isolates within the PIO group, as observed in Malawi. These and other results indicate that the introduction of BBTV in Mozambique from Malawi through the spread of infected material for planting is highly probable.

These results emphasize the need for intensive research to assess the extent of the geographical spread and severity of BBTV in Mozambique and for the implementation of quarantine and phytosanitary measures to prevent the internal and transboundary spread of this virus.

Author Contributions

All authors contributed to the study conception and design. The first draft of the manuscript was written by Sandra Carvalho I. Mussa Barros and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Funding

The authors would like to thank the National Council for Scientific and Technological Development (CNPq), the Coordination for the Improvement of Higher Education Personnel (CAPES), Foundation for Research Support of the State of Minas Gerais (FAPEMIG) for their support in this project.

Data availability

Data are contained within the article.

Declaration

The authors declare that they do not have any potential technical or financial conflict of interest in relation to the entire content of the study presented.

Conflicts of Interest

The authors declare no conflict of interest.

References

Dotto J, Matemu A, Ndakidemi P (2018) Potential of cooking bananas in addressing food security in East Africa. Int J Biosci 13(4):278-294. [CrossRef]
Perrier X, De Langhe E, Donohue M, Lentfer C, Vrydaghs L, Bakry F, Carreel F, Hippolyte I, Horry J-P, Christophe J et al. (2011) Multidisciplinary perspectives on banana (Musa spp.) domestication. Proc Natl Acad Sci 108(28):11311-11318.
Sharrock S, Frison EA (1998) Musa production around the world-trends, varieties and regional importance. In: Network for the Improvement of Banana and Plantain. Montpellier: France.
Karamura EB, Frison E, Karamura DA, Sharrock S (1999) Banana production Systems in Eastern and Southern Africa. In: Picq C, Fouré E, Frison EA (eds.). Bananas and Food security. International Network for the improvement of Banana an Plantain. France: Montpellier.
Viljoen A, Kunert K, Kiggundu A, Escalant JV, Bornman CH (2004) Biotechnology for sustainable banana and plantain production in Africa: the South African contribution. S Afr J Bot 70(1):67–74. [CrossRef]
Chabi MC, Dassou AG, Dossou-Aminon I, Ogouchoro D, Aman BO, Dansi A (2018) Banana and plantain production systems in Benin: ethnobotanical investigation, varietal diversity, pests, and implications for better production. J Ethnobiol Ethnomed 14(1):78. [CrossRef]
MADER - Ministério de Agricultura e Desenvolvimento Rural (2018) Indústria da banana em Moçambique - situação de Foc TR4 e BBTV. Moçambique: MADER.
Kumar PL, Hanna R, Alabi OJ, Soko MM, Oben TT, Vangu GHP et al. (2011) Banana bunchy top virus in sub-Saharan Africa: Investigations on virus distribution and diversity. Virus Res 159(2):171–182. [CrossRef]
Jooste AEC, Wessels N, van der Merwe M (2016) First Report of Banana bunchy top virus in Banana (Musa spp.) from South Africa. Plant Dis 100(6):1251. [CrossRef]
ICTV - International Committee on Taxonomy of Viruses (2021) International Committee on Taxonomy of Viruses 9th Report. https://talk.ictvonline.org/ictv-reports/ictv_9th_report/ssdna-viruses 2011/w/ssdna_viruses/150/nanoviridae-figures Accessed 19 July 2022.
Bell KE, Dale JL, Ha CV, Vu MT, Revill PA (2002) Characterisation of Rep-encoding components associated with banana bunchy top nanovirus in Vietnam. Arch Virol 147(4):695–707. [CrossRef]
Hu JS, Wang M, Sether D, Xie W, Leonhardt KW (1996) Use of polymerase chain reaction (PCR) to study transmission of banana bunchy top virus by the banana aphid (Pentalonia nigronervosa). Ann Appl Biol 128(1):55–64.
Burns TM, Harding RM, Dale JL (1995) The genome organization of banana bunchy top virus: analysis of six ssDNA components. J Gen Virol 76(6):1471–1482. [CrossRef]
Herrera-Valencia VA, Dugdale B, Harding RM, Dale JL (2006) An iterated sequence in the genome of Banana bunchy top virus is essential for efficient replication. J Gen Virol 87(11):3409–3412. [CrossRef]
Horser CL, Karan M, Harding RM, Dale JL (2001) Additional Rep-encoding DNAs associated with Banana bunchy top virus. Arch Virol 146(1):71–86. [CrossRef]
Ji X, Yu N, Qu L, Li B, Liu Z (2019) Banana bunchy top virus (BBTV) nuclear shuttle protein interacts and re-distributes BBTV coat protein in Nicotiana benthamiana. 3 Biotech 9(4). [CrossRef]
Amin I, Ilyas M, Qazi J, Bashir R, Yadav JS, Mansoor S et al. (2010) Identification of a major pathogenicity determinant and suppressors of RNA silencing encoded by a South Pacific isolate of Banana bunchy top virus originating from Pakistan. Virus Genes 42(2):272–281. [CrossRef]
Yu N, Xie H, Yu Liang Zhang, Wang J, Zhou X, Zhi Xin Liu (2019) Independent modulation of individual genomic component transcription and a cis-acting element related to high transcriptional activity in a multipartite DNA virus. BMC Genomics 20(1). [CrossRef]
Zhuang J, Lin W, Coates CJ, Shang P, Tàiyún Wèi, Wu Z et al. (2019) Cleavage of the Babuvirus Movement Protein B4 into Functional Peptides Capable of Host Factor Conjugation is Required for Virulence. Virol Sin 34(3):295–305. [CrossRef]
Yu NT, Zhang YL, Feng TC, Wang JH, Kulye M, Yang WJ et al. (2012) Cloning and sequence analysis of two Banana bunchy top virus genomes in Hainan. Virus Genes 44(3):488–494. [CrossRef]
Karan M, Harding RM, Dale JL (1994) Evidence for two groups of Banana bunchy top virus isolates. J Gen Virol 75(12):3541–3546. [CrossRef]
UNDP - United Nations Development Programme (2007) MDG Achievement Fund – Environment and Climate Change Thematic Window. MDG-F Terms of Reference. Madrid: MDG-F https://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.462.1402&rep=rep1&type=pdf Accessed 18 July 2021.
Lodhi MA, Ye GN, Weeden NF, Reisch BI (1994) A simple and efficient method for DNA extraction from grapevine cultivars andVitis species. Plant Mol Biol Rep 12(1):6–13. [CrossRef]
Kumar S, Stecher G, Tamura K. (2016) MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol 33(7):1870–1874. [CrossRef]
Islam, MN, Naqvi A, Jan A, Haq Q (2010) Genetic diversity and possible evidence of recombination among Banana Bunchy Top Virus (BBTV) isolates. Int Res J Microbiol 1(1):1-12.
James AP (2011) Viruses of Banana in East Africa: Queensland University of Technology.
Amin I, Qazi J, Mansoor S, Ilyas M, Briddon RW (2008) Molecular characterisation of Banana bunchy top virus (BBTV) from Pakistan. Virus Genes 36(1):191–198. [CrossRef]
Stainton D, Martin DP, Muhire BM, Lolohea S, Halafihi M, Lepoint P et al. (2015) The global distribution of Banana bunchy top virus reveals little evidence for frequent recent, human-mediated long distance dispersal events. Virus Evol 1(1):1-15. [CrossRef]
Furuya N, Kawano S, Natsuaki KT (2005) Characterization and genetic status of Banana bunchy top virus isolated from Okinawa, Japan. J Gen Plant Pathol 71(1):68–73. [CrossRef]
Sharman M, Thomas JE, Skabo S, Holton TA (20078) Abacá bunchy top virus, a new member of the genus Babuvirus (family Nanoviridae). Arch Virol 153(1):135–147. [CrossRef]
James AP, Soko MM, Mugini JA, Geijskes RJ, Dale JL, Harding RM (2012) Molecular characterisation of Banana bunchy top virus isolates from Malawi and Rwanda. Centre for Tropical Crops and Biocommodities. https://www.ncbi.nlm.nih.gov/nuccore/394326320 Accessed 15 July 2021.
Bashir S, Naqvi SMS, Muhammad A, Hussain İ, Ali K, Khan MR et al. (2022) Banana bunchy top virus genetic diversity in Pakistan and association of diversity with recombination in its genomes. Plos One 17(3). [CrossRef]
He Z et al. (2001) Cloning and sequencing of DNA component 1 of two BBTV strains. J Agric Biotechnol 9(2):145-148.
Kawabe K, Onuki M (2002) Banana bunchy top virus coat protein (CP) gene. https://www.ncbi.nlm.nih.gov/nuccore/394326320 Accessed 10 Jan 2023.
Rahayuniati RF, Subandiyah S, Hartono S et al. (2021) Recent distribution and diversity analysis on banana bunchy top virus of banana and alternative host in Indonesia. Trop Plant Pathol 46(5):506–517. [CrossRef]
Furuya N, Natsuaki KT (2008) Identification of banana viruses in Myanmar. https://www.ncbi.nlm.nih.gov/nuccore/394326320. Accessed 21 Jan 2023.
Adegbola RO, Ayodeji O, Awosusi OO, Atiri GI, Kumar PL (2013) First Report of Banana bunchy top virus in Banana and Plantain (Musa spp.) in Nigeria. Plant Dis 97(2):290–290. [CrossRef]
Oben TT, Hanna R, Ngeve J, Alabi OJ, Naidu RA, Kumar PL (2009) Occurrence of Banana Bunchy Top Disease Caused by the Banana bunchy top virus on Banana and Plantain (Musa sp.) in Cameroon. Plant Dis 93(10):1076. [CrossRef]
Karan M (1995) Sequence diversity of DNA fragments associated with Banana bunchy top virus. PhD thesis - Queensland University of Technology, Australia.
Abdel-Salam SM, Dahot MU, Sadik AS (2012) Molecular comparative analysis of component 1 (DNA-R) of an Egyptian isolate of banana bunchy top nanovirus isolated from banana aphid (Pentalonia nigronervosa). J Genet Eng Biotechnol 10(1):55–65. [CrossRef]
Karan M, Harding RM, Dale JL (1997) Association of Banana bunchy top virus DNA components 2 to 6 with bunchy top disease. Mol Plant Pathol 1-7.
Niyongere C, Lepoint P, Losenge T, Blomme G (2015) Towards understanding the diversity of banana bunchy top virus in the Great Lakes region of Africa. Afr J Agric Res 10(7):702–709.
Wanitchakorn R, Harding RM, Dale JL (2000) Sequence variability in the coat protein gene of two groups of banana bunchy top isolates. Arch Virol 145(3):593–602. [CrossRef]
Dale JL (1987) Banana Bunchy Top: An Economically Important Tropical Plant Virus Disease. Adv Virus Res 33:301–325. [CrossRef]
Fouré E, Manser PD (1982) Note sur l’apparition au Gabon d’une grave maladie virale des bananiers et plantains: le bunchy top. Fruits 37(6): 409-414.
Wardlaw CW (1961) Banana diseases including plantain and abaca. London: Longmans, Green.
Jeger MJ, Eden-Green SJ, Thresh JM, Johanson A, Waller JM, Brown AE (1995) Banana Diseases. In: Gowen S (ed.). Bananas and Plantains. London: Chapman and Hall.
EPPO – European and Mediterranean Plant Protection Organization (2021) Global Database. Banana bunchy top vírus (BBTV00). https://gd.eppo.int/taxon/BBTV00/distribution Accessed 15 July 2021.
Ximba SPF, Tshabalala J, Gubba A, Jooste AEC (2022) Monitoring the distribution of Banana bunchy top virus in South Africa: a country-wide survey. Arch Virol 167(6):1433–1441. [CrossRef]

Figure 1. Geographic distribution of BBTV isolates collected in the Chókwè district, Gaza Province, Mozambique (indicated with red dots).

Figure 2. Electrophoretic analysis of PCR products generated by the DNA S gene of BBTV isolates in fourteen different samples: 1: 1Kb marker; line of 2 to 15 samples that tested positive among those analyzed; 16: negative control; 17: positive control.

Figure 3. Electrophoretic analysis of PCR products amplified from the R DNA of BBTV isolates in fourteen different samples. 1: 1Kb marker; 2 to 15: band pattern of the amplified isolates; 16: negative control; 17: positive control.

Figure 4. Phylogenetic tree constructed based on the nucleotide sequence of the S DNA of BBTV isolates collected in the district of Chókwè – Mozambique and isolates from GenBank. Bootstrap values were obtained using the MEGA 7.0 and Neighbor Joining program with 2,000 repetitions.

Figure 5. Phylogenetic tree constructed based on the amino acid sequence of the S DNA of BBTV isolates collected in the district of Chókwè – Mozambique and isolates from GenBank. Bootstrap values were obtained using the MEGA 7.0 and Neighbor Joining program with 2,000 repetitions.

Figure 6. Phylogenetic tree constructed based on the nucleotide sequence of the R DNA of BBTV isolates collected in the district of Chókwè – Mozambique and isolates from GenBank. Bootstrap values were obtained using the MEGA 7.0 and Neighbor Joining program with 2,000 repetitions.

Figure 7. Phylogenetic tree constructed based on the amino acid sequence of the R DNA of BBTV isolates collected in the district of Chókwè – Mozambique and isolates from GenBank. Bootstrap values were obtained using the MEGA 7.0 and Neighbor Joining program with 2,000 repetitions.

Table 1. Primers used for amplification of the R and S Genes with the respective annealing temperatures and amplified fragment.

Primer	Sequence (5′-3′)	Gene	Temp (°C)	Size (bp)
BBTV 67SF	CGTTTAGATGGGTTTGGG	DNA S	49	718
BBTV 784SR	CGTTTAGATGGGTTTGGG	DNA S	49	718
BBTV 32RF	CGGGACGGGACATTTGC	DNA R	52	1007
BBTV 1038SR	CGTTAGGATAGCACAATC	DNA R	52	1007

Table 2. BBTV isolates of R DNA and S DNA and the number of accessions available in GenBank, used for comparison with the Mozambican isolates.

BBTV Group *	DNA R			DNA S
BBTV Group *	Origin	GenBank	Reference	Origin	GenBank	Reference
PIO	Egypt	HQ259074	lslam et al. [25]	Egypt	KM607471	Stainton et al. [28]
PIO	Egypt	AF416465	Horser et al. [15]	Egypt	KM607470	Stainton et al. [28]
PIO	Fiji	AF416466	Horser et al. [15]	Tonga	KM607584	Stainton et al. [28]
PIO	Malawi	JQ820453	James et al. [26]	Malawi	JF755980	Kumar et al. [8]
PIO	Rwanda	JQ820465	James et al. [26]	Rwanda	JQ820467	James et al. [31]
PIO	Tonga	AF416467	Horser et al. [15]	Tonga	KM607550	Stainton et al. [28]
PIO	Pakistan	AM418538	Amin et al. [27]	Pakistan	MK140621	Bashir et al. [32]
PIO	DRC	KM607635	Stainton et al. [28]	DRC	KM607505	Stainton et al. [28]
PIO	Australia	KM607623	Stainton et al. [28]	Australia	KM607489	Stainton et al. [28]
PIO	Burundi	KM607598	Stainton et al. [28]	Burundi	KM607456	Stainton et al. [28]
SEA	China	HQ616074	Yu et al. [20]	China	KM607536	Stainton et al. [28]
SEA	China	HQ378190	Yu et al. [20]	China	AF238876	He et al. [33]
SEA	Japan	AB108458	Furuya, Kawano and Natsuaki [29]	Japan	AB078023	Kawabe and Onuki [34]
SEA	Indonesia	AB186924	Furuya, Kawano and Natsuaki [29]	Indonesia	MT433375	Rahayuniati et al. [35]
SEA	Vietnam	AF416477	Bell et al. [11]	Myanmar	AB252644	Furuya and Natsuaki [36]
SEA	Vietnam	AF416472	Bell et al. [11]	Myanmar	AB252642	Furuya, Kawano and Natsuaki [29]
SEA	Japan	AB108455	Furuya, Kawano and Natsuaki [29]	Malawi	JQ820455	James et al. [31]
SEA	Taiwan	AF416468	Horser et al. [15]	Taiwan	KM607468	Stainton et al. [28]
SEA	Filipinas	AF416469	Horser et al. [15]	Filipinas	KM607469	Stainton et al. [28]
SEA	Indonesia	AB186925	Furuya, Kawano and Natsuaki [29]	Indonesia	MT433376	Rahayuniati et al. [35]
OuGR	Malaysia	EF546813	Sharman et al. [30]	Malaysia	NC_010316	Sharman et al. [30]

*PIO Group according to R DNA analysis, SSA countries in bold.

Table 3. Names of the BSV isolates collected Chokwe district, Gaza Province, Mozambique, with their respective accession number in GEnBAnk.

BSV Isolate Name	GenBank Accession Number
BSV Isolate Name	DNA R	DNA S
BBTV 252 F1	PP097007	PP096947
BBTV 260 F1	PP097009	PP096949
BBTV 261 F1	PP097010	PP096950
BBTV 262 F1	PP097011	PP096951
BBTV 271 F1	PP097012	PP096952
BBTV 273 F1	PP097013	PP096953
BBTV 280 F2	PP096991	PP096931
BBTV 281 F2	PP096992	PP096932
BBTV 285 F2	PP096995	PP096935
BBTV 286 F2	PP096996	PP096936
BBTV 291 F2	PP096999	PP096939
BBTV 292 F2	PP097000	PP096940
BBTV 299 F2	PP097001	PP096941
BBTV 301 F2	PP097003	PP096943
BBTV 308 F2	PP097005	PP096945
BBTV 232 F3	PP096986	PP096926
BBTV 241 F3	PP096988	PP096928
BBTV 243 F3	PP096989	PP096929
BBTV 340 F4	PP096983	PP096923
BBTV 344 F4	PP096984	PP096924
BBTV 329 F5	PP096980	PP096920
BBTV 335 F5	PP096982	PP096922
BBTV 311 F7	PP096979	PP096919
BBTV 201 F8	PP096977	PP096917
BBTV 204 F8	PP096978	PP096918
BBTV 205 F9	PP096975	PP096915
BBTV 210 F9	PP096976	PP096916
BBTV 217 F10	PP096971	PP096911
BBTV 222 F10	PP096972	PP096912
BBTV 224 F10	PP096973	PP096913
BBTV 225 F10	PP096974	PP096914
BBTV 171 F11	PP096954	PP096894
BBTV 172 F11	PP096955	PP096895
BBTV 177 F11	PP096959	PP096899
BBTV 179 F11	PP096961	PP096901
BBTV 180 F11	PP096962	PP096902
BBTV 183 F11	PP096963	PP096903
BBTV 187 F11	PP096967	PP096907
BBTV 191 F11	PP096968	PP096908
BBTV 193 F11	PP096969	PP096909

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

Downloads

152

Views

Comments

Subscription

Notify me about updates to this article or when a peer-reviewed version is published.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

Molecular Characterization of Isolates of the Banana bunchy top virus from the District of Chókwè, Mozambique

Abstract