Genome Resource for the Huanglongbing Causal Agent ‘Candidatus Liberibacter asiaticus’ Strain AHCA17 from Citrus Root Tissue in California, USA
- Weili Cai1 2
- Schyler Nunziata1
- Stefano Costanzo1
- Lucita Kumagai3
- John Rascoe1
- Michael J. Stulberg1 †
- 1United States Department of Agriculture–Animal and Plant Health Inspection Service, Plant Protection and Quarantine, Science and Technology, Beltsville, MD
- 2Department of Entomology and Plant Pathology, North Carolina State University, Raleigh, NC
- 3Plant Pests Diagnostics Branch, California Department of Food and Agriculture, Sacramento, CA
‘Candidatus Liberibacter asiaticus’ is the unculturable causative agent of citrus huanglongbing disease. Here, we report the first citrus root metagenome sequence containing the draft genome of ‘Ca. L. asiaticus’ strain AHCA17, obtained from a pummelo tree in California. The assembled genome was 1.2 Mbp and resulted in 37 contigs (N50 = 158.7 kbp) containing 1,057 predicted open reading frames and 45 RNA-coding genes. This draft genome will provide a valuable resource in further study of ‘Ca. L. asiaticus’ genome diversity and pathogen epidemiology.
‘Candidatus Liberibacter asiaticus’ is an unculturable bacteria and the causal agent of a serious citrus disease, huanglongbing (HLB) (Bové 2006; Jagoueix et al. 1994). HLB causes great economic losses annually in citrus production worldwide (Gottwald et al. 2007). Detection of ‘Ca. L. asiaticus’ is difficult due to uneven distribution in citrus tissues, low pathogen titer, and latency between infections and symptoms. A previous study showed that roots were one of the first tissues to show disease symptoms (Johnson et al. 2014). Currently, publicly available ‘Ca. L. asiaticus’ genomes from citrus samples were all obtained from foliar tissue (Cai et al. 2018; Dai et al. 2019; Wu et al. 2015; Zheng et al. 2014, 2018). Sequencing a ‘Ca. L. asiaticus’ genome from root tissue presents new challenges, considering the complex microbiota associated with this tissue (Zhang et al. 2017). Making additional ‘Ca. L. asiaticus’ genomes public provides a valuable supporting resource to understanding ‘Ca. L. asiaticus’ genome diversity, pathogen movement, and host interaction. Thus, here, we report the first draft whole-genome sequence of ‘Ca. L. asiaticus’ strain AHCA17, obtained directly from the root tissue of a pummelo tree (Citrus maxima) in Orange County, CA.
Genomic DNA of strain AHCA17 was extracted from root tissue using the DNeasy PowerPlant Pro kit (Qiagen, Valencia, CA, U.S.A.). The pathogen titer level was estimated using HLB nrdB-based real-time quantitative PCR, producing a threshold cycle value of 14.36 (Zheng et al. 2016), which signifies relatively high pathogen levels that are critical for successful sequencing. The 2 × 300-bp Illumina paired-end sequencing was performed on an Illumina MiSeq platform (Illumina, Inc., San Diego, CA, U.S.A.). The paired-end sequencing libraries were prepared using an Illumina Nextera DNA Flex library prep kit.
In total, 3.19 × 107 reads with a mean of 269 bp/read were generated from MiSeq sequencing. Kraken (version 1.1.1) (Wood and Salzberg 2014) analysis revealed that 10.28% of the reads mapped to bacteria and not plant host (Table 1; mapped to complete bacterial, archeal, and viral Refseq genomes as of 18 July 2019). This is in contrast to a citrus leaf midrib sample we sequenced (5.16 × 107 reads with the same mean read length), which only had 4.19% of the reads map to bacteria (Cai et al. 2018). Interestingly, within bacteria, the Proteobacteria phylum was represented similarly between the root and leaf sample (41.48 and 42.88%, respectively). Within proteobacteria, there is an increase in diversity in the root sample compared with the leaf, exemplified by Liberibacter spp. representing 54.51% of the Rhizobiaceae family in the root sample and representing 91.68% in the leaf. The larger diversity of closely related organisms in the root led us to a map-to-pathogen, then assemble approach to analyzing the data.
Using the ‘Ca. L. asiaticus’ genome strain psy62 as a reference, 104,637 reads with 25,536,511 bases were aligned using Bowtie2 (version 2.3.0) covering 99.8% of the psy62 genome (Duan et al. 2009; Langmead and Salzberg 2012). Aligned reads were assembled using the SPAdes assembler (version 3.10.0) (Nurk et al. 2013) generating 37 contigs ≥500 bp, ranging from 500 to 266,672 bp (N50 = 158.7 kbp), with 24.56 average coverage. The contigs comprised a total length of 1,208,622 bp, with a G+C content of 36.39%, and aligned with 98.08% of the complete psy62 genome. Annotation was performed by prokka (version 1.13.3) (Seemann 2014). The genome included 1,057 predicted open reading frames and 45 RNA-coding genes.
This Whole-Genome Shotgun project has been deposited at DNA Data Bank of Japan/European Nucleotide Archive/GenBank under the accession VNFL00000000. The version described in this article is version VNFL01000000. Data has also been deposited in the Sequence Read Archive database (PRJNA552562).
The author(s) declare no conflict of interest.
- 2006. Huanglongbing: A destructive, newly emerging, century-old disease of citrus. J. Plant Pathol. 88:7-37. ISI, Google Scholar
- 2018. Draft whole-genome sequence of “Candidatus Liberibacter asiaticus” strain TX1712 from citrus in Texas. Genome Announce. 6:e00554-18. https://doi.org/10.1128/genomeA.00554-18 Crossref, Google Scholar
- 2019. Prophage diversity of ‘Candidatus Liberibacter asiaticus’ strains in California. Phytopathology 109:551-559. https://doi.org/10.1094/PHYTO-06-18-0185-R Link, ISI, Google Scholar
- 2009. Complete genome sequence of citrus huanglongbing bacterium, ‘Candidatus Liberibacter asiaticus’ obtained through metagenomics. Mol. Plant-Microbe Interact. 22:1011-1020. https://doi.org/10.1094/MPMI-22-8-1011 Link, ISI, Google Scholar
- 2007. Citrus huanglongbing: The pathogen and its impact. Plant Health Prog. 8. doi.org/10.1094/PHP-2007-0906-01-RV Link, Google Scholar
- 1994. The phloem-limited bacterium of greening disease of citrus is a member of the α subdivision of the Proteobacteria. Int. J. Syst. Bacteriol. 44:379-386. https://doi.org/10.1099/00207713-44-3-379 Crossref, Google Scholar
- 2014. Association of ‘Candidatus Liberibacter asiaticus’ root infection, but not phloem plugging with root loss on huanglongbing‐affected trees prior to appearance of foliar symptoms. Plant Pathol. 63:290-298. https://doi.org/10.1111/ppa.12109 Crossref, ISI, Google Scholar
- 2012. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9:357-359. https://doi.org/10.1038/nmeth.1923 Crossref, ISI, Google Scholar
- 2013. Assembling genomes and mini-metagenomes from highly chimeric reads. Pages 158-170 in: Research in Computational Molecular Biology. M. Deng, R. Jiang, F. Sun, and X. Zhang, eds. RECOMB 2013. Lecture Notes in Computer Science, vol 7821. Springer, Berlin, Heidelberg, Germany. Crossref, Google Scholar
- 2014. Prokka: Rapid prokaryotic genome annotation. Bioinformatics 30:2068-2069. https://doi.org/10.1093/bioinformatics/btu153 Crossref, ISI, Google Scholar
- 2014. Kraken: Ultrafast metagenomic sequencing classification using exact alignments. Genome Biol. 15:R46. https://doi.org/10.1186/gb-2014-15-3-r46 Crossref, ISI, Google Scholar
- 2015. Draft genome sequence of “Candidatus Liberibacter asiaticus” from Diaphorina citri in Guangdong, China. Genome Announce. 3:e01316-15. https://doi.org/10.1128/genomeA.01316-15 Crossref, Google Scholar
- 2017. Huanglongbing impairs the rhizosphere-to-rhizoplane enrichment process of the citrus root-associated microbiome. Microbiome 5:97. https://doi.org/10.1186/s40168-017-0304-4 Crossref, ISI, Google Scholar
- 2018. A type 3 prophage of ‘Candidatus Liberibacter asiaticus’ carrying a restriction-modification system. Phytopathology 108:454-461. https://doi.org/10.1094/PHYTO-08-17-0282-R Link, ISI, Google Scholar
- 2014. Draft genome sequence of “Candidatus Liberibacter asiaticus” from California. Genome Announce. 2:e00999-14. Crossref, Google Scholar
- 2016. Unusual five copies and dual forms of nrdB in “Candidatus Liberibacter asiaticus”: Biological implications and PCR detection application. Sci. Rep. 6:39020. https://doi.org/10.1038/srep39020 Crossref, ISI, Google Scholar
The findings and conclusions in this publication are those of the authors and should not be construed to represent any official USDA or U.S. Government determination or policy.
The author(s) declare no conflict of interest.
Funding: This research was supported by the United States Department of Agriculture (USDA) Animal and Plant Health Inspection Service.