Global mtDNA genetic structure and hypothesized invasion history of a major pest of citrus, Diaphorina citri (Hemiptera: Liviidae)

Identifying a potentially invasive population in the native range of a species: The enlightenment from the phylogeography of the yellow spotted stink bug, Erthesina fullo (Hemiptera: Pentatomidae)

The yellow spotted stink bug (YSSB), Erthesina fullo (Thunberg, 1783) is an important Asian pest that has recently successfully invaded Europe and an excellent material for research on the initial stage of biological invasion. Here, we reported the native evolutionary history, recent invasion history, and potential invasion threats of YSSB for the first time based on population genetic methods [using double digest restriction-site associated DNA (ddRAD) data and mitochondrial COI and CYTB] and ecological niche modelling. The results showed that four lineages (east, west, southwest, and Hainan Island) were established in the native range with a strong east-west differentiation phylogeographical structure, and the violent climate fluctuation might cause population divergence during the Middle and Upper Pleistocene. In addition, land bridges and monsoon promote dispersal and directional genetic exchanging between island populations and neighboring continental populations. The east lineage (EA) was identified as the source of invasion in Albania. EA had the widest geographical distribution among all other lineages, with a star-like haplotype network with the main haplotype as the core. It also had a rapid population expansion history, indicating that the source lineage might have stronger diffusion ability and adaptability. Our findings provided a significant biological basis for fine tracking of invasive source at the lineage or population level and promote early invasion warning of potential invasive species on a much subtler lineage level.

Risk of spread of the Asian citrus psyllid Diaphorina citri Kuwayama (Hemiptera: Liviidae) in Ghana

The impact of invasive species on biodiversity, food security and economy is increasingly noticeable in various regions of the globe as a consequence of climate change. Yet, there is limited research on how climate change affects the distribution of the invasive Asian citrus psyllid Diaphorina citri Kuwayama (Hemiptera:Liviidae) in Ghana. Using maxnet package to fit the Maxent model in R software, we answered the following questions; (i) what are the main drivers for D. citri distribution, (ii) what are the D. citri-specific habitat requirements and (iii) how well do the risk maps fit with what we know to be correctly based on the available evidence?. We found that temperature seasonality (Bio04), mean temperature of warmest quarter (Bio10), precipitation of driest quarter (Bio17), moderate resolution imaging spectroradiometer land cover and precipitation seasonality (Bio15), were the most important drivers of D. citri distribution. The results follow the known distribution records of the pest with potential expansion of habitat suitability in the future. Because many invasive species, including D. citri, can adapt to the changing climates, our findings can serve as a guide for surveillance, tracking and prevention of D. citri spread in Ghana.

Characterization of the Asian Citrus Psyllid-‘Candidatus Liberibacter Asiaticus’ Pathosystem in Saudi Arabia Reveals Two Predominant CLas Lineages and One Asian Citrus Psyllid Vector Haplotype

In Saudi Arabia (SA), the citrus greening disease is caused by ‘Candidatus Liberibacter asiaticus’ (CLas) transmitted by the Asian citrus psyllid (ACP) Diaphorina citri. The origin and route(s) of the ACP-CLas pathosystem invasion in SA have not been studied. Adult ACP were collected from citrus trees in SA and differentiated by analysis of the mitochondrial cytochrome oxidase I (mtCOI) and nuclear copper transporting protein (atox1) genes. A phylogenetic analysis of the Wolbachia spp. surface protein (wsp) gene was used to identify the ACP-associated Wolbachia spp. A phylogenetic analysis of the atox1 and mtCOI gene sequences revealed one predominant ACP haplotype most closely related to the Indian subcontinent founder populations. The detection and identification of CLas in citrus trees were carried out by polymerase chain reaction (PCR) amplification and sequencing of the 16S rDNA gene. The CLas-integrated prophage genomes were sequenced, annotated, and used to differentiate CLas populations. The ML and ASTRAL trees reconstructed with prophages type 1 and 2 genome sequences, separately and concatenated, resolved two major lineages, CLas-1 and -2. The CLas-1 clade, reported here for the first time, consisted of isolates from SA isolates and Pakistan. The CLas-2 sequences formed two groups, CLas-2-1 and -2-2, previously the ‘Asiatic’ and ‘Floridian’ strains, respectively. Members of CLas-2-1 originated from Southeast Asia, the USA, and other worldwide locations, while CLas-2-2 was identified only in Florida. This study provides the first snapshot into the status of the ACP-CLas pathosystem in SA. In addition, the results provide new insights into the pathosystem coevolution and global invasion histories of two ACP-CLas lineages with a predicted center of origin in South and Southeast Asia, respectively.

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.

Two Liberibacter Proteins Combine to Suppress Critical Innate Immune Defenses in Citrus

We adopted a systems-based approach to determine the role of two Candidatus Liberibacter asiaticus (CLas) proteins, LasP ₂₃₅ and Effector 3, in Huanglongbing (HLB) pathogenesis. While a published work suggests the involvement of these CLas proteins HLB pathogenesis, the exact structure-based mechanism of their action has not been elucidated. We conducted the following experiments to determine the structure-based mechanisms of action. First, we immunoprecipitated the interacting citrus protein partners of LasP ₂₃₅ and Effector 3 from the healthy and CLas-infected Hamlin extracts and identified them by Liquid Chromatography with tandem mass spectrometry (LC-MS/MS). Second, we performed a split green fluorescent protein (GFP) assay in tobacco to validate that the interactions observed in vitro are also retained in planta. The notable in planta citrus targets of LasP ₂₃₅ and Effector 3 include citrus innate immune proteins. Third, in vitro and in planta studies were performed to show that LasP ₂₃₅ and Effector 3 interact with and inhibit the functions of multiple citrus proteins belonging to the innate immune pathways. These inhibitory interactions led to a high level of reactive oxygen species, blocking of bactericidal lipid transfer protein (LTP), and induction of premature programed cell death (PCD), all of which are beneficial to CLas lifecycle and HLB pathogenesis. Finally, we performed molecular dynamics simulations to visualize the interactions of LasP ₂₃₅ and Effector 3, respectively, with LTP and Kunitz protease inhibitor. This led to the design of an LTP mimic, which sequestered and blocked LasP ₂₃₅ and rescued the bactericidal activity of LTP thereby proving that LasP ₂₃₅ , indeed, participates in HLB pathogenesis.

Ancestral Haplotype Retention and Population Expansion Determine the Complicated Population Genetic Structure of the Hilly Lineage of Neolucanus swinhoei Complex (Coleoptera, Lucanidae) on the Subtropical Taiwan Island

The present study demonstrates that the complicated genetic structure of the hilly lineage of the Neolucanus swinhoei complex was driven by its biological features and habitat requirements as well as hindrance by the CMR during periodical Pleistocene glaciations. The results revealed a tendency of geographical differentiation and major and sub- lineage divergences before and after the Riss glaciation, followed by stable population growth during Würm glaciation. At least four refugia were inferred for N. swinhoei during the Riss-Würm glaciations. The ancestral haplotype retention in the cytochrome oxidase subunit I (COI) gene and compensated substitution in 16S rRNA gene is a possible evolutionary scenario resulting in the inconsistent evolution pattern between COI and 16S rRNA gene coupled with the long-distance dispersal of N. swinhoei. Although the CMR did hinder the dispersal of N. swinhoei, its ancestors may have dispersed to eastern Taiwan through the northern and southern low mountains of the CMR before the Riss glaciation. Our finding suggests that the population growth in the Würm glaciation led a dispersal back to western Taiwan, which is contrast to the more common dispersal scenario from western Taiwan to eastern populations proposed in other studies.

First Account of Phylogeographic Variation, Larval Characters, and Laboratory Rearing of the Endangered Cobblestone Tiger Beetle Cicindelidia marginipennis, Dejean, 1831 with Observations of Their Natural History

Simple Summary

Tiger beetles (Coleoptera: Cicindelidae) are highly predatory and colorful insects of long-standing fascination by entomologists. Most species in this group are in decline from range-wide habitat loss, including those with wide ranges such as the cobblestone tiger beetle (Cicindelidia marginipennis (Dejean, 1831)). This species is considered threatened and/or endangered range-wide, and conservation activities are hindered by a lack of basic information about this species’ biology. Here, we describe new aspects of this species’ basic biology in the lab and the field and use the mitochondrial locus cytochrome oxidase I (COI) to examine biogeographic patterns. We present larval descriptions with detailed drawings and find genetic evidence for geographically structured populations, suggesting that further conservation genetic research is warranted for this species. We expect the new tools and information presented in this paper will allow specialists to test further hypotheses about this species, advance survey methods, and guide proactive tiger beetle conservation.

Abstract

The cobblestone tiger beetle, Cicindelidia marginipennis (Dejean, 1831) is a North American species specializing in riparian habitats from New Brunswick, Canada, to Alabama in the United States. In the United States, this species is state-listed as threatened or endangered range-wide and periodically receives consideration for federal listing, mostly due to habitat decline. Despite its conservation status, intraspecific genetic diversity for this species has not been explored and little is known about its natural history. To support further inquiry into the biology of C. marginipennis, this study provides the first look at range-wide genetic diversity using mitochondrial DNA (mtDNA), describes all three larval instars, and describes natural history characteristics from captive rearing and field observation. Based on mtDNA analyses, our results suggest that geographically based population structure may exist throughout the range, with individuals from Alabama possessing haplotypes not found elsewhere in our sampling. Further genetic analyses, particularly multi-locus analyses, are needed to determine whether the Alabama population represents a separate cryptic species. Our morphological analysis and descriptions of larval instars reveal a combination of characteristics that can be used to differentiate C. marginipennis from closely related and co-occurring species. Based on our field observations, we find that the larval “throw pile” of soil excavated from burrows is a key search image for locating larvae, and we provide descriptions and detailed photographs to aid surveys. Lastly, we find that this species can be successfully reared in captivity and provide guidelines to aid future recovery efforts.

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 densovirus is a persistently infecting virus with a hybrid genome organization and unique transcription strategy

Diaphorina citri densovirus (DcDV) is an ambisense densovirus with a 5071 nt genome. Phylogenetic analysis places DcDV in an intermediate position between those in the Ambidensovirus and Iteradensovirus genera, a finding that is consistent with the observation that DcDV possesses an Iteradensoviris-like non-structural (NS) protein-gene cassette, but a capsid-protein (VP) gene cassette resembling those of other ambisense densoviruses. DcDV is maternally transmitted to 100 % of the progeny of infected female Diaphorina citri, and the progeny of infected females carry DcDV as a persistent infection without outward phenotypic effects. We were unable to infect naïve individuals by oral inoculation, however low levels of transient viral replication are detected following intrathoracic injection of DcDV virions into uninfected D. citri insects. Transcript mapping indicates that DcDV produces one transcript each from the NS and VP gene cassettes and that these transcripts are polyadenylated at internal sites to produce a ~2.2 kb transcript encoding the NS proteins and a ~2.4 kb transcript encoding the VP proteins. Additionally, we found that transcriptional readthrough leads to the production of longer non-canonical transcripts from both genomic strands.

Diaphorina citri (Hemiptera: Psylloidea) in China: Two Invasion Routes and Three Transmission Paths

Diaphorina citri Kuwayama (Hemiptera: Liviidae) is one of the most common pests impacting citrus orchards in southern China. Samples of D. citri were collected in southern China in order to systematically explore the genetic architecture of the species. Mitochondrial cytochrome b (Cytb) and cytochrome coxidase subunit I (COI) were amplified by polymerase chain reaction (PCR) which allowed highlighting low haplotype and nucleotide diversities among the population. Two clades could be observed in the haplotype network. Moreover, Bayesian and maximum parsimony phylogenetic trees were constructed based on the sequences of Cytb and COI. Here, we report on the significant genetic variation of the species when comparing southwestern China with other regions of southern China (southern and southeastern). This analysis also suggested that the genetic structure of D. citri in China originates may from long-term climate fluctuations concomitant with recent disturbances resulting from human activity. Combined with previous data, the present work indicates that D. citri potentially entered China through two distinct invasion routes and spread within the country via three transmission paths.

Identification and Distribution of Haplotypes of Diaphorina citri (Hemiptera: Liviidae) in Jamaica and the Caribbean

Diaphorina citri Kuwayama (Hemiptera: Psyllidae), also known as the Asian citrus psyllid, is a vector of the citrus huanglongbing (HLB) disease. D. citri transmits all three known strains of the HLB pathogen: Candidatus Liberibacter africanus, Candidatus Liberibacter americanus, and Candidatus Liberibacter asiaticus. The study involved 92 psyllids representing the Caribbean Basin and reference samples representing countries within Asia, North America, and South America. This study was aimed at characterizing D. citri on a molecular level in order to determine the haplotype diversity and uniqueness within Jamaica and the Caribbean Basin. D. citri-specific primers were used to amplify an 821 bp gene fragment from the mitochondrial cytochrome c oxidase subunit I gene (mtCOI). The statistical parsimony program, TCS, was used to determine the 12 haplotypes found within the Caribbean, with haplotypes H2 and H7 being the most prominent. The H2 haplotype was found to belong to the South Western Asia group originating from India. H2 represented 54% of the sequenced samples and dominated the Greater Antilles, 22% were grouped as H7, dominating the Lesser Antilles, while the remaining 24% of the sequences were grouped in the remaining 10 haplotypes and were variants seen within the Greater and Lesser Antilles.

Global mtDNA genetic structure and hypothesized invasion history of a major pest of citrus, Diaphorina citri (Hemiptera: Liviidae)

The Asian citrus psyllid Diaphorina citri Kuwayama is a key pest of citrus as the vector of the bacterium causing the “huanglongbing” disease (HLB). To assess the global mtDNA population genetic structure, and possible dispersal history of the pest, we investigated genetic variation at the COI gene collating newly collected samples with all previously published data. Our dataset consists of 356 colonies from 106 geographic sites worldwide. High haplotype diversity (H‐mean = 0.702 ± 0.017), low nucleotide diversity (π‐mean = 0.003), and significant positive selection (Ka/Ks = 32.92) were observed. Forty‐four haplotypes (Hap) were identified, clustered into two matrilines: Both occur in southeastern and southern Asia, North and South America, and Africa; lineages A and B also occur in eastern and western Asia, respectively. The most abundant haplotypes were Hap4 in lineage A (35.67%), and Hap9 in lineage B (41.29%). The haplotype network identified them as the ancestral haplotypes within their respective lineages. Analysis of molecular variance showed significant genetic structure (F_ST = 0.62, p < .0001) between the lineages, and population genetic analysis suggests geographic structuring. We hypothesize a southern and/or southeastern Asia origin, three dispersal routes, and parallel expansions of two lineages. The hypothesized first route involved the expansion of lineage B from southern Asia into North America via West Asia. The second, the expansion of some lineage A individuals from Southeast Asia into East Asia, and the third involved both lineages from Southeast Asia spreading westward into Africa and subsequently into South America. To test these hypotheses and gain a deeper understanding of the global history of D. citri, more data‐rich approaches will be necessary from the ample toolkit of next‐generation sequencing (NGS). However, this study may serve to guide such sampling and in the development of biological control programs against the global pest D. citri.