Complete genome sequences of two new virus isolates associated with cotton blue disease resistance breaking in Brazil

Genetic Diversity of Cotton Leafroll Dwarf Virus from the Southwestern United States and Its Implications for the Multi-Introduction Event Hypothesis and Future Evolution

Cotton leafroll dwarf virus (CLRDV) is a viral agent recently identified in the United States in 2017 in Alabama. Since its identification, CLRDV has spread to every cotton-growing state east of New Mexico. Oklahoma, Kansas, and Texas make up the westernmost border of reported CLRDV incidence, making monitoring of these states vital for proper control. Additionally, as the virus evolves, mutations that alter symptomology, such as mutations in the F-box-like motif in ORF0/P0, may occur and need to be monitored thoroughly during the growing seasons. Using high-throughput sequencing and PCR-derived Sanger sequencing, 4 CLRDV genomes and 21 P0 gene isolates were sequenced from Oklahoma, Kansas, and Texas from 2019 to 2021 to determine the genetic diversity among CLRDV isolates. Phylogenetic analyses of the complete genomes revealed seven clades, whereas ORF0 gene analyses resulted in large polytomic clusters. BEAST analyses of the 114 total P0 sequences from GenBank, downloaded before 2024, revealed a lower mean substitution rate than previously reported as well as an earlier root year (1914). In addition, using all available CLRDV genome sequences, 11 likely recombination events were determined. Examination of the P0 amino acid sequences revealed 13 mutations unique to the isolates collected in this study. Based on the phylogenetic and amino acid analyses, the CLRDV isolates from Texas (TX clade) may represent evidence for the multi-introduction event hypothesis into the United States. Additionally, based on our analyses in this study, we propose the Asian CLRDV isolates should be constituted as a potentially separate strain of CLRDV.

Genetic Variability Among U.S.-Sentinel Cotton Plot Cotton Leafroll Dwarf Virus and Globally Available Reference Isolates Based on ORF0 Diversity

The aphid-transmitted polerovirus, cotton leafroll dwarf virus (CLRDV), first characterized from symptomatic cotton plants in South America, has been identified in commercial cotton plantings in the United States. Here, the CLRDV intraspecific diversity was investigated by comparative sequence analysis of the most divergent CLRDV coding region, ORF0/P0. Bayesian analysis of ORF0 sequences for U.S. and reference populations resolved three well-supported sister clades comprising one U.S. and two South American lineages. Principal component analysis (PCA) identified seven statistically supported intraspecific populations. The Bayesian phylogeny and PCA dendrogram-inferred relationships were congruent. Population analysis of ORF0 sequences indicated most lineages have evolved under negative selection, albeit certain sites/isolates evolved under positive selection. Both U.S. and South American isolates exhibited extensive ORF0 diversity. At least two U.S. invasion foci were associated with their founder populations in Alabama-Georgia and eastern Texas. The Alabama-Georgia founder is implicated as the source of recent widespread expansion and establishment of secondary disease foci throughout the southeastern-central United States. Based on the geographically restricted distribution, spread of another extant Texas population appeared impeded by a population bottleneck. Extant CLRDV isolates represent several putative introductions potentially associated with catastrophic weather events dispersing viruliferous cotton aphids of unknown origin(s).

Aphid gene expression following polerovirus acquisition is host species dependent

Upon acquisition of persistent circulative viruses such as poleroviruses, the virus particles transcytose through membrane barriers of aphids at the midgut and salivary glands via hemolymph. Such intricate interactions can influence aphid behavior and fitness and induce associated gene expression in viruliferous aphids. Differential gene expression can be evaluated by omics approaches such as transcriptomics. Previously conducted aphid transcriptome studies used only one host species as the source of virus inoculum. Viruses typically have alternate hosts. Hence, it is not clear how alternate hosts infected with the same virus isolate alter gene expression in viruliferous vectors. To address the question, this study conducted a transcriptome analysis of viruliferous aphids that acquired the virus from different host species. A polerovirus, cotton leafroll dwarf virus (CLRDV), which induced gene expression in the cotton aphid, Aphis gossypii Glover, was assessed using four alternate hosts, viz., cotton, hibiscus, okra, and prickly sida. Among a total of 2,942 differentially expressed genes (DEGs), 750, 310, 1,193, and 689 genes were identified in A. gossypii that acquired CLRDV from infected cotton, hibiscus, okra, and prickly sida, respectively, compared with non-viruliferous aphids that developed on non-infected hosts. A higher proportion of aphid genes were overexpressed than underexpressed following CLRDV acquisition from cotton, hibiscus, and prickly sida. In contrast, more aphid genes were underexpressed than overexpressed following CLRDV acquisition from okra plants. Only four common DEGs (heat shock protein, juvenile hormone acid O-methyltransferase, and two unannotated genes) were identified among viruliferous aphids from four alternate hosts. Gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations indicated that the acquisition of CLRDV induced DEGs in aphids associated with virus infection, signal transduction, immune systems, and fitness. However, these induced changes were not consistent across four alternate hosts. These data indicate that alternate hosts could differentially influence gene expression in aphids and presumably aphid behavior and fitness despite being infected with the same virus isolate.

Analysis of Cotton Leafroll Dwarf Virus P0 Gene Sequences from South Carolina Reveals Low Variability Among Isolates

Cotton leafroll dwarf virus (CLRDV) is emerging across the major cotton-producing states of the southern United States. Because it was detected in nearly all cotton-producing states within a few years of its initial detection in the United States, the spread of the virus has apparently occurred rapidly. In this study spanning three growing seasons in South Carolina, we collected CLRDV isolates from symptomatic and asymptomatic cotton plants in 10 counties. The genomic region encoding P0, the viral suppressor of RNA silencing, was sequenced and compared among CLRDV isolates. Low variability among CLRDV P0 sequences from South Carolina isolates with similarities to other United States isolates was revealed by amino acid sequence alignment and phylogenetic analysis. Low variability among South Carolina isolates was also confirmed by sequencing a subset of eight near-complete genomes of CLRDV isolates. Although sequence variability was low among South Carolina isolates, this data should be taken in the context of all United States isolates, for which diversity may be higher than initially expected. Sequences gathered in this study add to the body of knowledge on CLRDV diversity in the United States.

Cotton leafroll dwarf disease: An enigmatic viral disease in cotton

TAXONOMY

Cotton leafroll dwarf virus (CLRDV) is a member of the genus Polerovirus, family Solemoviridae. Geographical Distribution: CLRDV is present in most cotton-producing regions worldwide, prominently in North and South America.

PHYSICAL PROPERTIES

The virion is a nonenveloped icosahedron with T = 3 icosahedral lattice symmetry that has a diameter of 26-34 nm and comprises 180 molecules of the capsid protein. The CsCl buoyant density of the virion is 1.39-1.42 g/cm3 and S20w is 115-127S. Genome: CLRDV shares genomic features with other poleroviruses; its genome consists of monopartite, single-stranded, positive-sense RNA, is approximately 5.7-5.8 kb in length, and is composed of seven open reading frames (ORFs) with an intergenic region between ORF2 and ORF3a.

TRANSMISSION

CLRDV is transmitted efficiently by the cotton aphid (Aphis gossypii Glover) in a circulative and nonpropagative manner. Host: CLRDV has a limited host range. Cotton is the primary host, and it has also been detected in different weeds in and around commercial cotton fields in Georgia, USA.

SYMPTOMS

Cotton plants infected early in the growth stage exhibit reddening or bronzing of foliage, maroon stems and petioles, and drooping. Plants infected in later growth stages exhibit intense green foliage with leaf rugosity, moderate to severe stunting, shortened internodes, and increased boll shedding/abortion, resulting in poor boll retention. These symptoms are variable and are probably influenced by the time of infection, plant growth stage, varieties, soil health, and geographical location. CLRDV is also often detected in symptomless plants.

CONTROL

Vector management with the application of chemical insecticides is ineffective. Some host plant varieties grown in South America are resistant, but all varieties grown in the United States are susceptible. Integrated disease management strategies, including weed management and removal of volunteer stalks, could reduce the abundance of virus inoculum in the field.

The Viral Threat in Cotton: How New and Emerging Technologies Accelerate Virus Identification and Virus Resistance Breeding

Cotton (Gossypium spp. L., Malvaceae) is the world's largest source of natural fibers. Virus outbreaks are fast and economically devasting regarding cotton. Identifying new viruses is challenging as virus symptoms usually mimic nutrient deficiency, insect damage, and auxin herbicide injury. Traditional viral identification methods are costly and time-consuming. Developing new resistant cotton lines to face viral threats has been slow until the recent use of molecular virology, genomics, new breeding techniques (NBT), remote sensing, and artificial intelligence (AI). This perspective article demonstrates rapid, sensitive, and cheap technologies to identify viral diseases and propose their use for virus resistance breeding.

Cotton Leafroll Dwarf Virus US Genomes Comprise Divergent Subpopulations and Harbor Extensive Variability

Cotton leafroll dwarf virus (CLRDV) was first reported in the United States (US) in 2017 from cotton plants in Alabama (AL) and has become widespread in cotton-growing states of the southern US. To investigate the genomic variability among CLRDV isolates in the US, complete genomes of the virus were obtained from infected cotton plants displaying mild to severe symptoms from AL, Florida, and Texas. Eight CLRDV genomes were determined, ranging in size from 5865 to 5867 bp, and shared highest nucleotide identity with other CLRDV isolates in the US, at 95.9–98.7%. Open reading frame (ORF) 0, encoding the P0 silencing suppressor, was the most variable gene, sharing 88.5–99.6% and 81.2–89.3% amino acid similarity with CLRDV isolates reported in cotton growing states in the US and in Argentina and Brazil in South America, respectively. Based on Bayesian analysis, the complete CLRDV genomes from cotton in the US formed a monophyletic group comprising three relatively divergent sister clades, whereas CLRDV genotypes from South America clustered as closely related sister-groups, separate from US isolates, patterns reminiscent of phylogeographical structuring. The CLRDV isolates exhibited a complex pattern of recombination, with most breakpoints evident in ORFs 2 and 3, and ORF5. Despite extensive nucleotide diversity among all available CLRDV genomes, purifying selection (dN/dS < 1) was implicated as the primary selective force acting on viral protein evolution.

Genome analysis of cotton leafroll dwarf virus reveals variability in the silencing suppressor protein, genotypes and genomic recombinants in the USA

Cotton leafroll dwarf virus (CLRDV) is an emerging virus in cotton production in Georgia and several other Southeastern states in the USA. To better understand the genetic diversity of the virus population, the near complete genome sequences of six isolates from Georgia and one from Alabama were determined. The isolates sequenced were 5,866 nucleotides with seven open reading frames (ORFs). The isolates from Georgia were >94% identical with other isolates from the USA and South America. In the silencing suppressor protein (P0), at amino acid position 72, the isolates from Georgia and Alabama had a valine (V), similar to resistant-breaking 'atypical' genotypes in South America, while the Texas isolate had isoleucine (I), similar to the more aggressive 'typical' genotypes of CLRDV. At position 120, arginine (R) is unique to Georgia and China isolates, but absent in Alabama, Texas and South American isolates. Ten potential recombinant events were detected in the isolates sequenced. An increased understanding of CLRDV population structure and genetic diversity will help develop management strategies for CLRDV in the USA cotton belt.

Genome Sequence of Cotton Leafroll Dwarf Virus Infecting Cotton in Georgia, USA

Cotton leafroll dwarf disease (CLRDD), caused by the aphid-borne Cotton leafroll dwarf virus (CLRDV; genus, Polerovirus; family, Luteoviridae), has been recently reported from the major cotton-growing regions of the United States. Here, we present the nearly complete genome sequence of a CLRDV isolate from cotton in Georgia.

Characterization of the Complete Genome and P0 Protein for a Previously Unreported Genotype of Cotton Leafroll Dwarf Virus, an Introduced Polerovirus in the United States

Virus-like disease symptoms consisting of leaf cupping, shortened internodes, and overall stunting were observed in commercial cotton fields in Alabama in 2017 to 2018. To determine the complete genome sequence of the suspected causal polerovirus, symptomatic leaf samples were collected in Macon County, Alabama, and subjected to Illumina RNA sequencing. Based on BLASTn analysis, the Illumina contig of 5,771 nt shared the highest nucleotide identity (approximately 95%) with members of the species Cotton leafroll dwarf virus (CLRDV) (genus Polerovirus; family Luteoviridae) from Argentina and Brazil. The full-length viral genome sequence was verified by reverse transcription (RT)-PCR amplification, cloning, and Sanger sequencing. The complete CLRDV genome of 5,865 nt in length shared 94.8 to 95.2% nucleotide identity with six previously reported CLRDV isolates. The genome of the CLRDV isolate amplified from Alabama samples (CLRDV-AL) has seven predicted open reading frames (ORFs). Viral proteins 1 to 5 (P1 to P5) shared 91.9 to 99.5% amino acid identity with the six CLRDV isolates from Argentina and Brazil. However, P0, the suppressor of host gene silencing, shared 82.4 to 88.5% pairwise amino acid identity with the latter CLRDV isolates. Phylogenetic analysis of the seven full-length CLRDV genomes resolved three sister clades: CLRDV-AL, CLRDV-typical, and CLRDV-atypical, respectively. Three recombination events were detected by the recombination detection program among the seven CLRDV isolates with breakpoints occurring along the genome. Pairwise nucleotide identity comparisons of ORF0 sequences for the three CLRDV-AL field isolates indicated that they were >99% identical, suggesting that this previously unknown CLRDV genotype represents a single introduction to Alabama.

Tomato Yellow Leaf Curl Virus-Resistant and -Susceptible Tomato Genotypes Similarly Impact the Virus Population Genetics

Tomato yellow leaf curl virus is a species in the genus Begomovirus and family Geminiviridae. Tomato yellow leaf curl virus (TYLCV) infection induces severe symptoms on tomato plants and causes serious yield losses worldwide. TYLCV is persistently transmitted by the sweetpotato whitefly, Bemisia tabaci (Gennadius). Cultivars and hybrids with a single or few genes conferring resistance against TYLCV are often planted to mitigate TYLCV-induced losses. These resistant genotypes (cultivars or hybrids) are not immune to TYLCV. They typically develop systemic infection, display mild symptoms, and produce more marketable tomatoes than susceptible genotypes under TYLCV pressure. In several pathosystems, extensive use of resistant cultivars with single dominant resistance-conferring gene has led to intense selection pressure on the virus, development of highly virulent strains, and resistance breakdown. This study assessed differences in TYLCV genomes isolated from susceptible and resistant genotypes in Florida and Georgia. Phylogenetic analyses indicated that Florida and Georgia isolates were distinct from each other. Population genetics analyses with genomes field-collected from resistant and susceptible genotypes from Florida and/or Georgia provided no evidence of a genetic structure between the resistant and susceptible genotypes. No codons in TYLCV genomes from TYLCV-resistant or susceptible genotypes were under positive selection, suggesting that highly virulent or resistance-breaking TYLCV strains might not be common in tomato farmscapes in Florida and Georgia. With TYLCV-resistant genotypes usage increasing recently and multiple tomato crops being planted during a calendar year, host resistance-induced selection pressure on the virus remains a critical issue. To address the same, a greenhouse selection experiment with one TYLCV-resistant and susceptible genotype was conducted. Each genotype was challenged with TYLCV through whitefly-mediated transmission serially 10 times (T₁-T₁₀). Population genetics parameters at the genome level were assessed at T₁, T₅, and T₁₀. Results indicated that genomes from resistant and susceptible genotypes did not differentiate with increasing transmission number, no specific mutations were repeatedly observed, and no positive selection was detected. These results reiterate that resistance in tomato might not be exerting selection pressure against TYLCV to facilitate development of resistance-breaking strains. TYLCV populations rather seem to be shaped by purifying selection and/or population expansion.

Identification of a New Cotton Disease Caused by an Atypical Cotton Leafroll Dwarf Virus in Argentina

An outbreak of a new disease occurred in cotton (Gossypium hirsutum) fields in northwest Argentina starting in the 2009-10 growing season and is still spreading steadily. The characteristic symptoms of the disease included slight leaf rolling and a bushy phenotype in the upper part of the plant. In this study, we determined the complete nucleotide sequences of two independent virus genomes isolated from cotton blue disease (CBD)-resistant and -susceptible cotton varieties. This virus genome comprised 5,866 nucleotides with an organization similar to that of the genus Polerovirus and was closely related to cotton leafroll dwarf virus, with protein identity ranging from 88 to 98%. The virus was subsequently transmitted to a CBD-resistant cotton variety using Aphis gossypii and symptoms were successfully reproduced. To study the persistence of the virus, we analyzed symptomatic plants from CBD-resistant varieties from different cotton-growing fields between 2013 and 2015 and showed the presence of the same virus strain. In addition, a constructed full-length infectious cDNA clone from the virus caused disease symptoms in systemic leaves of CBD-resistant cotton plants. Altogether, the new leafroll disease in CBD-resistant cotton plants is caused by an atypical cotton leafroll dwarf virus.

Function and diversity of P0 proteins among cotton leafroll dwarf virus isolates

Background

The RNA silencing pathway is an important anti-viral defense mechanism in plants. As a counter defense, some members of the viral family Luteoviridae are able to evade host immunity by encoding the P0 RNA silencing suppressor protein. Here we explored the functional diversity of P0 proteins among eight cotton leafroll dwarf virus (CLRDV) isolates, a virus associated with a worldwide cotton disease known as cotton blue disease (CBD).

Methods

CLRDV-infected cotton plants of different varieties were collected from five growing fields in Brazil and their P0 sequences compared to three previously obtained isolates. P0’s silencing suppression activities were scored based on transient expression experiments in Nicotiana benthamiana leaves.

Results

High sequence diversity was observed among CLRDV P0 proteins, indicating that some isolates found in cotton varieties formerly resistant to CLRDV should be regarded as new genotypes within the species. All tested proteins were able to suppress local and systemic silencing, but with significantly variable degrees. All P0 proteins were able to mediate the decay of ARGONAUTE proteins, a key component of the RNA silencing machinery.

Conclusions

The sequence diversity observed in CLRDV P0s is also reflected in their silencing suppression capabilities. However, the strength of local and systemic silencing suppression was not correlated for some proteins.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-015-0356-7) contains supplementary material, which is available to authorized users.