Genetic diversity of <Emphasis Type="Italic">avenin-like b</Emphasis> genes in <Emphasis Type="Italic">Aegilops tauschii</Emphasis> Coss

Avenin-like storage proteins influence the rheological properties and processing quality in common wheat, and the discovery of new alleles will benefit wheat quality improvement. In this study, 13 avenin-like b alleles (TaALPb7D-A–M) were discovered in 108 Aegilops tauschii Coss. accessions. Ten alleles were reported for the first time, while the remaining three alleles were the same as alleles in other species. A total of 15 nucleotide changes were detected in the 13 alleles, resulting in only 11 amino acid changes because of synonymous mutations. Alleles TaALPb7D-E, TaALPb7D-G, and TaALPb7D-J encoded the same protein. These polymorphic sites existed in the N-terminus, Repetitive region (Left), Repetitive region (Right) and C-terminus domains, with no polymorphisms in the signal peptide sequence nor in those encoding the 18 conserved cysteine residues. Phylogenetic analysis divided the TaALPb7Ds into four clades. The Ae. tauschii alleles were distributed in all four clades, while the alleles derived from common wheat, TaALPb7D-G and TaALPb7D-C, belonged to clade III and IV, respectively. Alleles TaALPb7D-G and TaALPb7D-C were the most widely distributed, being present in nine and six countries, respectively. Iran and Turkey exhibited the highest genetic diversity with respect to TaALPb7D alleles, accessions from these countries carrying seven and six alleles, respectively, which implied that these countries were the centers of origin of the avenin-like b gene. The new alleles discovered and the phylogenetic analysis of avenin-like b genes will provide breeding materials and a theoretical basis for wheat quality improvement.     

Development and validation of genomic simple sequence repeat markers in <Emphasis Type="Italic">Erianthus arundinaceus</Emphasis>

Erianthus arundinaceus, a member of the Saccharum complex, is of interest as a potential resource for sugarcane improvement and as a bioenergy crop. Genetic analyses of germplasm collections of E. arundinaceus are being used increasingly. To expand the genomic resources in E. arundinaceus, we aimed at developing simple sequence repeat markers. Using pyrosequencing on the 454 GS FLX system, we sequenced genomic DNA from “JW630” collected in Japan. A total of 1682 candidate loci were used to design the primers, and 1234 primer pairs amplified fragments of the expected size in the primer screening with three wild E. arundinaceus accessions (JW630, “JW4,” and “IJ76-349”). The efficiency of genotyping was validated with a subset of 174 primer pairs and 8 E. arundinaceus accessions. Of these primer pairs, 171 amplified fragments in all accessions tested and 162 detected polymorphic loci. The average values of genetic parameters were estimated as 0.30 (range, 0.09–0.49) for polymorphic information content, 1.65 (0.00–5.87) for marker index, and 2.78 (0.00–8.75) for resolving power. Using these parameters, we selected 61 primer pairs with large discriminatory power for the analyzed loci. Of the 174 primer pairs, 45 (25.9%) were also applicable to Saccharum and 33 (19.0%) to Miscanthus species. These markers would provide a valuable tool for estimating genetic diversity and constructing linkage maps in E. arundinaceus, which would be useful for genetic study and breeding.     

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.     

Prolonged culture of <Emphasis Type="Italic">Boesenbergia rotunda</Emphasis> cells reveals decreased growth and shoot regeneration capacity

We evaluated the ploidy levels and tissue culture responses of 16 Japanese Miscanthus accessions, which are registered and vegetatively maintained in the National Agriculture and Food Research Organization GeneBank, Japan, to screen suitable genotypes for the molecular breeding of Miscanthus species. A ploidy analysis showed that most M. sinensis and M. sinensis var. condensatus (var. condensatus) were putative diploids, but one accession identified as M. sinensis was unexpectedly a putative tetraploid. Additionally, M. sacchariflorus and its hybrid accessions were putative tetraploids. The deoxyribonucleic acid levels in var. condensatus were significantly higher than those in the diploid M. sinensis. Of the accessions, 10, including M. sinensis and var. condensatus, could induce plant regenerable embryogenic calli from apical meristems. We selected three of these M. sinensis accessions for further experiments because their calli growth rates were faster than those of the var. condensatus accessions. Tissue culture experiments with the selected accessions indicated that the frequencies of callus and green shoot formation strongly correlated with genotype. The broad-sense heritabilities of the embryogenic callus and green shoot formation frequencies in the selected accessions were 0.75 and 0.65, respectively, indicating that the cultures’ responses were mainly controlled by genetic factors. Thus, we further selected one accession that had the highest efficiencies in callus and green shoot formation, and we observed that light during callus culturing significantly inhibited calli growth, but promoted plant regeneration from calli in the selected accession.     

<Emphasis Type="Italic">Otidea</Emphasis> species from China,three new species with comments on some previously described species

The genus Otidea was recently monographed and studied phylogenetically, but knowledge of the diversity and distribution of Otidea species in China is fragmentary. In this study, collections from China were examined morphologically and included in phylogenetic analyses. Using LSU, TEF1-α, and RPB2 new species were placed within previously recognized clades in the genus. The results agree with both Genealogical Concordance Phylogenetic Species Recognition (GCPSR) and genetic divergence as previously reported. Three new species, Otidea hanseniae, Otidea korfii and Otidea purpureogrisea are recognized based on phylogenetic reconstruction using ITS, LSU, TEF1- α and RPB2. Comments on some incompletely known species are added. With the discovery of these three new species, the genus Otidea in China proves to be more diverse than previously recognized.     

Genetic structure of cultivated <Emphasis Type="Italic">Zanthoxylum</Emphasis> species investigated with SSR markers

Zanthoxylum is an economically and ecologically important genus of the Rutaceae family, of which Z. bungeanum and Z. armatum have a long history of cultivation in China. However, how the natural processes such as selection and drift and agriculture practices have influenced the genetic variation of cultivated Zanthoxylum species during long-term domestication remains elusive. Herein, we determined the population genetic structure of current widely cultivated Zanthoxylum species, Z. bungeanum and Z. armatum. Microsatellite markers revealed a high level of genetic variation and significant genetic differentiation for both species despite Z. bungeanum showed higher genetic diversity than Z. armatum. AMOVA indicated that most of the genetic variation exists within individuals rather than among provenances for both species. Population structure analyses generated three distinct groups within the entire accessions. All Z. bungeanum accessions were distinguished into two major geographic groups, north and south groups, with Qinling Mountains as the main geographic barrier to gene flow while a significant genetic differentiation was observed between cultivated and wild Z. armatum accessions. Mantel test of Z. bungeanum displayed a significant correlation between genetic and geographic distances within each inferred group but no correlation between genetic and geographic distance was observed when comparing genetic and geographic distances focusing only on pairwise of north vs. south provenances, ruling out the hypothesis that gene flow between north and south provenances followed an isolation-by-distance model. Our research provided a fundamental genetic profile that will improve the conservation and responsible exploitation of the extant germplasm of Zanthoxylum.     

A molecular phylogeny of the Old World cypresses (<Emphasis Type="Italic">Cupressus</Emphasis>: Cupressaceae): evidence from nuclear and chloroplast DNA sequences

Studies of phylogenetic relationships among cypresses of the Old World (Cupressus; Cupressaceae) have been plagued by unresolved relationships, poor branch support, and conflict between data sets and methods of analysis. In this study, we combined 5.4 kb of aligned DNA sequence and 157 binary characters with previously published data in examining phylogenetic relationships among Cupressus species. Bayesian and parsimony analysis of the combined data or of the nuclear data alone always recovered three principal clades of Cupressus; however, tests of phylogenetic incongruence could not distinguish between competing relationships among the three principal Cupressus lineages. In contrast, incongruence tests often found statistically significant conflict between the nuclear and plastid data, particularly with respect to the placement of C. chengiana. Consistent with previous studies and prevailing taxonomic opinion, we find C. darjeelingensis more closely related to cypresses of the New World (Hesperocyparis). In contrast, we placed accessions of C. assamica and C. tonkinensis, two putatively Old World species suggested to be misidentified New World taxa by some authors, within well-supported Old World clades. Statistical analysis of genetic distances suggests instances in which taxa recognized as distinct species by some authors are identical or nearly so and may best be considered a single taxon. Conversely, we identify instances in which infraspecific taxa are more distantly related to one another than those traditionally recognized as distinct species. Factors confounding cypress taxonomies, including poor morphological differentiation, misidentification, and the use of accessions of questionable provenance, are discussed.     

Phylogenetic and morphological analysis of a new cliff-dwelling species reveals a remnant ancestral diversity and evolutionary parallelism in <Emphasis Type="Italic">Sonchus</Emphasis> (Asteraceae)

We describe a new cliff-dwelling species within Sonchus (Asteraceae): Sonchus boulosii and analyze its systematic position and evolutionary significance; in addition, we provide a key to the species of Sonchus in Morocco. Both morphological and ecological characteristics suggest a close relationship of S. boulosii with taxa of section Pustulati. However, ITS nrDNA and cpDNA matK markers indicate its uncertain position within the genus, but clear genetic differentiation from the remaining major clades. ITS phylogenetic trees show that likely evolutionary shifts to rocky habitat took place at least five times within genus Sonchus and that sect. Pustulati and S. boulosii clades have a clearly independent evolutionary origin. We postulate that the strong resemblance of S. boulosii to other rocky species reflects a phenomenon of homoplasy, probably driven by parallel evolutionary adaptations to the severe ecological constraints of its cliff face habitat. Therefore, a new section is also described, which includes S. boulosii as its sole representative: section Pulvinati. According to phylogenetic trees, the new clade may share its common ancestor with the clade comprising sections Maritimi and Arvenses, from which it is widely divergent in morphology and ecology, with the exception of Sonchus novae-zelandiae. However, the latter is a derived taxon, with high level of polyploidy unlike S. boulosii that shows 2n?=?18, basal chromosome number of the genus. Since sections Pulvinati and Pustulati seem to be quite old in Sonchus, we also hypothesize that some similarities, such as fruit morphology, may reflect the persistence of some primitive characteristics.     

Genetic Diversity and Phylogenetic Relationships of Two Closely Related Northeast China <Emphasis Type="Italic">Vicia</Emphasis> Species Revealed with RAPD and ISSR Markers

RAPD and ISSR analyses revealed genetic diversity and relationships among 11 populations of two closely related northeast China Vicia species, Vicia ramuliflora and V. unijuga. Both methods yielded similar and complementary results, showing high genetic diversity. Vicia ramuliflora had 100% polymorphic loci in both RAPD and ISSR, and V. unijuga had 100% polymorphic loci for RAPD and 98.96% for ISSR. Genetic differentiation was moderate among populations of each species. Genetic variation was distributed mainly within populations for the two species. The high level of gene flow was important for the allocation of genetic variation. The UPGMA dendrogram and principal coordinates analysis at the level of individuals and populations showed that V. ramuliflora and V. unijuga were more closely related than either of them was to the outgroup species, V. cracca. The small molecular variance of V. ramuliflora and V. unijuga supports the conclusion that these two species had a common ancestor.     

Assessment of genetic diversity in <Emphasis Type="Italic">Mucuna</Emphasis> species of India using randomly amplified polymorphic DNA and inter simple sequence repeat markers

Genus Mucuna which is native to China and Eastern India comprises of perennial climbing legume with long slender branches, trifoliate leaves and bear green or brown pod covered with soft or rigid hairs that cause intense irritation. The plants of this genus are agronomically and economically important and commercially cultivated in India, China and other regions of the world. The high degrees of taxonomical confusions exist in Mucuna species that make authentic identification and classification difficult. In the present study, the genetic diversity among the 59 accessions of six species and three varieties of M. pruriens has been assessed using DNA fingerprinting based molecular markers techniques namely randomly amplified polymorphic DNA (RAPD), inter simple sequence repeats (ISSR) and combined dataset of RAPD and ISSR. Also, genetic relationship among two endemic species of Mucuna namely M. imbricata and M. macrocarpa and two varieties namely IIHR hybrid (MHR) and Dhanwantari (MD) with other species under study was investigated by using cluster analysis and principal coordinate analysis. The cluster analysis of RAPD, ISSR and combined dataset of RAPD and ISSR clearly demonstrated the existence of high interspecific variation than intra-specific variation in genus Mucuna. The utility and efficacy of RAPD and ISSR for the study of intra species and interspecies genetic diversity was evident from AMOVA and PCoA analysis. This study demonstrates the genetic diversity in Mucuna species and indicates that these markers could be successfully used to assess genetic variation among the accessions of Mucuna species.     

Study on chloroplast DNA diversity of cultivated and wild pears (<Emphasis Type="Italic">Pyrus</Emphasis> L.) in Northern China

Eight pairs of chloroplast DNA (cpDNA) universal primers selected from 34 pairs were used to assess the genetic diversity of 132 pear accessions in Northern China. Among them, six amplified cpDNA fragments showed genetic diversity. A total of 24 variable sites, including 1 singleton variable site and 23 parsimony informative sites, as well as 21 insertion-deletion fragments, were obtained from the combined cpDNA sequences (5309–5535 bp). Two trnL-trnF-487 haplotypes, five trnL-trnF-413 haplotypes, five rbcL haplotypes, six trnS-psbC haplotypes, eight accD-psaI haplotypes and 12 rps16-trnQ haplotypes were identified among the individuals. Twenty-one haplotypes were identified based on the combined fragments. The values of nucleotide diversity (P_i), average number of nucleotide differences (k) and haplotype diversity (H_d) were 0.00070, 3.56408 and 0.7960, respectively. No statistical significance was detected in Tajima’s D test. Remarkably, the important cpDNA haplotypes and their representing accessions were identified clearly in this study. H_19 was considered as one of the ancient haplotypes and was a divergent centre. H_16 was the most common haplotype of the wild accessions. H_2 was the haplotype representing the most pear germplasm resources (46 cultivars and two wild Ussurian Pear accessions), followed by haplotype H_5 (30 cultivars, two wild Ussurian Pear accessions and four sand pears in outgroups) representing the cultivars ‘Dangshan Suli’ and ‘Yali’, which harbour the largest and the second largest cultivation areas in China. More importantly, this study demonstrated, for the first time, the supposed evolution routes of Pyrus based on cpDNA divergence in the background of pear phylogeny in Northern China.     

<Emphasis Type="Italic">In vitro</Emphasis> propagation of endangered <Emphasis Type="Italic">Mammillaria</Emphasis> genus (Cactaceae) species and genetic stability assessment using SSR markers

In vitro propagation protocols were established for endangered species of cacti Mammillaria hernandezii, M. dixanthocentron, and M. lanata. In vitro-germinated seedlings were used as the explant source. Three explant types were evaluated as apical, basal, and lateral stem sections. Shoot multiplication was achieved using Murashige and Skoog (MS) medium supplemented with benzyladenine, kinetin, meta-topolin, and thidiazuron in equimolar concentrations (0.0, 0.4, 1.1, 2.2, 4.4, and 8.9 μM). Shoot regeneration was obtained primarily in the lateral stem section explants. In M. hernandezii, an average of 7.4 shoots was regenerated in MS medium with 2.2 μM meta-topolin. M. dixanthocentron and M. lanata averaged 16.7 and 17.9 shoots/explant, respectively, in MS medium supplemented with 1.1 μM meta-topolin. Rooting occurred in MS medium without growth regulators. Three in vitro culture cycles were performed to validate the propagation protocols and to verify genetic stability. Shoots were collected in each cycle and genomic DNA was extracted. Amplified microsatellites were used to compare each genotype with its respective donor plant. Polymorphic information content analysis showed low levels of intra-clonal polymorphisms—M. hernandezii 0.04 and M. dixanthocentron and M. lanata both 0.12. More than 95% of the plants were successfully acclimatized in the greenhouse. After 12 months, plants of M. hernandezii reached the flowering stage; M. dixanthocentron and M. lanata flowered at 24 mo.     

AFLP,RAPD, and ISSR analysis of intraspecific polymorphism and interspecific differences of allotetraploid species <Emphasis Type="Italic">Aegilops kotschyi</Emphasis> Boiss. and <Emphasis Type="Italic">Aegilops variabilis</Emphasis> Eig

To evaluate genetic variation, 27 accessions of allotetraploid species Aegilops kotschyi and Ae. variabilis with the US genome were analyzed using the AFLP, RAPD, and ISSR methods. A total of 316 polymorphic RAPD fragments, 750 polymorphic AFLP fragments, and 234 polymorphic ISSR fragments were obtained. It was demonstrated that the analyzed species were characterized by a considerable level of nuclear genome variation. According to the data of ISSR and RAPD analysis, the average value of the Jaccard similarity coefficient for the accessions of Ae. variabilis from different geographical regions was slightly lower than that for the accessions of Ae. kotschyi. At the same time, AFLP analysis showed no considerable differences in the levels of intraspecific variation of the studied species. Analysis of the summarized RAPD, ISSR, and AFLP marking data in the Structure software program showed that most of the analyzed accessions with high degree of probability could be assigned to one of two groups, the first of which corresponded to Ae. kotschyi and the second corresponded to Ae. variabilis, thereby confirming the species independence of Ae. kotschyi and Ae. variabilis. Accessions k900, k907, k908, and v90 could not with a sufficiently high degree of probability be assigned to one of the species, which possibly was the result of interspecific hybridization. Analysis of the species diversity using different molecular markers made it possible to identify the accessions that were notably different from other accessions of its species.     

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.     

Assessment of genetic diversity in Ethiopian cowpea [<Emphasis Type="Italic">Vigna unguiculata</Emphasis> (L.) Walp.] germplasm using simple sequence repeat markers

The genetic diversity of cowpea (Vigna unguiculata L. Walp.) in Ethiopia was analyzed using 19 uniform accessions, 62 variable accessions (yielding 185 sub-types), and two mungbean (Vigna radiata) accessions (four subtypes) as outgroup. A set of 23 polymorphic simple sequence repeat (SSR) markers was identified, and polymorphism in the various accessions was scored by determining amplicon variability. Allele frequency, genetic diversity, and polymorphism information content (PIC) were determined for each SSR marker, and a neighbor joining dendrogram was generated to show the genetic relationship among the individual accessions. A total of 75 allelic variants was defined, with the average number of alleles per locus calculated to be three. The average genetic diversity (D) was 0.47, and PIC was 0.4. Three main clusters were identified by phylogenetic analysis, and the clusters and sub-grouping were supported by STRUCTURE and principal component analysis. This grouping had a moderate fixation index value of 0.075 and gene flow (Nm) of 3.176, indicating that the accessions possess wide diversity within individuals and populations. The accessions showed no clustering by geographical origins. Three well-characterized molecular markers (SSR1, C42-2B, and 61RM2) for race specific resistance to Striga gesnerioides in the cowpea cultivar B301 were used to evaluate the accessions for their potential for use in genetic improvement against this pest. Based on this analysis, only two accessions, 222890–2 from Gambela and 286–2 from the Southern Nations, Nationalities, and Peoples (SNNP) region, were found to cluster with B301 and contain the SSR1 resistance allele. These findings will assist in germplasm conservation efforts by the Institute of Biodiversity and Conservation of Ethiopia, and contribute to future studies aimed at the genetic improvement of local germplasm for improved overall agronomic performance as well as Striga resistance in particular.     

What can the geographic distribution of mtDNA haplotypes tell us about the invasion of New Zealand by house mice <Emphasis Type="Italic">Mus musculus</Emphasis>?

We mapped the distribution and diversity of mitochondrial D-loop haplotypes among 502 New Zealand house mice (Mus musculus). By widespread sampling from 74 sites, we identified 14 new haplotypes. We used Bayesian phylogenetic reconstructions to estimate the genetic relationships between the New Zealand representatives of Mus musculus domesticus (all six known clades) and M. m. castaneus (clade HG2), and mice from other locales. We defined four distinct geographic regions of New Zealand with differing haplotype diversity indices. Our Results suggest (a) two independent pre-1840 invasions by mice of different origin (domesticus clade E and castaneus clade HG2) at opposite ends of the country; (b) multiple later invasions by domesticus clades E and F accompanying the post-1840 development of New Zealand port facilities in the central regions, plus limited local incursions by domesticus clades A, B, C and D1; (c) a separate invasion of Chatham I. by castaneus clade HG2; (d) previously undescribed New Zealand haplotypes, potentially the products of localised indigenous mutation, and (e) hybridisation between different lineages.     

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.     

Genetic diversity and population structure of pummelo (<Emphasis Type="Italic">Citrus maxima</Emphasis>) germplasm in China

Pummelo (Citrus maxima) is one of the basic species of Citrus. It has been cultivated for about 4000 years in China, and therefore, there are abundant germplasm during the long time of culture. However, there is still a lack of a detailed study of the genetic characteristics of pummelo population. In this study, genetic diversity and population structure among 274 pummelo accessions collected in China were analyzed using 31 nuclear simple sequence repeat (nSSR) markers. The observed heterozygosity was calculated as 0.325 and genetic differentiation Fst as 0.077. Genetic structure analysis divided the whole germplasm into three subpopulations, Pop-a, Pop-b, and Pop-c. Pop-a was composed of accessions mostly from Southeast China, Pop-b was composed of accessions from the central region of South China, and Pop-c was composed of accessions from Southwest China. Meanwhile, the analysis of principal coordinate analysis and neighbor-join tree supported the viewpoint of three subpopulations, and then the possible dispersal routes of pummelos in China were proposed. This study provides an insight into the genetic diversity, facilitates future genome-wide association studies, and promotes the breeding program of pummelo as well.