Annual and geographical variations in the specific composition of jassids and their damage on cotton in Ivory Coast

In recent years, jassids have become a real problem for cotton growing in Ivory Coast. It is important to investigate the causes of this problem. The aim of this study was to highlight the diversity of jassid species and their impact on cotton growing in Ivory Coast. The collections carried out in 2021 identified three species. Jacobiasca lybica (Bergevin & Zanon, 1922) is the most abundant, with proportions ranging from 73.3 to 93.3% depending on the site. The other two species are Empoasca papayae (Oman, 1937) (8.3%) and Empoasca facialis (Jacobi, 1912) (5%). In 2022, collections revealed the invasion of a new species, Amrasca biguttula (Shiraki, 1913), which became dominant with proportions of 90 to 100% depending on the site. Two other species, Jacobiasca lybica (Bergevin & Zanon, 1922) (2.7%) and Empoasca facialis (Jacobi, 1912) (1.3%), cohabit with Amrasca biguttula (Shiraki, 1913). Thus, while the damage noted in 2021 was attributable to Jacobiasca lybica (Bergevin & Zanon, 1922), that observed in 2022 is essentially due to Amrasca biguttula (Shiraki, 1913), with incidences exceeding the economic impact threshold. The North-East of the cotton basin was the area most affected by attacks by these two species. The results of the study reveal significant changes in jassid species composition and climatic conditions in the cotton-growing areas of Ivory Coast, from 1 year to the next. This situation also explains the variations in damage levels.

Ivory Coast is based on a calendar program of six (6) insecticide treatments, every 14 days, from the 45th to the 115th day after cotton emergence.Before starting insecticide treatments, beneficial insects are used to control cotton pests.Thus, starting spraying on the 45th day after cotton emergence enables beneficial insects to naturally regulate cotton pests.One restriction imposed by this strategy is the ban on the use of pyrethroid-based products.This reduces the resistance observed in Helicoverpa armigera Hübner to this chemical family 4 .Pyrethroids are banned at the start of the cropping season until August 10 in northern cotton-growing areas (above the 9th parallel) and until August 20 in southern cotton-growing areas (below the 9th parallel).
Despite these actions, during the 2022-2023 crop year, more severe attacks than in previous years were observed throughout the cotton production zone in Ivory Coast and in several neighboring countries.This damage was estimated at around 34 billion franc of the African Financial Community (f CFA).
According to 1 have shown that pockets of drought during the July-August period favor jassid outbreaks.However, other factors such as substitutions or introductions of new species could also explain these constantly observed changes in behavior.The jassid species known in Ivory Coast are Empoasca facialis (Jacobi, 1912) and Orosius cellulosus (Lindberg, 1927) 5,6 .The first species is characterized by the injection of toxic saliva into the host plant, resulting in leaf discoloration.The second species is known as the vector of a phytoplasma responsible for floral virescence 1,[6][7][8] .However, the main species known to cause the most damage is Empoasca facialis (Jacobi,  1912).Today, we don't know whether this species is still the only or the main cause of the damage observed.Could the emergence of new species in the field be the consequence of the invasion of new species?
The general aim of the present work was to contribute to understanding the causes of regular outbreaks of jassids in cotton crops.Specifically, it aimed to identify wasp species and their dynamics in cotton crops.

Materials
The study focused on leafhopper specimens collected on cotton (Gossypium hirsutum Linné, 1763) in a production environment.The main cotton varieties collected were CI 128, CI 123, Gouassou and Sicama.They were developed by the Centre National de Recherche Agronomique (CNRA).They are hybrid varieties.
The protection program popularized by growers is the six-treatment insecticide window program.This program reduces the use of pyrethroid-based insecticides by positioning alternatives in the first window.The main active ingredients used for crop protection are listed in Table 1.

Sample collections
Specimens were collected for identification during the 2021 and 2022 seasons, from June to November each year.The jassids were caught plant by plant, on a sample of 30 plants taken in groups of five consecutive plants per row, using the sequential "diagonal" method 9 .Adult jassids were collected using a mouth aspirator early in the morning from 7:30 a.m. to 1:00 p.m., then preserved on naphthalene, over a 0.25 ha area marked out in the plot center.

Identification of individuals
Individuals were identified based on taxonomic keys [10][11][12] consolidated by 13 from various resources.They were dissected and then meticulously observed using a Panthera L Motic microscope at 4× and 10× magnifications respectively, and a Motic SMZ-171-TLED magnifying glass.

Observation of jassid damage on cotton plants
Data as collected in collaboration with the cotton companies' Research and Development (R&D) departments.Each year, the monitoring system covered around 320 plots in farming areas.These plots were selected in 32 localities, with ten (10) plots per locality.The plots were chosen according to the most representative sowing dates in each locality (Table 2).An area of 0.25 ha was delimited in the Centre of the plot for data collection.
Each plot was subjected to the technical practices recommended for cotton production in Ivory Coast (variety, sowing density, fertilization, weeding, etc.).On each plot, weekly observations were made from the 30th to the 128th JAL (Jour Après Levée), i.e., for 14 weeks.Observations of plant damage and adult jassids were made plant by plant, on a sample of 30 plants in six groups of five consecutive plants using the diagonal method 9 .A plant is considered attacked when one of its five terminal leaves shows damage or symptoms of jassid attack.Data for 2021 to 2022 were collected using the SySParCotCI data collection application developed by the National Center for Agronomic Research 14 .Data from 2012 to 2020 were collected using the same method, based on the data collection sheets that were used prior to the development of the SySParCotCI data collection application.

Climatic data collection
Rainfall data were provided by the central soil, water, and plant laboratory of the Centre National de Recherche Agronomique (CNRA).These (in-situ) data are taken every day of the year, from weather stations installed not far from the farmers' plots.As for temperature and relative humidity data, we mobilized them from the Nasa power platform (https:// power.larc.nasa.gov/ data-access-viewer/), which provides access to a range of climatic data available at daily, monthly, and annual time steps.For this study, we focused on data from June to July, during the 2021 and 2022 season.

Data analysis
After identification, the percentages of individuals belonging to the different species were calculated according to crop year and locality, and an overall average was determined for both crop years simultaneously.The same analyses were carried out according to locality.Geographical and seasonal damage levels were determined from data collected in 2021 and 2022.Average annual damage levels for the last decade (2012 to 2022) were also determined by combining data collected from 2012 to 2022.Average annual damage levels for the last decade (2012 to 2022) have also been determined by combining data collected from 2012 to 2022.Ecological diversity indices were determined using the Vegan package for R. 15 , based on data from the 2 years combined.These were the species richness index (S), the Shannon index (H), the dominance index (D), the Pielou equitability index (J) and the Jaccard index (I) [16][17][18] .The species richness index (S) represents the number of species in a community, regardless of their relative abundance.The Shannon index (H) is a measure of community stability.It varies according to the number of species and the relative proportion of different species.Simpson index (D) measures the probability that two individuals drawn at random from a sample belong to the same species.Pielou equitability index (J) measures the distribution of individuals within species.It is useful for comparing potential dominance between samples.Jaccard index (I) enables a comparison between two sites, as it evaluates similarity by calculating the ratio between species common to both sites and those specific to each survey.Seasonal variation curves were plotted using Excel software from the Windows Office 2016 package.Analyses and graphical representations were carried out using the free Rstudio development environment of the R 4.2.2 software package 19 using the questionr and ggplot2 packages.The questionr package was used to determine species percentages according to locality and crop year.The ggplot2 package was used for graphical representations.Spatialization of damage and rainfall levels was carried out using Arc gis 10.2.2 software 20 .

Ecological community structure
A total of 600 individuals were identified, including 300 during the 2021 season and 300 during the 2022 season..At the 10 sites sampled, species richness (S) ranged from three (3) to five (5) species (Table 3).Species diversity determined using the Shannon index (H) ranged from 0.84 (Korhogo) to 1.11 (Boundiali) (Table 3).As the H index is greater than 0, this indicates a high level of species diversity.Piélou's equitability index (J) varied from 0.61 in Korhogo to 0.90 in Bouaké (Table 3).This variation, which differs from the 0 value, shows that the distribution of species in the different localities sampled appears to be identical.As for Jaccard index (I), values varied

Annual variation in overall species composition of leafhopper species
The proportion of species identified varied from one season to the next (Table 6).In 2021, the majority of species identified were Jacobiasca lybica

Geographic variation in overall species composition
The proportions of different species varied among localities (Table 7).Two species were found to be in the majority: Amrasca biguttula (Shiraki, 1913) and Jacobiasca lybica (Bergevin & Zanon, 1922).Their proportions varied respectively from 45 to 50% and from 40 to 50% according to the localities.However, Amrasca biguttula (Shiraki, 1913) was the most abundant.

Dynamics of geographic variation in species composition between 2021 and 2022
From 2021 to 2022, leafhopper species encountered on cotton varied by locality and crop year (Table 8).However, in both crop years, three main species were identified.In 2021, Jacobiasca lybica (Bergevin & Zanon, 1922) was the most dominant species with a proportion ranging from 73.3 to 96.7%.However, in 2022 it was reduced to a proportion ranging from 0 to 6.7% depending on the locality, making it the second most dominant species during this season.Empoasca papayae (Oman, 1937) increased from 0 to 23.3% in 2021, then decreased to 0% in 2022, depending on locality, making it an extinct species in this season.Empoasca facialis (Jacobi, 1912) was the minority species in 2021, with a proportion ranging from 0 to 16.7% and decreased in 2022 to a proportion ranging from 0 to 10% depending on the locality.It remained the minority species.However, this survey identified it in only two localities, including Bouaké and Boundiali.Amrasca biguttula (Shiraki, 1913) went from being absent from the samples collected in 2021 (0%) to the most abundant species in 2022, with a proportion ranging from 90 to 100%, depending on the locality.

Geographic variation in infestation levels
The levels of leafhopper damage varied among localities in the two crop years (Fig. 2).In the 2021 season, attacks were significant in the northern part of the cotton basin.However, they were particularly pronounced in the northeast of the cotton basin.During the 2022 season, attacks were widespread throughout the cotton basin.However, the highest infestations were observed in the northern part of the cotton basin, with higher levels in the northeastern part.

Seasonal and annual variation in infestation levels
Seasonal incidence of leafhoppers varied in both crop years (Fig. 3).In both years, attacks increased throughout the cropping season, but with higher levels in 2022 (more than 30 adult jassids/30 plants).Leafhopper attacks varied by year from 2012 to 2022 (Fig. 4).From 2012 to 2021, average attack levels ranged from 2.3 to 6.04 plants attacked per 30 plants observed.These levels were relatively low compared to those  www.nature.com/scientificreports/observed in 2022.Indeed, the average level of attacks was more than nine attacked plants per 30 plants observed in 2022.This level had never been reached in the previous 10 years (Fig. 4).

Geographical variation in rainfall
Cumulative rainfall from June to July varied downwards from 2021 to 2022 (Fig. 5).During the 2021 season, cumulative rainfall varied between 290 and 530 mm, following the North-East/South-West gradient.In 2022, however, the cumulative rainfall varied between 260 and 410 mm, a drop ranging from 30 to 120 mm.In addition   The two main species identified in this study were also identified in 2022 by 22 on okra and eggplant in similar proportions.These species therefore parasitize several host plants, notably eggplant and okra 22,23 .Amrasca biguttula (Shiraki, 1913) was first reported in Africa in Ghana on okra as a minor pest 24 .Climatic conditions could favour its expansion in the sub-region, as there are no geographical boundaries for these pests.
Damage from jassids also varied from one locality to another.However, damage was most pronounced in the north-east of the cotton-growing basin.The greatest damage was recorded in this part of the cotton-growing zone.They also varied from 2012 to 2022.However, the greatest damage was observed during the 2022 season.This could be due to the diversity observed, as well as lower rainfall and higher temperatures.Indeed 25,26 , have shown that low rainfall combined with high temperatures and low relative humidity contribute to the healthy development of jassids.Climatic factors therefore play an essential role in the composition and outbreak of jassid populations and in the manifestation of damage 1,26 .These conditions could create a favorable environment for the development of new species not previously present in the cotton basin pest complex.Indeed, the rise and fall in temperatures in July play a role in jassid outbreaks.
The appearance of these species, new to the Ivorian parasite facies, would explain the emergence of jassids in recent years.It would also explain the inadequacy of control programs observed in the field.

Conclusion
The present study revealed a variation in the diversity of jassid species.Three (3) to five (5) jassid species were identified in the cotton basin.This abundance follows a geographical and annual dynamic according to the crop year, resulting in the appearance of more invasive and virulent species with an adaptability that is favoured by climatic factors.
The most abundant species observed during this study period were Jacobiasca lybica (Bergevin & Zanon, 1922), identified during the 2021 campaign, and Amrasca biguttula (Shiraki, 1913), identified during the 2022 campaign.Both appear to be invasive.However, Amrasca biguttula (Shiraki, 1913) appears to be the most invasive and virulent.Its appearance has led to the virtual extinction of certain species in certain localities.It has thus disrupted the pest complex in cotton-growing areas, with damage levels exceeding the economic impact threshold throughout the cotton-growing basin.It has had a significant impact on cotton growing.
It is planned to continue identifying species using molecular tools in order to extend knowledge of any species that could not be identified using conventional taxonomy.Bio-ecological studies on the life cycle of Amrasca biguttula (Shiraki, 1913), its natural enemies and the susceptibility of cotton varieties and vegetable crops to jassids are also envisaged, in order to contribute to the implementation of an ecological control strategy.

Figure 3 .
Figure 3. Seasonal variation in adult jassid infestation levels observed during the 2021 and 2022 crop year.
2013 Year 2014 Year 2015 Year 2016 Year 2017 Year 2018 Year 2019 Year 2020 Year 2021 Year 2022 Crop year Number of attacked plants /

Figure 4 .
Figure 4. Annual variation in levels of damage caused by jassid attacks observed over the last 10 years in the Ivorian cotton basin.

Table 1 :
. Characteristics of the main active ingredients used in the phytosanitary protection of cotton in Côte d'Ivoire.

Table 3 .
Three 0.50 to 1 (Table4).Although the I-indices varied from one site to another, values close to 1 showed that the composition of the species found at the different sites was almost identical.The dominance index or simpson index (D) varied from species to species.It was 0.90, 0.90, 0.83, 0.78 and 0.76 respectively for Amrasca biguttula(Shiraki, 1913), Jacobiasca lybica(Bergevin & Zanon, 1922), Empoasca papayae (Oman, 1937), Cicadellidae sp and Empoasca facialis(Jacobi, 1912).The first two species had the highest probability of occurrence.They reveal their abundance in the samples collected.

Table 5 .
Relative proportion of the main leafhopper species identified on cotton in Ivory Coast.

Table 6 .
Annual variation in the overall specific composition of leafhoppers during the 2021 and 2022 cropping seasons in Ivory Coast.

Table 7 .
Geographical variation in the proportions of leafhopper species identified in the 10 localities sampled in Ivory Coast.

Table 8 .
Geographical dynamics of leafhopper species proportions during the 2021 and 2022 cropping seasons in 10 localities in Ivory Coast.