Karyotype diversity and 2C DNA content in species of the Caesalpinia group

Polliana Silva Rodrigues
Margarete Magalhães Souza
Cláusio Antônio
Ferreira Melo
Telma Nair
Santana Pereira
Ronan Xavier Corrêa
Karyotype diversity and 2C DNA content in species of the Caesalpinia group
10.1186/s12863-018-0610-2R E S E A R C H A R T I C L E Open Access
Background: The Leguminosae family is the third-largest family of angiosperms, and Caesalpinioideae is its secondlargest subfamily. A great number of species (approximately 205) are found in the Caesalpinia group within this subfamily; together with these species' phenotypic plasticity and the similarities in their morphological descriptors, make this a complex group for taxonomic and phylogenetic studies. The objective of the present work was to evaluate the karyotypic diversity and the 2C DNA content variation in 10 species of the Caesalpinia group, representing six genera: Paubrasilia, Caesalpinia, Cenostigma, Poincianella, Erythrostemon and Libidibia. The GC-rich heterochromatin and 45S rDNA sites (which are used as chromosome markers) were located to evaluate the karyotype diversity in the clade. The variation in the 2C DNA content was determined through flow cytometry. Results: The fluorochrome banding indicated that the chromomycin A 3 + /4′,6-diamidino-2-phenylindole − blocks were exclusively in the terminal regions of the chromosomes, coinciding with 45S rDNA sites in all analyzed species. Physical mapping of the species (through fluorescence in situ hybridization) revealed variation in the size of the hybridization signals and in the number and distribution of the 45S rDNA sites. All hybridization sites were in the terminal regions of the chromosomes. In addition, all species had a hybridization site in the fourth chromosome pair. The 2C DNA content ranged from 1.54 pg in Erythrostemon calycina to 2.82 pg in the Paubrasilia echinata large-leaf variant. The Pa. echinata small-leaf variant was isolated from the other leaf variants through Scoot-Knott clustering. Conclusions: The chromosome diversity and the variation in the 2C DNA content reinforce that the actual taxonomy and clustering of the analyzed taxa requires more genera that were previously proposed. This fact indicates that taxonomy, phylogeny and cytoevolutionary inference related to the complex Caesalpinia group have to be done through integrative evaluation.

Background

Leguminosae is the third-largest family among angiosperms [1], and Caesalpinioideae, its second-largest subfamily, is represented by about 170 genera, many of them with complex and confused taxonomies. This subfamily's phenotypic plasticity is a challenge for taxonomies that are based on morphology [2]. This group commonly occurs in Brazil, which is home to about 790 described species [3,4]. The Caesalpinia group within the Caesalpinioideae subfamily is a pantropical clade that presents about 205 pantropical species [5], including important Brazilian species that are threatened with extinction [6].

Taxonomic and phylogenetic changes have been done for some species and genera from Leguminosae, including the clustering of taxa in a new genetic system for the Caesalpinia group [5]. The major problem in the Caesalpinia clade taxonomy and its phylogenetic classification relates to morphological similarities, as there is little variation for some descriptors [7]. Solving this problem requires a broad mode of analysis for the taxa characterization [8,9], which has been helpful for the systematic distribution and taxonomy of the Caesalpinia group [10,11].

The banding obtained from the application of chromomycin A 3 (CMA 3 ) and 4′,6-diamidino-2-phenylindole (DAPI) fluorochromes, and from the localization of 45S rDNA sites using the fluorescence in situ hybridization (FISH) technique, has been used to identify specific sites; the positions and sizes of such tags may be useful as cytological markers. These data allow us to define the location and quantity of the chromosome markers that are commonly observed in a group of species, as well as the specific chromosome pattern of the markers for each species [12].

Flow cytometry has been used in biosystematics analyses, mainly to provide results regarding nuclear DNA content and, consequently, the level of ploidy. This allows for better species detection and delimitation, which is helpful in the study of a particular genus's phylogenetic relationships and evolutionary characteristics [13].

Previous studies in which fluorochrome staining was applied to the Caesalpinioideae subfamily revealed interand intraspecific differences. The heterochromatic blocks observed (CMA 3 + /DAPI − ) were distributed in regions proximal to the nuclear organizer regions, but the presence of CMA 3 + /DAPI − blocks was also observed in previous studies [14]. The DNA content indicated the existence of intra-and interspecific variability in some genera within Fabaceae [15][16][17]. However, these analyses included only one DNA-content analysis for a species of the genus Caesalpinia (Caesalpinia crista) [17].

This study aimed to evaluate the karyotype diversity in 10 species (representing six genera) of the pantropical Caesalpinia clade, using the location of GC-rich heterochromatin and the number and position of 45S rDNA sites. In addition, 2C DNA content was quantified using flow cytometry.

Methods

Botanical material and pretreatment

Seeds of 10 species of the Caesalpinia group were collected from several locations in the state of Bahia in Brazil ( Table 1). The seeds were randomly collected, with the name of obtaining as many species from the state as possible. After field collection, the seeds were treated with Captan (Fersol®) fungicide and germinated on moistened filter paper in a humid chamber at room temperature. Root tips of approximately 3 mm in length were collected shortly after germination and pretreated with an antimitotic solution of 0.002 M 8-hydroxyquinoline for 6 h; the root tips were then washed twice in distilled water, dried on filter paper fixed in Carnoy I (3:1 glacial acetic acid to absolute ethanol, v/v) [18] for 2 h at room The averages for the 2C DNA followed by the same letter did not differ statistically from one another by the Scott-Knott group at 0.05% with DMS: 0.0647 b During the duration of this study seeds of the medium and large variants have not been found and cytometry analyses were performed on seedlings obtained by donation temperature, and maintained at − 20°C until the time of use. After radicles were collected, the seedlings were planted in 2 kg bags with organic soil and monitored in a greenhouse, resulting in matrix plants for the cytogenetic characterization.

Preparation of slides and banding with CMA 3

and DAPI fluorochromes

For the localization of base regions specific to GC and AT, the fluorochromes CMA 3 and DAPI were used in a double-staining process. Distamycin A solution was added to the cytological preparation. This protocol followed the one proposed by Guerra and Souza [19], with some modifications. The slides were prepared through enzymatic digestion with 2% cellulase and 20% pectinase for 1 h; this was followed by maceration in a drop of 45% acetic acid and then by freezing in liquid nitrogen to remove the cover slip. The slides containing the cytological preparations were aged for 3 days at room temperature, after which 0.25 mg − 1 of CMA 3 was added for 1 h; this was followed by washing with distilled water and air-drying. Next, 0.1 mg − 1 distamycin A was added for 30 min, followed by another round of washing with distilled water and air-drying. Finally, DAPI was added for 30 min, followed by a last round of washing in distilled water and air-drying. The slides were assembled with 20 × 20 mm cover slips with 1:1 glycerol/McIlvaine medium (v/v), plus 2.5 mM MgCl 2 . After the application of the double staining, the slides were aged for another three days before analysis with epifluorescence microscopy.

Fluorescent in situ hybridization

The application of FISH was performed following the protocol developed by Souza et al. [20], with some modifications, such as eliminating the pretreatment of the slides and adding digestion with pepsin (for better interference of the cytoplasm and cellular walls). The cytological preparations were digested with RNAse (100 μg/ml) and washed twice in 2xSSC (salt, sodium citrate) for 5 min. Next, 50 μl HCl (10 mM) was added and incubated for 5 min at room temperature. After removal of the cover slip, 50 μL of pepsin solution was added (0.75 μL of pepsin and 49.25 μL HCl, 10 mM), and the slide was kept in a humid chamber at 37°C for 20 min. Next, the following steps were carried out: two washes with 2xSSC (5 min each); incubation in 4% paraformaldehyde for 10 min; two washes with 2xSSC (5 min each); and dehydration in an alcoholic series (70% ethanol and 95% ethanol; 5 min each). The slide was then air-dried for at least 30 min. The hybridization mixture was composed of 100% formamide (7.5 μL), 50% dextran (3.0 μL), 20xSSC (1.5 μL), 10% sodium dodecyl sulfate (0.2 μL) and the 45S probe (2.8 μL). This mix was heated in a thermocycler at 75°C for 10 min, transferred to ice for at least 2 min and then placed on a slide, which was then denatured in a thermocycler at 75°C for 10 min and placed in a humid chamber at 37°C overnight. For the post-hybridization baths, the slide was washed with 2xSSC at room temperature, followed by two washes in 2xSSC at 42°C (5 min each), two washes in 0.1xSSC at 42°C (5 min each), two washes in 2xSSC at 42°C (5 min each) and one wash in 4xSSC/0.2% Tween20 at room temperature (5 min). For detection, 50 μl of 5% BSA was applied to the slide for 10 min at room temperature, an antibody solution containing 0.7 μL of avidin-fluorescein isothiocyanate and 19.3 μL of 5% BSA was then added to the slide, which was kept in a dark, humid chamber for 1 h at 37°C. Three washes were performed in 4xSSC/0.2% Tween 20 at room temperature, while still in the dark. Excess 4xSSC/0.2% Tween 20 was removed with a blade rinse in 2xSSC, and slide assembly was completed with 15 μL of DAPI-conjugated Vectashield® (Vector® Laboratories). The slides were refrigerated in a dark container for at least 24 h. Blade analysis was performed using the Olympus® BX41 fluorescence microscope; the images were captured with a DP25 digital camera and DP2-BSW software from Olympus®. The overlap of the images and the drawing of the boards were completed using Adobe Photoshop® software.

Analysis of the 2C DNA

Five plants of each analyzed species were used in the analysis of the 2C DNA. For this, five leaves of each species were sampled for the analysis. The species Zea mays CV Kukurice (with 2C = 5.43 pg of DNA) and Glycine max L. (with 2C = 2.50 pg of DNA) [21] were used as internal standards to estimate the species' genome sizes. The Zea mays species was used as an internal standard for the DNA content of all species except Poincianella pluviosa, for which Glycine max was used as the standard. Suspensions of intact nuclei were prepared using the Cystain PI Absolut P kit (Partec®). About 17 mg of leaf tissue from the target species and 20 mg of leaf tissue from the standard were minced simultaneously on a slide in Petri dish with 1 mL of extraction buffer. The suspension material was filtered through a nylon mesh screen of 50 μm. Then, 2 mL of solution containing RNAse and propidium iodide was added, and the material was incubated in a light-protected vessel for at least 30 min at room temperature. The evaluation of the 2C nuclear DNA was conducted using the Partec® PAII flow cytometer. The gain parameter was adjusted so that the peak for the nuclei of target species G1 was positioned over channel 50. At least 10,000 nuclei were analyzed for each sample. The fluorescence intensity of the nuclei, after staining with propidium iodide, was analyzed at rates of 20-50 nuclei/s. The positions of the peaks, their areas and their coefficients of variation were obtained from the cytometer. The size of the nuclear genome was calculated according to Dolezel [22]:
2CDNA ¼ Average peak G0=G1 for Caesalpinia Average peak G0=G1 for standard x 2C DNA standard pg ð Þ
Analysis of variance (ANOVA) was used for the evaluation of significant differences in the flow cytometry data using a completely randomized design and five repetitions for each species. Additionally, the mean of the 2C DNA was clustered using Scott-Knott clustering. The ANOVA and the mean clustering were done using Sisvar software [23]. The editing of the Partec® flow cytometer histograms was carried out using Corel Draw® X7 software.

Results

Localization of GC-rich heterochromatin

The location of GC-rich heterochromatin using base-specific fluorochromes revealed CMA 3 + /DAPI − terminal blocks in all analyzed species. However,CMA 3 + pericentromeric blocks were also observed in metaphase chromosomes in Libidibia ferrea and Po. microphylla (Figs. 1 and 2  The application of FISH to localize the 45S rDNA sites allowed for the visualization of 45S rDNA in pairs of chromosomes (three, four or five, depending on the species), along with terminal 45S rDNA hybridization sites (Figs. 1 and 2). The karyotypes of the species C. macrophyllum, Po. pluviosa and L. ferrea each had three chromosome pairs with 45S rDNA hybridization sites  (Fig. 1). The species Po. bracteosa, Po. laxiflora, Po. microphylla, Ca. pulcherrima and Po. pyramidalis each had four chromosome pairs with 45S rDNA sites (Fig. 2). Five chromosome pairs with hybridization sites for 45S rDNA were observed in E. calycina (Fig. 1f ) and Pa. echinata (Fig. 2r). All species demonstrated 45S rDNA hybridization sites in the fourth chromosome pair, and only Ca. pulcherrima did not show a marking on the seventh chromosome pair. The 45S rDNA hybridization sites in the eighth and tenth chromosome pairs were limited to Po. microphylla and Po. pyramidalis, respectively. The most frequent markers were located in the second chromosome pair (E. calycina, L. ferrea and Po. pyramidalis), in the fifth chromosome pair (E. calycina, Po. bracteosa and Po. microphylla) and in the eleventh chromosome pair (E. calycina, Pa. echinata and Ca. pulcherrima).

Analysis of the 2C DNA

The histograms in Fig. 3 show the fluorescence distribution as a function of the number of nuclei in the examined sample (Fig. 3) The ANOVA for the 2C DNA of the species revealed a highly significant difference, with a low coefficient of variation: 1.45% ( Table 2). The Scott-Knott test clustered the taxa into six groups based on the average 2C values, with a minimum significant difference of 0.0647 (Table  1). Species from the genus Poincianella were arranged into group C. Two species remained in isolated groups: E. calycina and C. macrophyllum. The species L. ferrea and Ca. pulcherrima were placed in group E.

An ANOVA was carried out to estimate the variation of the DNA content for only the three morphotypes of Pa. echinata, thus demonstrating the existence of a significant difference at p-value = < 0.05 (Table 3). The Scott-Knott test with only the morphotypes of Pa. echinata showed a separation between the morphotypes, as the SV type was isolated in group B with a statistically significant difference in relation to the MV and LV types in group A.

Discussion

Prior studies have been carried out to locate heterochromatin rich in AT and GC among species from the Caesalpinia group using the fluorochromes CMA 3 and  [14]. In Copaifera, only CMA 3 + /DAPI − blocks were observed [24], similar to the terminal pattern observed in the present study, for which the analyzed species showed heterochromatin blocks that were rich in GC and poor in AT. These CMA 3 + /DAPI − bands have been observed only in the terminal regions of chromosomes and have ranged from two to four pairs of chromosomes in the analyzed species. Additionally, CMA 3 + blocks have been observed in the metaphase chromosomes of L. ferrea and Po. microphylla, indicating either the existence of GC-rich pericentromeric sequences in related genera (suggesting a shared trait with a common ancestor) or changes related to the composition of centromeric satellite DNA, which is qualitatively rich in GC in these species.

Fluorochrome CMA 3 − /DAPI + blocks have already been observed in Senna obtusifolia (L.) H.S. Irwin & Barneby and in one population of Chamaecrista nictitans Moench [14]. However, these bands were not reported in the present study. The absence of this type of heterochromatin may be a typical characteristic of species in the analyzed genera of Poincianella, Libidia, Erytrostemon, Paubrasilia, Caesalpinia and Cenostigma. However, the karyotype characterization on more population can relieve the interspecific variation in cytogenetic markers. In the Leguminosae family, various patterns of AT-rich and GC-poor heterochromatin were also observed for the Mimosa [25] and Erythrina [26] genera. These variations reinforce the need to characterize heterochromatin in other species so as to understand the distribution pattern and evolution of this class of DNA, which is variously colored within the Caesalpinia group.

In the present study, CMA 3 + /DAPI − terminal blocks coincided with certain 45S rDNA hybridization sites, reinforcing the fact that these CMA 3 + blocks relate to these 45S rDNA, which in turn are rich in GC bases [27][28][29]. However, the application of base-specific fluorochromes did not reveal rDNA sites with few repetitions, as the small heterochromatic block made detection and photographic documentation unviable for the epifluorescence microscope [30].

Molecular cytogenetics techniques have been widely used to localize specific in situ DNA sequences [31]. Genetically related species tend to have karyotypes with similar characteristics in terms of sequence localization and are useful in studies of plant systematics, taxonomy and evolution mainly contributing to groupings of species or cytotypes that share common characteristics, thus suggesting primitiveness or deactivation for a given cytological marker that is shared among a group of plants [32][33][34][35][36][37][38].

The hybridization sites of 45S rDNA probes in species of the Caesalpinia group have shown variations in both the number and the location of these sequences. The differences in the number of such sequences generally occur due to chromosomal rearrangements such as translocations, inversions, duplications and deletions [39], whereas variations in the signal intensity of hybridization are observed between sites with different numbers of rDNA replicates. Any changes in these sites' patterns of distribution are levels of speciation; this may assist in determining how evolution has occurred within a group of taxonomically complex plants [31,40]. Lower quantitative variation has been observed in the Poincianella genus, suggesting greater stability in the number of 45S rDNA sites (a total of eight). Conservation in the location of the rDNA genes (as revealed using FISH) was observed for species from the genus Trifolium (Leguminosae: Papilionoideae), which may indicate that some Leguminosae have great stability in this region [41]. This stability, which is based on the number and location of a chromosomal marker, is a good characteristic for species identification and delimitation through karyotype analysis.

The species in this study showed significant variations in 2C DNA and, consequently, in the size of the genomes for the evaluated species and genera. In addition, the low coefficient of variation among the replicates indicates the precision of the sample, as analyzed using flow cytometry. The variation in the amount of DNA across species can be attributed to the loss or gain of DNA sequences, which usually consist of repetitive DNA; this may occur due to evolutionary changes in accumulation and/or loss of repeating monomers in the micro and macro environments during the species' evolution [42,43]. This suggests that such losses or additions to the genome become stabilized during microevolution and selection [43].  In this work, only Ca. pulcherrima had a lower estimated amount of 2C DNA (1.63 pg) than its previous estimate (1.80 pg) [16]. This may be due to the variation in number of chromosomes for the two analyzed populations of Ca. pulcherrima, as the population evaluated in this study presented 2n = 24, but the population evaluated by OHRI et al. [16] presented 2n = 28. The estimated DNA nuclear content for diploid Ca. crista (2n = 24) indicated a too-high value of 0.707 pg per chromosome, which leads to 17.67 pg for the 2C DNA. This value, which is much higher than our result, can be attributed to the different cytophotometric methods, as Allium cepa L.'s DNA value was computed as a DNA size pattern [17]. The 2C value that we found for Ca. crista was also considerably higher than those that OHRI et al. [16] found for the species in the genus Caesalpinioideae. Analyzing the DNA content and chromosomal differences observed for taxa from the Caesalpinia genus, together with the results from the literature, requires an interdisciplinary mode in order to indicate the species' taxonomy, delimitation and clustering.

Many species in the Caesalpinia group exhibit a high degree of phenotypic plasticity, especially in foliage and leaflets. This has resulted in multiple nomenclature for the species, with each leaf-size variant having a specific condition, thus resulting in taxonomic problems [1,5,7,9]. This fact can be observed in Pa. echinata, which was previously arranged in the Caesalpinia genus and which has three morphotypes that were previously characterized using chloroplast DNA sequences [44]. The three morphotypes (leaf-size variants) presented small variations in 2C DNA, with values of 2.76, 2.81 and 2.82 pg for the SV, MV and LV types, respectively. The Scott-Knott test separated these morphotypes into two groups, one with only Pa. echinata SV and one composed of Pa. echinata MV and LV; this shows that, although the variations among Pa. echinata leaf-size morphotypes are not large, the values are sufficient to separate the SV morphotype from the other two variants, with this type's low DNA content acting as a differentiating feature.

In legume species, a positive correlation has been observed between leaf size and nuclear DNA content [45]. This relationship was also observed in this study, wherein the variants with the relatively large leaves (MV and LV) had more DNA than the SV variants. Therefore, diversification of genome size results from speciation, which, along with phenotypic changes in quantitative descriptors, is an adaptation response such as the ones observed in polyploid plants [46]. Thus, plants' 2C DNA can be used to estimate the taxonomic differentiation between species, as seen here for the variants of Pa. echinata.

The data obtained in our studies corroborate the new classification of the species that were initially placed in the Caesalpinia group [5,[9][10][11], showing that the species that had at least 1.87 pg of 2C DNA in this study should actually be grouped in the Poincianella genus. Among the species analyzed in the present study, the only representative of the Erythrostemon genus was E. calycina, which had the least amount of 2C DNA (1.54 pg) and which was also the only species to present five pairs of chromosomes with the presence of 45S rDNA sites; these were unique characteristics of this genus. In previous analyses, similarities in chromosome morphology have been reported for the six species studied herein, and karyotype formulas have shown the predominance of metacentric chromosomes [47].

Only one species was evaluated for the Cenostigma genus, C. macrophyllum; this species which was initially collected under the belief that it belonged to the Caesalpinia genus. The similarities between species of these two genera have also been visualized previously, as the species Cenostigma sclerophyllum Tul. was later described as a synonym for Caesalpinia marginata Tul. [8,48]. The analyses of C. macrophyllum enabled us to observe that it was the only species to present only six CMA 3 +/ DAPI − bands, showing that the amount of GC-rich heterochromatin was lower than that of the other species evaluated herein; this could be a feature exclusive to the Cenostigma genus.

In this work, the distribution pattern for heterochromatin, the physical location of 45S rDNA regions and the amount of DNA were all useful to corroborate studies of systematics and of evolution in Caesalpinia-group species. Although the quantity of species evaluated herein is only a small fraction of the diversity already described as belonging to the Caesalpinia group, and although some species were relocated within new genera, it was possible to observe a distinctive pattern for individual cytogenetic characteristics in the genera currently specified as Poncianella. This shows that the karyotypic analysis and the quantification of 2C DNA are valid methods to support taxonomic and biosystematics studies.

Conclusion

The quantitative variation in GC-rich heterochromatin among species in the Caesalpinia group indicates not only the variable number of satellite-related rDNA sites but also the existence of chromosomes with pericentromeric repetitive DNA with GC-rich heterochromatin in L. ferrea and Po. microphylla species. The intra-and interspecific variations in size of the GC-rich chromosomal blocks related to rDNA and satellites (relative to the location of these regions), as determined using the FISH technique with a probe for 45S rDNA. This relation suggests that the use of fluorochromes for the localization of 45S rDNA is not indicated for loci numbers identification in species of the Caesalpinia group. This fact is attributed to the minor size replications of the 45S rDNA genes in some chromosomes, which make it difficult to observe variations using CMA 3 fluorochromes. The 2C DNA may not be even related to morphological leaf size for Paubrasilia echinata, so this relationship must be evaluated again with other populations to increase the number of analyzed plants. On the other hand, the 2C DNA helped us to see all the variation in the Caesalpinia group within this trait. All the data indicate that the actual taxonomy appropriated by the Caesalpinia group is due to the larger chromosomal and genome-size variations, which could be clustered in a specific genus. This information shows the group from point of view that differs from the old taxonomy, which grouped all species analyzed herein as just one genus.

Fig. 1
1Application of fluorochromes in the Caesalpinia group. The fluorochromes DAPI (a, d, g, j) and CMA 3 (b, e, h, k), as well as the FISH (c, g, i, l) with a probe for 45S rDNA on metaphase chromosomes. ac Cynostigma macrophyllum, (d -f) Erythrostemon calycina, (g -i) Poincianella pluviosa and (j -l) Libidibia ferrea. White arrows indicate CMA 3 + /DAPI − blocks, orange arrows indicate CMA 3 + blocks, and red arrows indicate 45S rDNA sites; a bar is 10 μm Location of 45S rDNA sites

Fig. 2
2Application of fluorochromes in the Caesalpinia group. The fluorochromes DAPI (a, d, g, j, m, p) and CMA 3 (b, e, h, k, n, q), as well as the FISH (c, g, i, l, o, r) with a probe for 45S rDNA on metaphase chromosomes. ac Poincianella bracteosa, (d -f) Po. laxiflora, (g -i) Po. microphylla, (j -l) Po. pyramidalis (m -o), Caesalpinia pulcherrima and (p -r) Paubrasilia echinata (SV). White arrows indicate CMA 3 + /DAPI − blocks, orange arrows indicate CMA 3 + blocks, and red arrows indicate 45S rDNA sites

Fig. 3
3Histograms with 2C DNA content for species of the Caesalpinia group. Internal standards used: Glycine max (L.) Merr and Zea mays L DAPI. Previous studies in the Senna and Chamaecrista genera revealed the presence of

2C DNA: 2C DNA content; 2n: Diploid number; ANOVA: Analysis of variance; CMA 3 : Chromomycin A 3 ; DAPI: 4′-6-diamidino-2-phenylindole; FISH: Fluorescent in situ hybridization; LV: Large-leaf variant; MV: Medium-leaf variant; rDNA: Ribosomal DNA; SSC: Salt sodium citrate; SV: Small-leaf variant

Table 1
1Estimates of nuclear genome size (2C DNA content) for Caesalpinia groupSpecies

). In general, the terminal heterocyclic blocks were of distinct sizes. CMA 3 + / DAPI − terminal blocks were observed in two chromosome pairs of Cenostigma macrophyllum, Po. pluviosa, Caesalpinia pulcherrima and Pa. echinata. Three chromosome pairs with CMA 3 + /DAPI − terminal blocks were observed in Po. bracteosa, Po.laxiflora, Po. microphylla and L. ferrea. DAPI − terminal blocks was observed in Po. pyramidalis, which had four chromosome pairs with GC-rich heterochromatic blocks.Erythrostemon calycina had seven chromosomes with
CMA 3
+ /DAPI − terminal blocks. The highest number of
CMA 3
+ /

. The Pa. echinata large-leaf variant (LV) was the taxon with the highest 2C value of DNA (2.82 pg). The Pa. echinata LV was then used as the reference for high DNA content; compared to this reference value, the DNA content of the other species was lower by 46.3% for L. ferrea, 45.4% for E. calycina, 42.2% for Ca. pulcherrima, 35.1% for C. macrophyllum, 33.7% for Po. microphylla, 33.3% for Po. pluviosa, 32.6% % for Po. laxiflora, 31.9% for Po. bracteosa and Po. pyramidalis, 2.1% for the Pa. echinata small-leaf variant (SV), and 0.4% for the Pa. echinata medium-leaf variant (MV).

Table 2
2Summary of the ANOVA for the characteristic 2C DNA content among the analyzed speciesSources of variation
GL
QM

Taxon
11
1.136674 a

Error
48
0.000888

CV (%)
1.46

a

significant at 1% probability by F test
GL Level of freedom, QM Average square, CV Coefficient of variation

Table 3
3Summary of the ANOVA for 2C DNA content among the three morphological leaf variants of Paubrasilia echinataSources of variation
GL
QM

Taxon
2
0.012653 a

Error
12
0.012640

CV (%)
1.16

a

significant at 5% probability by F test
GL Level of freedom, QM Average square, CV Coefficient of variation

AcknowledgmentsThe authors would like to thank José Lima Paixão, a technician who assisted in the collection and identification of the botanical material, as well as the collaborators at the herbarium of Universidade Estadual de Feira de Santana for confirming the identification of the Cenostigma species.FundingThis research received financial support from Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) (Grant number 473393/ 2007-7) and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) for the collection of germoplasm and for the acquisition of laboratory supplies for the analysis. The cytogenetics techniques were performed at the Plant Breeding Laboratory at Universidade Estadual de Santa Cruz, which provided the financial support for the equipment and physical structures used in this study.Availability of data and materialsAll data sets that support this article's conclusions are included in the article.Authors' contributions PSR performed the cytogenetic laboratory procedures and prepared the manuscript text. RXC helped write and review the text and made several contributions on the relations between plants, as well as obtaining financial support. CAFM helped with the cytogenetic procedures, analyzed the results and reviewed the text. TNSP helped with the flow cytometry procedure and the data interpretation. MMS participated in the molecular cytogenetic analysis and helped with the data interpretation and the protocol adjustment. All authors have read and approved the manuscript and consent to its publication.Ethics approval and consent to participate Permission for seeds collection was not necessary because the collections were made in public areas. We also found some species near roads and in rural areas. Taxonomists at the herbarium museums of Universidade Estadual de Santa Cruz (Prof. Dr. André Márcio Araújo Amorim) and Universidade Estadual de Feira de Santana (Prof. Dr. Luciano Paganucci de Queiroz) certified all species in Bahia, Brazil.Consent for publicationNot applicable.Competing interestsThe authors declare that they have no competing interests.Publisher's NoteSpringer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Author details• We accept pre-submission inquiries • Our selector tool helps you to find the most relevant journal Submit your next manuscript to BioMed Central and we will help you at every step:
Legumes of the world. G P Lewis, B Schrire, B Mackinder, M Lock, 577Kew: Roy Bot GardLewis GP, Schrire B, Mackinder B, Lock M. Legumes of the world. Kew: Roy Bot Gard; 2005. p. 577.

Structural evolution in the Caesalpinioideae (Leguminosae). P S Herendeen, Advances in Legume Systematics. Roy Bot Gardpart 9Herendeen PS. Structural evolution in the Caesalpinioideae (Leguminosae). In: Herendeen PS, Bruneau A, editors. Advances in Legume Systematics, part 9. Kew: Roy Bot Gard, 2000. p. 45-64.

Sistemática de angiospermas do Brasil. G M Barroso, A L Peixoto, G C Costa, Clf Ichasso, E F Guimarães, H C Lima, 337Viçosa: UFVBarroso GM, Peixoto AL, Costa GC, Ichasso CLF, Guimarães EF, Lima HC. Sistemática de angiospermas do Brasil, vol. 337. Viçosa: UFV; 1984.

Abordagens taxonômicas e citogenéticas em Leguminosae-Caesalpinioideae na região sul do Brasil. Rlc Bortoluzzi, E Biondo, Sts Miotto, M T Schifino-Witmann, Revista Brasileira de Biociências. 5Bortoluzzi RLC, Biondo E, Miotto STS, Schifino-Witmann MT. Abordagens taxonômicas e citogenéticas em Leguminosae-Caesalpinioideae na região sul do Brasil. Revista Brasileira de Biociências. 2007;5:339-41.

A new generic system for the pantropical Caesalpinia group (Leguminosae). E Gagnon, A Bruneau, C E Hughes, L P Queiroz, G P Lewis, PhytoKeys. 71Gagnon E, Bruneau A, Hughes CE, Queiroz LP, Lewis GP. A new generic system for the pantropical Caesalpinia group (Leguminosae). PhytoKeys. 2016;71:1-160.

IUCN red list of threatened species. Iucn, IUCN. IUCN red list of threatened species; 2013.

R Polhill, J Vidal, Caesalpinieae, Advances in legume systematics. 1. Richmond: Roy Bot GardPolhill R, Vidal J. Caesalpinieae. In: Polhill R, Raven PH, editors. Advances in legume systematics, part 1. Richmond: Roy Bot Gard; 1981. p. 81-95.

Caesalpinia, a revision of the Poincianella -Erythrostemon group. G P Lewis, Kew: Roy Bot GardLewis GP. Caesalpinia, a revision of the Poincianella -Erythrostemon group. Kew: Roy Bot Gard; 1998.

Molecular phylogeny of Caesalpinia sensu lato: increased sampling reveals new insights and more genera than expected. E Gagnon, G P Lewis, J S Sotuyo, C E Hughes, A A Bruneau, S Afr J Bot. 89Gagnon E, Lewis GP, Sotuyo JS, Hughes CE, Bruneau AA. Molecular phylogeny of Caesalpinia sensu lato: increased sampling reveals new insights and more genera than expected. S Afr J Bot. 2013;89:111-27.

. L P Queiroz, Leguminosas Da Caatinga, Santana Feira De, Universidade Estadual de Feira de Santana. Associação Plantas do NordesteQueiroz LP. Leguminosas da Caatinga. Feira de Santana: Universidade Estadual de Feira de Santana. Associação Plantas do Nordeste; 2009.

New combinations in Libidibia (DC.) Schltdl. And Poincianella Britton & Rose (Leguminosae, Caesalpinioideae). L P Queiroz, Neodiversity. 5Queiroz LP. New combinations in Libidibia (DC.) Schltdl. And Poincianella Britton & Rose (Leguminosae, Caesalpinioideae). Neodiversity. 2010;5:11-2.

Karyotype analysis for diploid and polyploid species of the Solanum L. Caf Melo, Mig Martins, Mbm Oliveira, Amr Benko-Iseppon, R Carvalho, Plant Syst Evol. 293Melo CAF, Martins MIG, Oliveira MBM, Benko-Iseppon AMR, Carvalho R. Karyotype analysis for diploid and polyploid species of the Solanum L. Plant Syst Evol. 2011;293:227-35.

Applications of flow cytometry to evolutionary and population biology. P Kron, J Suda, B C Husband, Annu Rev Ecol Evol Syst. 38Kron P, Suda J, Husband BC. Applications of flow cytometry to evolutionary and population biology. Annu Rev Ecol Evol Syst. 2007;38:847-76.

Benko-Iseppon AM. Cytogenetics and chromosome banding patterns in Caesalpinioideae and Papilionioideae species of Pará. M D Soares-Scott, Lmm Meletti, L C Bernacci, Irs Passos, Faleiro FG, Souza MGC. Amazonas, Brazil144Citogenética clássica e molecular em passiflorasSoares-Scott MD, Meletti LMM, Bernacci LC, Passos IRS. Citogenética clássica e molecular em passifloras. In: Faleiro FG, Souza MGC, Benko- Iseppon AM. Cytogenetics and chromosome banding patterns in Caesalpinioideae and Papilionioideae species of Pará, Amazonas, Brazil. Bot J Linn Soc. 2004;144:181-91.

Variability in the pod and seed parameters and nuclear DNA content of Tunisian populations of Ceratonia siliqua L. Efo Héla, S Naghmouchi, D J Walker, E Correal, M Boussaïd, M L Khouja, Caryologia. 614Héla EFO, Naghmouchi S, Walker DJ, Correal E, Boussaïd M, Khouja ML. Variability in the pod and seed parameters and nuclear DNA content of Tunisian populations of Ceratonia siliqua L. Caryologia. 2008;61(4):354-62.

Correlations between 2C DNA values and habit in Cassia (Leguminosae: Caesalpinioideae). D Ohri, A Kumar, M Pal, Plant Syst Evol. 153Ohri D, Kumar A, Pal M. Correlations between 2C DNA values and habit in Cassia (Leguminosae: Caesalpinioideae). Plant Syst Evol. 1986;153:223-7.

Identification of RAPD markers, in situ DNA content and structural chromosomal diversity in some legumes of the mangrove flora of Orissa. S Jena, P Sahoo, S Mohanty, A B Das, Genetica. 122Jena S, Sahoo P, Mohanty S, Das AB. Identification of RAPD markers, in situ DNA content and structural chromosomal diversity in some legumes of the mangrove flora of Orissa. Genetica. 2004;122:217-26.

. D A Johansen, Plant Microtechique, McGraw-Hill Book CompanyNew YorkJohansen DA. Plant Microtechique. New York: McGraw-Hill Book Company; 1940.

Como observar cromossomos: Um guia de práticas em citogenética vegetal, animal e humana. Ribeirão Preto: Editora Funpec. M Guerra, M J Souza, Guerra M, Souza MJ. Como observar cromossomos: Um guia de práticas em citogenética vegetal, animal e humana. Ribeirão Preto: Editora Funpec; 2002.

Improvements in cytological preparations for fluorescent in situ hybridization in Passiflora. M M Souza, J D Urdampilleta, E R Forni-Martins, Genet Mol Res. 94Souza MM, Urdampilleta JD, Forni-Martins ER. Improvements in cytological preparations for fluorescent in situ hybridization in Passiflora. Genet Mol Res. 2010;9(4):2148-55.

Plant genome size estimation by flow cytometry: inter-laboratory comparison. J Dolezel, J Greilhuber, S Lucretti, A Meister, Kma Lysa, L Nardi, R Obermayer, Ann Bot. 82Dolezel J, Greilhuber J, Lucretti S, Meister A, Lysa KMA, Nardi L, Obermayer R. Plant genome size estimation by flow cytometry: inter-laboratory comparison. Ann Bot. 1998;82:17-26.

Sex determination in dioecious plants Melandrium album and M. Rubrum using high-resolution flow cytometry. J Dolezel, W Gohde, Cytometry. 19Dolezel J, Gohde W. Sex determination in dioecious plants Melandrium album and M. Rubrum using high-resolution flow cytometry. Cytometry. 1995;19:103-6.

Programa Sisvar. Software 5.0: UFLA. D F Ferreira, Ferreira DF. Programa Sisvar. Software 5.0: UFLA; 2003.

Occurrence and chromosome distribution of retroelements and NUPT sequences in Copaifera langsdorffii Desf. (Caesalpinioideae). M L Gaeta, P M Yuyama, D Sartori, Mhp F, All V, Chromosom Res. 18Gaeta ML, Yuyama PM, Sartori D, MHP F, ALL V. Occurrence and chromosome distribution of retroelements and NUPT sequences in Copaifera langsdorffii Desf. (Caesalpinioideae). Chromosom Res. 2010;18:515-24.

Polyploidy, B chromosomes, and heterochromatin characterization of Mimosa caesalpiniifolia Benth. S M Sousa, A C Reis, Viccini, Fabaceae-MimosoideaeSousa SM, Reis AC. Viccini. Polyploidy, B chromosomes, and heterochromatin characterization of Mimosa caesalpiniifolia Benth. (Fabaceae-Mimosoideae).

. Tree Genet Genomes. 9Tree Genet Genomes. 2013;9:613-9.

Karyological features and banding patterns in Arachis species belonging to the Heteranthae section. S C Silva, Mig Martins, R C Santos, Aps Peñaloza, P A Melo-Filho, A Benko-Iseppon, Jfm Valls, R Carvalho, Plant Syst Evol. 285Silva SC, Martins MIG, Santos RC, Peñaloza APS, Melo-Filho PA, Benko-Iseppon A, Valls JFM, Carvalho R. Karyological features and banding patterns in Arachis species belonging to the Heteranthae section. Plant Syst Evol. 2010;285:201-7.

Fluorescent chromosome banding in plants: application, mechanisms, and implications for chromosome structure. D Schweizer, Proceedings of the fourth John Innes symposium. the fourth John Innes symposiumSchweizer D. Fluorescent chromosome banding in plants: application, mechanisms, and implications for chromosome structure, Proceedings of the fourth John Innes symposium, held September; 1979. p. 61-72.

Chromosome banding. London: Unwin and Hyman. A T Sumner, Sumner AT. Chromosome banding. London: Unwin and Hyman; 1990.

Species relationships among the wild B genome of Arachis species (section Arachis) based on FISH mapping of rDNA loci and heterochromatin detection: a new proposal for genome arrangement. G Robledo, G Seijo, Theor Appl Genet. 1216Robledo G, Seijo G. Species relationships among the wild B genome of Arachis species (section Arachis) based on FISH mapping of rDNA loci and heterochromatin detection: a new proposal for genome arrangement. Theor Appl Genet. 2010;121(6):1033-46.

Karyomorphology and GC-rich heterochromatin pattern in Passiflora (Passifloraceae) wild species from Decaloba and Passiflora subgenera. Caf Melo, M M Souza, P P Abreu, Ajc Viana, Flora. 209Melo CAF, Souza MM, Abreu PP, Viana AJC. Karyomorphology and GC-rich heterochromatin pattern in Passiflora (Passifloraceae) wild species from Decaloba and Passiflora subgenera. Flora. 2014;209:620-31.

DNA : DNA in situ hybridization and methods for light microscopy. T Schwarzacher, A R Leitch, J S Heslop-Harrison, Plant cell biology: a practical approach. Harris N, Oparka KJOxfordOxford University PressSchwarzacher T, Leitch AR, Heslop-Harrison JS. DNA : DNA in situ hybridization and methods for light microscopy. In: Harris N, Oparka KJ, editors. Plant cell biology: a practical approach. Oxford: Oxford University Press; 1994. p. 127-55.

Physical mapping of rDNA loci in Brassica species. J Maluszynska, J S Heslop-Harrison, Genome. 36Maluszynska J, Heslop-Harrison JS. Physical mapping of rDNA loci in Brassica species. Genome. 1993;36:774-81.

Chromatin characterization by banding techniques, in situ hybridization, and nuclear DNA content in Cicer L. (Leguminosae). I Galasso, D Pignone, M Frediani, M Maggiani, Genome. 39Galasso I, Pignone D, Frediani M, Maggiani M. Chromatin characterization by banding techniques, in situ hybridization, and nuclear DNA content in Cicer L. (Leguminosae). Genome. 1996;39:258-65.

Comparison of ribosomal DNA sites in Lolium species by fluorescente in situ hibridization. H M Thomas, J Á Harper, M R Meredith, W G Morgan, I D Thomas, E Timms, I P King, Chromosom Res. 4Thomas HM, Harper JÁ, Meredith MR, Morgan WG, Thomas ID, Timms E, King IP. Comparison of ribosomal DNA sites in Lolium species by fluorescente in situ hibridization. Chromosom Res. 1996;4:486-90.

Physical mapping of ribosomal RNA genes in peonies (Paeonia, Paeoniaceae) by fluorescent in situ hybridization and concert evolution. D Zhang, T Sang, Am J Bot. 85Zhang D, Sang T. Physical mapping of ribosomal RNA genes in peonies (Paeonia, Paeoniaceae) by fluorescent in situ hybridization and concert evolution. Am J Bot. 1999;85:735-40.

Karyotype evolution in the genus Boronia (Rutaceae). J G Seijo, G I Lavia, A Fernández, A Krapovickas, D Ducasse, F Shan, G Yan, J A Plummer, Bot J Linn Soc. 142Seijo JG, Lavia GI, Fernández A, Krapovickas A, Ducasse D, Shan F, Yan G, Plummer JA. Karyotype evolution in the genus Boronia (Rutaceae). Bot J Linn Soc. 2003;142:309-20.

Chromosome evaluation of southern south American species of Camptosema and allied genera (Diocleinae-Phaseolae-Papilionoideae-Leguminosae). S M Sede, R H Fortunato, L Poggio, Bot J Linn Soc. 52Sede SM, Fortunato RH, Poggio L. Chromosome evaluation of southern south American species of Camptosema and allied genera (Diocleinae- Phaseolae-Papilionoideae-Leguminosae). Bot J Linn Soc. 2006;52:235-43.

Characterization of the Arachis (Leguminosae) D genome using fluorescence in situ hybridization (FISH) chromosome markers and total genome DNA hybridization. G Robledo, G Seijo, Genet Mol Biol. 313Robledo G, Seijo G. Characterization of the Arachis (Leguminosae) D genome using fluorescence in situ hybridization (FISH) chromosome markers and total genome DNA hybridization. Genet Mol Biol. 2008;31(3):717-24.

Quantitative karyotyping and dual-color FISH mapping of 5S and 18S-25S rDNA probes in the cultivated Phaseolus species (Leguminosae). E A Moscone, F Klein, M Lambrou, J Fuchs, D Schweizer, Moscone EA, Klein F, Lambrou M, Fuchs J, Schweizer D. Quantitative karyotyping and dual-color FISH mapping of 5S and 18S-25S rDNA probes in the cultivated Phaseolus species (Leguminosae).

. Genome. 426Genome. 1999;42(6):1224-33.

Comparative genome organization in plants: from sequence and markers to chromatin and chromosomes. Js H-H, Plant Cell. 12JS H-H. Comparative genome organization in plants: from sequence and markers to chromatin and chromosomes. Plant Cell. 2000;12:617-35.

Molecular cytogenetic of 5S and 18S-26S rDNA loci in white clover (Trifolium repens L.) and related species. H A Ansari, N W Ellison, S M Reader, E D Badaevas, B Friebe, T E Miller, W M Williams, Ann Bot. 83Ansari HA, Ellison NW, Reader SM, Badaevas ED, Friebe B, Miller TE, Williams WM. Molecular cytogenetic of 5S and 18S-26S rDNA loci in white clover (Trifolium repens L.) and related species. Ann Bot. 1999;83:199-206.

Nuclear DNA content variation within angiosperm species. H J Price, Evolution Trends. 12Price HJ. Nuclear DNA content variation within angiosperm species. Evolution Trends. 1998;1(2):53-60.

Interespecific genetic in 15 species of Cassia L. evident by chromosome and 4C nuclear DNA analysis. S Mohanty, A B Das, J Biolo Sci. 64Mohanty S, Das AB. Interespecific genetic in 15 species of Cassia L. evident by chromosome and 4C nuclear DNA analysis. J Biolo Sci. 2006;6(4):664-70.

Phylogenetic relationships among morphotypes of Caesalpinia echinata lam. (Caesalpinioideae: Leguminosae) evidenced by trnL intron sequences. F S Juchum, Naturwissenschaften. 95Juchum FS. Phylogenetic relationships among morphotypes of Caesalpinia echinata lam. (Caesalpinioideae: Leguminosae) evidenced by trnL intron sequences. Naturwissenschaften. 2008;95:1085-91.

Relationships between nuclear DNA content and seed and leaf size in soybean. Ljh Chungj, K Arumuganathan, G L Graef, J E Specht, Theo Appl Genet. 96ChungJ LJH, Arumuganathan K, Graef GL, Specht JE. Relationships between nuclear DNA content and seed and leaf size in soybean. Theo Appl Genet. 1998;96:1064-8.

Polyploidy and cellular mechanisms changing leaf size: comparison of diploid and autotetraploid populations in two species of Lolium. S-I Sugiyama, Ann Bot. 96Sugiyama S-I. Polyploidy and cellular mechanisms changing leaf size: comparison of diploid and autotetraploid populations in two species of Lolium. Ann Bot. 2005;96:931-8.

Karyomorphology and karyotype asymmetry in the south American Caesalpinia species (Leguminosae and Caesalpinioideae). P S Rodrigues, M M Souza, R X Correa, Genet Mol Biol. 13Rodrigues PS, Souza MM, Correa RX. Karyomorphology and karyotype asymmetry in the south American Caesalpinia species (Leguminosae and Caesalpinioideae). Genet Mol Biol. 2014;13:8278-93.

A revision of Cenostigma (Leguminosae -Caesalpinioideae -Caesalpinieae), a genus endemic to Brazil. M C Warwick, G P Lewis, Kew Bull. 64Warwick MC, Lewis GP. A revision of Cenostigma (Leguminosae - Caesalpinioideae -Caesalpinieae), a genus endemic to Brazil. Kew Bull. 2009;64:135-46.


Molecular Characterization of Rotavirus Gastroenteritis Strains, Iraqi Kurdistan

Herish M Ahmed
Liverpool School of Tropical Medicine
†University of Liverpool
Liverpool, LiverpoolUK;, UK

‡Erbil Paediatric Hospital
Erbil

Iraqi Kurdistan; §Nagasaki University
NagasakiJapan

‡ J Brian
S Coulter
Liverpool School of Tropical Medicine
†University of Liverpool
Liverpool, LiverpoolUK;, UK

‡Erbil Paediatric Hospital
Erbil

Iraqi Kurdistan; §Nagasaki University
NagasakiJapan

Osamu Nakagomi
C A Hart
Jamal M Zaki
Abas A Al-Rabaty
Winifred Dove
Nigel A Cunliffe
Molecular Characterization of Rotavirus Gastroenteritis Strains, Iraqi Kurdistan

Of 260 children with acute diarrhea in Erbil, Iraqi Kurdistan, 96 (37%) were infected with rotavirus. Reverse transcription-polymerase chain reaction identified G1, G4, G2, G9, P[8], P[6], and P[4] as the most common genotypes. Eight G/P combinations were found, but P[8]G1 and P[4]G2 accounted for >50% of the strains. DISPATCHES 824 Emerging Infectious Diseases • www.cdc.gov/eid • Vol. 12, No. 5,

R otavirus is the single most important cause of severe gastroenteritis in young children throughout the world. Globally, an estimated 702,000 children die each year due to rotavirus diarrhea (1). This large impact of rotavirus disease has speeded the development of rotavirus vaccines, and 2 live, attenuated rotavirus vaccines are expected to be available for global use within the next few years (1). Therefore, determining the prevalence and types of rotaviruses within regions is essential to prepare for introducing a vaccine.

Rotavirus, a member of the family Reoviridae, has a triple-layered capsid that contains 11 segments of doublestranded genomic RNA. While protective immunity against rotavirus infection is not completely understood, serotype-specific immunity is believed to play a major role (1). Rotavirus serotypes are defined by genome segment 4 for the P (protease-sensitive protein) type and by genome segment 9 (or 7 or 8, depending on the strain) for the G (glycoprotein) type. Fourteen G types exist, of which G1-G4 are commonly found in children with diarrhea, but a recent increase in the detection of serotype G8 and G9 strains has captured considerable attention (2)(3)(4). While >24 P types have been reported in the literature, only P [4], P [6], and P [8] are commonly found among human rotaviruses (1)(2)(3).

In Iraq, the death rate in children <5 years of age was reported to be 130/1,000 for boys and 120/1,000 for girls in 2003 (5). Diarrhea is a major cause of illness and death in Iraqi children; however, little information exists about the origin of childhood diarrhea. Only a single study showed that rotavirus accounted for 24% of acute diarrhea in hospitalized children in Basrah (6).

The Study

This descriptive, cross-sectional study of 6 weeks' duration was undertaken at Erbil Paediatric Hospital in Iraqi Kurdistan between March and May 2005. The study recruited 260 children from 1 month to 5 years of age who were admitted with acute diarrhea (defined as the passage of watery or loose stools >3 times per day for <2 weeks' duration). Basic demographic, epidemiologic, and clinical information were collected prospectively, according to a pro forma. Ethical approval for the research was obtained from the review boards of the Liverpool School of Tropical Medicine and Erbil Paediatric Hospital. The hospital serves a population of ≈1.5 million, and ≈3,116 births per month occur in this population.

A commercial enzyme-linked immunosorbent assay (ELISA) was used to detect rotavirus antigen (Rotaclone, Meridian Diagnostics, Cincinnati, OH, USA). Stool samples were then stored frozen in the laboratory of the study hospital until they were transported to Liverpool for rotavirus genotyping and electropherotyping. All samples (66) with an absorbance equal to or greater than the positive control for the ELISA were subjected to genotyping. Rotavirus genomic RNA was extracted with guanidine isothiocyanate, followed by adsorption to and elution from silica particles according to the method described by Gentsch et al. (7). The purified RNA was then used to determine the P type and G type of rotavirus present in the stool specimens by reverse transcription-polymerase chain reaction as described by Gentsch et al. (7) and by Gouvea et al. (8). Rotavirus electropherotypes were determined by polyacrylamide gel electrophoresis according to the method described by Koshimura et al. (9), with some modifications.

Of 260 stool specimens tested by ELISA, 96 (37%) were positive for rotavirus. Rotavirus-positive patients had a mean age (SD) of 9.3 (8.5) months compared to 11.1 (10.1) months in the rotavirus-negative patients. These results suggest that rotavirus positive cases were slightly younger, although the difference was not statistically significant (p = 0.14). Rotavirus-positive patients were similar to rotavirus-negative patients in most of the epidemiologic and clinical characteristics (data not shown). However, rotavirus-positive patients were more likely to exhibit vomiting and have a shorter duration of diarrhea (p<0.01 for both analyses).

Of the 66 rotavirus strains that underwent molecular characterization, 25 (38%) were G1, 11 (17%) were G2, 13 (20%) were G4, and 7 (11%) were G9. Four (6%) were mixed infections (3 G1/G2, 1 G2/G4), and 6 (9%) were G nontypeable. A total of 7 (11%) were P [4], 10 (15%) were P [6], and 45 (68%) were P [8]. One showed mixed P [4] and P [8] genotypes (mixed with G1/G2), and 3 (5%) were P nontypeable. None of the rotaviruses was both G and P nontypeable.

A total of 8 different P and G genotype combinations were detected (Table). The most common combinations were P An electropherotype was obtained for 50 of the 66 genotyped strains. Of these, 11 (22%) had a short electropherotype, and 39 (78%) had a long electropherotype (Table). Most of the short electropherotypes were the expected G2 strains; however, 1 strain (P[8]G9) also had a short electropherotype.

Conclusions

The only other study of viral gastroenteritis from Iraq (Basrah in the south) demonstrated that 24% of children with acute gastroenteritis were infected with rotavirus (6). This figure is somewhat lower than the 37% detection rate in our study. Moreover, the prevalence we found is similar to those reported from neighboring countries such as Iran (35%) (10), Jordan (33%) (11), Kuwait (40%) (12), and Turkey (37%) (13). However, our study was undertaken over a 6-week period from the end of March to the beginning of May 2005. No information is available on the seasonal prevalence of rotavirus infection in Iraq, and a longer study is warranted to determine the true prevalence of rotavirus infection and its seasonality in northern Iraq. However, the peaks of rotavirus infection in Iran, Kuwait, and Turkey were February-March, March-May, and December, respectively (10,12,13). More than 75% of our cases of rotavirus diarrhea occurred in children <1 year of age, with an overall mean age of slightly more than 9 months. This pattern is similar to that in many developing countries. In Jordan the mean age of children with rotavirus diarrhea was 7.2 months (10). However, in other countries in the region the distribution was different; 30% of the infants with rotavirus in Iran were <1 year of age (10,12), 50% in Kuwait were <1 year of age, and 63% in Turkey were <2 years of age (13).

Although this study period was brief, we detected a variety of rotavirus strains. Four of the major global human rotavirus genotypes (G1, G2, G4, G9) were detected, as were each of the major P genotypes (P [4], P [6], P [8]). In Iran, in a study undertaken in 2001 and 2002, only G1 and G2 rotaviruses were detected, and the only P types were P [4] and P [8] (10), and in Turkey over a 2-year period (2000)(2001)(2002), G types G1-G4 and G9, as well as each of the 3 major human P types were found (14). In Iraq, the combinations P[8]G1 and P [8]G4 accounted for >50% of the strains of rotavirus. In Iran, P [8]G1 accounted for 95% of the strains, but P[8]G4 was not detected (10). In Turkey, P[8]G4 (42%) and P[8]G1 (27%) accounted for more than two thirds of the strains (14). G3 rotaviruses were not detected in Iraq or Iran, and in Turkey only 1 of the 65 strains was of genotype G3. Genotype G9 was detected in 13% of the Iraqi strains, a similar finding to results in Turkey (14). We also detected mixed rotavirus infections in 6% of our patients, again similar to the findings in Turkey (14). The presence of mixed rotavirus infections indicates that new rotavirus strains may evolve by reassortment (1)(2)(3).

Finally, among the G9 strains, one P[6]G9 had a long electropherotype, and one P [8]G9 had a short electropherotype. The P[6]G9 and P [8]G9 strains were both cultured and subgrouped by ELISA with monoclonal antibodies and found to be of subgroup II. Partial sequences (831 bp) were obtained for their VP7 genes (AB247941 and AB247943; available from the DNA Data Bank of Japan: www.ddbj.nig.ac.jp). They showed 99.4% similarity to each other and >99% similarity to strains from Australia (AY307087), Belgium (AY487858, AY487856), and India (RG9491165). A strain similar to our P[6]G9, called variant 3, was first detected in India, and strains similar to our P [8]G9, called variant 2, have been described in Bangladesh and in the United States (15).

Although the major global genotypes (except for G3 strains) were detected, clearly, rotavirus strains are continuing to diversify in Iraq and other parts of the region. This circumstance may pose challenges to the efficacy of rotavirus vaccines.
 Rotavirus Gastroenteritis Strains, Iraqi Kurdistan Emerging Infectious Diseases • www.cdc.gov/eid • Vol. 12, No. 5, May 2006
AcknowledgmentsWe thank the Ministry of Health of the Kurdistan Regional Government in Iraq for giving permission to undertake the study and the management and staff of Erbil Paediatric Hospital for their cooperation and support.Dr Ahmed is a pediatrician who conducted this research as part of the requirements for his master's degree in tropical pediatrics. His research interests are in viral gastroenteritis and respiratory tract infections.
A critical time for rotavirus vaccines. N A Cunliffe, O Nakagomi, Expert Rev Vaccines. 4Cunliffe NA, Nakagomi O. A critical time for rotavirus vaccines. Expert Rev Vaccines. 2005;4:521-34.

Global distribution of rotavirus serotypes/genotypes and its implication for the development and implementation of an effective rotavirus vaccine. N Santos, Y Hoshino, Rev Med Virol. 15Santos N, Hoshino Y. Global distribution of rotavirus serotypes/geno- types and its implication for the development and implementation of an effective rotavirus vaccine. Rev Med Virol. 2005;15:29-56.

Serotype diversity and reassortment between human and animal rotavirus strains: implications for rotavirus vaccine programs. J R Gentsch, A R Laird, B Bielfelt, D D Griffin, K Banyai, M Ramachandran, J Infect Dis. 1921SupplGentsch JR, Laird AR, Bielfelt B, Griffin DD, Banyai K, Ramachandran M, et al. Serotype diversity and reassortment between human and animal rotavirus strains: implications for rotavirus vac- cine programs. J Infect Dis. 2005;192(Suppl 1):S146-59.

Expanding global distribution of rotavirus G9 serotype: detection in Libya. N A Cunliffe, W Dove, J E Bunn, M Ben Ramadan, J W Nyangao, R L Riveron, Emerg Infect Dis. 7Cunliffe NA, Dove W, Bunn JE, Ben Ramadan M, Nyangao JW, Riveron RL, et al. Expanding global distribution of rotavirus G9 serotype: detection in Libya, Kenya, and Cuba. Emerg Infect Dis. 2001;7:890-2.

World Health Organization-Iraq. Statistics by country or region. World Health Organization-Iraq. Statistics by country or region. 2003 [cited 2005 Oct 19]. Available from http://www3.who.int/who- sis/country/indicators.cfm?country=IRQ&language=English.

Clinical and epidemiological characteristics of rotavirus-and EPEC-associated hospitalized infantile diarrhoea in Basrah. D A Mahmood, R G Feachem, Iraq. J Trop Pediatr. 33Mahmood DA, Feachem RG. Clinical and epidemiological character- istics of rotavirus-and EPEC-associated hospitalized infantile diar- rhoea in Basrah, Iraq. J Trop Pediatr. 1987;33:319-25.

Identification of group A rotavirus gene 4 types by polymerase chain reaction. J R Gentsch, R I Glass, P Woods, V Gouvea, M Gorziglia, J Flores, J Clin Microbiol. 30Gentsch JR, Glass RI, Woods P, Gouvea V, Gorziglia M, Flores J, et al. Identification of group A rotavirus gene 4 types by polymerase chain reaction. J Clin Microbiol. 1992;30:1365-73.

Polymerase chain reaction amplification and typing of rotavirus nucleic acid from stool specimens. V Gouvea, R I Glass, P Woods, K Taniguchi, H F Clark, B Forrester, J Clin Microbiol. 28Gouvea V, Glass RI, Woods P, Taniguchi K, Clark HF, Forrester B, et al. Polymerase chain reaction amplification and typing of rotavirus nucleic acid from stool specimens. J Clin Microbiol. 1990;28: 276-82.

The relative frequencies of G serotypes of rotaviruses recovered from hospitalized children with diarrhoea: a 10 year survey (1987-1996) in Japan with a review of globally collected data. Y Koshimura, T Nakagomi, O Nakagomi, Microbiol Immunol. 44Koshimura Y, Nakagomi T, Nakagomi O. The relative frequencies of G serotypes of rotaviruses recovered from hospitalized children with diarrhoea: a 10 year survey (1987-1996) in Japan with a review of globally collected data. Microbiol Immunol. 2000;44:499-510.

Epidemiology of rotavirus diarrhoea in Iranian children. B Khalili, L E Cuevas, N Reisi, W Dove, N A Cunliffe, C A Hart, J Med Virol. 73Khalili B, Cuevas LE, Reisi N, Dove W, Cunliffe NA, Hart CA. Epidemiology of rotavirus diarrhoea in Iranian children. J Med Virol. 2004;73:309-12.

Bacterial, viral and parasitic enteric pathogens associated with acute diarrhea in hospitalized children from northern Jordan. M Youssef, A Shurman, M-E Bougnoux, M Rawashdeh, S Bretagne, N Strockbine, FEMS Immunol Med Microbiol. 28Youssef M, Shurman A, Bougnoux M-E, Rawashdeh M, Bretagne S, Strockbine N. Bacterial, viral and parasitic enteric pathogens associ- ated with acute diarrhea in hospitalized children from northern Jordan. FEMS Immunol Med Microbiol. 2000;28:257-63.

Acute diarrhea and rotavirus infections in young children in Kuwait. S K Sethi, W Al-Nakib, F A Khuffash, H A Majeed, Ann Trop Paediatr. 4Sethi SK, Al-Nakib W, Khuffash FA, Majeed HA. Acute diarrhea and rotavirus infections in young children in Kuwait. Ann Trop Paediatr. 1984;4:117-21.

Childhood diarrhoea in Ankara, Turkey: epidemiological and clinical features of rotavirus-positive versus rotavirus negative cases. A Karadag, Z C Acikgoz, Z Avci, F Catal, S Gocer, S Gamberzade, Scand J Infect Dis. 37Karadag A, Acikgoz ZC, Avci Z, Catal F, Gocer S, Gamberzade S, et al. Childhood diarrhoea in Ankara, Turkey: epidemiological and clin- ical features of rotavirus-positive versus rotavirus negative cases. Scand J Infect Dis. 2005;37:269-75.

Molecular characterization of rotaviruses circulating in the population in Turkey. O Cataloluk, M Iturriza, J Gray, Epidemiol Infect. 133Cataloluk O, Iturriza M, Gray J. Molecular characterization of rotaviruses circulating in the population in Turkey. Epidemiol Infect. 2005;133:673-8.

Molecular characterization of serotype G9 rotavirus strains from a global collection. M Ramachandran, C D Kirkwood, L Unicomb, N A Cunliffe, R L Ward, M K Bhan, Virology. 278Ramachandran M, Kirkwood CD, Unicomb L, Cunliffe NA, Ward RL, Bhan MK, et al. Molecular characterization of serotype G9 rotavirus strains from a global collection. Virology. 2000; 278:436-44.