A gene-rich linkage map in the dioecious species Actinidia chinensis (kiwifruit) reveals putative X/Y sex-determining chromosomes
BMC Genomics volume 10, Article number: 102 (2009)
The genus Actinidia (kiwifruit) consists of woody, scrambling vines, native to China, and only recently propagated as a commercial crop. All species described are dioecious, but the genetic mechanism for sex-determination is unknown, as is the genetic basis for many of the cluster of characteristics making up the unique fruit. It is, however, an important crop in the New Zealand economy, and a classical breeding program would benefit greatly by knowledge of the trait alleles carried by both female and male parents. The application of marker assisted selection (MAS) in seedling populations would also aid the accurate and efficient development of novel fruit types for the market.
Gene-rich female, male and consensus linkage maps of the diploid species A. chinensis have been constructed with 644 microsatellite markers. The maps consist of twenty-nine linkage groups corresponding to the haploid number n = 29. We found that sex-linked sequence characterized amplified region (SCAR) markers and the 'Flower-sex' phenotype consistently mapped to a single linkage group, in a subtelomeric region, in a section of inconsistent marker order. The region also contained markers of expressed genes, some of unknown function. Recombination, assessed by allelic distribution and marker order stability, was, in the remainder of the linkage group, in accordance with other linkage groups. Fully informative markers to other genes in this linkage group identified the comparative linkage group in the female map, where recombination ratios determining marker order were similar to the autosomes.
We have created genetic linkage maps that define the 29 linkage groups of the haploid genome, and have revealed the position and extent of the sex-determining locus in A. chinensis. As all Actinidia species are dioecious, we suggest that the sex-determining loci of other Actinidia species will be similar to that region defined in our maps. As the extent of the non-recombining region is limited, our result supports the suggestion that the subtelomeric region of an autosome is in the early stages of developing the characteristics of a sex chromosome. The maps provide a reference of genetic information in Actinidia for use in genetic analysis and breeding programs.
New Zealand has a long history of interest in the genus Actinidia, being the country which commercialized the fruit of Actinidia deliciosa (A. Chev.) C.F. Liang et A.R. Ferguson var. deliciosa, as kiwifruit, and which recently released Actinidia chinensis Planch. var. chinensis 'Hort16A', the gold-fleshed kiwifruit, as an alternative cultivar. While New Zealand was instrumental in bringing these fruits to commercial attention, the genus is native to China and neighbouring countries where more than 60 species are known. This germplasm is relatively unexplored in terms of horticultural development of new and novel cultivars and offers a huge range of fruit characters and 'eating attributes', and plants suited to a wide range of climatic conditions. The diversity of flavours, fragrances, colours, and health factors are also of interest in genomic studies, offering the possibility of defining chemical pathways and identifying gene function.
Actinidia species present challenges to research and breeding. All known species in the genus are dioecious. Female plants bear flowers that are hermaphroditic in appearance but produce only empty pollen grains, while male plants have flowers that are unisexual with numerous stamens surrounding a rudimentary pistil whose growth is suppressed before style elongation or ovule initiation. Full dioecism is shown by about 4% of seed plants, and a second group display a variety of sub-dioecious conditions . Genetic studies have shown that dioecy has evolved many times in plants, and have demonstrated a variety of sex-determining systems [2, 3]. In Actinidia, bulk segregant analysis with random amplified polymorphic DNA (RAPD) markers supported the hypothesis that sex-determining genes were localized in a pair of chromosomes that function like an XX/XY system with male heterogamety [4–6]. The small size (<1 μm) of the chromosomes has made cytological studies difficult with the techniques available, and sex-determining chromosomes have not been positively identified. He et al.  using an improved chromosome binding technique, analyzed the karyotypes of diploid A. chinensis at the primary differentiation stage and reported that the sex chromosomes could not be identified from karyotypes of somatic cells. However, when they examined the pachytene stage of pollen mother cell meiosis, all 29 pairs of homologous chromosomes of pistillate and staminate plants paired tightly, except for a pair of nucleolar (SAT-) chromosomes in staminate plants. The two SAT-chromosomes were similar in length and shape, but in staminate plants the SAT region, about 15% of the total nucleolar chromosomal length, did not pair. He et al.  suggested the SAT region of nucleolar chromosomes may be the region of sex determination. They also suggested that sex chromosomes were probably at an early stage of differentiation in Actinidia.
The DNA content of the 2C genome of A. chinensis measured by flow cytometry was reported to be 1.3 – 1.4 pg , which corresponds to about 1.3 × 109 bp per genome. The genus contains species that form a polyploid series from diploid to octoploid .
As kiwifruit is a relatively new crop, knowledge of its genetic make-up is limited, so the development of a comprehensive genetic map and the use of molecular markers have the potential to improve efficiency in breeding new cultivars. A map will also help to simplify genomic studies to identify and isolate genes. Genetic linkage maps based on the recombination values of molecular markers have been constructed in an increasing number of plants (tomato [10, 11], rice , barley, , lotus , Brassica , cotton , grape ) and are proving valuable tools for plant breeding. The construction of a genetic map in an obligate outbreeding species, such as A. chinensis, is more complex than one derived from inbred or homozygous parents. Maps in outbreeding species have been developed by utilising the two-way pseudo-testcross procedure [18, 19], where the mapping population is the F1 progeny of a cross between unrelated, highly heterozygous individuals. Constructing the linkage map is complicated, as the two-way pseudo-test cross may segregate for up to four alleles at any locus, with one or both parents heterozygous at any given locus. The linkage phase of the markers will often be unknown, and can be different for the two parents, which can lead to inaccuracies in the estimation of recombination frequencies [20, 21]. The recombination frequencies can, however, be separately estimated for each parent so that two maps are developed, and these maps can be integrated using markers that are heterozygous in both parents. Two low density linkage maps have been reported in Actinidia. Testolin et al.  used the progeny of an interspecific cross between A. chinensis and A. callosa to construct, at a LOD score ≥ 2.0, a female map of 203 loci over 38 linkage groups and a male map of 143 loci over 30 linkage groups.
While marker systems such as restriction fragment length polymorphisms (RFLPs), amplified fragment length polymorphisms (AFLPs), random amplified polymorphic DNA (RAPDs), or single nucleotide polymorphisms (SNPs) have been developed to facilitate genetic mapping and gene discovery, the marker system of choice in many plant species is microsatellites (simple sequence repeats or SSRs). Microsatellites are arrays of short tandem repeat motifs of 1 to 5 base pairs in length which are characterized by their abundance, their distribution in both non-coding and coding regions of eukaryotic genomes, reproducibility, Mendelian mode of inheritance and co-dominant nature . They are recognised as highly informative genetic markers because of their inherent variability. This hypervariability is due to the high mutation rate within the nucleotide sequences of the microsatellites, and increases with increasing number of tandem repeats. In humans, heterozygosities generally exceed 0.5 and range as high as 0.9, with as many as 50 alleles per locus , and mutation rates, though variable among loci, exceed rates for non-microsatellite loci by up to four orders of magnitude [25, 26]. Similar hypervariability within microsatellite loci has been reported for birds, insects and plants, and loci may be polymorphic even in species where low levels of genetic diversity make alternate marker systems less useful [27, 28].
The time-consuming and expensive process of developing enriched genomic libraries and the subsequent sequencing and seeking of the simple sequence repeats is now often replaced with data mining of expressed sequence tag (EST) libraries to give a rapid, efficient and low-cost alternative for identifying microsatellites in plant species. Microsatellites have been found to occur regularly in ESTs . The frequency of occurrence of microsatellites of suitable length (20 nucleotides or more) varied in five cereals examined from 1.5% for maize to 4.7% for rice . This percentage would be sufficient to yield numerous markers from plant species in which large numbers of ESTs have been developed. In Actinidia the frequency of occurrence and level of polymorphism of EST-derived di-nucleotide microsatellites were sampled and found to be numerous in both the 5' and 3' ends of the genes represented, and highly polymorphic (93.5%) in the mapping population .
The construction of a single map for a cross in an outbreeding species, rather than two separate maps for the parents, depends on the availability of markers that are heterozygous in both parents. These markers form allelic bridges . Dominant markers such as RAPDs or AFLPs are generally of very limited use in combining parental maps, therefore, when the homologous linkage groups of the parents of a mapping population are required to be integrated, co-dominant markers such as microsatellites or RFLPs are the markers of choice, and allow the construction of either separate parental maps, or an integrated map for the cross .
Here we present comprehensive genetic linkage maps of female and male informative markers mapped in a cross in the outbreeding species A. chinensis, and also an integrated map of the cross, achieved through the use of co-dominant microsatellite markers. The twenty nine linkage groups are defined, and the position of sex-determining loci identified. Genetic linkage maps in Actinidia have been developed to supply markers for breeding novel cultivars, to provide tools for comparative and quantitative trait mapping, and to investigate the evolution and function of genetic control mechanisms.
Linkage map construction
We have created three gene-rich genetic maps, female, male and consensus, identifying the 29 linkage groups of the haploid genome and incipient X and Y sex chromosomes (Figures 1, 2, 3, 4 and 5, Table 1). The female linkage map constructed at LOD 4 and higher, was composed of 464 markers clearly defining 29 linkage groups and covered 2266 cM in the Kosambi function (see Additional file 1). The male map, composed of 365 markers, was shorter than the female map at 2078 cM in length (see Additional file 2). Estimated genome lengths of the female and male parents, using Method 3 of Chakravarti, were 3090 and 2782 cM respectively. The method based on average marker spacing adjustment  gave a genome length of 2562 cM in the female and 2402 cM in the male.
The same genome lengths with the adjustment for chromosome ends as per Remington et al. were 2820 and 2518 cM. A statistical programme that assumes markers are randomly distributed gave an estimate of intra-marker distance. Markers were estimated to be within 10 cM of each other in over 96% and 94% of the female and male genomes respectively (see Additional file 3). These estimates could be somewhat biased due to the non-randomness of marker distribution as noted below. The estimates were based on the adjusted genome lengths. The same genome coverage estimates based on the ratio of observed to estimated genome lengths were 80% and 83% respectively.
In the construction of the male map, Linkage Groups LG11 and LG17 were formatted in two separate sections. However, when the consensus map was compiled, the relatedness of the sections was confirmed through fully informative markers associating them with the corresponding linkage group in the female map. The consensus map, of length 2341 cM, and composed of 636 markers, showed some regions of altered linear order of markers from the sex-related maps (Figures 1, 2, 3, 4 and 5).
The Chi-square value for goodness-of-fit for the female map was 46.5 with 28 df. The p-value for a larger Chi-square value than observed was 0.015 which indicated a statistically significant departure from a Poisson random process. Linkage groups 1, 10, 14, 18 and 28 had significantly higher number of markers than expected from a random distribution. Similarly, linkage groups 4, 16 and 26 had significantly fewer. For the male map the Chi-square value was 36.6 with 28 df and a p-value of 0.128. When related to randomness of marker distribution, this value was non-significant.
A total of 799 primer pairs defining potential microsatellite markers were trialled for map construction. Polymorphism was established between the parents before a marker was evaluated for segregation across the genotypes of 272 siblings in the F1 mapping population. A total of 793 SSR markers from the 799 primer pairs tested, the two sex-linked SCAR markers and eight indel/SNP markers were considered to be of sufficient quality to use for map construction. Flower sex phenotype (FlowerSex), with aa alleles in females and ab alleles in males, was mapped as a male informative marker. A number of the markers amplified more than one locus in the mapping population and these are identified by a number in brackets (Figures 1, 2, 3, 4 and 5). Markers that were non-polymorphic in the parents, or non-segregating in the progeny, together made up 39% of the discarded potential markers, 18% of the marker results were difficult to read when one amplifying locus was overlying another and the alleles could not be unequivocally allocated to an individual locus, and 43% showed low information content, non-conforming ratios or poor PCR amplification. In some instances the PCR product was larger than the standard, and unable to be read. This was due to introns amplified from the genomic DNA that were not seen in the EST sequences. The allelic information content of the markers was of four types (Table 2). On the consensus map (Figures 1, 2, 3, 4 and 5) fully informative markers segregated 1:1:1:1 and made up 34% of the markers mapped, female informative and male informative markers segregating 1:1 were 37% and 21% respectively, and 8% were only partly informative segregating 1:2:1 or 3:1. Null alleles featured in all informative groups. Of the 636 markers in the consensus map, 587 markers were EST-derived and representative of expressed genes.
Identification of the X and Y chromosomes in the female and male linkage maps
The male-sex-linked marker SmY and the phenotype 'FlowerSex' both mapped to a subtelometic region on LG17 in the male genetic map. The fully informative marker udkac096, which mapped in the vicinity of the two markers, was used to identify the corresponding sex linkage group in the female map. While the markers over a large part of the length of these linkage groups were reliably ordered, in a section around the sex markers, some markers were found to alter order after consecutive runs of the map. This unreliably ordered portion of the linkage group contained seven markers including five EST-derived SSRs indicative of expressed genes (Table 3). Primer pairs of these five markers amplified products from both X and Y chromosomes, suggesting conservation of DNA sequence at primer sites on X and Y chromosomes.
Genes mapped to the sex chromosomes
Seventeen markers that mapped to the linkage group containing the sex-determining locus were derived from SSRs associated with ESTs and therefore could be matched to expressed genes. BLAST searches of the GenBANK database were used to try and identify potential functions for these ESTs (Table 3). Nine ESTs showed homology to sequences of known function. However, with the exception of Ke587 and Ke470, with homology to transcription factors, the function of these ESTs appeared unlikely to influence flower sex. The seven remaining genes were all listed in the 'unknown function' category.
Genetic linkage maps
Gene-rich genetic maps of the Actinidia chinensis female and male genomes, together with a consensus map of the two genomes, have been constructed. EST-derived microsatellites proved to be extremely polymorphic and a high proportion were mapped (Table 1). The loci were widely distributed over the 29 linkage groups we have identified. Some linkage groups carried fewer markers than others (Figures 1, 2, 3, 4 and 5) and this may be a reflection of the number of genes contained on a particular linkage group, or may be due to sampling bias in the ESTs selected for mapping. We sampled expressed genes that will be found in euchromatic rather than heterochromatic regions of the genome, not all tissues of the plant were represented, for example, there were no root libraries sampled, and we only selected ESTs that contained a large microsatellite in the transcribed sequence. While our sampling methods will have influenced the non-randomness of the markers on our genetic map, the abundance and relative distribution of microsatellites between transcribed and non-transcribed regions of the genome has been reviewed [33, 34], and has been reported to be non-random. Morgante et al.  report that microsatellite frequency was higher in transcribed regions, particularly in the untranslated portions, than in genomic DNA, and suggest that most microsatellites are found in regions pre-dating the genome expansion in many plants.
Many of the EST-derived primer pairs amplified more than one locus. The loci from a single primer pair could be on separate linkage groups, at a distance from one another on the same linkage group, or, with the resolving power of JoinMap software, unable to be separated on the map. Various authors [7, 22, 35, 36], have considered the possibility that diploid A. chinensis is a paleopolyploid, and some evidence they present would support this view. The haploid number of 29 is high, and would suggest that polyploidization may have occurred more than once, and may also have involved hybridization. Alternatively, duplication of a DNA segment may have occurred. This duplication event may have been physically separated by cross-over events or inversions, or tandem repeats may still be in evidence and markers that were unable to be separated could reflect this condition. As one of the primers of a pair was located within the translated portion of the EST to facilitate marker transfer to other species, it is reasonable to suggest that the genes were either duplicated, or members of a gene family with strong sequence homology in the transcribed region. However, a random priming event cannot be ruled out entirely.
Theoretical putative sex chromosomes
In the genus Actinidia all known species are dioecious, so it is reasonable to suggest that dioecy preceded speciation. For the sexes to remain separate over the period of differentiation to speciation, the genes responsible for the male and female characteristics would need to be tightly linked on the two haplotypes of one chromosome, and suppression of recombination would be essential to prevent the recurrence of hermaphroditism. It is not known how many genes are responsible for sex determination. In Actinidia we believe there must be at least two genes involved in the development of dioecy. One possible model that has been proposed has a dominant allele for pistil suppression closely linked to a dominant allele for pollen development on the putative Y chromosome, while the equivalent differential segment on the X chromosome has two alleles that function as recessives, one allowing pistil development and one leading to programmed pollen death [4, 37].
The genetic structure of the sex-determining region in the genus Actinidia has not previously been described. The chromosomes of Actinidia are small and of a fairly uniform size, and no definite sex chromosomes have been identified, though He et al.  described physical characteristics of a theoretical Y chromosome based on cytological studies. The genetic structure of LG 17 would also suggest that, in the male genotype, in a subtelomeric region of the chromosome, recombination was suppressed and marker order was difficult to establish. Lack of recombination is typical of a sex-determining region, so our data support the observation of He et al. .
There is evidence that sex chromosomes originate from autosomes [38–40], and it is thought that translocation of genetic material to sex chromosomes has occurred . A mutation that produced female or male sterility could, in theory, be found in any part of the chromosome. However, for full dioecy to develop, more than one mutation would need to have occurred, and the sexual differences would need to be fixed in the genome for gender-specific chromosomes to result. Suppression of recombination in the region of the mutations, allowing multiple loci to remain linked, would be required. Such regions of suppressed recombination are known, especially in the heterochromatic regions of the chromosomes such as the pericentromeric and subtelomeric regions of autosomes [10, 42]. Pericentromeric locations have been identified as incipient sex-determining loci in asparagus, Asparagus officinalis L. [43, 44], and papaya, Carica papaya , both species showing severe suppression of recombination around the sex-determining locus. Again, the recently reported incipient sex chromosomes in the genus Populus showed recombination suppression in the vicinity of the gender-linked locus. In Populus, like Actinidia chinensis, the sex-determining region was identified in the subtelomeric portion of a single chromosome pair . In the sex-determining locus of Actinidia several genes were associated with the sex-linked SmX and SmY markers and the phenotype 'FlowerSex'. The suppression of recombination around the sex-determining locus prevents recombination mapping from determining the accurate linear ordering of these genes, and the estimation of their distances from each other.
The genes of unknown function in the sex-determining locus may be similar to unidentified genes in other genera, or they may be specific to floral development in Actinidia. The markers have been used to identify bacterial artificial chromosomes (BACs) specific to the region and this is the first step in isolating and characterizing genes of the sex chromosomes in this genus. In addition, the gene-rich linkage map we have constructed will be a valuable resource for quantitative trait loci (QTL) analyses to identify markers related to traits of importance in breeding new and novel kiwifruits for the markets of the world. It will also considerably advance the development of a physical map for map-based cloning of genes for characterization.
We have described the genetic structure of all 29 linkage groups corresponding to n = 29 in a diploid species of the genus Actinidia. Map construction has been robust at LODs between 4 and 10. As the majority of the markers represent expressed genes, we anticipate that this resource will be useful in understanding genetic diversity in the genus. We have also identified, through sex-linked markers, putative X and Y chromosomes. These would appear to be in the early stages of evolution, as the subtelomeric sex-determining locus occupies only a small portion of the chromosome, as assessed through evidence of non-recombination, while the remainder of the chromosomes have the character of autosomes.
Plant material and DNA extraction
An intraspecific mapping population of 272 plants was created in the diploid species A. chinensis. Seedlings were screened with the sex marker SmY , to ensure equal numbers of female and male plants were propagated. The parents were chosen for their geographic separation in China, the female parent originating from seed from Henan province, Central China, and the male parent from a seed accession from Guangxi province, South China, and for the diversity of fruiting characters which they exhibited, those of the male being inferred from the attributes of female siblings. The mapping population was grown in the Plant and Food research orchard in Te Puke, Bay of Plenty, New Zealand. At budbreak, leaf tissue was taken from each genotype, held at 4° for 24 h, then stored at -80° until required. Seedlings were screened by flow cytometry to ascertain ploidy, and fingerprinted, using a kit of seven previously identified variable microsatellites which occurred in the parents. Only seedlings with the expected genotypes were included in the mapping analysis, all rogue plants were removed. DNA was extracted from young leaf tissue of both parental genotypes and every individual in the mapping population. A sample was ground to powder in liquid nitrogen before being processed through a DNeasy Plant Mini Kit (Qiagen™) according to the manufacturer's instructions. The final eluate was 200 μl in volume. 5 μl of a one in ten dilution of this eluate was used in each PCR reaction.
Microsatellite identification and primer design
Microsatellites (SSRs) suitable for use as markers were obtained from three sources. The first source was cDNA libraries constructed at Plant and Food Research from floral tissues of A. chinensis. The libraries were constructed in the pSport 1 (Not 1-Sal 1) vector and transformed into MAX Efficiency DH5α™ Competent Cells (Life Technologies). Clones were screened for the presence of microsatellite repeats by plaque hybridisation with short (16 bp), 32P-labelled probes and positive clones were identified. The second source of SSRs was two genomic libraries enriched either in (AC/GT)n or (AG/CT)n microsatellite repeats respectively which were constructed at the University of Udine, Italy , in Lambda Zap II vector (Stratagene). Transformation into XL1-Blue MRF Escherichia coli cells was followed by efficient excision of the plasmid from the Lambda Zap vector by the use of ExAssist® interference-resistant helper phage with SOLR™ cells (Stratagene). DNA from positive clones from both cDNA and genomic libraries was prepared, sequenced, and microsatellites identified. The third source of SSRs was from Plant and Food Research EST databases of A. chinensis and A. deliciosa sequences [31, 47]. Microsatellites from the EST libraries were identified in silico as described in Fraser et al. . Sequence data from this article have been deposited in the GenBank Data Libraries under accession nos. FG396279–FG528563.
Primer pairs were designed for non-duplicated sequences using the software programme Primer3 (©1996, 1997, 1998 ). Primer pair sequences were chosen which gave a theoretical PCR product size between 200 and 450 bp, with an annealing temperature between 55° and 60°, and with a GC content of approximately 50%. One of the primers of each pair was located before the microsatellite in the transcribed region, and the other was designed within the translated portion of the EST to facilitate marker transfer to other species. The primer pairs were synthesised and fluorescently-labelled (Dye Sets DS-31 or DS-34) by Applied Biosystems, Australia.
When a microsatellite was not found in the EST database, insertions and deletions (indels) and single nucleotide polymorphisms (SNPs) that were present were used to map particular genes of interest.
Polymerase chain reaction and electrophoresis
Primer pairs were screened for PCR amplification and length polymorphism with DNA samples of both parents of the mapping population, and the 272 progeny. A reaction mix of 15 μl containing 1 × PCR buffer (20 mM Tris-HCl, 50 mM KCl), MgCl2 5 mM (the buffer and MgCl2 were those supplied with the polymerase), 0.2 mM each of dNTPs, 4.5 pmol of each primer, and 1.25 units of Platinum Taq polymerase (Invitrogen), was prepared for each DNA sample. About 12.5 ng of genomic DNA was added in 5 μl to bring the total PCR volume to 20 μl. PCRs were performed in a Techne™ TC-412 thermal cycler with a single cycle of 94° for 3 min preceding 35 cycles of denaturing at 94° for 30 sec, annealing for 30 sec, and elongation at 72° for 1 min. PCR reactions were carried out individually before three colour multiplexes of products labelled with 6FAM, TET or HEX (Filter Set C), or 6FAM, VIC or NED (Filter Set D) were prepared for analysis. The allelic content of each genotype was determined by either gel electrophoresis in an ABI PRISM 377 DNA Sequencer (Filter Set C, TAMRA™ size standard), and analyzed with GeneScan Analysis and Genotyper software (Applied Biosystems), or capillary electrophoresis in an ABI Prism® 3100 Genetic Analyzer (Filter Set D, ROX™ size standard), and analyzed with GeneMapper™ Software Version 3.0 (Applied Biosystems). All markers were scored by at least two people independently for verification.
Data analysis and map construction
Chi-square tests of goodness-of-fit to expected segregation ratios of 1:1:1:1, 1:2:1, 3:1 or 1:1 were carried out for all markers segregating in the F1 progeny. The two sex-linked SCAR markers SmX and SmY previously developed in A. chinensis , were scored in the same fashion. The sex of the flowers of all genotypes in the mapping population was scored, and this phenotype was mapped as the 'FlowerSex' locus. Female, male and consensus linkage maps were constructed for the F1 progeny using JoinMap® 3.0 . For grouping of marker loci Joinmap uses a minimum LOD score only. A range of LOD scores, 4 – 10, were used so that linkage groups that were consistent across the range could be identified. In JoinMap, linkage is considered transitive, i.e. if A is linked to B which is linked to C, then A, B and C are linked. This gives one grouping of markers for each of the given threshold values. JoinMap uses only two-point (pairwise) analysis for ordering of markers. There are other software packages that use multi-point likelihood methods. One such package OutMap , which is also designed for out-crossing species, was used to validate a sample linkage group for consistency in marker order and distances. For ordering of markers within a linkage group, JoinMap uses a sequential build-up starting with the most informative, and 'rippling' for local optima. Map distances are estimated by minimising the sum of weighted (LOD scores as weights), squared difference of observed and expected pairwise distances .
The observed total map length was calculated simply as the sum of map distances of the terminal markers of each linkage group. Several methods were then used to provide estimates of the total genetic length of the genome. Firstly, the genome length was estimated using Method 3 of Chakravarti et al. , which is a modification of the method-of-moments estimator proposed by Hulbert et al. . This method estimates the genome length by multiplying, the ratio of total marker pairs to the number of marker pairs that equals or exceeds a specified LOD threshold value, by the map distance corresponding to the largest observed recombination fraction among the latter marker pairs. A threshold LOD value of 7 and the Kosambi genetic distance were used . The second method used was an adjustment to the Hulbert's method proposed by Remington et al.  to correct for the upward bias due to ends of linkage groups. The adjustment is given as
where n is the total number of markers, k is the number of marker pairs having a LOD equal to, or greater than, the specified threshold value, d is the map distance corresponding to the above LOD threshold and C is the haploid chromosome number. A third method used estimated individual linkage group length as the observed length plus twice the average marker spacing. This is based on the assumption that marker position is uniformly distributed, and the expected distance from a terminal marker to chromosome end is then equal to the average marker spacing . Adding these individual linkage group lengths gave an estimate of the total genome length.
Genome coverage was estimated first as the proportion of observed total map length to the estimated genome length. A second method used the following equation proposed by Lange and Boehnke , which makes the assumption markers are randomly distributed throughout the genome
c = 1 - e-2dn/L
where c is the proportion of genome within d cM of a marker, n and L are the number of markers and the estimated genome length.
Distribution of markers among linkage groups was assessed by comparing observed marker numbers in linkage groups with expectations under a random Poisson process. The expected number of markers for linkage group i under the null hypothesis is given by the Poisson parameter λ i = nL i /∑ L i , where L i s are lengths estimated by the third method described above. A chi-square value for goodness-of-fit was computed as ∑ [(O i - E i )2/E i ], where O and E are the observed and expected values, and the corresponding degrees of freedom is equal to one less than the number of linkage groups.
Lloyd DG, Webb CJ: Secondary sex characters in plants. Bot Rev. 1977, 43: 177-216. 10.1007/BF02860717.
Polley A, Seigner E, Ganal MW: Identification of sex in hop (Humulus lupulus) using molecular markers. Genome. 1997, 40 (3): 357-10.1139/g97-048.
Hamdi S, Teller G, Louis J-P: Master regulatory genes, auxin levels, and sexual organogeneses in the dioecious plant Mercurialis annua. Plant Physiol. 1987, 85: 393-399. 10.1104/pp.85.2.393.
Harvey CF, Gill GP, Fraser LG, McNeilage MA: Sex determination in Actinidia .1. Sex-linked markers and progeny sex ratio in diploid A. chinensis. Sex Plant Reprod. 1997, 10 (3): 149-154. 10.1007/s004970050082.
Gill GP, Harvey CF, Gardner RC, Fraser LG: Development of sex-linked PCR markers for gender identification in Actinidia. Theor Appl Genet. 1998, 97 (3): 439-10.1007/s001220050914.
Testolin R, Cipriani G, Costa G: Sex segregation ratio and gender expression in the genus Actinidia. Sex Plant Reprod. 1995, 8 (3): 129-10.1007/BF00242255.
He ZC, Huang HW, Zhong Y: Cytogenetic study of diploid Actinidia chinensis – Karyotype, morphology of sex chromosomes at primary differentiation stage and evolutionary significance. Acta Hort. 2003, 379-385. 610
Ferguson AR, O'Brien IEW, Yan GJ: Ploidy in Actinidia. Acta Hort. 1997, 444: 67-71.
Ferguson AR: 1904 – the year that kiwifruit (Actinidia deliciosa) came to New Zealand. N Z J Crop Hortic Sci. 2004, 32 (1): 3-
Tanksley SD, Ganal MW, Prince JP, Devicente MC, Bonierbale MW, Broun P, Fulton TM, Giovannoni JJ, Grandillo S, Martin GB, et al: High density molecular linkage maps of the tomato and potato genomes. Genetics. 1992, 132 (4): 1141-1160.
Chetelat RT, Meglic V, Cisneros P: A genetic map of tomato based on BC1 Lycopersicon esculentum × Solanum lycopersicoides reveals overall synteny but suppressed recombination between these homeologous genomes. Genetics. 2000, 154 (2): 857-
Harushima Y, Yano M, Shomura P, Sato M, Shimano T, Kuboki Y, Yamamoto T, Lin SY, Antonio BA, Parco A, et al: A high-density rice genetic linkage map with 2275 markers using a single F2 population. Genetics. 1998, 148 (1): 479-
Ramsay L, Macaulay M, Ivanissevich SD, MacLean K, Cardle L, Fuller J, Edwards KJ, Tuvesson S, Morgante M, Massari A, et al: A simple sequence repeat-based linkage map of barley. Genetics. 2000, 156 (4): 1997-
Sandal N, Krusell L, Radutoiu S, Olbryt M, Pedrosa A, Stracke S, Sato S, Kato T, Tabata S, Parniske M, et al: A genetic linkage map of the model legume Lotus japonicus and strategies for fast mapping of new loci. Genetics. 2002, 161 (4): 1673-
Kim JS, Chung TY, King GJ, Jin M, Yang TJ, Jin YM, Kim HI, Park BS: A sequence-tagged linkage map of Brassica rapa. Genetics. 2006, 174 (1): 29-10.1534/genetics.106.060152.
Guo WZ, Cai CP, Wang CB, Han ZG, Song XL, Wang K, Niu XW, Wang C, Lu KY, Shi B, et al: A microsatellite-based, gene-rich linkage map reveals genome structure, function and evolution in gossypium. Genetics. 2007, 176 (1): 527-10.1534/genetics.107.070375.
Troggio M, Malacarne G, Coppola G, Segala C, Cartwright DA, Pindo M, Stefanini M, Mank R, Moroldo M, Morgante M, et al: A dense single-nucleotide polymorphism-based genetic linkage map of grapevine (Vitis vinifera L.) anchoring pinot noir bacterial artificial chromosome contigs. Genetics. 2007, 176 (4): 2637-10.1534/genetics.106.067462.
Grattapaglia D, Sederoff R: Genetic linkage maps of Eucalyptus grandis and Eucalyptus urophylla using a pseudo-testcross: Mapping strategy and RAPD markers. Genetics. 1994, 137 (4): 1121-
Ritter E, Gebhardt C, Salamini F: Estimation of recombination frequencies and construction of RFLP linkage maps in plants from crosses between heterozygous parents. Genetics. 1990, 125 (3): 645-
Alm V, Fang C, Busso CS, Devos KM, Vollan K, Grieg Z, Rognli OA: A linkage map of meadow fescue (Festuca pratensis Huds.) and comparative mapping with other Poaceae species. Theor Appl Genet. 2003, 108 (1): 25-10.1007/s00122-003-1399-5.
Maliepaard C, Jansen J, Van Ooijen JW: Linkage analysis in a full-sib family of an outbreeding plant species: overview and consequences for applications. Genetical Research. 1997, 70 (3): 237-10.1017/S0016672397003005.
Testolin R, Huang WG, Lain O, Messina R, Vecchione A, Cipriani G: A kiwifruit (Actinidia spp.) linkage map based on microsatellites and integrated with AFLP markers. Theor Appl Genet. 2001, 103 (1): 30-10.1007/s00122-001-0555-z.
Scott KD, Eggler P, Seaton G, Rossetto M, Ablett EM, Lee LS, Henry RJ: Analysis of SSRs derived from grape ESTs. Theor Appl Genet. 2000, 100 (5): 723-10.1007/s001220051344.
Weber JL: Human DNA polymorphisms and methods of analysis. Curr Opin Biotechnol. 1990, 1: 166-171. 10.1016/0958-1669(90)90026-H.
Peakall R, Gilmore S, Keys W, Morgante M, Rafalski A: Cross-species amplification of soybean (Glycine max) simple sequence repeats (SSRs) within the genus and other legume genera: Implications for the transferability of SSRs in plants. Mol Biol Evol. 1998, 15 (10): 1275-
Weber JL, Wong C: Mutation of human short tandem repeats. Hum Mol Genet. 1993, 2 (8): 1123-10.1093/hmg/2.8.1123.
Hughes CR, Queller DC: Detection of highly polymorphic microsatellite loci in a species with little allozyme polymorphism. Mol Ecol. 1993, 2 (3): 131-10.1111/j.1365-294X.1993.tb00102.x.
Rongwen J, Akkaya MS, Bhagwat AA, Lavi U, Cregan PB: The use of microsatellite DNA markers for soybean genotype identification. Theor Appl Genet. 1995, 90 (1): 43-10.1007/BF00220994.
Varshney RK, Graner A, Sorrells ME: Genic microsatellite markers in plants: features and applications. Trends Biotechnol. 2005, 23 (1): 48-55. 10.1016/j.tibtech.2004.11.005.
Kantety RV, La Rota M, Matthews DE, Sorrells ME: Data mining for simple sequence repeats in expressed sequence tags from barley, maize, rice, sorghum and wheat. Plant Mol Biol. 2002, 48 (5): 501-10.1023/A:1014875206165.
Fraser LG, Harvey CF, Crowhurst RN, De Silva HN: EST-derived microsatellites from Actinidia species and their potential for mapping. Theor Appl Genet. 2004, 108 (6): 1010-10.1007/s00122-003-1517-4.
Remington DL, Whetten RW, Liu BH, O'Malley DM: Construction of an AFLP genetic map with nearly complete genome coverage in Pinus taeda. Theor Appl Genet. 1999, 98 (8): 1279-1292. 10.1007/s001220051194.
Li YC, Korol AB, Fahima T, Beiles A, Nevo E: Microsatellites: genomic distribution, putative functions and mutational mechanisms: a review. Mol Ecol. 2002, 11 (12): 2453-2465. 10.1046/j.1365-294X.2002.01643.x.
Morgante M, Hanafey M, Powell W: Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes. Nature Genet. 2002, 30 (2): 194-200. 10.1038/ng822.
Huang WG, Cipriani G, Morgante M, Testolin R: Microsatellite DNA in Actinidia chinensis : isolation, characterisation, and homology in related species. Theor Appl Genet. 1998, 97 (8): 1269-10.1007/s001220051019.
McNeilage MA, Considine JA: Chromosome studies in some Actinidia taxa and implications for breeding. N Z J Bot. 1989, 27 (1): 71-
McNeilage MA: Gender variation in Actinidia deliciosa, the kiwifruit. Sex Plant Reprod. 1991, 4 (4): 267-
Charlesworth B, Charlesworth D: Model for evolution of dioecy and gynodioecy. Am Nat. 1978, 112: 975-997. 10.1086/283342.
Charlesworth D, Charlesworth B, Marais G: Steps in the evolution of heteromorphic sex chromosomes. Heredity. 2005, 95 (2): 118-10.1038/sj.hdy.6800697.
Filatov DA: Evolutionary history of Silene latifolia sex chromosomes revealed by genetic mapping of four genes. Genetics. 2005, 170: 975-979. 10.1534/genetics.104.037069.
Ainsworth C: Boys and girls come out to play: The molecular biology of dioecious plants. Ann Bot. 2000, 86 (2): 211-10.1006/anbo.2000.1201.
Liu ZY, Moore PH, Ma H, Ackerman CM, Ragiba M, Yu QY, Pearl HM, Kim MS, Charlton JW, Stiles JI, et al: A primitive Y chromosome in papaya marks incipient sex chromosome evolution. Nature. 2004, 427 (6972): 348-10.1038/nature02228.
Telgmann-Rauber A, Jamsari A, Kinney MS, Pires JC, Jung C: Genetic and physical maps around the sex-determining M-locus of the dioecious plant asparagus. Mol Genet Genomics. 2007, 278 (3): 221-234. 10.1007/s00438-007-0235-z.
Jamsari A, Nitz I, Reamon-Buttner SM, Jung C: BAC-derived diagnostic markers for sex determination in asparagus. Theor Appl Genet. 2004, 108 (6): 1140-1146. 10.1007/s00122-003-1529-0.
Ma H, Moore PH, Liu ZY, Kim MS, Yu QY, Fitch MMM, Sekioka T, Paterson AH, Ming R: High-density linkage mapping revealed suppression of recombination at the sex determination locus in papaya. Genetics. 2004, 166 (1): 419-10.1534/genetics.166.1.419.
Yin T, DiFazio SP, Gunter LE, Zhang X, Sewell MM, Woolbright SA, Allan GJ, Kelleher CT, Douglas CJ, Wang M, et al: Genome structure and emerging evidence of an incipient sex chromosome in Populus. Genome Res. 2008, 18: 422-430. 10.1101/gr.7076308.
Crowhurst RN, Gleave AP, MacRae EA, Ampomah-Dwamena C, Atkinson RG, Beuning LL, Bulley SM, Chagne D, Marsh KB, Matich AJ, et al: Analysis of expressed sequence tags from Actinidia : applications of a cross species EST database for gene discovery in the areas of flavor, health, color and ripening. BMC Genomics. 2008, 9: 351-10.1186/1471-2164-9-351.
Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics methods and protocols: methods in molecular biology. Edited by: Krawetz S, Misener S. 2000, Humana: Totowa, 365-386.
Van Ooijen JW, Voorrips RE: JoinMap® 3.0, Software for the calculation of genetic linkage maps. 2001, Wageningen
Butcher PA, Williams ER, Whitaker D, Ling S, Speed TP, Moran GF: Improving linkage analysis in outcrossed forest trees – an example from Acacia mangium. Theor Appl Genet. 2002, 104 (6–7): 1185-1191.
Stam P: Construction of integrated genetic linkage maps by means of a new computer package: Joinmap. Plant J. 1993, 3 (5): 739-
Chakravarti A, Lasher LK, Reefer JE: A maximum likelihood method for estimating genome length using genetic linkage data. Genetics. 1991, 128 (1): 175-182.
Hulbert SH, Ilott TW, Legg EJ, Lincoln SE, Lander ES, Michelmore RW: Genetic analysis of the fungus, Bremia lactucae, using Restriction Fragment Length Polymorphisms. Genetics. 1988, 120 (4): 947-958.
Kosambi DD: The estimation of map distances from recombination values. Ann Eugen. 1944, 12: 172-175.
Lange K, Boehnke M: How many polymorphic genes will it take to span the human genome?. Am J Hum Genet. 1982, 34 (6): 842-845.
We thank Professor Raffaele Testolin and Dr Guido Cipriani of the University of Udine, Italy, for their valued collaboration during the early stages of library construction, and for sharing their enriched genomic libraries with us. We also thank Professor Richard Gardner, University of Auckland, for his input during the initial stages of the project, and Dr Elspeth McRae for continuing support. We acknowledge excellent technical help from students Erica Svenson, Hanna-Mari Tervo, Zhening Zhang, Reshma Pujari, Bhavneet Kaur, Natasha Austin, Heesuk Shin, and Brendan Swan. This work was partially funded by the New Zealand Foundation for Research, Science and Technology (CO6X0213) and by the HortResearch Kiwifruit Royalty Investment Programme.
LGF conceived the study, carried out SSR discovery, marker polymorphism discovery and population genotyping, and drafted the manuscript. GKT, CFH, and GPG performed SSR and marker polymorphism discovery and population genotyping, GKT also assisted in preparing the manuscript. PMD performed genotyping analyses and map construction. MAM recommended the most suitable mapping population and performed genotyping analyses. HNDS assisted with map construction, carried out all statistical analyses and wrote the statistical methods section for the manuscript. RNC developed the bioinformatics system used for SSR discovery. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Genetic linkage map (female) of Actinidia chinensis. The markers prefixed 'Ke' were from the kiwifruit EST database and represent expressed genes. Those prefixed 'udk' were from enriched genomic libraries, while all other prefixes relate to the bud libraries, and various markers as described in materials and methods. A number in brackets following a marker name indicates that a single primer pair amplified more than one locus. In the female map 29 linkage groups were defined. Incipient sex chromosomes were identified in Linkage Group 17 where the sex-determining locus was located in the subtelomeric region. (DOC 404 KB)
Additional file 2: Genetic linkage map (male) of Actinidia chinensis. The markers prefixed 'Ke' were from the kiwifruit EST database and represent expressed genes. Those prefixed 'udk' were from enriched genomic libraries, while all other prefixes relate to the bud libraries, and various markers as described in materials and methods. A number in brackets following a marker name indicates that a single primer pair amplified more than one locus. In the male map 29 linkage groups were defined. Incipient sex chromosomes were identified in Linkage Group 17 where the sex-determining locus was located in the subtelomeric region. (DOC 324 KB)
Additional file 3: Marker distribution in the female and male maps. A statistical programme that assumes markers are randomly distributed gave an estimate of intra-marker distance. Markers were estimated to be within 10 cM of each other in over 96% and 94% of the female and male genomes respectively. (DOC 64 KB)
About this article
Cite this article
Fraser, L.G., Tsang, G.K., Datson, P.M. et al. A gene-rich linkage map in the dioecious species Actinidia chinensis (kiwifruit) reveals putative X/Y sex-determining chromosomes. BMC Genomics 10, 102 (2009). https://doi.org/10.1186/1471-2164-10-102