Improving Nelumbo nucifera genome assemblies using high‐resolution genetic maps and BioNano genome mapping reveals ancient chromosome rearrangements

Genetic and physical maps are powerful tools to anchor fragmented draft genome assemblies generated from next‐generation sequencing. Currently, two draft assemblies of Nelumbo nucifera, the genomes of ‘China Antique’ and ‘Chinese Tai‐zi’, have been released. However, there is presently no information on how the sequences are assembled into chromosomes in N. nucifera. The lack of physical maps and inadequate resolution of available genetic maps hindered the assembly of N. nucifera chromosomes. Here, a linkage map of N. nucifera containing 2371 bin markers [217 577 single nucleotide polymorphisms (SNPs)] was constructed using restriction‐site associated DNA sequencing data of 181 F₂ individuals and validated by adding 197 simple sequence repeat (SSR) markers. Additionally, a BioNano optical map covering 86.20% of the ‘Chinese Tai‐zi’ genome was constructed. The draft assembly of ‘Chinese Tai‐zi’ was improved based on the BioNano optical map, showing an increase of the scaffold N50 from 0.989 to 1.48 Mb. Using a combination of multiple maps, 97.9% of the scaffolds in the ‘Chinese Tai‐zi’ draft assembly and 97.6% of the scaffolds in the ‘China Antique’ draft assembly were anchored into pseudo‐chromosomes, and the centromere regions along the pseudo‐chromosomes were identified. An evolutionary scenario was proposed to reach the modern N. nucifera karyotype from the seven ancestral eudicot chromosomes. The present study provides the highest‐resolution linkage map, the optical map and chromosome level genome assemblies for N. nucifera, which are valuable for the breeding and cultivation of N. nucifera and future studies of comparative and evolutionary genomics in angiosperms.     

The last bastion? X chromosome genotyping of Anopheles gambiae species pair males from a hybrid zone reveals complex recombination within the major candidate ‘genomic island of speciation’

Speciation with gene flow may be aided by reduced recombination helping to build linkage between genes involved in the early stages of reproductive isolation. Reduced recombination on chromosome X has been implicated in speciation within the Anopheles gambiae complex, species of which represent the major Afrotropical malaria vectors. The most recently diverged, morphologically indistinguishable, species pair, A. gambiae and Anopheles coluzzii, ubiquitously displays a ‘genomic island of divergence’ spanning over 4 Mb from chromosome X centromere, which represents a particularly promising candidate region for reproductive isolation genes, in addition to containing the diagnostic markers used to distinguish the species. Very low recombination makes the island intractable for experimental recombination studies, but an extreme hybrid zone in Guinea Bissau offers the opportunity for natural investigation of X‐island recombination. SNP analysis of chromosome X hemizygous males revealed: (i) strong divergence in the X‐island despite a lack of autosomal divergence; (ii) individuals with multiple‐recombinant genotypes, including likely double crossovers and localized gene conversion; (iii) recombination‐driven discontinuity both within and between the molecular species markers, suggesting that the utility of the diagnostics is undermined under high hybridization. The largely, but incompletely protected nature of the X centromeric genomic island is consistent with a primary candidate area for accumulation of adaptive variants driving speciation with gene flow, while permitting some selective shuffling and removal of genetic variation.     

A second‐generation genetic linkage map for bighead carp (Aristichthys nobilis) based on microsatellite markers

Bighead carp (Aristichthys nobilis) is an important aquaculture fish worldwide. Genetic linkage maps for the species were previously reported, but map resolution remained to be improved. In this study, a second‐generation genetic linkage map was constructed for bighead carp through a pseudo‐testcross strategy using interspecific hybrids between bighead carp and silver carp. Of the 754 microsatellites genotyped in two interspecific mapping families (with 77 progenies for each family), 659 markers were assigned to 24 linkage groups, which were equal to the chromosome numbers of the haploid genome. The consensus map spanned 1917.3 cM covering 92.8% of the estimated bighead carp genome with an average marker interval of 2.9 cM. The length of linkage groups ranged from 52.2 to 133.5 cM with an average of 79.9 cM. The number of markers per linkage group varied from 11 to 55 with an average of 27.5 per linkage group. Normality tests on interval distances of the map showed a non‐normal marker distribution; however, significant correlation was found between the length of linkage group and the number of markers below the 0.01 significance level (two‐tailed). The length of the female map was 1.12 times that of the male map, and the average recombination ratio of female to male was 1.10:1. Visual inspection showed that distorted markers gathered in some linkage groups and in certain regions of the male and female maps. This well‐defined genetic linkage map will provide a basic framework for further genome mapping of quantitative traits, comparative mapping and marker‐assisted breeding in bighead carp.     

Allotetraploidy of Zoysia species with 2n=40 based on a RFLP genetic map

 A RFLP linkage map of Zoysia spp. (2n=40), a warm-season turfgrass, was constructed by using the self-pollinated progenies obtained from an interspecific hybrid. Out of 115 DNA clones tested, 100 (87.0%), including 55 genomic clones, 38 cDNA clones, and seven gene clones encoding photosynthetic enzymes showed allelic-RFLP banding patterns among the parental accessions. We found that 26 probes detected two or more loci segregating in the self-pollinated progenies independently. The RFLP linkage map of Zoysia spp. consists of 115 loci in 22 linkage groups ranging in size from 12.5 cM to 141.3 cM with a total map distance of 1506 cM. Six RFLP loci (5.4%) showed significant segregation distortion (P<0.01). Two loci out of six were mapped to linkage group II, and another two loci were mapped to group VII. In the RFLP linkage map of zoysiagrass, five pairs of linkage groups sharing a series of duplicated loci with approximately the same order were identified. Therefore, we conclude that Zoysia spp. with 2n=40 should be considered as allotetraploids, which might have evolved from progenitors with a basic chromosome number of ten (x=10). Received: 20 March 1998 / Accepted: 17 September 1998     

Structural variation and rates of genome evolution in the grass family seen through comparison of sequences of genomes greatly differing in size

Homology was searched with genes annotated in the Aegilops tauschii pseudomolecules against genes annotated in the pseudomolecules of tetraploid wild emmer wheat, Brachypodium distachyon, sorghum and rice. Similar searches were performed with genes annotated in the rice pseudomolecules. Matrices of collinear genes and rearrangements in their order were constructed. Optical BioNano genome maps were constructed and used to validate rearrangements unique to the wild emmer and Ae. tauschii genomes. Most common rearrangements were short paracentric inversions and short intrachromosomal translocations. Intrachromosomal translocations outnumbered segmental intrachromosomal duplications. The densities of paracentric inversion lengths were approximated by exponential distributions in all six genomes. Densities of collinear genes along the Ae. tauschii chromosomes were highly correlated with meiotic recombination rates but those of rearrangements were not, suggesting different causes of the erosion of gene collinearity and evolution of major chromosome rearrangements. Frequent rearrangements sharing breakpoints suggested that chromosomes have been rearranged recurrently at some sites. The distal 4 Mb of the short arms of rice chromosomes Os11 and Os12 and corresponding regions in the sorghum, B. distachyon and Triticeae genomes contain clusters of interstitial translocations including from 1 to 7 collinear genes. The rates of acquisition of major rearrangements were greater in the large wild emmer wheat and Ae. tauschii genomes than in the lineage preceding their divergence or in the B. distachyon, rice and sorghum lineages. It is suggested that synergy between large quantities of dynamic transposable elements and annual growth habit have been the primary causes of the fast evolution of the Triticeae genomes.     

The high‐quality genome of Brassica napus cultivar ‘ZS11’ reveals the introgression history in semi‐winter morphotype

Allotetraploid oilseed rape (Brassica napus L.) is an agriculturally important crop. Cultivation and breeding of B. napus by humans has resulted in numerous genetically diverse morphotypes with optimized agronomic traits and ecophysiological adaptation. To further understand the genetic basis of diversification and adaptation, we report a draft genome of an Asian semi‐winter oilseed rape cultivar ‘ZS11’ and its comprehensive genomic comparison with the genomes of the winter‐type cultivar ‘Darmor‐bzh’ as well as two progenitors. The integrated BAC‐to‐BAC and whole‐genome shotgun sequencing strategies were effective in the assembly of repetitive regions (especially young long terminal repeats) and resulted in a high‐quality genome assembly of B. napus ‘ZS11’. Within a short evolutionary period (~6700 years ago), semi‐winter‐type ‘ZS11’ and the winter‐type ‘Darmor‐bzh’ maintained highly genomic collinearity. Even so, certain genetic differences were also detected in two morphotypes. Relative to ‘Darmor‐bzh’, both two subgenomes of ‘ZS11’ are closely related to its progenitors, and the ‘ZS11’ genome harbored several specific segmental homoeologous exchanges (HEs). Furthermore, the semi‐winter‐type ‘ZS11’ underwent potential genomic introgressions with B. rapa (A_r). Some of these genetic differences were associated with key agronomic traits. A key gene of A03.FLC3 regulating vernalization‐responsive flowering time in ‘ZS11’ was first experienced HE, and then underwent genomic introgression event with A_r, which potentially has led to genetic differences in controlling vernalization in the semi‐winter types. Our observations improved our understanding of the genetic diversity of different B. napus morphotypes and the cultivation history of semi‐winter oilseed rape in Asia.     

Identification of sex‐specific molecular markers using restriction site‐associated DNA sequencing

A major barrier to evolutionary studies of sex determination and sex chromosomes has been a lack of information on the types of sex‐determining mechanisms that occur among different species. This is particularly problematic in groups where most species lack visually heteromorphic sex chromosomes, such as fish, amphibians and reptiles, because cytogenetic analyses will fail to identify the sex chromosomes in these species. We describe the use of restriction site‐associated DNA (RAD) sequencing, or RAD‐seq, to identify sex‐specific molecular markers and subsequently determine whether a species has male or female heterogamety. To test the accuracy of this technique, we examined the lizard Anolis carolinensis. We performed RAD‐seq on seven male and ten female A. carolinensis and found one male‐specific molecular marker. Anolis carolinensis has previously been shown to possess male heterogamety and the recently published A. carolinensis genome facilitated the characterization of the sex‐specific RAD‐seq marker. We validated the male specificity of the new marker using PCR on additional individuals and also found that it is conserved in some other Anolis species. We discuss the utility of using RAD‐seq to identify sex‐determining mechanisms in other species with cryptic or homomorphic sex chromosomes and the implications for the evolution of male heterogamety in Anolis.     

A high-resolution comparative radiation hybrid map of ovine chromosomal regions that are homologous to human chromosome 6 (HSA6)

In this study, we constructed high-resolution radiation hybrid (RH) and comparative maps of ovine chromosomes or chromosomal segments that are homologous to human chromosome 6 (HSA6). A total of 251 markers were successfully genotyped across the recently developed USUoRH5000 whole-genome panel; 208 of these markers were assigned to five RH linkage groups distributed on three ovine chromosomes (OAR8, 9 and 20). The RH maps have good correspondence with previous chromosome painting data, although a small centromeric region on OAR9 that is homologous to HSA6 had not been previously detected using human chromosome paints on ovine chromosomal spreads. High percentages of the ovine markers were identified as orthologues in the bovine (86.3%), dog (85.8%), horse (69.3%) and human (88.7%) genomes. These maps contribute to investigations in mammalian chromosome evolution and the search for economic trait loci in sheep.     

An integrated SSR based linkage map for <Emphasis Type="Italic">Zoysia matrella</Emphasis> L. and <Emphasis Type="Italic">Z. japonica</Emphasis> Steud.

Japanese lawngrass (Zoysia japonica) and Manila grass (Z. matrella) are the two most important and commonly used Zoysia species. A consensus based SSR linkage map was developed for the genus by combining maps from each species. This used previously constructed maps for two Z. japonica populations and a new map from Z. matrella. The new SSR linkage map for Z. matrella was based on 86 F₂ individuals and contained 213 loci and covered a map distance of 1,351.2 cM in 32 linkage groups. Comparison of the three linkage maps constructed from populations with different genetic backgrounds indicated that most markers exhibited a consensus order, although some intervals or regions displayed discrepancy in marker orders or positions. The integrated map comprises 507 loci with a mean interval of 4.1 cM, covering a map distance of 2,066.6 cM in 22 linkage groups. The SSR-based map will allow marker-assisted selection and be useful for the mapping and cloning of economically important genes or quantitative trait loci.     

Construction of High-Density Linkage Maps of Populus deltoides × P. simonii Using Restriction-Site Associated DNA Sequencing

Although numerous linkage maps have been constructed in the genus Populus, they are typically sparse and thus have limited applications due to low throughput of traditional molecular markers. Restriction-site associated DNA sequencing (RADSeq) technology allows us to identify a large number of single nucleotide polymorphisms (SNP) across genomes of many individuals in a fast and cost-effective way, and makes it possible to construct high-density genetic linkage maps. We performed RADSeq for 299 progeny and their two parents in an F₁ hybrid population generated by crossing the female Populus deltoides ‘I-69’ and male Populus simonii ‘L3’. A total of 2,545 high quality SNP markers were obtained and two parent-specific linkage maps were constructed. The female genetic map contained 1601 SNPs and 20 linkage groups, spanning 4,249.12 cM of the genome with an average distance of 2.69 cM between adjacent markers, while the male map consisted of 940 SNPs and also 20 linkage groups with a total length of 3,816.24 cM and an average marker interval distance of 4.15 cM. Finally, our analysis revealed that synteny and collinearity are highly conserved between the parental linkage maps and the reference genome of P. trichocarpa. We demonstrated that RAD sequencing is a powerful technique capable of rapidly generating a large number of SNPs for constructing genetic maps in outbred forest trees. The high-quality linkage maps constructed here provided reliable genetic resources to facilitate locating quantitative trait loci (QTLs) that control growth and wood quality traits in the hybrid population.     

A sugar beet (Beta vulgaris L.) reference FISH karyotype for chromosome and chromosome‐arm identification,integration of genetic linkage groups and analysis of major repeat family distribution

We developed a reference karyotype for B. vulgaris which is applicable to all beet cultivars and provides a consistent numbering of chromosomes and genetic linkage groups. Linkage groups of sugar beet were assigned to physical chromosome arms by FISH (fluorescent in situ hybridization) using a set of 18 genetically anchored BAC (bacterial artificial chromosome) markers. Genetic maps of sugar beet were correlated to chromosome arms, and North–South orientation of linkage groups was established. The FISH karyotype provides a technical platform for genome studies and can be applied for numbering and identification of chromosomes in related wild beet species. The discrimination of all nine chromosomes by BAC probes enabled the study of chromosome‐specific distribution of the major repetitive components of sugar beet genome comprising pericentromeric, intercalary and subtelomeric satellites and 18S‐5.8S‐25S and 5S rRNA gene arrays. We developed a multicolor FISH procedure allowing the identification of all nine sugar beet chromosome pairs in a single hybridization using a pool of satellite DNA probes. Fiber‐FISH was applied to analyse five chromosome arms in which the furthermost genetic marker of the linkage group was mapped adjacently to terminal repetitive sequences on pachytene chromosomes. Only on two arms telomere arrays and the markers are physically linked, hence these linkage groups can be considered as terminally closed making the further identification of distal informative markers difficult. The results support genetic mapping by marker localization, the anchoring of contigs and scaffolds for the annotation of the sugar beet genome sequence and the analysis of the chromosomal distribution patterns of major families of repetitive DNA.     

家蚕细胞遗传学及其应用

由于家蚕染色体数目较多、着丝粒弥散, 在较长时期内, 家蚕染色体识别、核型分析、染色体结构和功能的研究都受到很大限制。近年来, 应用比较基因组杂交、基因组原位杂交、基于细菌人工染色体克隆的原位杂交技术建立了家蚕的细胞学图, 综合分子连锁图构建高密度的细胞遗传学图已成为可能。分子细胞遗传学的应用正在推动家蚕染色体结构和功能的研究, 揭示出家蚕W染色体密集地分布着嵌套结构的逆转座子, 染色体端粒由重复序列(TTAGG)n以及端粒特异的非长末端逆转座子TRAS1和SART1组成, TRAS1、SART1具有较高的转录活性, 可能与维持染色体的稳定性有关。     

Genetic Linkage Map Construction and QTL Mapping of Salt Tolerance Traits in Zoysiagrass (Zoysia japonica)

Zoysiagrass (Zoysia Willd.) is an important warm season turfgrass that is grown in many parts of the world. Salt tolerance is an important trait in zoysiagrass breeding programs. In this study, a genetic linkage map was constructed using sequence-related amplified polymorphism markers and random amplified polymorphic DNA markers based on an F₁ population comprising 120 progeny derived from a cross between Zoysia japonica Z105 (salt-tolerant accession) and Z061 (salt-sensitive accession). The linkage map covered 1211 cM with an average marker distance of 5.0 cM and contained 24 linkage groups with 242 marker loci (217 sequence-related amplified polymorphism markers and 25 random amplified polymorphic DNA markers). Quantitative trait loci affecting the salt tolerance of zoysiagrass were identified using the constructed genetic linkage map. Two significant quantitative trait loci (qLF-1 and qLF-2) for leaf firing percentage were detected; qLF-1 at 36.3 cM on linkage group LG4 with a logarithm of odds value of 3.27, which explained 13.1% of the total variation of leaf firing and qLF-2 at 42.3 cM on LG5 with a logarithm of odds value of 2.88, which explained 29.7% of the total variation of leaf firing. A significant quantitative trait locus (qSCW-1) for reduced percentage of dry shoot clipping weight was detected at 44.1 cM on LG5 with a logarithm of odds value of 4.0, which explained 65.6% of the total variation. This study provides important information for further functional analysis of salt-tolerance genes in zoysiagrass. Molecular markers linked with quantitative trait loci for salt tolerance will be useful in zoysiagrass breeding programs using marker-assisted selection.     

Cytogenetic mapping with centromeric bacterial artificial chromosomes contigs shows that this recombination‐poor region comprises more than half of barley chromosome 3H

Genetic maps are based on the frequency of recombination and often show different positions of molecular markers in comparison to physical maps, particularly in the centromere that is generally poor in meiotic recombinations. To decipher the position and order of DNA sequences genetically mapped to the centromere of barley (Hordeum vulgare) chromosome 3H, fluorescence in situ hybridization with mitotic metaphase and meiotic pachytene chromosomes was performed with 70 genomic single‐copy probes derived from 65 fingerprinted bacterial artificial chromosomes (BAC) contigs genetically assigned to this recombination cold spot. The total physical distribution of the centromeric 5.5 cM bin of 3H comprises 58% of the mitotic metaphase chromosome length. Mitotic and meiotic chromatin of this recombination‐poor region is preferentially marked by a heterochromatin‐typical histone mark (H3K9me2), while recombination enriched subterminal chromosome regions are enriched in euchromatin‐typical histone marks (H3K4me2, H3K4me3, H3K27me3) suggesting that the meiotic recombination rate could be influenced by the chromatin landscape.     

A transgenomic cytogenetic sorghum (Sorghum propinquum) bacterial artificial chromosome fluorescence in situ hybridization map of maize (Zea mays L.) pachytene chromosome 9, evidence for regions of genome hyperexpansion

A cytogenetic FISH map of maize pachytene-stage chromosome 9 was produced with 32 maize marker-selected sorghum BACs as probes. The genetically mapped markers used are distributed along the linkage maps at an average spacing of 5 cM. Each locus was mapped by means of multicolor direct FISH with a fluorescently labeled probe mix containing a whole-chromosome paint, a single sorghum BAC clone, and the centromeric sequence, CentC. A maize-chromosome-addition line of oat was used for bright unambiguous identification of the maize 9 fiber within pachytene chromosome spreads. The locations of the sorghum BAC-FISH signals were determined, and each new cytogenetic locus was assigned a centiMcClintock position on the short (9S) or long (9L) arm. Nearly all of the markers appeared in the same order on linkage and cytogenetic maps but at different relative positions on the two. The CentC FISH signal was localized between cdo17 (at 9L.03) and tda66 (at 9S.03). Several regions of genome hyperexpansion on maize chromosome 9 were found by comparative analysis of relative marker spacing in maize and sorghum. This transgenomic cytogenetic FISH map creates anchors between various maps of maize and sorghum and creates additional tools and information for understanding the structure and evolution of the maize genome.     

Developing a flexible,high‐efficiency Agrobacterium‐mediated sorghum transformation system with broad application

Sorghum is the fifth most widely planted cereal crop in the world and is commonly cultivated in arid and semi‐arid regions such as Africa. Despite its importance as a food source, sorghum genetic improvement through transgenic approaches has been limited because of an inefficient transformation system. Here, we report a ternary vector (also known as cohabitating vector) system using a recently described pVIR accessory plasmid that facilitates efficient Agrobacterium‐mediated transformation of sorghum. We report regeneration frequencies ranging from 6% to 29% in Tx430 using different selectable markers and single copy, backbone free ‘quality events’ ranging from 45% to 66% of the total events produced. Furthermore, we successfully applied this ternary system to develop transformation protocols for popular but recalcitrant African varieties including Macia, Malisor 84‐7 and Tegemeo. In addition, we report the use of this technology to develop the first stable CRISPR/Cas9‐mediated gene knockouts in Tx430.     

Recombination patterns reveal information about centromere location on linkage maps

Linkage mapping is often used to identify genes associated with phenotypic traits and for aiding genome assemblies. Still, many emerging maps do not locate centromeres – an essential component of the genomic landscape. Here, we demonstrate that for genomes with strong chiasma interference, approximate centromere placement is possible by phasing the same data used to generate linkage maps. Assuming one obligate crossover per chromosome arm, information about centromere location can be revealed by tracking the accumulated recombination frequency along linkage groups, similar to half‐tetrad analyses. We validate the method on a linkage map for sockeye salmon (Oncorhynchus nerka) with known centromeric regions. Further tests suggest that the method will work well in other salmonids and other eukaryotes. However, the method performed weakly when applied to a male linkage map (rainbow trout; O. mykiss) characterized by low and unevenly distributed recombination – a general feature of male meiosis in many species. Further, a high frequency of double crossovers along chromosome arms in barley reduced resolution for locating centromeric regions on most linkage groups. Despite these limitations, our method should work well for high‐density maps in species with strong recombination interference and will enrich many existing and future mapping resources.     

RAD‐seq linkage mapping and patterns of segregation distortion in sedges: meiosis as a driver of karyotypic evolution in organisms with holocentric chromosomes

Meiotic drive, the class of meiotic mechanisms that drive unequal segregation of alleles among gametes, may be an important force in karyotype evolution. Its role in holocentric organisms, whose chromosomes lack localized centromeres, is poorly understood. We crossed two individuals of Carex scoparia (Cyperaceae) with different chromosome numbers (2n = 33^II = 66 × 2n = 32^II = 64) to obtain F1 individuals, which we then self‐pollinated to obtain second‐generation (F2) crosses. RAD‐seq was performed for 191 individuals (including the parents, five F1 individuals and 184 F2 individuals). Our F2 linkage map based on stringent editing of the RAD‐seq data set yielded 32 linkage groups. In the final map, 865 loci were located on a linkage map of 3966.99 cM (linkage groups ranged from 24.39 to 193.31 cM in length and contained 5–51 loci each). Three linkage groups exhibit more loci under segregation distortion than expected by chance; within linkage groups, loci exhibiting segregation distortion are clustered. This finding implicates meiotic drive in the segregation of chromosome variants, suggesting that selection of chromosome variants in meiosis may contribute to the establishment and fixation of chromosome variants in Carex, which is renowned for high chromosomal and species diversity. This is an important finding as previous studies demonstrate that chromosome divergence may play a key role in differentiation and speciation in Carex.     

Isolation, characterization and mapping of simple sequence repeat markers in zoysiagrass (Zoysia spp.)

The genus Zoysia consists of 16 species that are naturally distributed on sea coasts and grasslands around the Pacific. Of these, Zoysia japonica, Zoysia matrella, and Zoysia tenuifolia are grown extensively as turfgrasses, and Z. japonica is also used as forage grass in Japan and other countries in East Asia. To develop simple sequence repeat (SSR) markers for zoysiagrass (Zoysia spp.), we used four SSR-enriched genomic libraries to isolate 1,163 unique SSR clones. All four libraries contained a high percentage of perfect clones, ranging from 67.1 to 96.0%, and compound clones occurred with higher frequencies in libraries A (28.6%) and D (11.6%). From these clones, we developed 1,044 SSR markers when we tested all 1,163 SSR primer pairs. Using all 1,044 SSR markers, we tested one screening panel consisting of eight Zoysia clones for testing PCR amplifications, from which five unrelated clones, among the eight, were used for polymorphism assessment, and found that the polymorphic information content ranged from 0 (monomorphic loci) to 0.88. Of the 1,044 SSR markers, 170 were segregated in our mapping population and we mapped 161 on existing amplified fragment length polymorphism-based linkage groups, using this mapping population. These SSR markers will provide an ideal marker system to assist with gene targeting, quantitative trait locus mapping, variety or species identification, and marker-assisted selection in Zoysia species.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.