Journal of Heredity Advance Access first published online on May 22, 2008
This version published online on July 22, 2008
Journal of Heredity, doi:10.1093/jhered/esn036
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
The First Sex-Specific Molecular Marker Discovered in the Moss Pseudocalliergon trifarium
From the Department of Applied Biology, University of Helsinki, PO Box 27, FI-00014 Helsinki, Finland (Korpelainen and Kolehmainen); and the Swedish Museum of Natural History, Department of Cryptogamic Botany, Box 50007, SE-104 05 Stockholm, Sweden (Bisang and Hedenäs)
Address correspondence to H. Korpelainen at the address above, or e-mail: helena.korpelainen{at}helsinki.fi.
Most dioecious plants do not exhibit discernible sexual dimorphism before sexual maturity. Therefore, it is impossible to address any sex-related questions during the prereproductive phase unless a genetic sex marker is available for gender determination. The aim of the present study was to develop a genetic sex marker for the moss Pseudocalliergon trifarium to allow gender and sex ratio determination at any stage in the life cycle. A high proportion of P. trifarium populations do not express sex. The screening of genomic DNA with inter simple sequence repeat (ISSR) primers was used to discover sex-specific polymerase chain reaction (PCR) amplification products. A presumably female-specific band was found, excised from the gel, cloned, and sequenced. A sequence-walking method was used to characterize the same region in males. A primer pair was designed to allow the amplification of a 159-bp portion of the female-specific DNA region. All tested material, up to 16-year-old herbarium specimens, provided unambiguous amplification products. This study successfully provides, for the first time in a moss, a sex-specific DNA marker. It allows reliable determination of gender and sex ratios. The short length of the amplification product is an advantage as satisfactory PCR products are more likely when the targeted sequence is short. The amount of variation in the DNA region shared by both sexes was relatively high. If the male sequence can be better characterized, the sex-specific regions could possibly be used to evaluate sex-specific phylogeographic patterns.
The genetic basis of sex determination in dioecious plants appears diverse (for "dioicy" vs. "dioecy," see e.g., Wyatt 1985), including heteromorphic sex chromosomes and also one or several autosomal nuclear sex-determining loci, possibly influenced by cytoplasmic genes and environmental cues (Korpelainen 1998; Ainsworth 2000). Previously, sex markers have been primarily developed for plant species that have been shown or proposed to possess sex chromosomes. They included sex-linked randomly amplified polymorphic DNA (RAPD)-, amplified fragment length polymorphism (AFLP)-, and RAPD-derived SCAR markers and simple sequence repeats in Actinidia chinensis (Harvey et al. 1997), Asparagus officinalis (Jiang and Sink 1997), Cannabis sativa (e.g., Sakamoto et al. 1995; Törjék et al. 2002; Rode et al. 2005), Carica papaya (Urasaki et al. 2002), Dioscorea tokoro (Terauchi and Kahl 1999), Ginkgo biloba (Jiang et al. 2003), Hippophae rhamnoides (Persson and Nybom 1998), Humulus lupulus (Polley et al. 1997), and Pistacia vera (Hormaza et al. 1994; Yakubov et al. 2005). Also, sex chromosome–specific DNA sequences have been characterized further and used as sex markers in Rumex acetosa (Shibata et al. 1999; Korpelainen 2002), Rumex nivalis (Stehlik and Blattner 2004), and Silene latifolia (e.g., Mulcahy et al. 1992; Lyons et al. 1995; Zhang et al. 1998). In bryophytes, characterized by a haploid-dominant life cycle, sex-specific markers have so far been discovered in the liverworts Marchantia polymorpha (Okada et al. 2000, 2001; Yamato et al. 2007) and Sphaerocarpos texanus (McLetchie and Collins 2001), both of which possess heteromorphic sex chromosomes. As to our knowledge, mosses have not been studied with respect to molecular sexual markers to date. However, Newton (1971) used cytological evidence for gender identification in Plagiomnium undulatum.
In general, sex markers have been developed only rarely for plants with presumably autosomal sex determination, including Atriplex garrettii (Ruas et al. 1998), Salix viminalis (Alström-Rapaport et al. 1998; Gunter et al. 2003), and Uapaca kirkiana (Mwase et al. 2007). It is clearly more likely to find sex-linked markers in cases where a sex chromosome system operates. Also, most of the known sex-linked markers are male associated, which reflects males being the heterogametic sex (Ainsworth 2000). Mainly, the genetic methods to identify gender have been developed for the purpose of genetic mapping or breeding work (Ainsworth 2000), and only rarely have such methods been applied to investigations on the sex ratios of natural populations (Lyons et al. 1995; Korpelainen 2002).
Most dioecious plants are not known to be sexually dimorphic before sexual maturity. Therefore, it is impossible to address any sex-related questions during the prereproductive phase unless a genetic sex marker is available for gender determination. Examples of investigations based on the utilization of genetic sex markers include a study by Korpelainen (2002), which showed that the sex ratios of R. acetosa populations are about 1:1 in the whole seed pool but that a significant female bias develops by the time of flowering. Comparably, Mulcahy et al. (1992) have examined sex-specific germination patterns in S. latifolia. Dioecious species with infrequent sex expression are even more challenging: The gender of most individuals cannot be determined at any stage in the life cycle, and sex-related issues are not possible to investigate in-depth in species with rare sexuality without a genetic sex marker.
A dioecious breeding system is widespread among bryophytes: More than half of all moss taxa and roughly two-thirds of the liverworts worldwide are unisexual (Wyatt 1982) Among these, there is a considerable number of species and many populations which consist of nonexpressing gametophytes only, bearing neither male nor female inflorescences (Longton and Schuster 1983; Bisang and Hedenäs 2005 and references therein). Biased expressed sex ratios are common, and females outnumber males in most species studied so far (reviewed by Bisang and Hedenäs 2005). The underlying causes for this, however, are poorly understood (Bisang et al. 2006).
The pleurocarpous P. trifarium (Web. & Mohr) Loeske is a dioecious moss of the order Hypnales. It is relatively common in the northern temperate zone and occurs also in the mountains of South America. The species mainly grows in deep fens or in sloping fens with slowly moving mineral-rich water (Hedenäs 1992; Hedenäs et al. 2003). A high proportion of P. trifarium populations are nonexpressing (67%), sporophyte production is rare (4.2%), and the expressed sex ratio is skewed toward the female gender (F:M = 2:7; N = 215; Bisang and Hedenäs, 2005; Bisang I, Hedenäs L, unpublished data). It is therefore a well-suited target species for the development of a molecular method to determine sex in nonexpressing moss individuals. Although detailed mechanisms for sex determination are not fully understood in many bryophyte species, there is evidence of sex chromosomes in a number of taxa (e.g., Ramsay and Berry 1982; McDaniel et al. 2006). Pseudocalliergon trifarium, however, has not been studied in this respect, and even its chromosome number is unknown (Fritsch 1991).
The aim of the present study is to develop a genetic sex marker for P. trifarium to allow gender and sex ratio determination at any stage in the life cycle and, eventually, the understanding of sex ratio dynamics in natural moss populations. In addition, we attempt to provide insights into the genetic basis of sex determination in P. trifarium.
| Materials and Methods |
|---|
|
|
|---|
Marker Development
The process for developing a sex-specific molecular marker was initiated with the screening of genomic DNA, extracted from dry material representing 6 female and 3 male samples of P. trifarium originating from Sweden (PT1–PT9; Table 1), with a set of ISSR primers (tandem arrays of simple nucleotide motifs interrupted in one end by out-of-phase bases). Due to a very small amount of sex-expressing plant material and no possibility to obtain sex-expressing progeny, there was no possibility to conduct a bulk or bulk segregant screening. The used screening method involved amplification by polymerase chain reaction (PCR), in separate reactions (no pooling of DNA or primers) with arbitrary ISSR primers, which had been assembled by the Nucleic Acid–Protein Service Unit, University of British Columbia, Canada. The primers used were the following: 807–812, 814–817, 819–826, 828, 834, 836, 843, 848–853, 858–862, 866, 881, 884, 887–889, and 891–892. The reaction mixture contained about 10–20 ng of genomic DNA, 1.2 units of DyNAzyme II DNA polymerase (Finnzymes, Espoo, Finland), 1x PCR buffer, 0.4 µl of 10 mM dNTP mix, and 1 µl of a single 5 µM primer. The thermocyler (MJ Research, Inc., model PTC-200) was programed for 4-min denaturation at 94 °C, followed by 45 cycles of denaturation at 94 °C for 45 s, annealing at 50–60 °C for 45 s, and elongation at 72 °C for 90 s. An additional 8-min elongation followed the last cycle. All amplification products were electrophoretically separated on 1.4% agarose gels. The only presumably sex-specific band resulting from the amplification with primer 807 was excised from the gel. The band was present in all females but in none of the males. The excised fragment representing sample PT1 was purified using an E.Z.N.A. Gel Extraction Kit (Omega Bio-tek Inc., Norcross, GA). Then, the purified DNA was concentrated to 3–5 µl and cloned using the TOPO TA Cloning kit by Invitrogen (Invitrogen, Carlsbad, CA). Positive clones were selected from the plates and transferred into 50 µl of ddH2O. DNA amplification was conducted in 20-µl reaction volumes containing 4 µl clone solution, 0.5 units of Phusion High-Fidelity DNA polymerase (Finnzymes), 1x PCR buffer, 0.4 µl of 10 mM dNTP mix, and 1 µl each of primers M13-f and M13-r (at the concentration of 5 µM) included in the cloning kit. The thermocyler was programed for 30 s denaturation at 98 °C, followed by 30 cycles of denaturation at 98 °C for 10 s, annealing at 55 °C for 20 s, and elongation at 72 °C for 50 s. An additional 8-min elongation followed the last cycle. The amplification product was separated electrophoretically, excised from the gel and purified. The following sequencing was conducted using the M13-f primer. The size of the sequenced female product was 670 bp (GenBank accession number EU368956; Figure 1). Specific primers were developed based on the cloned female sequence. Three primer pairs were tested, PT-1f (5'-TTC CTA GTG GGG AAC AGA AAA A-3') and PT-1r (5'-GAA TGA TGT TCT TGC ACG TCT C-3'), PT-1f and PT-2r (5'-TAA GCG GAA AAA TGG GAT TAG A-3'), and PT-3f (5'-GGA TTG ATA TTG GCA TTG AGT T-3') and PT-3r (5'-TGG AAT GTC ACA TTG TTT AGG A-3') with expected product sizes of 468 bp, 182 bp and 159 bp, respectively.
|
|
Marker Testing
The plant material used for testing the functioning of the sex marker included all available samples of sex-expressing P. trifarium specimen: 9 Swedish samples (6 females, 3 males) used in the screening for a sex-specific marker and 24 additional samples (PT10–PT33; 15 females, 9 males) originating from Sweden, Norway, Switzerland, and the United States (Table 1). When testing the functioning of the primer pairs, the following PCR conditions were used: The 20-µl reaction volumes contained about 10–20 ng of genomic DNA, 1.2 units of DyNAzyme II DNA polymerase (Finnzymes), 1x PCR buffer, 0.4 µl of 10 mM dNTP mix, and 1 µl of both 5 µM primers. The thermocyler was programed for denaturation at 94 °C for 4 min, followed by 30 cycles of denaturation at 94 °C for 45 s, annealing at 53–55 °C 45 s, and elongation at 72 °C for 60 s. An additional 8-min elongation followed the last cycle. All amplification products were electrophoretically separated on 1.4% agarose gels.
Characterization of Marker Region
To further characterize the sex-specific DNA region, genomic DNA of a female (sample PT23) and a male (sample PT24) were subjected to sequence walking using a method modified from Korpelainen et al. (2007) (Figure 2). Concentrated genomic DNA (about 15–20 ng/µl) was treated with two 4-cutter restriction enzymes leaving blunt ends after the protocol: restriction for 4 h at +37 °C utilizing HaeIII and RsaI restriction enzymes (both 4-base cutters) in separate reactions each in a volume of 8 µl containing 6 µl of genomic DNA, 4 units of either restriction enzyme, 1x restriction buffer, and 8 µg BSA. Then, a 2-stranded adaptor, composed of a 44-bp sequence (5'-GAA CTA GTC TCG ACT CCA GTC AGA GAT TCC ACC GCC GTG ACC GC-3') and an 8-bp sequence with 5'-phosphorylation and an amino modifier at the 3' end (5'-pGCG GTC AC-NH2-3'), was ligated to the restricted DNA. The ligation reaction was conducted for 2 h at +22 °C in a volume of 20 µl containing the previous 8-µl mix of restricted DNA, 0.5 µl of double-stranded 50 µM adaptor, 5 units of T4 ligase, 1x ligation buffer, 2 µl of 50% PEG 4000, and followed by inactivation at +65 °C for 15 min. After ligation, PCR amplifications were conducted with the primer PT-2r (found to target both male and female DNA template) and the adaptor primer ADAP3 (5'-TCG ACT CCA GTC AGA GAT TCC AC-3') and with the primer PT-2f (complementary to PT-2r) and the adaptor primer ADAP3 in 20-µl reaction volumes containing 2 µl restriction-ligation DNA, 0.5 units of DyNAzyme II Hot Start DNA polymerase (Finnzymes), 1x PCR buffer, 0.4 µl of 10 mM dNTP mix, and 1 µl of both 5 µM primers. The thermocyler was programed for 10-min denaturation at 94 °C, followed by 30 cycles of denaturation at 94 °C for 30 s, annealing at 59 °C for 25 s, and elongation at 72 °C for 40 s. An additional 8-min elongation followed the last cycle. The PCRs were followed by electrophoretic separation of the amplification products, fragment excision and purification, and sequencing with the adaptor primer. The resulting male or female sequence was presumed to represent the sequence either upstream or downstream from the primer region PT-2r, depending whether the primer PT-2f or PT-2r was used in combination with the primer ADAP3 (Figures 1 and 2).
|
| Results |
|---|
|
|
|---|
Out of more than 100 clear, scorable bands, only one band produced using the ISSR primer 807 showed a consistent sex-specific pattern present in all females but in none of the males. The fragment was cloned and sequenced, and primer pairs were designed. When testing the suitability of the 3 primer combinations to be used to amplify the female-specific marker, the first primer pair (PT-1f and PT-1r) resulted in a 468-bp amplification product only in females, as expected, but the amplification suffered from an uneven quality, varying from strong to weak. Surprisingly, the second primer pair combination (PT-1f and PT-2r) amplified well in both males and females (product size 182 bp). The third combination (PT-3f and PT-3r) amplified very well in females, but did not amplify in males (size 159 bp, Figure 3). Therefore, PT-3f and PT-3r were chosen to be the adequate primers for the sex marker of P. trifarium.
|
The clearly readable portions of the sequence-walked female (sample PT23) and male (sample PT24) sequences equalled 319 bp and 510 bp, respectively (GenBank accession numbers EU368957 [GenBank] and EU368958 [GenBank] ). When the female PT23 sequence was aligned with the cloned female PT1 sequence, the differences between the sequences included 17 substitutions (5.3%). On the other hand, when the male sequence PT24 was aligned with the cloned female sequence PT1, the first 290 bp aligned well, including 33 substitutions (11.4%), whereas the rest of the male sequence (220 bp) could not be aligned with the female sequence (Figure 1). When comparing the sequences PT23 and PT24 within the aligned 290-bp sequences, the number of substitutions equalled 16 (5.5%). All 3 sequenced samples represent different populations, with pairwise geographic distances equalling 24 km (PT1–PT24), 135 km (PT1–PT23), and 150 km (PT23–PT24).
A BLAST search (a nucleotide query vs. nucleotide database) of the cloned female-specific sequence and the sequence-walked male sequence against the sequences in GenBank failed to identify any match with the available sequence resources.
| Discussion |
|---|
|
|
|---|
This study successfully provides, for the first time in a moss, a sex-specific DNA marker. The marker is specific to females. The marker region was partially characterized in males as well, revealing the beginning of the DNA region unique to each sex. In addition, the ISSR primer 807, which resulted in the amplification of 3–5 clear bands, including the female-specific band (670 bp), can be used as a gender determination tool. However, the longer, specific primers PT-3f and PT-3r provide a single, unambiguous, easy to score amplification product. They are therefore the recommended way to conduct genetic gender determination in P. trifarium. When applying the sex marker to determine the gender of P. trifarium samples, it is important to use another PCR as an internal control for successful PCRs to avoid false negatives.
When the examined DNA is degraded, the short length (159 bp) of the female-specific amplification product is an advantage as obtaining satisfactory PCR products is more likely when the targeted sequence is short. The uneven quality of the larger female-specific product resulting from the amplification with the primers PT-1f and PT-1r may have been due to problems with degraded DNA in some herbarium specimens. In the tested material, amplifications with the primers PT-3f and PT-3r were very good in all samples, ranging from DNA extracted within weeks after field collection of the material to DNA extracted from herbarium specimens up to 16 years old.
Investigations of the sex ratios of plant populations have mostly relied on information obtained at sexual maturity. If uneven sex ratios are observed at maturity, it is of interest to reveal the processes leading to deviations from 1:1, which is the expected primary sex ratio at the end of parental investment (Fisher 1930). Biased expressed adult sex ratios may result from mechanisms acting at different ontogenetic stages, such as differences in life-history traits (e.g., germination rate, mortality, and vegetative vigor), sex-specific adult reproductive costs, or environmental responses (Bisang et al. 2006 and references therein). They may also be due to genetic mechanisms distorting the sex ratios (Taylor 1999). In organisms with infrequent sexuality, as in our study species P. trifarium, an observed skewed sex ratio might be the consequence of different degrees of sex expression (i.e., gametangia formation) between genetic sexes. Such species are a major challenge for studies of the sex ratio dynamics as the actual sex ratios cannot be properly revealed at any stage of life unless a molecular, sex-specific marker system is available.
It is not known whether the sequenced male and female regions of P. trifarium represent autosomal or partly homologous sex-chromosomal sequences. Because the males and females share a large part of the sequenced 670-bp female-specific product, it is possible that this region does not represent DNA in a heterogametic sex chromosome. However, most other plants, for which sex markers are available, have sex chromosomes, including 2 bryophyte species, M. polymorpha (Okada et al. 2000, 2001) and S. texanus (McLetchie and Collins 2001). On the other hand, it is conceivable that P. trifarium possesses heterogametic sex chromosomes and that the X and Y chromosomes contain partially shared sequences. Yamato et al. (2007) described the gene organization of the Y chromosome in M. polymorpha in detail. They identified 64 genes, of which 14 were detected only in the male genome and were expressed in the male reproductive organs suggesting their participation in male reproductive functions. At least 6 of the 40 male genes expressed both in thalli and male sexual organs have divergent X-linked female counterparts that are expressed in the thalli and sexual organs of female plants. These findings suggest that the Y and X chromosomes of M. polymorpha share the same ancestral autosome (Yamato et al. 2007).
Because the aligned part of the DNA region sequenced in P. trifarium was highly variable, it is likely to represent noncoding DNA which has accumulated mutations in both sexes. Previously, high mutation rates have been detected in Y-chromosomal heterochromatin in the males of R. acetosa (Navajas-Pérez et al. 2006; Korpelainen and Kostamo 2007). Evidently, the region sequenced in the moss P. trifarium, a part of which is the basis of the female-specific marker, is unlikely to represent a gene region involved in sex determination. Most likely, it is linked to a sex-determining locus, which is either autosomal or represents a partly homologous X- and Y-chromosomal regions. There may be one or, more likely, several loci which are involved in the sex determination of P. trifarium.
| Future Perspectives |
|---|
|
|
|---|
The characterization of a sex marker for P. trifarium opens up several fields of potential interest. We can now, for the first time, detect actual sex ratios in this species, from the species or wide geographical levels down to the population level, and compare these with the ratios of expressing plants. Gametangia formation in P. trifarium may not occur in all habitats where one sex grows. The new marker allows the investigation of expression rates in various habitats for each sex separately. Thus, it is possible to study both whether sex-specific habitat requirements occur or whether specific conditions are necessary for sex expression.
Depending on when the sex-specific regions evolved, a smaller or larger monophyletic group, including P. trifarium, could be expected to share the same sex marker. This should be tested at first for members of Pseudocalliergon and Drepanocladus (cf., Vanderpoorten et al. 2002, 2003). The amount of variation among the specimens in the region shared by both sexes was relatively high. If the male sequence can be better characterized, the specific regions could possibly be used to evaluate sex-specific phylogeographic patterns.
| Funding |
|---|
|
|
|---|
Swedish Research Council (621-2003-3338 to L.H.); Maj and Tor Nessling Foundation (2007012 to H.K.).
| Acknowledgments |
|---|
We thank H. Karttunen, K. Mirbakhsh, and M. Pietiläinen for help in the laboratory; C. Persson for field assistance; and Z. Schnyder and N. Schnyder for loan of specimens.
| Footnotes |
|---|
Corresponding Editor: Brian Murray
Received February 4, 2008
Accepted April 14, 2008
| References |
|---|
|
|
|---|
-
Ainsworth C. Boys and girls come out to play: the molecular biology of dioecious plants. Ann Bot (2000) 86:211–221.
Alström-Rapaport C, Lascoux M, Wang YC, Roberts G, Tuskan GA. Identification of a RAPD marker linked to sex determination in the basket willow (Salix viminalis L.). J Hered (1998) 89:44–49.
Bisang I, Ehrlén J, Hedenäs L. Reproductive effort and costs of reproduction do not explain female-biased sex ratios in the moss Pseudocalliergon trifarium (Amblystegiaceae). Am J Bot (2006) 93:1313–1319.
Bisang I, Hedenäs L. Sex ratio patterns in dioicous bryophytes re-visited. J Bryol (2005) 27:207–219.[CrossRef][Web of Science]
Fisher RA. The genetical theory of natural selection (1930) Oxford: Oxford University Press.
Fritsch R. Index to bryophyte chromosome counts. Bryophyt Bibl (1991) 40:1–352.
Gunter LE, Riberts GT, Lee K, Larimer FW, Tuskan GA. The development of two flanking SCAR markers linked to a sex determination locus in Salix viminalis L. J Hered (2003) 94:185–189.
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:149–154.[CrossRef]
Hedenäs L. The genus Pseudocalliergon in northern Europe. Lindbergia (1992) 16:80–99.[Web of Science]
Hedenäs L, Bisang I, Schnyder N. The distribution of bryophytes in Switzerland and in Liechtenstein IV. Hamatocaulis and Pseudocalliergon. Bot Helv (2003) 113:111–123.
Hormaza JI, Dollo L, Polito VS. Identification of a RAPD marker linked to sex determination in Pistacia vera using bulk segregant analysis. Theor Appl Genet (1994) 89:9–13.[Web of Science]
Jiang C, Sink KC. RAPD and SCAR markers linked to the sex expression locus M in asparagus. Euphytica (1997) 94:329–333.[CrossRef][Web of Science]
Jiang L, You RL, Li MX, Shi C. Identification of a sex-associated RAPD marker in Ginkgo biloba. Acta Bot Sin (2003) 45:742–747.
Korpelainen H. Labile sex expression in plants. Biol Rev (1998) 73:157–180.
Korpelainen H. A genetic method to resolve gender complement investigations on sex ratios in Rumex acetosa. Mol Ecol (2002) 11:2151–2156.[CrossRef][Medline]
Korpelainen H, Kostamo K. Excessive variation in Y-chromosomal DNA in Rumex acetosa (Polygonaceae). Plant Biol (2007) 9:383–389.[CrossRef][Medline]
Korpelainen H, Kostamo K, Virtanen V. Microsatellite marker identification using genome screening and restriction-ligation. Biotechniques (2007) 42:479–486.[Web of Science][Medline]
Longton RE, Schuster RM. Reproductive biology. In: New manual of bryology—Schuster RM, ed. (1983) Nichinan (Japan): The Hattori Botanical Laboratory. 386–462.
Lyons EE, Shah-Mahoney N, Lombard LA. Evolutionary dynamics of sex ratio and gender dimorphism in Silene latifolia: II. Sex ratio and flowering status in a potentially male-biased population. J Hered (1995) 86:107–113.
McDaniel SF, Quatrano RS, Shaw J. Sex specific evolutionary processes in Ceratodon purpureus inferred from variation in the sex linked des6 gene family [abstract], Moss 2006. (2006) Available from: http://ucjeps.berkeley.edu/bryolab/Moss2006/abstracts.htm. Accessed February 3, 2008.
McLetchie DN, Collins AL. Identification of DNA regions specific to the X and Y chromosomes in Sphaerocarpos texanus. Bryologist (2001) 104:543–547.[CrossRef]
Mulcahy DL, Weeden NF, Kesseli R, Carroll SB. DNA probes for the Y-chromosome of Silene latifolia, a dioecious angiosperm. Sex Plant Reprod (1992) 5:86–88.[CrossRef]
Mwase WF, Lid E, Björnstad A, Stedje B, Kwapata MB, Bokosi JM. Application of amplified fragment length polymorphism (AFLPs) for detection of sex-specific markers in dioecious Uapaca kirkiana Muell. Årg. Afr J Biotechnol (2007) 6:137–142.
Navajas-Pérez R, Schwarzacher T, de la Herrán R, Ruiz Rejón C, Ruiz Rejón M, Garrido-Ramos MA. The origin and evolution of the variability in a Y-specific satellite-DNA of Rumex acetosa and its relatives. Gene (2006) 368:61–71.[CrossRef][Web of Science][Medline]
Newton ME. A cytological distinction between male and female Mnium undulatum Hedw. Trans Br Bryol Soc (1971) 6:230–243.
Okada S, Fujisawa M, Sone T, Nakayama S, Nishiyama R, Takenaka M, Yamaoka S, Sakaida M, Kono K, Takahama M, et al. Construction of male and female PAC genomic libraries suitable for identification of Y-chromosome-specific clones from the liverwort, Marchantia polymorpha. Plant J (2000) 24:421–428.[CrossRef][Web of Science][Medline]
Okada S, Sone T, Fujisawa M, Nakayama S, Takenaka M, Ishizaki K, Kono K, Shimizu-Ueda Y, Hanajiri T, Yamato KT, et al. The Y chromosome in the liverwort Marchantia polymorpha has accumulated unique repeat sequences harboring a male-specific gene. Proc Natl Acad Sci USA (2001) 98:9454–9459.
Persson HA, Nybom H. Genetic sex determination and RAPD marker segregation in the dioecious species sea buckthorn (Hippophae rhamnoides L.). Hereditas (1998) 129:45–51.[CrossRef][Web of Science]
Polley A, Seigner E, Ganal MW. Identification of sex in hop (Humulus lupulus) using molecular markers. Genome (1997) 40:357–361.[Medline]
Ramsay HP, Berry GK. Sex determination in bryophytes. J Hattori Bot Lab (1982) 52:255–274.
Rode J, In-Chol K, Saal B, Flachowsky H, Kriese U, Weber WE. Sex-linked SSR markers in hemp. Plant Breed (2005) 124:167–170.[CrossRef]
Ruas CF, Fairbanks DJ, Evans RP, Stutz HC, Andersen WR, Ruas PM. Male-specific DNA in the dioecious species Atriplex garrettii (Chenopodiaceae). Am J Bot (1998) 85:162–167.[Abstract]
Sakamoto K, Shimomura K, Komeda Y, Kamada H, Satoh S. A male-associated DNA sequence in a dioecious plant, Cannabis sativa L. Plant Cell Physiol (1995) 36:1549–1554.
Shibata F, Hizume M, Kuroki Y. Chromosome painting of Y chromosomes and isolation of a Y-chromosome-specific repetitive sequence in the dioecious plant Rumex acetosa. Chromosome (1999) 108:266–270.[CrossRef]
Stehlik I, Blattner FR. Sex-specific SCAR markers in the dioecious plant Rumex nivalis (Polygonaceae) and implications for the evolution of sex chromosomes. Theor Appl Genet (2004) 108:238–242.[CrossRef][Web of Science][Medline]
Taylor DR. Genetics of sex ratio variation among natural populations of a dioecious plant. Evolution (1999) 53:55–62.[CrossRef][Web of Science]
Terauchi R, Kahl G. Mapping of the Dioscorea tokoro genome: AFLP markers linked to sex. Genome (1999) 97:1083–1089.
Törjék O, Bucherna N, Kiss E, Homoki H. Novel male-specific markers (MADC5, MADC6) in hemp. Euphytica (2002) 127:209–218.[CrossRef][Web of Science]
Urasaki N, Tokumoto M, Tarora K, Ban Y, Kayano T, Tanaka H, Oku H, Chinen I, Terauchi R. A male and hermaphrodite specific RAPD marker for papaya (Carica papaya L.). Theor Appl Genet (2002) 104:281–285.[CrossRef][Web of Science][Medline]
Vanderpoorten A, Goffinet B, Hedenäs L, Cox CJ, Shaw AJ. A taxonomic reassessment of the Vittiaceae (Hypnales, Bryopsida): evidence from phylogenetic analyses of combined chloroplast and nuclear sequence data. Plant Syst Evol (2003) 241:1–12.[CrossRef]
Vanderpoorten A, Hedenäs L, Cox C, Shaw AJ. Phylogeny and morphological evolution of the Amblystegiaceae (Bryopsida). Mol Phylogenet Evol (2002) 23:1–21.[CrossRef][Web of Science][Medline]
Wyatt R. Population ecology of bryophytes. J Hattori Bot Lab (1982) 52:179–198.
Wyatt R. Terminology for bryophyte sexuality: toward an unified system. Taxon (1985) 34:420–425.[CrossRef][Web of Science]
Yakubov B, Barazani O, Golan-Goldhirsh A. Combination of SCAR primers and touchdown-PCR for sex identification in Pistacia vera L. Sci Hortic (2005) 103:473–478.[CrossRef]
Yamato KT, Ishizaki K, Fujisawa M, Okada S, Nakayama S, Fujishita M, Bando H, Yodoya K, Hayashi K, Bando T, et al. Gene organization of the liverwort Y chromosome revelas distinct sex chromosome evolution in a haploid system. Proc Natl Acad Sci USA (2007) 104:6472–6477.
Zhang YH, Di Stilio VS, Rehman F, Avery A, Mulcahy D, Kesseli R. Y chromosome specific markers and the evolution of dioecy in the genus Silene. Genome (1998) 41:141–147.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||


