Detectability of Hibiscus Mealybug, Nipaecoccus viridis (Hemiptera: Pseudoccocidae), DNA in the Mealybug Destroyer, Cryptolaemus montrouzieri (Coleoptera: Coccinellidae), and Survey of Its Predators in Florida Citrus Groves.

The Hibiscus mealybug, Nipaecoccus viridis (Newstead), has recently established in Florida citrus and become a pest of concern given secondary pest outbreaks associated with management of citrus greening disease. Chemical controls used to manage other citrus arthropod pests are not as effective against N. viridis due to its waxy secretions, clumping behavior, and induced cellular changes to host plant tissue which increase microhabitats. Populations of this mealybug pest are regulated by natural enemies in its native region, but it remains unclear if resident natural enemies in Florida citrus could similarly suppress N. viridis populations. This investigation: 1) established species-specific primers for N. viridis based on the mitochondrial gene Cytochrome-oxidase 1 (COI), 2) determined duration of N. viridis DNA detectability in a known predator, the mealybug destroyer (Cryptolaemus montrouzieri Mulsant), by using identified primers in molecular gut content analysis, and 3) screened field-collected predators for the presence of N. viridis DNA. The detection rate of N. viridis DNA was >50% at 36 h after adult C. montrouzieri feeding but DNA was no longer detectable by 72 h after feeding. Field-collected predators were largely comprised of spiders, lacewings, and C. montrouzieri. Spiders, beetles (primarily C. montrouzieri), and juvenile lacewings were the most abundant predators of N. viridis, with 17.8, 43.5, and 58.3 of field-collected samples testing positive for N. viridis DNA, respectively. Our results indicate that Florida citrus groves are hosts to abundant predators of N. viridis and encourage the incorporation of conservation or augmentative biological control for management of this pest.

Complete de novo assembly of Wolbachia endosymbiont of Diaphorina citri Kuwayama (Hemiptera: Liviidae) using long-read genome sequencing.

Wolbachia, a gram-negative [Formula: see text]-proteobacterium, is an endosymbiont found in some arthropods and nematodes. Diaphorina citri Kuwayama, the vector of 'Candidatus Liberibacter asiaticus' (CLas), are naturally infected with a strain of Wolbachia (wDi), which has been shown to colocalize with the bacteria pathogens CLas, the pathogen associated with huanglongbing (HLB) disease of citrus. The relationship between wDi and CLas is poorly understood in part because the complete genome of wDi has not been available. Using high-quality long-read PacBio circular consensus sequences, we present the largest complete circular wDi genome among supergroup-B members. The assembled circular chromosome is 1.52 megabases with 95.7% genome completeness with contamination of 1.45%, as assessed by checkM. We identified Insertion Sequences (ISs) and prophage genes scattered throughout the genomes. The proteins were annotated using Pfam, eggNOG, and COG that assigned unique domains and functions. The wDi genome was compared with previously sequenced Wolbachia genomes using pangenome and phylogenetic analyses. The availability of a complete circular chromosome of wDi will facilitate understanding of its role within the insect vector, which may assist in developing tools for disease management. This information also provides a baseline for understanding phylogenetic relationships among Wolbachia of other insect vectors.

Comparative transcriptome analysis of thiamethoxam susceptible and resistant Asian citrus psyllid, Diaphorina citri (Hemiptera: Liviidae), using RNA-sequencing.

Asian citrus psyllid, Diaphorina citri Kuwayama (Hemiptera: Liviidae), transmits the causal pathogen of huanglongbing and is a global pest of citrus. D. citri populations exhibit resistance to multiple insecticide modes of action in areas where these chemicals have been overused. We performed genome-wide transcriptional analysis for two field populations of D. citri (Wauchula and Lake Alfred, Florida, USA) that exhibit 1300-fold resistance to the neonicotinoid insecticide, thiamethoxam, and compared it to that of susceptible psyllids collected from the same area and without imposed selection. The Lake Alfred population responded to insecticide resistance by up-regulation of 240 genes and down-regulation of 148 others. The Wauchula population exhibited similar patterns to the Lake Alfred population with up-regulation of 253 genes and down-regulation of 115 others. Gene Ontology annotation associated with cellular processes, cell, and catalytic activity were assigned to differentially expressed genes (DEGs). The DEGs from Lake Alfred and Wauchula populations were mapped to Kyoto Encyclopedia of Gene and Genomes pathways and implicated enrichment of metabolic pathways, oxidative phosphorylation, extracellular matrix-receptor interaction, terpenoid backbone biosynthesis, and insect hormone biosynthesis in the resistant populations. Up-regulation of 60s ribosomal proteins, UDP-gluscoyltransferases, cytochrome c oxidases, and CYP and ABC transporters among thiamethoxam-resistant D. citri implicates a broad array of novel and conventionally understood resistance mechanisms.

Insecticide rotation scheme restores insecticide susceptibility in thiamethoxam-resistant field populations of Asian citrus psyllid, Diaphorina citri Kuwayama (Hemiptera: Liviidae), in Florida.

BACKGROUND: We investigated rotation using insecticides with multiple modes of action as a resistance management strategy for Asian citrus psyllid, Diaphorina citri, Kuwayama (Hemiptera: Liviidae), in Florida. The stability of thiamethoxam resistance was investigated in the laboratory by establishing populations of field-collected, resistant D. citri and rearing them under no insecticide exposure. Furthermore, recovery of susceptibility was investigated in the field by initiating rotation to insecticides in plots that previously were treated with consecutive thiamethoxam applications. RESULTS: The resistance ratio (RR) for thiamethoxam reached between 1266.29- and 1395.00-fold after three and four consecutive applications of thiamethoxam, respectively. However, the RR for thiamethoxam remained low (1.71-5.28-fold) under both rotations at both Lake Alfred and Wauchula. Thiamethoxam was cross-resistant with imidacloprid (RR site 1 = 1059.65-fold, RR site 2 = 1595.43-fold) and clothianidin (RR site 1 = 1798.78-fold, RR site 2 = 1270.57-fold) in the nonrotated treatment at both sites. There was very low cross-resistance to other insecticides with different modes of action. Both laboratory and field investigations indicated that susceptibility to thiamethoxam fully recovered after five D. citri generations. Expression of CYP4C67 was significantly increased in resistant populations. CONCLUSION: Our results revealed that D. citri populations develop a high level of resistance following only three or four consecutive neonicotinoid sprays; this was associated with subsequent product failure. Our data suggest that metabolic detoxification by cytochrome P450s contributes to thiamethoxam resistance in D. citri. Overall, the investigation demonstrated that resistance to thiamethoxam can be managed readily in populations of D. citri by rotating modes of action.

Near-Complete Genome Sequences of a Wolbachia Strain Isolated from Diaphorina citri Kuwayama (Hemiptera: Liviidae).

Wolbachia strains are one of three endosymbionts associated with the insect vector of "Candidatus Liberibacter asiaticus," Diaphorina citri Kuwayama (Hemiptera: Liviidae). We report three near-complete genome sequences of samples of Wolbachia from D. citri (wDi), with sizes of 1,518,595, 1,542,468, and 1,538,523 bp.

Molecular Basis of Soybean Resistance to Soybean Aphids and Soybean Cyst Nematodes.

Soybean aphid (SBA; Aphis glycines Matsumura) and soybean cyst nematode (SCN; Heterodera glycines Ichninohe) are major pests of the soybean (Glycine max [L.] Merr.). Substantial progress has been made in identifying the genetic basis of limiting these pests in both model and non-model plant systems. Classical linkage mapping and genome-wide association studies (GWAS) have identified major and minor quantitative trait loci (QTLs) in soybean. Studies on interactions of SBA and SCN effectors with host proteins have identified molecular cues in various signaling pathways, including those involved in plant disease resistance and phytohormone regulations. In this paper, we review the molecular basis of soybean resistance to SBA and SCN, and we provide a synthesis of recent studies of soybean QTLs/genes that could mitigate the effects of virulent SBA and SCN populations. We also review relevant studies of aphid-nematode interactions, particularly in the soybean-SBA-SCN system.

Transcriptome profiling of interaction effects of soybean cyst nematodes and soybean aphids on soybean.

Soybean aphid (Aphis glycines; SBA) and soybean cyst nematode (Heterodera glycines; SCN) are two major pests of soybean (Glycine max) in the United States of America. This study aims to characterize three-way interactions among soybean, SBA, and SCN using both demographic and genetic datasets. SCN-resistant and SCN-susceptible soybean cultivars with a combination of soybean aphids (biotype 1) and SCN (HG type 0) in a randomized complete block design (RCBD) with six blocks were used to evaluate the three-way interactions in a greenhouse setup. Treatments receiving SCN were infested at planting with 2000 nematode eggs, and the treatments with soybean aphids were infested at second trifoliate growth stage (V2) with 15 soybean aphids. The whole roots were sampled from plants at 5 and 30 days post SBA infestation for RNA sequencing using Illumina Hiseq. 3000. The data comprises of 47 libraries that are useful for further analyses of important genes, which are involved in interaction effects of SBA and SCN on soybean.

Transcriptome profiling of induced susceptibility effects on soybean-soybean aphid (Hemiptera: Aphididae) interaction.

OBJECTIVES: Soybean aphid (Aphis glycines Matsumura; SBA) is the most economically damaging insect of soybean (Glycine max) in the United States. One previous study demonstrated that avirulent (biotype 1) and virulent (biotype 2) biotypes could co-occur and interact on resistant (i.e., Rag1) and susceptible soybean resulting in induced susceptibility after 11 days of feeding. The main objective of this research was to employ RNA sequencing (RNA-seq) technique to compare the induced susceptibility effect of biotype 2 on susceptible and resistant soybean at day 1 and day 11 (i.e., both susceptible and resistant soybean were initially challenged by biotype 2 and the effect was monitored through biotype 1 populations). DATA DESCRIPTION: We investigated susceptible and Rag1 transcriptome response to SBA feeding in soybean plants colonized by biotype 1 in the presence or absence of an inducer population (i.e., biotype 2). Ten RNA datasets are reported with 266,535,654 sequence reads (55.2 GB) obtained from pooled samples derived from the leaves collected at day 1 and day 11 post SBA infestation. A comprehensive understanding of these transcriptome data will enhance our understanding of interactions among soybean and two different biotypes of soybean aphids at the molecular level.

Identification and Characterization of Mitogen-Activated Protein Kinase (MAPK) Genes in Sunflower (Helianthus annuus L.).

Mitogen-Activated Protein Kinase (MAPK) genes encode proteins that regulate biotic and abiotic stresses in plants through signaling cascades comprised of three major subfamilies: MAP Kinase (MPK), MAPK Kinase (MKK), and MAPKK Kinase (MKKK). The main objectives of this research were to conduct genome-wide identification of MAPK genes in Helianthus annuus and examine functional divergence of these genes in relation to those in nine other plant species (Amborella trichopoda, Aquilegia coerulea, Arabidopsis thaliana, Daucus carota, Glycine max, Oryza sativa, Solanum lycopersicum, Sphagnum fallax, and Vitis vinifera), representing diverse taxonomic groups of the Plant Kingdom. A Hidden Markov Model (HMM) profile of the MAPK genes utilized reference sequences from A. thaliana and G. max, yielding a total of 96 MPKs and 37 MKKs in the genomes of A. trichopoda, A. coerulea, C. reinhardtii, D. carota, H. annuus, S. lycopersicum, and S. fallax. Among them, 28 MPKs and eight MKKs were confirmed in H. annuus. Phylogenetic analyses revealed four clades within each subfamily. Transcriptomic analyses showed that at least 19 HaMPK and seven HaMKK genes were induced in response to salicylic acid (SA), sodium chloride (NaCl), and polyethylene glycol (Peg) in leaves and roots. Of the seven published sunflower microRNAs, five microRNA families are involved in targeting eight MPKs. Additionally, we discussed the need for using MAP Kinase nomenclature guidelines across plant species. Our identification and characterization of MAP Kinase genes would have implications in sunflower crop improvement, and in advancing our knowledge of the diversity and evolution of MAPK genes in the Plant Kingdom.

Genome-Wide Identification of NBS-Encoding Resistance Genes in Sunflower (Helianthus annuus L.).

Nucleotide Binding Site-Leucine-Rich Repeat (NBS-LRR) genes encode disease resistance proteins involved in plants' defense against their pathogens. Although sunflower is affected by many diseases, only a few molecular details have been uncovered regarding pathogenesis and resistance mechanisms. Recent availability of sunflower whole genome sequences in publicly accessible databases allowed us to accomplish a genome-wide identification of Toll-interleukin-1 receptor-like Nucleotide-binding site Leucine-rich repeat (TNL), Coiled Coil (CC)-NBS-LRR (CNL), Resistance to powdery mildew 8 (RPW8)-NBS-LRR (RNL) and NBS-LRR (NL) protein encoding genes. Hidden Markov Model (HMM) profiling of 52,243 putative protein sequences from sunflower resulted in 352 NBS-encoding genes, among which 100 genes belong to CNL group including 64 genes with RX_CC like domain, 77 to TNL, 13 to RNL, and 162 belong to NL group. We also identified signal peptides and nuclear localization signals present in the identified genes and their homologs. We found that NBS genes were located on all chromosomes and formed 75 gene clusters, one-third of which were located on chromosome 13. Phylogenetic analyses between sunflower and Arabidopsis NBS genes revealed a clade-specific nesting pattern in CNLs, with RNLs nested in the CNL-A clade, and species-specific nesting pattern for TNLs. Surprisingly, we found a moderate bootstrap support (BS = 50%) for CNL-A clade being nested within TNL clade making both the CNL and TNL clades paraphyletic. Arabidopsis and sunflower showed 87 syntenic blocks with 1049 high synteny hits between chromosome 5 of Arabidopsis and chromosome 6 of sunflower. Expression data revealed functional divergence of the NBS genes with basal level tissue-specific expression. This study represents the first genome-wide identification of NBS genes in sunflower paving avenues for functional characterization and potential crop improvement.

Evolutionary Divergence of TNL Disease-Resistant Proteins in Soybean (Glycine max) and Common Bean (Phaseolus vulgaris).

Disease-resistant genes (R genes) encode proteins that are involved in protecting plants from their pathogens and pests. Availability of complete genome sequences from soybean and common bean allowed us to perform a genome-wide identification and analysis of the Toll interleukin-1 receptor-like nucleotide-binding site leucine-rich repeat (TNL) proteins. Hidden Markov model (HMM) profiling of all protein sequences resulted in the identification of 117 and 77 regular TNL genes in soybean and common bean, respectively. We also identified TNL gene homologs with unique domains, and signal peptides as well as nuclear localization signals. The TNL genes in soybean formed 28 clusters located on 10 of the 20 chromosomes, with the majority found on chromosome 3, 6 and 16. Similarly, the TNL genes in common bean formed 14 clusters located on five of the 11 chromosomes, with the majority found on chromosome 10. Phylogenetic analyses of the TNL genes from Arabidopsis, soybean and common bean revealed less divergence within legumes relative to the divergence between legumes and Arabidopsis. Syntenic blocks were found between chromosomes Pv10 and Gm03, Pv07 and Gm10, as well as Pv01 and Gm14. The gene expression data revealed basal level expression and tissue specificity, while analysis of available microRNA data showed 37 predicted microRNA families involved in targeting the identified TNL genes in soybean and common bean.

Comparative Genomics of Non-TNL Disease Resistance Genes from Six Plant Species.

Disease resistance genes (R genes), as part of the plant defense system, have coevolved with corresponding pathogen molecules. The main objectives of this project were to identify non-Toll interleukin receptor, nucleotide-binding site, leucine-rich repeat (nTNL) genes and elucidate their evolutionary divergence across six plant genomes. Using reference sequences from Arabidopsis, we investigated nTNL orthologs in the genomes of common bean, Medicago, soybean, poplar, and rice. We used Hidden Markov Models for sequence identification, performed model-based phylogenetic analyses, visualized chromosomal positioning, inferred gene clustering, and assessed gene expression profiles. We analyzed 908 nTNL R genes in the genomes of the six plant species, and classified them into 12 subgroups based on the presence of coiled-coil (CC), nucleotide binding site (NBS), leucine rich repeat (LRR), resistance to Powdery mildew 8 (RPW8), and BED type zinc finger domains. Traditionally classified CC-NBS-LRR (CNL) genes were nested into four clades (CNL A-D) often with abundant, well-supported homogeneous subclades of Type-II R genes. CNL-D members were absent in rice, indicating a unique R gene retention pattern in the rice genome. Genomes from Arabidopsis, common bean, poplar and soybean had one chromosome without any CNL R genes. Medicago and Arabidopsis had the highest and lowest number of gene clusters, respectively. Gene expression analyses suggested unique patterns of expression for each of the CNL clades. Differential gene expression patterns of the nTNL genes were often found to correlate with number of introns and GC content, suggesting structural and functional divergence.

Diversity and Evolution of Disease Resistance Genes in Barley (Hordeum vulgare L.).

Plant disease resistance genes (R-genes) play a critical role in the defense response to pathogens. Barley is one of the most important cereal crops, having a genome recently made available, for which the diversity and evolution of R-genes are not well understood. The main objectives of this research were to conduct a genome-wide identification of barley Coiled-coil, Nucleotide-binding site, Leucine-rich repeat (CNL) genes and elucidate their evolutionary history. We employed a Hidden Markov Model using 52 Arabidopsis thaliana CNL reference sequences and analyzed for phylogenetic relationships, structural variation, and gene clustering. We identified 175 barley CNL genes nested into three clades, showing (a) evidence of an expansion of the CNL-C clade, primarily due to tandem duplications; (b) very few members of clade CNL-A and CNL-B; and (c) a complete absence of clade CNL-D. Our results also showed that several of the previously identified mildew locus A (MLA) genes may be allelic variants of two barley CNL genes, MLOC_66581 and MLOC_10425, which respond to powdery mildew. Approximately 23% of the barley CNL genes formed 15 gene clusters located in the extra-pericentromeric regions on six of the seven chromosomes; more than half of the clustered genes were located on chromosomes 1H and 7H. Higher average numbers of exons and multiple splice variants in barley relative to those in Arabidopsis and rice may have contributed to a diversification of the CNL-C members. These results will help us understand the evolution of R-genes with potential implications for developing durable resistance in barley cultivars.