Simple Sequence Repeat and <Emphasis Type="Italic">S</Emphasis>-locus Genotyping to Explore Genetic Variability in Polyploid <Emphasis Type="Italic">Prunus spinosa</Emphasis> and <Emphasis Type="Italic">P. insititia</Emphasis>

Polyploid Prunus spinosa (2n = 4×) and P. insititia (2n = 6×) represent enormous genetic potential in Central Europe, which can be exploited in breeding programmes. In Hungary, 17 cultivar candidates were selected from wild-growing populations including 10 P. spinosa, 4 P. insititia and three P. spinosa × P. domestica hybrids (2n = 5×). Their taxonomic classification was based on their phenotypic characteristics. Six simple sequence repeats (SSRs) and the multiallelic S-locus genotyping were used to characterize genetic variability and reliable identification of the tested accessions. A total of 98 SSR alleles were identified, which presents 19.5 average allele number per locus, and each of the 17 genotypes could be discriminated based on unique SSR fingerprints. A total of 23 S-RNase alleles were identified. The complete and partial S-genotype was determined for 8 and 9 accessions, respectively. The identification of a cross-incompatible pair of cultivar candidates and several semi-compatible combinations help maximize fruit set in commercial orchards. Our results indicate that the S-allele pools of wild-growing P. spinosa and P. insititia are overlapping in Hungary. A phylogenetic and principal component analysis confirmed the high level of diversity and genetic differentiation present within the analysed genotypes and helped clarify doubtful taxonomic identities. Our data confirm that S-locus genotyping is suitable for diversity studies in polyploid Prunus species. The analysed accessions represent huge genetic potential that can be exploited in commercial cultivation.     

Inheritance and gene mapping of the white flower trait in <Emphasis Type="Italic">Brassica juncea</Emphasis>

Despite being a unique marker trait, white flower inheritance in Brassica juncea remains poorly understood at the gene level. In this study, we investigated a B. juncea landrace with white petal in China. The white petal phenotype possessed defective chromoplasts with less plastoglobuli than the yellow petal phenotype. Genetic analysis confirmed that two independent recessive genes (Bjpc1 and Bjpc2) controlled the white flower trait. We then mapped the BjPC1 gene in a BC₄ population comprising 2295 individuals. We identified seven AFLP (amplified fragment length polymorphism) markers closely linked to the white flower gene. BLAST search revealed the sequence of AFLP fragments were highly homologous with the Scaffold000085 and Scaffold000031 sequences on the A02 chromosome in the Brassica rapa genome. Based on this sequence homology, we developed simple sequence repeat (SSR) primer pairs and identified 13 SSRs linked to the BjPC1 gene, including two that were co-segregated (SSR9 and SSR10). The two closest markers (SSR4 and SSR11) were respectively 0.9 and 0.4 cM on either side of BjPC1. BLAST analysis revealed that these marker sequences corresponded highly to A02 in B. juncea. They were mapped within a 33 kb genomic region on B. rapa A02 (corresponds to a 40 kb genomic region on B. juncea A02) that included three genes. Sequence BjuA008406, homologous to AtPES2 in Arabidopsis thaliana and Bra032956 in B. rapa, was the most likely candidate for BjPC1. These results should accelerate BjPC1 cloning and facilitate our understanding of the molecular mechanisms controlling B. juncea petal color.     

Natural genetic variation in <Emphasis Type="Italic">Brassica</Emphasis> homologs of <Emphasis Type="Italic">FLOWERING LOCUS T</Emphasis> and characterization of its expression domains

FLOWERING LOCUS T (FT), a major effect gene, regulates flowering time in Arabidopsis. We analyzed evolutionary changes distinguishing two FT homeologous loci in B. rapa, described genetic variation in homologs isolated and reported expression pattern of FT in B. juncea. Synteny analysis confirmed presence of two FT genomic copies in B. rapa ssp. pekinensis and resolved pre-existing anomalies regarding copy number in “AA” genome. Synteny analysis of B. rapa homeologous regions CR1 (129 kb) and CR2 (232 kb) revealed differential gene fractionation and wide-spread re-arrangements. Seven genomic DNA (gDNA) variants (2.1–2.2 kb) and 10 complementary DNA (cDNA) variants (528 bp) were isolated from 6 Brassica species. The gDNA variants shared 72–99 % similarity within Brassica and 58–60 % between Arabidopsis and Brassica. FT cDNA variants shared 92–100 % similarity within Brassica and 87 % between Arabidopsis and Brassica. Phylogenetic analysis of FT gDNA, cDNA and protein sequences revealed two major clades, differentiating homologs derived from species containing shared “BB” and “CC” genomes. Phylogram based on Brassica FT gDNA differentiated homeologs derived from AA-LF (Least fractioned) and AA-MF1 (Moderately fractioned) sub-genomes. Analysis of FT expression pattern in B. juncea revealed increasing levels correlating with attainment of physiological maturity; highest levels were detected in older leaves implying conservation in spatio-temporal expression pattern vis-à-vis Arabidopsis. In conclusion, our study reveals that polyploidy in Brassicas resulted in expansion of FT gene copies with homologs charting independent evolutionary course through accumulation of mutations. However, expression domains of FT remained conserved across Brassicaceae to preserve the critical function of FT in controlling flowering time.     

Development of cost-effective single nucleotide polymorphism marker assays for genetic diversity analysis in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Competitive allele-specific PCR (KASPar) assay is a user-friendly system that provides flexibility in the numbers of single nucleotide polymorphisms (SNPs) and genotypes. Based on Illumina-GA-IIx genomic data from 10 genotypes with a broad genetic background, 3183 SNPs were selected for KASPar assays development, and 568 were finally converted and selected for Brassica rapa germplasm characterization (17.8%) on the basis of reproducibility, missing data rate, and uniform genetic distribution. High levels of polymorphism of these markers across 231 B. rapa genotypes were verified, illustrating by high polymorphic information content (averaged 0.34), minor allele frequency (0.37), genetic diversity (0.45), and the low observed heterozygosity (0.10). Based on the SNP dataset, structure and principal coordinates analysis, and neighbor-joining phylogenetic methods were used to examine the population structure and gave highly consistent results. The 231 accessions were divided into the four primary subspecies, representing 99 accessions from B. rapa ssp. pekinensis, 85 from B. rapa ssp. chinensis, 30 from B. rapa ssp. rapifera, and 17 from B. rapa ssp. oleifera and were further subdivided into 12 lower-order clusters according to different morphotypes. The genetic variability and pairwise fixation index analysis revealed that the ssp. pekinensis accessions possess the most extensive genetic variation among the four subspecies. The KASPar system is highly useful for validating SNPs and will be valuable for genetics research and breeding applications in B. rapa.     

Development of <Emphasis Type="Italic">Brassica oleracea</Emphasis>-<Emphasis Type="Italic">nigra</Emphasis> monosomic alien addition lines: genotypic,cytological and morphological analyses

Key message

We report the development and characterization of Brassica oleracea - nigra monosomic alien addition lines (MAALs) to dissect the Brassica B genome.

Abstract

Brassica nigra (2n = 16, BB) represents the diploid Brassica B genome which carries many useful genes and traits for breeding but received limited studies. To dissect the B genome from B. nigra, the triploid F₁ hybrid (2n = 26, CCB) obtained previously from the cross B. oleracea var. alboglabra (2n = 18, CC) × B. nigra was used as the maternal parent and backcrossed successively to parental B. oleracea. The progenies in BC₁ to BC₃ generations were analyzed by the methods of FISH and SSR markers to screen the monosomic alien addition lines (MAALs) with each of eight different B-genome chromosomes added to C genome (2n = 19, CC + 1B_1?8), and seven different MAALs were established, except for the one with chromosome B2 which existed in one triple addition. Most of these MAALs were distinguishable morphologically from each other, as they expressed the characters from B. nigra differently and at variable extents. The alien chromosome remained unpaired as a univalent in 86.24% pollen mother cells at diakinesis or metaphase I, and formed a trivalent with two C-genome chromosomes in 13.76% cells. Transmission frequency of all the added chromosomes was far higher through the ovules (averagely 14.40%) than the pollen (2.64%). The B1, B4 and B5 chromosomes were transmitted by female at much higher rates (22.38–30.00%) than the other four (B3, B6, B7, B8) (5.04–8.42%). The MAALs should be valuable for exploiting the genome structure and evolution of B. nigra.

     

Mapping of genetic locus for leaf trichome in <Emphasis Type="Italic">Brassica oleracea</Emphasis>

Key message

The genetic locus for leaf trichome was identified via marker-based mapping and SNP microarray assay, and a functional marker was developed to facilitate the breeding for hairiness in Brassica oleracea.

Abstract

Plant trichomes are involved in various functions particularly in protecting plants against some biotic and abiotic damages. In the present study, an F2 segregating population was developed from the cross between a glabrous cultivated B. oleracea (CC, 2n = 18) and a hairy wild relative, B. incana (CC, 2n = 18). A 1:3 segregation pattern between glabrous and hairy plants was detected among 1063 F2 genotypes, and the locus for hairiness was mapped in a 4.3-cM genetic region using 267 SSR markers among 149 F2 genotypes, corresponding to a 17.6-Mb genomic region on chromosome C01. To narrow the genetic region for hairiness, the Brassica 60 K SNP Bead Chip Arrays were applied to genotype 64 glabrous and 30 hairy F2 plants, resulting in a 1.04-Mb single peak region located in the 17.6-Mb interval. A candidate gene, BoTRY, was identified by qRT-PCR which revealed significant higher expression in glabrous F2 genotypes as compared with that in hairy plants. A cleaved amplified polymorphic site marker was successfully developed to distinguish the sequence variations of BoTRY between hairy and glabrous plants. Our study will be helpful for molecular breeding for hairiness in B. oleracea.

     

Identification,characterization, validation and cross-species amplification of genic-SSRs in Indian Mustard (<Emphasis Type="Italic">Brassica juncea</Emphasis>)

Brassica juncea is an economically important oilseed crop worldwide. It has limited genomic resources at present. We generated 47,962,057 expressed sequence reads which were assembled into 45,280 unigenes. A total of 4108 SSR loci (≥10 bp) were identified in these unigenes. Trinucleotide was the most frequent repeat unit (59.91 %) followed by di- (38.66 %), tetra - (0.71 %), hexa - (0.49 %) and pentanucleotide repeats (0.24 %). Primers were designed for 2863 SSR loci among which 460 were selected for primer synthesis. A total of 339 loci amplified successfully of which 134 (39.5 %) exhibited polymorphism among six B. juncea genotypes with PIC values ranging from 0.18 to 0.81. Further, 25 polymorphic SSRs were used for analysis of genetic variability in 25 genotypes of Brassicas and their wild relatives. Two to five alleles with PIC values 0.22–0.66 were detected at these loci. The dendrogram grouped the genotypes according to their known pedigree/systematic position.     

Genotyping and Persistence of <Emphasis Type="Italic">Candida albicans</Emphasis> from Pregnant Women with Vulvovaginal Candidiasis

Objective

To study Candida albicans genotypes using RAPD and their susceptibility to fluconazole in healthy pregnant women and in vulvovaginal candidiasis (VVC) patients after topical treatment with clotrimazole.

Methods

Vaginal swabs were collected at t = 0 and t = 1 (1 month later) in pregnant women (control group, n = 33), and before (t = 0), at 1 month (t = 1) and at 2 months (t = 2) after clotrimazole treatment in pregnant women with VVC.

Results

Candida albicans was isolated in 30% of healthy pregnant women and 80% of patients with VVC. A high genetic heterogeneity was observed in C. albicans genotypes between individuals. In patients with VVC, topical antifungal treatment with clotrimazole was clinically effective, but only in a 62% C. albicans was eradicated. In patients in which C. albicans was not eradicated, this microorganism persisted for 1 or 2 months after the antifungal treatment. The persistent colonies were not associated with a specific genotype, but they were associated with higher MICs in comparison with colonies isolated from the control group.

Conclusions

Therapy with topical clotrimazole, despite a good clinical outcome, could not eradicate completely C. albicans allowing the persistence of genotypes, with higher MICs to fluconazole. More studies with higher number of patients are needed to validate this preliminary finding.

     

Genome-wide association mapping of resistance to eyespot disease (<Emphasis Type="Italic">Pseudocercosporella herpotrichoides</Emphasis>) in European winter wheat (<Emphasis Type="Italic">Triticum aestivum</Emphasis> L.) and fine-mapping of <Emphasis Type="Italic">Pch1</Emphasis>

Key message

Genotypes with recombination events in the Triticum ventricosum introgression on chromosome 7D allowed to fine-map resistance gene Pch1, the main source of eyespot resistance in European winter wheat cultivars.

Abstract

Eyespot (also called Strawbreaker) is a common and serious fungal disease of winter wheat caused by the necrotrophic fungi Oculimacula yallundae and Oculimacula acuformis (former name Pseudocercosporella herpotrichoides). A genome-wide association study (GWAS) for eyespot was performed with 732 microsatellite markers (SSR) and 7761 mapped SNP markers derived from the 90 K iSELECT wheat array using a panel of 168 European winter wheat varieties as well as three spring wheat varieties and phenotypic evaluation of eyespot in field tests in three environments. Best linear unbiased estimations (BLUEs) were calculated across all trials and ranged from 1.20 (most resistant) to 5.73 (most susceptible) with an average value of 4.24 and a heritability of H ² = 0.91. A total of 108 SSR and 235 SNP marker–trait associations (MTAs) were identified by considering associations with a ?log₁₀ (P value) ≥3.0. Significant MTAs for eyespot-score BLUEs were found on chromosomes 1D, 2A, 2D, 3D, 5A, 5D, 6A, 7A and 7D for the SSR markers and chromosomes 1B, 2A, 2B, 2D, 3B and 7D for the SNP markers. For 18 varieties (10.5%), a highly resistant phenotype was detected that was linked to the presence of the resistance gene Pch1 on chromosome 7D. The identification of genotypes with recombination events in the introgressed genomic segment from Triticum ventricosum harboring the Pch1 resistance gene on chromosome 7DL allowed the fine-mapping of this gene using additional SNP markers and a potential candidate gene Traes_7DL_973A33763 coding for a CC-NBS-LRR class protein was identified.

     

Feed-backs between genetic structure and perturbation-driven decline in seagrass (<Emphasis Type="Italic">Posidonia oceanica</Emphasis>) meadows

We explored the relationships between perturbation-driven population decline and genetic/genotypic structure in the clonal seagrass Posidonia oceanica, subject to intensive meadow regression around four Mediterranean fish-farms, using seven specific microsatellites. Two meadows were randomly sampled (40 shoots) within 1,600 m² at each site: the “impacted” station, 5–200 m from fish cages, and the “control” station, around 1,000 m downstream further away (considered a proxy of the pre-impact genetic structure at the site). Clonal richness (R), Simpson genotypic diversity (D*) and clonal sub-range (CR) were highly variable among sites. Nevertheless, the maximum distance at which clonal dispersal was detected, indicated by CR, was higher at impacted stations than at the respective control station (paired t-test: P < 0.05, N = 4). The mean number of alleles (Â) and the presence of rare alleles ( _r) decreased at impacted stations (paired t-test: P < 0.05, and P < 0.02, respectively, N = 4). At a given perturbation level (quantified by the organic and nutrient loads), shoot mortality at the impacted stations significantly decreased with CR at control stations (R ² = 0.86, P < 0.05). Seagrass mortality also increased with  (R ² = 0.81, P < 0.10), R (R ² = 0.96, P < 0.05) and D* (R ² = 0.99, P < 0.01) at the control stations, probably because of the negative correlation between those parameters and CR. Therefore, the effects of clonal size structure on meadow resistance could play an important role on meadow survival. Large genotypes of P. oceanica meadows thus seem to resist better to fish farm-derived impacts than little ones. Clonal integration, foraging advantage or other size-related fitness traits could account for this effect.     

Identification of candidate genes of QTLs for seed weight in <Emphasis Type="Italic">Brassica napus</Emphasis> through comparative mapping among <Emphasis Type="Italic">Arabidopsis</Emphasis> and <Emphasis Type="Italic">Brassica</Emphasis>species

Background

Map-based cloning of quantitative trait loci (QTLs) in polyploidy crop species remains a challenge due to the complexity of their genome structures. QTLs for seed weight in B. napus have been identified, but information on candidate genes for identified QTLs of this important trait is still rare.

Results

In this study, a whole genome genetic linkage map for B. napus was constructed using simple sequence repeat (SSR) markers that covered a genetic distance of 2,126.4 cM with an average distance of 5.36 cM between markers. A procedure was developed to establish colinearity of SSR loci on B. napus with its two progenitor diploid species B. rapa and B. oleracea through extensive bioinformatics analysis. With the aid of B. rapa and B. oleracea genome sequences, the 421 homologous colinear loci deduced from the SSR loci of B. napus were shown to correspond to 398 homologous loci in Arabidopsis thaliana. Through comparative mapping of Arabidopsis and the three Brassica species, 227 homologous genes for seed size/weight were mapped on the B. napus genetic map, establishing the genetic bases for the important agronomic trait in this amphidiploid species. Furthermore, 12 candidate genes underlying 8 QTLs for seed weight were identified, and a gene-specific marker for BnAP2 was developed through molecular cloning using the seed weight/size gene distribution map in B. napus.

Conclusions

Our study showed that it is feasible to identify candidate genes of QTLs using a SSR-based B. napus genetic map through comparative mapping among Arabidopsis and B. napus and its two progenitor species B. rapa and B. oleracea. Identification of candidate genes for seed weight in amphidiploid B. napus will accelerate the process of isolating the mapped QTLs for this important trait, and this approach may be useful for QTL identification of other traits of agronomic significance.

     

Epigenetic and genetic variation in <Emphasis Type="Italic">GATA5</Emphasis> is associated with gastric disease risk

Gastric cancer incidence varies considerably among populations, even those with comparable rates of Helicobacter pylori infection. To test the hypothesis that genetic variation plays a role in gastric disease, we assessed the relationship between genotypes and gastric histopathology in a Colombian study population, using a genotyping array of immune-related single nucleotide polymorphisms (SNPs). Two synonymous SNPs (rs6061243 and rs6587239) were associated with progression of premalignant gastric lesions in a dominant-effects model after correction for multiple comparisons (p = 2.63E?07 and p = 7.97E?07, respectively); effect sizes were β = ?0.863 and β = ?0.815, respectively, where β is an estimate of effect on histopathology scores, which ranged from 1 (normal) to 5 (dysplasia). In our replication cohort, a second Colombian population, both SNPs were associated with histopathology when additively modeled (β = ?0.256, 95 % CI = ?0.47, ?0.039; and β = ?0.239, 95 % CI = ?0.45, ?0.024), and rs6587239 was significantly associated in a dominant-effects model (β = ?0.330, 95 % CI = ?0.66, 0.00). Because promoter methylation of GATA5 has previously been associated with gastric cancer, we also tested for the association of methylation status with more advanced histopathology scores in our samples and found a significant relationship (p = 0.001). A multivariate regression model revealed that the effects of both the promoter methylation and the exonic SNPs in GATA5 were independent. A SNP-by-methylation interaction term was also significant. This interaction between GATA5 variants and GATA5 promoter methylation indicates that the association of either factor with gastric disease progression is modified by the other.     

Association between <Emphasis Type="Italic">cagA</Emphasis>, <Emphasis Type="Italic">vacAi</Emphasis>, and <Emphasis Type="Italic">dupA</Emphasis> genes of <Emphasis Type="Italic">Helicobacter pylori</Emphasis> and gastroduodenal pathologies in Chilean patients

In addition to the already known cagA gene, novel genetic markers have been associated with Helicobacter pylori (H. pylori) virulence: the dupA and vacAi genes. These genes might play an important role as specific markers to determine the clinical outcome of the disease, especially the vacAi gene, which has been expected to be a good marker of severe pathologies like gastric adenocarcinoma. In the present study, the association of cagA, dupA, and vacAi genes with gastroduodenal pathologies in Chilean patients was studied. One hundred and thirty-two patients positive for H. pylori were divided into two groups—non-severe and severe gastric pathologies—and investigated for the presence of cagA, dupA, and vacAi H. pylori virulence genes by PCR. The cagA gene was detected in 20/132 patients (15.2%), the vacAi1 gene was detected in 54/132 patients (40.9%), the vacAi2 gene was detected in 26/132 patients (19.7%), and the dupA gene was detected in 50/132 (37.9%) patients. Logistic regression model analysis showed that the vacAi1 isoform gene in the infected strains and the severity of the diseases outcome were highly associated, causing severe gastric damage that may lead to gastric cancer (p < 0.0001; OR = 8.75; 95% CI 3.54–21.64). Conversely, cagA (p = 0.3507; OR = 1.62; 95% CI 0.59–4.45) and vacAi2 (p = 0.0114; OR = 3.09; 95% CI 1.26–7.60) genes were not associated with damage, while the dupA gene was associated significantly with non-severe clinical outcome (p = 0.0032; OR = 0.25; 95% CI 0.09–0.65). In addition, dupA gene exerts protection against severe gastric pathologies induced by vacAi1 by delaying the outcome of the disease by approximately 20 years.     

Identification of a nuclear-recessive gene locus for male sterility on A2 chromosome using the Brassica 60 K SNP array in non-heading Chinese cabbage

WS24-3 is a newly bred recessive genic male sterility line of the non-heading Chinese cabbage (Brassica rapa ssp. chinensis). Here, an F₂ population was produced from the cross between WS24-3 and a male-fertile breeding line (WS135). The Illumina Brassica 60 K single nucleotide polymorphism (SNP) array was used for SNPs detecting between sterile and fertile bulks from the F₂ population, and 62 SNPs were identified. BLAST analysis of the 62 SNPs revealed that the A2 chromosome of Brassica rapa genome contained 22 SNPs, whereas the other chromosomes did not contain more than 6 SNPs each. These data indicated that the potential target gene locus, named Bra2Ms, might be located on A2. Based on 10 of the 22 SNPs, allele-specific-polymerase chain reaction (AS-PCR) primers and single sequence repeat (SSR) primers were designed, 5 AS-PCR primers and 9 SSR primers showed difference between the bulks in electrophoretic determination. Analysis of these markers in F₂ population revealed that Bra2Ms was genetically delimited to a region of 1.2 cM. We also detected two co-segregated markers SSRa2-951 and SSRa2-960 in this region. The markers identified in our study might facilitate the transfer of recessive genic male sterility alleles to other favorable genetic backgrounds. Furthermore, these markers will support a map-based clone of Bra2Ms.     

Microsatellite-based genetic diversity analyses of <Emphasis Type="Italic">sugary1</Emphasis>-, <Emphasis Type="Italic">shrunken2</Emphasis>- and double mutant- sweet corn inbreds for their utilization in breeding programme

Sweet corn has recently experienced sharp rise in demand worldwide. Recessive sugary1 (su1) and shrunken2 (sh2) that enhances kernel sweetness have been abundantly used in sweet corn breeding. Analyses of genetic diversity among sweet corn inbreds assume great significance for their effective utilization in hybrid breeding. A set of 48 diverse sweet corn genotypes encompassing su1su1, sh2sh2 and su1su1/sh2sh2 types were analyzed using 56 microsatellite markers. A total of 213 alleles with mean of 3.8 alleles per locus were generated. Two unique- and 12 rare- alleles were identified. The average PIC and genetic dissimilarity was 0.50 and 0.73, respectively. Cluster analysis grouped the inbreds into three major clusters, with each of the su1su1-, sh2sh2- and su1su1/sh2sh2-types were broadly clustered together. Principal coordinate analyses also depicted the diverse origin of the genotypes. The study identified inbreds for synthesis of pools and pedigree populations to develop novel inbreds. The study led to the identification of prospective heterotic combinations in various genetic backgrounds (sh2sh2 × sh2sh2, su1su1 × su1su1, su1su1/sh2sh2 × su1su1/sh2sh2, sh2sh2 × su1su1/sh2sh2 and su1su1 × su1su1/sh2sh2).     

Gene loss/retention and evolutionary pattern of ascorbic acid biosynthesis and recycling genes in <Emphasis Type="Italic">Brassica rapa</Emphasis> following whole genome triplication

Ascorbic acid (AsA) is an inevitable antioxidant found abundantly in higher plants. Despite the importance of AsA in plants, how AsA biosynthesis (ABGs; d-mannose/l-galactose pathway) and AsA recycling genes (ARGs) evolved through polyploidization has not been addressed to date. Here, we evaluated the impacts of whole genome triplication (WGT) on ABGs and ARGs in Chinese cabbage (Brassica rapa ssp. pekinensis), which diverged from Arabidopsis thaliana before the WGT event. Twenty-three ABGs coded in 13 loci representing nine different enzyme classes and 29 ARGs coded in 19 loci representing five different enzyme classes were identified in the B. rapa genome by whole-genome screening through comparative genomic analyses. Five of these loci maintained three gene copies, 10 loci maintained two gene copies and the majority of the loci (n = 17) maintained single gene copies. Segmental (62 %) and tandem duplication (6 %), and fragment (21 %) and large-scale recombination (10 %) events accelerated the diversification of ABGs and ARGs. Thirteen of the 52 (25 %) identified genes experienced intron losses and two (4 %) experienced intron gains implying that intron losses outnumbered intron gains. The expansion and the retention of ABGs and ARGs were similar to the whole genome gene expansion and retention (P > 0.05). These findings provide new insights into the structural characteristics and evolutionary trends of ABGs and ARGs. In addition, our data could become a useful resource to further the functional characterization of these genes.     

Genetic Diversity in Various Accessions of Pineapple [<Emphasis Type="Italic">Ananas comosus</Emphasis> (L.) Merr.] Using ISSR and SSR Markers

Inter simple sequence repeat (ISSR) and simple sequence repeat (SSR) markers were used to assess the genetic diversity of 36 pineapple accessions that were introduced from 10 countries/regions. Thirteen ISSR primers amplified 96 bands, of which 91 (93.65%) were polymorphic, whereas 20 SSR primers amplified 73 bands, of which 70 (96.50%) were polymorphic. Nei’s gene diversity (h = 0.28), Shannon’s information index (I = 0.43), and polymorphism information content (PIC = 0.29) generated using the SSR primers were higher than that with ISSR primers (h =  0.23, I = 0.37, PIC = 0.24), thereby suggesting that the SSR system is more efficient than the ISSR system in assessing genetic diversity in various pineapple accessions. Mean genetic similarities were 0.74, 0.61, and 0.69, as determined using ISSR, SSR, and combined ISSR/SSR, respectively. These results suggest that the genetic diversity among pineapple accessions is very high. We clustered the 36 pineapple accessions into three or five groups on the basis of the phylogenetic trees constructed based on the results of ISSR, SSR, and combined ISSR/SSR analyses using the unweighted pair-group with arithmetic averaging (UPGMA) method. The results of principal components analysis (PCA) also supported the UPGMA clustering. These results will be useful not only for the scientific conservation and management of pineapple germplasm but also for the improvement of the current pineapple breeding strategies.