Assessment of Genetic Diversity and Population Genetic Structure of Corylus mandshurica in China Using SSR Markers

Molecular evidence for multiple origins and high genetic differentiation of non-native winter crane fly, Trichocera maculipennis (Diptera: Trichoceridae), in the maritime Antarctic

Native biodiversity and ecosystems of Antarctica safeguarded from biological invasion face recent threats from non-native species, accelerated by increasing human activities and climate changes. Over two decades ago, the winter crane fly, Trichocera maculipennis, was first detected on King George Island. It has now successfully colonized several research stations across King George Island. To understand the origin, genetic diversity, and population structure of this Holarctic species, we conducted mitochondrial DNA cytochrome c oxidase subunit I (COI) sequence analysis across both its native and invasive ranges. In parallel, we performed microsatellite loci analysis within the invasive ranges, utilizing 12 polymorphic microsatellite markers. Furthermore, we compared body sizes among adult males and females collected from three different locations of King George Island. Our COI sequence analysis exhibited two different lineages present on King George Island. Lineage I was linked to Arctic Svalbard and Polish cave populations and Lineage II was related to Canadian Terra Nova National Park populations, implying multiple origins. Microsatellite analysis further exhibited high levels of genetic diversity and significant levels of genetic differentiation among invasive populations. Body sizes of adult T. maculipennis were significantly different among invasive populations but were not attributed to genetics. This significant genetic diversity likely facilitated the rapid colonization and establishment of T. maculipennis on King George Island, contributing to their successful invasion. Molecular analysis results revealed a substantial amount of genetic variation within invasive populations, which can serve as management units for invasive species control. Furthermore, the genetic markers we developed in the study will be invaluable tools for tracking impending invasion events and the travel routes of new individuals. Taken together, these findings illustrate the highly invasive and adaptable characteristics of T. maculipennis. Therefore, immediate action is necessary to mitigate their ongoing invasion and facilitate their eradication.

Microsatellites reveal divergence in population genetic diversity, and structure of osyris lanceolata (santalaceae) in Uganda and Kenya

BACKGROUND

Osyris lanceolata (Hochst. & Steud.) (Santalaceae) is a multipurpose plant highly valued culturally and economically in Africa. However, O. lanceolata populations have rapidly dwindled in East Africa due to overexploitation and this is believed to cause further consequences on the species' genetic diversity and structure within the region. Information regarding a species' genetic diversity and structure is necessary for conservation but this is currently lacking for O. lanceolata in Uganda and Kenya. Lack of adequate scientific data hinders conservation efforts hence threatening the species survival and livelihoods. This study investigated patterns in genetic diversity and structure of O. lanceolata in Uganda and Kenya. Ten polymorphic microsatellite loci were used to genotype 210 individuals: 96 from Ugandan and 114 from Kenyan populations.

RESULTS

All populations were highly polymorphic (80-100% polymorphism). A genetic differentiation was found between Kenyan and Ugandan populations. The highest genetic differentiation was among individuals and the least among populations. The Kenyan populations showed higher genetic diversity than Ugandan populations. The Ugandan populations showed more marker deviations from Hardy-Weinberg equilibrium and inbreeding coefficient. Two populations showed evidence of going through a recent bottleneck. There was significant genetic differentiation and structuring at higher K values into larger clusters and observed admixture between populations. The populations were significantly isolated by altitude as opposed to distance and climatic variables. Main barriers were associated with altitude differences. The data supports the idea of long-distance gene-flow between high altitude populations in both countries.

CONCLUSION

The divergence in genetic structure suggests unrecognised taxonomic units within O. lanceolata which are characteristic to lower altitudes and higher altitudes including most Kenyan populations with divergent evolutionary patterns. Geographical barriers and environmental gradients could have influenced this genetic divergence, and such patterns may escalate the species microevolutionary processes into full allopatric speciation. Further investigations into the species' genetic admixture and emerging taxonomic units are necessary to guide conservation strategies in the region.

Different rates of pollen and seed gene flow cause branch-length and geographic cytonuclear discordance within Asian butternuts

Topological cytonuclear discordance is commonly observed in plant phylogenetic and phylogeographic studies, yet few studies have attempted to detect two other forms of cytonuclear discordance (branch length and geographical) and to uncover the causes of the discordance. We used the whole nuclear and chloroplast genome data from 80 individual Asian butternuts to reveal the pattern and processes of cytonuclear discordance. Our findings indicate that the chloroplast genome had substantially deeper divergence (branch-length discordance) and a steeper cline in the contact zone (geographic discordance) compared with the nuclear genome. After various hypothesis have been tested, the results suggest that incomplete lineage sorting, positive selection and cytonuclear incompatibility are probably insufficient to explain this pattern. However, isolation-by-distance analysis and gene flow estimation point to a much higher level of gene flow by pollen compared with by seeds, which may have slowed down lineage divergence and mediated wider contact for nuclear genome compared with the chloroplast genome. Altogether, this study highlights a critical role of sex-biased dispersal in causing discordance between the nuclear and plastid genome of Asian butternuts. Given its ubiquity among plants, asymmetric gene flow should be given a high priority in future studies of cytonuclear discordance.

Species delimitation and hybridization history of a hazel species complex

BACKGROUND AND AIMS

Hybridization increases species adaptation and biodiversity but also obscures species boundaries. In this study, species delimitation and hybridization history were examined within one Chinese hazel species complex (Corylus chinensis-Corylus fargesii). Two species including four varieties have already been described for this complex, with overlapping distributions.

METHODS

A total of 322 trees from 44 populations of these four varieties across their ranges were sampled for morphological and molecular analyses. Climatic datasets based on 108 geographical locations were used to evaluate their niche differentiations. Flowering phenology was also observed for two co-occurring species or varieties.

KEY RESULTS

Four statistically different phenotypic clusters were revealed, but these clusters were highly inconsistent with the traditional taxonomic groups. All the clusters showed statistically distinct niches, with complete or partial geographical isolation. Only two clusters displayed a distributional overlap, but they had distinct flowering phenologies at the site where they co-occurred. Population-level evidence based on the genotypes of ten simple sequence repeat loci supported four phenotypic clusters. In addition, one cluster was shown to have an admixed genetic composition derived from the other three clusters through repeated historical hybridizations.

CONCLUSIONS

Based on our new evidence, it is better to treat the four clusters identified here as four independent species. One of them was shown to have an admixed genetic composition derived from the other three through repeated historical hybridizations. This study highlights the importance of applying integrative and statistical methods to infer species delimitations and hybridization history. Such a protocol should be adopted widely for future taxonomic studies.

A chromosome-level reference genome of the hazelnut, Corylus heterophylla Fisch

BACKGROUND

Corylus heterophylla Fisch. is a species of the Betulaceae family native to China. As an economically and ecologically important nut tree, C. heterophylla can survive in extremely low temperatures (-30 to -40 °C). To deepen our knowledge of the Betulaceae species and facilitate the use of C. heterophylla for breeding and its genetic improvement, we have sequenced the whole genome of C. heterophylla.

FINDINGS

Based on >64.99 Gb (∼175.30×) of Nanopore long reads, we assembled a 370.75-Mb C. heterophylla genome with contig N50 and scaffold N50 sizes of 2.07 and 31.33  Mb, respectively, accounting for 99.23% of the estimated genome size (373.61 Mb). Furthermore, 361.90 Mb contigs were anchored to 11 chromosomes using Hi-C link data, representing 97.61% of the assembled genome sequences. Transcriptomes representing 4 different tissues were sequenced to assist protein-coding gene prediction. A total of 27,591 protein-coding genes were identified, of which 92.02% (25,389) were functionally annotated. The phylogenetic analysis showed that C. heterophylla is close to Ostrya japonica, and they diverged from their common ancestor ∼52.79 million years ago.

CONCLUSIONS

We generated a high-quality chromosome-level genome of C. heterophylla. This genome resource will promote research on the molecular mechanisms of how the hazelnut responds to environmental stresses and serves as an important resource for genome-assisted improvement in cold and drought resistance of the Corylus genus.

Development of SSR markers based on transcriptome data and association mapping analysis for fruit shell thickness associated traits in oil palm (Elaeis guineensis Jacq.)

Present study mainly aimed to ascertain the distribution characteristics of gene-based microsatellite loci and to develop polymorphic SSR markers from the already available transcriptome data of Elaeis guineensis Jacq, an important oil crop. From this study, we identified the sum of 5791 SSRs across 51,425 unigenes from the transcripts of oil palm. We were able to evaluate 331primer pairs and characterized 183 polymorphic gene-based SSR markers. We identified a total of 506 alleles from the 183 polymorphic SSR loci, with an average of 2.77 alleles per locus. The characterized gene-based SSR markers from the transcriptome data of oil palm exhibited moderate levels of polymorphism with a significant level of heterozygosity ranges from 0.096 to 0.594 (mean = 0.336 ± 0.11). Among the identified SSR markers, sixty polymorphic markers were used to analyze genotypes of 55 oil palm accessions selected from three different provinces of China. Association mapping analysis provided the information of four markers that are associated with fruit shell thickness trait of oil palm. Among the four markers identified from association analysis, one SSR marker obtained from Unigene17150 is strictly associated with the oil palm fruit shell thickness trait.

Genetic Variation May Have Promoted the Successful Colonization of the Invasive Gall Midge, Obolodiplosis robiniae, in China

Invasive species often cause serious economic and ecological damage. Despite decades of extensive impacts of invasives on bio-diversity and agroforestry, the mechanisms underlying the genetic adaptation and rapid evolution of invading populations remain poorly understood. The black locust gall midge, Obolodiplosis robiniae, a highly invasive species that originated in North America, spread widely throughout Asia and Europe in the past decade. Here, we used 11 microsatellite DNA markers to analyze the genetic variation of 22 O. robiniae populations in China (the introduced region) and two additional US populations (the native region). A relatively high level of genetic diversity was detected among the introduced populations, even though they exhibited lower diversity than the native US populations. Evidence for genetic differentiation among the introduced Chinese populations was also found based on the high Fst value compared to the relatively low among the native US populations. Phylogenetic trees, structure graphical output, and principal coordinate analysis plots suggested that the Chinese O. robiniae populations (separated by up to 2,540 km) cluster into two main groups independent of geographical distance. Genetic variation has been observed to increase rapidly during adaptation to a new environment, possibly contributing to population establishment and spread. Our results provide insights into the genetic mechanisms underlying successful invasion, and identify factors that have contributed to colonization by an economically important pest species in China. In addition, the findings improve our understanding of the role that genetic structure plays during invasion by O. robiniae.

Comparison of different annotation tools for characterization of the complete chloroplast genome of Corylus avellana cv Tombul

BACKGROUND

Several bioinformatics tools have been designed for assembly and annotation of chloroplast (cp) genomes, making it difficult to decide which is most useful and applicable to a specific case. The increasing number of plant genomes provide an opportunity to accurately obtain cp genomes from whole genome shotgun (WGS) sequences. Due to the limited genetic information available for European hazelnut (Corylus avellana L.) and as part of a genome sequencing project, we analyzed the complete chloroplast genome of the cultivar 'Tombul' with multiple annotation tools.

RESULTS

Three different annotation strategies were tested, and the complete cp genome of C. avellana cv Tombul was constructed, which was 161,667 bp in length, and had a typical quadripartite structure. A large single copy (LSC) region of 90,198 bp and a small single copy (SSC) region of 18,733 bp were separated by a pair of inverted repeat (IR) regions of 26,368 bp. In total, 125 predicted functional genes were annotated, including 76 protein-coding, 25 tRNA, and 4 rRNA unique genes. Comparative genomics indicated that the cp genome sequences were relatively highly conserved in species belonging to the same order. However, there were still some variations, especially in intergenic regions, that could be used as molecular markers for analyses of phylogeny and plant identification. Simple sequence repeat (SSR) analysis showed that there were 83 SSRs in the cp genome of cv Tombul. Phylogenetic analysis suggested that C. avellana cv Tombul had a close affinity to the sister group of C. fargesii and C. chinensis, and then a closer evolutionary relationship with Betulaceae family than other species of Fagales.

CONCLUSION

In this study, the complete cp genome of Corylus avellana cv Tombul, the most widely cultivated variety in Turkey, was obtained and annotated, and additionally phylogenetic relationships were predicted among Fagales species. Our results suggest a very accurate assembly of chloroplast genome from next generation whole genome shotgun (WGS) sequences. Enhancement of taxon sampling in Corylus species provide genomic insights into phylogenetic analyses. The nucleotide sequences of cv Tombul cp genomes can provide comprehensive genetic insight into the evolution of genus Corylus.

Narrow gene pool can threaten the survival of Calamus nagbettai R. R. Fernald & Dey: a highly, endemic dioecious rattan species in the Western Ghats of India

Rattans, the spiny climbing palms of Arecaceae (Palmae) family exhibit high endemism to the biodiversity hot spots in India. Of the five rattan genera, Calamus is the only genus found in peninsular India with 15 of 21 species, endemic to the Western Ghats. The extensive utilization of rattans owing to their strength, durability and huge demand has resulted in depletion of their natural resources. Of the 15 endemic species, C. nagbettai is the most affected species on account of endemism, low population size and restricted distribution with fragmented populations. The present study revealed high amount of genetic diversity in the surviving scattered populations of the species using microsatellite markers. High gene flow (N_m = 1.498) observed across the populations resulted in low genetic differentiation (14%). A clear genetic admixture could be seen in Kerala as well as one of the Karnataka's populations while the remaining two populations were genetically distinct. UPGMA, PCoA and STRUCTURE analyses showed significantly different genetic composition in Kerala population compared to other populations. Kerala and Karnataka populations of C. nagbettai were also unique in their genetic structure and allelic composition. Therefore, effective management and conservation strategies have to be implemented to preserve the rarealleles with adaptive potential to protect this economically valuable Calamus species from endangerment. Over exploitation, low seed set and poor regeneration, as well as habitat fragmentation can further threaten the survival of this endemic, narrowly distributed dioecious rattan species in the Western Ghats region.

Genetic Diversity and Population Structure of Bermudagrass [Cynodon dactylon (L.) Pers.] along Latitudinal Gradients and the Relationship with Polyploidy Level

Understanding the population genetic pattern and process of gene flow requires a detailed knowledge of how landscape characteristics structure populations. Although Cynodon dactylon (L.) Pers. (common bermudagrass) is widely distributed in the world, information on its genetic pattern and population structure along latitudinal gradients is limited. We tried to estimate the genetic diversity and genetic structure of C. dactylon along a latitudinal gradient across China. Genetic diversity among different ploidy levels was also compared in the study. The material used consisted of 296 C. dactylon individuals sampled from 16 geographic sites from 22°35′ N to 36°18′ N. Genetic diversity was estimated using 153 expressed sequence tag-derived simple sequence repeat (EST-SSR) loci. Higher within-population genetic diversity appeared at low-latitude, as well as having positive correlation with temperature and precipitation. The genetic diversity increased with the ploidy level of C. dactylon, suggesting polyploidy creates higher genetic diversity. No isolation by distance and notable admixture structure existed among populations along latitudes. Both seed dispersal (or vegetative organs) and extrinsic pollen played important roles for gene flow in shaping the spatial admixture population structure of C. dactylon along latitudes. In addition, populations were separated into three clusters according to ploidy levels. C. dactylon has many such biological characters of perennial growth, wind-pollination, polyploidy, low genetic differentiation among populations, sexual and asexual reproduction leading to higher genetic diversity, which gives it strong adaptability with its genetic patterns being very complex across all the sampled latitudes. The findings of this study are related to landscape population evolution, polyploidy speciation, preservation, and use of bermudagrass breeding.

Genetic assessment of three Fagaceae species in forest restoration trials

Restoring isolated patches of forest ecosystems in degraded landscapes could potentially lead to genetic loss and inbreeding. Therefore, this study determined the occurrence of genetic diversity among the tree species Castanopsis tribuloides, C. calathiformis, and Lithocarpus polystachyus all of which were proven previously to be effective native tree species in the restoration of upland evergreen forests in northern Thailand when using the seed sample collection method. We tested our hypothesis as to whether the genetic diversity of a plant population that had been planted from the seeds of 4–6 adult trees would be lower and whether incidences of fixation index (Fis) would be higher among the second generation seedlings of these three Fagaceae species in isolated forest restoration trial plots. Microsatellite primers were selected from the entire genome sequence of C. tribuloides and the genetic sequences of C. tribuloides, L. polystachyus, and C. calathiformis were analyzed. Our results indicated a high degree of genetic diversity (He) in C. tribuloides (0.736) and C. calathiformis (0.481); however, a low level of genetic diversity was observed in L. polystachyus (0.281) within the restored forest. The fixation index for the second generation of L. polystachyus and C. calathiformis in the restored forest showed evidence of inbreeding. These results imply the efficiency of the seed sample collection method and verify that it does not reduce the level of genetic diversity in C. tribuloides and C. calathiformis. However, it may result in incidences of an inbreeding phenomena, suggesting the need to increase the number of adult trees used at the seed collection stage.

Genetic diversity and structure in the northern populations of European hazelnut (Corylus avellana L.)

European hazelnut (Corylus avellana L.) is a strictly cross-pollinated diploid tree species, which has its northernmost populations in Fennoscandia, and it was one of the first species to recolonize northern Europe after the last ice age. Hazelnut produces edible nuts in Finland but nowadays they are underutilized as food, and currently no breeding programmes exist. In the present study, 300 hazelnut specimens were collected from 20 different locations (= populations) in Finland, and they were genetically analyzed using nine simple sequence repeat (SSR) markers. Most of the genetic diversity existed within populations (83%). According to different genetic analyses (STRUCTURE, principal coordinates analysis, and clustering), a general lack of structure was observed, suggesting extensive gene flow among hazelnuts between 17 investigated populations. However, genetic structuring was clearly observed in three populations: Hakavuori, Mustiala, and Pähkinämäki, which might have become isolated due to geographical barriers that kept them separate, diminishing gene flow from other populations. Studying the diversity of European hazelnut is of great interest for understanding population genetics of a species distributed in its marginal areas in the north, and the results are also valuable for further uses in plant conservation, selection, and possible future breeding actions in Finland.

Assessment of Genetic Diversity and Population Genetic Structure of Norway Spruce (Picea abies (L.) Karsten) at Its Southern Lineage in Europe. Implications for Conservation of Forest Genetic Resources

In the present paper we studied the genetic diversity and genetic structure of five Norway spruce (Picea abies (L.) Karsten) natural populations situated in Serbia, belonging to the southern lineage of the species at the southern margin of the species distribution range. Four populations occur as disjunct populations on the outskirts of the Dinaric Alps mountain chain, whereas one is located at the edge of Balkan Mountain range and, therefore, can be considered as ecologically marginal due to drier climatic conditions occurring in this region. Due to the negative effect of biotic and abiotic stress factors, the sustainability of these populations is endangered, making conservation of their genetic resources one of the key measures of Norway spruce persistence in Serbia under climatic changes. The insight on genetic diversity and genetic structure of the studied spruce populations can provide the information required for the initiation of programs aimed at the conservation and utilization of spruce genetic resources at the rear edge of species environmental limits. Norway spruce genetic variation and population genetic structure were estimated using eight EST-SSR markers. The results showed that mean expected heterozygosity was 0.616 and allelic richness 10.22. Genetic differentiation among populations was low (Fst = 0.007). No recent bottleneck effect or isolation by distance were detected. Bayesian clustering, obtained with STRUCTURE, grouped the populations into two genetic clusters, whereas UPGMA analysis distinguished three main groups approximately in line with the geographic area of occurrence. Based on the study results and the EUFORGEN Pan-European strategy for genetic conservation of forest trees, the establishment of additional dynamic gene conservation units must be considered in Serbia in order to protect the adaptive and neutral genetic diversity of the species.

SSR identification and marker development for sago palm based on NGS genome data

Sago palm (Metroxylon sagu Rottb.) is one of the most productive carbohydrate-producing crops. Unfortunately, only limited information regarding sago palm genetics is available. This study aimed to develop simple sequence repeat (SSR) markers using sago palm NGS genomic data and use these markers to evaluate the genetic diversity of sago palm from Indonesia. De novo assembly of partial sago palm genomic data and subsequent SSR mining identified 29,953 contigs containing 31,659 perfect SSR loci and 31,578 contigs with 33,576 imperfect SSR loci. The perfect SSR loci density was 132.57/Mb, and AG, AAG and AAAT were the most frequent SSR motifs. Five hundred perfect SSR loci were randomly selected and used for designing SSR primers; 93 SSR primer pairs were identified. After synteny analysis using rice genome sequences, 20 primer pairs were validated using 11 sago palm accessions, and seven primers generated polymorphic alleles. Genetic diversity analysis of 41 sago palm accessions from across Indonesia using polymorphic SSR loci indicated the presence of three clusters. These results demonstrated the success of SSR identification and marker development for sago palm based on NGS genome data, which can be further used for assisting sago palm breeding in the future.

Next Generation Sequencing-Based Molecular Marker Development: A Case Study in Betula Alnoides

Betula alnoides is a fast-growing valuable indigenous tree species with multiple uses in the tropical and warm subtropical regions in South-East Asia and southern China. It has been proved to be tetraploid in most parts of its distribution in China. In the present study, next generation sequencing (NGS) technology was applied to develop numerous SSR markers for B. alnoides, and 64,376 contig sequences of 106,452 clean reads containing 164,357 candidate SSR loci were obtained. Among the derived SSR repeats, mono-nucleotide was the main type (77.05%), followed by di- (10.18%), tetra- (6.12%), tri- (3.56%), penta- (2.14%) and hexa-nucleotide (0.95%). The short nucleotide sequence repeats accounted for 90.79%. Among the 291 repeat motifs, AG/CT (46.33%) and AT/AT (44.15%) were the most common di-nucleotide repeats, while AAT/ATT (48.98%) was the most common tri-nucleotide repeats. A total of 2549 primer sets were designed from the identified putative SSR regions of which 900 were randomly selected for evaluation of amplification successfulness and detection of polymorphism if amplified successfully. Three hundred and ten polymorphic markers were obtained through testing with 24 individuals from B. alnoides natural forest in Jingxi County, Guangxi, China. The number of alleles (NA) of each marker ranged from 2 to 19 with a mean of 5.14. The observed (HO) and expected (HE) heterozygosities varied from 0.04 to 1.00 and 0.04 to 0.92 with their means being 0.64 and 0.57, respectively. Shannon-Wiener diversity index (I) ranged from 0.10 to 2.68 with a mean of 1.12. Cross-species transferability was further examined for 96 pairs of SSR primers randomly selected, and it was found that 48.96⁻84.38% of the primer pairs could successfully amplify each of six related Betula species. The obtained SSR markers can be used to study population genetics and molecular marker assisted breeding, particularly genome-wide association study of these species in the future.

Comparative Genomics and Phylogenetic Analysis Revealed the Chloroplast Genome Variation and Interspecific Relationships of Corylus (Betulaceae) Species

Corylus L. is an economically and phylogenetically important genus in the family Betulaceae. Taxonomic and phylogenetic relationships of Corylus species have long been controversial for lack of effective molecular markers. In this study, the complete chloroplast (cp) genomes of six Corylus species were assembled and characterized using next-generation sequencing. We compared the genome features, repeat sequences, sequence divergence, and constructed the phylogenetic relationships of the six Corylus species. The results indicated that Corylus cp genomes were typical of the standard double-stranded DNA molecule, ranging from 160,445 base pairs (bp) (C. ferox var. thibetca) to 161,621 bp (C. yunnanensis) in length. Each genome contained a pair of inverted repeats (IRs), a large single-copy (LSC) region and a small single-copy (SSC) region. Each of the six cp genomes possessed 113 unique genes arranged in the same order, including 80 protein-coding, 29 tRNA, and 4 rRNA genes. C. yunnanensis contained the highest number of repeat sequences, and the richest SSRs in six cp genomes were A/T mononucleotides. Comparative analyses of six Corylus cp genomes revealed four hotspot regions (trnH-psbA, rpoB-trnC, trnF-ndhJ, and rpl32-trnL) that could be used as potential molecular markers. Phylogenetic analyses of the complete chloroplast genomes and 80 protein-coding genes exhibited nearly identical topologies that strongly supported the monophyly of Corylus and simultaneously revealed the generic relationships among Betulaceae. The availability of these genomes can offer valuable genetic information for further taxonomy, phylogeny, and species delimitation in Corylus or even Betulaceae plants.

Genetic diversity and population structure of Chinese natural bermudagrass [Cynodon dactylon (L.) Pers.] germplasm based on SRAP markers

Bermudagrass [Cynodon dactylon (L.) Pers.], an important turfgrass used in public parks, home lawns, golf courses and sports fields, is widely distributed in China. In the present study, sequence-related amplified polymorphism (SRAP) markers were used to assess genetic diversity and population structure among 157 indigenous bermudagrass genotypes from 20 provinces in China. The application of 26 SRAP primer pairs produced 340 bands, of which 328 (96.58%) were polymorphic. The polymorphic information content (PIC) ranged from 0.36 to 0.49 with a mean of 0.44. Genetic distance coefficients among accessions ranged from 0.04 to 0.61, with an average of 0.32. The results of STRUCTURE analysis suggested that 157 bermudagrass accessions can be grouped into three subpopulations. Moreover, according to clustering based on the unweighted pair-group method of arithmetic averages (UPGMA), accessions were divided into three major clusters. The UPGMA dendrogram revealed that accessions from identical or adjacent areas were generally, but not entirely, clustered into the same cluster. Comparison of the UPGMA dendrogram and the Bayesian STRUCTURE analysis showed general agreement between the population subdivisions and the genetic relationships among accessions. Principal coordinate analysis (PCoA) with SRAP markers revealed a similar grouping of accessions to the UPGMA dendrogram and STRUCTUE analysis. Analysis of molecular variance (AMOVA) indicated that 18% of total molecular variance was attributed to diversity among subpopulations, while 82% of variance was associated with differences within subpopulations. Our study represents the most comprehensive investigation of the genetic diversity and population structure of bermudagrass in China to date, and provides valuable information for the germplasm collection, genetic improvement, and systematic utilization of bermudagrass.

Male and Female Subpopulations of Salix viminalis Present High Genetic Diversity and High Long-Term Migration Rates between Them

Dioecy distributed in 157 flowering plant families and 959 flowering plant genera. Morphological and physiological differences between male and female plants have been studied extensively, but studies of sex-specific genetic diversity are relatively scarce in dioecious plants. In this study, 20 SSR loci were employed to examine the genetic variance of male subpopulations and female subpopulations in Salix viminalis. The results showed that all of the markers were polymorphic (Na = 14.15, He = 0.7566) and workable to reveal the genetic diversity of S. viminalis. No statistically significant difference was detected between male and female subpopulations, but the average genetic diversity of male subpopulations (Na = 7.12, He = 0.7071) and female subpopulations (Na = 7.31, He = 0.7226) were high. Under unfavorable environments (West Liao basin), the genetic diversity between male and female subpopulations was still not significantly different, but the genetic diversity of sexual subpopulations were lower. The differentiation of the ten subpopulations in S. viminalis was moderate (FST = 0.0858), which was conformed by AMOVA that most of genetic variance (94%) existed within subpopulations. Pairwise FST indicated no differentiation between sexual subpopulations, which was accompanied by high long-term migrate between them (M = 0.73~1.26). However, little recent migration was found between sexual subpopulations. Therefore, artificial crossing or/and transplantation by cutting propagation should be carried out so as to increase the migration during the process of ex situ conservation.