Complete Genome Sequence of the Largest Known Flavi-Like Virus, Diaphorina citri flavi-like virus, a Novel Virus of the Asian Citrus Psyllid, Diaphorina citri

An insect virus differentially alters gene expression among life stages of an insect vector and enhances bacterial phytopathogen transmission

Diaphorina citri transmits Candidatus Liberibacter asiaticus (CLas) between citrus plants which causes the expression of huanglongbing disease in citrus. D. citri flavi-like virus (DcFLV) co-occurs intracellularly with CLas in D. citri populations in the field. However, the impact(s) of DcFLV presence on the insect vector and its interaction with the CLas phytopathogen remain unclear. We compared CLas acquisition and transmission efficiencies as well as transcriptomic expression between viruliferous and non-viruliferous psyllids at multiple life stages. Viruliferous nymphs acquired higher titers of CLas than non-viruliferous nymphs, whereas viruliferous adults acquired less CLas than those without virus. The presence of DcFLV increased the transmission of CLas by both nymphs and adults. Furthermore, RNA-seq and functional gene expression analyses revealed that endoplasmic reticulum stress-, autophagy-, and defense-related genes were significantly upregulated in viruliferous adult psyllids, whereas most of these genes were downregulated in viruliferous nymphs. Our work demonstrates that DcFLV differentially modulates various cellular and physiological functions in D. citri in a life stage-dependent manner and promotes the acquisition of CLas at the nymphal stage and transmission of the pathogen at the adult stage of the vector. Collectively, our results suggest that D. citri vectors with DcFLV exhibit greater pathogen transmission efficiency than those without virus.

IMPORTANCE

Huanglongbing (HLB), caused by fastidious bacteria from three Candidatus Liberibacter species, is the most damaging disease impacting the citrus industry worldwide. Spread by the Asian citrus psyllid (Diaphorina citri) in Asia and the Americas, HLB causes substantial financial losses, and has reduced citrus production in Florida by more than 90%. Although there are ongoing efforts to limit spread of the disease, effective HLB management remains elusive. Suppressing vector populations and decreasing CLas transmission are the two strategies that need to be urgently improved. Recently, a D. citri flavi-like virus (DcFLV) was characterized within its D. citri host, and it co-occurs intracellularly with CLas in psyllid populations. Here, we show that viruliferous nymphs exhibit higher CLas acquisition than non-viruliferous nymphs. Furthermore, both viruliferous adults and nymphs exhibit increased CLas transmission efficiency. We suggest the possibility of manipulating DcFLV in D. citri populations to reduce CLas transmission for HLB disease management.

Integrated analysis of protein sequence and structure redefines viral diversity and the taxonomy of the Flaviviridae

The Flaviviridae are a family of non-segmented positive-sense enveloped RNA viruses containing significant pathogens including hepatitis C virus and yellow fever virus. Recent large-scale metagenomic surveys have identified many diverse RNA viruses related to classical orthoflaviviruses and pestiviruses but quite different genome lengths and configurations, and with a hugely expanded host range that spans multiple animal phyla, including molluscs, cnidarians and stramenopiles,, and plants. Grouping of RNA-directed RNA polymerase (RdRP) hallmark gene sequences of flavivirus and 'flavi-like' viruses into four divergent clades and multiple lineages within them was congruent with helicase gene phylogeny, PPHMM profile comparisons, and comparison of RdRP protein structure predicted by AlphFold2. These results support their classification into the established order, Amarillovirales, in three families (Flaviviridae, Pestiviridae, and Hepaciviridae), and 14 genera. This taxonomic framework informed by RdRP hallmark gene evolutionary relationships provides a stable reference from which major genome re-organisational events can be understood.

A ~40-kb flavi-like virus does not encode a known error-correcting mechanism

It is commonly held that there is a fundamental relationship between genome size and error rate, manifest as a notional "error threshold" that sets an upper limit on genome sizes. The genome sizes of RNA viruses, which have intrinsically high mutation rates due to a lack of mechanisms for error correction, must therefore be small to avoid accumulating an excessive number of deleterious mutations that will ultimately lead to population extinction. The proposed exceptions to this evolutionary rule are RNA viruses from the order Nidovirales (such as coronaviruses) that encode error-correcting exonucleases, enabling them to reach genome lengths greater than 40 kb. The recent discovery of large-genome flavi-like viruses (Flaviviridae), which comprise genomes up to 27 kb in length yet seemingly do not encode exonuclease domains, has led to the proposal that a proofreading mechanism is required to facilitate the expansion of nonsegmented RNA virus genomes above 30 kb. Herein, we describe a ~40 kb flavi-like virus identified in a Haliclona sponge metatranscriptome that does not encode a known exonuclease. Structural analysis revealed that this virus may have instead captured cellular domains associated with nucleic acid metabolism that have not been previously found in RNA viruses. Phylogenetic inference placed this virus as a divergent pesti-like lineage, such that we have provisionally termed it "Maximus pesti-like virus." This virus represents an instance of a flavi-like virus achieving a genome size comparable to that of the Nidovirales and demonstrates that RNA viruses have evolved multiple solutions to overcome the error threshold.

Discovery of Brassica Yellows Virus and Porcine Reproductive and Respiratory Syndrome Virus in Diaphorina citri and Changes in Virome Due to Infection with 'Ca. L. asiaticus'

Detection of new viruses or new virus hosts is essential for the protection of economically important agroecosystems and human health. Increasingly, metatranscriptomic data are being used to facilitate this process. Such data were obtained from adult Asian citrus psyllids (ACP) (Diaphorina citri Kuwayama) that fed solely on mandarin (Citrus ×aurantium L.) plants grafted with buds infected with 'Candidatus Liberibacter asiaticus' (CLas), a phloem-limited bacterium associated with the severe Asian variant of huanglongbing (HLB), the most destructive disease of citrus. Brassica yellows virus (BrYV), the causative agent of yellowing or leafroll symptoms in brassicaceous plants, and its associated RNA (named as BrYVaRNA) were detected in ACP. In addition, the porcine reproductive and respiratory syndrome virus (PRRSV), which affects pigs and is economically important to pig production, was also found in ACP. These viruses were not detected in insects feeding on plants grafted with CLas-free buds. Changes in the concentrations of insect-specific viruses within the psyllid were caused by coinfection with CLas. IMPORTANCE The cross transmission of pathogenic viruses between different farming systems or plant communities is a major threat to plants and animals and, potentially, human health. The use of metagenomics is an effective approach to discover viruses and vectors. Here, we collected buds from the CLas-infected and CLas-free mandarin (Citrus ×aurantium L. [Rutaceae: Aurantioideae: Aurantieae]) trees from a commercial orchard and grafted them onto CLas-free mandarin plants under laboratory conditions. Through metatranscriptome sequencing, we first identified the Asian citrus psyllids feeding on plants grafted with CLas-infected buds carried the plant pathogen, brassica yellows virus and its associated RNA, and the swine pathogen, porcine reproductive and respiratory syndrome virus. These discoveries indicate that both viruses can be transmitted by grafting and acquired by ACP from CLas+ mandarin seedlings.

High-quality, chromosome-scale genome assemblies: comparisons of three Diaphorina citri (Asian citrus psyllid) geographic populations

The Asian citrus psyllid, Diaphorina citri, is the insect vector of the causal agent of huanglongbing (HLB), a devastating bacterial disease of commercial citrus. Presently, few genomic resources exist for D. citri. In this study, we utilized PacBio HiFi and chromatin confirmation contact (Hi-C) sequencing to sequence, assemble, and compare three high-quality, chromosome-scale genome assemblies of D. citri collected from California, Taiwan, and Uruguay. Our assemblies had final sizes of 282.67 Mb (California), 282.89 Mb (Taiwan), and 266.67 Mb (Uruguay) assembled into 13 pseudomolecules—a reduction in assembly size of 41–45% compared with previous assemblies which we validated using flow cytometry. We identified the X chromosome in D. citri and annotated each assembly for repetitive elements, protein-coding genes, transfer RNAs, ribosomal RNAs, piwi-interacting RNA clusters, and endogenous viral elements. Between 19,083 and 20,357 protein-coding genes were predicted. Repetitive DNA accounts for 36.87–38.26% of each assembly. Comparative analyses and mitochondrial haplotype networks suggest that Taiwan and Uruguay D. citri are more closely related, while California D. citri are closely related to Florida D. citri. These high-quality, chromosome-scale assemblies provide new genomic resources to researchers to further D. citri and HLB research.

Analysis of Citrus Tristeza Virus Incidences within Asian Citrus Psyllid (Diaphorina citri) Populations in Florida via High-Throughput Sequencing

The destructive citrus disease, Huanglongbing (HLB) or citrus greening, continues to devastate Florida's citrus industry. A hemipteran insect, the Asian citrus psyllid (ACP), disperses Candidatus Liberibacter asiaticus, one of the putative bacterial pathogens of HLB. This study builds upon ongoing research utilizing high-throughput sequencing to analyze the virome of ACP populations collected from citrus groves throughout Florida. Following the widespread detection of sequences aligning to the genome of citrus tristeza virus (CTV) across consecutive years in the Florida ACP virome, we continued to detect a pervasive amount of CTV in Florida ACPs during subsequent years. Simultaneously, we also detected mixed infections of CTV strains in pooled ACPs from different Florida regions. Predating the HLB epidemic, CTV has been present in Florida for many years and our results confirm its widespread and diverse persistence in Florida citrus groves through a unique lens, the ACP. CTV presence in the ACP likely results from feeding on CTV-infected citrus trees in Florida citrus groves, which may help to understand an overlapping presence of CTV and HLB, both endemic citrus pathosystems in the state, and their role in future integrated pest management strategies.

Diaphorina citri flavi-like virus localization, transmission, and association with Candidatus Liberibacter asiaticus in its psyllid host

Huanglongbing is caused by Candidatus Liberibacter asiaticus (CLas) and transmitted by Diaphorina citri. D. citri harbors various insect-specific viruses, including the Diaphorina citri flavi-like virus (DcFLV). The distribution and biological role of DcFLV in its host and the relationship with CLas are unknown. DcFLV was found in various organs of D. citri, including the midgut and salivary glands, where it co-localized with CLas. CLas-infected nymphs had the highest DcFLV titers compared to the infected adults and CLas-free adults and nymphs. DcFLV was vertically transmitted to offspring from female D. citri and was temporarily detected in Citrus macrophylla and grapefruit leaves from greenhouse and field. The incidences of DcFLV and CLas were positively correlated in field-collected D. citri samples, suggesting that DcFLV might be associated with CLas in the vector. These results provide new insights on the interactions between DcFLV, the D. citri, and CLas.

Complexity and local specificity of the virome associated with tospovirus-transmitting thrips species

Frankliniella occidentalis (western flower thrips=WFT) and Thrips tabaci (onion thrips=OT) are insect species that greatly impact horticultural crops through direct damage and their efficient vectoring of tomato spotted wilt virus and iris yellow spot virus. In this study we collected thrips of these species from 12 field populations in various regions in Italy. We also included one field population of Neohydatothrips variabilis (soybean thrips=ST) from the U.S.A. Total RNAseq from high-throughput sequencing (HTS) was used to assemble the virome and then we assigned putative viral contigs to each thrips sample by qRT-PCR. Excluding plant and fungal viruses, we were able to identify 61 viral segments, corresponding to 41 viruses: 14 were assigned to WFT, 17 to OT, one from ST and 9 viruses could not be assigned to any species based on our stringent criteria. All these viruses are putative representative of new species (with only the exception of a sobemo-like virus that is 100% identical to a virus recently characterized in ST) and some belong to new higher-ranking taxa. These additions to the viral phylogeny suggest previously undescribed evolutionary niches. Most of the Baltimore's classes of RNA viruses were present (positive- and minus- strand and dsRNA viruses), but only one DNA virus was identified in our collection. Repeated sampling in a subset of locations in 2019 and 2020 and further virus characterization in a subset of four thrips populations maintained in laboratory allowed us to provide evidence of a locally persistent thrips core virome that characterizes each population. IMPORTANCE Harnessing the insect microbiome can result in new approaches to contain their populations or the damage they cause vectoring viruses of medical, veterinary, or agricultural importance. Persistent insect viruses are a neglected component of their microbiota. Here for the first time, we characterize the virome associated with the two model systems for tospovirus-transmitting thrips species, of utmost importance for the direct and indirect damage they cause to a number of different crops. The thrips virome here characterized includes several novel viruses, that in some cases reveal previously undescribed clades. More importantly, some of the viruses we describe are part of a core virome that is specific and consistently present in distinct geographical locations monitored over the years, hinting at a possible mutualistic symbiotic relationship with their host.

The enzymes for genome size increase and maintenance of large (+)RNA viruses

With sizes <50 kb, viral RNA genomes are at the crossroads of genetic, biophysical, and biochemical stability in their host cell. Here, we analyze the enzymatic assets accompanying large RNA genome viruses, mostly based on recent scientific advances in Coronaviridae. We argue that, in addition to the presence of an RNA exonuclease (ExoN), two markers for the large size of viral RNA genomes are (i) the presence of one or more RNA methyltransferases (MTases) and (ii) a specific architecture of the RNA-dependent RNA polymerase active site. We propose that RNA genome expansion and maintenance are driven by an evolutionary ménage-à-trois made of fast and processive RNA polymerases, RNA repair ExoNs, and RNA MTases that relates to the transition between RNA- to DNA-based life.

Endogenous Viral Element-Derived Piwi-Interacting RNAs (piRNAs) Are Not Required for Production of Ping-Pong-Dependent piRNAs from Diaphorina citri Densovirus

Piwi-interacting RNAs (piRNAs) are a class of small RNAs primarily responsible for silencing transposons in the animal germ line. The ping-pong cycle, the posttranscriptional silencing branch of the piRNA pathway, relies on piRNAs produced from endogenous transposon remnants to direct cleavage of transposon RNA via association with Piwi-family Argonaute proteins. In some mosquito species and mosquito-derived cell lines expressing a functionally expanded group of Piwi-family Argonaute proteins, both RNA and DNA viruses are targeted by piRNAs in a manner thought to involve direct processing of exogenous viral RNA into piRNAs. Whether viruses are targeted by piRNAs in nonmosquito species is unknown. Partial integrations of DNA and nonretroviral RNA virus genomes, termed endogenous viral elements (EVEs), are abundant in arthropod genomes and often produce piRNAs that are speculated to target cognate viruses through the ping-pong cycle. Here, we describe a Diaphorina citri densovirus (DcDV)-derived EVE in the genome of Diaphorina citri We found that this EVE gives rise to DcDV-specific primary piRNAs and is unevenly distributed among D. citri populations. Unexpectedly, we found that DcDV is targeted by ping-pong-dependent virus-derived piRNAs (vpiRNAs) in D. citri lacking the DcDV-derived EVE, while four naturally infecting RNA viruses of D. citri are not targeted by vpiRNAs. Furthermore, a recombinant Cricket paralysis virus containing a portion of the DcDV genome corresponding to the DcDV-derived EVE was not targeted by vpiRNAs during infection in D. citri harboring the EVE. These results demonstrate that viruses can be targeted by piRNAs in a nonmosquito species independently of endogenous piRNAs.IMPORTANCE Small RNAs serve as specificity determinants of antiviral responses in insects. Piwi-interacting RNAs (piRNAs) are a class of small RNAs found in animals, and their primary role is to direct antitransposon responses. These responses require endogenous piRNAs complementary to transposon RNA. Additionally, piRNAs have been shown to target RNA and DNA viruses in some mosquito species. In contrast to transposons, targeting of viruses by the piRNA pathway in these mosquito species does not require endogenous piRNAs. Here, we show that piRNAs target a DNA virus, but not RNA viruses, in an agricultural insect pest. We found that targeting of this DNA virus did not require endogenous piRNAs and that endogenous piRNAs did not mediate targeting of an RNA virus with which they shared complementary sequence. Our results highlight differences between mosquitoes and our experimental system and raise the possibility that DNA viruses may be targeted by piRNAs in other species.

Diaphorina citri reovirus is most closely related to fijiviruses

The Asian citrus psyllid, Diaphorina citri Kuwayama, is an important insect vector of Candidatus Liberibacter asiaticus, the causal agent of Huanglongbing, which is the most destructive disease of citrus worldwide. Sequences for putative Diaphorina citri reovirus (DcRV) were identified from some worldwide populations of D. citri. Here, field surveys indicated that the virus was common in D. citri populations from Hawaii and Fuzhou of PR China. Electron microscopy showed that DcRV virions possessed a typical reovirus-like morphology. The U. S. and Chinese DcRV isolates both showed 10 segments of double-stranded RNA sharing >96% nucleotide sequence identity, and encoding 11 deduced proteins. All genome segments contained conserved 5' and 3' terminal nucleotide sequences and inverted repeats that are hallmarks of reovirus sequence. Phylogenetic analysis showed that DcRV may be considered a new species of the genus Fijivirus sharing a most recent common ancestor with the insect-specific fijivirus Nilaparvata lugens reovirus.

Proteomics Computational Analyses Suggest That the Envelope Glycoproteins of Segmented Jingmen Flavi-Like Viruses are Class II Viral Fusion Proteins (b-Penetrenes) with Mucin-Like Domains

Jingmen viruses are newly described segmented flavi-like viruses that have a worldwide distribution in ticks and have been associated with febrile illnesses in humans. Computational analyses were used to predict that Jingmen flavi-like virus glycoproteins have structural features of class II viral fusion proteins, including an ectodomain consisting of beta-sheets and short alpha-helices, a fusion peptide with interfacial hydrophobicity and a three-domain architecture. Jingmen flavi-like virus glycoproteins have a sequence enriched in serine, threonine, and proline at the amino terminus, which is a feature of mucin-like domains. Several of the serines and threonines are predicted be modified by the addition of O-linked glycans. Some of the glycoproteins are predicted to have an additional mucin-like domain located prior to the transmembrane anchor, whereas others are predicted to have a stem consisting of two alpha-helices. The flavivirus envelope protein and Jingmen flavi-virus glycoproteins may have diverged from a common class II precursor glycoprotein with a mucin-like domain or domains acquired after divergence.

The Detection and Surveillance of Asian Citrus Psyllid (Diaphorina citri)-Associated Viruses in Florida Citrus Groves

The plant pathogenic bacterium Candidatus Liberibacter asiaticus (CLas), the causal agent of the citrus disease Huanglongbing (HLB), and its insect vector, the Asian citrus psyllid (ACP; Diaphorina citri), have been devastating the Florida citrus industry. To restore the competitive production presence of Florida in the worldwide citrus market, effective and sustainable control of HLB and the ACP needs to be identified. As alternatives for resistance-inducing insecticides, viruses are currently being considered for biological control of the ACP. To identify possible biological control candidates, we conducted one of the most comprehensive surveys of natural ACP populations in major citrus production regions spanning 21 counties in Florida. By optimizing PCRs and RT-PCRs, we were able to successfully detect and monitor the prevalence of five previously identified ACP-associated RNA and DNA viruses throughout Florida citrus groves, which include: Diaphorina citri-associated C virus (DcACV), Diaphorina citri flavi-like virus (DcFLV), Diaphorina citri densovirus (DcDNV), Diaphorina citri reovirus (DcRV), and Diaphorina citri picorna-like virus (DcPLV). Adult and nymph ACP populations from 21 of Florida's major citrus-producing counties were collected each month during approximately 18 consecutive months. RNA extracts used for these viral screens were also regionally combined and subjected to High Throughput Sequencing (HTS) to reveal a more comprehensive picture of known and unknown viruses in Florida ACP populations. We discovered that DcACV was the most prevalent ACP-associated virus throughout nymph and adult ACP populations in Florida, detected in more than 60% of all samples tested, followed by DcPLV and DcFLV. HTS allowed us to identify a novel ACP-associated reo-like virus and a picorna-like virus. The putative reo-like virus, tentatively named Diaphorina citri cimodo-like virus, was later surveyed and detected back in seasonal adult and nymph ACP samples collected in Florida during this study. HTS generated data also revealed that the most abundant virus in Florida ACP populations was Citrus tristeza virus (CTV), which is not an ACP-associated virus, suggesting persistent presence of CTV infection in citrus throughout Florida groves. Collectively, information obtained from our study may be able to help guide the direction of biotechnological pest control efforts involving a number of viruses that were detected for the first time in Florida ACP populations, including two newly identified ACP-associated viruses.