Construction of citrus gene coexpression networks from microarray data using random matrix theory
Identifieur interne : 000440 ( Pmc/Corpus ); précédent : 000439; suivant : 000441Construction of citrus gene coexpression networks from microarray data using random matrix theory
Auteurs : Dongliang Du ; Nidhi Rawat ; Zhanao Deng ; Fred G. GmitterSource :
- Horticulture Research [ 2052-7276 ] ; 2015.
Abstract
After the sequencing of citrus genomes, gene function annotation is becoming a new challenge. Gene coexpression analysis can be employed for function annotation using publicly available microarray data sets. In this study, 230 sweet orange (
Url:
DOI: 10.1038/hortres.2015.26
PubMed: 26504573
PubMed Central: 4595991
Links to Exploration step
PMC:4595991Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Construction of citrus gene coexpression networks from microarray data using random matrix theory</title>
<author><name sortKey="Du, Dongliang" sort="Du, Dongliang" uniqKey="Du D" first="Dongliang" last="Du">Dongliang Du</name>
<affiliation><nlm:aff id="aff1"><institution>Citrus Research and Education Center, Institute of Food and Agricultural Sciences, University of Florida</institution>
, Lake Alfred, FL 33850,<country>USA</country>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Rawat, Nidhi" sort="Rawat, Nidhi" uniqKey="Rawat N" first="Nidhi" last="Rawat">Nidhi Rawat</name>
<affiliation><nlm:aff id="aff2"><institution>Gulf Coast Research and Education Center, Institute of Food and Agricultural Sciences, University of Florida</institution>
, Wimauma, FL 33598,<country>USA</country>
</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="aff2"><institution>Gulf Coast Research and Education Center, Institute of Food and Agricultural Sciences, University of Florida</institution>
, Wimauma, FL 33598,<country>USA</country>
</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="aff1"><institution>Citrus Research and Education Center, Institute of Food and Agricultural Sciences, University of Florida</institution>
, Lake Alfred, FL 33850,<country>USA</country>
</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">26504573</idno>
<idno type="pmc">4595991</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4595991</idno>
<idno type="RBID">PMC:4595991</idno>
<idno type="doi">10.1038/hortres.2015.26</idno>
<date when="2015">2015</date>
<idno type="wicri:Area/Pmc/Corpus">000440</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Construction of citrus gene coexpression networks from microarray data using random matrix theory</title>
<author><name sortKey="Du, Dongliang" sort="Du, Dongliang" uniqKey="Du D" first="Dongliang" last="Du">Dongliang Du</name>
<affiliation><nlm:aff id="aff1"><institution>Citrus Research and Education Center, Institute of Food and Agricultural Sciences, University of Florida</institution>
, Lake Alfred, FL 33850,<country>USA</country>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Rawat, Nidhi" sort="Rawat, Nidhi" uniqKey="Rawat N" first="Nidhi" last="Rawat">Nidhi Rawat</name>
<affiliation><nlm:aff id="aff2"><institution>Gulf Coast Research and Education Center, Institute of Food and Agricultural Sciences, University of Florida</institution>
, Wimauma, FL 33598,<country>USA</country>
</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="aff2"><institution>Gulf Coast Research and Education Center, Institute of Food and Agricultural Sciences, University of Florida</institution>
, Wimauma, FL 33598,<country>USA</country>
</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="aff1"><institution>Citrus Research and Education Center, Institute of Food and Agricultural Sciences, University of Florida</institution>
, Lake Alfred, FL 33850,<country>USA</country>
</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">Horticulture Research</title>
<idno type="eISSN">2052-7276</idno>
<imprint><date when="2015">2015</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p>After the sequencing of citrus genomes, gene function annotation is becoming a new challenge. Gene coexpression analysis can be employed for function annotation using publicly available microarray data sets. In this study, 230 sweet orange (<italic>Citrus sinensis</italic>
) microarrays were used to construct seven coexpression networks, including one condition-independent and six condition-dependent (Citrus canker, Huanglongbing, leaves, flavedo, albedo, and flesh) networks. In total, these networks contain 37 633 edges among 6256 nodes (genes), which accounts for 52.11% measurable genes of the citrus microarray. Then, these networks were partitioned into functional modules using the Markov Cluster Algorithm. Significantly enriched Gene Ontology biological process terms and KEGG pathway terms were detected for 343 and 60 modules, respectively. Finally, independent verification of these networks was performed using another expression data of 371 genes. This study provides new targets for further functional analyses in citrus.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></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">Hortic Res</journal-id>
<journal-id journal-id-type="iso-abbrev">Hortic Res</journal-id>
<journal-title-group><journal-title>Horticulture Research</journal-title>
</journal-title-group>
<issn pub-type="epub">2052-7276</issn>
<publisher><publisher-name>Nature Publishing Group</publisher-name>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">26504573</article-id>
<article-id pub-id-type="pmc">4595991</article-id>
<article-id pub-id-type="pii">hortres201526</article-id>
<article-id pub-id-type="doi">10.1038/hortres.2015.26</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Construction of citrus gene coexpression networks from microarray data using random matrix theory</article-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Du</surname>
<given-names>Dongliang</given-names>
</name>
<xref ref-type="aff" rid="aff1">1</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Rawat</surname>
<given-names>Nidhi</given-names>
</name>
<xref ref-type="aff" rid="aff2">2</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Deng</surname>
<given-names>Zhanao</given-names>
</name>
<xref ref-type="aff" rid="aff2">2</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Gmitter</surname>
<given-names>Fred G.</given-names>
<suffix>Jr.</suffix>
</name>
<xref ref-type="aff" rid="aff1">1</xref>
<xref ref-type="corresp" rid="caf1">*</xref>
</contrib>
<aff id="aff1"><label>1</label>
<institution>Citrus Research and Education Center, Institute of Food and Agricultural Sciences, University of Florida</institution>
, Lake Alfred, FL 33850,<country>USA</country>
</aff>
<aff id="aff2"><label>2</label>
<institution>Gulf Coast Research and Education Center, Institute of Food and Agricultural Sciences, University of Florida</institution>
, Wimauma, FL 33598,<country>USA</country>
</aff>
</contrib-group>
<author-notes><corresp id="caf1"><label>*</label>
(<email>fgmitter@ufl.edu</email>
)</corresp>
</author-notes>
<pub-date pub-type="epub"><day>10</day>
<month>06</month>
<year>2015</year>
</pub-date>
<pub-date pub-type="collection"><year>2015</year>
</pub-date>
<volume>2</volume>
<fpage>15026</fpage>
<lpage></lpage>
<history><date date-type="received"><day>30</day>
<month>03</month>
<year>2015</year>
</date>
<date date-type="rev-recd"><day>29</day>
<month>04</month>
<year>2015</year>
</date>
<date date-type="accepted"><day>04</day>
<month>05</month>
<year>2015</year>
</date>
</history>
<permissions><copyright-statement>Copyright © 2015 Nanjing Agricultural University</copyright-statement>
<copyright-year>2015</copyright-year>
<copyright-holder>Nanjing Agricultural University</copyright-holder>
<license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by-nc-sa/3.0/"><pmc-comment>author-paid</pmc-comment>
<license-p>This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by-nc-sa/3.0/">http://creativecommons.org/licenses/by-nc-sa/3.0/</ext-link>
</license-p>
</license>
</permissions>
<abstract><p>After the sequencing of citrus genomes, gene function annotation is becoming a new challenge. Gene coexpression analysis can be employed for function annotation using publicly available microarray data sets. In this study, 230 sweet orange (<italic>Citrus sinensis</italic>
) microarrays were used to construct seven coexpression networks, including one condition-independent and six condition-dependent (Citrus canker, Huanglongbing, leaves, flavedo, albedo, and flesh) networks. In total, these networks contain 37 633 edges among 6256 nodes (genes), which accounts for 52.11% measurable genes of the citrus microarray. Then, these networks were partitioned into functional modules using the Markov Cluster Algorithm. Significantly enriched Gene Ontology biological process terms and KEGG pathway terms were detected for 343 and 60 modules, respectively. Finally, independent verification of these networks was performed using another expression data of 371 genes. This study provides new targets for further functional analyses in citrus.</p>
</abstract>
</article-meta>
</front>
<body><sec sec-type="intro"><title>Introduction</title>
<p>Gene coexpression network analysis is an attractive method for gene function annotation, which has been used in many model organisms, including yeast, mouse, human, Arabidopsis, and grapevine<sup><xref ref-type="bibr" rid="bib1 bib2 bib3 bib4 bib5">1–5</xref>
</sup>
. In the gene coexpression network, nodes represent genes and edges represent significant correlations between the expression patterns of connected genes<sup><xref ref-type="bibr" rid="bib6">6</xref>
</sup>
. After network construction, highly connected genes are clustered into modules. Genes within one module tend to participate in similar biological processes. Therefore, the function of unannotated genes could be hypothesized based on “guilt-by-association” principle<sup><xref ref-type="bibr" rid="bib7">7</xref>
</sup>
.</p>
<p>After the sequencing of citrus genomes<sup><xref ref-type="bibr" rid="bib8">8</xref>
</sup>
, gene function annotation is becoming a new challenge. For citrus, large amounts of data from microarray and RNA-seq experiments are available in public databases<sup><xref ref-type="bibr" rid="bib9 bib10 bib11 bib12 bib13">9–13</xref>
</sup>
. These data make it possible to construct gene coexpression networks for citrus. Several papers on citrus gene coexpression networks have been published<sup><xref ref-type="bibr" rid="bib14 bib15 bib16 bib17">14–17</xref>
</sup>
. Most of these studies focused on specific areas and used small data sets. Only one study used 297 citrus microarrays, and covered the general area and several specific areas<sup><xref ref-type="bibr" rid="bib17">17</xref>
</sup>
. However, a limitation of this study was that probe sets were used to construct gene coexpression networks, not the genes, which were used in many coexpression studies<sup><xref ref-type="bibr" rid="bib18 bib19 bib20">18–20</xref>
</sup>
. There are also some protein–protein interaction (PPI) networks, but these networks were inferred based on PPI networks of Arabidopsis<sup><xref ref-type="bibr" rid="bib21 bib22 bib23">21–23</xref>
</sup>
.</p>
<p>In this study, we first made a customized Chip Definition File (CDF) by AffyProbeMiner to transfer probes to gene locus. Then, seven gene coexpression networks were constructed by RMTGeneNet using all or part of 230 citrus microarrays. These networks were partitioned into modules, and the functional coherence of modules was assessed by Gene Ontology (GO) and KEGG pathway enrichment analyses. Finally, RNA-seq data of 371 genes were used to test the validity of these networks.</p>
</sec>
<sec sec-type="materials|methods"><title>Materials and methods</title>
<sec><title>Data collection and preprocessing</title>
<p>The sweet orange (<italic>Citrus sinensis</italic>
) microarray data used in this study were downloaded from National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO)<sup><xref ref-type="bibr" rid="bib24">24</xref>
</sup>
. A total of 231 CEL files were obtained from the platform GPL5731. The raw CEL data were preprocessed with RMA normalization using affy package of R 3.1.0<sup><xref ref-type="bibr" rid="bib25">25</xref>
</sup>
. One sample (GSM825502) that failed more than one test of arrayQualityMetrics<sup><xref ref-type="bibr" rid="bib26">26</xref>
</sup>
was removed and a total of 230 samples remained for network construction (<xref ref-type="supplementary-material" rid="xob2">Table S1</xref>
). Based on hierarchical cluster analysis, these 230 microarrays (called “all data”) were classified into “citrus canker” (30 arrays) and “HLB” (36 arrays), or “leaves” (63 arrays), “flavedo” (40 arrays), “albedo” (31 arrays), and “flesh” (43 arrays). These sub-data sets were also used for networks construction. To map the microarray probes to citrus genes, a customized CDF was generated by AffyProbeMiner<sup><xref ref-type="bibr" rid="bib27">27</xref>
</sup>
using <italic>C. x clementina</italic>
v1.0 annotation as the reference<sup><xref ref-type="bibr" rid="bib8">8</xref>
</sup>
. Probes mapping to multiple citrus loci and probe sets containing less than five members were discarded. Data about the NICCE networks were downloaded from <ext-link ext-link-type="uri" xlink:href="http://citrus.adelaide.edu.au/nicce/home.aspx">http://citrus.adelaide.edu.au/nicce/home.aspx</ext-link>
.</p>
</sec>
<sec><title>Coexpression network construction and topological analysis</title>
<p>The coexpression networks were constructed using the RMTGeneNet package<sup><xref ref-type="bibr" rid="bib28">28</xref>
</sup>
. A minimum of 25 input microarrays is required for this application. First, a gene expression correlation matrix was constructed using pair-wise Pearson correlation coefficients (PCC). Then, a threshold was determined according to the transition of nearest neighbor spacing distribution from the Gaussian distribution to Poisson distribution (<italic>p</italic>
= 0.001). Coexpression networks were visualized using the Cytoscape 2.8.3<sup><xref ref-type="bibr" rid="bib29">29</xref>
</sup>
, and all topological analyses were performed using the NetworkAnalyzer package<sup><xref ref-type="bibr" rid="bib30">30</xref>
</sup>
for Cytoscape 2.8.3.</p>
</sec>
<sec><title>Module clustering and functional enrichment analysis</title>
<p>The Markov Cluster (MCL) Algorithm<sup><xref ref-type="bibr" rid="bib31">31</xref>
</sup>
, an efficient graph clustering algorithm based on the simulation of random walk, was used to partition the network into modules. The inflation parameter (I) was scanned from 1.2 to 5.0 with increments of 0.2. Area fraction, mass fraction, and efficiency were used to determine the inflation parameter for MCL. The GO terms and Arabidopsis homologs for <italic>C. x clementina</italic>
genes were downloaded from Phytozome v10<sup><xref ref-type="bibr" rid="bib32">32</xref>
</sup>
. The KEGG annotations of Arabidopsis genes were obtained through KEGGREST of Bioconductor<sup><xref ref-type="bibr" rid="bib33">33</xref>
</sup>
. GO biological process term enrichment analysis was carried out using topGO package of Bioconductor<sup><xref ref-type="bibr" rid="bib33">33</xref>
</sup>
. KEGG pathway enrichment was performed in R 3.1.0. Terms enriched with a Fisher’s test <italic>p</italic>
-value <0.05 were considered.</p>
</sec>
<sec><title>Genome synteny analysis and network validation</title>
<p>To use the RNA-seq data of the <italic>C. sinensis</italic>
annotation project (CAP)<sup><xref ref-type="bibr" rid="bib34">34</xref>
</sup>
, synteny analysis between <italic>C. sinensis</italic>
genome and <italic>C. x clementina</italic>
genome was conducted locally using the similar method developed for the Plant Genome Duplication Database<sup><xref ref-type="bibr" rid="bib35">35</xref>
,<xref ref-type="bibr" rid="bib36">36</xref>
</sup>
. First, BLASTP<sup><xref ref-type="bibr" rid="bib37">37</xref>
</sup>
was conducted using all <italic>C. sinensis</italic>
proteins to search for potential anchors (<italic>E</italic>
< 1e−10, top 5 matches) in the <italic>C. x clementina</italic>
genome. Afterwards, MCscan was employed to identify homologous regions<sup><xref ref-type="bibr" rid="bib38">38</xref>
</sup>
. Finally, syntenic blocks were evaluated by ColinearScan<sup><xref ref-type="bibr" rid="bib39">39</xref>
</sup>
. Alignments with an <italic>E</italic>
value < 1e−10 were considered as significant matches. The expression data of 371 <italic>C. sinensis</italic>
genes were downloaded from CAP and the correlation coefficients between them were calculated using R 3.1.0.</p>
</sec>
</sec>
<sec sec-type="results"><title>Results</title>
<sec><title>Network construction</title>
<p>As shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>
, 231 Affymetrix citrus microarrays were downloaded from the NCBI GEO. After quality check, 230 high-quality microarrays (<xref ref-type="supplementary-material" rid="xob2">Table S1</xref>
) were chosen for downstream analyses. Based on the hierarchical cluster analysis (<xref ref-type="supplementary-material" rid="xob1">Figure S1</xref>
), these microarrays were first distributed into seven organ groups: flower, stem, leaves, fruit, seed, roots, and epicotyls. In the fruit group, they were further divided into flavedo, albedo, flesh, and vascular core (also called central core) subgroup. The data from albedo and vascular core were first clustered together and then clustered with data from other parts of fruit. This is reasonable considering that albedo and vascular core are composed of a colorless, spongy network of parenchymatous cells. These data sets were combined into one group, and labeled “albedo”, because neither was large enough for RMTGeneNet analysis. Five groups (flower, stem, seed, roots, and epicotyls), which had fewer microarrays than the minimum requirement of RMTGeneNet (<xref rid="tbl1" ref-type="table">Table 1</xref>
), were not included for condition-dependent coexpression analysis. Within groups, microarrays of the same treatment were clustered together. Two major diseases of citrus<sup><xref ref-type="bibr" rid="bib40">40</xref>
</sup>
, citrus canker and HLB, constituted 38.3% of the experiments (citrus canker: 30, HLB: 58), or 81.5% of the experiments if controls were not included. Other treatments were not included for network construction because of insufficient numbers of microarrays. For citrus canker, all microarrays are included in the leaves group. However, HLB data covered five groups (stem, leaves, fruit, seed, and roots). Only 36 microarrays in the fruit group were used for constructing “HLB” coexpression network. Therefore, these 230 microarrays (called “all data”) were divided into sub-data sets of “citrus canker”, “HLB”, “leaves”, “flavedo”, “albedo”, and “flesh” based on their experimental conditions or organ types. Data from these seven groups were analyzed individually to construct coexpression networks.</p>
<p>The Affymetrix citrus microarray contains 30 217 probe sets and 341 730 probes. In order to map the probes to citrus gene loci, a customized CDF was generated by AffyProbeMiner<sup><xref ref-type="bibr" rid="bib27">27</xref>
</sup>
using <italic>C. x clementina</italic>
v1.0 annotation as the reference<sup><xref ref-type="bibr" rid="bib8">8</xref>
</sup>
. After removing ambiguous probes mapping to multiple gene loci and probe sets with less than 5 probes, 158 557 probes belonging to 12 005 gene loci were kept in the customized CDF. Therefore, the expression of 48.9% citrus genes (12 005/24 533) can be tested accurately using the Affymetrix citrus microarray. According to the study of NICCE network, 47.6% (14 020/29 445) <italic>C. sinensis</italic>
genes can be tested by this citrus microarray<sup><xref ref-type="bibr" rid="bib17">17</xref>
</sup>
. Similar results were also found in maize that only 56.5% genes could be detected by maize microarrays<sup><xref ref-type="bibr" rid="bib41">41</xref>
</sup>
.</p>
<p>The coexpression networks were constructed using the RMTGeneNet<sup><xref ref-type="bibr" rid="bib28">28</xref>
</sup>
, which first calculated pair-wise Pearson Correlation Coefficients (PCC) for all genes and then identified a threshold for cutting PCC values using Random Matrix Theory. The PCC thresholds for these coexpression networks are shown in <xref rid="tbl2" ref-type="table">Table 2</xref>
, ranging from 0.882 for “all data” to 0.968 for “HLB”. At these relatively stringent thresholds, only the top 0.24% to 1.06% of all possible edges was retained. The nodes of these networks range from 1137 to 2263, accounting for 9.47%–18.85% measurable genes of the citrus microarray.</p>
</sec>
<sec><title>Network topology</title>
<p><xref ref-type="fig" rid="fig2">Figure 2</xref>
displays the coexpression network of “all data” using Cytoscape 2.8.3<sup><xref ref-type="bibr" rid="bib29">29</xref>
</sup>
. Although these networks have different numbers of nodes and edges, they have similar topological characteristics (<xref rid="tbl2" ref-type="table">Table 2</xref>
). All these networks are composed of a major component and other small components. All nodes within one component are directly or indirectly connected. Nodes in the major components account for 62.56% to 82.28% of that in corresponding networks. The average path length of these networks ranges from 6.66 to 10.75, implying the small-world properties. The nodes degree of these networks fits to a power law distribution with the degree exponent (<italic>r</italic>
) ranging from 1.13 to 2.10, indicating that these networks are scale free. These networks demonstrate modular and hierarchical characteristics with the average clustering coefficient ranging from 0.20 to 0.40, which is more than 36 times higher than that of same size random networks (data not shown). <xref ref-type="supplementary-material" rid="xob3">Table S2</xref>
provides a list of all edges in these networks. Taken together, these seven networks contain 37 633 edges among 6256 nodes (genes, <xref ref-type="supplementary-material" rid="xob4">Table S3</xref>
), which account for 52.11% measurable genes of the microarray or 25.50% total genes in <italic>C. x clementina</italic>
v1.0 genome. <xref rid="tbl3" ref-type="table">Table 3</xref>
shows the intersections between nodes and edges of these networks. Generally, the intersections among them are relatively low. In total, 3304 nodes (52.81%) and 34 860 edges (92.63%) were found only in one network. The intersection among “leaves”, “all data”, and “citrus canker” network is relatively higher than that between other networks.</p>
</sec>
<sec><title>Network clustering and functional enrichment</title>
<p>MCL algorithm was used to identify sets of nodes (i.e. coexpression module) that are more densely connected with each other than with the remaining nodes of the network<sup><xref ref-type="bibr" rid="bib4">4</xref>
</sup>
. The inflation parameter (I), the most important parameter of MCL, was chosen according to area fraction, mass fraction, and efficiency. In the present study, more than 80% of the entire edge masses could be captured only using less than 3% of the network area (<xref rid="tbl4" ref-type="table">Table 4</xref>
). A total of 2338 modules were detected in these seven networks (<xref ref-type="supplementary-material" rid="xob5">Table S4</xref>
), with 525 of them containing five or more nodes. The size of biggest modules in these networks ranges from 47 to 200. Functional enrichment analyses of these 525 modules were performed using terms from the GO biological process and KEGG pathway (<xref ref-type="supplementary-material" rid="xob6">Tables S5</xref>
and <xref ref-type="supplementary-material" rid="xob7">S6</xref>
). Only terms enriched within a module with a Fisher’s <italic>p</italic>
-value of 0.05 or less were considered. Only 343 modules in these networks had some degree of GO enrichment. Some GO terms were commonly enriched in these networks, such as gene expression (GO: 0010467), translation (GO: 0006412), and photosynthesis (GO: 0015979). However, the gene numbers associated with these common GO terms varied among networks. For example, 33 and 28 genes were associated with photosynthesis (GO: 0015979) in the “all data” and “leaves” network, respectively. However, zero and five genes were related to photosynthesis in the “albedo” and “flesh” network, respectively. A total of 132 GO terms (28.5%) were enriched exclusively in one network, such as polysaccharide catabolic process (GO: 0000272) and trehalose metabolic process (GO: 0005991) in “citrus canker” network.</p>
<p>Because a small portion of all nodes (25.15%, 1574/6256) was annotated with KEGG orthology identifiers in the <italic>C. x clementina</italic>
annotation file, their homologs in Arabidopsis were used for KEGG enrichment. A total of 60 modules were detected with significantly enriched KEGG pathways, and 36 KEGG pathways were enriched in at least one module. Some pathways were commonly enriched in these networks, such as Ribosome (ath03010), and photosynthesis (ath00195). A clear correspondence was observed between GO and KEGG enrichment analyses.</p>
</sec>
<sec><title>Predominant function of selected modules</title>
<p>Four modules were presented below to illustrate the correspondence of these modules with defined biological functions and methods that can be used to explore functional modules from these gene coexpression networks.</p>
</sec>
<sec><title>(1) Citrus lateral organ boundaries 1 in “citrus canker“ network</title>
<p>The guide-gene approach is commonly used to explore functional modules from gene coexpression networks. A lateral organ boundaries 1 (<italic>CsLOB1</italic>
) gene has recently been identified as a citrus canker disease susceptibility gene in sweet orange<sup><xref ref-type="bibr" rid="bib42">42</xref>
</sup>
. The precise function of <italic>CsLOB1</italic>
is still not clear. Using its homolog in <italic>C. x clementina</italic>
(<italic>Ciclev10033956m</italic>
) as a guide, 25 coexpressed genes were identified in module 1 of the “citrus canker” network (<xref ref-type="fig" rid="fig3">Figure 3</xref>
). Six of them were involved in cell wall metabolism: <italic>Ciclev10005888m</italic>
(plant pectin methylesterase inhibitor superfamily protein), <italic>Ciclev10016123m</italic>
(xyloglucan endotransglucosylase/hydrolase 5), <italic>Ciclev10021623m</italic>
(expansin B2), <italic>Ciclev10007670m</italic>
(proline-rich extensin-like receptor kinase), <italic>Ciclev10014994m</italic>
(glycosyl hydrolase), and <italic>Ciclev10019941m</italic>
(pectin lyase-like superfamily protein). Similar results were reported in the NICCE networks<sup><xref ref-type="bibr" rid="bib17">17</xref>
</sup>
. Interestingly, three minichromosome maintenance family genes (<italic>Ciclev10007588m</italic>
, <italic>Ciclev10027769m</italic>
, and <italic>Ciclev10019324m</italic>
) were coexpressed with <italic>Ciclev10033956m</italic>
, implying the functions of LOB1 in DNA replication. Another candidate target of TAL effectors, <italic>CsSWEET1</italic>
(<italic>Ciclev10002276m</italic>
)<sup><xref ref-type="bibr" rid="bib42">42</xref>
</sup>
, was also included in module 1 of the “citrus canker” network. It encodes a sugar transporter for pathogen nutrition and is linked to <italic>Ciclev10033956m</italic>
through three nodes (the shortest path).</p>
</sec>
<sec><title>(2) Module 25 in “citrus canker network” (canker-module 25): plant hormone signal transduction</title>
<p>Canker-module 25 was selected based on functional enrichment analyses. It has 10 nodes, 14 edges and a density of 0.311 (<xref ref-type="fig" rid="fig4">Figure 4</xref>
). The highest ranked (lowest <italic>p</italic>
value) GO term of this module was response to oxidative stress (GO: 0006979, <italic>p</italic>
= 0.05). The highest ranked KEGG pathway of this module was plant hormone signal transduction (ath04075, <italic>p</italic>
= 0.00024). Increased ethylene production was reported in citrus leaves inoculated with <italic>Xanthomonas campestris</italic>
pv. <italic>citri</italic>
(Hasse) Dye (Xc), a strain of bacteria that causes citrus canker<sup><xref ref-type="bibr" rid="bib43">43</xref>
</sup>
. However, the ethylene signal transduction pathway is not clear in citrus. Three nodes in this module, <italic>Ciclev10019132m</italic>
(<italic>ERS1</italic>
, ethylene response sensor 1), <italic>Ciclev10021170m</italic>
(MAP kinase kinase) and <italic>Ciclev10005820m</italic>
(<italic>ERF1</italic>
, ethylene response factor 1), may be involved in the ethylene signal transduction. The hub gene of this module is <italic>Ciclev10019132m</italic>
(<italic>ERS1</italic>
). In Arabidopsis, ethylene signal is first perceived by endoplasmic reticulum localized receptor (including ERS1) and then transduced to ERF and downstream targets through MAPK cascades<sup><xref ref-type="bibr" rid="bib44">44</xref>
,<xref ref-type="bibr" rid="bib45">45</xref>
</sup>
. A jasmonic acid-amido synthetase gene (<italic>Ciclev10019459m</italic>
) and a protein phosphatase 2C gene (<italic>Ciclev10004981m</italic>
) were also included in this module, implying the cross-talk among ethylene, JA and ABA signaling pathways. Other genes may also be involved in plant hormone signal transduction, such as <italic>Ciclev10024032m</italic>
(cysteine-rich receptor-like protein kinase) and <italic>Ciclev10001726m</italic>
(peroxidase gene). Therefore, canker-module 25 is likely to carry on the functions of plant hormone signal transduction. Unannotated genes in this module would be hypothesized to be related to plant hormone signal transduction.</p>
</sec>
<sec><title>(3) Module 19 in “flesh” network (flesh-module 19): fruit ripening</title>
<p>Flesh-module 19 was also selected based on functional enrichment analyses. It has 11 nodes, 11 edges, and a density of 0.2 (<xref ref-type="fig" rid="fig5">Figure 5</xref>
). The highest ranked KEGG pathway of this module was the citrate cycle (TCA cycle) (ath00020, <italic>p</italic>
= 0.00026). The citrate cycle is the major pathway for the synthesis of citric acid, the most abundant organic acid in citrus<sup><xref ref-type="bibr" rid="bib46">46</xref>
</sup>
. At least three nodes of this module were related to the citrate cycle: <italic>Ciclev10008189m</italic>
(dihydrolipoamide succinyltransferase gene), <italic>Ciclev10025308m</italic>
(dihydrolipoamide acetyltransferase gene), and <italic>Ciclev10013692m</italic>
(acyl-activating enzyme 5 gene). Two nodes were involved in the biosynthesis of the polyphenol compounds: <italic>Ciclev10019346m</italic>
(UDP-glycosyltransferase gene) and <italic>Ciclev10011175m</italic>
(phenylalanine ammonia lyase gene). One node, <italic>Ciclev10028195m</italic>
(glucose-1-phosphate adenylyltransferase gene) was involved in glycogen biosynthesis. All these nodes were linked by <italic>Ciclev10006509m</italic>
, which encodes a subunit of a RUB (Related to Ubiquitin)-activating enzyme. The proteins encoded by these genes may be subject to similar post-translational modifications.</p>
</sec>
<sec><title>(4) Module 6 in “HLB” network (HLB-module 6): programmed cell death</title>
<p>HLB-module 6 has 28 nodes, 75 edges, and a density of 0.198 (<xref ref-type="fig" rid="fig6">Figure 6</xref>
). This module was selected because 13 of these 28 genes (46.43%) were only included in “HLB” network. Seventeen genes were assigned to specific GO terms. The highest ranked GO term of this module was programmed cell death (PCD, GO: 0012501, <italic>p</italic>
= 0.005). PCD is widely observed in plants in response to pathogenic infection. At least eight genes in this module were related to PCD. Bcl-2-associated athanogene gene (<italic>Ciclev10018596m</italic>
) plays a critical role in PCD<sup><xref ref-type="bibr" rid="bib47">47</xref>
</sup>
. It can suppress PCD via its interaction with Hsc70 and Hsp40 (<italic>Ciclev10000372m</italic>
)<sup><xref ref-type="bibr" rid="bib48">48</xref>
</sup>
. However, the up-regulation of genes involved in the ubiquitin-proteasome system can activate PCD<sup><xref ref-type="bibr" rid="bib49">49</xref>
</sup>
. <italic>Ciclev10008240m</italic>
(polyubiquitin 10) and <italic>Ciclev10005221m</italic>
(RING finger E3 ubiquitin ligases) are parts of the ubiquitin-proteasome system. Other genes related to PCD include: <italic>Ciclev10005800m</italic>
(myosin heavy chain-related), <italic>Ciclev10032432m</italic>
(sphingoid base hydroxylase), <italic>Ciclev10021281m</italic>
(LAG1 longevity assurance homolog 3), and <italic>Ciclev10032631m</italic>
(Glutaredoxin family protein). Their functions in HLB still need to be determined.</p>
</sec>
<sec><title>Comparison with NICCE network</title>
<p>When this manuscript was being prepared, a citrus gene coexpression network (called “NICCE network” in this study) based on publicly available microarray data sets was reported<sup><xref ref-type="bibr" rid="bib17">17</xref>
</sup>
. There are several differences between the NICCE networks and networks in this study.</p>
<p>First, probe sets, rather than genes, were used to construct the NICCE networks. In the 30 217 nodes of the NICCE networks, 5960 (19.7%) nodes were not mapped to any citrus transcripts; 9336 (30.9%) nodes belonged to the “one probe set per transcript” group. 5775 transcripts (38.2%) were represented by the remaining 14 921 (49.4%) probe sets (<xref ref-type="supplementary-material" rid="xob8">Table S7</xref>
). Therefore, 5.9% of the edges of NICCE networks were between probe sets of the same transcript/gene. Probe sets representing the same transcript were expected to have similar expression levels and appear in the same cluster of one network. However, this is not the case in the NICCE networks. One example (<italic>Cs1g07330.1</italic>
) was shown in <xref ref-type="supplementary-material" rid="xob9">Table S8</xref>
.</p>
<p>Second, when constructing the NICCE networks, PCC values between probe sets were transformed into highest reciprocal ranks (HRR), and the top 100 HRR for a given probe set was considered. This leads to most PCC values between nodes of NICCE networks being very low. <italic>Cs5g33560</italic>
was given as an example in the website of NICCE (<ext-link ext-link-type="uri" xlink:href="http://citrus.adelaide.edu.au/nicce/home.aspx">http://citrus.adelaide.edu.au/nicce/home.aspx</ext-link>
). However, the PCC values between <italic>Cs5g33560</italic>
and its coexpressed genes in condition-independent network range from 0.68 to 0.39. More attention should be paid to assess gene pairs with low PCC values.</p>
<p>Third, only sweet orange microarrays were used in this study, and they were classified into six condition-dependent data sets: citrus canker, HLB, leaves, flavedo, albedo, and flesh. In the NICCE networks, 297 microarrays from different species of citrus (including mandarin, sweet orange, lemon, and pummelo) were used, and they were classified into four condition-dependent data sets: sweet orange, fruit, leaf, and stress<sup><xref ref-type="bibr" rid="bib17">17</xref>
</sup>
.</p>
<p>In order to compare our networks with the NICCE networks, <italic>C. x clementina</italic>
gene IDs from our networks were transformed to <italic>C. sinensis</italic>
gene IDs. <italic>C. sinensis</italic>
orthologs were not identified for 1504 <italic>C. x clementina</italic>
genes in our networks. Therefore, only 26 191 edges in our networks were used in the comparison with the NICCE networks, whose nodes were also transformed to <italic>C. sinensis</italic>
gene IDs. Only 3868 edges were found in common between the two networks. About 85% of edges in our networks were not included in the NICCE networks. This may be due to different classification methods for microarray data sets. Most edges (72.84%) in our networks were exclusively found in condition-dependent networks.</p>
</sec>
<sec><title>Validation of coexpression networks using RNA-seq data</title>
<p>To confirm the coexpression networks in this study, 500 edges among 371 genes were randomly selected from the “all data” network. The expression of these genes was examined using another gene expression data set (<xref ref-type="supplementary-material" rid="xob10">Table S9</xref>
) in CAP<sup><xref ref-type="bibr" rid="bib34">34</xref>
</sup>
. The correlation coefficients (<italic>r</italic>
) between them were computed. The distribution of these correlation coefficients was highly skewed, as shown in <xref ref-type="fig" rid="fig7">Figure 7</xref>
. For 353 edges (70.6%), <italic>r</italic>
values were higher than the PCC thresholds that were used to construct the “all data” networks (0.882). <italic>R</italic>
values of 385 edges (77.0%) were higher than 0.8. These results suggest that the coexpression networks in this study are reliable.</p>
</sec>
</sec>
<sec sec-type="discussion"><title>Discussion</title>
<p>In this study, 230 citrus microarrays from a diverse collection of experiments were used to construct seven coexpression networks. The nodes of these networks range from 1137 to 2263, accounting for 9.47%–18.85% measurable genes of the citrus microarray. This is consistent with Ficklin’s work on rice<sup><xref ref-type="bibr" rid="bib20">20</xref>
</sup>
, which also employed the RMT method to select a threshold for rice coexpression network. 10% of the measurable genes on rice microarray were included in their network. The percentage is relatively low compared with other studies using empirical thresholds. For example, in the Arabidopsis coexpression network, the PCC cutoff value was set to 0.75 and 38% measurable genes were retained<sup><xref ref-type="bibr" rid="bib4">4</xref>
</sup>
. RMT method was taken from the field of particle physics and had been used to construct gene coexpression networks for <italic>Escherichia coli</italic>
, yeast, human, Arabidopsis, rice, and maize<sup><xref ref-type="bibr" rid="bib20">20</xref>
,<xref ref-type="bibr" rid="bib28">28</xref>
,<xref ref-type="bibr" rid="bib50">50</xref>
,<xref ref-type="bibr" rid="bib51">51</xref>
</sup>
. It has been demonstrated to be a reliable method for generating networks across a wide range of data sets<sup><xref ref-type="bibr" rid="bib50">50</xref>
</sup>
. It should be mentioned that after combining the seven coexpression networks, the nodes captured in our study reached 52.11% of the measurable genes of the microarray.</p>
<p>Both condition-independent and condition-dependent analyses were employed to ensure that coexpressed genes in special conditions were not lost. 77.77% nodes and 72.84% edges in our networks were exclusively found in condition-dependent networks. Function analysis of modules yielded similar results. 66.31% enriched GO terms were identified only in condition-dependent networks, such as programmed cell death in “citrus canker”, “HLB”, and “albedo” network. Condition-independent analysis was considered to be suitable for identifying globally coexpressed genes<sup><xref ref-type="bibr" rid="bib7">7</xref>
</sup>
, such as genes in photosynthesis, ribosome and DNA metabolism. In this study, we found that condition-independent analysis was not sufficient to identify all the genes in these pathways. For example, 159 ribosome genes could be measured in the citrus microarray (<xref ref-type="supplementary-material" rid="xob11">Table S10</xref>
). Thirty-seven ribosome genes were included in the condition-independent network (“all data”). This number is much smaller compared with 147 ribosome genes in “leaves” network. It has been demonstrated that gene coexpression analysis using too many microarray samples could result in the loss of information<sup><xref ref-type="bibr" rid="bib52">52</xref>
</sup>
. Therefore, condition-dependent analysis is necessary even for identifying globally coexpressed genes.</p>
<p>According to the present annotation of the <italic>C. x clementina</italic>
genome, 2485 and 4682 (39.72% and 74.84%) genes in these networks were not assigned to a specific GO and KEGG pathway term, respectively<sup><xref ref-type="bibr" rid="bib8">8</xref>
</sup>
. The function of these genes could be predicted based on well-annotated genes within the same module. For example, 28 genes were included in HLB-module 6. Eleven of them were not labeled with a specific GO term, and only four genes were assigned to a specific KEGG pathway. Based on the above analysis, HLB-module 6 is likely to carry on the functions of programmed cell death. Unannotated genes in this module could be hypothesized to be related to programmed cell death. In addition to gene function prediction, gene coexpression analysis is also helpful for hypothesis generation and testing<sup><xref ref-type="bibr" rid="bib7">7</xref>
</sup>
. For example, several genes encoding transcription factors were also included in HLB-module 6, such as ERF and KH domain-containing putative RNA-binding protein. It has been demonstrated in Arabidopsis that a KH domain-containing putative RNA-binding protein is critical for HSF and HSP regulation<sup><xref ref-type="bibr" rid="bib53">53</xref>
</sup>
. Therefore, it would be reasonable to hypothesize that those transcription factors can regulate the expression of other genes within the same module.</p>
</sec>
</body>
<back><ack><p>This work was supported by a grant from the Citrus Research and Development Foundation (CRDF-724).</p>
</ack>
<notes><p>The authors declare no conflict of interest.</p>
</notes>
<sec sec-type="supplementary-material"><supplementary-material content-type="local-data" id="xob1"><label>Figure S1.</label>
<media xlink:href="hortres201526-s1.jpg"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="xob2"><label>Table S1.</label>
<media xlink:href="hortres201526-s2.xls"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="xob3"><label>Table S2.</label>
<media xlink:href="hortres201526-s3.xls"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="xob4"><label>Table S3.</label>
<media xlink:href="hortres201526-s4.xls"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="xob5"><label>Table S4.</label>
<media xlink:href="hortres201526-s5.xls"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="xob6"><label>Table S5.</label>
<media xlink:href="hortres201526-s6.xls"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="xob7"><label>Table S6.</label>
<media xlink:href="hortres201526-s7.xls"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="xob8"><label>Table S7.</label>
<media xlink:href="hortres201526-s8.xls"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="xob9"><label>Table S8.</label>
<media xlink:href="hortres201526-s9.xls"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="xob10"><label>Table S9.</label>
<media xlink:href="hortres201526-s10.xls"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="xob11"><label>Table S10.</label>
<media xlink:href="hortres201526-s11.xls"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
<ref-list><ref id="bib1"><mixed-citation publication-type="journal">Wu LF, Hughes TR, Davierwala AP, Robinson MD, Stoughton R, Altschuler SJ. <article-title>Large-scale prediction of <italic>Saccharomyces cerevisiae</italic>
gene function using overlapping transcriptional clusters.</article-title>
<source>Nat Genet</source>
<year>2002</year>
; <volume>31</volume>
: <fpage>255</fpage>
–265.<pub-id pub-id-type="pmid">12089522</pub-id>
</mixed-citation>
</ref>
<ref id="bib2"><mixed-citation publication-type="journal">Lee HK, Hsu AK, Sajdak J, Qin J, Pavlidis P. <article-title>Coexpression analysis of human genes across many microarray data sets.</article-title>
<source>Genome Res</source>
<year>2004</year>
; <volume>14</volume>
: <fpage>1085</fpage>
–1094.<pub-id pub-id-type="pmid">15173114</pub-id>
</mixed-citation>
</ref>
<ref id="bib3"><mixed-citation publication-type="journal">Fuller TF, Ghazalpour A, Aten JE, Drake TA, Lusis AJ, Horvath S. <article-title>Weighted gene coexpression network analysis strategies applied to mouse weight.</article-title>
<source>Mamm Genome</source>
<year>2007</year>
; <volume>18</volume>
: <fpage>463</fpage>
–472.<pub-id pub-id-type="pmid">17668265</pub-id>
</mixed-citation>
</ref>
<ref id="bib4"><mixed-citation publication-type="journal">Mao LY, Hemert JLv, Dash S, Dickerson JA. <article-title>Arabidopsis gene co-expression network and its functional modules.</article-title>
<source>BMC Bioinformatics</source>
<year>2009</year>
; <volume>10</volume>
: <fpage>346</fpage>
.<pub-id pub-id-type="pmid">19845953</pub-id>
</mixed-citation>
</ref>
<ref id="bib5"><mixed-citation publication-type="journal">Liang Y-H, Cai B, Chen F et al. <article-title>Construction and validation of a gene co-expression network in grapevine (<italic>Vitis vinifera</italic>
. L.).</article-title>
<source>Hort Res</source>
<year>2014</year>
; <volume>1</volume>
: <fpage>14040</fpage>
.</mixed-citation>
</ref>
<ref id="bib6"><mixed-citation publication-type="journal">Aoki K, Ogata Y, Shibata D. <article-title>Approaches for extracting practical information from gene co-expression networks in plant biology.</article-title>
<source>Plant Cell Physiol</source>
<year>2007</year>
; <volume>48</volume>
: <fpage>381</fpage>
–390.<pub-id pub-id-type="pmid">17251202</pub-id>
</mixed-citation>
</ref>
<ref id="bib7"><mixed-citation publication-type="journal">Usadel B, Obayashi T, Mutwil M et al. <article-title>Co-expression tools for plant biology: opportunities for hypothesis generation and caveats.</article-title>
<source>Plant Cell Environ</source>
<year>2009</year>
; <volume>32</volume>
: <fpage>1633</fpage>
–1651.<pub-id pub-id-type="pmid">19712066</pub-id>
</mixed-citation>
</ref>
<ref id="bib8"><mixed-citation publication-type="journal">Wu GA, Prochnik S, Jenkins J et al. <article-title>Sequencing of diverse mandarin, pummelo and orange genomes reveals complex history of admixture during citrus domestication.</article-title>
<source>Nat Biotechnol</source>
<year>2014</year>
; <volume>32</volume>
: <fpage>656</fpage>
–662.<pub-id pub-id-type="pmid">24908277</pub-id>
</mixed-citation>
</ref>
<ref id="bib9"><mixed-citation publication-type="journal">Aritua V, Achor D, Gmitter FG, Albrigo G, Wang N. <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>
: <fpage>e73742</fpage>
.<pub-id pub-id-type="pmid">24058486</pub-id>
</mixed-citation>
</ref>
<ref id="bib10"><mixed-citation publication-type="journal">Febres VJ, Khalaf A, Gmitter FG, Moore GA. <article-title>Evaluating gene expression responses of citrus to two types of defense inducers using a newly developed citrus agilent microarray.</article-title>
<source>Acta Hortic</source>
<year>2012</year>
; <volume>929</volume>
: <fpage>59</fpage>
–64.</mixed-citation>
</ref>
<ref id="bib11"><mixed-citation publication-type="journal">Fan J, Chen C, Yu Q et al. <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>
: <fpage>1396</fpage>
–1407.<pub-id pub-id-type="pmid">22809274</pub-id>
</mixed-citation>
</ref>
<ref id="bib12"><mixed-citation publication-type="journal">Khalaf AA, Gmitter FG Jr., Conesa A, Dopazo J, Moore GA. <article-title><italic>Fortunella margarita</italic>
transcriptional reprogramming triggered by <italic>Xanthomonas citri</italic>
subsp. <italic>citri</italic>
.</article-title>
<source>BMC Plant Biol</source>
<year>2011</year>
; <volume>11</volume>
: <fpage>159</fpage>
.<pub-id pub-id-type="pmid">22078099</pub-id>
</mixed-citation>
</ref>
<ref id="bib13"><mixed-citation publication-type="journal">Fan J, Chen C, Yu Q, Brlansky RH, Li Z-G, Gmitter FG Jr. <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>
: <fpage>235</fpage>
–245.</mixed-citation>
</ref>
<ref id="bib14"><mixed-citation publication-type="journal">Khalaf AA, Gmitter FG, Brlansky RH, Fan J, Moore GA. <article-title>Gene expression profiling of two citrus cultivars in response to huanglongbing (HLB) using the Agilent Citrus custom microarray Chip.</article-title>
<source>Phytopathology</source>
<year>2010</year>
; <volume>100</volume>
: <fpage>S61</fpage>
.</mixed-citation>
</ref>
<ref id="bib15"><mixed-citation publication-type="journal">Zheng Z, Zhao Y. <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="pmid">23324561</pub-id>
</mixed-citation>
</ref>
<ref id="bib16"><mixed-citation publication-type="journal">Martinelli F, Uratsu SL, Albrecht U et al. <article-title>Transcriptome profiling of citrus fruit response to Huanglongbing disease.</article-title>
<source>PloS One</source>
<year>2012</year>
; <volume>7</volume>
: <fpage>16</fpage>
.</mixed-citation>
</ref>
<ref id="bib17"><mixed-citation publication-type="journal">Wong DCJ, Sweetman C, Ford CM. <article-title>Annotation of gene function in citrus using gene expression information and co-expression networks.</article-title>
<source>BMC Plant Biol</source>
<year>2014</year>
; <volume>14</volume>
: <fpage>186</fpage>
.<pub-id pub-id-type="pmid">25023870</pub-id>
</mixed-citation>
</ref>
<ref id="bib18"><mixed-citation publication-type="journal">Atias O, Chor B, Chamovitz DA. <article-title>Large-scale analysis of Arabidopsis transcription reveals a basal co-regulation network.</article-title>
<source>BMC Syst Biol</source>
<year>2009</year>
; <volume>3</volume>
: <fpage>22</fpage>
.<pub-id pub-id-type="pmid">19228404</pub-id>
</mixed-citation>
</ref>
<ref id="bib19"><mixed-citation publication-type="journal">Childs KL, Davidson RM, Buell CR. <article-title>Gene coexpression network analysis as a source of functional annotation for rice genes.</article-title>
<source>PLoS One</source>
<year>2011</year>
; <volume>6</volume>
: <fpage>e22196</fpage>
.<pub-id pub-id-type="pmid">21799793</pub-id>
</mixed-citation>
</ref>
<ref id="bib20"><mixed-citation publication-type="journal">Ficklin SP, Luo F, Feltus FA. <article-title>The association of multiple interacting genes with specific phenotypes in rice using gene coexpression networks.</article-title>
<source>Plant Physiol</source>
<year>2010</year>
; <volume>154</volume>
: <fpage>13</fpage>
–24.<pub-id pub-id-type="pmid">20668062</pub-id>
</mixed-citation>
</ref>
<ref id="bib21"><mixed-citation publication-type="journal">Ding Y, Chang J, Guo J et al. <article-title>Prediction and functional analysis of the sweet orange protein-protein interaction network.</article-title>
<source>BMC Plant Biol</source>
<year>2014</year>
; <volume>14</volume>
: <fpage>213</fpage>
.<pub-id pub-id-type="pmid">25091279</pub-id>
</mixed-citation>
</ref>
<ref id="bib22"><mixed-citation publication-type="journal">Ibanez AM, Martinelli F, Reagan RL et al. <article-title>Transcriptome and metabolome analysis of citrus fruit to elucidate puffing disorder.</article-title>
<source>Plant Sci</source>
<year>2014</year>
; <volume>217</volume>
: <fpage>87</fpage>
–98.<pub-id pub-id-type="pmid">24467900</pub-id>
</mixed-citation>
</ref>
<ref id="bib23"><mixed-citation publication-type="journal">Martinelli F, Reagan RL, Uratsu SL et al. <article-title>Gene regulatory networks elucidating Huanglongbing disease mechanisms.</article-title>
<source>PLoS One</source>
<year>2013</year>
; <volume>8</volume>
: <fpage>12</fpage>
.</mixed-citation>
</ref>
<ref id="bib24"><mixed-citation publication-type="journal">Barrett T, Wilhite SE, Ledoux P et al. <article-title>NCBI GEO: archive for functional genomics data sets-update.</article-title>
<source>Nucleic Acids Res</source>
<year>2013</year>
; <volume>41</volume>
: <fpage>D991</fpage>
–D995.<pub-id pub-id-type="pmid">23193258</pub-id>
</mixed-citation>
</ref>
<ref id="bib25"><mixed-citation publication-type="journal">Gautier L, Cope L, Bolstad BM, Irizarry RA. <article-title>Affy – analysis of Affymetrix GeneChip data at the probe level.</article-title>
<source>Bioinformatics</source>
<year>2004</year>
; <volume>20</volume>
: <fpage>307</fpage>
–315.<pub-id pub-id-type="pmid">14960456</pub-id>
</mixed-citation>
</ref>
<ref id="bib26"><mixed-citation publication-type="journal">Kauffmann A, Gentleman R, Huber W. <article-title>ArrayQualityMetrics – a bioconductor package for quality assessment of microarray data.</article-title>
<source>Bioinformatics</source>
<year>2009</year>
; <volume>25</volume>
: <fpage>415</fpage>
–416.<pub-id pub-id-type="pmid">19106121</pub-id>
</mixed-citation>
</ref>
<ref id="bib27"><mixed-citation publication-type="journal">Liu H, Zeeberg BR, Qu G et al. <article-title>AffyProbeMiner: a web resource for computing or retrieving accurately redefined Affymetrix probe sets.</article-title>
<source>Bioinformatics</source>
<year>2007</year>
; <volume>23</volume>
: <fpage>2385</fpage>
–2390.<pub-id pub-id-type="pmid">17660211</pub-id>
</mixed-citation>
</ref>
<ref id="bib28"><mixed-citation publication-type="journal">Gibson SM, Ficklin SP, Isaacson S, Luo F, Feltus FA, Smith MC. <article-title>Massive-scale gene co-expression network construction and robustness testing using random matrix theory.</article-title>
<source>PLoS One</source>
<year>2013</year>
; <volume>8</volume>
:<fpage>e55871</fpage>
.<pub-id pub-id-type="pmid">23409071</pub-id>
</mixed-citation>
</ref>
<ref id="bib29"><mixed-citation publication-type="journal">Smoot ME, Ono K, Ruscheinski J, Wang P-L, Ideker T. <article-title>Cytoscape 2.8: new features for data integration and network visualization.</article-title>
<source>Bioinformatics</source>
<year>2011</year>
; <volume>27</volume>
: <fpage>431</fpage>
–432.<pub-id pub-id-type="pmid">21149340</pub-id>
</mixed-citation>
</ref>
<ref id="bib30"><mixed-citation publication-type="journal">Doncheva NT, Assenov Y, Domingues FS, Albrecht M. <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>
: <fpage>670</fpage>
–685.<pub-id pub-id-type="pmid">22422314</pub-id>
</mixed-citation>
</ref>
<ref id="bib31"><mixed-citation publication-type="journal">Enright AJ, Van Dongen S, Ouzounis CA. <article-title>An efficient algorithm for large-scale detection of protein families.</article-title>
<source>Nucleic Acids Res</source>
<year>2002</year>
; <volume>30</volume>
: <fpage>1575</fpage>
–1584.<pub-id pub-id-type="pmid">11917018</pub-id>
</mixed-citation>
</ref>
<ref id="bib32"><mixed-citation publication-type="journal">Goodstein DM, Shu S, Howson R et al. <article-title>Phytozome: a comparative platform for green plant genomics.</article-title>
<source>Nucleic Acids Res</source>
<year>2012</year>
; <volume>40</volume>
: <fpage>D1178</fpage>
–D1186.<pub-id pub-id-type="pmid">22110026</pub-id>
</mixed-citation>
</ref>
<ref id="bib33"><mixed-citation publication-type="journal">Gentleman RC, Carey VJ, Bates DM et al. <article-title>Bioconductor: open software development for computational biology and bioinformatics.</article-title>
<source>Genome Biol</source>
<year>2004</year>
; <volume>5</volume>
: <fpage>R80</fpage>
.<pub-id pub-id-type="pmid">15461798</pub-id>
</mixed-citation>
</ref>
<ref id="bib34"><mixed-citation publication-type="journal">Wang J, Chen D, Lei Y et al. <article-title><italic>Citrus sinensis</italic>
annotation project (CAP): a comprehensive database for sweet orange genome.</article-title>
<source>PLoS One</source>
<year>2014</year>
; <volume>9</volume>
: <fpage>e87723</fpage>
.<pub-id pub-id-type="pmid">24489955</pub-id>
</mixed-citation>
</ref>
<ref id="bib35"><mixed-citation publication-type="journal">Lee T-H, Tang H, Wang X, Paterson AH. <article-title>PGDD: a database of gene and genome duplication in plants.</article-title>
<source>Nucleic Acids Res</source>
<year>2013</year>
; <volume>41</volume>
: <fpage>D1152</fpage>
–D1158.<pub-id pub-id-type="pmid">23180799</pub-id>
</mixed-citation>
</ref>
<ref id="bib36"><mixed-citation publication-type="journal">Du D, Cheng T, Pan H, Yang W, Wang J, Zhang Q. <article-title>Genome-wide identification, molecular evolution and expression analyses of the phospholipase D gene family in three Rosaceae species.</article-title>
<source>Sci Hortic (Amsterdam)</source>
<year>2013</year>
; <volume>153</volume>
: <fpage>13</fpage>
–21.</mixed-citation>
</ref>
<ref id="bib37"><mixed-citation publication-type="journal">McGinnis S, Madden TL. <article-title>BLAST: at the core of a powerful and diverse set of sequence analysis tools.</article-title>
<source>Nucleic Acids Res</source>
<year>2004</year>
; <volume>32</volume>
: <fpage>W20</fpage>
–W25.<pub-id pub-id-type="pmid">15215342</pub-id>
</mixed-citation>
</ref>
<ref id="bib38"><mixed-citation publication-type="journal">Wang Y, Tang H, DeBarry JD et al. <article-title>MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity.</article-title>
<source>Nucleic Acids Res</source>
<year>2012</year>
; <volume>40</volume>
: <fpage>e49</fpage>
.<pub-id pub-id-type="pmid">22217600</pub-id>
</mixed-citation>
</ref>
<ref id="bib39"><mixed-citation publication-type="journal">Wang X, Shi X, Li Z et al. <article-title>Statistical inference of chromosomal homology based on gene colinearity and applications to Arabidopsis and rice.</article-title>
<source>BMC Bioinformatics</source>
<year>2006</year>
; <volume>7</volume>
: <fpage>447</fpage>
.<pub-id pub-id-type="pmid">17038171</pub-id>
</mixed-citation>
</ref>
<ref id="bib40"><mixed-citation publication-type="journal">Stover E, Stange RR Jr., McCollum TG, Jaynes J, Irey M, Mirkov E. <article-title>Screening antimicrobial peptides in vitro for use in developing transgenic citrus resistant to Huanglongbing and citrus canker.</article-title>
<source>J Am Soc Hortic Sci</source>
<year>2013</year>
; <volume>138</volume>
: <fpage>142</fpage>
–148.</mixed-citation>
</ref>
<ref id="bib41"><mixed-citation publication-type="journal">Sekhon RS, Briskine R, Hirsch CN et al. <article-title>Maize gene atlas developed by RNA sequencing and comparative evaluation of transcriptomes based on RNA sequencing and microarrays.</article-title>
<source>PLoS One</source>
<year>2013</year>
; <volume>8</volume>
: <fpage>11</fpage>
.</mixed-citation>
</ref>
<ref id="bib42"><mixed-citation publication-type="journal">Hu Y, Zhang JL, Jia HG et al. <article-title>Lateral organ boundaries 1 is a disease susceptibility gene for citrus bacterial canker disease.</article-title>
<source>Proc Natl Acad Sci U S A</source>
<year>2014</year>
; <volume>111</volume>
: <fpage>E521</fpage>
–E529.<pub-id pub-id-type="pmid">24474801</pub-id>
</mixed-citation>
</ref>
<ref id="bib43"><mixed-citation publication-type="journal">Dutta S, Biggs RH. <article-title>Regulation of ethylene biosynthesis in citrus leaves infected with <italic>Xanthomonas campestris</italic>
pv. <italic>citri</italic>
.</article-title>
<source>Physiol Plantarum</source>
<year>1991</year>
; <volume>82</volume>
: <fpage>225</fpage>
–230.</mixed-citation>
</ref>
<ref id="bib44"><mixed-citation publication-type="journal">Ouaked F, Rozhon W, Lecourieux D, Hirt H. <article-title>A MAPK pathway mediates ethylene signaling in plants.</article-title>
<source>EMBO J</source>
<year>2003</year>
; <volume>22</volume>
: <fpage>1282</fpage>
–1288.<pub-id pub-id-type="pmid">12628921</pub-id>
</mixed-citation>
</ref>
<ref id="bib45"><mixed-citation publication-type="journal">Hall AE, Findell JL, Schaller GE, Sisler EC, Bleecker AB. <article-title>Ethylene perception by the ERS1 protein in Arabidopsis.</article-title>
<source>Plant Physiol</source>
<year>2000</year>
; <volume>123</volume>
: <fpage>1449</fpage>
–1458.<pub-id pub-id-type="pmid">10938361</pub-id>
</mixed-citation>
</ref>
<ref id="bib46"><mixed-citation publication-type="journal">Yu KQ, Xu Q, Da XL, Guo F, Ding YD, Deng XX. <article-title>Transcriptome changes during fruit development and ripening of sweet orange (<italic>Citrus sinensis</italic>
).</article-title>
<source>BMC Genomics</source>
<year>2012</year>
; <volume>13</volume>
: <fpage>10</fpage>
.<pub-id pub-id-type="pmid">22230690</pub-id>
</mixed-citation>
</ref>
<ref id="bib47"><mixed-citation publication-type="journal">Fang S, Li L, Cui B, Men S, Shen Y, Yang X. <article-title>Structural insight into plant programmed cell death mediated by BAG proteins in <italic>Arabidopsis thaliana</italic>
.</article-title>
<source>Acta Crystallogr D</source>
<year>2013</year>
; <volume>69</volume>
: <fpage>934</fpage>
–945.<pub-id pub-id-type="pmid">23695238</pub-id>
</mixed-citation>
</ref>
<ref id="bib48"><mixed-citation publication-type="journal">Luders J, Demand J, Papp O, Hohfeld J. <article-title>Distinct isoforms of the cofactor BAG-1 differentially affect Hsc70 chaperone function.</article-title>
<source>J Biol Chem</source>
<year>2000</year>
; <volume>275</volume>
: <fpage>14817</fpage>
–14823.<pub-id pub-id-type="pmid">10809723</pub-id>
</mixed-citation>
</ref>
<ref id="bib49"><mixed-citation publication-type="journal">Serrano I, Gu Y, Qi D, Dubiella U, Innes RW. <article-title>The Arabidopsis EDR1 protein kinase negatively regulates the ATL1 E3 ubiquitin ligase to suppress cell death.</article-title>
<source>Plant Cell</source>
<year>2014</year>
; <volume>26</volume>
: <fpage>4532</fpage>
–4546.<pub-id pub-id-type="pmid">25398498</pub-id>
</mixed-citation>
</ref>
<ref id="bib50"><mixed-citation publication-type="journal">Luo F, Yang Y, Zhong J et al. <article-title>Constructing gene co-expression networks and predicting functions of unknown genes by random matrix theory.</article-title>
<source>BMC Bioinformatics</source>
<year>2007</year>
; <volume>8</volume>
: <fpage>299</fpage>
.<pub-id pub-id-type="pmid">17697349</pub-id>
</mixed-citation>
</ref>
<ref id="bib51"><mixed-citation publication-type="journal">Ficklin SP, Feltus FA. <article-title>Gene coexpression network alignment and conservation of gene modules between two grass species: maize and rice.</article-title>
<source>Plant Physiol</source>
<year>2011</year>
; <volume>156</volume>
: <fpage>1244</fpage>
–1256.<pub-id pub-id-type="pmid">21606319</pub-id>
</mixed-citation>
</ref>
<ref id="bib52"><mixed-citation publication-type="journal">Cosgrove EJ, Gardner TS, Kolaczyk ED. <article-title>On the choice and number of microarrays for transcriptional regulatory network inference.</article-title>
<source>BMC Bioinformatics</source>
<year>2010</year>
; <volume>11</volume>
: <fpage>454</fpage>
.<pub-id pub-id-type="pmid">20825684</pub-id>
</mixed-citation>
</ref>
<ref id="bib53"><mixed-citation publication-type="journal">Guan QM, Wen CL, Zeng HT, Zhu JH. <article-title>A KH domain-containing putative RNA-binding protein is critical for heat stress-responsive gene regulation and thermotolerance in Arabidopsis.</article-title>
<source>Mol Plant</source>
<year>2013</year>
; <volume>6</volume>
: <fpage>386</fpage>
–395.<pub-id pub-id-type="pmid">23087326</pub-id>
</mixed-citation>
</ref>
</ref-list>
</back>
<floats-group><fig id="fig1"><label>Figure 1</label>
<caption><p>Work flow used for networks construction and clustering in the present study.</p>
</caption>
<graphic xlink:href="hortres201526-f1"></graphic>
</fig>
<fig id="fig2"><label>Figure 2</label>
<caption><p>Layout of the citrus “all data” coexpression network. The most overrepresented GO terms were shown for the 12 largest color-coded modules.</p>
</caption>
<graphic xlink:href="hortres201526-f2"></graphic>
</fig>
<fig id="fig3"><label>Figure 3</label>
<caption><p>Graph showing coexpressed genes of the C. clementina homolog of citrus LOB1 (Ciclev10033956m) and SWEET1 (Ciclev10002276m) in canker-module 1.</p>
</caption>
<graphic xlink:href="hortres201526-f3"></graphic>
</fig>
<fig id="fig4"><label>Figure 4</label>
<caption><p>Genes and edges in canker-module 25.</p>
</caption>
<graphic xlink:href="hortres201526-f4"></graphic>
</fig>
<fig id="fig5"><label>Figure 5</label>
<caption><p>Genes and edges in flesh-module 19.</p>
</caption>
<graphic xlink:href="hortres201526-f5"></graphic>
</fig>
<fig id="fig6"><label>Figure 6</label>
<caption><p>Genes and edges in HLB-module 6.</p>
</caption>
<graphic xlink:href="hortres201526-f6"></graphic>
</fig>
<fig id="fig7"><label>Figure 7</label>
<caption><p>Distribution of absolute value of correlation coefficients.</p>
</caption>
<graphic xlink:href="hortres201526-f7"></graphic>
</fig>
<table-wrap id="tbl1"><label>Table 1</label>
<caption><title>Composition of the 230 microarrays according to the experiment conditions and organs.</title>
</caption>
<table frame="hsides" rules="groups" border="1"><colgroup><col align="left"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
</colgroup>
<thead valign="bottom"><tr><th align="left" valign="top" charoff="50"> </th>
<th align="center" valign="top" charoff="50">Epicotyls</th>
<th align="center" valign="top" charoff="50">Root</th>
<th align="center" valign="top" charoff="50">Stem</th>
<th align="center" valign="top" charoff="50">Leaves</th>
<th align="center" valign="top" charoff="50">Flower</th>
<th align="center" valign="top" charoff="50">Peel</th>
<th align="center" valign="top" charoff="50">Flavedo</th>
<th align="center" valign="top" charoff="50">Albedo</th>
<th align="center" valign="top" charoff="50">Vascular core</th>
<th align="center" valign="top" charoff="50">Flesh</th>
<th align="center" valign="top" charoff="50">Seed</th>
<th align="center" valign="top" charoff="50">Total</th>
</tr>
</thead>
<tbody valign="top"><tr><td align="left" valign="top" charoff="50">Citrus canker</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">30</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">30</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Huanglongbing</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">3</td>
<td align="center" valign="top" charoff="50">3</td>
<td align="center" valign="top" charoff="50">12</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">12</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">12</td>
<td align="center" valign="top" charoff="50">12</td>
<td align="center" valign="top" charoff="50">4</td>
<td align="center" valign="top" charoff="50">58</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Other treatments</td>
<td align="center" valign="top" charoff="50">12</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">2</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">6</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">20</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Control</td>
<td align="center" valign="top" charoff="50">0</td>
<td align="center" valign="top" charoff="50">3</td>
<td align="center" valign="top" charoff="50">3</td>
<td align="center" valign="top" charoff="50">19</td>
<td align="center" valign="top" charoff="50">9</td>
<td align="center" valign="top" charoff="50">12</td>
<td align="center" valign="top" charoff="50">22</td>
<td align="center" valign="top" charoff="50">12</td>
<td align="center" valign="top" charoff="50">7</td>
<td align="center" valign="top" charoff="50">31</td>
<td align="center" valign="top" charoff="50">4</td>
<td align="center" valign="top" charoff="50">122</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Total</td>
<td align="center" valign="top" charoff="50">12</td>
<td align="center" valign="top" charoff="50">6</td>
<td align="center" valign="top" charoff="50">6</td>
<td align="center" valign="top" charoff="50">63</td>
<td align="center" valign="top" charoff="50">9</td>
<td align="center" valign="top" charoff="50">12</td>
<td align="center" valign="top" charoff="50">40</td>
<td align="center" valign="top" charoff="50">12</td>
<td align="center" valign="top" charoff="50">19</td>
<td align="center" valign="top" charoff="50">43</td>
<td align="center" valign="top" charoff="50">8</td>
<td align="center" valign="top" charoff="50">230</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap id="tbl2"><label>Table 2</label>
<caption><title>Topological characteristics of seven coexpression networks</title>
</caption>
<table frame="hsides" rules="groups" border="1"><colgroup><col align="left"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
</colgroup>
<thead valign="bottom"><tr><th align="left" valign="top" charoff="50"> </th>
<th align="center" valign="top" charoff="50">Arrays</th>
<th align="center" valign="top" charoff="50">RT</th>
<th align="center" valign="top" charoff="50">Nodes</th>
<th align="center" valign="top" charoff="50">Edges</th>
<th align="center" valign="top" charoff="50">AD</th>
<th align="center" valign="top" charoff="50">ND</th>
<th align="center" valign="top" charoff="50">NCC</th>
<th align="center" valign="top" charoff="50">NBC</th>
<th align="center" valign="top" charoff="50">EBC</th>
<th align="center" valign="top" charoff="50">APL</th>
<th align="center" valign="top" charoff="50">DE</th>
<th align="center" valign="top" charoff="50">CC</th>
</tr>
</thead>
<tbody valign="top"><tr><td align="left" valign="top" charoff="50">All data</td>
<td align="center" valign="top" charoff="50">230</td>
<td align="center" valign="top" charoff="50">0.882</td>
<td align="center" valign="top" charoff="50">1391</td>
<td align="center" valign="top" charoff="50">10 220</td>
<td align="center" valign="top" charoff="50">14.69</td>
<td align="center" valign="top" charoff="50">0.011</td>
<td align="center" valign="top" charoff="50">104</td>
<td align="center" valign="top" charoff="50">1046</td>
<td align="center" valign="top" charoff="50">9878</td>
<td align="center" valign="top" charoff="50">7.62</td>
<td align="center" valign="top" charoff="50">1.13</td>
<td align="center" valign="top" charoff="50">0.40</td>
</tr>
<tr><td align="left" valign="top" charoff="50">HLB</td>
<td align="center" valign="top" charoff="50">36</td>
<td align="center" valign="top" charoff="50">0.968</td>
<td align="center" valign="top" charoff="50">1400</td>
<td align="center" valign="top" charoff="50">5036</td>
<td align="center" valign="top" charoff="50">7.19</td>
<td align="center" valign="top" charoff="50">0.005</td>
<td align="center" valign="top" charoff="50">91</td>
<td align="center" valign="top" charoff="50">1040</td>
<td align="center" valign="top" charoff="50">4379</td>
<td align="center" valign="top" charoff="50">6.76</td>
<td align="center" valign="top" charoff="50">1.58</td>
<td align="center" valign="top" charoff="50">0.29</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Citrus canker</td>
<td align="center" valign="top" charoff="50">30</td>
<td align="center" valign="top" charoff="50">0.938</td>
<td align="center" valign="top" charoff="50">1841</td>
<td align="center" valign="top" charoff="50">5565</td>
<td align="center" valign="top" charoff="50">6.05</td>
<td align="center" valign="top" charoff="50">0.003</td>
<td align="center" valign="top" charoff="50">211</td>
<td align="center" valign="top" charoff="50">1265</td>
<td align="center" valign="top" charoff="50">5174</td>
<td align="center" valign="top" charoff="50">8.79</td>
<td align="center" valign="top" charoff="50">1.50</td>
<td align="center" valign="top" charoff="50">0.23</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Leaves</td>
<td align="center" valign="top" charoff="50">63</td>
<td align="center" valign="top" charoff="50">0.911</td>
<td align="center" valign="top" charoff="50">2263</td>
<td align="center" valign="top" charoff="50">11 535</td>
<td align="center" valign="top" charoff="50">10.19</td>
<td align="center" valign="top" charoff="50">0.005</td>
<td align="center" valign="top" charoff="50">158</td>
<td align="center" valign="top" charoff="50">1862</td>
<td align="center" valign="top" charoff="50">11 247</td>
<td align="center" valign="top" charoff="50">10.75</td>
<td align="center" valign="top" charoff="50">1.47</td>
<td align="center" valign="top" charoff="50">0.33</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Flavedo</td>
<td align="center" valign="top" charoff="50">40</td>
<td align="center" valign="top" charoff="50">0.964</td>
<td align="center" valign="top" charoff="50">1425</td>
<td align="center" valign="top" charoff="50">2467</td>
<td align="center" valign="top" charoff="50">3.46</td>
<td align="center" valign="top" charoff="50">0.002</td>
<td align="center" valign="top" charoff="50">123</td>
<td align="center" valign="top" charoff="50">1133</td>
<td align="center" valign="top" charoff="50">2289</td>
<td align="center" valign="top" charoff="50">9.63</td>
<td align="center" valign="top" charoff="50">2.10</td>
<td align="center" valign="top" charoff="50">0.20</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Albedo</td>
<td align="center" valign="top" charoff="50">31</td>
<td align="center" valign="top" charoff="50">0.947</td>
<td align="center" valign="top" charoff="50">1592</td>
<td align="center" valign="top" charoff="50">3725</td>
<td align="center" valign="top" charoff="50">4.68</td>
<td align="center" valign="top" charoff="50">0.003</td>
<td align="center" valign="top" charoff="50">162</td>
<td align="center" valign="top" charoff="50">996</td>
<td align="center" valign="top" charoff="50">2644</td>
<td align="center" valign="top" charoff="50">6.66</td>
<td align="center" valign="top" charoff="50">1.80</td>
<td align="center" valign="top" charoff="50">0.20</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Flesh</td>
<td align="center" valign="top" charoff="50">43</td>
<td align="center" valign="top" charoff="50">0.948</td>
<td align="center" valign="top" charoff="50">1137</td>
<td align="center" valign="top" charoff="50">2431</td>
<td align="center" valign="top" charoff="50">4.28</td>
<td align="center" valign="top" charoff="50">0.004</td>
<td align="center" valign="top" charoff="50">90</td>
<td align="center" valign="top" charoff="50">911</td>
<td align="center" valign="top" charoff="50">2290</td>
<td align="center" valign="top" charoff="50">7.59</td>
<td align="center" valign="top" charoff="50">1.86</td>
<td align="center" valign="top" charoff="50">0.20</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Total</td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50">6256</td>
<td align="center" valign="top" charoff="50">37 633</td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50"> </td>
</tr>
</tbody>
</table>
<table-wrap-foot><fn id="t2-tfoot1"><p>RT, RMT threshold; AD, average degree; ND, network density; NCC, number of connected components; NBC, nodes of biggest components; EBC, edges of biggest components; APL, average path length; DE, degree exponent; CC, clustering coefficient.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap id="tbl3"><label>Table 3</label>
<caption><title>Intersection between edges/nodes (upper/lower triangular) of networks</title>
</caption>
<table frame="hsides" rules="groups" border="1"><colgroup><col align="left"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
</colgroup>
<thead valign="bottom"><tr><th align="left" valign="top" charoff="50"> </th>
<th align="center" valign="top" charoff="50">All data</th>
<th align="center" valign="top" charoff="50">HLB</th>
<th align="center" valign="top" charoff="50">Citrus canker</th>
<th align="center" valign="top" charoff="50">Leaves</th>
<th align="center" valign="top" charoff="50">Flavedo</th>
<th align="center" valign="top" charoff="50">Albedo</th>
<th align="center" valign="top" charoff="50">Flesh</th>
</tr>
</thead>
<tbody valign="top"><tr><td align="left" valign="top" charoff="50">All data</td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50">206 (1.35%)</td>
<td align="center" valign="top" charoff="50">574 (3.64%)</td>
<td align="center" valign="top" charoff="50">1305 (6%)</td>
<td align="center" valign="top" charoff="50">158 (1.25%)</td>
<td align="center" valign="top" charoff="50">45 (0.32%)</td>
<td align="center" valign="top" charoff="50">28 (0.22%)</td>
</tr>
<tr><td align="left" valign="top" charoff="50">HLB</td>
<td align="center" valign="top" charoff="50">334 (11.97%)</td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50">26 (0.25%)</td>
<td align="center" valign="top" charoff="50">29 (0.18%)</td>
<td align="center" valign="top" charoff="50">14 (0.19%)</td>
<td align="center" valign="top" charoff="50">28 (0.32%)</td>
<td align="center" valign="top" charoff="50">15 (0.2%)</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Citrus canker</td>
<td align="center" valign="top" charoff="50">482 (14.91%)</td>
<td align="center" valign="top" charoff="50">308 (9.5%)</td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50">1271 (7.43%)</td>
<td align="center" valign="top" charoff="50">15 (0.19%)</td>
<td align="center" valign="top" charoff="50">10 (0.11%)</td>
<td align="center" valign="top" charoff="50">3 (0.04%)</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Leaves</td>
<td align="center" valign="top" charoff="50">665 (18.2%)</td>
<td align="center" valign="top" charoff="50">274 (7.48%)</td>
<td align="center" valign="top" charoff="50">925 (22.54%)</td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50">59 (0.42%)</td>
<td align="center" valign="top" charoff="50">21 (0.14%)</td>
<td align="center" valign="top" charoff="50">9 (0.06%)</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Flavedo</td>
<td align="center" valign="top" charoff="50">278 (9.87%)</td>
<td align="center" valign="top" charoff="50">257 (9.1%)</td>
<td align="center" valign="top" charoff="50">276 (8.45%)</td>
<td align="center" valign="top" charoff="50">361 (9.79%)</td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50">70 (1.13%)</td>
<td align="center" valign="top" charoff="50">32 (0.65%)</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Albedo</td>
<td align="center" valign="top" charoff="50">262 (8.78%)</td>
<td align="center" valign="top" charoff="50">267 (8.92%)</td>
<td align="center" valign="top" charoff="50">269 (7.84%)</td>
<td align="center" valign="top" charoff="50">317 (8.22%)</td>
<td align="center" valign="top" charoff="50">565 (18.73%)</td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50">40 (0.65%)</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Flesh</td>
<td align="center" valign="top" charoff="50">144 (5.7%)</td>
<td align="center" valign="top" charoff="50">238 (9.38%)</td>
<td align="center" valign="top" charoff="50">178 (5.98%)</td>
<td align="center" valign="top" charoff="50">217 (6.38%)</td>
<td align="center" valign="top" charoff="50">410 (16%)</td>
<td align="center" valign="top" charoff="50">382 (14%)</td>
<td align="center" valign="top" charoff="50"> </td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap id="tbl4"><label>Table 4</label>
<caption><title>Network clustering and functional enrichment of modules</title>
</caption>
<table frame="hsides" rules="groups" border="1"><colgroup><col align="left"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
<col align="center"></col>
</colgroup>
<tbody valign="top"><tr><td align="left" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50">Inflation</td>
<td align="center" valign="top" charoff="50">Efficiency</td>
<td align="center" valign="top" charoff="50">Mass fraction (%)</td>
<td align="center" valign="top" charoff="50">Area fraction (%)</td>
<td align="center" valign="top" charoff="50">Modules</td>
<td align="center" valign="top" charoff="50">M5</td>
<td align="center" valign="top" charoff="50">SBM</td>
<td align="center" valign="top" charoff="50">GO</td>
<td align="center" valign="top" charoff="50">KEGG</td>
</tr>
<tr><td align="left" valign="top" charoff="50">All data</td>
<td align="center" valign="top" charoff="50">2</td>
<td align="center" valign="top" charoff="50">0.52</td>
<td align="center" valign="top" charoff="50">83.63</td>
<td align="center" valign="top" charoff="50">2.85</td>
<td align="center" valign="top" charoff="50">285</td>
<td align="center" valign="top" charoff="50">55</td>
<td align="center" valign="top" charoff="50">200</td>
<td align="center" valign="top" charoff="50">39</td>
<td align="center" valign="top" charoff="50">7</td>
</tr>
<tr><td align="left" valign="top" charoff="50">HLB</td>
<td align="center" valign="top" charoff="50">1.8</td>
<td align="center" valign="top" charoff="50">0.45</td>
<td align="center" valign="top" charoff="50">80.93</td>
<td align="center" valign="top" charoff="50">1.38</td>
<td align="center" valign="top" charoff="50">267</td>
<td align="center" valign="top" charoff="50">72</td>
<td align="center" valign="top" charoff="50">90</td>
<td align="center" valign="top" charoff="50">47</td>
<td align="center" valign="top" charoff="50">9</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Citrus canker</td>
<td align="center" valign="top" charoff="50">1.8</td>
<td align="center" valign="top" charoff="50">0.56</td>
<td align="center" valign="top" charoff="50">85.90</td>
<td align="center" valign="top" charoff="50">0.81</td>
<td align="center" valign="top" charoff="50">447</td>
<td align="center" valign="top" charoff="50">88</td>
<td align="center" valign="top" charoff="50">103</td>
<td align="center" valign="top" charoff="50">58</td>
<td align="center" valign="top" charoff="50">9</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Leaves</td>
<td align="center" valign="top" charoff="50">1.8</td>
<td align="center" valign="top" charoff="50">0.46</td>
<td align="center" valign="top" charoff="50">82.00</td>
<td align="center" valign="top" charoff="50">1.40</td>
<td align="center" valign="top" charoff="50">432</td>
<td align="center" valign="top" charoff="50">87</td>
<td align="center" valign="top" charoff="50">144</td>
<td align="center" valign="top" charoff="50">57</td>
<td align="center" valign="top" charoff="50">17</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Flavedo</td>
<td align="center" valign="top" charoff="50">1.6</td>
<td align="center" valign="top" charoff="50">0.48</td>
<td align="center" valign="top" charoff="50">86.92</td>
<td align="center" valign="top" charoff="50">0.75</td>
<td align="center" valign="top" charoff="50">287</td>
<td align="center" valign="top" charoff="50">86</td>
<td align="center" valign="top" charoff="50">47</td>
<td align="center" valign="top" charoff="50">58</td>
<td align="center" valign="top" charoff="50">12</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Albedo</td>
<td align="center" valign="top" charoff="50">1.8</td>
<td align="center" valign="top" charoff="50">0.54</td>
<td align="center" valign="top" charoff="50">83.03</td>
<td align="center" valign="top" charoff="50">0.84</td>
<td align="center" valign="top" charoff="50">402</td>
<td align="center" valign="top" charoff="50">73</td>
<td align="center" valign="top" charoff="50">102</td>
<td align="center" valign="top" charoff="50">37</td>
<td align="center" valign="top" charoff="50">3</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Flesh</td>
<td align="center" valign="top" charoff="50">1.6</td>
<td align="center" valign="top" charoff="50">0.44</td>
<td align="center" valign="top" charoff="50">84.90</td>
<td align="center" valign="top" charoff="50">1.67</td>
<td align="center" valign="top" charoff="50">218</td>
<td align="center" valign="top" charoff="50">64</td>
<td align="center" valign="top" charoff="50">108</td>
<td align="center" valign="top" charoff="50">47</td>
<td align="center" valign="top" charoff="50">3</td>
</tr>
<tr><td align="left" valign="top" charoff="50">Total</td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50">2338</td>
<td align="center" valign="top" charoff="50">525</td>
<td align="center" valign="top" charoff="50"> </td>
<td align="center" valign="top" charoff="50">343</td>
<td align="center" valign="top" charoff="50">60</td>
</tr>
</tbody>
</table>
<table-wrap-foot><fn id="t4-tfoot1"><p>M5, number of modules containing more than five nodes; SBM, size of biggest module; GO, KEGG: number of modules with significantly enriched GO biological process terms, KEGG pathways.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</floats-group>
</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 000440 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000440 | 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:4595991 |texte= Construction of citrus gene coexpression networks from microarray data using random matrix theory }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:26504573" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \ | NlmPubMed2Wicri -a OrangerV1
![]() | This area was generated with Dilib version V0.6.25. | ![]() |