Comprehensive meta-analysis, co-expression, and miRNA nested network analysis identifies gene candidates in citrus against Huanglongbing disease
Identifieur interne : 000631 ( Pmc/Corpus ); précédent : 000630; suivant : 000632Comprehensive meta-analysis, co-expression, and miRNA nested network analysis identifies gene candidates in citrus against Huanglongbing disease
Auteurs : Nidhi Rawat ; Sandhya P. Kiran ; Dongliang Du ; Fred G. Gmitter ; Zhanao DengSource :
- BMC Plant Biology [ 1471-2229 ] ; 2015.
Abstract
Huanglongbing (HLB), the most devastating disease of citrus, is associated with infection by
Twenty-two publically available transcriptome studies on citrus–HLB interactions, comprising 18 susceptible (S) datasets and four resistant (R) datasets, were investigated using Limma and RankProd methods of meta-analysis. A combined list of 7,412 differentially expressed probe sets was generated using a Teradata in-house Structured Query Language (SQL) script. We identified the 65 most common probe sets modulated in HLB disease among different tissues from the S and R datasets. Gene ontology analysis of these probe sets suggested that carbohydrate metabolism, nutrient transport, and biotic stress were the core pathways that were modulated in citrus by
Most of the potential candidate hub probe sets were co-expressed with gibberellin pathway (GA)-related probe sets, implying the role of GA signalling in HLB resistance. Our findings contribute to the integration of existing citrus–HLB transcriptome data that will help to elucidate the holistic picture of the citrus–HLB interaction. The citrus probe sets identified in this analysis signify a robust set of HLB-responsive candidates that are useful for further validation.
The online version of this article (doi:10.1186/s12870-015-0568-4) contains supplementary material, which is available to authorized users.
Url:
DOI: 10.1186/s12870-015-0568-4
PubMed: 26215595
PubMed Central: 4517500
Links to Exploration step
PMC:4517500Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Comprehensive meta-analysis, co-expression, and miRNA nested network analysis identifies gene candidates in citrus against Huanglongbing disease</title><author><name sortKey="Rawat, Nidhi" sort="Rawat, Nidhi" uniqKey="Rawat N" first="Nidhi" last="Rawat">Nidhi Rawat</name><affiliation><nlm:aff id="Aff1">University of Florida, Institute of Food and Agricultural Sciences, Gulf Coast Research and Education Center, Wimauma, FL 33598 USA</nlm:aff></affiliation></author><author><name sortKey="Kiran, Sandhya P" sort="Kiran, Sandhya P" uniqKey="Kiran S" first="Sandhya P." last="Kiran">Sandhya P. Kiran</name><affiliation><nlm:aff id="Aff2">Ocimum BioSolutions, Banjara Hills Road No. 1, VI Floor Reliance Classic, Hyderabad, 500039 India</nlm:aff></affiliation></author><author><name sortKey="Du, Dongliang" sort="Du, Dongliang" uniqKey="Du D" first="Dongliang" last="Du">Dongliang Du</name><affiliation><nlm:aff id="Aff3">University of Florida, Institute of Food and Agricultural Sciences, Citrus Research and Education Center, Lake Alfred, FL 33850 USA</nlm:aff></affiliation></author><author><name sortKey="Gmitter, Fred G" sort="Gmitter, Fred G" uniqKey="Gmitter F" first="Fred G." last="Gmitter">Fred G. Gmitter</name><affiliation><nlm:aff id="Aff3">University of Florida, Institute of Food and Agricultural Sciences, Citrus Research and Education Center, Lake Alfred, FL 33850 USA</nlm:aff></affiliation></author><author><name sortKey="Deng, Zhanao" sort="Deng, Zhanao" uniqKey="Deng Z" first="Zhanao" last="Deng">Zhanao Deng</name><affiliation><nlm:aff id="Aff1">University of Florida, Institute of Food and Agricultural Sciences, Gulf Coast Research and Education Center, Wimauma, FL 33598 USA</nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">26215595</idno><idno type="pmc">4517500</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4517500</idno><idno type="RBID">PMC:4517500</idno><idno type="doi">10.1186/s12870-015-0568-4</idno><date when="2015">2015</date><idno type="wicri:Area/Pmc/Corpus">000631</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Comprehensive meta-analysis, co-expression, and miRNA nested network analysis identifies gene candidates in citrus against Huanglongbing disease</title><author><name sortKey="Rawat, Nidhi" sort="Rawat, Nidhi" uniqKey="Rawat N" first="Nidhi" last="Rawat">Nidhi Rawat</name><affiliation><nlm:aff id="Aff1">University of Florida, Institute of Food and Agricultural Sciences, Gulf Coast Research and Education Center, Wimauma, FL 33598 USA</nlm:aff></affiliation></author><author><name sortKey="Kiran, Sandhya P" sort="Kiran, Sandhya P" uniqKey="Kiran S" first="Sandhya P." last="Kiran">Sandhya P. Kiran</name><affiliation><nlm:aff id="Aff2">Ocimum BioSolutions, Banjara Hills Road No. 1, VI Floor Reliance Classic, Hyderabad, 500039 India</nlm:aff></affiliation></author><author><name sortKey="Du, Dongliang" sort="Du, Dongliang" uniqKey="Du D" first="Dongliang" last="Du">Dongliang Du</name><affiliation><nlm:aff id="Aff3">University of Florida, Institute of Food and Agricultural Sciences, Citrus Research and Education Center, Lake Alfred, FL 33850 USA</nlm:aff></affiliation></author><author><name sortKey="Gmitter, Fred G" sort="Gmitter, Fred G" uniqKey="Gmitter F" first="Fred G." last="Gmitter">Fred G. Gmitter</name><affiliation><nlm:aff id="Aff3">University of Florida, Institute of Food and Agricultural Sciences, Citrus Research and Education Center, Lake Alfred, FL 33850 USA</nlm:aff></affiliation></author><author><name sortKey="Deng, Zhanao" sort="Deng, Zhanao" uniqKey="Deng Z" first="Zhanao" last="Deng">Zhanao Deng</name><affiliation><nlm:aff id="Aff1">University of Florida, Institute of Food and Agricultural Sciences, Gulf Coast Research and Education Center, Wimauma, FL 33598 USA</nlm:aff></affiliation></author></analytic><series><title level="j">BMC Plant Biology</title><idno type="eISSN">1471-2229</idno><imprint><date when="2015">2015</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><sec><title>Background</title><p>Huanglongbing (HLB), the most devastating disease of citrus, is associated with infection by <italic>Candidatus</italic> Liberibacter asiaticus (<italic>Ca</italic>Las) and is vectored by the Asian citrus psyllid (ACP). Recently, the molecular basis of citrus–HLB interactions has been examined using transcriptome analyses, and these analyses have identified many probe sets and pathways modulated by <italic>Ca</italic>Las infection among different citrus cultivars. However, lack of consistency among reported findings indicates that an integrative approach is needed. This study was designed to identify the candidate probe sets in citrus–HLB interactions using meta-analysis and gene co-expression network modelling.</p></sec><sec><title>Results</title><p>Twenty-two publically available transcriptome studies on citrus–HLB interactions, comprising 18 susceptible (S) datasets and four resistant (R) datasets, were investigated using Limma and RankProd methods of meta-analysis. A combined list of 7,412 differentially expressed probe sets was generated using a Teradata in-house Structured Query Language (SQL) script. We identified the 65 most common probe sets modulated in HLB disease among different tissues from the S and R datasets. Gene ontology analysis of these probe sets suggested that carbohydrate metabolism, nutrient transport, and biotic stress were the core pathways that were modulated in citrus by <italic>Ca</italic>Las infection and HLB development. We also identified R-specific probe sets, which encoded leucine-rich repeat proteins, chitinase, constitutive disease resistance (CDR), miraculins, and lectins. Weighted gene co-expression network analysis (WGCNA) was conducted on 3,499 probe sets, and 21 modules with major hub probe sets were identified. Further, a miRNA nested network was created to examine gene regulation of the 3,499 target probe sets. Results suggest that csi-miR167 and csi-miR396 could affect ion transporters and defence response pathways, respectively.</p></sec><sec><title>Conclusion</title><p>Most of the potential candidate hub probe sets were co-expressed with gibberellin pathway (GA)-related probe sets, implying the role of GA signalling in HLB resistance. Our findings contribute to the integration of existing citrus–HLB transcriptome data that will help to elucidate the holistic picture of the citrus–HLB interaction. The citrus probe sets identified in this analysis signify a robust set of HLB-responsive candidates that are useful for further validation.</p></sec><sec><title>Electronic supplementary material</title><p>The online version of this article (doi:10.1186/s12870-015-0568-4) contains supplementary material, which is available to authorized users.</p></sec></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Albrecht, U" uniqKey="Albrecht U">U Albrecht</name></author><author><name sortKey="Bowman, Kd" uniqKey="Bowman K">KD Bowman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Albrecht, U" uniqKey="Albrecht U">U Albrecht</name></author><author><name sortKey="Bowman, Kd" uniqKey="Bowman K">KD Bowman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Aritua, V" uniqKey="Aritua V">V Aritua</name></author><author><name sortKey="Achor, D" uniqKey="Achor D">D Achor</name></author><author><name sortKey="Gmitter, Fg" uniqKey="Gmitter F">FG Gmitter</name></author><author><name sortKey="Albrigo, G" uniqKey="Albrigo G">G Albrigo</name></author><author><name sortKey="Wang, N" uniqKey="Wang N">N Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bhargava, A" uniqKey="Bhargava A">A Bhargava</name></author><author><name sortKey="Clabaugh, I" uniqKey="Clabaugh I">I Clabaugh</name></author><author><name sortKey="To, Jp" uniqKey="To J">JP To</name></author><author><name sortKey="Maxwell, Bb" uniqKey="Maxwell B">BB Maxwell</name></author><author><name sortKey="Chiang, Yh" uniqKey="Chiang Y">YH Chiang</name></author><author><name sortKey="Schaller, Ge" uniqKey="Schaller G">GE Schaller</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bove, Jm" uniqKey="Bove J">JM Bové</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Breitling, R" uniqKey="Breitling R">R Breitling</name></author><author><name sortKey="Amtmann, A" uniqKey="Amtmann A">A Amtmann</name></author><author><name sortKey="Herzyk, P" uniqKey="Herzyk P">P Herzyk</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Butte, Aj" uniqKey="Butte A">AJ Butte</name></author><author><name sortKey="Tamayo, P" uniqKey="Tamayo P">P Tamayo</name></author><author><name sortKey="Slonim, D" uniqKey="Slonim D">D Slonim</name></author><author><name sortKey="Golub, Tr" uniqKey="Golub T">TR Golub</name></author><author><name sortKey="Kohane, Is" uniqKey="Kohane I">IS Kohane</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Buttner, D" uniqKey="Buttner D">D Buttner</name></author><author><name sortKey="He, Sy" uniqKey="He S">SY He</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Canales, J" uniqKey="Canales J">J Canales</name></author><author><name sortKey="Moyano, Tc" uniqKey="Moyano T">TC Moyano</name></author><author><name sortKey="Villarroel, E" uniqKey="Villarroel E">E Villarroel</name></author><author><name sortKey="Gutierrez, Ra" uniqKey="Gutierrez R">RA Gutiérrez</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chang, C" uniqKey="Chang C">C Chang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Da Graca, J" uniqKey="Da Graca J">J Da Graca</name></author><author><name sortKey="Korsten, L" uniqKey="Korsten L">L Korsten</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ding, Y" uniqKey="Ding Y">Y Ding</name></author><author><name sortKey="Wu, W" uniqKey="Wu W">W Wu</name></author><author><name sortKey="Wei, W" uniqKey="Wei W">W Wei</name></author><author><name sortKey="Davis, R" uniqKey="Davis R">R Davis</name></author><author><name sortKey="Lee, I" uniqKey="Lee I">I Lee</name></author><author><name sortKey="Hammond, R" uniqKey="Hammond R">R Hammond</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Doncheva, Nt" uniqKey="Doncheva N">NT Doncheva</name></author><author><name sortKey="Assenov, Y" uniqKey="Assenov Y">Y Assenov</name></author><author><name sortKey="Domingues, Fs" uniqKey="Domingues F">FS Domingues</name></author><author><name sortKey="Albrecht, M" uniqKey="Albrecht M">M Albrecht</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Du, Z" uniqKey="Du Z">Z Du</name></author><author><name sortKey="Zhou, X" uniqKey="Zhou X">X Zhou</name></author><author><name sortKey="Ling, Y" uniqKey="Ling Y">Y Ling</name></author><author><name sortKey="Zhang, Z" uniqKey="Zhang Z">Z Zhang</name></author><author><name sortKey="Su, Z" uniqKey="Su Z">Z Su</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Duan, Y" uniqKey="Duan Y">Y Duan</name></author><author><name sortKey="Zhou, L" uniqKey="Zhou L">L Zhou</name></author><author><name sortKey="Hall, Dg" uniqKey="Hall D">DG Hall</name></author><author><name sortKey="Li, W" uniqKey="Li W">W Li</name></author><author><name sortKey="Doddapaneni, H" uniqKey="Doddapaneni H">H Doddapaneni</name></author><author><name sortKey="Lin, H" uniqKey="Lin H">H Lin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Eisen, Mb" uniqKey="Eisen M">MB Eisen</name></author><author><name sortKey="Spellman, Pt" uniqKey="Spellman P">PT Spellman</name></author><author><name sortKey="Brown, Po" uniqKey="Brown P">PO Brown</name></author><author><name sortKey="Botstein, D" uniqKey="Botstein D">D Botstein</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Fagen, Jr" uniqKey="Fagen J">JR Fagen</name></author><author><name sortKey="Leonard, Mt" uniqKey="Leonard M">MT Leonard</name></author><author><name sortKey="Mccullough, Cm" uniqKey="Mccullough C">CM McCullough</name></author><author><name sortKey="Edirisinghe, Jn" uniqKey="Edirisinghe J">JN Edirisinghe</name></author><author><name sortKey="Henry, Cs" uniqKey="Henry C">CS Henry</name></author><author><name sortKey="Davis, Mj" uniqKey="Davis M">MJ Davis</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Fan, J" uniqKey="Fan J">J Fan</name></author><author><name sortKey="Chen, C" uniqKey="Chen C">C Chen</name></author><author><name sortKey="Yu, Q" uniqKey="Yu Q">Q Yu</name></author><author><name sortKey="Brlansky, Rh" uniqKey="Brlansky R">RH Brlansky</name></author><author><name sortKey="Li, Z" uniqKey="Li Z">Z Li</name></author><author><name sortKey="Gmitter, Fg" uniqKey="Gmitter F">FG Gmitter</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Fan, J" uniqKey="Fan J">J Fan</name></author><author><name sortKey="Chen, C" uniqKey="Chen C">C Chen</name></author><author><name sortKey="Yu, Q" uniqKey="Yu Q">Q Yu</name></author><author><name sortKey="Khalaf, A" uniqKey="Khalaf A">A Khalaf</name></author><author><name sortKey="Achor, Ds" uniqKey="Achor D">DS Achor</name></author><author><name sortKey="Brlansky, Rh" uniqKey="Brlansky R">RH Brlansky</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Folimonova, Sy" uniqKey="Folimonova S">SY Folimonova</name></author><author><name sortKey="Robertson, Cj" uniqKey="Robertson C">CJ Robertson</name></author><author><name sortKey="Garnsey, Sm" uniqKey="Garnsey S">SM Garnsey</name></author><author><name sortKey="Gowda, S" uniqKey="Gowda S">S Gowda</name></author><author><name sortKey="Dawson, Wo" uniqKey="Dawson W">WO Dawson</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Gautier, L" uniqKey="Gautier L">L Gautier</name></author><author><name sortKey="Cope, L" uniqKey="Cope L">L Cope</name></author><author><name sortKey="Bolstad, Bm" uniqKey="Bolstad B">BM Bolstad</name></author><author><name sortKey="Irizarry, Ra" uniqKey="Irizarry R">RA Irizarry</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gottwald, Tr" uniqKey="Gottwald T">TR Gottwald</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gottwald, T" uniqKey="Gottwald T">T Gottwald</name></author><author><name sortKey="Bassanezi, R" uniqKey="Bassanezi R">R Bassanezi</name></author><author><name sortKey="Amorim, L" uniqKey="Amorim L">L Amorim</name></author><author><name sortKey="Bergamin Filho, A" uniqKey="Bergamin Filho A">A Bergamin-Filho</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Halbert, Se" uniqKey="Halbert S">SE Halbert</name></author><author><name sortKey="Manjunath, Kl" uniqKey="Manjunath K">KL Manjunath</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hammond, Jp" uniqKey="Hammond J">JP Hammond</name></author><author><name sortKey="White, Pj" uniqKey="White P">PJ White</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hodges, Aw" uniqKey="Hodges A">AW Hodges</name></author><author><name sortKey="Spreen, Th" uniqKey="Spreen T">TH Spreen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Holton, N" uniqKey="Holton N">N Holton</name></author><author><name sortKey="Nekrasov, V" uniqKey="Nekrasov V">V Nekrasov</name></author><author><name sortKey="Ronald, Pc" uniqKey="Ronald P">PC Ronald</name></author><author><name sortKey="Zipfel, C" uniqKey="Zipfel C">C Zipfel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hong, F" uniqKey="Hong F">F Hong</name></author><author><name sortKey="Breitling, R" uniqKey="Breitling R">R Breitling</name></author><author><name sortKey="Mcentee, Cw" uniqKey="Mcentee C">CW McEntee</name></author><author><name sortKey="Wittner, Bs" uniqKey="Wittner B">BS Wittner</name></author><author><name sortKey="Nemhauser, Jl" uniqKey="Nemhauser J">JL Nemhauser</name></author><author><name sortKey="Chory, J" uniqKey="Chory J">J Chory</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jagoueix, S" uniqKey="Jagoueix S">S Jagoueix</name></author><author><name sortKey="Bove, Jm" uniqKey="Bove J">JM Bove</name></author><author><name sortKey="Garnier, M" uniqKey="Garnier M">M Garnier</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kim, J" uniqKey="Kim J">J Kim</name></author><author><name sortKey="Sagaram, Us" uniqKey="Sagaram U">US Sagaram</name></author><author><name sortKey="Burns, Jk" uniqKey="Burns J">JK Burns</name></author><author><name sortKey="Li, J" uniqKey="Li J">J Li</name></author><author><name sortKey="Wang, N" uniqKey="Wang N">N Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Landi, L" uniqKey="Landi L">L Landi</name></author><author><name sortKey="Romanazzi, G" uniqKey="Romanazzi G">G Romanazzi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Langfelder, P" uniqKey="Langfelder P">P Langfelder</name></author><author><name sortKey="Horvath, S" uniqKey="Horvath S">S Horvath</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liao, Hl" uniqKey="Liao H">HL Liao</name></author><author><name sortKey="Burns, Jk" uniqKey="Burns J">JK Burns</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liu, Y" uniqKey="Liu Y">Y Liu</name></author><author><name sortKey="Heying, E" uniqKey="Heying E">E Heying</name></author><author><name sortKey="Tanumihardjo, Sa" uniqKey="Tanumihardjo S">SA Tanumihardjo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lopes, S" uniqKey="Lopes S">S Lopes</name></author><author><name sortKey="Bertolini, E" uniqKey="Bertolini E">E Bertolini</name></author><author><name sortKey="Frare, G" uniqKey="Frare G">G Frare</name></author><author><name sortKey="Martins, E" uniqKey="Martins E">E Martins</name></author><author><name sortKey="Wulff, N" uniqKey="Wulff N">N Wulff</name></author><author><name sortKey="Teixeira, D" uniqKey="Teixeira D">D Teixeira</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lu, Ca" uniqKey="Lu C">CA Lu</name></author><author><name sortKey="Ho, Th" uniqKey="Ho T">TH Ho</name></author><author><name sortKey="Ho, Sl" uniqKey="Ho S">SL Ho</name></author><author><name sortKey="Yu, Sm" uniqKey="Yu S">SM Yu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lu, Yt" uniqKey="Lu Y">YT Lu</name></author><author><name sortKey="Li, My" uniqKey="Li M">MY Li</name></author><author><name sortKey="Cheng, Kt" uniqKey="Cheng K">KT Cheng</name></author><author><name sortKey="Tan, Cm" uniqKey="Tan C">CM Tan</name></author><author><name sortKey="Su, Lw" uniqKey="Su L">LW Su</name></author><author><name sortKey="Lin, Wy" uniqKey="Lin W">WY Lin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mafra, V" uniqKey="Mafra V">V Mafra</name></author><author><name sortKey="Martins, Pk" uniqKey="Martins P">PK Martins</name></author><author><name sortKey="Francisco, Cs" uniqKey="Francisco C">CS Francisco</name></author><author><name sortKey="Ribeiro Alves, M" uniqKey="Ribeiro Alves M">M Ribeiro-Alves</name></author><author><name sortKey="Freitas Astua, J" uniqKey="Freitas Astua J">J Freitas-Astua</name></author><author><name sortKey="Machado, Ma" uniqKey="Machado M">MA Machado</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Martinelli, F" uniqKey="Martinelli F">F Martinelli</name></author><author><name sortKey="Reagan, Rl" uniqKey="Reagan R">RL Reagan</name></author><author><name sortKey="Uratsu, Sl" uniqKey="Uratsu S">SL Uratsu</name></author><author><name sortKey="Phu, Ml" uniqKey="Phu M">ML Phu</name></author><author><name sortKey="Albrecht, U" uniqKey="Albrecht U">U Albrecht</name></author><author><name sortKey="Zhao, W" uniqKey="Zhao W">W Zhao</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Martinelli, F" uniqKey="Martinelli F">F Martinelli</name></author><author><name sortKey="Uratsu, Sl" uniqKey="Uratsu S">SL Uratsu</name></author><author><name sortKey="Albrecht, U" uniqKey="Albrecht U">U Albrecht</name></author><author><name sortKey="Reagan, Rl" uniqKey="Reagan R">RL Reagan</name></author><author><name sortKey="Phu, Ml" uniqKey="Phu M">ML Phu</name></author><author><name sortKey="Britton, M" uniqKey="Britton M">M Britton</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nazri, A" uniqKey="Nazri A">A Nazri</name></author><author><name sortKey="Lio, P" uniqKey="Lio P">P Lio</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nwugo, Cc" uniqKey="Nwugo C">CC Nwugo</name></author><author><name sortKey="Duan, Y" uniqKey="Duan Y">Y Duan</name></author><author><name sortKey="Lin, H" uniqKey="Lin H">H Lin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Peumans, Wj" uniqKey="Peumans W">WJ Peumans</name></author><author><name sortKey="Van Damme, Ej" uniqKey="Van Damme E">EJ Van Damme</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Plank, J" uniqKey="Plank J">J Plank</name></author><author><name sortKey="Siebenhofer, A" uniqKey="Siebenhofer A">A Siebenhofer</name></author><author><name sortKey="Berghold, A" uniqKey="Berghold A">A Berghold</name></author><author><name sortKey="Jeitler, K" uniqKey="Jeitler K">K Jeitler</name></author><author><name sortKey="Horvath, K" uniqKey="Horvath K">K Horvath</name></author><author><name sortKey="Mrak, P" uniqKey="Mrak P">P Mrak</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Puniya, Bl" uniqKey="Puniya B">BL Puniya</name></author><author><name sortKey="Kulshreshtha, D" uniqKey="Kulshreshtha D">D Kulshreshtha</name></author><author><name sortKey="Verma, Sp" uniqKey="Verma S">SP Verma</name></author><author><name sortKey="Kumar, S" uniqKey="Kumar S">S Kumar</name></author><author><name sortKey="Ramachandran, S" uniqKey="Ramachandran S">S Ramachandran</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pustika, A" uniqKey="Pustika A">A Pustika</name></author><author><name sortKey="Subandiyah, S" uniqKey="Subandiyah S">S Subandiyah</name></author><author><name sortKey="Holford, P" uniqKey="Holford P">P Holford</name></author><author><name sortKey="Beattie, G" uniqKey="Beattie G">G Beattie</name></author><author><name sortKey="Iwanami, T" uniqKey="Iwanami T">T Iwanami</name></author><author><name sortKey="Masaoka, Y" uniqKey="Masaoka Y">Y Masaoka</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rawat, N" uniqKey="Rawat N">N Rawat</name></author><author><name sortKey="Neeraja, Cn" uniqKey="Neeraja C">CN Neeraja</name></author><author><name sortKey="Nair, S" uniqKey="Nair S">S Nair</name></author><author><name sortKey="Bentur, Js" uniqKey="Bentur J">JS Bentur</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sarwar, M" uniqKey="Sarwar M">M Sarwar</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shaik, R" uniqKey="Shaik R">R Shaik</name></author><author><name sortKey="Ramakrishna, W" uniqKey="Ramakrishna W">W Ramakrishna</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Smyth, Gk" uniqKey="Smyth G">GK Smyth</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sugio, A" uniqKey="Sugio A">A Sugio</name></author><author><name sortKey="Maclean, Am" uniqKey="Maclean A">AM MacLean</name></author><author><name sortKey="Hogenhout, Sa" uniqKey="Hogenhout S">SA Hogenhout</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Taheri, F" uniqKey="Taheri F">F Taheri</name></author><author><name sortKey="Nematzadeh, G" uniqKey="Nematzadeh G">G Nematzadeh</name></author><author><name sortKey="Zamharir, Mg" uniqKey="Zamharir M">MG Zamharir</name></author><author><name sortKey="Nekouei, Mk" uniqKey="Nekouei M">MK Nekouei</name></author><author><name sortKey="Naghavi, M" uniqKey="Naghavi M">M Naghavi</name></author><author><name sortKey="Mardi, M" uniqKey="Mardi M">M Mardi</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Thirunavukkarasu, N" uniqKey="Thirunavukkarasu N">N Thirunavukkarasu</name></author><author><name sortKey="Hossain, F" uniqKey="Hossain F">F Hossain</name></author><author><name sortKey="Mohan, S" uniqKey="Mohan S">S Mohan</name></author><author><name sortKey="Shiriga, K" uniqKey="Shiriga K">K Shiriga</name></author><author><name sortKey="Mittal, S" uniqKey="Mittal S">S Mittal</name></author><author><name sortKey="Sharma, R" uniqKey="Sharma R">R Sharma</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ting, Sv" uniqKey="Ting S">SV Ting</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tsukuda, S" uniqKey="Tsukuda S">S Tsukuda</name></author><author><name sortKey="Gomi, K" uniqKey="Gomi K">K Gomi</name></author><author><name sortKey="Yamamoto, H" uniqKey="Yamamoto H">H Yamamoto</name></author><author><name sortKey="Akimitsu, K" uniqKey="Akimitsu K">K Akimitsu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Usadel, B" uniqKey="Usadel B">B Usadel</name></author><author><name sortKey="Nagel, A" uniqKey="Nagel A">A Nagel</name></author><author><name sortKey="Steinhauser, D" uniqKey="Steinhauser D">D Steinhauser</name></author><author><name sortKey="Gibon, Y" uniqKey="Gibon Y">Y Gibon</name></author><author><name sortKey="Blasing, Oe" uniqKey="Blasing O">OE Blasing</name></author><author><name sortKey="Redestig, H" uniqKey="Redestig H">H Redestig</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Usadel, B" uniqKey="Usadel B">B Usadel</name></author><author><name sortKey="Nagel, A" uniqKey="Nagel A">A Nagel</name></author><author><name sortKey="Thimm, O" uniqKey="Thimm O">O Thimm</name></author><author><name sortKey="Redestig, H" uniqKey="Redestig H">H Redestig</name></author><author><name sortKey="Blaesing, Oe" uniqKey="Blaesing O">OE Blaesing</name></author><author><name sortKey="Palacios Rojas, N" uniqKey="Palacios Rojas N">N Palacios-Rojas</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Van Loon, L" uniqKey="Van Loon L">L Van Loon</name></author><author><name sortKey="Bakker, P" uniqKey="Bakker P">P Bakker</name></author><author><name sortKey="Pieterse, C" uniqKey="Pieterse C">C Pieterse</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, N" uniqKey="Wang N">N Wang</name></author><author><name sortKey="Trivedi, P" uniqKey="Trivedi P">P Trivedi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, M" uniqKey="Wang M">M Wang</name></author><author><name sortKey="Zheng, Q" uniqKey="Zheng Q">Q Zheng</name></author><author><name sortKey="Shen, Q" uniqKey="Shen Q">Q Shen</name></author><author><name sortKey="Guo, S" uniqKey="Guo S">S Guo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Weintraub, Pg" uniqKey="Weintraub P">PG Weintraub</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Xia, Y" uniqKey="Xia Y">Y Xia</name></author><author><name sortKey="Suzuki, H" uniqKey="Suzuki H">H Suzuki</name></author><author><name sortKey="Borevitz, J" uniqKey="Borevitz J">J Borevitz</name></author><author><name sortKey="Blount, J" uniqKey="Blount J">J Blount</name></author><author><name sortKey="Guo, Z" uniqKey="Guo Z">Z Guo</name></author><author><name sortKey="Patel, K" uniqKey="Patel K">K Patel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yip, Am" uniqKey="Yip A">AM Yip</name></author><author><name sortKey="Horvath, S" uniqKey="Horvath S">S Horvath</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, B" uniqKey="Zhang B">B Zhang</name></author><author><name sortKey="Horvath, S" uniqKey="Horvath S">S Horvath</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, Z" uniqKey="Zhang Z">Z Zhang</name></author><author><name sortKey="Yu, J" uniqKey="Yu J">J Yu</name></author><author><name sortKey="Li, D" uniqKey="Li D">D Li</name></author><author><name sortKey="Zhang, Z" uniqKey="Zhang Z">Z Zhang</name></author><author><name sortKey="Liu, F" uniqKey="Liu F">F Liu</name></author><author><name sortKey="Zhou, X" uniqKey="Zhou X">X Zhou</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhao, H" uniqKey="Zhao H">H Zhao</name></author><author><name sortKey="Sun, R" uniqKey="Sun R">R Sun</name></author><author><name sortKey="Albrecht, U" uniqKey="Albrecht U">U Albrecht</name></author><author><name sortKey="Padmanabhan, C" uniqKey="Padmanabhan C">C Padmanabhan</name></author><author><name sortKey="Wang, A" uniqKey="Wang A">A Wang</name></author><author><name sortKey="Coffey, Md" uniqKey="Coffey M">MD Coffey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zheng, Zl" uniqKey="Zheng Z">ZL Zheng</name></author><author><name sortKey="Zhao, Y" uniqKey="Zhao Y">Y Zhao</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhong, R" uniqKey="Zhong R">R Zhong</name></author><author><name sortKey="Demura, T" uniqKey="Demura T">T Demura</name></author><author><name sortKey="Ye, Zh" uniqKey="Ye Z">ZH Ye</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zou, H" uniqKey="Zou H">H Zou</name></author><author><name sortKey="Gowda, S" uniqKey="Gowda S">S Gowda</name></author><author><name sortKey="Zhou, L" uniqKey="Zhou L">L Zhou</name></author><author><name sortKey="Hajeri, S" uniqKey="Hajeri S">S Hajeri</name></author><author><name sortKey="Chen, G" uniqKey="Chen G">G Chen</name></author><author><name sortKey="Duan, Y" uniqKey="Duan Y">Y Duan</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">BMC Plant Biol</journal-id><journal-id journal-id-type="iso-abbrev">BMC Plant Biol</journal-id><journal-title-group><journal-title>BMC Plant Biology</journal-title></journal-title-group><issn pub-type="epub">1471-2229</issn><publisher><publisher-name>BioMed Central</publisher-name><publisher-loc>London</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">26215595</article-id><article-id pub-id-type="pmc">4517500</article-id><article-id pub-id-type="publisher-id">568</article-id><article-id pub-id-type="doi">10.1186/s12870-015-0568-4</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group></article-categories><title-group><article-title>Comprehensive meta-analysis, co-expression, and miRNA nested network analysis identifies gene candidates in citrus against Huanglongbing disease</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Rawat</surname><given-names>Nidhi</given-names></name><address><email>nrawat@ufl.edu</email></address><xref ref-type="aff" rid="Aff1"></xref></contrib><contrib contrib-type="author"><name><surname>Kiran</surname><given-names>Sandhya P.</given-names></name><address><email>p_sandhyakiran@yahoo.com</email></address><xref ref-type="aff" rid="Aff2"></xref></contrib><contrib contrib-type="author"><name><surname>Du</surname><given-names>Dongliang</given-names></name><address><email>dldu@ufl.edu</email></address><xref ref-type="aff" rid="Aff3"></xref></contrib><contrib contrib-type="author"><name><surname>Gmitter</surname><given-names>Fred G.</given-names><suffix>Jr</suffix></name><address><email>fgmitter@ufl.edu</email></address><xref ref-type="aff" rid="Aff3"></xref></contrib><contrib contrib-type="author" corresp="yes"><name><surname>Deng</surname><given-names>Zhanao</given-names></name><address><email>zdeng@ufl.edu</email></address><xref ref-type="aff" rid="Aff1"></xref></contrib><aff id="Aff1"><label></label>University of Florida, Institute of Food and Agricultural Sciences, Gulf Coast Research and Education Center, Wimauma, FL 33598 USA</aff><aff id="Aff2"><label></label>Ocimum BioSolutions, Banjara Hills Road No. 1, VI Floor Reliance Classic, Hyderabad, 500039 India</aff><aff id="Aff3"><label></label>University of Florida, Institute of Food and Agricultural Sciences, Citrus Research and Education Center, Lake Alfred, FL 33850 USA</aff></contrib-group><pub-date pub-type="epub"><day>28</day><month>7</month><year>2015</year></pub-date><pub-date pub-type="pmc-release"><day>28</day><month>7</month><year>2015</year></pub-date><pub-date pub-type="collection"><year>2015</year></pub-date><volume>15</volume><elocation-id>184</elocation-id><history><date date-type="received"><day>9</day><month>4</month><year>2015</year></date><date date-type="accepted"><day>7</day><month>7</month><year>2015</year></date></history><permissions><copyright-statement>© Rawat et al. 2015</copyright-statement><license license-type="open-access"><license-p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0">http://creativecommons.org/licenses/by/4.0</ext-link>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/publicdomain/zero/1.0/">http://creativecommons.org/publicdomain/zero/1.0/</ext-link>) applies to the data made available in this article, unless otherwise stated.</license-p></license></permissions><abstract id="Abs1"><sec><title>Background</title><p>Huanglongbing (HLB), the most devastating disease of citrus, is associated with infection by <italic>Candidatus</italic> Liberibacter asiaticus (<italic>Ca</italic>Las) and is vectored by the Asian citrus psyllid (ACP). Recently, the molecular basis of citrus–HLB interactions has been examined using transcriptome analyses, and these analyses have identified many probe sets and pathways modulated by <italic>Ca</italic>Las infection among different citrus cultivars. However, lack of consistency among reported findings indicates that an integrative approach is needed. This study was designed to identify the candidate probe sets in citrus–HLB interactions using meta-analysis and gene co-expression network modelling.</p></sec><sec><title>Results</title><p>Twenty-two publically available transcriptome studies on citrus–HLB interactions, comprising 18 susceptible (S) datasets and four resistant (R) datasets, were investigated using Limma and RankProd methods of meta-analysis. A combined list of 7,412 differentially expressed probe sets was generated using a Teradata in-house Structured Query Language (SQL) script. We identified the 65 most common probe sets modulated in HLB disease among different tissues from the S and R datasets. Gene ontology analysis of these probe sets suggested that carbohydrate metabolism, nutrient transport, and biotic stress were the core pathways that were modulated in citrus by <italic>Ca</italic>Las infection and HLB development. We also identified R-specific probe sets, which encoded leucine-rich repeat proteins, chitinase, constitutive disease resistance (CDR), miraculins, and lectins. Weighted gene co-expression network analysis (WGCNA) was conducted on 3,499 probe sets, and 21 modules with major hub probe sets were identified. Further, a miRNA nested network was created to examine gene regulation of the 3,499 target probe sets. Results suggest that csi-miR167 and csi-miR396 could affect ion transporters and defence response pathways, respectively.</p></sec><sec><title>Conclusion</title><p>Most of the potential candidate hub probe sets were co-expressed with gibberellin pathway (GA)-related probe sets, implying the role of GA signalling in HLB resistance. Our findings contribute to the integration of existing citrus–HLB transcriptome data that will help to elucidate the holistic picture of the citrus–HLB interaction. The citrus probe sets identified in this analysis signify a robust set of HLB-responsive candidates that are useful for further validation.</p></sec><sec><title>Electronic supplementary material</title><p>The online version of this article (doi:10.1186/s12870-015-0568-4) contains supplementary material, which is available to authorized users.</p></sec></abstract><kwd-group xml:lang="en"><title>Keywords</title><kwd>Citrus</kwd><kwd>Gene co-expression analysis</kwd><kwd>Gene ontology</kwd><kwd>Huanglongbing</kwd><kwd>HLB resistance</kwd><kwd>Meta-analysis</kwd><kwd>miRNA network analysis</kwd><kwd>Common and R-specific probe sets</kwd></kwd-group><custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name><meta-value>© The Author(s) 2015</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec id="Sec1" sec-type="introduction"><title>Background</title><p>Citrus are among the most popular fruit crops in the world, providing energy (carbohydrates) and nutrients, and are an important component of the daily diet in many parts of the world [<xref ref-type="bibr" rid="CR35">35</xref>]. Fresh citrus is also a good source of dietary fibre and vitamins B and C. Citrus fruits and fruit products are globally important from nutritional and economic perspectives [<xref ref-type="bibr" rid="CR56">56</xref>]. However, there are various biotic and abiotic challenges to citrus production, among which Huanglongbing (HLB), or citrus greening disease, is the most devastating [<xref ref-type="bibr" rid="CR5">5</xref>, <xref ref-type="bibr" rid="CR25">25</xref>, <xref ref-type="bibr" rid="CR61">61</xref>]. HLB was first reported in China in the early 1900s and is now well established in many citrus-producing regions, including India, China, the United States, Indonesia, the Philippines, the Arabian Peninsula, Brazil, and Africa [<xref ref-type="bibr" rid="CR24">24</xref>, <xref ref-type="bibr" rid="CR25">25</xref>]. Brazil and the United States produce more than 90 % of the world’s supply of orange juice, and HLB is a current threat to the U.S. and the Brazilian citrus industry. HLB was first found in 2005 in Florida, the second-largest orange producing region in the world. Since then, HLB has reached epidemic proportions in Florida and has caused more than $4 billion in economic losses between 2005 and 2011 [<xref ref-type="bibr" rid="CR27">27</xref>]. HLB attacks all important commercial citrus, including oranges, grapefruit, and tangerines [<xref ref-type="bibr" rid="CR5">5</xref>, <xref ref-type="bibr" rid="CR20">20</xref>]. Sweet oranges and mandarins are considered highly susceptible, and sour oranges and grapefruits are moderately susceptible. It seems that some lemons are tolerant to HLB and some trifoliate orange (a close relative of citrus) are resistant of HLB disease [<xref ref-type="bibr" rid="CR20">20</xref>].</p><p>HLB is a bacterial disease caused by gram negative and phloem-restricted Liberibacter species, which are vectored by the Asian citrus psyllid (<italic>Diaphorina citri</italic> Kuwayama) [<xref ref-type="bibr" rid="CR21">21</xref>, <xref ref-type="bibr" rid="CR30">30</xref>]. <italic>D. citri</italic> feeds on new leaf growth, rendering twisted and curled leaves [<xref ref-type="bibr" rid="CR36">36</xref>]. The disease can also be transmitted to healthy trees by grafting of diseased budwood [<xref ref-type="bibr" rid="CR36">36</xref>]. Among the three known liberibacters that cause HLB disease, <italic>Candidatus</italic> Liberibacter asiaticus (<italic>Ca</italic>Las) is the most widespread. The other two species are more geographically constrained; <italic>Ca.</italic> L. africanus is present primarily in Africa [<xref ref-type="bibr" rid="CR30">30</xref>], and <italic>Ca.</italic> L. americanus has only been found in Brazil and China [<xref ref-type="bibr" rid="CR54">54</xref>]. The bacterium resides in phloem tissues and causes phloem collapse, which leads to decreased productivity. The HLB disease causes a rapid tree decline with blotchy mottling of leaves, and small, misshapen, irregularly coloured, bitter fruit with aborted seeds [<xref ref-type="bibr" rid="CR11">11</xref>, <xref ref-type="bibr" rid="CR23">23</xref>]. At present, there are no effective control methods for HLB, except for the use of HLB-free budwood for plant propagation, removal of infected trees to minimize inoculum, and insect vector control [<xref ref-type="bibr" rid="CR5">5</xref>]. The full genome sequencing of <italic>Ca</italic>Las (1.23 Mb) has made genome-based identification of virulence factors associated with HLB symptoms possible [<xref ref-type="bibr" rid="CR15">15</xref>]. However, <italic>Ca</italic>Las has not been cultured and a full understanding of the molecular basis of citrus–HLB interactions is lacking [<xref ref-type="bibr" rid="CR17">17</xref>]. It is difficult to create HLB-tolerant citrus cultivars via conventional breeding because of a lack of any known resistance (R) genes against HLB and because of the complex biology of the citrus host. However, examination of citrus genomics through transcriptome and proteome analysis can elucidate the differentially expressed genes/proteins as potential candidate to study citrus–HLB interactions.</p><p>Several transcriptomic studies of citrus–HLB interactions have been published so far [<xref ref-type="bibr" rid="CR1">1</xref>–<xref ref-type="bibr" rid="CR3">3</xref>, <xref ref-type="bibr" rid="CR18">18</xref>, <xref ref-type="bibr" rid="CR19">19</xref>, <xref ref-type="bibr" rid="CR31">31</xref>, <xref ref-type="bibr" rid="CR34">34</xref>, <xref ref-type="bibr" rid="CR39">39</xref>–<xref ref-type="bibr" rid="CR41">41</xref>]. These studies showed that the expression of genes in carbohydrate metabolism, nutrient transport, cell wall synthesis, and defence- and hormone-response pathways were reprogrammed by <italic>Ca</italic>Las infection and HLB disease development. These studies also suggested that disruption in carbohydrate source and sink pathway and phloem plugging were likely the main cause of HLB symptoms [<xref ref-type="bibr" rid="CR34">34</xref>, <xref ref-type="bibr" rid="CR39">39</xref>, <xref ref-type="bibr" rid="CR40">40</xref>]. Numerous differentially expressed probe sets have been identified in these studies, but there was only a limited level of overlapping among the studies in terms of differentially expressed gene probe sets, possibly because of technique-, tissue-, or genotype-specific effects. Furthermore, some studies focused only on fully symptomatic plants, while others followed a time course of events during HLB development. Collective interrogation of this collection of transcriptomic studies should provide a better understanding than could be gained from analysis of individual studies. However, there are significant challenges to meta-analyses with respect to data comparability, normalization, and analysis tools. Meta-analysis of microarrays has been used extensively in animal systems to define robust, regulated probe sets [<xref ref-type="bibr" rid="CR42">42</xref>, <xref ref-type="bibr" rid="CR45">45</xref>]. Recently, meta-analyses have also been used to identify differentially expressed probe sets in plants [<xref ref-type="bibr" rid="CR4">4</xref>, <xref ref-type="bibr" rid="CR9">9</xref>, <xref ref-type="bibr" rid="CR50">50</xref>, <xref ref-type="bibr" rid="CR55">55</xref>]. A recent study presented a meta-analysis of HLB-responsive probe sets in citrus, but this study used a small number of datasets (six), and only one R genotype (US-897) was included for comparison [<xref ref-type="bibr" rid="CR69">69</xref>]. To fill this gap, we performed a comprehensive meta-analysis using a larger dataset (22 studies) that included tissue-specific, susceptible, and resistance-specific HLB-responsive probe sets, and we further analysed these probe sets using co-expression and miRNA nested network approaches.</p><p>Our objective was to investigate the citrus–HLB interactions in detail using available transcriptome data and to identify robust probe sets in citrus that are regulated by <italic>Ca</italic>Las infection and HLB development. To do this, we analysed 46 publicly available citrus–HLB transcriptome microarray datasets from the Affymetrix GeneChip Citrus Genome Array using the Limma and RankProd methods. We identified the most statistically important tissue-specific probe sets and pathways, which we further analysed for enrichment of Gene Ontology (GO) terms using MapMan and AgriGO. In addition, we identified common probe sets that provided information on core citrus pathways and resistance-specific probe sets that gave clues about defence responses in citrus against <italic>Ca</italic>Las infection. Correlations among common probe sets identified with both Limma and RankProd were examined using WGCNA co-expression network analysis. We identified potential hub probe sets that showed the greatest number of interconnections and contributed most to citrus–HLB interactions. We also characterised the miRNA-mediated nested network in response to <italic>Ca</italic>Las infection. These findings provide potential candidate probe sets that can be used to dissect citrus–HLB interactions in detail.</p></sec><sec id="Sec2" sec-type="results"><title>Results</title><sec id="Sec3"><title>Meta-analysis of HLB-responsive transcriptome data</title><p>We used 22 datasets, including 18 datasets from susceptible citrus plants (S datasets, one each for root and stem tissue, eight for leaf tissues, and eight for fruit) and four from leaf tissues of HLB-resistant citrus plants (R datasets) (Table <xref rid="Tab1" ref-type="table">1</xref>). We retrieved the Gene Expression Omnibus (GEO) CEL files for the above-mentioned datasets that include 46 array samples and performed Robust MultiChip Analysis (RMA) normalisation using the affy package in the R program. Differentially expressed probe sets were identified using two statistical approaches, Limma and RankProd, in the Bioconductor package. In the Limma method, probe sets were identified based on a moderated <italic>t</italic> statistic, <italic>P</italic> < 0.05 and a false discovery rate (FDR) < 0.1. We did not use FDR < 0.05 because we could have missed many differentially expressed probe sets from some datasets (especially those for roots). The RankProd method was used under more stringent conditions (FDR < 0.01) to detect probe sets at a higher level of statistical significance. Hence, fewer probe sets were identified with RankProd than with Limma. More probe sets were identified with the Limma method also because we added common probe sets from Fan et al. [<xref ref-type="bibr" rid="CR19">19</xref>] and fruit data [<xref ref-type="bibr" rid="CR34">34</xref>, <xref ref-type="bibr" rid="CR40">40</xref>] using an in-house Teradata SQL script. We identified 7,412 probe sets using Limma and 4,221 using RankProd (Additional file <xref rid="MOESM1" ref-type="media">1</xref>). The above-mentioned studies reported a total of 7,326 differentially expressed probe sets. Of these, 6,634 were common in the Limma results and 3,184 were common in RankProd results (Fig. <xref rid="Fig1" ref-type="fig">1</xref>); 3,499 probe sets were differentially expressed in both Limma and RankProd.<table-wrap id="Tab1"><label>Table 1</label><caption><p>Details of microarray studies used for meta-analysis</p></caption><table frame="hsides" rules="groups"><thead><tr><th>Dataset</th><th>Designated in present analysis</th><th>Reference</th><th>Genotype</th><th>Number of samples</th><th>Tissue</th><th>Time after <italic>Ca</italic>Las infection</th><th>GSE Dataset</th><th>Analysis method</th><th>Analysis statistic</th><th>Up-regulated Probesets</th><th>Down- regulated Probesets</th></tr></thead><tbody><tr><td>Root dataset</td><td>Valencia root</td><td>[<xref ref-type="bibr" rid="CR3">3</xref>]</td><td>Valencia</td><td>6 (3C<sup>f</sup> + 3I<sup>f</sup>)</td><td>Root</td><td>16 months</td><td>GSE33004</td><td>Limma R Package</td><td>Unadjusted <italic>P</italic> < 0.05 FC > 2</td><td>56</td><td>55</td></tr><tr><td>Stem dataset</td><td>Valencia stem</td><td>[<xref ref-type="bibr" rid="CR3">3</xref>]</td><td>Valencia</td><td>6 (3C + 3I)</td><td>Stem</td><td>16 months</td><td>GSE33004</td><td>Limma R Package</td><td>Unadjusted <italic>P</italic> < 0.05 FC > 2</td><td>551</td><td>334</td></tr><tr><td rowspan="8">Susceptible (S) leaf dataset</td><td>Valencia 5–9 wai<sup>a</sup></td><td>[<xref ref-type="bibr" rid="CR1">1</xref>]</td><td>Valencia</td><td>5 (2C + 3I)</td><td>Leaf</td><td>5-9 week</td><td>GSE33459</td><td>local pooled Error (LPE)</td><td>Adjusted <italic>P</italic> (FDR < 0.05) FC > 2</td><td>195</td><td>84</td></tr><tr><td>Valencia 13–17 wai</td><td>[<xref ref-type="bibr" rid="CR1">1</xref>]</td><td>Valencia</td><td>5 (2C + 3I)</td><td>Leaf</td><td>13-17 week</td><td>GSE33459</td><td>local pooled Error (LPE)</td><td>Adjusted <italic>P</italic> (FDR < 0.05) FC > 2</td><td>310</td><td>205</td></tr><tr><td>Sweet orange</td><td>[<xref ref-type="bibr" rid="CR31">31</xref>]</td><td>sweet orange</td><td>6 (3C + 3I)</td><td>Leaf</td><td>8 months</td><td>GSE33003</td><td>Limma R Package</td><td>Unadjusted <italic>P</italic> < 0.05 FC > 2</td><td>307</td><td>317</td></tr><tr><td>Madam Vinous</td><td>[<xref ref-type="bibr" rid="CR18">18</xref>]</td><td>Madam Vinous</td><td>6 (3C + 3I)</td><td>Leaf</td><td>7 months</td><td>GSE29633</td><td>Limma R Package</td><td>Unadjusted <italic>P</italic> < 0.05 FC > 2</td><td>943</td><td>1060</td></tr><tr><td>Cleopatra</td><td>[<xref ref-type="bibr" rid="CR2">2</xref>]</td><td>Cleopatra</td><td>6 (3C + 3I)</td><td>Leaf</td><td>32 week</td><td>GSE30502</td><td>Partek suit two way ANOVA</td><td>Adjusted <italic>P</italic> (FDR < 0.05) FC > 4</td><td>326</td><td>68</td></tr><tr><td>Madam Vinous 5 wai</td><td>[<xref ref-type="bibr" rid="CR19">19</xref>]</td><td>Madam Vinous</td><td>6 (3C + 3I)</td><td>Leaf</td><td>5 week</td><td>Not available</td><td>Limma R Package</td><td>Unadjusted <italic>P</italic> < 0.05 FC > 2</td><td>83</td><td>101</td></tr><tr><td>Madam Vinous 17 wai</td><td>[<xref ref-type="bibr" rid="CR19">19</xref>]</td><td>Madam Vinous</td><td>6 (3C + 3I)</td><td>Leaf</td><td>17 week</td><td>Not available</td><td>Limma R Package</td><td>Unadjusted <italic>P</italic> < 0.05 FC > 2</td><td>614</td><td>414</td></tr><tr><td>Madam Vinous 27 wai</td><td>[<xref ref-type="bibr" rid="CR19">19</xref>]</td><td>Madam Vinous</td><td>6 (3C + 3I)</td><td>Leaf</td><td>27 week</td><td>Not available</td><td>Limma R Package</td><td>Unadjusted <italic>P</italic> < 0.05 FC > 2</td><td>1955</td><td>1523</td></tr><tr><td rowspan="4">Resistant (R) leaf dataset</td><td>RL<sup>b</sup> 5 wai</td><td>[<xref ref-type="bibr" rid="CR19">19</xref>]</td><td>Rough Lemon</td><td>6 (3C + 3I)</td><td>Leaf</td><td>5 week</td><td>Not available</td><td>Limma R Package</td><td>Unadjusted <italic>P</italic> < 0.05 FC > 2</td><td>239</td><td>447</td></tr><tr><td>RL 17 wai</td><td>[<xref ref-type="bibr" rid="CR19">19</xref>]</td><td>Rough Lemon</td><td>6 (3C + 3I)</td><td>Leaf</td><td>17 week</td><td>Not available</td><td>Limma R Package</td><td>Unadjusted <italic>P</italic> < 0.05 FC > 2</td><td>260</td><td>154</td></tr><tr><td>RL 27 wai</td><td>[<xref ref-type="bibr" rid="CR19">19</xref>]</td><td>Rough Lemon</td><td>6 (3C + 3I)</td><td>Leaf</td><td>27 week</td><td>Not available</td><td>Limma R Package</td><td>Unadjusted <italic>P</italic> < 0.05 FC > 2</td><td>867</td><td>190</td></tr><tr><td>US-897</td><td>[<xref ref-type="bibr" rid="CR2">2</xref>]</td><td>US-897</td><td>6 (3C + 3I)</td><td>Leaf</td><td>32 week</td><td>GSE30502</td><td>Partek suit two way ANOVA</td><td>Adjusted <italic>P</italic> (FDR < 0.05) FC > 4</td><td>17</td><td>none</td></tr><tr><td rowspan="8">Fruit dataset</td><td>Immature fruit</td><td>[<xref ref-type="bibr" rid="CR40">40</xref>]</td><td>Valencia</td><td>10 (5C + 5I)</td><td>Fruit</td><td>_</td><td>SRP022979</td><td>DESeq package</td><td>_</td><td>45</td><td>24</td></tr><tr><td>Mature fruit</td><td>[<xref ref-type="bibr" rid="CR40">40</xref>]</td><td>Valencia</td><td>10 (5C + 5I)</td><td>Fruit</td><td>_</td><td>SRP022979</td><td>DESeq package</td><td>_</td><td>234</td><td>73</td></tr><tr><td>Hamlin_JV<sup>c</sup></td><td>[<xref ref-type="bibr" rid="CR34">34</xref>]</td><td>Hamlin</td><td>8 (4C + 4I)</td><td>Fruit</td><td>_</td><td>GSE33373</td><td>Limma R Package</td><td>Adjusted <italic>P</italic> (FDR < 0.01) FC > 2</td><td>174</td><td>116</td></tr><tr><td>Hamlin_VT<sup>d</sup></td><td>[<xref ref-type="bibr" rid="CR34">34</xref>]</td><td>Hamlin</td><td>8 (4C + 4I)</td><td>Fruit</td><td>_</td><td>GSE33373</td><td>Limma R Package</td><td>Adjusted <italic>P</italic> (FDR < 0.01) FC > 2</td><td>398</td><td>471</td></tr><tr><td>Hamlin_FF<sup>e</sup></td><td>[<xref ref-type="bibr" rid="CR34">34</xref>]</td><td>Hamlin</td><td>8 (4C + 4I)</td><td>Fruit</td><td>_</td><td>GSE33373</td><td>Limma R Package</td><td>Adjusted <italic>P</italic> (FDR < 0.01) FC > 2</td><td>696</td><td>706</td></tr><tr><td>Valencia_JV</td><td>[<xref ref-type="bibr" rid="CR34">34</xref>]</td><td>Valencia</td><td>8 (4C + 4I)</td><td>Fruit</td><td>_</td><td>GSE33373</td><td>Limma R Package</td><td>Adjusted <italic>P</italic> (FDR < 0.01) FC > 2</td><td>291</td><td>110</td></tr><tr><td>Valencia_VT</td><td>[<xref ref-type="bibr" rid="CR34">34</xref>]</td><td>Valencia</td><td>8 (4C + 4I)</td><td>Fruit</td><td>_</td><td>GSE33373</td><td>Limma R Package</td><td>Adjusted <italic>P</italic> (FDR < 0.01) FC > 2</td><td>288</td><td>257</td></tr><tr><td>Valencia_FF</td><td>[<xref ref-type="bibr" rid="CR34">34</xref>]</td><td>Valencia</td><td>8 (4C + 4I)</td><td>Fruit</td><td>_</td><td>GSE33373</td><td>Limma R Package</td><td>Adjusted <italic>P</italic> (FDR < 0.01) FC > 2</td><td>210</td><td>324</td></tr></tbody></table><table-wrap-foot><p><sup>a</sup>weeks after inoculation, <sup>b</sup>Rough Lemon, <sup>c</sup>Juice Vesicle, <sup>d</sup>Vascular Tissue, <sup>e</sup>Flavedo, <sup>f</sup>Control and Infected</p></table-wrap-foot></table-wrap><fig id="Fig1"><label>Fig. 1</label><caption><p>Venn diagram showing distribution of HLB-responsive differentially expressed probe sets in citrus. Figure is displaying comparison of differentially expressed probe sets of citrus-HLB interactions identified from previous published reports (blue circle) and, from the present meta-analysis by ‘Limma’ (yellow circle) and ‘RankProd’ (green circle) methods</p></caption><graphic xlink:href="12870_2015_568_Fig1_HTML" id="MO1"></graphic></fig></p></sec><sec id="Sec4"><title>Gene enrichment and pathway analysis of differentially expressed probe sets</title><p>Gene enrichment and pathway analysis of differentially expressed probe sets were conducted using MapMan and AgriGO (Additional file <xref rid="MOESM2" ref-type="media">2</xref>). The MapMan ontology was designed specifically for plants [<xref ref-type="bibr" rid="CR59">59</xref>], where probe sets were assigned into largely non-redundant and hierarchically organized BINs based on their annotation. Of the 7,412 probe sets identified with Limma, 5,003 were mapped to 34 different BINs by MapMan, 2,083 fell under a ‘not assigned’ category, and the remaining 2,920 probe sets were mapped to different metabolic pathways. Signalling (<italic>P </italic>< 7.2e-7, cell wall (<italic>P</italic> < 4.9e<sup>−4</sup>), and carbohydrate metabolism (<italic>P</italic> < 0.002) were among the top significant MapMan categories. Of the 4,221 probe sets identified with RankProd, 1,703 were ‘not assigned’, and the rest 2,518 were mapped to different pathways. The RankProd analysis suggested that carbohydrate metabolism was the most significant class modulated by HLB disease in citrus.</p><p>AgriGO is a web-based tool designed specifically for GO analysis of agricultural species [<xref ref-type="bibr" rid="CR14">14</xref>]. GO terms for differentially expressed probe sets from Limma and RankProd were assigned to biological processes, molecular functions, or cellular components. Among biological processes, ‘response to stimulus’ (Limma, FDR = 9.4e<sup>−05</sup>; RankProd, FDR = 0.0013) was the most significant class. Among molecular functions and cellular component terms, ‘catalytic activity’ (Limma, FDR = 0.00023; RankProd, FDR = 1.2e<sup>−05</sup>) and ‘extracellular region’ (Limma, FDR = 8.7e<sup>−08</sup>; RankProd, FDR = 9.8e<sup>−08</sup>) were the most significant classes. The GO analysis of differentially expressed probe sets, using hypergeometric statistics with the Yekutieli FDR multi-test adjustment (FDR < 0.01), indicated that carbohydrate metabolism, secondary metabolism, response to stimulus, cell wall metabolism, nutrient transport, and response to stress were among the most significant metabolic pathways modulated in citrus<bold>–</bold>HLB interactions.</p><p>To study the distribution of these six metabolic pathways (carbohydrate metabolism, cell wall metabolism, hormone metabolism, signalling, secondary metabolism, and transport) among the 22 datasets, we analysed them using the Wilcoxon test in PageMan [<xref ref-type="bibr" rid="CR58">58</xref>]. The results were divided into an S dataset (one each from roots and stems and eight from leaves), an R dataset (resistant leaves), and a fruit dataset (Additional file <xref rid="MOESM3" ref-type="media">3</xref>: Figure S1, Additional file <xref rid="MOESM4" ref-type="media">4</xref>: Figure S2, Additional file <xref rid="MOESM5" ref-type="media">5</xref>: Figure S3, Additional file <xref rid="MOESM6" ref-type="media">6</xref>: Figure S4, Additional file <xref rid="MOESM7" ref-type="media">7</xref>: Figure S5, Additional file <xref rid="MOESM8" ref-type="media">8</xref>: Figure S6). Carbohydrate metabolism generally showed down-regulation in the fruit dataset and during early stages of <italic>Ca</italic>Las infection in both the R and S datasets. Probe sets related to starch degradation were up-regulated in the fruit dataset but down-regulated in the R dataset. However, starch synthesis-related probe sets were up-regulated in the resistant citrus variety US-897 and S datasets but down-regulated in the fruit dataset (Additional file <xref rid="MOESM3" ref-type="media">3</xref>: Figure S1). In total, 125 probe sets related to the cell wall were induced. Most cell wall-related probe sets were induced at 17 weeks after infection (wai) with <italic>Ca</italic>Las and were up-regulated in the R dataset and down-regulated in the S dataset. Probe sets related to cell wall synthesis, such as α-expansin 3 and pectinesterase were induced in tolerant citrus variety Rough Lemon (RL), but cell wall-modification-related probe sets were up-regulated in resistant US-897 (Additional file <xref rid="MOESM4" ref-type="media">4</xref>: Figure S2). In hormone metabolism, brassinosteroid and GA pathway-related probe sets were up-regulated in the R dataset. Few salicylic acid (SA)-responsive probe sets were up-regulated with ethylene-signalling probes in the S dataset (Additional file <xref rid="MOESM5" ref-type="media">5</xref>: Figure S3). In the signalling, all of the probe sets were up-regulated in the R dataset except RL 27 wai and coded for sugar and nutrient signalling, leucine-rich receptor, or wall-associated kinase receptor. Calcium signalling related probe sets were down-regulated in the S but not in the R dataset (Additional file <xref rid="MOESM6" ref-type="media">6</xref>: Figure S4). We found 217 probe sets that were related to transport. Among these, six probe sets encoding aquaporin were up-regulated only in RL 17 wai and were down-regulated in the S dataset. ABC transporters coding probe sets were up-regulated only in the S dataset (Additional file <xref rid="MOESM7" ref-type="media">7</xref>: Figure S5). In the defence and stress pathway, 19 probe sets coding for PR proteins or chitinase were induced in 12 datasets and showed up-regulation in R datasets (US-897 and RL 17 wai). Seventeen probe sets encoding miraculin, a member of the Kunitz serine and trypsin protease inhibitor family proteins, was also induced. Five miraculin-encoding probe sets were up-regulated in RL 17 wai, and six were down-regulated in RL 27 wai (Additional file <xref rid="MOESM8" ref-type="media">8</xref>: Figure S6).</p></sec><sec id="Sec5"><title>Tissue-specific gene expression in citrus–HLB interactions</title><p>An in-house Teradata SQL script was used to extract fruit data corresponding to the identified 7,412 probe sets. Differentially expressed probe sets were further classified according to tissue specificity (root, stem, leaf, or fruit datasets) (Fig. <xref rid="Fig2" ref-type="fig">2</xref>, Additional file <xref rid="MOESM9" ref-type="media">9</xref>). The root dataset included 92 probe sets (Fold change >= 2 or <= -2), of which 45 were common to other datasets and the remaining 47 probe sets were specific to root data. Only 17 of these 47 probe sets were annotated in citrus or <italic>Arabidopsis</italic> genome databases, and they coded for Cu/Zn superoxidase dismutase, ent-kaurenoic oxidase, MYB transcription factor (two probe sets), and polygalacturonase. Stem dataset showed significant expression changes for 672 probe sets, of which 542 probe sets were common to other datasets and 130 probe sets were specific to the stem dataset. These probe sets mainly coded for transporter (<italic>n</italic> = 9), protein signalling (<italic>n</italic> = 8), and the cell wall (<italic>n</italic> = 6). All transport-related probe sets that coded for protein transporters (<italic>n</italic> = 6), peptide and oligopeptide transporters, and ABC transporters were up-regulated in HLB-diseased citrus. Cell wall-related probe sets coding for xyloglucan endotransglycosylase and expansin were down-regulated in the stem dataset. Among signalling-related probe sets, those involved in calcium signalling were down-regulated, but those related to sugar and nutrient physiology were up-regulated in stems. Most of the transcriptome studies for citrus–HLB interactions were performed using leaf tissues, and we retrieved 12 leaf datasets. A total of 7,182 probe sets showed differential expression in at least one leaf dataset, and 5,327 probe sets were unique to leaf datasets, showing no expression changes in other tissues. These 5,327 probe sets were further grouped into four R and eight S datasets; 318 probe sets were specific to the R datasets, 4,003 were specific to the S datasets, and 1,006 were common to both datasets. There were 288 probe sets present in at least one of the four R datasets, and 30 were present in more than two R datasets. From 4,003 probe sets in S datasets, 3,923 were from at least one dataset, and 80 were present in more than four datasets. For the eight fruit datasets, 1,345 probe sets showed significant fold changes in at least one dataset. Because we used fruit datasets only for comparison with our meta-analysis, no probe sets were specific to these datasets.<fig id="Fig2"><label>Fig. 2</label><caption><p>Classification of tissue-specific probe sets from 22 sample datasets. Figure represent the number of differentially expressed probe sets that are unique and common among different citrus tissues datasets (root, stem, leaf and fruit). R and S represent resistant and susceptible datasets of leaf, respectively</p></caption><graphic xlink:href="12870_2015_568_Fig2_HTML" id="MO2"></graphic></fig></p></sec><sec id="Sec6"><title>Conserved expression of top common probe sets from both meta-analysis methods</title><p>Using the TopGene function in RankProd, we identified the top 100 common probe sets (50 up-regulated and 50 down-regulated). We searched the status of these 100 probe sets with Limma results and taken into account only those showed differential expression at least in four datasets. Thus, the 65 most common probe sets from both methods were identified and analysed. Because these 65 probe sets are expressed in most datasets, whether or not in a tissue-specific or response-specific manner, these probe sets should signify core pathways in citrus–HLB interactions. We performed hierarchal clustering of these common probe sets with the 22 datasets using Karl Pearson correlations with average linkage [<xref ref-type="bibr" rid="CR16">16</xref>]. The clusters revealed relatedness among the datasets for expression of the common probe sets. Using a distance threshold of 0.9, we calculated groups/nodes among differentially expressed probe sets and datasets and found five vertical clusters for probe sets and four horizontal clusters for datasets from the main hierarchical clusters (Fig. <xref rid="Fig3" ref-type="fig">3</xref>, Additional file <xref rid="MOESM10" ref-type="media">10</xref>). Cluster I comprised 21 probe sets; of these, 11 were of unknown functions, and others mainly coded for orcinol-<italic>O</italic>-methyltransferase (Cit.12172.1.S1_at and Cit.10244.1.S1_s_at), phloem protein 2(PP2)-B15 (Cit.35955.1.S1_at), NAC domain-containing protein (Cit.12214.1.S1_s_at), and glucose-6-phosphate/phosphate translocator 2 (GPT2) (Cit.9625.1.S1_s_at and Cit.22602.1.S1_at). Out of 21 probe sets, only four were induced in fruit datasets, and all four showed down-regulation. All probe sets were up-regulated in the R and S leaf datasets, except orcinol-<italic>O</italic>-methyltransferase and PP2-B15, which were not induced in any of the R datasets. Cluster II consisted of 26 probe sets that coded for Zn transporter 5 (Cit.11459.1.S1_s_at, Cit.11460.1.S1_at, and Cit.28305.1.S1_at), GASA1 (Cit.10032.1.S1_x_at), ADP-glucose pyrophosphorylase large subunit (APL3) (Cit.13437.1.S1_s_at), starch synthase (Cit.9504.1.S1_s_at), AtWRKY40 (Cit.10816.1.S1_at), CDR (Cit.28117.1.S1_s_at), Cu/Zn superoxide dismutase (Cit.28102.1.S1_s_at), and 2OG-Fe(II) oxygenase family protein (Cit.15355.1.S1_at). Cluster III consisted of six probe sets: β-amylase (Cit.39675.1.S1_at and Cit.36677.1.S1_at), nodulin (Cit.4972.1.S1_s_at), ripening-related proteins (Cit.18669.1.S1_at), and membrane protein-encoding probe sets (Cit.16855.1.S1_at). Cluster IV was made up of two probe sets that coded for chalcone synthase (Cit.13366.1.S1_at) and amino acid transport (Cit.18023.1.S1_at). Cluster V consisted of 10 probe sets mainly coding for ABC transporter (Cit.6534.1.S1_at), peroxidase (Cit.6534.1.S1_at), and 9-cis-epoxycarotenoid dioxygenase 2 (Cit.17235.1.S1_s_at). In cluster V, all of the probe sets were down-regulated in S datasets, but four were up-regulated in one R dataset (RL 17 wai).<fig id="Fig3"><label>Fig. 3</label><caption><p>Hierarchical clustering of 65 common probe sets identified by both meta-analysis methods. The expression datasets were gene-wise normalized and, the clusters for 65 probe sets and 22 datasets were made using Pearson correlation coefficients. Heatmap showed four vertical clusters and four single nodes for 22 datasets. Five horizontal clusters were made for 65 probe sets</p></caption><graphic xlink:href="12870_2015_568_Fig3_HTML" id="MO3"></graphic></fig></p><p>The 22 datasets (Table <xref rid="Tab1" ref-type="table">1</xref>) were grouped into four clustered and four single nodes. Madam Vinous 5 wai, Madam Vinous 27 wai, Valencia root, and RL 27 wai were not part of any cluster and remained as single nodes; US-897 and Valencia 13–17 wai were clustered together. All of the fruit data clustered together except for mature and immature fruit datasets from Martinelli et al. [<xref ref-type="bibr" rid="CR40">40</xref>], which clustered with RL 5 wai, RL 17 wai, and the stem dataset. The remaining S datasets were clustered together. RL 5 wai and RL 17 wai showed gene responses similar to those of the stem dataset, but the US-897 response for common probe sets was similar to that for Valencia 13–17 wai.</p></sec><sec id="Sec7"><title>Differentially expressed probe sets specific to the R datasets</title><p>Thirty differentially expressed probe sets were identified in at least two of the four R datasets (Fig. <xref rid="Fig2" ref-type="fig">2</xref>). The hierarchical cluster analysis using Pearson correlation at a distance threshold of 0.9 showed three vertical clusters for 30 probe sets (7, 10, and 13 probe sets) and one horizontal cluster for two R datasets (US-897 and RL 17 wai). The other two R datasets (RL 5 wai and RL 27 wai) were present as single nodes (Fig. <xref rid="Fig4" ref-type="fig">4a</xref>, Additional file <xref rid="MOESM11" ref-type="media">11</xref>). The seven probe sets comprising cluster I encoded GA-responsive protein (Cit.18244.1.S1_at), leucine-rich repeat (LRR) VIII-2 (Cit.5979.1.S1_at), terpene synthase (Cit.38319.1.S1_s_at), aquaporin (Cit.8763.1.S1_s_at), lectin-related protein precursor (Cit.35004.1.S1_s_at and Cit.8609.1.S1_x_at), and SNARE 11 (Cit.37265.1.S1_at). These seven probe sets were up-regulated in RL 17 wai and down-regulated in RL 27 wai. However, both of the lectin-related protein-coding probe sets were up-regulated in US-897 and RL 17 wai. In cluster II, eight of 10 probe sets were down-regulated in RL 5 wai and RL 17 wai encoding cinnamoyl-CoA reductase (Cit.1800.1.S1_at and Cit.28372.1.S1_at), C<sub>2</sub>H<sub>2</sub> Zn family protein (Cit.6562.1.S1_at), DNAJ heat shock family protein (Cit.2424.1.S1_s_at), putative ABC transporter protein (Cit.23189.1.S1_at), and unknown function proteins (Cit.14457.1.S1_at, Cit.17300.1.S1_s_at, and Cit.28939.1.S1_at). One probe set from Cluster II that encoded an oxidoreductase 2OG-Fe(II) oxygenase family protein (Cit.940.1.S1_s_at) was up-regulated in US-897 and RL 17 wai. Four probe sets from Cluster III that encoded 5-AMP-activated protein kinase beta-1 subunit-related (Cit.24484.1.S1_at and Cit.14610.1.S1_s_at), abscisic stress ripening-like protein (Cit.8661.1.S1_x_at), and unknown protein (Cit.21210.1.S1_s_at) were up-regulated at RL 5 wai and RL 27 wai. Four probe sets that coded for B-box Zn family protein (Cit.18262.1.S1_at), jumonji transcription factor (Cit.7734.1.S1_at), mutT domain protein-like (Cit.23890.1.S1_at), and unknown protein (Cit.7734.1.S1_at) were up-regulated at RL 5 wai and RL 17 wai. From the remaining probe sets of Cluster III, two that encoded an unknown protein (Cit.23443.1.S1_at) were up-regulated in rough lemon at all three time points, and one that encoded NADH dehydrogenase subunit 4 (Cit.29311.1.S1_at) was down-regulated in US-897 and RL 17 wai. Cluster of US-897 and RL 17 wai suggested that both had similar expression changes against <italic>Ca</italic>Las infection for these 30 R-specific probe sets. However, RL 5 wai and RL 27 wai showed distinct regulation patterns of R-specific probe sets.<fig id="Fig4"><label>Fig. 4</label><caption><p><bold>a</bold> Hierarchical clustering of 30 R-specific probe sets. The expression datasets were gene-wise normalized and, the clusters for 30 probe sets and four datasets were made using Pearson correlation coefficients. Heatmap showed one vertical cluster and two single nodes for four datasets. Three horizontal clusters were made for 30 probe sets. <bold>b</bold> Relevance network of 30 R-specific probe sets. Six relevance sub-networks were generated using 30 R-specific probe sets. Sub-networks showed that the expression of one probe set affected the expression of the other probe set either positively (shown in red line) or negatively (shown in blue line). The expression profiles of these probesets among resistant datasets (US-897, RL 5, 17 and 27 wai) are shown in rectangle boxes</p></caption><graphic xlink:href="12870_2015_568_Fig4_HTML" id="MO4"></graphic></fig></p><p>We further analysed the relevance network of these 30 probe sets (Fig. <xref rid="Fig4" ref-type="fig">4b</xref>, Additional file <xref rid="MOESM11" ref-type="media">11</xref>). A relevance network is a group of probe sets whose expression profiles are highly predictive of one another. Each pair of probe sets is related by a correlation coefficient greater than a minimum threshold (0.97) and less than a maximum threshold (1.0) [<xref ref-type="bibr" rid="CR7">7</xref>]. The relevance network analysis resulted in six tightly correlated sub-networks that showed either positive or negative effects of one probe set on the expression of the other member of the pair (red and blue lines, respectively, Fig. <xref rid="Fig4" ref-type="fig">4b</xref>). Sub-network 1 consisted of seven probe sets that regulated one another’s gene expression in coding for unknown proteins (three probe sets), cinnamoyl CoA reductase (two probe sets), heat shock protein, and C<sub>2</sub>H<sub>2</sub> Zn finger protein. Except for unknown protein Cit.14457.1.S1_at, which negatively regulated the expression of the other six probe sets, these probe sets positively regulated one another’s expression. Sub-network 2 consisted of five probe sets (<italic>LRR-VIII-2</italic>, <italic>SNARE11</italic>, terpene synthase, aquaporin, and abscisic stress ripening coding protein). <italic>LRR-VIII-2</italic> positively regulated the cell vessel transport <italic>SNARE 11</italic>. However, both of these probe sets were negatively regulated by abscisic stress ripening coding protein. Sub-network 3 consisted of four probe sets, two of which encoded unknown proteins, and one each for mutT domain protein-like and myo-inositol-1-phosphate synthase. Sub-network 4 included three probe sets that encoded 5-AMP-activated protein kinase beta-1 subunit-related (two probe sets) and an unknown protein. Sub-network 5 contained two probe sets, one each for diacylglycerol kinase and an unknown protein, and sub-network 6 included two probe sets for lectin-related proteins.</p></sec><sec id="Sec8"><title>Gene co-expression network analysis</title><p>Using WGCNA package [<xref ref-type="bibr" rid="CR33">33</xref>], we constructed a weighted co-expression network of the 3,499 probe sets identified by both Limma and RankProd. This network has helped to elucidate tightly co-expressed modules and to identify hub probe sets in the respective modules [<xref ref-type="bibr" rid="CR33">33</xref>, <xref ref-type="bibr" rid="CR46">46</xref>]. Our network construction yielded 2,071 nodes and 134,217 edges with a Pearson correlation coefficient of 0.85 and scale-free topological matrix. The global network was further clustered into 27 co-expressed modules using topological overlap-based average hierarchical clustering and a dissection threshold of 0.2 (Fig. <xref rid="Fig5" ref-type="fig">5</xref>). We computed the module eigengene (ME), which is the first principal component of each module, and we merged highly correlated MEs to yield 21 modules. A dendrogram showing gene modules before and after merging is illustrated in Additional file <xref rid="MOESM12" ref-type="media">12</xref>: Figure S7. A heat map was generated to visualise topological overlap values between pairs of co-expressed probe sets in the modules (Additional file <xref rid="MOESM13" ref-type="media">13</xref>: Figure S8). Functional analysis showed that six of the 21 modules had probe sets with overrepresented biological functions from pathways such as carbohydrate metabolism, stress, cell wall, signalling, and secondary metabolism. These six modules together contained 2,043 probe sets (Additional file <xref rid="MOESM13" ref-type="media">13</xref>: Figure S8: blue = 715, red = 439, black = 315, yellow = 245, magenta = 165, and purple = 164), which represented 58 % of all probe sets used for co-expressed network construction.<fig id="Fig5"><label>Fig. 5</label><caption><p>Co-expression network of 3,499 common probe sets identified in the meta-analysis. The edges with r ≥ 0.95 are shown in red and r ≤ −0.9 are show in black color. Different modules are color-coded. Circle represent a node (probe set) and edges represent connections. The network was drawn and analyzed using Cytoscape 3.1.1</p></caption><graphic xlink:href="12870_2015_568_Fig5_HTML" id="MO5"></graphic></fig></p></sec><sec id="Sec9"><title>Functional enrichment analysis of HLB-responsive modules and identification of hub probe sets</title><p>To understand the biological relevance of the co-expressed modules, we conducted GO and pathway analyses using AgriGO tool and MapMan. Hub probe sets having maximum connections were identified from each module (Additional file <xref rid="MOESM14" ref-type="media">14</xref>). Here, we describe the ‘blue’ module in detail because it contained the largest number of nodes (probe sets) and connections. Functional annotation by AgriGO revealed that carbohydrate metabolic processes (GO: 0005975, <italic>P</italic> = 1.7e<sup>−10</sup>), response to endogenous stimulus (GO: 0009719, <italic>P</italic> = 0.002), hydrolase activity (GO: 0016787, <italic>P</italic> = 6.9e<sup>−06</sup>), and cell wall (GO: 0005618, <italic>P</italic> = 9.3e<sup>−08</sup>) were among the most significant GO terms in the ‘blue’ module. Using MapMan analysis, nodes were classified into different metabolic pathways, including cell wall (51 nodes), RNA and transcription factors (46 nodes), signalling (45 nodes), secondary metabolism (42 nodes), hormone metabolism (33 nodes), and transport (19 nodes) (Fig. <xref rid="Fig6" ref-type="fig">6</xref>). Next, the top hub probe sets were identified in each category based on K<sub>within</sub>, and K<sub>extracted</sub> >50 and node degree >500 from the 21 modules. Hub probe sets in the secondary metabolism pathway coded for alkaloids, transferase, and <italic>O</italic>-methyltransferase. Top probe sets in the signalling pathway encoded for LRR-XI. In relation to hormone signalling, hub probe sets encoded GA-induced proteins, and bHLH and pepsin A protease were the hub probe sets for transcription factors group. Pathway details, K<sub>within</sub>, and K<sub>extracted</sub> for other modules are provided in Additional file <xref rid="MOESM14" ref-type="media">14</xref>.<fig id="Fig6"><label>Fig. 6</label><caption><p>Detailed analysis of the ‘blue’ module. Figure displays secondary metabolism, signaling, RNA and transcription factor, and, hormone metabolism related sub-network with their respective hub probe sets. The edges with r ≥ 0.95 are shown in red and r ≤ −0.9 are show in black color</p></caption><graphic xlink:href="12870_2015_568_Fig6_HTML" id="MO6"></graphic></fig></p><p>We also analysed the co-expression network to identify the hub probe sets among the top 65 common and 30 R-specific probe sets. All of the top 65 common probe sets were connected among different modules in the network. For example, probe sets for lipoxygenase and <italic>O</italic>-methyltransferase were co-expressed in the ‘black’ module; probe sets coding for glucose transporter (GPT2) was co-expressed with <italic>WRKY40</italic> in the ‘grey’ module; probe sets for phloem protein (PP2-B8) and three probe sets encoding Zn transporters (ZIP5) were co-expressed in the ‘paleturquoise’ module. Probe sets for Cu/Zn superoxide dismutase, WBC11, miraculin protein, and CDR were co-expressed in the ‘midnightblue’ module (Additional file <xref rid="MOESM14" ref-type="media">14</xref>). Interestingly, the ‘midnightblue’ module also contained six probe sets that encoded miraculin and others that encoded disease resistance proteins such as NBS-LRR, chitinases, and LRR transmembrane protein kinase. In contrast to the top 65 common probe sets, only 11 of the 30 R-specific probe sets were co-expressed among modules in the network. Lectin protein kinase (two probe sets) and terpene synthase encoding probe sets were co-expressed in the ‘blue’ module. Three probe sets coding for nucleotide binding, Zn finger protein, and NADH dehydrogenase subunit were co-expressed in the ‘greenyellow’ module. The remaining five R-specific probe sets were distributed among different modules (Additional file <xref rid="MOESM14" ref-type="media">14</xref>).</p></sec><sec id="Sec10"><title>MicroRNA–probe set nested networks</title><p>Recently, Zhao et al. [<xref ref-type="bibr" rid="CR68">68</xref>] identified HLB-responsive miRNAs in susceptible citrus cultivars after 10 and 14 weeks post inoculation (wpi) with <italic>Ca</italic>Las. We analysed these functionally validated miRNAs and found their targets among the 3,499 common probe sets identified here. Interestingly, 10 miRNA classes were identified as targeting 24 probe sets (Additional file <xref rid="MOESM15" ref-type="media">15</xref>). We searched these 24 target probe sets in the co-expressed networks and used them as seed nodes to create a nested network by connecting them with first-degree neighbouring nodes. Thus, we identified nested networks for five miRNA classes (csi-miR396, csi-miR166, csi-miR156, csi-miR167, and csi-miR172) (Fig. <xref rid="Fig7" ref-type="fig">7</xref>).<fig id="Fig7"><label>Fig. 7</label><caption><p>miRNA nested networks. Five miRNA families (csi-miR156, csi-miR166, csi-miR167, csi-miR172, and csi-miR396) forming nested network with their respective target probe sets. The targets are shown in red color node while miRNA classes are represented in yellow color. Connected nested networks with target probe sets are shown in green color. The edges with r ≥ 0.95 are shown in red and r ≤ −0.9 are show in black color</p></caption><graphic xlink:href="12870_2015_568_Fig7_HTML" id="MO7"></graphic></fig></p><p>The nested network of csi-miR396 targeted four probe sets: cytosolic ascorbate peroxidase (Cit.8235.1.S1_at), papain-like cysteine protease (Cit.26437.1.S1_s_at), nuclear matrix constituent protein (Cit.16271.1.S1_at), and cysteine protease (Cit.18551.1.S1_at). Among these, cysteine protease was linked to two large hub probe sets coding for LRR transmembrane protein kinase (Cit13717.1.S1_at) and cellulose synthase (Cit.3390.1.S1_at), which in turn were connected to the 36 stress- and secondary metabolism-related probe sets with a clustering coefficient of 0.6. The expression of csi-miR396 was reduced at both time points (10 and 14 wpi) after <italic>Ca</italic>Las infection [<xref ref-type="bibr" rid="CR68">68</xref>]; therefore, its targeted probe sets should be up-regulated. Except for cytosolic ascorbate peroxidase, the other target probe sets (<italic>n</italic> = 3) were up-regulated in at least one S dataset (Additional file <xref rid="MOESM15" ref-type="media">15</xref>).</p><p>The nested network of csi-miR166 targeted five probe sets coding for unknown protein (Cit.12380.1.S1_at), HD-ZIPIII protein (Cit.7481.1.S1_at), spermine synthase (Cit.14666.1.S1_at), HD-Zip protein (Cit.15346.1.S1_s_at), and homeodomain-leucine zipper protein (Cit.11989.1.S1_at). Of these, the first three targets were connected to stress, secondary metabolism, and hormone-related hubs with 213 nodes at a clustering coefficient of 0.65. The strongest connectivity was shown by four hub probe sets coding for oxidoreductase (Cit.3648.2.S1_at), RHM1/UDP-glucose-4,6-dehydratase (Cit.39620.1.S1_s_at), ferredoxin (Cit.6529.1.S1_at), and GA-regulated protein (Cit.30950.1.S1_at). The expression of csi-miR166 was down-regulated at 14 wpi [<xref ref-type="bibr" rid="CR68">68</xref>]. All five target probe sets were down-regulated in at least two S datasets, while unknown protein (Cit.12380.1.S1_at) was up-regulated in R datasets (RL17 wai).</p><p>The nested network of csi-miR156 targeted two probe sets, an unknown protein (Cit.14355.1.S1_at) and the SQUAMOSA promoter-binding protein (Cit.15466.1.S1_at), which were connected to the developmental pathway by three hub probe sets. These hub probe sets coded for NAC domain-containing protein (Cit.26317.1.S1_at), unknown protein (Cit.6138.1.S1_at), and G protein signalling (Cit.13023.1.S1_at). The expression of csi-miR156 was up-regulated at 14 wpi, and our data suggested that both the connected target probe sets were down-regulated in S datasets.</p><p>The nested network of csi-miR167 targeted two probe sets: an unknown protein (Cit.22035.1.S1_a_at) and potassium transporter (Cit.22029.1.S1_at). The unknown protein was connected to oxidoreductase and 2OG-Fe(II) oxygenase family protein. A potassium transporter encoding probe set was co-expressed with many transport-related probe sets, such as sulphate, phosphate (P), metal, peptide, metabolite, and ABC transporters. The expression of csi-miR167 was down-regulated at 10 wpi and up-regulated at 14 wpi in susceptible citrus [<xref ref-type="bibr" rid="CR68">68</xref>]. Our data suggested that an unknown protein-encoding probe set was up-regulated while the potassium transport probe set was down-regulated in S datasets.</p><p>The nested network of csi-miR172 targeted three probe sets: MYB transcription factor (Cit.13801.1.S1_at), GDSL-motif lipase (Cit.30073.1.S1_at), and an unknown protein (Cit.19822.1.S1_s_at). The MYB transcription factor was connected to three hub probe sets coding for gluconeogenesis (Cit.18718.1.S1_s_at), Rubisco (Cit.20550.1.S1_at), and senescence-associated protein-related (Cit.17917.1.S1_at) at a clustering coefficient of 0.77. GDSL-motif lipase was connected with two hub probe sets coding for serine-type endopeptidase and cell wall pectinesterase family protein. Down-regulation of csi-miR172 was reported at both time points (10 and 14 wpi) in citrus [<xref ref-type="bibr" rid="CR68">68</xref>]. Among the target probe sets, MYB transcription factor and unknown protein showed up-regulation in two S datasets, and GDSL-motif lipase showed down-regulation in one S dataset.</p></sec></sec><sec id="Sec11" sec-type="discussion"><title>Discussion</title><p>Development of HLB-resistant or tolerant citrus cultivars would be the most effective, economic, and ecologically friendly approach for managing HLB. To date, no R genes against HLB disease have been reported, and identification of the best candidate genes for functional dissection of the <italic>Ca</italic>Las response in citrus remains a daunting challenge. Several transcriptome studies have been conducted to identify potential citrus candidates against HLB and that have yielded more than 8,000 differentially expressed genes/probe sets (DEGs). However, information about interactions between the identified DEGs is lacking. We conducted a thorough meta-analysis and network co-expression analysis of publically available citrus transcriptome data and identified 65 common and 30 resistance-specific HLB-responsive candidate probe sets. Our study aimed to identify specific features that define the host response and targets manipulated by <italic>Ca</italic>Las infection and HLB development. Prior to this study, only one meta-analysis of citrus–HLB interactions was conducted [<xref ref-type="bibr" rid="CR69">69</xref>], which included four datasets. The authors evaluated early and late HLB-responsive probe sets and suggested that hormone and defence pathways play vital roles in citrus–HLB interactions. Here, we performed a more comprehensive analysis using 22 datasets, a miRNA network, and a co-expression network to expand our understanding of citrus–HLB interactions.</p><p>We used two statistical approaches (Limma and RankProd) to identify probe sets in this meta-analysis, and we restricted the study to experiments conducted on the Affymetrix GeneChip to reduce potential technical bias in transcriptome measurements. The identified differentially expressed probe sets were subjected to Gene enrichment analysis. Our results suggested that carbohydrate metabolism, phloem plugging, transport, cell wall, hormone, defence, and stress-related pathways were the most important classes of probe sets regulated by <italic>Ca</italic>Las infection in citrus. We also found that very few probe sets were modulated in the citrus root dataset and those were not from any specific pathway. Further, transport-related probe sets were abundant in the stem dataset, while leaf datasets included the probe sets from many metabolic pathways. R datasets for leaves consisted mostly of cell wall-related probe sets, gibberellin signalling, and aquaporin transporters; S datasets for leaves were characterized by down-regulation of Ca signalling and up-regulation of ethylene signalling and ABC transporters. We also found that most of the HLB-responsive probe sets were modulated at intermediate stages (17 wai) rather than at early (5 wai) or late (27 wai) stages of infection. Fruit and leaf datasets were characterised by opposing directions of gene regulation; most probe sets that were up-regulated in the leaf dataset were down-regulated in fruit datasets or vice-versa. These results indicate that different tissues of a citrus plant respond in distinct ways to <italic>Ca</italic>Las infection and further susceptible and resistant responses of citrus add up the complexity of the interactions. To explore this further, we identified the 65 most common probe sets that were present in a maximum number of datasets and that had higher levels of expression change in their respective fold values. The 65 most common probe sets indicated that three core pathways were modulated by <italic>Ca</italic>Las infection. These pathways coded for carbohydrate metabolism, transporters, and hormone- and stress-related probe sets and were abundant in at least 10 of 22 citrus datasets. Hierarchical clustering of the 65 probe sets showed that six fruit datasets and five S datasets clustered separately, suggesting tight regulation of these core pathways among different tissues. Because the stem dataset was clustered with two R datasets, we speculate that stem tissues of a susceptible citrus can withstand better than its leaf tissue against <italic>Ca</italic>Las infection. Further, US-897 clustered with Valencia (13 and 17 wai), suggesting that 17 wai could be a crucial time point for the divergence of gene expression among susceptible and resistant citrus. We unveiled connections among common probe sets through WGCNA co-expression and miRNA nested networks. The integrated connections among these probe sets have not been characterised in earlier reports, and therefore our study suggests novel candidates for citrus–HLB interactions.</p><sec id="Sec12"><title>Carbohydrate metabolism</title><p>Carbohydrate metabolism is a critical pathway in citrus–HLB interactions. HLB-diseased citrus trees are characterised by the accumulation of excessive amounts of starch and sucrose in leaves, which disrupts the starch–sucrose pathway between source (leaves) and sink (fruit) tissues. Several reports have shown that genes coding for ADP-glucose pyrophosphorylase (AGPase) large subunit, the rate-limiting enzyme in starch synthesis, was up-regulated and that β-amylase, the rate-limiting enzyme in starch degradation, was down-regulated in HLB-diseased citrus plants [<xref ref-type="bibr" rid="CR1">1</xref>, <xref ref-type="bibr" rid="CR3">3</xref>, <xref ref-type="bibr" rid="CR31">31</xref>]. Another gene coding for glucose-6-phosphate/phosphate translocator (GPT) mediates the import of glucose-6-phosphate, an essential substrate for starch synthesis in plastids. GPT, like AGPase, is associated with the characteristic starch accumulation in HLB-affected plants [<xref ref-type="bibr" rid="CR31">31</xref>]. In this study, these three genes were among the 65 most common probe sets. Co-expression analysis revealed that probe sets for starch synthase and AGPase were tightly connected with a hub probe set encoding NAC domain protein in the ‘plum 1’ module (Additional file <xref rid="MOESM14" ref-type="media">14</xref>). NAC domain transcription factors (TF) are related to plant development and have been implicated as transcriptional switches that regulate secondary wall synthesis [<xref ref-type="bibr" rid="CR70">70</xref>]. Over-expression of NAC TFs in <italic>Arabidopsis</italic> leaves causes activation of secondary wall biosynthetic genes that leads to extensive deposition of secondary walls in cells [<xref ref-type="bibr" rid="CR70">70</xref>]. Our results suggest that up-regulation of this NAC TF could control carbon deposition in citrus secondary cell walls and subsequent thickening in the leaves due to <italic>Ca</italic>Las infection. This study also revealed that <italic>GPT</italic> was co-expressed with <italic>WRKY40</italic>, <italic>MYB15</italic>, and ethylene signalling genes in the ‘grey60’ module. Previous reports also suggested the role of these genes (<italic>WRKY40</italic> and <italic>MYB</italic>) in HLB-symptomatic citrus leaves [<xref ref-type="bibr" rid="CR1">1</xref>, <xref ref-type="bibr" rid="CR3">3</xref>, <xref ref-type="bibr" rid="CR31">31</xref>]. Interestingly, the same <italic>MYB</italic>-like gene was induced nearly 200-fold in symptomatic leaves of susceptible plants infected with <italic>Ca</italic>Las, but not in the resistant (US-897) genotype [<xref ref-type="bibr" rid="CR1">1</xref>, <xref ref-type="bibr" rid="CR3">3</xref>, <xref ref-type="bibr" rid="CR31">31</xref>]. <italic>MYB</italic> genes have been reported as key regulators of sugar-responsive genes during sugar starvation in rice [<xref ref-type="bibr" rid="CR37">37</xref>]. We suggest that the NAC domain, MYB15, and WRKY40 TFs may play important roles in the transcriptional regulation of carbohydrate metabolism in citrus–HLB interactions.</p></sec><sec id="Sec13"><title>Phloem plugging and nutrient transport</title><p>Excessive starch accumulation in citrus phloem cells induced by <italic>Ca</italic>Las infection affects nutrient transport and other transporter genes due to vascular blockage, caused by callose deposition and phloem plugging. Phloem protein encoding probe set (<italic>PP2-B15</italic>), from the 65 most common probe sets, is involved in phloem blockage [<xref ref-type="bibr" rid="CR69">69</xref>]. Three probe sets encoding <italic>bZIP5</italic> were co-expressed with <italic>PP2-B15</italic> and ankyrin-repeat family encoding protein in the ‘paleturquoise’ module. Our findings were in agreement with Zheng and Zhao [<xref ref-type="bibr" rid="CR69">69</xref>], who suggested that Zn transporters and phloem protein were linked in the co-expression network. <italic>PP2-B15</italic> was induced only in the S datasets; further, HLB symptoms resemble those caused by Zn deficiency. Hence, phloem plugging by <italic>PP2-B15</italic> could disrupt Zn transport in infected citrus leaves and may cause symptoms of Zn deficiency. However, the role of ankyrin-repeat protein in phloem plugging needs to be clarified. Cationic plant nutrients such as Ca, K, Mg, Fe, Cu, Mn, and Zn play important roles in citrus–HLB interactions [<xref ref-type="bibr" rid="CR43">43</xref>]. HLB-responsive miRNA (miR399) is up-regulated during phosphorous (P) starvation in susceptible plants [<xref ref-type="bibr" rid="CR68">68</xref>]. Our integrated miRNA nested network showed that an HLB-responsive miRNA, csi-miR167, is linked to a hub probe set that codes for potassium (K) transport, which in turn is connected to the nodes of sulphate, phosphate, metal, peptide, and metabolic transport-related probe sets. K-deficient plants tend to be more susceptible to infection than those with an adequate supply of K [<xref ref-type="bibr" rid="CR49">49</xref>, <xref ref-type="bibr" rid="CR62">62</xref>]. Furthermore, it has been suggested that foliar fertilization with several mineral elements (Zn, Fe, Ca, K, and Mn) reduced HLB symptoms on infected trees [<xref ref-type="bibr" rid="CR47">47</xref>]. Our study suggests that the expression of csi-miR167 could regulate the hub gene for K-transporters, which in turn may modulate the expression of other transport-regulated genes such as P transporters.</p></sec><sec id="Sec14"><title>Hormone and stress response</title><p>Gene reprogramming in plants against pathogen is initiated by rapid respiratory burst and redox as a part of their basal defences [<xref ref-type="bibr" rid="CR48">48</xref>]. The 65 most common probe sets showed enrichment of oxidoreductase, 2OG-Fe(II) oxygenase, Cu/Zn superoxide dismutase, peroxidases, and catalase. Previous reports suggested that these probe sets could sequester free radicals and maintain cell homeostasis in citrus–HLB interactions [<xref ref-type="bibr" rid="CR18">18</xref>, <xref ref-type="bibr" rid="CR43">43</xref>]. Present study showed that 2OG-Fe(II) oxygenase was co-expressed in the ‘Orangered4’ module connected to a hub probe set that encoded a CCCH-type Zn finger protein. Cu/Zn superoxide dismutase was co-expressed with many PR-encoding probe sets in the ‘Midnightblue’ module, which in turn was connected to a hub probe set that encoded a protease inhibitor. We suspect that up-regulation of protease inhibitor after insect vector feeding in turn may up-regulate these free radical-sequestering systems related genes in order to restrict the HLB disease-related membrane damage after <italic>Ca</italic>Las infection.</p></sec><sec id="Sec15"><title>R-specific probe sets in citrus–HLB interactions</title><p>We identified 30 R-specific probe sets that coded primarily for chitinase, lectins, terpene synthase, miraculins, aquaporin, GA-responsive protein, and Zn finger (C<sub>2</sub>H<sub>2</sub> and C<sub>2</sub>C<sub>2</sub> type) family proteins. <italic>Ca</italic>Las is characterised as being more parasitic than pathogenic [<xref ref-type="bibr" rid="CR15">15</xref>] and as lacking a type III secretion system (T3SS). T3SS is required for bacterial effector proteins to be secreted into the host plant, and these effector proteins interact with host R proteins to induce a defence pathway cascade [<xref ref-type="bibr" rid="CR8">8</xref>]. The most common group of R genes in plants is that of the nucleotide binding site-leucine-rich repeat (NBS-LRR) class. <italic>Ca</italic>Las is also known to suppress plant defences [<xref ref-type="bibr" rid="CR61">61</xref>]. Here, none of the eight NBS-LRR-coding probe sets showed a specific pattern of expression in the R datasets. Similar results were reported by Aritua et al. [<xref ref-type="bibr" rid="CR3">3</xref>], who showed that NBS-LRR and other resistance gene probe sets were down-regulated after <italic>Ca</italic>Las infection. The <italic>Ca</italic>Las genome contains most flagellin genes (<italic>fla</italic>), including <italic>flg22</italic> peptide, which could act as a pathogen-associated molecular pattern (PAMP) [<xref ref-type="bibr" rid="CR71">71</xref>]. Flagellin, a well-characterized PAMP, is recognized by the LRR receptor kinase XII (<italic>FLS2)</italic> in <italic>Arabidopsis</italic> [<xref ref-type="bibr" rid="CR32">32</xref>]. Here, among various LRR kinases modulated by HLB, <italic>LRR-II</italic> and <italic>LRR-XII</italic> were up-regulated in the R datasets (RL17 wai) and down-regulated in at least three S datasets. These LRR-family probe sets were also among the hub probe sets in the ‘blue’ module and were connected to a csi-miR166 target probe set encoding HD-Zip III protein. The role of these LRR probe sets in recognizing <italic>Ca</italic>Las effector molecules needs to be investigated.</p><p>A large number of pathogenesis-related probe sets were more abundant in the R than in the S datasets. Two of these probe sets encode1,3-β-glucanase and chitinase, and they were also implicated in <italic>Bois noir</italic> phytoplasma infection of grapevines [<xref ref-type="bibr" rid="CR28">28</xref>]. Like <italic>Ca</italic>Las, phytoplasma (<italic>Candidatus Phytoplasma asteris</italic>) is a cell wall-less bacteria that parasitises plant phloem sieve cells and lacks bacterial T3SS. Co-expression network analysis showed that five probe sets coding for chitinase were co-expressed in the ‘black’ module and connected to the WRKY23-encoding hub probe set. Other R-specific probe sets coded for protease inhibitors such as miraculin and CDR. A large number of miraculin-encoding probe sets (18) were induced in citrus in response to <italic>Ca</italic>Las. Miraculin-like proteins have been reported in the leaves of rough lemon [<xref ref-type="bibr" rid="CR57">57</xref>]; proteomics studies have shown higher accumulation of miraculin proteins after HLB disease [<xref ref-type="bibr" rid="CR18">18</xref>], and miraculin-like proteins were up-regulated in lime trees by phytoplasma [<xref ref-type="bibr" rid="CR53">53</xref>]. The present study showed that probe set coding for miraculin was co-expressed in the ‘midnightblue’ module along with PR4, MLO-like protein, and chitinase encoding probe sets. The roles of miraculins in HLB resistance require investigation. Another protease gene, <italic>CDR</italic>, was thought to be involved in the US-897 resistance mechanism against <italic>Ca</italic>Las and it encodes aspartic protease, which can release endogenous peptides for defence responses [<xref ref-type="bibr" rid="CR64">64</xref>]. Our results showed that <italic>Ca</italic>Las-induced csi-miR396 targets protease probe sets that are co-expressed with LRR receptor hubs, indicating a role of proteases and LRR receptor kinase in HLB defence responses. It should be noted that expression of csi-miR396 decreased with <italic>Ca</italic>Las infection. Reduced expression of miR396 was also documented in transgenic <italic>Arabidopsis</italic> expressing phytoplasma effector protein SAP11 [<xref ref-type="bibr" rid="CR38">38</xref>]. It has been suggested that miR396 is required for cell proliferation during root and leaf development [<xref ref-type="bibr" rid="CR38">38</xref>]. The role of miR396 in HLB defence mechanisms needs to be clarified in future studies.</p><p>Citrus defence against HLB disease seems not to be characterised by jasmonic acid (JA) and/or salicylic acid (SA) signalling [<xref ref-type="bibr" rid="CR31">31</xref>, <xref ref-type="bibr" rid="CR61">61</xref>]. <italic>Ca</italic>Las encodes salicylate hydroxylase, which converts SA into catechol, a product that does not induce resistance [<xref ref-type="bibr" rid="CR60">60</xref>]. The effects of different plant growth regulators have been studied in phytoplasma-affected plants [<xref ref-type="bibr" rid="CR10">10</xref>]. Phytoplasma are also known to suppresses JA signalling responses in plant against insect vectors and also down-regulate the SA-mediated defence responses against bacterial effector protein SAP11 [<xref ref-type="bibr" rid="CR38">38</xref>]. Phytoplasma manipulates development and defence hormone biosynthesis by suppressing the JA pathway and by modulating phosphate homeostasis through triggering the P starvation response [<xref ref-type="bibr" rid="CR52">52</xref>]. Down-regulation of the JA pathway and P starvation has also been documented in <italic>Ca</italic>Las infection [<xref ref-type="bibr" rid="CR31">31</xref>, <xref ref-type="bibr" rid="CR68">68</xref>]. Zhao et al. [<xref ref-type="bibr" rid="CR68">68</xref>] reported that phosphate deficiency worsened HLB symptoms, while application of exogenous P reduced the symptoms. It has also been suggested that under P starvation, plants accumulate sugars and starch in their leaves [<xref ref-type="bibr" rid="CR26">26</xref>]. In the present study, LRR receptor kinase and PR encoding probe sets were co-expressed with the GA signalling pathway related probe sets. GA signalling-related genes such as <italic>GASA</italic> and <italic>GAST</italic> were up-regulated in the R datasets but down-regulated in the S datasets. Our results were in concordance with Martinelli et al. [<xref ref-type="bibr" rid="CR41">41</xref>] showing an opposite pattern of GA-related genes expression; down regulated in fruits but up regulated in leaves. The authors [<xref ref-type="bibr" rid="CR41">41</xref>] also suspected that change in fruit sugar metabolism might be linked with the down-regulation of GA pathway which in turn regulate energy and carbohydrate metabolism. It has also been documented that <italic>GASA5</italic> was constitutively more expressed in resistant US-897, and expression levels were significantly induced in response to infection [<xref ref-type="bibr" rid="CR2">2</xref>]. GA deficiency has been implicated in phytoplasma diseases [<xref ref-type="bibr" rid="CR12">12</xref>]. Ding et al. [<xref ref-type="bibr" rid="CR12">12</xref>] suggested that in potato purple-top phytoplasma infection in tomato, increased GA signalling is required to activate the defence-related enzymes β-1,3-glucanase and chitinase. We propose that GA signalling could co-ordinate defence responses in citrus–HLB interactions. Effect of different plant growth regulators on phytoplasma affected plants has been studied [<xref ref-type="bibr" rid="CR10">10</xref>].</p><p>We found that lectin precursor probe sets were up-regulated in US-897 and RL 17 wai but not in the S datasets. Lectins are a unique group of plant proteins that can recognise and bind to glycol conjugates present on the external bacterial surface [<xref ref-type="bibr" rid="CR44">44</xref>]. Weintraub [<xref ref-type="bibr" rid="CR63">63</xref>] reported that lectins might have direct effects on the insect vector during phytoplasma infection. Snowdrop lectins bind to the insect midgut and are not fully degraded by midgut proteases. Because <italic>Ca</italic>Las is mediated through the insect vector citrus psyllid, we propose that lectins may have direct roles by inhibiting insect feeding on citrus plants. Engineering plants with these lectin genes may deter citrus psyllid, and subsequently HLB incidence.</p></sec></sec><sec id="Sec16" sec-type="conclusion"><title>Conclusion</title><p>Here, we have assembled currently available transcriptome data on differential gene expression in citrus after <italic>Ca</italic>Las infection and upgraded it from mere a ‘list’ to ‘interconnections’ through co-expression and miRNA network analyses. By integrating transcriptomic information, we identified major hub genes for susceptible and resistant response of citrus–HLB interaction. Our data showed that carbohydrate metabolism is regulated by hub genes coding for MYB, NAC and WRKY transcription factor. Further reduced ion transport in susceptible plants could be restore by regulating csi-miR167. Our analysis identified 30 R-specific probe sets as potential candidates against HLB disease. We hypothesise that basal defences in citrus against HLB could be mediated by degradation of <italic>Ca</italic>Las -derived PAMPs such as flagellin or by other unknown effectors which could be recognised by LRR family receptor proteins and chitinases. R gene-mediated defences against HLB may be mediated by plant proteases such as miraculin and CDR. Analysis also suggest that citrus proteases are regulated by csi-miR396 after <italic>Ca</italic>Las infection and csi-miR396 may have a role in defence pathway. Our study indicate that GA signalling could play a vital role in HLB resistance.</p></sec><sec id="Sec17" sec-type="materials|methods"><title>Methods</title><sec id="Sec18"><title>Microarray data acquisition and pre-processing</title><p>We included 14 studies from six reports that are available from the NCBI Gene Expression Omnibus (GEO) database (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/geo">http://www.ncbi.nlm.nih.gov/geo</ext-link>) [<xref ref-type="bibr" rid="CR1">1</xref>–<xref ref-type="bibr" rid="CR3">3</xref>, <xref ref-type="bibr" rid="CR18">18</xref>, <xref ref-type="bibr" rid="CR19">19</xref>, <xref ref-type="bibr" rid="CR31">31</xref>]. All of these studies were based on the Affymetrix system. The GEO series datasets were available for eight studies from five reports (GenBank: GSE33004, GSE33003, GSE33459, GSE30502, and GSE29633). The remaining six studies were from the data of Fan et al. [<xref ref-type="bibr" rid="CR19">19</xref>]. The 14 studies included one dataset for roots (GenBank: GSE33004), one dataset for stems (GenBank: GSE33004), and 12 datasets for leaf tissues (GenBank: GSE29633, GSE30502, GSE33003, GSE33459, and Fan et al. [<xref ref-type="bibr" rid="CR19">19</xref>]). Our study did not include meta-analysis of fruit tissues, but we included fruit transcriptome data from two published studies [<xref ref-type="bibr" rid="CR34">34</xref>, <xref ref-type="bibr" rid="CR40">40</xref>] and compared it with our dataset using an in-house Teradata database SQL script (<ext-link ext-link-type="uri" xlink:href="http://www.teradata.com/products-and-services/database">http://www.teradata.com/products-and-services/database</ext-link>). From different combinations of fruit studies mentioned in these reports, we included HLB-symptomatic fruit probe sets against their respective healthy controls in order to reduce the complexity of the fruit data. We included six fruit datasets from Liao and Burns [<xref ref-type="bibr" rid="CR34">34</xref>] and two datasets from Martinelli et al. [<xref ref-type="bibr" rid="CR40">40</xref>]. Therefore, this study analyse transcriptional changes in HLB-responsive probe sets from 22 datasets after including eight studies on fruit datasets. These 22 datasets include different time points (early and late time points after <italic>Ca</italic>Las infection), different tissues (roots, stems, leaves, and fruit), and different host responses against HLB disease (susceptibility and tolerance/resistance). The susceptible group consisted of 18 datasets, and the resistant group included four datasets. The accession numbers for the data sets are indicated in Table <xref rid="Tab1" ref-type="table">1</xref>. Background correction and normalisation of the raw data sets were performed using RMA implemented in the affy package in R [<xref ref-type="bibr" rid="CR22">22</xref>].</p></sec><sec id="Sec19"><title>Statistical analysis and identification of differentially expressed probe sets</title><p>To identify differentially expressed probe sets, two meta-analysis approaches were applied to the normalized gene expression datasets. The first approach, the Limma method (linear modelling of microarray data) used the Limma package of Bioconductor [<xref ref-type="bibr" rid="CR51">51</xref>]. Limma computes moderated <italic>t</italic>-statistics and log-odds of differential expression by empirical Bayes shrinkage of standard errors towards a common value. Probe sets with <italic>P</italic> < 0.05 (FDR < 0.1) and fold change were considered. Limma generated separate files for each study, and these files were combined using an in-house SQL script on the Teradata database (<ext-link ext-link-type="uri" xlink:href="http://www.teradata.com/products-and-services/database">http://www.teradata.com/products-and-services/database</ext-link>). A second nonparametric meta-analysis, the RankProd method, combines <italic>P</italic> values for identifying probe sets from individual experiments [<xref ref-type="bibr" rid="CR6">6</xref>]. We used the RPadvance function in the Bioconductor package [<xref ref-type="bibr" rid="CR29">29</xref>], which is specifically designed for meta-analysis. The RankProd output was a single file in the form of Tables <xref rid="Tab1" ref-type="table">1</xref> and 2 (down-regulated and up-regulated HLB-responsive probe sets, respectively), gene expression, FC, <italic>P</italic>-values, and percentage of false predictions (PFP). The number of permutation tests was set to 250, and the TopGene function with a PFP cut-off value of <0.01 was used to identify the top differentially expressed probe sets among the studies.</p></sec><sec id="Sec20"><title>Functional enrichment analysis of differentially expressed probe sets</title><p>We used the singular enrichment analysis method AgriGO (<ext-link ext-link-type="uri" xlink:href="http://bioinfo.cau.edu.cn/agriGO/">http://bioinfo.cau.edu.cn/agriGO/</ext-link>) with default settings for Fisher’s <italic>t</italic>-test (<italic>P</italic> < 0.05) and FDR correction by the Hochberg method for the GO analysis. Probe sets were annotated against the <italic>Arabidopsis</italic> genome, and citrus orthologs were identified using the HarvEST database (<ext-link ext-link-type="uri" xlink:href="http://harvest-web.org">http://harvest-web.org</ext-link>). Functions and metabolic pathways of probe sets were visualized using MapMan with the <italic>Citrus sinensis</italic> mapping file (Csinensis_154.txt) (<ext-link ext-link-type="uri" xlink:href="http://mapman.gabipd.org/">http://mapman.gabipd.org/</ext-link>). The PageMan analysis plugin of MapMan was used to visualise differences among HLB-responsive tissue-specific metabolic pathways using Wilcoxon tests, no correction, and an over-representation analysis (ORA) cutoff value of 1.0. Hierarchical clustering of microarray data, using Pearson correlation and average linkage method with distance threshold of 0.90, was performed using multiarray viewer software from TIAR (<ext-link ext-link-type="uri" xlink:href="http://www.tm4.org/mev.html">http://www.tm4.org/mev.html</ext-link>).</p></sec><sec id="Sec21"><title>Weighted gene co-expression network construction</title><p>Differentially expressed probe sets were taken to infer a weighted gene co-expression network by using the WGCNA package in R version 1.27.1 [<xref ref-type="bibr" rid="CR33">33</xref>]. Pairwise Pearson correlations (<italic>r</italic><sup>2</sup>) were calculated for probe sets across all samples to find correlations and to generate similarity matrices. Similarity (S<sub><italic>ij</italic></sub>) in correlation between the <italic>i</italic><sup>th</sup> and <italic>j</italic><sup>th</sup> probe sets was calculated using the following equation [<xref ref-type="bibr" rid="CR6">6</xref>]: S<sub><italic>ij</italic></sub> = abs|cor(X<sub><italic>i</italic></sub> – X<sub><italic>j</italic></sub>)|, where, X<sub><italic>i</italic></sub> and X<sub><italic>j</italic></sub> = log (base 2) ratio of expression of the <italic>i</italic>th and <italic>j</italic>th genes across the samples, respectively; cor = Pearson correlation coefficient; and abs = absolute value. The absolute value of the Pearson correlation coefficient was used to generate an undirected, weighted network. The similarity matrix was weighted to adjacency matrices by raising it to a power (β). The PickSoftThreshold function of WGCNA was used to choose the appropriate power for the network topology from various soft-thresholding powers. The scale-free network was rendered by raising the soft thresholding power (β) to 14. At this threshold power, the model was fitted with <italic>r</italic><sup>2</sup> < 0.85. This similarity matrix was transformed into the adjacency matrix, which was transformed into a topological overlap matrix (TOM) similarity measure, a robust measure of pairwise interconnectedness [<xref ref-type="bibr" rid="CR65">65</xref>]. The TOM matrix of probe sets was coupled with average linkage hierarchical clustering to cluster the probe sets into distinct modules using the Dynamic Tree Cut algorithm (cutreeDynamic method; deepSplit = 3, cutheight = 0.993; minimal module size = 40) [<xref ref-type="bibr" rid="CR33">33</xref>]. Further, similar modules were merged using parameter cut height = 0.2. The co-expressed network was visualised using Cytoscape version 3.1.1 and analysed using the Network Analyser plugin [<xref ref-type="bibr" rid="CR13">13</xref>]. We considered edges only above a threshold of 0.2 to simplify and concentrate on relevant functions of the co-expressed network. Hub probe sets in the constructed network were identified as suggested by Puniya et al. [<xref ref-type="bibr" rid="CR46">46</xref>].</p></sec><sec id="Sec22"><title>miRNA target prediction and nested network construction</title><p>We obtained the list of miRNAs that targeted HLB-responsive probe sets from three sources: PMRD, the plant miRNA database [<xref ref-type="bibr" rid="CR67">67</xref>], and the <italic>Citrus sinensis</italic> annotation project (<ext-link ext-link-type="uri" xlink:href="http://citrus.hzau.edu.cn/orange/">http://citrus.hzau.edu.cn/orange/</ext-link>), and validated HLB-responsive miRNAs from Zhao et al. [<xref ref-type="bibr" rid="CR68">68</xref>]. Target genes interacting with each HLB responsive miRNA class from Zhao et al. [<xref ref-type="bibr" rid="CR68">68</xref>] were searched in the <italic>Citrus sinensis</italic> annotation project (<ext-link ext-link-type="uri" xlink:href="http://citrus.hzau.edu.cn/orange/">http://citrus.hzau.edu.cn/orange/</ext-link>). We identified 64 target genes for 14 miRNA classes. The corresponding probe sets for these 64 target genes were searched from our 7,412 differentially expressed probe sets based on common <italic>Arabidopsis</italic> annotation and NCBI BLASTX (<ext-link ext-link-type="uri" xlink:href="http://blast.ncbi.nlm.nih.gov/">http://blast.ncbi.nlm.nih.gov/</ext-link>), resulting 24 target probe sets for 10 miRNA classes. The nested networks of 10 miRNAs and their 24 target probe sets were inferred using the nested network plugin in Cytoscape 3.1.1 with default settings.</p></sec></sec></body><back><app-group><app id="App1"><sec id="Sec23"><title>Additional files</title><p><media position="anchor" xlink:href="12870_2015_568_MOESM1_ESM.xls" id="MOESM1"><label>Additional file 1:</label><caption><p><bold>Probe sets identified in the present analysis.</bold> List of the total number of probe sets identified with Limma and RankProd methods of meta-analysis.</p></caption></media><media position="anchor" xlink:href="12870_2015_568_MOESM2_ESM.xls" id="MOESM2"><label>Additional file 2:</label><caption><p><bold>Gene Ontology of identified probe sets.</bold> Gene enrichment analysis based on MapMan and AgriGO.</p></caption></media><media position="anchor" xlink:href="12870_2015_568_MOESM3_ESM.tif" id="MOESM3"><label>Additional file 3: Figure S1.</label><caption><p>PageMan display of coordinated changes of “major CHO metabolism” gene categories in 22 data set.</p></caption></media><media position="anchor" xlink:href="12870_2015_568_MOESM4_ESM.tif" id="MOESM4"><label>Additional file 4: Figure S2.</label><caption><p>PageMan display of coordinated changes of “cell wall” gene categories in 22 data set.</p></caption></media><media position="anchor" xlink:href="12870_2015_568_MOESM5_ESM.jpg" id="MOESM5"><label>Additional file 5: Figure S3.</label><caption><p>PageMan display of coordinated changes of “hormone metabolism” gene categories in 22 data set.</p></caption></media><media position="anchor" xlink:href="12870_2015_568_MOESM6_ESM.jpg" id="MOESM6"><label>Additional file 6: Figure S4.</label><caption><p>PageMan display of coordinated changes of “signaling” gene categories in 22 data set.</p></caption></media><media position="anchor" xlink:href="12870_2015_568_MOESM7_ESM.jpg" id="MOESM7"><label>Additional file 7: Figure S5.</label><caption><p>PageMan display of coordinated changes of “transport” gene categories in 22 data set.</p></caption></media><media position="anchor" xlink:href="12870_2015_568_MOESM8_ESM.jpg" id="MOESM8"><label>Additional file 8: Figure S6.</label><caption><p>PageMan display of coordinated changes of “stress” gene categories in 22 data set.</p></caption></media><media position="anchor" xlink:href="12870_2015_568_MOESM9_ESM.xls" id="MOESM9"><label>Additional file 9:</label><caption><p><bold>List of tissue-specific probe sets.</bold> Description of unique and common tissue-specific probe sets identified in root, stem, leaf and fruit datasets.</p></caption></media><media position="anchor" xlink:href="12870_2015_568_MOESM10_ESM.xls" id="MOESM10"><label>Additional file 10:</label><caption><p><bold>List of 65 top common probe sets and their expression values.</bold> Description of hierarchical clusters identified in the top 65 common probe sets and their fold values among 22 datasets.</p></caption></media><media position="anchor" xlink:href="12870_2015_568_MOESM11_ESM.xls" id="MOESM11"><label>Additional file 11:</label><caption><p><bold>List of 30 R-specific probe sets and their expression values.</bold> Description of hierarchical clusters identified in 30 R-specific probe sets and their fold values among four HLB-resistant datasets (Sheet 1). List of probe sets present in relevance network (Sheet 2).</p></caption></media><media position="anchor" xlink:href="12870_2015_568_MOESM12_ESM.jpg" id="MOESM12"><label>Additional file 12: Figure S7.</label><caption><p>Dendrogram for co-expressed modules. Dendrogram showing co-expressed modules before and after merging using the WGCNA package.</p></caption></media><media position="anchor" xlink:href="12870_2015_568_MOESM13_ESM.jpg" id="MOESM13"><label>Additional file 13: Figure S8.</label><caption><p>Network Heatmap plot for co-expressed modules. Heatmap view of topological overlap values in different modules. Red color shows highly co-expressed probe sets forming a module.</p></caption></media><media position="anchor" xlink:href="12870_2015_568_MOESM14_ESM.xls" id="MOESM14"><label>Additional file 14:</label><caption><p><bold>Description of WGCNA networks and hub genes.</bold> Co-expression module connectivity and description for the identified 21 modules. Description of common and R-specific probe sets among different co-expression modules.</p></caption></media><media position="anchor" xlink:href="12870_2015_568_MOESM15_ESM.xls" id="MOESM15"><label>Additional file 15:</label><caption><p><bold>Detail of miRNA- probe set interactions.</bold> Description of 10 miRNA classes and their 24 targeted probe sets with respective expression profile among 22 datasets.</p></caption></media></p></sec></app></app-group><glossary><title>Abbreviations</title><def-list><def-item><term>ACP</term><def><p>Asian citrus psyllid</p></def></def-item><def-item><term>AGPase</term><def><p>ADP-glucose pyrophosphorylase</p></def></def-item><def-item><term><italic>Ca</italic>Las</term><def><p><italic>Candidatus</italic> Liberibacter asiaticus</p></def></def-item><def-item><term>CDR</term><def><p>Constitutive disease resistance</p></def></def-item><def-item><term>FDR</term><def><p>False discovery rate</p></def></def-item><def-item><term>GA</term><def><p>Gibberellin</p></def></def-item><def-item><term>GEO</term><def><p>Gene expression omnibus</p></def></def-item><def-item><term>GO</term><def><p>Gene ontology</p></def></def-item><def-item><term>GPT</term><def><p>Glucose-6-phosphate/phosphate translocator</p></def></def-item><def-item><term>HLB</term><def><p>Huanglongbing</p></def></def-item><def-item><term>JA</term><def><p>Jasmonic acid</p></def></def-item><def-item><term>ME</term><def><p>Module eigengene</p></def></def-item><def-item><term>NBS-LRR</term><def><p>Nucleotide-binding site leucine-rich repeat</p></def></def-item><def-item><term>P</term><def><p>Phosphorous</p></def></def-item><def-item><term><italic>P</italic></term><def><p>Probability value</p></def></def-item><def-item><term>PP2</term><def><p>Phloem protein</p></def></def-item><def-item><term>PR</term><def><p>Pathogenesis-related</p></def></def-item><def-item><term>R</term><def><p>Resistant</p></def></def-item><def-item><term>RL</term><def><p>Rough lemon</p></def></def-item><def-item><term>RMA</term><def><p>Robust multiChip analysis</p></def></def-item><def-item><term>S</term><def><p>Susceptible</p></def></def-item><def-item><term>SA</term><def><p>Salicylic acid</p></def></def-item><def-item><term>SQL</term><def><p>Structured query language</p></def></def-item><def-item><term>WGCNA</term><def><p>Weighted gene co-expression network analysis</p></def></def-item></def-list></glossary><fn-group><fn><p><bold>Competing interests</bold></p><p>The authors declare that they have no competing interests.</p></fn><fn><p><bold>Authors’ contributions</bold></p><p>NR and DD collected citrus-HLB transcriptome data. NR and SPK executed meta-analysis and co-expression network analysis of the retrieved data. NR conceived, designed and wrote the paper. SPK wrote the custom scripts, helped in the troubleshooting during the analysis. ZD and FG contributed by critically evaluating the study and improving the manuscript. All authors read and approved the final manuscript.</p></fn></fn-group><ack><title>Acknowledgement</title><p>This research was financially supported by grants from the Citrus Research and Development Foundation, Inc. (CDRF) (project number #108766 and #105077). We are thankful for Ravi Pratap Singh (database administrator) for writing custom design scripts and providing technical help with Teradata database.</p></ack><ref-list id="Bib1"><title>References</title><ref id="CR1"><label>1.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Albrecht</surname><given-names>U</given-names></name><name><surname>Bowman</surname><given-names>KD</given-names></name></person-group><article-title>Gene expression in <italic>Citrus sinensis</italic> (L.) Osbeck following infection with the bacterial pathogen ‘<italic>Candidatus Liberibacter asiaticus’</italic> causing Huanglongbing in Florida</article-title><source>Plant Sci</source><year>2008</year><volume>175</volume><issue>3</issue><fpage>291</fpage><lpage>306</lpage><pub-id pub-id-type="doi">10.1016/j.plantsci.2008.05.001</pub-id></element-citation></ref><ref id="CR2"><label>2.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Albrecht</surname><given-names>U</given-names></name><name><surname>Bowman</surname><given-names>KD</given-names></name></person-group><article-title>Transcriptional response of susceptible and tolerant citrus to infection with ‘<italic>Candidatus Liberibacter asiaticus’</italic></article-title><source>Plant Sci</source><year>2012</year><volume>185–186</volume><fpage>118</fpage><lpage>30</lpage><pub-id pub-id-type="doi">10.1016/j.plantsci.2011.09.008</pub-id><pub-id pub-id-type="pmid">22325873</pub-id></element-citation></ref><ref id="CR3"><label>3.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Aritua</surname><given-names>V</given-names></name><name><surname>Achor</surname><given-names>D</given-names></name><name><surname>Gmitter</surname><given-names>FG</given-names></name><name><surname>Albrigo</surname><given-names>G</given-names></name><name><surname>Wang</surname><given-names>N</given-names></name></person-group><article-title>Transcriptional and microscopic analyses of citrus stem and root responses to ‘<italic>Candidatus Liberibacter asiaticus’</italic> infection</article-title><source>PloS One</source><year>2013</year><volume>8</volume><issue>9</issue><pub-id pub-id-type="doi">10.1371/journal.pone.0073742</pub-id><pub-id pub-id-type="pmid">24058486</pub-id></element-citation></ref><ref id="CR4"><label>4.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bhargava</surname><given-names>A</given-names></name><name><surname>Clabaugh</surname><given-names>I</given-names></name><name><surname>To</surname><given-names>JP</given-names></name><name><surname>Maxwell</surname><given-names>BB</given-names></name><name><surname>Chiang</surname><given-names>YH</given-names></name><name><surname>Schaller</surname><given-names>GE</given-names></name><etal></etal></person-group><article-title>Identification of cytokinin-responsive genes using microarray meta-analysis and RNA-seq in <italic>Arabidopsis</italic></article-title><source>Plant Physiol</source><year>2013</year><volume>162</volume><issue>1</issue><fpage>272</fpage><lpage>94</lpage><pub-id pub-id-type="doi">10.1104/pp.113.217026</pub-id><pub-id pub-id-type="pmid">23524861</pub-id></element-citation></ref><ref id="CR5"><label>5.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bové</surname><given-names>JM</given-names></name></person-group><article-title>Huanglongbing: a destructive, newly-emerging, century-old disease of citrus</article-title><source>J Plant Pathol</source><year>2006</year><volume>88</volume><fpage>7</fpage><lpage>37</lpage></element-citation></ref><ref id="CR6"><label>6.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Breitling</surname><given-names>R</given-names></name><name><surname>Amtmann</surname><given-names>A</given-names></name><name><surname>Herzyk</surname><given-names>P</given-names></name></person-group><article-title>Graph-based iterative group analysis enhances microarray interpretation</article-title><source>BMC Bioinformatics</source><year>2004</year><volume>5</volume><fpage>100</fpage><pub-id pub-id-type="doi">10.1186/1471-2105-5-100</pub-id><pub-id pub-id-type="pmid">15272936</pub-id></element-citation></ref><ref id="CR7"><label>7.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Butte</surname><given-names>AJ</given-names></name><name><surname>Tamayo</surname><given-names>P</given-names></name><name><surname>Slonim</surname><given-names>D</given-names></name><name><surname>Golub</surname><given-names>TR</given-names></name><name><surname>Kohane</surname><given-names>IS</given-names></name></person-group><article-title>Discovering functional relationships between RNA expression and chemotherapeutic susceptibility using relevance networks</article-title><source>Proc Natl Acad Sci USA</source><year>2000</year><volume>97</volume><issue>22</issue><fpage>12182</fpage><lpage>6</lpage><pub-id pub-id-type="doi">10.1073/pnas.220392197</pub-id><pub-id pub-id-type="pmid">11027309</pub-id></element-citation></ref><ref id="CR8"><label>8.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Buttner</surname><given-names>D</given-names></name><name><surname>He</surname><given-names>SY</given-names></name></person-group><article-title>Type III protein secretion in plant pathogenic bacteria</article-title><source>Plant Physiol</source><year>2009</year><volume>150</volume><issue>4</issue><fpage>1656</fpage><lpage>64</lpage><pub-id pub-id-type="doi">10.1104/pp.109.139089</pub-id><pub-id pub-id-type="pmid">19458111</pub-id></element-citation></ref><ref id="CR9"><label>9.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Canales</surname><given-names>J</given-names></name><name><surname>Moyano</surname><given-names>TC</given-names></name><name><surname>Villarroel</surname><given-names>E</given-names></name><name><surname>Gutiérrez</surname><given-names>RA</given-names></name></person-group><article-title>Systems analysis of transcriptome data provides new hypotheses about <italic>Arabidopsis</italic> root response to nitrate treatments</article-title><source>Front Plant Sci</source><year>2014</year><volume>5</volume><fpage>22</fpage><pub-id pub-id-type="doi">10.3389/fpls.2014.00022</pub-id><pub-id pub-id-type="pmid">24570678</pub-id></element-citation></ref><ref id="CR10"><label>10.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chang</surname><given-names>C</given-names></name></person-group><article-title>Pathogenicity of aster yellows phytoplasma and <italic>Spiroplasma citri</italic> on periwinkle</article-title><source>Phytopathology</source><year>1998</year><volume>88</volume><issue>12</issue><fpage>1347</fpage><lpage>50</lpage><pub-id pub-id-type="doi">10.1094/PHYTO.1998.88.12.1347</pub-id><pub-id pub-id-type="pmid">18944838</pub-id></element-citation></ref><ref id="CR11"><label>11.</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Da Graca</surname><given-names>J</given-names></name><name><surname>Korsten</surname><given-names>L</given-names></name></person-group><person-group person-group-type="editor"><name><surname>Naqvi</surname><given-names>S</given-names></name></person-group><article-title>Citrus Huanglongbing: Review, present status and future strategies</article-title><source>Diseases of Fruits and Vegetables</source><year>2004</year><publisher-loc>The Netherlands</publisher-loc><publisher-name>Kluwer</publisher-name><fpage>229</fpage><lpage>45</lpage></element-citation></ref><ref id="CR12"><label>12.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ding</surname><given-names>Y</given-names></name><name><surname>Wu</surname><given-names>W</given-names></name><name><surname>Wei</surname><given-names>W</given-names></name><name><surname>Davis</surname><given-names>R</given-names></name><name><surname>Lee</surname><given-names>I</given-names></name><name><surname>Hammond</surname><given-names>R</given-names></name><etal></etal></person-group><article-title>Potato purple top phytoplasma‐induced disruption of gibberellin homeostasis in tomato plants</article-title><source>Ann Appl Biol</source><year>2013</year><volume>162</volume><issue>1</issue><fpage>131</fpage><lpage>9</lpage><pub-id pub-id-type="doi">10.1111/aab.12008</pub-id></element-citation></ref><ref id="CR13"><label>13.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Doncheva</surname><given-names>NT</given-names></name><name><surname>Assenov</surname><given-names>Y</given-names></name><name><surname>Domingues</surname><given-names>FS</given-names></name><name><surname>Albrecht</surname><given-names>M</given-names></name></person-group><article-title>Topological analysis and interactive visualization of biological networks and protein structures</article-title><source>Nat Protoc</source><year>2012</year><volume>7</volume><issue>4</issue><fpage>670</fpage><lpage>85</lpage><pub-id pub-id-type="doi">10.1038/nprot.2012.004</pub-id><pub-id pub-id-type="pmid">22422314</pub-id></element-citation></ref><ref id="CR14"><label>14.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Du</surname><given-names>Z</given-names></name><name><surname>Zhou</surname><given-names>X</given-names></name><name><surname>Ling</surname><given-names>Y</given-names></name><name><surname>Zhang</surname><given-names>Z</given-names></name><name><surname>Su</surname><given-names>Z</given-names></name></person-group><article-title>AgriGO: a GO analysis toolkit for the agricultural community</article-title><source>Nucleic Acids Res</source><year>2010</year><volume>38</volume><issue>Web Server issue</issue><fpage>W64</fpage><lpage>70</lpage><pub-id pub-id-type="doi">10.1093/nar/gkq310</pub-id><pub-id pub-id-type="pmid">20435677</pub-id></element-citation></ref><ref id="CR15"><label>15.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Duan</surname><given-names>Y</given-names></name><name><surname>Zhou</surname><given-names>L</given-names></name><name><surname>Hall</surname><given-names>DG</given-names></name><name><surname>Li</surname><given-names>W</given-names></name><name><surname>Doddapaneni</surname><given-names>H</given-names></name><name><surname>Lin</surname><given-names>H</given-names></name><etal></etal></person-group><article-title>Complete genome sequence of citrus Huanglongbing bacterium, <italic>‘Candidatus Liberibacter asiaticus’</italic> obtained through metagenomics</article-title><source>Mol Plant-Microbe Interact</source><year>2009</year><volume>22</volume><issue>8</issue><fpage>1011</fpage><lpage>20</lpage><pub-id pub-id-type="doi">10.1094/MPMI-22-8-1011</pub-id><pub-id pub-id-type="pmid">19589076</pub-id></element-citation></ref><ref id="CR16"><label>16.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Eisen</surname><given-names>MB</given-names></name><name><surname>Spellman</surname><given-names>PT</given-names></name><name><surname>Brown</surname><given-names>PO</given-names></name><name><surname>Botstein</surname><given-names>D</given-names></name></person-group><article-title>Cluster analysis and display of genome-wide expression patterns</article-title><source>Proc Natl Acad Sci USA</source><year>1998</year><volume>95</volume><issue>25</issue><fpage>14863</fpage><lpage>8</lpage><pub-id pub-id-type="doi">10.1073/pnas.95.25.14863</pub-id><pub-id pub-id-type="pmid">9843981</pub-id></element-citation></ref><ref id="CR17"><label>17.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fagen</surname><given-names>JR</given-names></name><name><surname>Leonard</surname><given-names>MT</given-names></name><name><surname>McCullough</surname><given-names>CM</given-names></name><name><surname>Edirisinghe</surname><given-names>JN</given-names></name><name><surname>Henry</surname><given-names>CS</given-names></name><name><surname>Davis</surname><given-names>MJ</given-names></name><etal></etal></person-group><article-title>Comparative genomics of cultured and uncultured strains suggests genes essential for free-living growth of Liberibacter</article-title><source>PloS One</source><year>2014</year><volume>9</volume><issue>1</issue><pub-id pub-id-type="doi">10.1371/journal.pone.0084469</pub-id><pub-id pub-id-type="pmid">24416233</pub-id></element-citation></ref><ref id="CR18"><label>18.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fan</surname><given-names>J</given-names></name><name><surname>Chen</surname><given-names>C</given-names></name><name><surname>Yu</surname><given-names>Q</given-names></name><name><surname>Brlansky</surname><given-names>RH</given-names></name><name><surname>Li</surname><given-names>Z</given-names></name><name><surname>Gmitter</surname><given-names>FG</given-names></name></person-group><article-title>Comparative iTRAQ proteome and transcriptome analyses of sweet orange infected by <italic>‘Candidatus Liberibacter asiaticus’</italic></article-title><source>Physiol Plantarum</source><year>2011</year><volume>143</volume><issue>3</issue><fpage>235</fpage><lpage>45</lpage><pub-id pub-id-type="doi">10.1111/j.1399-3054.2011.01502.x</pub-id></element-citation></ref><ref id="CR19"><label>19.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fan</surname><given-names>J</given-names></name><name><surname>Chen</surname><given-names>C</given-names></name><name><surname>Yu</surname><given-names>Q</given-names></name><name><surname>Khalaf</surname><given-names>A</given-names></name><name><surname>Achor</surname><given-names>DS</given-names></name><name><surname>Brlansky</surname><given-names>RH</given-names></name><etal></etal></person-group><article-title>Comparative transcriptional and anatomical analyses of tolerant rough lemon and susceptible sweet orange in response to <italic>‘Candidatus Liberibacter asiaticus’</italic> Infection</article-title><source>Mol Plant-Microbe Interact</source><year>2012</year><volume>25</volume><issue>11</issue><fpage>1396</fpage><lpage>407</lpage><pub-id pub-id-type="doi">10.1094/MPMI-06-12-0150-R</pub-id><pub-id pub-id-type="pmid">22809274</pub-id></element-citation></ref><ref id="CR20"><label>20.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Folimonova</surname><given-names>SY</given-names></name><name><surname>Robertson</surname><given-names>CJ</given-names></name><name><surname>Garnsey</surname><given-names>SM</given-names></name><name><surname>Gowda</surname><given-names>S</given-names></name><name><surname>Dawson</surname><given-names>WO</given-names></name></person-group><article-title>Examination of the responses of different genotypes of citrus to Huanglongbing (citrus greening) under different conditions</article-title><source>Phytopathology</source><year>2009</year><volume>99</volume><issue>12</issue><fpage>1346</fpage><lpage>54</lpage><pub-id pub-id-type="doi">10.1094/PHYTO-99-12-1346</pub-id><pub-id pub-id-type="pmid">19900000</pub-id></element-citation></ref><ref id="CR21"><label>21.</label><mixed-citation publication-type="other">Garnier M, Danel N, Bove´ J. The greening organism is a gram negative bacterium. In: Garnsey SM, Timmer LW, Dodds JA, editors. Proceedings of the Ninth Conference of the International Organization of Citrus Virologists. Riverside CA; 1984. p. 115–124.</mixed-citation></ref><ref id="CR22"><label>22.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gautier</surname><given-names>L</given-names></name><name><surname>Cope</surname><given-names>L</given-names></name><name><surname>Bolstad</surname><given-names>BM</given-names></name><name><surname>Irizarry</surname><given-names>RA</given-names></name></person-group><article-title>affy--analysis of Affymetrix GeneChip data at the probe level</article-title><source>Bioinformatics</source><year>2004</year><volume>20</volume><issue>3</issue><fpage>307</fpage><lpage>15</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/btg405</pub-id><pub-id pub-id-type="pmid">14960456</pub-id></element-citation></ref><ref id="CR23"><label>23.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gottwald</surname><given-names>TR</given-names></name></person-group><article-title>Current epidemiological understanding of citrus Huanglongbing</article-title><source>Annu Rev Phytopathol</source><year>2010</year><volume>48</volume><fpage>119</fpage><lpage>39</lpage><pub-id pub-id-type="doi">10.1146/annurev-phyto-073009-114418</pub-id><pub-id pub-id-type="pmid">20415578</pub-id></element-citation></ref><ref id="CR24"><label>24.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gottwald</surname><given-names>T</given-names></name><name><surname>Bassanezi</surname><given-names>R</given-names></name><name><surname>Amorim</surname><given-names>L</given-names></name><name><surname>Bergamin-Filho</surname><given-names>A</given-names></name></person-group><article-title>Spatial pattern analysis of citrus canker-infected plantings in São Paulo, Brazil, and augmentation of infection elicited by the Asian leafminer</article-title><source>Phytopathology</source><year>2007</year><volume>97</volume><issue>6</issue><fpage>674</fpage><lpage>83</lpage><pub-id pub-id-type="doi">10.1094/PHYTO-97-6-0674</pub-id><pub-id pub-id-type="pmid">18943598</pub-id></element-citation></ref><ref id="CR25"><label>25.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Halbert</surname><given-names>SE</given-names></name><name><surname>Manjunath</surname><given-names>KL</given-names></name></person-group><article-title>Asian citrus psyllids (<italic>Sternorrhyncha: Psyllidae</italic>) and greening disease of citrus: a literature review and assessment of risk in Florida</article-title><source>Fla Entomol</source><year>2004</year><volume>87</volume><issue>3</issue><fpage>330</fpage><lpage>53</lpage><pub-id pub-id-type="doi">10.1653/0015-4040(2004)087[0330:ACPSPA]2.0.CO;2</pub-id></element-citation></ref><ref id="CR26"><label>26.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hammond</surname><given-names>JP</given-names></name><name><surname>White</surname><given-names>PJ</given-names></name></person-group><article-title>Sugar signaling in root responses to low phosphorus availability</article-title><source>Plant Physiol</source><year>2011</year><volume>156</volume><issue>3</issue><fpage>1033</fpage><lpage>40</lpage><pub-id pub-id-type="doi">10.1104/pp.111.175380</pub-id><pub-id pub-id-type="pmid">21487049</pub-id></element-citation></ref><ref id="CR27"><label>27.</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Hodges</surname><given-names>AW</given-names></name><name><surname>Spreen</surname><given-names>TH</given-names></name></person-group><source>Economic impacts of citrus greening (HLB) in Florida, 2006/07–2010/11</source><year>2012</year><publisher-loc>Gainesville, FL, USA</publisher-loc><publisher-name>University of Florida, Institute of Food and Agricultural Science</publisher-name></element-citation></ref><ref id="CR28"><label>28.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Holton</surname><given-names>N</given-names></name><name><surname>Nekrasov</surname><given-names>V</given-names></name><name><surname>Ronald</surname><given-names>PC</given-names></name><name><surname>Zipfel</surname><given-names>C</given-names></name></person-group><article-title>The phylogenetically-related pattern recognition receptors EFR and XA21 recruit similar immune signaling components in monocots and dicots</article-title><source>PLoS Pathog</source><year>2015</year><volume>11</volume><issue>1</issue><pub-id pub-id-type="doi">10.1371/journal.ppat.1004602</pub-id><pub-id pub-id-type="pmid">25607985</pub-id></element-citation></ref><ref id="CR29"><label>29.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hong</surname><given-names>F</given-names></name><name><surname>Breitling</surname><given-names>R</given-names></name><name><surname>McEntee</surname><given-names>CW</given-names></name><name><surname>Wittner</surname><given-names>BS</given-names></name><name><surname>Nemhauser</surname><given-names>JL</given-names></name><name><surname>Chory</surname><given-names>J</given-names></name></person-group><article-title>RankProd: a bioconductor package for detecting differentially expressed genes in meta-analysis</article-title><source>Bioinformatics</source><year>2006</year><volume>22</volume><issue>22</issue><fpage>2825</fpage><lpage>7</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/btl476</pub-id><pub-id pub-id-type="pmid">16982708</pub-id></element-citation></ref><ref id="CR30"><label>30.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jagoueix</surname><given-names>S</given-names></name><name><surname>Bove</surname><given-names>JM</given-names></name><name><surname>Garnier</surname><given-names>M</given-names></name></person-group><article-title>The phloem-limited bacterium of greening disease of citrus is a member of the alpha subdivision of the Proteobacteria</article-title><source>Int J Syst Bacteriol</source><year>1994</year><volume>44</volume><issue>3</issue><fpage>379</fpage><lpage>86</lpage><pub-id pub-id-type="doi">10.1099/00207713-44-3-379</pub-id><pub-id pub-id-type="pmid">7520729</pub-id></element-citation></ref><ref id="CR31"><label>31.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname><given-names>J</given-names></name><name><surname>Sagaram</surname><given-names>US</given-names></name><name><surname>Burns</surname><given-names>JK</given-names></name><name><surname>Li</surname><given-names>J</given-names></name><name><surname>Wang</surname><given-names>N</given-names></name></person-group><article-title>Response of sweet orange (<italic>Citrus sinensis</italic>) to <italic>‘Candidatus Liberibacter asiaticus’</italic> infection: microscopy and microarray analyses</article-title><source>Phytopathology</source><year>2009</year><volume>99</volume><issue>1</issue><fpage>50</fpage><lpage>7</lpage><pub-id pub-id-type="doi">10.1094/PHYTO-99-1-0050</pub-id><pub-id pub-id-type="pmid">19055434</pub-id></element-citation></ref><ref id="CR32"><label>32.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Landi</surname><given-names>L</given-names></name><name><surname>Romanazzi</surname><given-names>G</given-names></name></person-group><article-title>Seasonal variation of defense-related gene expression in leaves from Bois noir affected and recovered grapevines</article-title><source>J Agric Food Chem</source><year>2011</year><volume>59</volume><issue>12</issue><fpage>6628</fpage><lpage>37</lpage><pub-id pub-id-type="doi">10.1021/jf104297n</pub-id><pub-id pub-id-type="pmid">21627069</pub-id></element-citation></ref><ref id="CR33"><label>33.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Langfelder</surname><given-names>P</given-names></name><name><surname>Horvath</surname><given-names>S</given-names></name></person-group><article-title>WGCNA: an R package for weighted correlation network analysis</article-title><source>BMC Bioinformatics</source><year>2008</year><volume>9</volume><fpage>559</fpage><pub-id pub-id-type="doi">10.1186/1471-2105-9-559</pub-id><pub-id pub-id-type="pmid">19114008</pub-id></element-citation></ref><ref id="CR34"><label>34.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liao</surname><given-names>HL</given-names></name><name><surname>Burns</surname><given-names>JK</given-names></name></person-group><article-title>Gene expression in <italic>Citrus sinensis</italic> fruit tissues harvested from Huanglongbing-infected trees: comparison with girdled fruit</article-title><source>J Exp Bot</source><year>2012</year><volume>63</volume><issue>8</issue><fpage>3307</fpage><lpage>19</lpage><pub-id pub-id-type="doi">10.1093/jxb/ers070</pub-id><pub-id pub-id-type="pmid">22407645</pub-id></element-citation></ref><ref id="CR35"><label>35.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname><given-names>Y</given-names></name><name><surname>Heying</surname><given-names>E</given-names></name><name><surname>Tanumihardjo</surname><given-names>SA</given-names></name></person-group><article-title>History, global distribution, and nutritional importance of citrus fruits</article-title><source>Compr Rev Food Sci F</source><year>2012</year><volume>11</volume><issue>6</issue><fpage>530</fpage><lpage>45</lpage><pub-id pub-id-type="doi">10.1111/j.1541-4337.2012.00201.x</pub-id></element-citation></ref><ref id="CR36"><label>36.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lopes</surname><given-names>S</given-names></name><name><surname>Bertolini</surname><given-names>E</given-names></name><name><surname>Frare</surname><given-names>G</given-names></name><name><surname>Martins</surname><given-names>E</given-names></name><name><surname>Wulff</surname><given-names>N</given-names></name><name><surname>Teixeira</surname><given-names>D</given-names></name><etal></etal></person-group><article-title>Graft transmission efficiencies and multiplication of <italic>‘Candidatus Liberibacter americanus’</italic> and <italic>‘Ca. Liberibacter asiaticus’</italic> in citrus plants</article-title><source>Phytopathology</source><year>2009</year><volume>99</volume><issue>3</issue><fpage>301</fpage><lpage>6</lpage><pub-id pub-id-type="doi">10.1094/PHYTO-99-3-0301</pub-id><pub-id pub-id-type="pmid">19203283</pub-id></element-citation></ref><ref id="CR37"><label>37.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lu</surname><given-names>CA</given-names></name><name><surname>Ho</surname><given-names>TH</given-names></name><name><surname>Ho</surname><given-names>SL</given-names></name><name><surname>Yu</surname><given-names>SM</given-names></name></person-group><article-title>Three novel MYB proteins with one DNA binding repeat mediate sugar and hormone regulation of alpha-amylase gene expression</article-title><source>Plant Cell</source><year>2002</year><volume>14</volume><issue>8</issue><fpage>1963</fpage><lpage>80</lpage><pub-id pub-id-type="doi">10.1105/tpc.001735</pub-id><pub-id pub-id-type="pmid">12172034</pub-id></element-citation></ref><ref id="CR38"><label>38.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lu</surname><given-names>YT</given-names></name><name><surname>Li</surname><given-names>MY</given-names></name><name><surname>Cheng</surname><given-names>KT</given-names></name><name><surname>Tan</surname><given-names>CM</given-names></name><name><surname>Su</surname><given-names>LW</given-names></name><name><surname>Lin</surname><given-names>WY</given-names></name><etal></etal></person-group><article-title>Transgenic plants that express the phytoplasma effector SAP11 show altered phosphate starvation and defense responses</article-title><source>Plant Physiol</source><year>2014</year><volume>164</volume><issue>3</issue><fpage>1456</fpage><lpage>69</lpage><pub-id pub-id-type="doi">10.1104/pp.113.229740</pub-id><pub-id pub-id-type="pmid">24464367</pub-id></element-citation></ref><ref id="CR39"><label>39.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mafra</surname><given-names>V</given-names></name><name><surname>Martins</surname><given-names>PK</given-names></name><name><surname>Francisco</surname><given-names>CS</given-names></name><name><surname>Ribeiro-Alves</surname><given-names>M</given-names></name><name><surname>Freitas-Astua</surname><given-names>J</given-names></name><name><surname>Machado</surname><given-names>MA</given-names></name></person-group><article-title><italic>Candidatus Liberibacter americanus’</italic> induces significant reprogramming of the transcriptome of the susceptible citrus genotype</article-title><source>BMC Genomics</source><year>2013</year><volume>14</volume><fpage>247</fpage><pub-id pub-id-type="doi">10.1186/1471-2164-14-247</pub-id><pub-id pub-id-type="pmid">23586643</pub-id></element-citation></ref><ref id="CR40"><label>40.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Martinelli</surname><given-names>F</given-names></name><name><surname>Reagan</surname><given-names>RL</given-names></name><name><surname>Uratsu</surname><given-names>SL</given-names></name><name><surname>Phu</surname><given-names>ML</given-names></name><name><surname>Albrecht</surname><given-names>U</given-names></name><name><surname>Zhao</surname><given-names>W</given-names></name><etal></etal></person-group><article-title>Gene regulatory networks elucidating Huanglongbing disease mechanisms</article-title><source>PloS One</source><year>2013</year><volume>8</volume><issue>9</issue><pub-id pub-id-type="doi">10.1371/journal.pone.0074256</pub-id><pub-id pub-id-type="pmid">24086326</pub-id></element-citation></ref><ref id="CR41"><label>41.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Martinelli</surname><given-names>F</given-names></name><name><surname>Uratsu</surname><given-names>SL</given-names></name><name><surname>Albrecht</surname><given-names>U</given-names></name><name><surname>Reagan</surname><given-names>RL</given-names></name><name><surname>Phu</surname><given-names>ML</given-names></name><name><surname>Britton</surname><given-names>M</given-names></name><etal></etal></person-group><article-title>Transcriptome profiling of citrus fruit response to Huanglongbing disease</article-title><source>PloS One</source><year>2012</year><volume>7</volume><issue>5</issue><pub-id pub-id-type="doi">10.1371/journal.pone.0038039</pub-id><pub-id pub-id-type="pmid">22675433</pub-id></element-citation></ref><ref id="CR42"><label>42.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nazri</surname><given-names>A</given-names></name><name><surname>Lio</surname><given-names>P</given-names></name></person-group><article-title>Investigating meta-approaches for reconstructing gene networks in a mammalian cellular context</article-title><source>PloS One</source><year>2012</year><volume>7</volume><issue>1</issue><pub-id pub-id-type="doi">10.1371/journal.pone.0028713</pub-id><pub-id pub-id-type="pmid">22253694</pub-id></element-citation></ref><ref id="CR43"><label>43.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nwugo</surname><given-names>CC</given-names></name><name><surname>Duan</surname><given-names>Y</given-names></name><name><surname>Lin</surname><given-names>H</given-names></name></person-group><article-title>Study on citrus response to Huanglongbing highlights a down-regulation of defense-related proteins in lemon plants upon ‘<italic>Ca. Liberibacter asiaticus’</italic> infection</article-title><source>PloS One</source><year>2013</year><volume>8</volume><issue>6</issue><fpage>e67442</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0067442</pub-id><pub-id pub-id-type="pmid">23922636</pub-id></element-citation></ref><ref id="CR44"><label>44.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Peumans</surname><given-names>WJ</given-names></name><name><surname>Van Damme</surname><given-names>EJ</given-names></name></person-group><article-title>Lectins as plant defense proteins</article-title><source>Plant Physiol</source><year>1995</year><volume>109</volume><issue>2</issue><fpage>347</fpage><lpage>52</lpage><pub-id pub-id-type="doi">10.1104/pp.109.2.347</pub-id><pub-id pub-id-type="pmid">7480335</pub-id></element-citation></ref><ref id="CR45"><label>45.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Plank</surname><given-names>J</given-names></name><name><surname>Siebenhofer</surname><given-names>A</given-names></name><name><surname>Berghold</surname><given-names>A</given-names></name><name><surname>Jeitler</surname><given-names>K</given-names></name><name><surname>Horvath</surname><given-names>K</given-names></name><name><surname>Mrak</surname><given-names>P</given-names></name><etal></etal></person-group><article-title>Systematic review and meta-analysis of short-acting insulin analogues in patients with diabetes mellitus</article-title><source>Arch Intern Med</source><year>2005</year><volume>165</volume><issue>12</issue><fpage>1337</fpage><lpage>44</lpage><pub-id pub-id-type="doi">10.1001/archinte.165.12.1337</pub-id><pub-id pub-id-type="pmid">15983281</pub-id></element-citation></ref><ref id="CR46"><label>46.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Puniya</surname><given-names>BL</given-names></name><name><surname>Kulshreshtha</surname><given-names>D</given-names></name><name><surname>Verma</surname><given-names>SP</given-names></name><name><surname>Kumar</surname><given-names>S</given-names></name><name><surname>Ramachandran</surname><given-names>S</given-names></name></person-group><article-title>Integrated gene co-expression network analysis in the growth phase of <italic>Mycobacterium tuberculosis</italic> reveals new potential drug targets</article-title><source>Mol BioSyst</source><year>2013</year><volume>9</volume><issue>11</issue><fpage>2798</fpage><lpage>815</lpage><pub-id pub-id-type="doi">10.1039/c3mb70278b</pub-id><pub-id pub-id-type="pmid">24056838</pub-id></element-citation></ref><ref id="CR47"><label>47.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pustika</surname><given-names>A</given-names></name><name><surname>Subandiyah</surname><given-names>S</given-names></name><name><surname>Holford</surname><given-names>P</given-names></name><name><surname>Beattie</surname><given-names>G</given-names></name><name><surname>Iwanami</surname><given-names>T</given-names></name><name><surname>Masaoka</surname><given-names>Y</given-names></name></person-group><article-title>Interactions between plant nutrition and symptom expression in mandarin trees infected with the disease Huanglongbing</article-title><source>Australas Plant Dis Notes</source><year>2008</year><volume>3</volume><issue>1</issue><fpage>112</fpage><lpage>5</lpage><pub-id pub-id-type="doi">10.1071/DN08045</pub-id></element-citation></ref><ref id="CR48"><label>48.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rawat</surname><given-names>N</given-names></name><name><surname>Neeraja</surname><given-names>CN</given-names></name><name><surname>Nair</surname><given-names>S</given-names></name><name><surname>Bentur</surname><given-names>JS</given-names></name></person-group><article-title>Differential gene expression in gall midge- susceptible rice genotypes revealed by suppressive subtraction hybridization (SSH) cDNA libraries and microarray analysis</article-title><source>Rice</source><year>2012</year><volume>5</volume><issue>1</issue><fpage>1</fpage><lpage>15</lpage><pub-id pub-id-type="doi">10.1186/1939-8433-5-8</pub-id><pub-id pub-id-type="pmid">24764501</pub-id></element-citation></ref><ref id="CR49"><label>49.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sarwar</surname><given-names>M</given-names></name></person-group><article-title>Management of rice stem borers (<italic>Lepidoptera: Pyralidae</italic>) through host plant resistance in early, medium and late plantings of rice (<italic>Oryza sativa</italic> L.)</article-title><source>J Cereals Oilseeds</source><year>2012</year><volume>3</volume><issue>1</issue><fpage>10</fpage><lpage>4</lpage></element-citation></ref><ref id="CR50"><label>50.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shaik</surname><given-names>R</given-names></name><name><surname>Ramakrishna</surname><given-names>W</given-names></name></person-group><article-title>Machine learning approaches distinguish multiple stress conditions using stress-responsive genes and identify candidate genes for broad resistance in rice</article-title><source>Plant Physiol</source><year>2014</year><volume>164</volume><issue>1</issue><fpage>481</fpage><lpage>95</lpage><pub-id pub-id-type="doi">10.1104/pp.113.225862</pub-id><pub-id pub-id-type="pmid">24235132</pub-id></element-citation></ref><ref id="CR51"><label>51.</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Smyth</surname><given-names>GK</given-names></name></person-group><article-title>Limma: linear models for microarray data</article-title><source>Bioinformatics and computational biology solutions using R and Bioconductor</source><year>2005</year><fpage>397</fpage><lpage>420</lpage></element-citation></ref><ref id="CR52"><label>52.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sugio</surname><given-names>A</given-names></name><name><surname>MacLean</surname><given-names>AM</given-names></name><name><surname>Hogenhout</surname><given-names>SA</given-names></name></person-group><article-title>The small phytoplasma virulence effector SAP11 contains distinct domains required for nuclear targeting and CIN‐TCP binding and destabilization</article-title><source>New Phytol</source><year>2014</year><volume>202</volume><issue>3</issue><fpage>838</fpage><lpage>48</lpage><pub-id pub-id-type="doi">10.1111/nph.12721</pub-id><pub-id pub-id-type="pmid">24552625</pub-id></element-citation></ref><ref id="CR53"><label>53.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Taheri</surname><given-names>F</given-names></name><name><surname>Nematzadeh</surname><given-names>G</given-names></name><name><surname>Zamharir</surname><given-names>MG</given-names></name><name><surname>Nekouei</surname><given-names>MK</given-names></name><name><surname>Naghavi</surname><given-names>M</given-names></name><name><surname>Mardi</surname><given-names>M</given-names></name><etal></etal></person-group><article-title>Proteomic analysis of the Mexican lime tree response to <italic>‘Candidatus Phytoplasma aurantifolia’</italic> infection</article-title><source>Mol BioSyst</source><year>2011</year><volume>7</volume><issue>11</issue><fpage>3028</fpage><lpage>35</lpage><pub-id pub-id-type="doi">10.1039/c1mb05268c</pub-id><pub-id pub-id-type="pmid">21853195</pub-id></element-citation></ref><ref id="CR54"><label>54.</label><mixed-citation publication-type="other">Texeira DC, Ayres J, Kitajima E, Danet L, Jagoueix-Eveillard S, Saillard C, et al. First report of a Huanglongbing-like disease of citrus in Sao Paulo State, Brazil and association of a new Liberibacter species, <italic>‘Candidatus Liberibacter americanus’</italic>, with the disease. Plant Dis. 2005;89(1):107–7.</mixed-citation></ref><ref id="CR55"><label>55.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Thirunavukkarasu</surname><given-names>N</given-names></name><name><surname>Hossain</surname><given-names>F</given-names></name><name><surname>Mohan</surname><given-names>S</given-names></name><name><surname>Shiriga</surname><given-names>K</given-names></name><name><surname>Mittal</surname><given-names>S</given-names></name><name><surname>Sharma</surname><given-names>R</given-names></name><etal></etal></person-group><article-title>Genome-wide expression of transcriptomes and their co-expression pattern in subtropical maize (<italic>Zea mays</italic> L.) under waterlogging stress</article-title><source>PloS One</source><year>2013</year><volume>8</volume><issue>8</issue><fpage>e70433</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0070433</pub-id><pub-id pub-id-type="pmid">23936429</pub-id></element-citation></ref><ref id="CR56"><label>56.</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Ting</surname><given-names>SV</given-names></name></person-group><person-group person-group-type="editor"><name><surname>Nagy</surname><given-names>S</given-names></name><name><surname>Attaway</surname><given-names>JA</given-names></name></person-group><article-title>Nutrients and nutrition of citrus fruits</article-title><source>Citrus Nutrition and Quality</source><year>1980</year><publisher-loc>Washington DC</publisher-loc><publisher-name>American Chemical Society</publisher-name><fpage>3</fpage><lpage>24</lpage></element-citation></ref><ref id="CR57"><label>57.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tsukuda</surname><given-names>S</given-names></name><name><surname>Gomi</surname><given-names>K</given-names></name><name><surname>Yamamoto</surname><given-names>H</given-names></name><name><surname>Akimitsu</surname><given-names>K</given-names></name></person-group><article-title>Characterization of cDNAs encoding two distinct miraculin-like proteins and stress-related modulation of the corresponding mRNAs in <italic>Citrus jambhiri</italic> Lush</article-title><source>Plant Mol Biol</source><year>2006</year><volume>60</volume><issue>1</issue><fpage>125</fpage><lpage>36</lpage><pub-id pub-id-type="doi">10.1007/s11103-005-2941-4</pub-id><pub-id pub-id-type="pmid">16463104</pub-id></element-citation></ref><ref id="CR58"><label>58.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Usadel</surname><given-names>B</given-names></name><name><surname>Nagel</surname><given-names>A</given-names></name><name><surname>Steinhauser</surname><given-names>D</given-names></name><name><surname>Gibon</surname><given-names>Y</given-names></name><name><surname>Blasing</surname><given-names>OE</given-names></name><name><surname>Redestig</surname><given-names>H</given-names></name><etal></etal></person-group><article-title>PageMan: an interactive ontology tool to generate, display, and annotate overview graphs for profiling experiments</article-title><source>BMC Bioinformatics</source><year>2006</year><volume>7</volume><fpage>535</fpage><pub-id pub-id-type="doi">10.1186/1471-2105-7-535</pub-id><pub-id pub-id-type="pmid">17176458</pub-id></element-citation></ref><ref id="CR59"><label>59.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Usadel</surname><given-names>B</given-names></name><name><surname>Nagel</surname><given-names>A</given-names></name><name><surname>Thimm</surname><given-names>O</given-names></name><name><surname>Redestig</surname><given-names>H</given-names></name><name><surname>Blaesing</surname><given-names>OE</given-names></name><name><surname>Palacios-Rojas</surname><given-names>N</given-names></name><etal></etal></person-group><article-title>Extension of the visualization tool MapMan to allow statistical analysis of arrays, display of corresponding genes, and comparison with known responses</article-title><source>Plant Physiol</source><year>2005</year><volume>138</volume><issue>3</issue><fpage>1195</fpage><lpage>204</lpage><pub-id pub-id-type="doi">10.1104/pp.105.060459</pub-id><pub-id pub-id-type="pmid">16009995</pub-id></element-citation></ref><ref id="CR60"><label>60.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Van Loon</surname><given-names>L</given-names></name><name><surname>Bakker</surname><given-names>P</given-names></name><name><surname>Pieterse</surname><given-names>C</given-names></name></person-group><article-title>Systemic resistance induced by rhizosphere bacteria</article-title><source>Annu Rev Phytopathol</source><year>1998</year><volume>36</volume><issue>1</issue><fpage>453</fpage><lpage>83</lpage><pub-id pub-id-type="doi">10.1146/annurev.phyto.36.1.453</pub-id><pub-id pub-id-type="pmid">15012509</pub-id></element-citation></ref><ref id="CR61"><label>61.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>N</given-names></name><name><surname>Trivedi</surname><given-names>P</given-names></name></person-group><article-title>Citrus Huanglongbing: a newly relevant disease presents unprecedented challenges</article-title><source>Phytopathology</source><year>2013</year><volume>103</volume><issue>7</issue><fpage>652</fpage><lpage>65</lpage><pub-id pub-id-type="doi">10.1094/PHYTO-12-12-0331-RVW</pub-id><pub-id pub-id-type="pmid">23441969</pub-id></element-citation></ref><ref id="CR62"><label>62.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>M</given-names></name><name><surname>Zheng</surname><given-names>Q</given-names></name><name><surname>Shen</surname><given-names>Q</given-names></name><name><surname>Guo</surname><given-names>S</given-names></name></person-group><article-title>The critical role of potassium in plant stress response</article-title><source>Int J Mol Sci</source><year>2013</year><volume>14</volume><issue>4</issue><fpage>7370</fpage><lpage>90</lpage><pub-id pub-id-type="doi">10.3390/ijms14047370</pub-id><pub-id pub-id-type="pmid">23549270</pub-id></element-citation></ref><ref id="CR63"><label>63.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Weintraub</surname><given-names>PG</given-names></name></person-group><article-title>Insect vectors of phytoplasmas and their control-an update</article-title><source>B Insectol</source><year>2007</year><volume>60</volume><issue>2</issue><fpage>169</fpage><lpage>73</lpage></element-citation></ref><ref id="CR64"><label>64.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Xia</surname><given-names>Y</given-names></name><name><surname>Suzuki</surname><given-names>H</given-names></name><name><surname>Borevitz</surname><given-names>J</given-names></name><name><surname>Blount</surname><given-names>J</given-names></name><name><surname>Guo</surname><given-names>Z</given-names></name><name><surname>Patel</surname><given-names>K</given-names></name><etal></etal></person-group><article-title>An extracellular aspartic protease functions in <italic>Arabidopsis</italic> disease resistance signaling</article-title><source>EMBO J</source><year>2004</year><volume>23</volume><issue>4</issue><fpage>980</fpage><lpage>8</lpage><pub-id pub-id-type="doi">10.1038/sj.emboj.7600086</pub-id><pub-id pub-id-type="pmid">14765119</pub-id></element-citation></ref><ref id="CR65"><label>65.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yip</surname><given-names>AM</given-names></name><name><surname>Horvath</surname><given-names>S</given-names></name></person-group><article-title>Gene network interconnectedness and the generalized topological overlap measure</article-title><source>BMC Bioinformatics</source><year>2007</year><volume>8</volume><fpage>22</fpage><pub-id pub-id-type="doi">10.1186/1471-2105-8-22</pub-id><pub-id pub-id-type="pmid">17250769</pub-id></element-citation></ref><ref id="CR66"><label>66.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname><given-names>B</given-names></name><name><surname>Horvath</surname><given-names>S</given-names></name></person-group><article-title>A general framework for weighted gene co-expression network analysis</article-title><source>Statistical applications in genetics and molecular biology</source><year>2005</year><volume>4</volume><fpage>1128</fpage><pub-id pub-id-type="doi">10.2202/1544-6115.1128</pub-id></element-citation></ref><ref id="CR67"><label>67.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname><given-names>Z</given-names></name><name><surname>Yu</surname><given-names>J</given-names></name><name><surname>Li</surname><given-names>D</given-names></name><name><surname>Zhang</surname><given-names>Z</given-names></name><name><surname>Liu</surname><given-names>F</given-names></name><name><surname>Zhou</surname><given-names>X</given-names></name><etal></etal></person-group><article-title>PMRD: plant microRNA database</article-title><source>Nucleic Acids Res</source><year>2010</year><volume>38</volume><fpage>D806</fpage><lpage>13</lpage><pub-id pub-id-type="doi">10.1093/nar/gkp818</pub-id><pub-id pub-id-type="pmid">19808935</pub-id></element-citation></ref><ref id="CR68"><label>68.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhao</surname><given-names>H</given-names></name><name><surname>Sun</surname><given-names>R</given-names></name><name><surname>Albrecht</surname><given-names>U</given-names></name><name><surname>Padmanabhan</surname><given-names>C</given-names></name><name><surname>Wang</surname><given-names>A</given-names></name><name><surname>Coffey</surname><given-names>MD</given-names></name><etal></etal></person-group><article-title>Small RNA profiling reveals phosphorus deficiency as a contributing factor in symptom expression for citrus Huanglongbing disease</article-title><source>Mol Plant</source><year>2013</year><volume>6</volume><issue>2</issue><fpage>301</fpage><lpage>10</lpage><pub-id pub-id-type="doi">10.1093/mp/sst002</pub-id><pub-id pub-id-type="pmid">23292880</pub-id></element-citation></ref><ref id="CR69"><label>69.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zheng</surname><given-names>ZL</given-names></name><name><surname>Zhao</surname><given-names>Y</given-names></name></person-group><article-title>Transcriptome comparison and gene coexpression network analysis provide a systems view of citrus response to <italic>‘Candidatus Liberibacter asiaticus’</italic> infection</article-title><source>BMC Genomics</source><year>2013</year><volume>14</volume><fpage>27</fpage><pub-id pub-id-type="doi">10.1186/1471-2164-14-27</pub-id><pub-id pub-id-type="pmid">23324561</pub-id></element-citation></ref><ref id="CR70"><label>70.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhong</surname><given-names>R</given-names></name><name><surname>Demura</surname><given-names>T</given-names></name><name><surname>Ye</surname><given-names>ZH</given-names></name></person-group><article-title>SND1, a NAC domain transcription factor, is a key regulator of secondary wall synthesis in fibers of <italic>Arabidopsis</italic></article-title><source>Plant Cell</source><year>2006</year><volume>18</volume><issue>11</issue><fpage>3158</fpage><lpage>70</lpage><pub-id pub-id-type="doi">10.1105/tpc.106.047399</pub-id><pub-id pub-id-type="pmid">17114348</pub-id></element-citation></ref><ref id="CR71"><label>71.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zou</surname><given-names>H</given-names></name><name><surname>Gowda</surname><given-names>S</given-names></name><name><surname>Zhou</surname><given-names>L</given-names></name><name><surname>Hajeri</surname><given-names>S</given-names></name><name><surname>Chen</surname><given-names>G</given-names></name><name><surname>Duan</surname><given-names>Y</given-names></name></person-group><article-title>The destructive citrus pathogen, ‘<italic>Candidatus Liberibacter asiaticus’</italic> encodes a functional flagellin characteristic of a pathogen-associated molecular pattern</article-title><source>PloS One</source><year>2012</year><volume>7</volume><issue>9</issue><pub-id pub-id-type="doi">10.1371/journal.pone.0046447</pub-id><pub-id pub-id-type="pmid">23029520</pub-id></element-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Bois/explor/OrangerV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000631 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000631 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Bois
|area= OrangerV1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé= PMC:4517500
|texte= Comprehensive meta-analysis, co-expression, and miRNA nested network analysis identifies gene candidates in citrus against Huanglongbing disease
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:26215595" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \
| NlmPubMed2Wicri -a OrangerV1
| This area was generated with Dilib version V0.6.25. |  |