<title>Background</title> <p>The rapid increase of world population, the field degradation by soil salinization and erosion, and the likely fluctuations in climate caused by global warming will pose new and known challenges to agriculture during this century [##REF##10938790##1##]. Crop improvements required to cope with these challenges could be attained through agronomic advances, leading to a better use of fertilizers, protection agents or soil rescue, and exploitation of recent technologies for plant breeding. Despite the outstanding importance of genetics-based breeding applied to spontaneous mutations and conventional hybrids, molecular and genomic tools are expected to develop their great potential for crop improvement through functional genetics analysis, involving gene and function discovery and genome modification.</p> <p>Citrus, some of the most important fruit crops worldwide, are perennial trees requiring a juvenility period of several years and frequently are parthenocarpic and sexually self-incompatible [##UREF##0##2##,##UREF##1##3##], which considerably impairs traditional breeding. Genomic technology, including methods to rapidly identify and manipulate genes of agricultural interest, holds promise of improvements that may be difficult through traditional approaches. In recent years, <italic>Citrus </italic>has been the target of several genomic developments including large EST collections [##REF##15830128##4##, ####REF##17254327##5##, ##UREF##2##6##, ##UREF##3##7####3##7##], cDNA and oligonucleotide-based microarrays [##REF##15830128##4##,##UREF##4##8##,##UREF##5##9##], BAC libraries and BAC end sequencing (BES) (to be published). However, functional studies, i.e. genetic transformation and the capability to perform reverse genetic analyses, are also considerably impaired. In citrus, high throughput transgenic programs such as the generation of RNA interference knockouts, activation tagging through enhancer elements, gene-trap T-DNA insertions, or transposon tagging systems have not been developed yet. Since no efficient tagging or insertional procedures are available in these species, other gene disruption methods including strategies based on genome-wide mutagenesis such as TILLING and fast neutron mutagenesis have been initiated. These approaches are non-transgenic and may have particular interest for the industry where the debate on genetically modified organisms has restricted application of these technologies to crop improvement. Both approaches, however, are of limited usefulness as strategies for reverse genetics because of the lack of knowledge on <italic>Citrus </italic>genomic sequence and the large amount of space required for the establishment of mutant populations. ECOTILLING on natural citrus variants and microarray-based detection of deletions in fast neutron citrus mutants are apparently very straightforward approaches. In this work we explore the potential of this last idea using two fast neutron <italic>Citrus clementina </italic>hemizygous mutants from the IVIA collection and a 20K cDNA citrus microarray.</p> <p>Physical mutagenesis through fast neutron irradiation has been reported to cause variable genomic deletions ranging in size from few base pairs to 12 kb in <italic>Arabidopsis thaliana </italic>[##REF##8893551##10##,##REF##11532169##11##]. Several approaches have been used to characterize plant genomic deletions at the molecular level. These mostly include positional cloning [##REF##17496113##12##], a method applicable to any kind of genetic lesion that, however, needs highly saturated genetic maps; PCR-based reverse genetics techniques [##REF##11532169##11##,##UREF##6##13##], requiring a previous considerable knowledge of genomic sequence; and genomic subtraction procedures [##REF##2408039##14##, ####REF##1354004##15##, ##REF##9490740##16####9490740##16##], which do not need sequence information but are strongly dependent on the gene dosage. Since very little is known about <italic>Citrus </italic>genome sequence and the <italic>Citrus </italic>induced deletions are in hemizygous gene dosage, an array-based procedure as the one employed for identifying homozygous gene deletions in <italic>Arabidopsis </italic>[##REF##15486089##17##] seems more suitable for our purpose than those methods. Although the main application of microarrays is transcriptome profiling analysis, microarrays can also be used to study DNA variation. Oligonucleotide arrays are particularly suited for the detection of single nucleotide mismatches during hybridization, and hence for the discovery of novel DNA variants or the determination of known variants. The origin of this technique relies on a cytogenetic method described 25 years ago named "Comparative Genomic Hybridization" (CGH) that used differential DNA hybridization on chromosome spreads for visualization of deleted or amplified genomic regions in tumour tissues [##REF##1359641##18##]. Subsequently, different laboratories mostly working on cancer research independently applied microarray technology to genomic DNA hybridization procedure, a technique consequently named array-CGH [##REF##9842716##19##, ####REF##10706083##20##, ##REF##9771718##21##, ##REF##10471496##22##, ##REF##9408757##23####9408757##23##]. Array-CGH was successfully utilized to detect gene duplications in <italic>Arabidopsis </italic>and rice [##UREF##7##24##], and to validate aneuploidy analysis performed by quantitative fluorescent PCR in <italic>Arabidopsis </italic>[##REF##16995901##25##]. Therefore, this method has proven to be suitable to study chromosomal imbalances in plants.</p> <p>For the characterization of the deleted regions we also leaned on comparative genomics with other dicots since available physical citrus maps are not yet integrated with known genetic maps. Comparative genomics takes advantage of available information on gene content and order in genomic DNA from different species to infer phylogenetic relationships and formulate hypotheses on DNA evolutionary dynamics. Whole genomes are preferentially compared when available, but more often relatively short stretches of DNA or polymorphic markers are used.</p> <p>The main objective of this work was to identify deleted genes on a heterozygous genetic <italic>Citrus </italic>background, provided by fast neutron generated mutants, through array-Comparative Genomic Hybridization. In addition, we also explored the possibility of using comparative genomics with annotated dicot genomes assisted by BAC end sequencing for the generation of partial physical maps of the deleted <italic>Citrus </italic>regions.</p>
<title>Methods</title> <title>Plant material</title> <p>Approximately 6 years-old clementine trees (<italic>Citrus clementina </italic>Hort. Ex Tan. cv. clemenules) grown under standard agricultural practices at the Instituto Valenciano de Investigaciones Agrarias (IVIA) were used in this study. Commercial highly heterozygous clementine cultivars are considered "wild type" material, while the 39B3 and 39E7 genotypes that belong to the IVIA mutant collection were obtained through bud irradiation with fast neutrons (5–6 Gy) at the Instituto Tecnologico e Nuclear (Sacavem, Portugal) in the frame of a much wider breeding program. Both mutants are expected to carry DNA deletion lesions in hemizygous dosage and showed altered patterns of colour change of fruit peel.</p> <title>Array-CGH</title> <p>The protocol was adapted from several published array-Comparative Genomic Hybridization (array-CGH) methods pursuing mainly the measurement of copy-number changes in human genomic DNA [##REF##16030317##48##, ####REF##15961730##49##, ##UREF##12##50####12##50##], and the study of large-scale genetic variation of the symbiotic bacteria <italic>Sinorhizobium meliloti </italic>[##REF##16283928##51##]. Genomic DNA was isolated from leaves of wild type and mutant plants, using DNeasy plant mini kit (Qiagen). Four Cy3 or Cy5-labelled independent biological samples from each mutant plant were co-profiled on four 20K <italic>Citrus </italic>cDNA microarrays containing 21240 EST, using Cy5 or Cy3-labelled control genomic DNA, respectively. Label probes were prepared as follow: Cy3- or Cy5-dCTP fluorescent nucleotides (Amersham Biosciences) were incorporated directly in control and mutant genomic DNA (2 μg) using BioPrime Array CGH Genomic Labelling System (Invitrogene). Purified Cy5 and Cy3 labelled probes (about 50 μl each) were combined and mixed with 30 μg Cot-1 DNA (Invitrogene), 100 μg yeast tRNA (Invitrogene), and 346 μl TE buffer pH 7.4. Cot-1 DNA and yeast tRNA were used to block non-specific hybridization. Samples were laid on a microcon YM-30 filter (Millipore), and subsequently centrifuged until sample volume was reduced to approximately 48 μl. Finally, 10.2 μl 20× SSC and 1.8 μl 10% SDS were added to the probe mixture to reach a final volume of 60 μl containing 3.4× SSC and 0.3% SDS. For microarray hybridization, the probe mixture was denatured by heating at 97°C for 5 minutes, and immediately incubated at 37°C during 30 minutes to block repetitive DNA sequences. Hybridization mixture was applied to a 37°C pre-warmed hybrid-slip (Sigma), and a pre-warmed array slide was lowered onto the mix. Microarrays were hybridized in darkness at 65°C overnight (16–20 hours) using a glass array cassette following manufacturer's instructions (Ambion, cat. n° AM10040). To prevent evaporation of hybridization solution during incubation, 5 μl of 3× SCC were poured into the reservoir inside the cassette chamber. Following hybridization, microarray slides were placed in a rack and the cover slip removed by 10 minutes immersion in a washing chamber containing 2× SSC and 0.03% SDS at room temperature (RT). Microarray slides were passed through a series of washes on a shaking platform. Wash series were as follow: 2× SSC, 0.03% SDS for 5 min at 65°C, followed by 1× SSC for 5 min at RT, and 3 × 15 min washes in 0.2× SSC at RT. After first wash slides were transferred to new racks to minimize transference of SDS to the next washing solution. Microarray slides were dried by centrifugation for 5 min at 300 rpm by using an Eppendorf 5804-R tabletop centrifuge. Arrays were immediately scanned at 5 μm. Cy3 and Cy5 fluorescence intensity was collected by using a ScanArray Gx (Perkin Elmer). The resulting images were overlaid and spots identified by the ScanArray Express program (Perkin Elmer). Spot quality was first measured by the signal-to-background method with parameters lower limit (200) and multiplier (2), and subsequently confirmed by visual test. Data analysis was performed using the Limma package from the R statistical computing software [##UREF##13##52##, ####REF##16646809##53##, ##REF##15461798##54####15461798##54##]. A mutant/wild type signal lower than 0.7, with a P-value not higher than 0.1 (39E7) or 0.2 (39B3) were the cut-off values for positive EST identification. The experimental design of microarray experiments has been loaded into the ArrayExpress database [##UREF##14##55##] under accessions E-MEXP-1432 and E-MEXP-1433.</p> <title>Gene dosage measurements</title> <p>Quantitative real-time PCR was performed on a LightCycler 2.0 instrument (Roche), using the LightCycler FastStart DNA MasterPLUS SYBR Green I kit (Roche). Reaction composition and conditions followed manufacturer's instructions. Each individual PCR reaction contained 2 ng of genomic DNA from wild type or mutant, obtained with the DNeasy plant mini kit (Qiagen). Cycling protocol consisted of 10 min at 95°C for pre-incubation, then 40 cycles of 10 sec at 95°C for denaturation, 10 sec at 60°C for annealing and 10–25 sec at 72°C for extension. Fluorescent intensity data were acquired during the extension time. Specificity of the PCR reaction was assessed by the presence of a single peak in the dissociation curve after the amplification and through size estimation of the amplified product. For gene dosage measurements, we used the relative quantification-monocolor analysis from the LightCycler Software 4.0 package (Roche). This program compares the ratio of a target sequence to a reference DNA sequence in the mutant sample with the ratio of these sequences in a wild type sample. PCR and normalized calculations were repeated in at least three independent samples from each mutant and wild type, rendering an estimation of target gene dosage in the mutant genotype. Primers for the reference sequence were obtained from CX293764.</p> <title>Similarity searches</title> <p>DNA sequences of <italic>Citrus </italic>unigenes containing positive array-CGH ESTs were used in online TBLASTX searches against genomic databases from the annotated genomes of <italic>Arabidopsis thaliana </italic>[##UREF##15##56##], <italic>Populus trichocarpa </italic>[##UREF##16##57##] and <italic>Vitis vinifera </italic>[##UREF##17##58##] at an E-value cut-off of 10<sup>-5</sup>. For each gene, the best hit was placed on a chromosomal map while the second and third hits were only positioned in the map if they were located closer than 250 kb to any other hit. Two 700 kb regions from chromosomes 12 and 15 from the <italic>Populus </italic>genome including homologous genes to 39B3 array-CGH positive unigenes, were used as queries in a BLASTN local search on a <italic>Citrus </italic>BAC end sequence database. Only hits corresponding to those BAC ends showing an E-value lower than 10<sup>-5 </sup>in both chromosome searches were considered for the building of a local physical map of the 39B3 deletion.</p> <title>BAC isolation and analysis</title> <p>DNA from <italic>Citrus </italic>BACs was isolated with the Rapid Plasmid Miniprep System (Marligen Biosciences). Purified BACs were used as templates in PCR reactions in a total volume of 15 μl, including 0.2 mM dNTP, 2 mM MgCl<sub>2</sub>, 0.5 μM of each primer, 0.38 units of Netzyme DNA polymerase (Molecular Netline Bioproducts) and 0.1 ng of BAC DNA. After an initial denaturing step for 5 min at 95°C, amplification was performed for 35 cycles of 30 sec at 95°C, 30 sec at 60°C and 30 sec at 72°C, followed by 5 min incubation at 72°C. The PCR product was subjected to 1.5% agarose DNA electrophoresis.</p> <title>Chlorophyll measurements</title> <p>At least, three developing and mature leaves and fruit exocarp sectors from standard and 39B3 mutant lines of clementine were randomly collected per sample. Fruit exocarp tissues from a wild type clementine tree showing fruit colour delay due to altered environmental conditions were also sampled for chlorophyll analyses. Chlorophylls <italic>a </italic>and <italic>b </italic>were extracted with N,N-dimethylformamide for 72 h in the dark at 4°C and quantified through the absorbance at 647 and 664 nm following a reported procedure [##REF##16662407##59##]. Absorbance was measured using a Varian Cary 50 UV-visible spectrophotometer (Varian).</p> <title>Gene expression measurements</title> <p>Total RNA was extracted from fruit exocarp of wild type and 39B3 mutant using the RNeasy Plant Mini Kit (Qiagen). RNA concentration was determined by a fluorometric assay with the RiboGreen dye (Molecular Probes) following the manufacturer's instructions. About 5 μg of total RNA were reverse transcribed with the SuperScript III First-Strand Synthesis System for RT-PCR (Invitrogen) in a total volume of 20 μl. Single strand cDNA corresponding to <italic>ClpC</italic>-like and <italic>ClpD</italic>-like genes was amplified by quantitative real-time PCR on a LightCycler 2.0 instrument (Roche), using the LightCycler FastStart DNA MasterPLUS SYBR Green I kit (Roche). One μl of a 20 times diluted first-strand cDNA was used for each amplification reaction. Cycling protocol consisted of 10 min at 95°C for pre-incubation, then 40 cycles of 10 sec at 95°C for denaturation, 10 sec at 60°C for annealing and 15 sec at 72°C for extension. Melting curve analysis by applying increasing temperature from 65°C to 95°C (0.1°C/s) and gel electrophoresis of final product confirmed single amplicons. For expression measurements, we used the absolute quantification analysis from the LightCycler Software 4.0 package (Roche), and calculated expression levels relative to wild type values. Three independent biological samples were analyzed for wild type and mutant genotypes. Primers sequences are provided in Additional file ##SUPPL##1##2##.</p> <title>ClpC-like genomic sequence</title> <p><italic>ClpC</italic>-like genomic sequence from very few base pairs after the ATG until few base pairs before the stop codon was divided in four PCR fragments: Amplicon 3/4 (1820 bp) was amplified and sequenced with primers CLPC3 and CLPC4, amplicon 5/8 (2168 bp) was amplified with primers CLPC5 and CLPC8 and sequenced with primers CLPC5, CLPC8, CLPC10 and CLPC11, amplicon 7/2 (1446 bp) was amplified and sequenced with primers CLPC7 and CLPC2, and amplicon 1/6 (1158 bp) was amplified and sequenced with primers CLPC1 and CLPC6. Each amplicon was obtained by combining the product of 6–8 independent reactions. Primers sequences are provided in Additional file ##SUPPL##1##2##.</p>
<title>Results and discussion</title> <title>Procedure for the characterization of hemizygous deletions in <italic>Citrus</italic></title> <p>The proposed procedure to identify deleted genes is illustrated in Figure ##FIG##0##1## and its potential to structurally characterize hemizygous deletions is exemplified below with <italic>Citrus </italic>mutants as starting plant material. Its usefulness to describe genomic deletions in other species might be dependent upon genome complexity and ploidy. This method uses cDNA microarrays to hybridize genomic DNA extracted from the deletion mutants to render a list of underrepresented genes. The putative deleted genes are then validated through gene dosage evaluation by real-time PCR using gene specific primers. Deleted genes could subsequently contribute to the identification of the molecular mechanisms underlying the observed phenotypes by means of a candidate gene approach, validated by physiological analyses or genetic transformation [##UREF##8##26##]. In non-sequenced genomes or in plants with poorly developed physical maps, further characterization of deletions at the structural level requires TBLASTX similarity searches against databases containing the sequence annotation of known eudicot genomes, such as <italic>Arabidopsis thaliana</italic>, <italic>Populus trichocarpa </italic>and <italic>Vitis vinifera</italic>. These searches yield putative orthologous genes and syntenic genomic regions between these four species. Local physical maps of deletions are built allocating the deleted gene sequences and the syntenic genomic fragments from these other eudicots into a BES database of the species of interest. Lastly, specific PCR on the array of BACs confirms gene content and order on the lineal structure of the deletions. The results may also be used in comparative genomics analyses to study evolutionary dynamics and phylogenetics.</p> <title>Identification of deleted alleles in 39B3 and 39E7 fast neutron mutants of <italic>Citrus clementina</italic></title> <p>For this study, two mutants obtained by fast neutron mutagenesis of wild type <italic>Citrus clementina </italic>were selected from the IVIA mutant collection. These mutants, named 39B3 and 39E7, were expected to carry DNA deletion lesions in hemizygous dosage and showed a delay in natural colour break in fruit peel. The 39B3 mutant exhibited a delay in colour change from green to orange while 39E7 was better characterized by an abnormal final yellowish colour instead of the natural orange coloration. Putative deleted genes in the mutants were first identified through an approach based on genomic hybridization (array-CGH) that exploited a recently developed <italic>Citrus </italic>microarray containing 21240 cDNAs [##REF##15830128##4##,##REF##17254327##5##]. To this end, total genomic DNA from four independent samples of mutants 39B3 and 39E7 were Cy3 or Cy5-labelled and cohybridized with wild type DNA labelled with the complementary Cy5 or Cy3 probe on four independent microarray slides. Fluorescence intensity data were normalized and single ESTs showing a mutant/wild type signal ratio lower than 0.7 fold, with a P-value lower than 0.2 (39B3) or 0.1 (39E7), were selected as putative candidates.</p> <p>The number of ESTs fulfilling these criteria was 24 and 78 for mutants 39B3 and 39E7, respectively. One of the 39B3 positives [GenBank: <ext-link ext-link-type="gen" xlink:href="CX299090">CX299090</ext-link>], composed of three unrelated sequences was discarded for subsequent analysis due to its chimerical nature. In order to validate the array-CGH results, gene dosage of several putative candidates was determined through real-time PCR quantification of mutant/wild type signals for candidate ESTs as related to a reference undeleted gene [GenBank: <ext-link ext-link-type="gen" xlink:href="CX293764">CX293764</ext-link>]. The results showed that gene dosage for 39B3 candidates ranged from 0.50 to 0.60 when genomic DNA from the 39B3 genotype was tested, while ranged from 0.96 to 1.15 when the assayed DNA originated from the 39E7 genotype (Table ##TAB##0##1##). Similar results, corroborating the presence of putative deleted genes at half dosage, were also obtained for the 39E7 mutant. Therefore, the developed array-CGH procedure proved to be an appropriate tool to identify genes in hemizygous content in the self-incompatible clementine.</p> <title>Clustering of homologues of <italic>Citrus </italic>deleted genes in the poplar genome</title> <p>Microsynteny comparisons with homologous stretches from the sequenced genomes of <italic>Arabidopsis thaliana</italic>, <italic>Populus trichocarpa </italic>and <italic>Vitis vinifera </italic>[##REF##11130711##27##, ####REF##16973872##28##, ##REF##17721507##29####17721507##29##] were performed in order to elucidate hypothetical clustering of <italic>Citrus </italic>deleted genes in the genome. TBLASTX, which searches for translations of a crude genome similar to a translated query, was utilized with an E-value cut-off of 10<sup>-5</sup>. The homologous regions produced by the best TBLASTX hit of each of the <italic>Citrus </italic>candidate genes were located on the chromosome maps of <italic>Arabidopsis</italic>, poplar and grapevine. Homologues of <italic>Citrus </italic>genes were then grouped into clusters in each species when the distance between them was shorter than 250 kb. The second and third TBLASTX best alignments were similarly placed in the respective maps when they were included in an existing cluster. In this case, a binding line was drawn linking the second and third hits to the best hit of the same <italic>Citrus </italic>query. Thus, two chromosomal maps, one for each mutant, in the three species was obtained. Figure ##FIG##1##2## represents in detail chromosome mappings of the 39B3 mutation, which was subjected to further analyses. The results indicate that the <italic>Populus </italic>mapping exhibited rather lower complexity than the <italic>Arabidopsis </italic>and grapevine ones since it included fewer chromosomes and only 3 clusters although the number of 39B3 candidate genes represented in the map was identical (21) for the three genomes. Note that the number of represented hits in these mappings is higher than 21 due to the inclusion of second or third homologues. In <italic>Populus</italic>, most of the candidate genes mapped to two different genome regions of approximately 700 kb long in chromosomes 12 and 15, two duplicated chromosomes that probably originated during the recent genome duplication event that occurred in this species [##REF##16973872##28##]. These two clusters contained 17 and 15 hits respectively while the third one placed in chromosome 16 had only one hit. In contrast, the number of clusters in <italic>Arabidopsis </italic>and <italic>Vitis </italic>were 9 and 11, respectively, and none of them contained more than 11 hits. Furthermore, cluster number (and clustering density) of the homologues of 39E7 putative deleted genes was also lower (and higher) in <italic>Populus </italic>than in <italic>Arabidopsis </italic>or <italic>Vitis</italic>, although the differences were smaller: 26, 30 and 30 clusters were obtained for poplar, <italic>Arabidopsis </italic>and grapevine respectively (Figure ##FIG##2##3##).</p> <p>Overall, these observations suggest that the <italic>Populus </italic>genomic regions homologous to the <italic>Citrus </italic>deletions were less fragmented than their counterparts in <italic>Arabidopsis </italic>and <italic>Vitis</italic>, and consequently microsynteny on the considered segments was higher with the <italic>Populus </italic>genome. These results are striking since <italic>Citrus </italic>and <italic>Arabidopsis </italic>belong to Sapindales and Brassicales orders (inside the same clade eurosids II) while <italic>Populus </italic>is included in the eurosids I clade, and <italic>Vitis </italic>is part of Vitaceae, a family outside of rosids [##UREF##9##30##].</p> <title>Gene arrangement and partial physical map of the 39B3 deletion</title> <p>The closer microsynteny observed between the 39B3 deletion and the two duplicated homologous regions in poplar enabled prediction of gene order by direct inference from the <italic>Populus </italic>sequences. This assumption led to the gene arrangement depicted in Figure ##FIG##3##4##. Twenty genes out of twenty-one having high similarity with <italic>Populus </italic>homologues were directly located on the <italic>Citrus </italic>deletion fragment by combining the two clusters found on <italic>Populus </italic>chromosomes 12 and 15, which shared 12 hits. Inclusion of the 21<sup>st </sup>gene, a homologue of a <italic>Populus </italic>gene placed on chromosome 16, in the 39B3 deletion was based on its location on the right end of the <italic>Citrus </italic>BAC CCER1019D04 (named B12, see below), whose left end shared identity with another deleted gene [GenBank: <ext-link ext-link-type="gen" xlink:href="CX295702">CX295702</ext-link>]. The accession number and protein similarity of these 21 genes, numbered according to the ordered position of their homologues on the poplar genome (Figure ##FIG##3##4##), are depicted in Table ##TAB##1##2## that also shows coding strand sense of poplar homologues. The coding strand was coincident for the <italic>Populus </italic>paralogous genes present in chromosomes 12 and 15, except for genes similar to <italic>Citrus </italic><ext-link ext-link-type="gen" xlink:href="CX308429">CX308429</ext-link>, located in position 8 in Figure ##FIG##3##4##.</p> <p>Furthermore, the recent sequencing of 46,000 <italic>Citrus clementina </italic>BAC ends (to be published) enabled the construction of a physical map of the 39B3 deleted region. To this end, two DNA sequences covering 700 kb along the <italic>Populus </italic>chromosomes 12 and 15, containing the genes homologous to the <italic>Citrus </italic>deleted candidates, were BLASTed against the <italic>Citrus </italic>BAC end database. The homology search identified 33 BACs with a BLASTN E value lower than 10<sup>-5 </sup>for both paralogous regions. In subsequent analyses, redundant BACs were discarded, while additional candidate BACs were obtained by comparing these previous ones with the BES database to yield overlapping BACs. Moreover, BACs with both ends showing similarity to repetitive DNA that may cause ambiguous positioning and inaccurate gene dosage measurement were also discarded. Finally, a partial physical map containing 13 BACs systematically named B1 to B13 (Table ##TAB##2##3##) was provided by standard PCR of BAC end amplicons against BAC templates and <italic>in silico </italic>search of overlapping antiparallel ends (Figures ##FIG##4##5A, B##).</p> <p>This mapping contained three gaps, one at the 5' deletion junction and two internal ones (Figure ##FIG##4##5B##) delimiting three main BAC clusters, composed of B1 to B4, B5 to B8, and B9 to B13. BACs B11 and B12 were connected by unigene aCL4690Contig1 coding for a putative subunit ClpD of an ATP-dependent Clp protease, whose sequence was shared by both BACs. Similarly B12 and B13 interaction is mediated by unigene aCL1915Contig2 (Table ##TAB##1##2##, ##TAB##2##3##). Real-time PCR quantification of gene dosage for some of the BAC ends (Figure ##FIG##4##5A##) confirmed the presence of these sequences at half dosage in the mutant genotype, indicating that the 39B3 mutation is a hemizygous deletion. Indeed, all analyzed BACs covered an internal segment of the deletion except B13 that exhibited haploid gene dosage on the left end and diploid dosage on the right one, suggesting that B13 contained the 3' border of the 39B3 deletion.</p> <p>The above results indicated that the microsynteny between <italic>Citrus </italic>and <italic>Populus </italic>genomes was high enough to predict gene arrangement and to build a partial physical map of a <italic>Citrus </italic>genomic segment of about 700 kb, as inferred from the length of poplar homologous regions. Nevertheless, the observation that a 700 kb <italic>Citrus </italic>fragment only contains 21 genes may result striking considering an average distance of 10 Kb between adjacent genes, as deduced from the estimations of <italic>Citrus </italic>genome size (367 Mb) and gene number (35,000–40,000). It should be noted, however, that the microarray used in these analyses contains between approximately 2/3 and 1/2 of the estimated gene content of the <italic>Citrus </italic>genome, which may account for a major part of the hypothetical "loss" of deleted candidates. While this is a weakness of the currently available Citrus arrays, non-attributable to the array-CGH procedure, more complete results are expected after the development of a more representative cDNA microarray. Other limitations of the method may be related to the differential hybridization potential of different cDNAs, including for instance cross-hybridizations. In this regard, oligonucleotide arrays are particularly suited for the detection of dissimilar DNA variants. Alternatively, synteny might be limited to several genes located on a bulk of non-conserved sequences inside this 700 Kb region, a possibility that may only be corroborated after genome sequencing.</p> <p>Overall, the data indicated that the <italic>Populus </italic>genome is a useful model for comparative genomics which may be used to characterize hemizygous deletions in <italic>Citrus</italic>.</p> <title>The <italic>Citrus </italic>39B3 deletion shows higher local gene colinearity with <italic>Populus </italic>than with <italic>Arabidopsis</italic></title> <p>Local gene colinearity between two genomic fragments is determined by the number of paralogous genes arranged in the same order. Therefore, not only permanence of genes in their original chromosomal location, but also conservation of gene order, affects local colinearity. In order to validate the gene arrangement postulated in Figure ##FIG##3##4## and consequently to estimate gene colinearity of the 39B3 <italic>Citrus </italic>deletion with <italic>Populus </italic>homologous fragments, we mapped by PCR the 21 genes listed in Table ##TAB##1##2## on the physical map of Figure ##FIG##4##5B##. All but three genes showed at least one PCR product on the array of 13 BACs, confirming that those genes were effectively included in the 39B3 deletion (Figure ##FIG##5##6A##). In addition, PCR reactions reproduced at the BAC size resolution the expected gene order outlined in Figure ##FIG##3##4##, corroborating the gene arrangement deduced by comparative genomics. In Figure ##FIG##5##6B##, the genes rendering a positive PCR signal were linked to the physical map position with an arrow. Moreover, genes 3, 4 and 9 corresponding to unigenes aCL3991Contig1, aC18005F10Rv_c and aC16014F08SK_c, respectively, did not show a detectable PCR signal on purified BACs, although their respective primers produced a band of the expected size when tested against genomic DNA from normal clementine cultivar (data not shown). These genes were most likely placed into the two reported internal gaps of the physical map, as suggested by the border situation of their neighbouring genes.</p> <p>These results confirm high local gene colinearity with poplar in the genomic region covered by 39B3 deletion. Taking together gene content and order conservation (Figures ##FIG##1##2##, ##FIG##2##3## and ##FIG##5##6##), it is inferred that in the studied DNA deleted segment there was higher gene colinearity with <italic>Populus</italic>, which diverged about 109 million years ago (Mya), than with <italic>Arabidopsis</italic>, splitting from the <italic>Citrus </italic>lineage about 87 Mya [##UREF##9##30##], despite gene colinearity generally being correlated with phylogenetic relatedness. A similar conclusion has been reached in our group, after comparing the whole collection of <italic>Citrus </italic>BES with the poplar and <italic>Arabidopsis </italic>genomes (to be published), and also in previous works in papaya and melon. In papaya, BES alignment to the annotated genomes rendered higher gene colinearity with <italic>Populus </italic>than with <italic>Arabidopsis</italic>, although both <italic>Arabidopsis </italic>and papaya belong to the order Brassicales [##REF##16703363##31##]. In melon, microsynteny studies based on the sequence of two BACs also concluded that melon was closer to <italic>Populus </italic>than to <italic>Arabidopsis </italic>or <italic>Medicago truncatula </italic>[##REF##17665215##32##]. These observations may be explained by a differential genome evolutionary dynamics in poplar and <italic>Arabidopsis </italic>lineages [##UREF##10##33##]. The more recent appraisals estimated that last whole genome duplications occurred not later than 60–65 Mya in <italic>Populus </italic>and around 24–40 Mya in <italic>Arabidopsis </italic>lineages [##REF##16973872##28##,##REF##12566392##34##, ####REF##12777046##35##, ##REF##12660784##36####12660784##36##]. Despite the older poplar event, genome rearrangements involving gene loss and translocation following these duplications were much more frequent in <italic>Arabidopsis </italic>ancestors [##REF##11118139##37##]. Such a highly active genome dynamics probably caused the dispersion of genes and the subsequent reduction in synteny and gene colinearity with even related species. The different behaviour of <italic>Populus </italic>and <italic>Arabidopsis </italic>ancestral genomes still deserves further explanation. It has been suggested that woody long-lived species like poplar trees may undergo a slower genome dynamics due to their juvenile period that delays sexual fecundation for several years and to the recurrent contribution of gametes from aged individuals of previous generations [##REF##16973872##28##]. In addition, species like <italic>Arabidopsis thaliana </italic>may have very active mechanisms for unequal or illegitimate recombination causing frequent chromosomal rearrangements such as translocations, insertions and deletions. In this context, it is notable that nearly all <italic>Citrus </italic>species and many related genera have 2n = 18, probably indicating slow chromosomal evolution in this group.</p> <title>Chlorophyll <italic>a/b </italic>ratio is modified in 39B3 mutant</title> <p>Structural studies describing gene arrangement on a particular deletion have outstanding importance for linking a specific mutant phenotype with an impaired gene. The 39B3 deletion removed at least a set of 21 genes and resulted in delayed chlorophyll catabolism. Although in principle, no obvious candidate genes could be unequivocally related to the exocarp colour break retardation, the 39B3 mutant certainly exhibited altered chlorophyll <italic>a </italic>and <italic>b </italic>content. Ratios of chlorophyll <italic>a </italic>to chlorophyll <italic>b </italic>contents in 39B3 mutant were about 15% to 23% lower than those found in wild type when three different green tissues were tested: fruit exocarp, old and young leaves (Figure ##FIG##6##7##). This distinct chlorophyll composition was not accompanied by alterations in the total content of chlorophylls in the leaves although pigment levels in 39B3 fruit exocarp, as expected, were clearly higher (0.48 mg/g fresh weight) than in the peel of control fruit (0.15 mg/g fresh weight) that has initiated chlorophyll degradation (Table ##TAB##3##4##). The chlorophyll accumulation observed in the 39B3 exocarp, however, is higher than the maximum reached in normal clementine fruits (0.35 mg/g f w) [##REF##16787044##38##], suggesting that the mutation also induced total chlorophyll build-up in the fruit peel. Indeed, fruit exocarps of a "wild type" clementine tree showing fruit colour delay due to altered environmental conditions showed chlorophyll <italic>a</italic>/<italic>b </italic>ratios equivalent to those found in the standard cultivar (Figure ##FIG##6##7##) while total pigments had an intermediate value (0.29 mg/g f w) between those of normal and 39B3 genotypes (Table ##TAB##3##4##).</p> <p>Unigene aCL766Contig1, one of the 39B3 hits validated by real-time quantitative PCR (Table ##TAB##0##1##) coding for a ClpC-like protein, may have certain relevance to the altered chlorophyll composition found in 39B3 mutant. Plant ClpCs are ATP-binding proteins located in the stroma of chloroplasts which have been found to be associated with the protein import machinery [##REF##9118955##39##] and with the Clp protease complex [##UREF##11##40##]. In fact, ClpC has been related to protein translocation across the chloroplast inner envelope membrane and to multiple processes requiring proteolytic cleavage, as protein turnover and regulation [##REF##16603408##41##,##REF##11750666##42##]. In <italic>Arabidopsis</italic>, insertional mutagenesis in the <italic>ClpC1 </italic>gene caused chlorosis, growth retardation, photosynthetic damage and defects in chloroplast protein import [##REF##15516497##43##, ####REF##15659100##44##, ##REF##15563614##45####15563614##45##] and no double knock-outs of <italic>ClpC1 </italic>and the less expressed <italic>ClpC2 </italic>genes were obtained, suggesting that <italic>ClpC </italic>function is essential in plants [##REF##17376159##46##]. In addition, a mutant impaired in <italic>ClpC1 </italic>mRNA processing accumulated chlorophyllide <italic>a </italic>oxygenase protein (CAO), a key enzyme for the synthesis of chlorophyll <italic>b </italic>from chlorophyll <italic>a</italic>, leading to a reduced chlorophyll <italic>a</italic>/<italic>b </italic>ratio [##REF##17291312##47##]. Interestingly, aCL4690Contig1 unigene coding for another subunit of Clp complexes (<italic>ClpD</italic>-like) with sequence similarity to aCL766Contig1 showed half gene dosage (Table ##TAB##0##1##) and was also included in the 39B3 deletion. Expression of both <italic>ClpC</italic>-like and <italic>ClpD</italic>-like genes was analyzed in fruit exocarps from wild type and 39B3 mutant at two different developmental stages: green immature peel (September) and shortly after the time of natural colour break in wild type peel (November). Both genes showed reduced expression in the 39B3 mutant, an observation that was well correlated with the alteration in chlorophyll composition since <italic>ClpC</italic>-like and <italic>ClpD-</italic>like alleles in the hemizygous 39B3 mutant reached about a half of the expression values found in the wild type (Figure ##FIG##7##8##). These results suggest that wild type alleles are similarly expressed in the peel of clementine fruit. Furthermore, sequencing of the <italic>ClpC</italic>-like gene that according to its <italic>Arabidopsis </italic>homologue plays a major role in chlorophyll composition, also revealed that there were no essential differences between wild type and 39B3 mutant coding regions (Additional file ##SUPPL##0##1##). This observation corroborated that not only expression but also protein sequence were identical in the analyzed <italic>ClpC</italic>-like alleles. In the wild type, three single nucleotide polymorphisms (SNP) were observed in two different introns, which were very likely unable to alter protein stability or function. These single base variants that were detected as sequence ambiguities (N) were due to the presence of overlapping base peaks contributed by both alleles, while hemizygous <italic>ClpC</italic>-like gene in 39B3 mutant produced an unambiguous signal in the same positions. Bases 2572, 4104 and 4119 located on the forth and sixth introns were identified as guanine, guanine and adenine in the mutant, while a mix of guanine and thymine, guanine and adenine and adenine and thymine were respectively found in wild type DNA (Figure ##FIG##8##9##). While the contribution of the <italic>ClpC</italic>-like gene dosage to the retardation of the natural exocarp degreening remains to be unequivocally demonstrated, the data presented above clearly shows that there is a strong correlation between chlorophyll composition and the presence of a single allele in the mutant. Additional analysis of the remaining genes inside 39B3 deletion should be performed in order to accomplish a complete candidate gene approach</p>
<title>Results and discussion</title> <title>Procedure for the characterization of hemizygous deletions in <italic>Citrus</italic></title> <p>The proposed procedure to identify deleted genes is illustrated in Figure ##FIG##0##1## and its potential to structurally characterize hemizygous deletions is exemplified below with <italic>Citrus </italic>mutants as starting plant material. Its usefulness to describe genomic deletions in other species might be dependent upon genome complexity and ploidy. This method uses cDNA microarrays to hybridize genomic DNA extracted from the deletion mutants to render a list of underrepresented genes. The putative deleted genes are then validated through gene dosage evaluation by real-time PCR using gene specific primers. Deleted genes could subsequently contribute to the identification of the molecular mechanisms underlying the observed phenotypes by means of a candidate gene approach, validated by physiological analyses or genetic transformation [##UREF##8##26##]. In non-sequenced genomes or in plants with poorly developed physical maps, further characterization of deletions at the structural level requires TBLASTX similarity searches against databases containing the sequence annotation of known eudicot genomes, such as <italic>Arabidopsis thaliana</italic>, <italic>Populus trichocarpa </italic>and <italic>Vitis vinifera</italic>. These searches yield putative orthologous genes and syntenic genomic regions between these four species. Local physical maps of deletions are built allocating the deleted gene sequences and the syntenic genomic fragments from these other eudicots into a BES database of the species of interest. Lastly, specific PCR on the array of BACs confirms gene content and order on the lineal structure of the deletions. The results may also be used in comparative genomics analyses to study evolutionary dynamics and phylogenetics.</p> <title>Identification of deleted alleles in 39B3 and 39E7 fast neutron mutants of <italic>Citrus clementina</italic></title> <p>For this study, two mutants obtained by fast neutron mutagenesis of wild type <italic>Citrus clementina </italic>were selected from the IVIA mutant collection. These mutants, named 39B3 and 39E7, were expected to carry DNA deletion lesions in hemizygous dosage and showed a delay in natural colour break in fruit peel. The 39B3 mutant exhibited a delay in colour change from green to orange while 39E7 was better characterized by an abnormal final yellowish colour instead of the natural orange coloration. Putative deleted genes in the mutants were first identified through an approach based on genomic hybridization (array-CGH) that exploited a recently developed <italic>Citrus </italic>microarray containing 21240 cDNAs [##REF##15830128##4##,##REF##17254327##5##]. To this end, total genomic DNA from four independent samples of mutants 39B3 and 39E7 were Cy3 or Cy5-labelled and cohybridized with wild type DNA labelled with the complementary Cy5 or Cy3 probe on four independent microarray slides. Fluorescence intensity data were normalized and single ESTs showing a mutant/wild type signal ratio lower than 0.7 fold, with a P-value lower than 0.2 (39B3) or 0.1 (39E7), were selected as putative candidates.</p> <p>The number of ESTs fulfilling these criteria was 24 and 78 for mutants 39B3 and 39E7, respectively. One of the 39B3 positives [GenBank: <ext-link ext-link-type="gen" xlink:href="CX299090">CX299090</ext-link>], composed of three unrelated sequences was discarded for subsequent analysis due to its chimerical nature. In order to validate the array-CGH results, gene dosage of several putative candidates was determined through real-time PCR quantification of mutant/wild type signals for candidate ESTs as related to a reference undeleted gene [GenBank: <ext-link ext-link-type="gen" xlink:href="CX293764">CX293764</ext-link>]. The results showed that gene dosage for 39B3 candidates ranged from 0.50 to 0.60 when genomic DNA from the 39B3 genotype was tested, while ranged from 0.96 to 1.15 when the assayed DNA originated from the 39E7 genotype (Table ##TAB##0##1##). Similar results, corroborating the presence of putative deleted genes at half dosage, were also obtained for the 39E7 mutant. Therefore, the developed array-CGH procedure proved to be an appropriate tool to identify genes in hemizygous content in the self-incompatible clementine.</p> <title>Clustering of homologues of <italic>Citrus </italic>deleted genes in the poplar genome</title> <p>Microsynteny comparisons with homologous stretches from the sequenced genomes of <italic>Arabidopsis thaliana</italic>, <italic>Populus trichocarpa </italic>and <italic>Vitis vinifera </italic>[##REF##11130711##27##, ####REF##16973872##28##, ##REF##17721507##29####17721507##29##] were performed in order to elucidate hypothetical clustering of <italic>Citrus </italic>deleted genes in the genome. TBLASTX, which searches for translations of a crude genome similar to a translated query, was utilized with an E-value cut-off of 10<sup>-5</sup>. The homologous regions produced by the best TBLASTX hit of each of the <italic>Citrus </italic>candidate genes were located on the chromosome maps of <italic>Arabidopsis</italic>, poplar and grapevine. Homologues of <italic>Citrus </italic>genes were then grouped into clusters in each species when the distance between them was shorter than 250 kb. The second and third TBLASTX best alignments were similarly placed in the respective maps when they were included in an existing cluster. In this case, a binding line was drawn linking the second and third hits to the best hit of the same <italic>Citrus </italic>query. Thus, two chromosomal maps, one for each mutant, in the three species was obtained. Figure ##FIG##1##2## represents in detail chromosome mappings of the 39B3 mutation, which was subjected to further analyses. The results indicate that the <italic>Populus </italic>mapping exhibited rather lower complexity than the <italic>Arabidopsis </italic>and grapevine ones since it included fewer chromosomes and only 3 clusters although the number of 39B3 candidate genes represented in the map was identical (21) for the three genomes. Note that the number of represented hits in these mappings is higher than 21 due to the inclusion of second or third homologues. In <italic>Populus</italic>, most of the candidate genes mapped to two different genome regions of approximately 700 kb long in chromosomes 12 and 15, two duplicated chromosomes that probably originated during the recent genome duplication event that occurred in this species [##REF##16973872##28##]. These two clusters contained 17 and 15 hits respectively while the third one placed in chromosome 16 had only one hit. In contrast, the number of clusters in <italic>Arabidopsis </italic>and <italic>Vitis </italic>were 9 and 11, respectively, and none of them contained more than 11 hits. Furthermore, cluster number (and clustering density) of the homologues of 39E7 putative deleted genes was also lower (and higher) in <italic>Populus </italic>than in <italic>Arabidopsis </italic>or <italic>Vitis</italic>, although the differences were smaller: 26, 30 and 30 clusters were obtained for poplar, <italic>Arabidopsis </italic>and grapevine respectively (Figure ##FIG##2##3##).</p> <p>Overall, these observations suggest that the <italic>Populus </italic>genomic regions homologous to the <italic>Citrus </italic>deletions were less fragmented than their counterparts in <italic>Arabidopsis </italic>and <italic>Vitis</italic>, and consequently microsynteny on the considered segments was higher with the <italic>Populus </italic>genome. These results are striking since <italic>Citrus </italic>and <italic>Arabidopsis </italic>belong to Sapindales and Brassicales orders (inside the same clade eurosids II) while <italic>Populus </italic>is included in the eurosids I clade, and <italic>Vitis </italic>is part of Vitaceae, a family outside of rosids [##UREF##9##30##].</p> <title>Gene arrangement and partial physical map of the 39B3 deletion</title> <p>The closer microsynteny observed between the 39B3 deletion and the two duplicated homologous regions in poplar enabled prediction of gene order by direct inference from the <italic>Populus </italic>sequences. This assumption led to the gene arrangement depicted in Figure ##FIG##3##4##. Twenty genes out of twenty-one having high similarity with <italic>Populus </italic>homologues were directly located on the <italic>Citrus </italic>deletion fragment by combining the two clusters found on <italic>Populus </italic>chromosomes 12 and 15, which shared 12 hits. Inclusion of the 21<sup>st </sup>gene, a homologue of a <italic>Populus </italic>gene placed on chromosome 16, in the 39B3 deletion was based on its location on the right end of the <italic>Citrus </italic>BAC CCER1019D04 (named B12, see below), whose left end shared identity with another deleted gene [GenBank: <ext-link ext-link-type="gen" xlink:href="CX295702">CX295702</ext-link>]. The accession number and protein similarity of these 21 genes, numbered according to the ordered position of their homologues on the poplar genome (Figure ##FIG##3##4##), are depicted in Table ##TAB##1##2## that also shows coding strand sense of poplar homologues. The coding strand was coincident for the <italic>Populus </italic>paralogous genes present in chromosomes 12 and 15, except for genes similar to <italic>Citrus </italic><ext-link ext-link-type="gen" xlink:href="CX308429">CX308429</ext-link>, located in position 8 in Figure ##FIG##3##4##.</p> <p>Furthermore, the recent sequencing of 46,000 <italic>Citrus clementina </italic>BAC ends (to be published) enabled the construction of a physical map of the 39B3 deleted region. To this end, two DNA sequences covering 700 kb along the <italic>Populus </italic>chromosomes 12 and 15, containing the genes homologous to the <italic>Citrus </italic>deleted candidates, were BLASTed against the <italic>Citrus </italic>BAC end database. The homology search identified 33 BACs with a BLASTN E value lower than 10<sup>-5 </sup>for both paralogous regions. In subsequent analyses, redundant BACs were discarded, while additional candidate BACs were obtained by comparing these previous ones with the BES database to yield overlapping BACs. Moreover, BACs with both ends showing similarity to repetitive DNA that may cause ambiguous positioning and inaccurate gene dosage measurement were also discarded. Finally, a partial physical map containing 13 BACs systematically named B1 to B13 (Table ##TAB##2##3##) was provided by standard PCR of BAC end amplicons against BAC templates and <italic>in silico </italic>search of overlapping antiparallel ends (Figures ##FIG##4##5A, B##).</p> <p>This mapping contained three gaps, one at the 5' deletion junction and two internal ones (Figure ##FIG##4##5B##) delimiting three main BAC clusters, composed of B1 to B4, B5 to B8, and B9 to B13. BACs B11 and B12 were connected by unigene aCL4690Contig1 coding for a putative subunit ClpD of an ATP-dependent Clp protease, whose sequence was shared by both BACs. Similarly B12 and B13 interaction is mediated by unigene aCL1915Contig2 (Table ##TAB##1##2##, ##TAB##2##3##). Real-time PCR quantification of gene dosage for some of the BAC ends (Figure ##FIG##4##5A##) confirmed the presence of these sequences at half dosage in the mutant genotype, indicating that the 39B3 mutation is a hemizygous deletion. Indeed, all analyzed BACs covered an internal segment of the deletion except B13 that exhibited haploid gene dosage on the left end and diploid dosage on the right one, suggesting that B13 contained the 3' border of the 39B3 deletion.</p> <p>The above results indicated that the microsynteny between <italic>Citrus </italic>and <italic>Populus </italic>genomes was high enough to predict gene arrangement and to build a partial physical map of a <italic>Citrus </italic>genomic segment of about 700 kb, as inferred from the length of poplar homologous regions. Nevertheless, the observation that a 700 kb <italic>Citrus </italic>fragment only contains 21 genes may result striking considering an average distance of 10 Kb between adjacent genes, as deduced from the estimations of <italic>Citrus </italic>genome size (367 Mb) and gene number (35,000–40,000). It should be noted, however, that the microarray used in these analyses contains between approximately 2/3 and 1/2 of the estimated gene content of the <italic>Citrus </italic>genome, which may account for a major part of the hypothetical "loss" of deleted candidates. While this is a weakness of the currently available Citrus arrays, non-attributable to the array-CGH procedure, more complete results are expected after the development of a more representative cDNA microarray. Other limitations of the method may be related to the differential hybridization potential of different cDNAs, including for instance cross-hybridizations. In this regard, oligonucleotide arrays are particularly suited for the detection of dissimilar DNA variants. Alternatively, synteny might be limited to several genes located on a bulk of non-conserved sequences inside this 700 Kb region, a possibility that may only be corroborated after genome sequencing.</p> <p>Overall, the data indicated that the <italic>Populus </italic>genome is a useful model for comparative genomics which may be used to characterize hemizygous deletions in <italic>Citrus</italic>.</p> <title>The <italic>Citrus </italic>39B3 deletion shows higher local gene colinearity with <italic>Populus </italic>than with <italic>Arabidopsis</italic></title> <p>Local gene colinearity between two genomic fragments is determined by the number of paralogous genes arranged in the same order. Therefore, not only permanence of genes in their original chromosomal location, but also conservation of gene order, affects local colinearity. In order to validate the gene arrangement postulated in Figure ##FIG##3##4## and consequently to estimate gene colinearity of the 39B3 <italic>Citrus </italic>deletion with <italic>Populus </italic>homologous fragments, we mapped by PCR the 21 genes listed in Table ##TAB##1##2## on the physical map of Figure ##FIG##4##5B##. All but three genes showed at least one PCR product on the array of 13 BACs, confirming that those genes were effectively included in the 39B3 deletion (Figure ##FIG##5##6A##). In addition, PCR reactions reproduced at the BAC size resolution the expected gene order outlined in Figure ##FIG##3##4##, corroborating the gene arrangement deduced by comparative genomics. In Figure ##FIG##5##6B##, the genes rendering a positive PCR signal were linked to the physical map position with an arrow. Moreover, genes 3, 4 and 9 corresponding to unigenes aCL3991Contig1, aC18005F10Rv_c and aC16014F08SK_c, respectively, did not show a detectable PCR signal on purified BACs, although their respective primers produced a band of the expected size when tested against genomic DNA from normal clementine cultivar (data not shown). These genes were most likely placed into the two reported internal gaps of the physical map, as suggested by the border situation of their neighbouring genes.</p> <p>These results confirm high local gene colinearity with poplar in the genomic region covered by 39B3 deletion. Taking together gene content and order conservation (Figures ##FIG##1##2##, ##FIG##2##3## and ##FIG##5##6##), it is inferred that in the studied DNA deleted segment there was higher gene colinearity with <italic>Populus</italic>, which diverged about 109 million years ago (Mya), than with <italic>Arabidopsis</italic>, splitting from the <italic>Citrus </italic>lineage about 87 Mya [##UREF##9##30##], despite gene colinearity generally being correlated with phylogenetic relatedness. A similar conclusion has been reached in our group, after comparing the whole collection of <italic>Citrus </italic>BES with the poplar and <italic>Arabidopsis </italic>genomes (to be published), and also in previous works in papaya and melon. In papaya, BES alignment to the annotated genomes rendered higher gene colinearity with <italic>Populus </italic>than with <italic>Arabidopsis</italic>, although both <italic>Arabidopsis </italic>and papaya belong to the order Brassicales [##REF##16703363##31##]. In melon, microsynteny studies based on the sequence of two BACs also concluded that melon was closer to <italic>Populus </italic>than to <italic>Arabidopsis </italic>or <italic>Medicago truncatula </italic>[##REF##17665215##32##]. These observations may be explained by a differential genome evolutionary dynamics in poplar and <italic>Arabidopsis </italic>lineages [##UREF##10##33##]. The more recent appraisals estimated that last whole genome duplications occurred not later than 60–65 Mya in <italic>Populus </italic>and around 24–40 Mya in <italic>Arabidopsis </italic>lineages [##REF##16973872##28##,##REF##12566392##34##, ####REF##12777046##35##, ##REF##12660784##36####12660784##36##]. Despite the older poplar event, genome rearrangements involving gene loss and translocation following these duplications were much more frequent in <italic>Arabidopsis </italic>ancestors [##REF##11118139##37##]. Such a highly active genome dynamics probably caused the dispersion of genes and the subsequent reduction in synteny and gene colinearity with even related species. The different behaviour of <italic>Populus </italic>and <italic>Arabidopsis </italic>ancestral genomes still deserves further explanation. It has been suggested that woody long-lived species like poplar trees may undergo a slower genome dynamics due to their juvenile period that delays sexual fecundation for several years and to the recurrent contribution of gametes from aged individuals of previous generations [##REF##16973872##28##]. In addition, species like <italic>Arabidopsis thaliana </italic>may have very active mechanisms for unequal or illegitimate recombination causing frequent chromosomal rearrangements such as translocations, insertions and deletions. In this context, it is notable that nearly all <italic>Citrus </italic>species and many related genera have 2n = 18, probably indicating slow chromosomal evolution in this group.</p> <title>Chlorophyll <italic>a/b </italic>ratio is modified in 39B3 mutant</title> <p>Structural studies describing gene arrangement on a particular deletion have outstanding importance for linking a specific mutant phenotype with an impaired gene. The 39B3 deletion removed at least a set of 21 genes and resulted in delayed chlorophyll catabolism. Although in principle, no obvious candidate genes could be unequivocally related to the exocarp colour break retardation, the 39B3 mutant certainly exhibited altered chlorophyll <italic>a </italic>and <italic>b </italic>content. Ratios of chlorophyll <italic>a </italic>to chlorophyll <italic>b </italic>contents in 39B3 mutant were about 15% to 23% lower than those found in wild type when three different green tissues were tested: fruit exocarp, old and young leaves (Figure ##FIG##6##7##). This distinct chlorophyll composition was not accompanied by alterations in the total content of chlorophylls in the leaves although pigment levels in 39B3 fruit exocarp, as expected, were clearly higher (0.48 mg/g fresh weight) than in the peel of control fruit (0.15 mg/g fresh weight) that has initiated chlorophyll degradation (Table ##TAB##3##4##). The chlorophyll accumulation observed in the 39B3 exocarp, however, is higher than the maximum reached in normal clementine fruits (0.35 mg/g f w) [##REF##16787044##38##], suggesting that the mutation also induced total chlorophyll build-up in the fruit peel. Indeed, fruit exocarps of a "wild type" clementine tree showing fruit colour delay due to altered environmental conditions showed chlorophyll <italic>a</italic>/<italic>b </italic>ratios equivalent to those found in the standard cultivar (Figure ##FIG##6##7##) while total pigments had an intermediate value (0.29 mg/g f w) between those of normal and 39B3 genotypes (Table ##TAB##3##4##).</p> <p>Unigene aCL766Contig1, one of the 39B3 hits validated by real-time quantitative PCR (Table ##TAB##0##1##) coding for a ClpC-like protein, may have certain relevance to the altered chlorophyll composition found in 39B3 mutant. Plant ClpCs are ATP-binding proteins located in the stroma of chloroplasts which have been found to be associated with the protein import machinery [##REF##9118955##39##] and with the Clp protease complex [##UREF##11##40##]. In fact, ClpC has been related to protein translocation across the chloroplast inner envelope membrane and to multiple processes requiring proteolytic cleavage, as protein turnover and regulation [##REF##16603408##41##,##REF##11750666##42##]. In <italic>Arabidopsis</italic>, insertional mutagenesis in the <italic>ClpC1 </italic>gene caused chlorosis, growth retardation, photosynthetic damage and defects in chloroplast protein import [##REF##15516497##43##, ####REF##15659100##44##, ##REF##15563614##45####15563614##45##] and no double knock-outs of <italic>ClpC1 </italic>and the less expressed <italic>ClpC2 </italic>genes were obtained, suggesting that <italic>ClpC </italic>function is essential in plants [##REF##17376159##46##]. In addition, a mutant impaired in <italic>ClpC1 </italic>mRNA processing accumulated chlorophyllide <italic>a </italic>oxygenase protein (CAO), a key enzyme for the synthesis of chlorophyll <italic>b </italic>from chlorophyll <italic>a</italic>, leading to a reduced chlorophyll <italic>a</italic>/<italic>b </italic>ratio [##REF##17291312##47##]. Interestingly, aCL4690Contig1 unigene coding for another subunit of Clp complexes (<italic>ClpD</italic>-like) with sequence similarity to aCL766Contig1 showed half gene dosage (Table ##TAB##0##1##) and was also included in the 39B3 deletion. Expression of both <italic>ClpC</italic>-like and <italic>ClpD</italic>-like genes was analyzed in fruit exocarps from wild type and 39B3 mutant at two different developmental stages: green immature peel (September) and shortly after the time of natural colour break in wild type peel (November). Both genes showed reduced expression in the 39B3 mutant, an observation that was well correlated with the alteration in chlorophyll composition since <italic>ClpC</italic>-like and <italic>ClpD-</italic>like alleles in the hemizygous 39B3 mutant reached about a half of the expression values found in the wild type (Figure ##FIG##7##8##). These results suggest that wild type alleles are similarly expressed in the peel of clementine fruit. Furthermore, sequencing of the <italic>ClpC</italic>-like gene that according to its <italic>Arabidopsis </italic>homologue plays a major role in chlorophyll composition, also revealed that there were no essential differences between wild type and 39B3 mutant coding regions (Additional file ##SUPPL##0##1##). This observation corroborated that not only expression but also protein sequence were identical in the analyzed <italic>ClpC</italic>-like alleles. In the wild type, three single nucleotide polymorphisms (SNP) were observed in two different introns, which were very likely unable to alter protein stability or function. These single base variants that were detected as sequence ambiguities (N) were due to the presence of overlapping base peaks contributed by both alleles, while hemizygous <italic>ClpC</italic>-like gene in 39B3 mutant produced an unambiguous signal in the same positions. Bases 2572, 4104 and 4119 located on the forth and sixth introns were identified as guanine, guanine and adenine in the mutant, while a mix of guanine and thymine, guanine and adenine and adenine and thymine were respectively found in wild type DNA (Figure ##FIG##8##9##). While the contribution of the <italic>ClpC</italic>-like gene dosage to the retardation of the natural exocarp degreening remains to be unequivocally demonstrated, the data presented above clearly shows that there is a strong correlation between chlorophyll composition and the presence of a single allele in the mutant. Additional analysis of the remaining genes inside 39B3 deletion should be performed in order to accomplish a complete candidate gene approach</p>
<title>Conclusion</title> <p>In this study, we propose a procedure for the genetic characterization of genomic hemizygous deletions in citrus mutants. The procedure that might be applied to other non-sequenced species of similar genome size and ploidy level is illustrated with the study of the 39B3 <italic>Citrus clementina </italic>deletion, generated by fast neutron bombardment. The proposed strategy utilizes several genomic resources such as array-Comparative Genomic Hybridization (array-CGH) technology, EST and BAC end sequencing databases and poplar genome annotation.</p> <p>The array-CGH results led to the conclusion that the 39B3 deletion removed at least 21 genes while a partial physical map of about 700 kb of the deleted region was inferred by comparison of two homologous genomic regions from poplar with a <italic>Citrus </italic>BES database.</p> <p>Structural data including gene content and order in the deletion was utilized for microsynteny and local gene colinearity studies concluding that in the studied region <italic>Citrus </italic>is more similar to <italic>Populus </italic>than to <italic>Arabidopsis</italic>, a phylogenetically closer species. This observation supports previous works on other species and suggests that the <italic>Arabidopsis </italic>lineage underwent a quicker genome evolutionary dynamics than the <italic>Populus </italic>one.</p> <p>Among the deleted alleles, the function of <italic>ClpC</italic>-like, coding for a putative subunit of a protease involved in chlorophyll <italic>b </italic>synthesis was directly related to the mutant phenotype since green mutant tissues had a lower chlorophyll <italic>a</italic>/<italic>b </italic>ratio.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Many fruit-tree species, including relevant <italic>Citrus </italic>spp varieties exhibit a reproductive biology that impairs breeding and strongly constrains genetic improvements. In citrus, juvenility increases the generation time while sexual sterility, inbreeding depression and self-incompatibility prevent the production of homozygous cultivars. Genomic technology may provide citrus researchers with a new set of tools to address these various restrictions. In this work, we report a valuable genomics-based protocol for the structural analysis of deletion mutations on an heterozygous background.</p> <title>Results</title> <p>Two independent fast neutron mutants of self-incompatible clementine (<italic>Citrus clementina </italic>Hort. Ex Tan. cv. Clemenules) were the subject of the study. Both mutants, named 39B3 and 39E7, were expected to carry DNA deletions in hemizygous dosage. Array-based Comparative Genomic Hybridization (array-CGH) using a <italic>Citrus </italic>cDNA microarray allowed the identification of underrepresented genes in these two mutants. Subsequent comparison of citrus deleted genes with annotated plant genomes, especially poplar, made possible to predict the presence of a large deletion in 39B3 of about 700 kb and at least two deletions of approximately 100 and 500 kb in 39E7. The deletion in 39B3 was further characterized by PCR on available <italic>Citrus </italic>BACs, which helped us to build a partial physical map of the deletion. Among the deleted genes, <italic>ClpC</italic>-like gene coding for a putative subunit of a multifunctional chloroplastic protease involved in the regulation of chlorophyll <italic>b </italic>synthesis was directly related to the mutated phenotype since the mutant showed a reduced chlorophyll <italic>a</italic>/<italic>b </italic>ratio in green tissues.</p> <title>Conclusion</title> <p>In this work, we report the use of array-CGH for the successful identification of genes included in a hemizygous deletion induced by fast neutron irradiation on <italic>Citrus clementina</italic>. The study of gene content and order into the 39B3 deletion also led to the unexpected conclusion that microsynteny and local gene colinearity in this species were higher with <italic>Populus trichocarpa </italic>than with the phylogenetically closer <italic>Arabidopsis thaliana</italic>. This work corroborates the potential of <italic>Citrus </italic>genomic resources to assist mutagenesis-based approaches for functional genetics, structural studies and comparative genomics, and hence to facilitate citrus variety improvement.</p>
<title>Authors' contributions</title> <p>GR carried out the microarray hybridizations, standard PCR reactions, expression measurements, gene sequencing, similarity searches and data analysis, and drafted the manuscript. MAN isolated DNA from <italic>Citrus </italic>BACs and carried out quantitative real-time PCR. DJI carried out mutant collection screenings, selected mutants and performed chlorophyll measurements. OR-R designed the array-CGH protocol. MG assisted in microarray hybridizations and data analysis. AU provided plant material, identified mutant genotypes and carried out relevant work on the field. MT conceived and coordinated the project and elaborated the final manuscript. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>Work at Centro de Genómica was supported by INIA grant RTA04-013, INCO contract 015453 and Ministerio de Educación y Ciencia grant AGL2007-65437-C04-01/AGR. We also thank Dr. José Marqués at the Instituto Tecnologico e Nuclear (Sacavem, Portugal) for irradiation with fast neutrons. Help and expertise of A. Almenar, E Blázquez, I. López, I. Sanchís and M. Sancho are gratefully acknowledged.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Schematic guidelines for the characterization of hemizygous deletion <italic>Citrus </italic>mutants</bold>. Arrows indicate successive steps of mutant characterization. Plant material and genomic resources are highlighted in grey boxes while gained knowledge of genetics and genomics are highlighted in white rectangles. Genes are shown in bold and methods and approaches in italics.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Chromosome mapping of poplar, <italic>Arabidopsis </italic>and grapevine homologues of the 39B3 <italic>Citrus </italic>deleted genes</bold>. The first TBLASTX hit for each <italic>Citrus </italic>deleted gene with an E value cut-off &lt; 10<sup>-5 </sup>is represented on linkage groups (LG) from <italic>Populus trichocarpa</italic>, <italic>Arabidopsis thaliana </italic>and <italic>Vitis vinifera</italic>. Homologues of <italic>Citrus </italic>genes were grouped into clusters in each species when the distance between them was shorter than 250 kb. Second and third hits are only represented when they are located in a previously identified cluster, and in this case are linked to the first hit by a line. The value on each cluster indicates the hit number of the cluster.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Chromosome mapping of poplar, <italic>Arabidopsis </italic>and grapevine homologues of the 39E7 deleted genes</bold>. The TBLASTX hits for the 78 39E7 candidates, identified as stated in figure 2, are represented on linkage groups (LG) from <italic>Populus trichocarpa</italic>, <italic>Arabidopsis thaliana </italic>and <italic>Vitis vinifera</italic>. The number of hits contained in each cluster is shown.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Gene composition of the <italic>Citrus </italic>39B3 deletion inferred from poplar homologous regions</bold>. The 39B3 deleted <italic>Citrus </italic>genes are arranged in the centre of the figure in the order inferred from the position of their <italic>Populus </italic>homologues found in linkage groups 12, 15 and 16. Genes are numbered following this order. Strand sense deduced from poplar counterparts is indicated by an arrow.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Local physical mapping of the <italic>Citrus </italic>39B3 deletion</bold>. (A) Electrophoretic analysis of PCR products showing overlapping BACs. Purified BAC templates are distributed horizontally and divided in two panels. Primer pairs were designed from BAC end sequences containing non-repetitive DNA and named with the number of the BAC plus "-L" for left end and "-R" for right end according to the drawing orientation. B7-I primers amplify an internal sequence from B7 instead of an end. Gene dosage measurements for some of the primer pairs in the 39B3 genotype are shown on the left side of the electrophoretic images. (B) Physical map of the <italic>Citrus </italic>39B3 deletion. Horizontal lines represent BACs, which are numbered from left to right. Vertical arrows show overlapping as inferred from PCR reactions and the head of the arrow indicates the BAC template. The vertical lines without arrow show connection of B11 with B12 by sequence of unigene aCL4690Contig1 and B12 with B13 by aCL1915Contig2 instead of a PCR reaction.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p><bold>Gene arrangement on the physical map of the <italic>Citrus </italic>39B3 deletion</bold>. (A) Primer pairs designed for the putative deleted genes included in the 39B3 deletion (Table 2) were utilized in PCR reactions on the BAC templates shown in Table 3. Genes are numbered and arranged vertically, on the left side of the electrophoretic image, and BAC templates are listed horizontally. (B) <italic>Citrus </italic>genes included in the 39B3 deletion and arranged as drawn in Figure 4 but without indication of strand sense, are connected with arrows to the deletion physical map according to PCR results.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p><bold>Chlorophyll <italic>a</italic>/<italic>b </italic>ratio in green tissues from 39B3 mutant and wild type cultivar of <italic>Citrus clementina</italic></bold>. Chlorophyll <italic>a </italic>and <italic>b </italic>content was measured in young and old leaves and fruit peel exocarp from wild type clementine cultivar (Wt) and 39B3 mutant (39B3). Measurements were also taken from exocarps of a normal wild type clementine tree showing fruit colour delay (Wt-d) due to altered environmental conditions. The relative content of chlorophyll <italic>a </italic>to chlorophyll <italic>b </italic>is represented as the Ca/Cb ratio. Data are average of 3 (exocarp) or 5 (leaves) independent determinations and error bars show standard deviation.</p></caption></fig> <fig position="float" id="F8"><label>Figure 8</label><caption><p><bold>Relative expression level of <italic>ClpC</italic>-like and <italic>ClpD</italic>-like genes in 39B3 mutant</bold>. Quantitative real-time PCR with specific primers for <italic>ClpC</italic>-like (A) and <italic>ClpD</italic>-like genes (B) was performed on reverse transcribed RNA from fruit exocarps at two different developmental stages (September and November) from wild type (wt) and 39B3 mutant. Specific first strand cDNA concentration in 39B3 mutant is related to wild type values. The results are average and standard deviation of three independent biological replicates that were assayed twice.</p></caption></fig> <fig position="float" id="F9"><label>Figure 9</label><caption><p><bold>Single nucleotide polymorphisms in <italic>ClpC</italic>-like alleles</bold>. Four colour representations of polymorphic stretches in the sequence of the <italic>ClpC</italic>-like gene in wild type (upper panels) and 39B3 mutant DNA (lower panel). Differential nucleotides are labelled in red. The residues are numbered from the A in the ATG codon.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Gene dosage measurement of deleted genes in 39B3 and 39E7 <italic>Citrus </italic>mutants.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">Unigene</td><td align="center">EST accession number (GenBank)</td><td align="center" colspan="2">Array-CGH</td><td align="center" colspan="2">Real-time PCR gene dosage</td></tr><tr><td/><td/><td colspan="4"><hr/></td></tr><tr><td/><td/><td align="center">Mutant/wt ratio</td><td align="center">P value</td><td align="center">Template 39B3</td><td align="center">Template 39E7</td></tr></thead><tbody><tr><td align="center">aCL4690Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="CX295702">CX295702</ext-link></td><td align="center">0.59</td><td align="center">0.10</td><td align="center">0.56 ± 0.02</td><td align="center">0.99 ± 0.09</td></tr><tr><td align="center">aCL1915Contig2</td><td align="center"><ext-link ext-link-type="gen" xlink:href="DY300024">DY300024</ext-link></td><td align="center">0.62</td><td align="center">0.12</td><td align="center">0.60 ± 0.08</td><td align="center">0.96 ± 0.04</td></tr><tr><td align="center">aCL3317Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="DY265056">DY265056</ext-link></td><td align="center">0.62</td><td align="center">0.10</td><td align="center">0.60 ± 0.03</td><td align="center">0.98 ± 0.03</td></tr><tr><td align="center">aC20009H03SK_c</td><td align="center"><ext-link ext-link-type="gen" xlink:href="CX308429">CX308429</ext-link></td><td align="center">0.63</td><td align="center">0.10</td><td align="center">0.50 ± 0.02</td><td align="center">0.98 ± 0.11</td></tr><tr><td align="center">aCL766Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="CX288964">CX288964</ext-link></td><td align="center">0.65</td><td align="center">0.13</td><td align="center">0.56 ± 0.06</td><td align="center">0.96 ± 0.11</td></tr><tr><td align="center">aCL7097Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="FC868864">FC868864</ext-link></td><td align="center">0.65</td><td align="center">0.12</td><td align="center">0.59 ± 0.04</td><td align="center">1.15 ± 0.07</td></tr><tr><td colspan="6"><hr/></td></tr><tr><td align="center">aCL2087Contig2</td><td align="center"><ext-link ext-link-type="gen" xlink:href="DY300006">DY300006</ext-link></td><td align="center">0.57</td><td align="center">0.05</td><td align="center">1.05 ± 0.13</td><td align="center">0.59 ± 0.11</td></tr><tr><td align="center">aCL6684Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="DY265447">DY265447</ext-link></td><td align="center">0.60</td><td align="center">0.05</td><td align="center">1.04 ± 0.12</td><td align="center">0.59 ± 0.12</td></tr><tr><td align="center">aCL6641Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="FC930062">FC930062</ext-link></td><td align="center">0.62</td><td align="center">0.05</td><td align="center">1.14 ± 0.19</td><td align="center">0.58 ± 0.08</td></tr><tr><td align="center">aCL3902Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="DY267778">DY267778</ext-link></td><td align="center">0.62</td><td align="center">0.05</td><td align="center">1.12 ± 0.10</td><td align="center">0.64 ± 0.04</td></tr><tr><td align="center">aC05139C12SK_c</td><td align="center"><ext-link ext-link-type="gen" xlink:href="CX296347">CX296347</ext-link></td><td align="center">0.66</td><td align="center">0.05</td><td align="center">1.24 ± 0.22</td><td align="center">0.61 ± 0.05</td></tr><tr><td align="center">aC01019E12SK_c</td><td align="center"><ext-link ext-link-type="gen" xlink:href="CX288357">CX288357</ext-link></td><td align="center">0.66</td><td align="center">0.05</td><td align="center">1.04 ± 0.07</td><td align="center">0.63 ± 0.11</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Gene components of the <italic>Citrus </italic>39B3 deletion.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">N°</td><td align="left">Citrus unigene</td><td align="center">EST accession number (GenBank)</td><td align="center">Strand</td><td align="center">Similarity</td></tr></thead><tbody><tr><td align="center">1</td><td align="left">aC01006D04SK_c</td><td align="center"><ext-link ext-link-type="gen" xlink:href="CX287243">CX287243</ext-link></td><td align="center">-</td><td align="left">Hypothetical protein</td></tr><tr><td align="center">2</td><td align="left">aC20006C06SK_c</td><td align="center"><ext-link ext-link-type="gen" xlink:href="CX308114">CX308114</ext-link></td><td align="center">-</td><td align="left">Ubiquitin conjugating enzyme</td></tr><tr><td align="center">3</td><td align="left">aCL3991Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="DY278065">DY278065</ext-link></td><td align="center">-</td><td align="left">Sterile alpha motif (SAM) domain-containing protein</td></tr><tr><td align="center">4</td><td align="left">aC18005F10Rv_c</td><td align="center"><ext-link ext-link-type="gen" xlink:href="CX305429">CX305429</ext-link></td><td align="center">-</td><td align="left">Sialyltransferase-like protein</td></tr><tr><td align="center">5</td><td align="left">aCL3317Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="DY265056">DY265056</ext-link></td><td align="center">+</td><td align="left">Hypothetical protein</td></tr><tr><td align="center">6</td><td align="left">aCL766Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="CX288964">CX288964</ext-link></td><td align="center">-</td><td align="left">ATP-dependent Clp protease, clpC homolog</td></tr><tr><td align="center">7</td><td align="left">aC01012C02SK_c</td><td align="center"><ext-link ext-link-type="gen" xlink:href="CX287682">CX287682</ext-link></td><td align="center">-</td><td align="left">Alpha-mannosidase</td></tr><tr><td align="center">8</td><td align="left">aC20009H03SK_c</td><td align="center"><ext-link ext-link-type="gen" xlink:href="CX308429">CX308429</ext-link></td><td align="center">?</td><td align="left">Mei2-like protein</td></tr><tr><td align="center">9</td><td align="left">aC16014F08SK_c</td><td align="center"><ext-link ext-link-type="gen" xlink:href="CX304691">CX304691</ext-link></td><td align="center">-</td><td align="left">Putative pol polyprotein</td></tr><tr><td align="center">10</td><td align="left">aCL6210Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="DY282423">DY282423</ext-link></td><td align="center">-</td><td align="left">Hypothetical protein</td></tr><tr><td align="center">11</td><td align="left">aCL8592Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="DY267639">DY267639</ext-link></td><td align="center">-</td><td align="left">Tudor domain-containing protein</td></tr><tr><td align="center">12</td><td align="left">aCL1065Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="DY282340">DY282340</ext-link></td><td align="center">-</td><td align="left">Putative amidase</td></tr><tr><td align="center">13</td><td align="left">aC32108G01EF_c</td><td align="center"><ext-link ext-link-type="gen" xlink:href="FC921733">FC921733</ext-link></td><td align="center">+</td><td align="left">Hypothetical protein</td></tr><tr><td align="center">14</td><td align="left">aCL503Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="CX292510">CX292510</ext-link></td><td align="center">-</td><td align="left">Respiratory burst oxidase homolog</td></tr><tr><td align="center">15</td><td align="left">aCL6269Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="DY263746">DY263746</ext-link></td><td align="center">+</td><td align="left">FHA domain-containing protein</td></tr><tr><td align="center">16</td><td align="left">aCL7097Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="FC868864">FC868864</ext-link></td><td align="center">+</td><td align="left">Putative pentatricopeptide (PPR) repeat protein</td></tr><tr><td align="center">17</td><td align="left">aCL4690Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="CX295702">CX295702</ext-link></td><td align="center">+</td><td align="left">ERD1 protein, chloroplast precursor</td></tr><tr><td align="center">18</td><td align="left">aCL8011Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="FC923875">FC923875</ext-link></td><td align="center">-</td><td align="left">Fe-superoxide dismutase</td></tr><tr><td align="center">19</td><td align="left">aIC0AAA60DF12RM1_c</td><td align="center"><ext-link ext-link-type="gen" xlink:href="DY284274">DY284274</ext-link></td><td align="center">+</td><td align="left">Poly(A)-binding protein II-like</td></tr><tr><td align="center">20</td><td align="left">aCL1848Contig1</td><td align="center"><ext-link ext-link-type="gen" xlink:href="FC875470">FC875470</ext-link></td><td align="center">+</td><td align="left">Hypothetical protein</td></tr><tr><td align="center">21</td><td align="left">aCL1915Contig2</td><td align="center"><ext-link ext-link-type="gen" xlink:href="DY300024">DY300024</ext-link></td><td align="center">-</td><td align="left">Tubulin-specific chaperone C-related</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Listing of BACs included in the <italic>Citrus </italic>39B3 deletion.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">N°</td><td align="left">BAC</td><td align="center">Ends</td><td align="center">BES ID (GenBank)</td><td align="center">BLASTX against plant proteins</td><td align="center">E value</td></tr></thead><tbody><tr><td align="center">B1</td><td align="left">CCL021E18</td><td align="center">B1-L</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET070583">ET070583</ext-link></td><td align="left">Nhf</td><td align="center">--</td></tr><tr><td/><td/><td align="center">B1-R</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET070584">ET070584</ext-link></td><td align="left">gi| 91805627| hypothetical protein</td><td align="center">9e-24</td></tr><tr><td align="center">B2</td><td align="left">CCL011O24</td><td align="center">B2-L</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET086992">ET086992</ext-link></td><td align="left">Nhf</td><td align="center">--</td></tr><tr><td/><td/><td align="center">B2-R</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET086991">ET086991</ext-link></td><td align="left">gi| 7576215| hypothetical protein</td><td align="center">7e-47</td></tr><tr><td align="center">B3</td><td align="left">CCER1037B12</td><td align="center">B3-L</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET077105">ET077105</ext-link></td><td align="left">gi| 7576215| hypothetical protein</td><td align="center">3e-90</td></tr><tr><td/><td/><td align="center">B3-R</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET077106">ET077106</ext-link></td><td align="left">Nhf</td><td align="center">--</td></tr><tr><td align="center">B4</td><td align="left">CCER1032N17</td><td align="center">B4-L</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET101817">ET101817</ext-link></td><td align="left">gi| 25411577| probable retroelement pol polyprotein</td><td align="center">2e-06</td></tr><tr><td/><td/><td align="center">B4-R</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET101816">ET101816</ext-link></td><td align="left">Nhf</td><td align="center">--</td></tr><tr><td align="center">B5</td><td align="left">CCL011N15</td><td align="center">B5-L</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET087145">ET087145</ext-link></td><td align="left">gi| 6469119| mitochondrial phosphate transporter</td><td align="center">5e-56</td></tr><tr><td/><td/><td align="center">B5-R</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET087144">ET087144</ext-link></td><td align="left">Nhf</td><td align="center">--</td></tr><tr><td align="center">B6</td><td align="left">CCER1045A09</td><td align="center">B6-L</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET077286">ET077286</ext-link></td><td align="left">gi| 92895029| Polynucleotidyl transferase (retrotransposon protein)</td><td align="center">8e-63</td></tr><tr><td/><td/><td align="center">B6-R</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET077285">ET077285</ext-link></td><td align="left">gi| 30027167| auxin response factor-like protein</td><td align="center">6e-85</td></tr><tr><td align="center">B7</td><td align="left">CCH3037D01</td><td align="center">B7-L</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET112059">ET112059</ext-link></td><td align="left">gi| 87240692| Helix-loop-helix DNA-binding</td><td align="center">1e-21</td></tr><tr><td align="center">B8</td><td align="left">CCER1005N09</td><td align="center">B8-L</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET079746">ET079746</ext-link></td><td align="left">Nhf</td><td align="center">--</td></tr><tr><td/><td/><td align="center">B8-R</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET079745">ET079745</ext-link></td><td align="left">gi| 79331867| AML1; RNA binding/nucleic acid binding</td><td align="center">2e-09</td></tr><tr><td align="center">B9</td><td align="left">CCH3005L04</td><td align="center">B9-L</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET081228">ET081228</ext-link></td><td align="left">gi| 33113977| putative copia-type pol polyprotein</td><td align="center">2e-85</td></tr><tr><td/><td/><td align="center">B9-R</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET081227">ET081227</ext-link></td><td align="left">gi| 51968598| peroxisomal Ca-dependent solute carrier-like protein</td><td align="center">2e-21</td></tr><tr><td align="center">B10</td><td align="left">CCER1033B14</td><td align="center">B10-L</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET102435">ET102435</ext-link></td><td align="left">gi| 51968598| peroxisomal Ca-dependent solute carrier-like protein</td><td align="center">2e-37</td></tr><tr><td/><td/><td align="center">B10-R</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET102434">ET102434</ext-link></td><td align="left">Nhf</td><td align="center">--</td></tr><tr><td align="center">B11</td><td align="left">CCL011K21</td><td align="center">B11-L</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET086761">ET086761</ext-link></td><td align="left">gi| 25402907| protein F5M15.26 (retrotransposon protein)</td><td align="center">4e-78</td></tr><tr><td/><td/><td align="center">B11-R</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET086760">ET086760</ext-link></td><td align="left">gi| 14334878| putative ATP-dependent Clp protease ClpD</td><td align="center">4e-57</td></tr><tr><td align="center">B12</td><td align="left">CCER1019D04</td><td align="center">B12-L</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET098996">ET098996</ext-link></td><td align="left">gi| 14334878| putative ATP-dependent Clp protease ClpD</td><td align="center">4e-35</td></tr><tr><td/><td/><td align="center">B12-R</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET098995">ET098995</ext-link></td><td align="left">gi| 6729532| putative protein</td><td align="center">3e-28</td></tr><tr><td align="center">B13</td><td align="left">CCL032E17</td><td align="center">B13-L</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET094320">ET094320</ext-link></td><td align="left">gi| 6729532| putative protein</td><td align="center">9e-34</td></tr><tr><td/><td/><td align="center">B13-R</td><td align="center"><ext-link ext-link-type="gen" xlink:href="ET094321">ET094321</ext-link></td><td align="left">Nhf</td><td align="center">--</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>Total chlorophyll content in green tissues from 39B3 mutant and "wild type" cultivar of <italic>Citrus clementina</italic>.</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="center">Young leaves</td><td align="center">Old leaves</td><td align="center">Exocarp</td></tr></thead><tbody><tr><td align="center">Wt</td><td align="center">0.76 ± 0.25</td><td align="center">2.23 ± 0.19</td><td align="center">0.15 ± 0.05</td></tr><tr><td align="center">39B3</td><td align="center">0.71 ± 0.24</td><td align="center">2.49 ± 0.17</td><td align="center">0.48 ± 0.01</td></tr><tr><td align="center">Wt-d</td><td align="center">--</td><td align="center">--</td><td align="center">0.29 ± 0.01</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Nucleotide sequence of <italic>ClpC</italic>-like coding region plus introns in wild type and 39B3 mutant.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Primers employed in real-time and standard PCR experiments.</p></caption></supplementary-material>
<table-wrap-foot><p>Accession numbers and the corresponding assembled unigenes for several candidate ESTs are shown. Array-CGH data are signal ratios of mutant with respect to wild type samples for each single EST. Gene dosage was calculated through real-time PCR on 39B3 and 39E7 genomic templates. A horizontal line separates array-CGH hits of 39B3 mutant (upper) from those of the 39E7 mutant (lower).</p></table-wrap-foot> <table-wrap-foot><p>Unigenes are numbered according to the ordered position of their homologues on the poplar genome. The accession numbers correspond to 39B3 candidate ESTs. The coding frame of <italic>Populus </italic>homologues follow the proposed strand (+) or the complementary reverse one (-).</p></table-wrap-foot> <table-wrap-foot><p>BACs are numbered according to the ordered position in the deletion from B1 to B13. BES are named with the number of the BAC plus "-L" for left end and "-R" for right end according to the drawing orientation in Figure ##FIG##4##5b##. Nhf: no hits found.</p></table-wrap-foot> <table-wrap-foot><p>Total chlorophyll content (mg/g fresh weight) was measured in young and old leaves and fruit peel exocarp from samples shown in Figure ##FIG##6##7##. Data are average of 3 (exocarp) or 5 (leaves) independent determinations. Standard deviation is shown.</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-381-1"/> <graphic xlink:href="1471-2164-9-381-2"/> <graphic xlink:href="1471-2164-9-381-3"/> <graphic xlink:href="1471-2164-9-381-4"/> <graphic xlink:href="1471-2164-9-381-5"/> <graphic xlink:href="1471-2164-9-381-6"/> <graphic xlink:href="1471-2164-9-381-7"/> <graphic xlink:href="1471-2164-9-381-8"/> <graphic xlink:href="1471-2164-9-381-9"/>
<media xlink:href="1471-2164-9-381-S1.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-381-S2.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Iglesias", "Cercos", "Colmenero-Flores", "Naranjo", "Rios", "Carrera", "Ruiz-Rivero", "Lliso", "Morillon", "Tadeo", "Talon"], "given-names": ["DJ", "M", "JM", "MA", "G", "E", "O", "I", "R", "FR", "M"], "article-title": ["Physiology of citrus fruiting"], "source": ["Braz J Plant Physiol"], "year": ["2008"]}, {"surname": ["Tadeo", "Cercos", "Colmenero-Flores", "Iglesias", "Naranjo", "Morillon", "Carrera", "Lliso", "Rios", "Ruiz-Rivero", "Talon"], "given-names": ["FR", "M", "JM", "DJ", "MA", "R", "E", "I", "G", "O", "M"], "article-title": ["Molecular physiology of development and quality of citrus"], "source": ["Adv Bot Res"], "year": ["2008"]}, {"surname": ["Bausher", "Shatters", "Chaparro", "Dang", "Hunter", "Niedz"], "given-names": ["M", "R", "J", "P", "W", "R"], "article-title": ["An expressed sequence tag (EST) set from Citrus sinensis L. Osbeck whole seedlings and the implications of further perennial source investigations"], "source": ["Plant Sci"], "year": ["2003"], "volume": ["165"], "fpage": ["415"], "lpage": ["422"], "pub-id": ["10.1016/S0168-9452(03)00202-4"]}, {"surname": ["Shimada", "Kita", "Endo", "Fujii", "Ueda", "Moriguchi", "Omura"], "given-names": ["T", "M", "T", "H", "T", "T", "M"], "article-title": ["Expressed sequence tags of ovary tissue cDNA library in Citrus unshiu Marc"], "source": ["Plant Sci"], "year": ["2003"], "volume": ["165"], "fpage": ["167"], "lpage": ["168"], "pub-id": ["10.1016/S0168-9452(03)00155-9"]}, {"surname": ["Fujii", "Shimada", "Sugiyama", "Nishikawa", "Endo", "Nakano", "Ikoma", "Shimizu", "Omura"], "given-names": ["H", "T", "A", "F", "T", "M", "Y", "T", "M"], "article-title": ["Profiling ethylene-responsive genes in mature mandarin fruit using a citrus 22K oligoarray"], "source": ["Plant Science"], "year": ["2007"], "volume": ["173"], "fpage": ["340"], "lpage": ["348"], "pub-id": ["10.1016/j.plantsci.2007.06.006"]}, {"surname": ["Shimada", "Fujii", "Endo", "Yazaki", "Kishimoto", "Shimbo", "Kikuchi", "Omura"], "given-names": ["T", "H", "T", "J", "N", "K", "S", "M"], "article-title": ["Toward comprehensive expression profiling by microarray analysis in citrus: monitoring the expression profiles of 2213 genes during fruit development"], "source": ["Plant Sci"], "year": ["2005"], "volume": ["168"], "fpage": ["1383"], "lpage": ["1385"], "pub-id": ["10.1016/j.plantsci.2005.01.013"]}, {"surname": ["Li", "Lassner", "Zhang"], "given-names": ["X", "M", "YL"], "article-title": ["Deleteagene: a fast neutron deletion mutagenesis-based gene knockout system for plants"], "source": ["Compar Funct Genom"], "year": ["2002"], "volume": ["3"], "fpage": ["158"], "lpage": ["160"], "pub-id": ["10.1002/cfg.148"]}, {"surname": ["Fan", "Vibranovski", "Chen", "Long"], "given-names": ["CZ", "MD", "Y", "MY"], "article-title": ["A microarray based genomic hybridization method for identification of new genes in plants: Case analyses of Arabidopsis and Oryza"], "source": ["J Integr Plant Biol"], "year": ["2007"], "volume": ["49"], "fpage": ["915"], "lpage": ["926"], "pub-id": ["10.1111/j.1744-7909.2007.00503.x"]}, {"surname": ["Pflieger", "Lefebvre", "Causse"], "given-names": ["S", "V", "M"], "article-title": ["The candidate gene approach in plant genetics: a review"], "source": ["Mol Breeding"], "year": ["2001"], "volume": ["7"], "fpage": ["275"], "lpage": ["291"], "pub-id": ["10.1023/A:1011605013259"]}, {"surname": ["Wikstrom", "Savolainen", "Chase"], "given-names": ["N", "V", "MW"], "article-title": ["Evolution of the angiosperms: calibrating the family tree"], "source": ["P Roy Soc Lond B Bio"], "year": ["2001"], "volume": ["268"], "fpage": ["2211"], "lpage": ["2220"], "pub-id": ["10.1098/rspb.2001.1782"]}, {"surname": ["Kellogg", "Bennetzen"], "given-names": ["EA", "JL"], "article-title": ["The evolution of nuclear genome structure in seed plants"], "source": ["Am J Bot"], "year": ["2004"], "volume": ["91"], "fpage": ["1709"], "lpage": ["1725"], "pub-id": ["10.3732/ajb.91.10.1709"]}, {"surname": ["Desimone", "WeissWichert", "Wagner", "Altenfeld", "Johanningmeier"], "given-names": ["M", "C", "E", "U", "U"], "article-title": ["Immunochemical studies on the Clp-protease in chloroplasts: Evidence for the formation of a ClpC/P complex"], "source": ["Bot Acta"], "year": ["1997"], "volume": ["110"], "fpage": ["234"], "lpage": ["239"]}, {"surname": ["PO"], "given-names": ["B"], "article-title": ["Preparation of fluorescent DNA probe from HUMAN mRNA or total RNA using direct incorporation"]}, {"article-title": ["R statistical computing software "]}, {"article-title": ["ArrayExpress database "]}, {"article-title": ["NCBI BLAST home "]}, {"article-title": ["Populus trichocarpa genome v1.1 BLAST page "]}, {"article-title": ["Grape genome BLAST server "]}]
BMC Genomics. 2008 Aug 9; 9:381