Graph pangenome captures missing heritability and empowers tomato breeding.

Missing heritability in genome-wide association studies defines a major problem in genetic analyses of complex biological traits1,2. The solution to this problem is to identify all causal genetic variants and to measure their individual contributions3,4. Here we report a graph pangenome of tomato constructed by precisely cataloguing more than 19 million variants from 838 genomes, including 32 new reference-level genome assemblies. This graph pangenome was used for genome-wide association study analyses and heritability estimation of 20,323 gene-expression and metabolite traits. The average estimated trait heritability is 0.41 compared with 0.33 when using the single linear reference genome. This 24% increase in estimated heritability is largely due to resolving incomplete linkage disequilibrium through the inclusion of additional causal structural variants identified using the graph pangenome. Moreover, by resolving allelic and locus heterogeneity, structural variants improve the power to identify genetic factors underlying agronomically important traits leading to, for example, the identification of two new genes potentially contributing to soluble solid content. The newly identified structural variants will facilitate genetic improvement of tomato through both marker-assisted selection and genomic selection. Our study advances the understanding of the heritability of complex traits and demonstrates the power of the graph pangenome in crop breeding.

Draft Genome Sequence Resource of the Citrus Stem-End Rot Fungal Pathogen Lasiodiplodia theobromae CITRA15.

Lasiodiplodia theobromae is a fungal pathogen associated with perennial tropical fruit plants worldwide. In citrus, L. theobromae causes stem-end rot (Diplodia stem-end rot), a damaging postharvest disease that is aggravated when trees are also infected with the citrus greening bacteria 'Candidatus Liberibacter asiaticus'. Due to the latent infection of L. theobromae during the preharvest stage, it becomes difficult to control the disease by chemical or physical treatment. In the current study, we sequenced and assembled strain CITRA15, the first genome of L. theobromae obtained from diseased Citrus paradise 'Flame' grapefruit in Florida, and thereby provided a genomic resource for future research on diagnostics, and postharvest and preharvest disease management of citrus and other fruit crops.

Longitudinal Transcriptomic, Proteomic, and Metabolomic Analysis of Citrus limon Response to Graft Inoculation by Candidatus Liberibacter asiaticus.

Presymptomatic detection of citrus trees infected with Candidatus Liberibacter asiaticus (CLas), the bacterial pathogen associated with Huanglongbing (HLB; citrus greening disease), is critical to controlling the spread of the disease. To test whether infected citrus trees produce systemic signals that may be used for indirect disease detection, lemon (Citrus limon) plants were graft-inoculated with either CLas-infected or control (CLas-) budwood, and leaf samples were longitudinally collected over 46 weeks and analyzed for plant changes associated with CLas infection. RNA, protein, and metabolite samples extracted from leaves were analyzed using RNA-Seq, mass spectrometry, and 1H NMR spectroscopy, respectively. Significant differences in specific transcripts, proteins, and metabolites were observed between CLas-infected and control plants as early as 2 weeks post graft (wpg). The most dramatic differences between the transcriptome and proteome of CLas-infected and control plants were observed at 10 wpg, including coordinated increases in transcripts and proteins of citrus orthologs of known plant defense genes. This integrated approach to quantifying plant molecular changes in leaves of CLas-infected plants supports the development of diagnostic technology for presymptomatic or early disease detection as part of efforts to control the spread of HLB into uninfected citrus groves.

Longitudinal Transcriptomic, Proteomic, and Metabolomic Analyses of Citrus sinensis (L.) Osbeck Graft-Inoculated with "Candidatus Liberibacter asiaticus".

"Candidatus Liberibacter asiaticus" (CLas) is the bacterium associated with the citrus disease Huanglongbing (HLB). Current CLas detection methods are unreliable during presymptomatic infection, and understanding CLas pathogenicity to help develop new detection techniques is challenging because CLas has yet to be isolated in pure culture. To understand how CLas affects citrus metabolism and whether infected plants produce systemic signals that can be used to develop improved detection techniques, leaves from Washington Navel orange (Citrus sinensis (L.) Osbeck) plants were graft-inoculated with CLas and longitudinally studied using transcriptomics (RNA sequencing), proteomics (liquid chromatography-tandem mass spectrometry), and metabolomics (proton nuclear magnetic resonance). Photosynthesis gene expression and protein levels were lower in infected plants compared to controls during late infection, and lower levels of photosynthesis proteins were identified as early as 8 weeks post-grafting. These changes coordinated with higher sugar concentrations, which have been shown to accumulate during HLB. Cell wall modification and degradation gene expression and proteins were higher in infected plants during late infection. Changes in gene expression and proteins related to plant defense were observed in infected plants as early as 8 weeks post-grafting. These results reveal coordinated changes in greenhouse navel leaves during CLas infection at the transcript, protein, and metabolite levels, which can inform of biomarkers of early infection.

A quick guide for student-driven community genome annotation.

High quality gene models are necessary to expand the molecular and genetic tools available for a target organism, but these are available for only a handful of model organisms that have undergone extensive curation and experimental validation over the course of many years. The majority of gene models present in biological databases today have been identified in draft genome assemblies using automated annotation pipelines that are frequently based on orthologs from distantly related model organisms and usually have minor or major errors. Manual curation is time consuming and often requires substantial expertise, but is instrumental in improving gene model structure and identification. Manual annotation may seem to be a daunting and cost-prohibitive task for small research communities but involving undergraduates in community genome annotation consortiums can be mutually beneficial for both education and improved genomic resources. We outline a workflow for efficient manual annotation driven by a team of primarily undergraduate annotators. This model can be scaled to large teams and includes quality control processes through incremental evaluation. Moreover, it gives students an opportunity to increase their understanding of genome biology and to participate in scientific research in collaboration with peers and senior researchers at multiple institutions.

The Ptr1 Locus of Solanum lycopersicoides Confers Resistance to Race 1 Strains of Pseudomonas syringae pv. tomato and to Ralstonia pseudosolanacearum by Recognizing the Type III Effectors AvrRpt2 and RipBN.

Race 1 strains of Pseudomonas syringae pv. tomato, which cause bacterial speck disease of tomato, are becoming increasingly common and no simply inherited genetic resistance to such strains is known. We discovered that a locus in Solanum lycopersicoides, termed Pseudomonas tomato race 1 (Ptr1), confers resistance to race 1 P. syringae pv. tomato strains by detecting the activity of type III effector AvrRpt2. In Arabidopsis, AvrRpt2 degrades the RIN4 protein, thereby activating RPS2-mediated immunity. Using site-directed mutagenesis of AvrRpt2, we found that, like RPS2, activation of Ptr1 requires AvrRpt2 proteolytic activity. Ptr1 also detected the activity of AvrRpt2 homologs from diverse bacteria, including one in Ralstonia pseudosolanacearum. The genome sequence of S. lycopersicoides revealed no RPS2 homolog in the Ptr1 region. Ptr1 could play an important role in controlling bacterial speck disease and its future cloning may shed light on an example of convergent evolution for recognition of a widespread type III effector.

The Sol Genomics Network (SGN)--from genotype to phenotype to breeding.

The Sol Genomics Network (SGN, http://solgenomics.net) is a web portal with genomic and phenotypic data, and analysis tools for the Solanaceae family and close relatives. SGN hosts whole genome data for an increasing number of Solanaceae family members including tomato, potato, pepper, eggplant, tobacco and Nicotiana benthamiana. The database also stores loci and phenotype data, which researchers can upload and edit with user-friendly web interfaces. Tools such as BLAST, GBrowse and JBrowse for browsing genomes, expression and map data viewers, a locus community annotation system and a QTL analysis tools are available. A new tool was recently implemented to improve Virus-Induced Gene Silencing (VIGS) constructs called the SGN VIGS tool. With the growing genomic and phenotypic data in the database, SGN is now advancing to develop new web-based breeding tools and implement the code and database structure for other species or clade-specific databases.

Comprehensive repeatome annotation reveals strong potential impact of repetitive elements on tomato ripening.

BACKGROUND: Plant genomes are populated by different types of repetitive elements including transposable elements (TEs) and simple sequence repeats (SSRs) that can have a strong impact on genome size and dynamic as well as on the regulation of gene transcription. At least two-thirds of the tomato genome is composed of repeats. While their bulk impact on genome organization has been recently revealed by whole genome assembly, their influence on tomato biology and phenotype remains largely unaddressed. More specifically, the effects and roles of DNA repeats on the maturation of fleshy fruits, which is a complex process of key agro-economic interest, still needs to be investigated comprehensively and tomato is arguably an excellent model for such study. RESULTS: We have performed a comprehensive annotation of the tomato repeatome to explore its potential impact on tomato genome composition and gene transcription. Our results show that the tomato genome can be fractioned into three compartments with different gene and repeat density, each compartment presenting contrasting repeat and gene composition, repeat-gene associations and different gene transcriptional levels. In the context of fruit ripening, we found that repeats are present in the majority of differentially methylated regions (DMRs) and thousands of repeat-associated DMRs are found in gene proximity including hundreds that are differentially regulated. Furthermore, we found that repeats are also present in the proximity of binding sites of the key ripening protein RIN. We also observed that some repeat families are present at unexpected high frequency in the proximity of genes that are differentially expressed during tomato ripening. CONCLUSION: Altogether, our study emphasizes the fractionation as defined by repeat content in the tomato genome and enables to further characterize the specificities of each genomic compartment. Additionally, our results present strong associations between differentially regulated genes, differentially methylated regions and repeats, suggesting a potential adaptive function of repeats in tomato ripening. Our work therefore provides significant perspectives for the understanding of the impact of repeats on the maturation of fleshy fruits.

Traits of selected Clostridium strains for syngas fermentation to ethanol.

Syngas fermentation is an anaerobic bioprocess that could become industrially relevant as a biorefinery platform for sustainable production of fuels and chemicals. An important prerequisite for commercialization is adequate performance of the biocatalyst (i.e., sufficiently high production rate, titer, selectivity, yield, and stability of the fermentation). Here, we compared the performance of three potential candidate Clostridium strains in syngas-to-ethanol conversion: Clostridium ljungdahlii PETC, C. ljungdahlii ERI-2, and Clostridium autoethanogenum JA1-1. Experiments were conducted in a two-stage, continuously fed syngas-fermentation system that had been optimized for stable ethanol production. The two C. ljungdahlii strains performed similar to each other but different from C. autoethanogenum. When the pH value was lowered from 5.5 to 4.5 to induce solventogenesis, the cell-specific carbon monoxide and hydrogen consumption (similar rate for all strains at pH 5.5), severely decreased in JA1-1, but hardly in PETC and ERI-2. Ethanol production in strains PETC and ERI-2 remained relatively stable while the rate of acetate production decreased, resulting in a high ethanol/acetate ratio, but lower overall productivities. With JA1-1, lowering the pH severely lowered rates of both ethanol and acetate production; and as a consequence, no pronounced shift to solventogenesis was observed. The highest overall ethanol production rate of 0.301 g · L(-1) · h(-1) was achieved with PETC at pH 4.5 with a corresponding 19 g/L (1.9% w/v) ethanol concentration and a 5.5:1 ethanol/acetate molar ratio. A comparison of the genes relevant for ethanol metabolism revealed differences between C. ljungdahlii and C. autoethanogenum that, however, did not conclusively explain the different phenotypes.

Genome-Assisted Development of a Diagnostic Protocol for Distinguishing High Virulence Pseudomonas syringae pv. tomato Strains.

A severe outbreak of bacterial speck of tomato, caused by Pseudomonas syringae pv. tomato, occurred in central New York in 2009. Isolate 09150, collected from this outbreak and subsequently named NYS-T1, was found to be highly virulent on tomato. To better understand the relationship of 09150 to other P. syringae strains and develop a diagnostic assay for aggressive strains of this pathogen, the 09150 genome was sequenced. Genome comparison revealed it to be highly similar to a previously sequenced isolate, T1. Genetic factors linked to host interaction including type III effectors, toxin biosynthetic genes, and elicitors of host innate immunity were identified. Type III effector repertoires were compared with other strains in the high virulence T1-like subgroup and lower virulence DC3000/P. syringae pv. maculicola subgroup within P. syringae phylogenetic Group I. Primers for conventional PCR were developed using sequences for avrA, hopW, conserved in the former subgroup and hopN, present in the latter. These were tested on isolates in the two subgroups, other pseudomonads, and other bacterial pathogens of tomato. Primers developed for avaA and hopW were diagnostic for more virulent strains of P. syringae pv. tomato while primers for hopN were diagnostic for P. syringae pv. tomato DC3000 and related P. syringe pv. maculicola strains. Primers designed against hopR distinguished both of these P. syringae subgroups from other P. syringae strains.

Bound to Succeed: transcription factor binding-site prediction and its contribution to understanding virulence and environmental adaptation in bacterial plant pathogens.

Bacterial plant pathogens rely on a battalion of transcription factors to fine-tune their response to changing environmental conditions and to marshal the genetic resources required for successful pathogenesis. Prediction of transcription factor binding sites (TFBS) represents an important tool for elucidating regulatory networks and has been conducted in multiple genera of plant-pathogenic bacteria for the purpose of better understanding mechanisms of survival and pathogenesis. The major categories of TFBS that have been characterized are reviewed here, with emphasis on in silico methods used for site identification and challenges therein, their applicability to different types of sequence datasets, and insights into mechanisms of virulence and survival that have been gained through binding-site mapping. An improved strategy for establishing E-value cutoffs when using existing models to screen uncharacterized genomes is also discussed.

Differential gene expression of Asian citrus psyllids infected with 'Ca. Liberibacter asiaticus' reveals hyper-susceptibility to invasion by instar fourth-fifth and teneral adult stages.

The bacterial pathogen Candidatus Liberibacter asiaticus (CLas) is the causal agent of citrus greening disease. This unusual plant pathogenic bacterium also infects its psyllid host, the Asian citrus psyllid (ACP). To investigate gene expression profiles with a focus on genes involved in infection and circulation within the psyllid host of CLas, RNA-seq libraries were constructed from CLas-infected and CLas-free ACP representing the five different developmental stages, namely, nymphal instars 1-2, 3, and 4-5, and teneral and mature adults. The Gbp paired-end reads (296) representing the transcriptional landscape of ACP across all life stages and the official gene set (OGSv3) were annotated based on the chromosomal-length v3 reference genome and used for de novo transcript discovery resulting in 25,410 genes with 124,177 isoforms. Differential expression analysis across all ACP developmental stages revealed instar-specific responses to CLas infection, with greater overall responses by nymphal instars, compared to mature adults. More genes were over-or under-expressed in the 4-5th nymphal instars and young (teneral) adults than in instars 1-3, or mature adults, indicating that late immature instars and young maturing adults were highly responsive to CLas infection. Genes identified with potential for direct or indirect involvement in the ACP-CLas circulative, propagative transmission pathway were predominantly responsive during early invasion and infection processes and included canonical cytoskeletal remodeling and endo-exocytosis pathway genes. Genes with predicted functions in defense, development, and immunity exhibited the greatest responsiveness to CLas infection. These results shed new light on ACP-CLas interactions essential for pathogenesis of the psyllid host, some that share striking similarities with effector protein-animal host mechanisms reported for other culturable and/or fastidious bacterial- or viral- host pathosystems.

A chromosome scale assembly of the tarnished plant bug, Lygus lineolaris (Palisot de Beauvois), genome.

OBJECTIVE: The tarnished plant bug (TPB), Lygus lineolaris (Palisot de Beauvois) (Hemiptera: Miridae), is a pest damaging many cultivated crops in North America. Although partial transcriptome data are available for this pest, a genome assembly was not available for this species. This assembly of a high-quality chromosome-length genome of TPB is aimed to develop the genetic resources that can provide the foundation required for advancing research on this species. RESULTS: The initial genome of TPB assembled with paired-end nucleotide sequences generated with Illumina technology was scaffolded with Illumina HiseqX reads generated from a proximity ligated (HiC) library to obtain a high-quality genome assembly. The final assembly contained 3963 scaffolds longer than 1 kbp to yield a genome of 599.96 Mbp. The N50 of the TPB genome assembly was 35.64 Mbp and 98.68% of the genome was assembled into 17 scaffolds larger than 1 Mbp. This megabase scaffold number is the same as the number of chromosomes observed in karyotyping of this insect. The TPB genome is known to have high repetitive DNA content, and the reduced assembled genome size compared to flowcytometric estimates of approximately 860 Mbp may be due to the collapsed assembly of highly similar regions.

Chromosome length genome assembly of the redbanded stink bug, Piezodorus guildinii (Westwood).

OBJECTIVE: The redbanded stink bug (RBSB), Piezodorus guildinii (Hemiptera: Pentatomidae), is native to the Caribbean Basin and is currently considered an invasive pest in Florida, Louisiana, Mississippi, and Texas in the southern United States. Although RBSB is an economically important invasive pest in the USA, relatively few studies have been conducted to understand molecular mechanisms, population genetic structure, and the genetic basis of resistance to insecticides. The objective of this work was to obtain a high-quality genome assembly to develop genomic resources to conduct population genetic, genomic, and physiological studies of the RBSB. RESULTS: The genome of RBSB was sequenced with Pacific Biosciences technology followed by two rounds of scaffolding using Chicago libraries and HiC proximity ligation to obtain a high-quality assembly. The genome assembly contained 800 scaffolds larger than 1 kbp and the N50 was 170.84 Mbp. The largest scaffold was 222.22 Mbp and 90% of the genome was included in the 7 scaffolds larger than 118 Mbp. The number of megabase scaffolds also matched the number of chromosomes in this insect. The genome sequence will facilitate the development of resources to conduct studies on genetics, transcriptomics, and physiology of RBSB.

Annotation of Hox cluster and Hox cofactor genes in the Asian citrus psyllid, Diaphorina citri, reveals novel features.

Hox genes and their cofactors are essential developmental genes specifying regional identity in animals. Hox genes have a conserved arrangement in clusters in the same order in which they specify identity along the anterior-posterior axis. A few insect species have breaks in the cluster, but these are exceptions. We annotated the 10 Hox genes of the Asian citrus psyllid Diaphorina citri, and found a split in its Hox cluster between the Deformed and Sex combs reduced genes - the first time a break at this position has been observed in an insect Hox cluster. We also annotated D. citri orthologs of the Hox cofactor genes homothorax, PKNOX and extradenticle and found an additional copy of extradenticle in D. citri that appears to be a retrogene. Expression data and sequence conservation suggest that the extradenticle retrogene may have retained the original extradenticle function and allowed divergence of the parental extradenticle gene.

Lessons learned about the biology and genomics of Diaphorina citri infection with "Candidatus Liberibacter asiaticus" by integrating new and archived organ-specific transcriptome data.

BACKGROUND: Huanglongbing, a devastating disease of citrus, is caused by the obligate, intracellular bacterium "Candidatus Liberibacter asiaticus" (CLas). CLas is transmitted by Diaphorina citri, the Asian citrus psyllid. Development of transmission-blocking strategies to manage huanglongbing relies on knowledge of CLas and D. citri interactions at the molecular level. Prior transcriptome analyses of D. citri point to changes in psyllid biology due to CLas infection but have been hampered by incomplete versions of the D. citri genome, proper host plant controls, and/or a lack of a uniform data analysis approach. In this work, we present lessons learned from a quantitative transcriptome analysis of excised heads, salivary glands, midguts, and bacteriomes from CLas-positive and CLas-negative D. citri using the chromosomal length D. citri genome assembly. RESULTS: Each organ had a unique transcriptome profile and response to CLas infection. Though most psyllids were infected with the bacterium, CLas-derived transcripts were not detected in all organs. By analyzing the midgut dataset using both the Diaci_v1.1 and v3.0 D. citri genomes, we showed that improved genome assembly led to significant and quantifiable differences in RNA-sequencing data interpretation. CONCLUSIONS: Our results support the hypothesis that future transcriptome studies on circulative, vector-borne pathogens should be conducted at the tissue-specific level using complete, chromosomal-length genome assemblies for the most accurate understanding of pathogen-induced changes in vector gene expression.

Genomic identification, annotation, and comparative analysis of Vacuolar-type ATP synthase subunits in Diaphorina citri.

The hemipteran insect Diaphorina citri, or Asian citrus psyllid, is a vector for Candidatus Liberibacter asiaticus (CLas), the bacterium causing citrus greening disease, or Huanglongbing (HLB). Millions of citrus trees have been destroyed, and every grove in Florida, USA, has been directly affected by this disease. In eukaryotes, vacuolar-type ATP synthase (V-ATPase) is an abundant heterodimeric enzyme that serves the cell with essential compartment acidification through the active processes that transport protons across the membrane. Fifteen putative V-ATPase genes in the D. citri genome were manually curated. Comparative genomic analysis revealed that D. citri V-ATPase subunits share domains and motifs with other insects, including the V-ATPase-A superfamily domain. Phylogenetic analysis separates D. citri V-ATPase subunits into expected clades with orthologous sequences. Annotation of the D. citri genome is a critical step towards developing directed pest management strategies to reduce the spread of HLB throughout the citrus industry.

Manual curation and phylogenetic analysis of chitinase family genes in the Asian citrus psyllid, Diaphorina citri.

Chitinases are enzymes that digest the polysaccharide polymer chitin. During insect development, breakdown of chitin is an essential step in molting of the exoskeleton. Knockdown of chitinases required for molting is lethal to insects, making chitinase genes an interesting target for RNAi-based pest control methods. The Asian citrus psyllid, Diaphorina citri, carries the bacterium causing Huanglongbing, or citrus greening disease, a devastating citrus disease. We identified and annotated 12 chitinase family genes from D. citri as part of a community effort to create high-quality gene models to aid the design of interdictory molecules for pest control. We categorized the D. citri chitinases according to an established classification scheme and re-evaluated the classification of chitinases in other hemipterans. In addition to chitinases from known groups, we identified a novel class of chitinases present in D. citri and several related hemipterans that appears to be the result of horizontal gene transfer.

Ubiquitin-proteasome pathway annotation in Diaphorina citri can reveal potential targets for RNAi-based pest management.

Ubiquitination is an ATP-dependent process that targets proteins for degradation by the proteasome. Here, we annotated 15 genes from the ubiquitin-proteasome pathway in the Asian citrus psyllid, Diaphorina citri. This psyllid vector has come to prominence in the last decade owing to its role in the transmission of the devastating bacterial pathogen, Candidatus Liberibacter asiaticus (CLas). Infection of citrus crops by this pathogen causes Huanglongbing (HLB), or citrus greening disease, and results in the eventual death of citrus trees. The identification and correct annotation of these genes in D. citri will be useful for functional genomic studies to aid the development of RNAi-based management strategies aimed at reducing the spread of HLB. Investigating the effects of CLas infection on the expression of ubiquitin-proteasome pathway genes may provide new information about the role these genes play in the acquisition and transmission of CLas by D. citri.

Annotation of glycolysis, gluconeogenesis, and trehaloneogenesis pathways provide insight into carbohydrate metabolism in the Asian citrus psyllid.

Citrus greening disease is caused by the pathogen Candidatus Liberibacter asiaticus and transmitted by the Asian citrus psyllid, Diaphorina citri. No curative treatment or significant prevention mechanism exists for this disease, which causes economic losses from reduced citrus production. A high-quality genome of D. citri is being manually annotated to provide accurate gene models to identify novel control targets and increase understanding of this pest. Here, we annotated 25 D. citri genes involved in glycolysis and gluconeogenesis, and seven in trehaloneogenesis. Comparative analysis showed that glycolysis genes in D. citri are highly conserved but copy numbers vary. Analysis of expression levels revealed upregulation of several enzymes in the glycolysis pathway in the thorax, consistent with the primary use of glucose by thoracic flight muscles. Manually annotating these core metabolic pathways provides accurate genomic foundation for developing gene-targeting therapeutics to control D. citri.