Comparative analyses of six complete chloroplast genomes from the genus Cupressus and Juniperus (Cupressaceae).

Cupressaceae is a conifer family distributed around the world. Cupressus and Juniperus are the main genera of the Cupressaceae family and have important medicinal value. This leads to confusion between Cupressus and Juniperus due to similar morphologies. Here, the complete cp genomes of two Cupressus (C. duclouxiana and C. funebri) and four Juniperus (J. chinensis, J. gaussenii J. pingii and J. procumbens) were sequenced. The results revealed that the length of the cp genomes ranged from 126,996 bp to 129,959 bp, with 119 genes comprising 82 protein-coding genes, 33 transfer RNAs and 4 ribosomal RNAs. All chloroplast genomes of Cupressus and Juniperus lost whole IR regions, which is consistent with gymnosperm cp genome studies. In addition, the number of SSRs per species ranged from 54 to 73 and was dominated by mononucleotide repeats. In the six cp genomes of Cupressus and Juniperus, five highly divergent regions, including accD, accD-rpl2, ycf1, ycf2 and rrn23-rrn4.5, can be used as DNA barcodes of interspecific relationships and potential genetic markers. We compared the gene selection pressures (C. chengiana as reference species), and 6 genes underwent positive selection, the majority of which were related to photosynthesis. Phylogenetic results showed that the monophyly of Cupressus and Juniperus supported most bootstrap support. Cupressus funebris and J. chinensis were resolved to be early diverging species within Cupressus and Juniperus, and the two genera were sister groups to each other. This research revealed a new understanding of the structural pluralism and phylogenetic relationships of Cupressaceae cp genomes. These results will facilitate comprehension of the complexity and diversity of conifer cp genomes.

Symbiosis and pathogenicity of Geosmithia and Talaromyces spp. associated with the cypress bark beetles Phloeosinus spp. and their parasitoids.

Fungi associated with cypress bark beetles are practically unknown in the Eastern Mediterranean. Our study focused on the fungi associated with the body parts and galleries of two indigenous cypress bark beetles, Phloeosinus armatus and P. bicolor, sampled from Cupressus sempervirens trees in different regions in Israel. Arbitrarily primed PCR, performed on genomic DNA of 302 isolates, clustered the fungal population into five distinct groups. Multilocus phylogeny, split-network analyses and morphological characterization identified the isolates as Geosmithia omnicola, Geosmithia langdonii, Geosmithia sp. 708b, Geosmithia cupressina sp. nov. CBS147103 and Talaromyces cupressi sp. nov. CBS147104. Of these fungal isolates, G. cupressina and T. cupressi are newly described, and their morphological features and phylogenetic designations are presented. Inoculation of intact cypress saplings in an outdoor net-house revealed that only the representative isolate T. cupressi sp. nov. CBS147104 causes 100% disease incidence, whereas Geosmithia spp. isolates are not pathogenic. A number of these fungi were isolated from parasitoids that emerged from branch and stem sections colonized by P. armatus. This study suggests a long and stable association between Phloeosinus and Geosmithia species, and a possible role for additional associated fungal species as pathogens or endophytes of C. sempervirens trees in Israel.

Polymorphisms in Cha o 1 and Cha o 2, major allergens of Japanese cypress (Chamaecyparis obtusa) pollen from a restricted region in Japan.

Japanese cedar pollinosis is a major seasonal allergy in Japan, and Japanese cypress pollinosis is a growing concern because the cypress pollen season follows the cedar pollen season and cross-reactivity among allergens occurs between these closely related species. Allergens purified from pollen under unspecified collecting conditions can potentially heterogenous allergens profiles and batch to batch variability, and amino acid sequence variants in allergens possibly exist among trees. Polymorphisms have not been investigated for the cypress pollen major allergens, Cha o 1 and Cha o 2. Our aim was to examine the homogeneity of allergen amino acid sequences. DNA sequences of Cha o 1 and Cha o 2 from pollen collected from Chiba and Ibaraki prefectures and from needles of 47 plus trees located at seed orchards in Chiba Prefecture were examined by amplicon sequencing and amino acid substitutions were deduced. Sequence analysis of the pollen samples revealed that eight and seven residues of Cha o 2 were polymorphic, respectively. Thirteen residues in Cha o 2, including those residues identified in pollen, were deduced to be polymorphic for the plus trees. Cha o 2 expressed by the 47 plus trees included amino acid differences when compared with that of isoallergen Cha o 2.0101. No substitution was deduced in Cha o 1 for pollen taken from the two prefectures. One conservative amino acid substitution was deduced in Cha o 1 for the plus trees. Of the 47 plus trees examined, 38 were deduced to express only the isoallergen Cha o 1.0101 isoform, whereas eight trees were heterozygous and a single tree was homozygous for the non-synonymous mutation, which indicates relative uniformity of Cha o 1. Cha o 2 was found to be a heterogeneous allergen which suggests that studies using pollen from different trees may not give the same results.

Genome skimming and exploration of DNA barcodes for Taiwan endemic cypresses.

Cypresses are characterized by their longevity and valuable timber. In Taiwan, two endemic cypress species, Chamaecyparis formosensis and C. obtusa var. formosana, are threatened by prevalent illegal logging. A DNA barcode system is urgently needed for reforestation and conservation of these two cypresses. In this study, both plastomes and 35S rDNAs from 16, 10, and 6 individuals of C. formosensis, C. obtusa var. formosana, and C. obtusa var. obtusa were sequenced, respectively. We show that the loss of plastid trnT-GGU readily distinguishes C. formosensis from its congeneric species. We demonstrate that entire sequences of plastomes or 35S rDNAs are capable of correctly identifying cypress species and varieties, suggesting that they are effective super-barcodes. We also discover three short hypervariable loci (i.e., 3'ETS, ITS1, and trnH-psbA) that are promising barcodes for identifying cypress species and varieties. Moreover, nine species-specific indels of > 100 bp were detected in the cypress plastomes. These indels, together with the three aforementioned short barcodes, constitute an alternative and powerful barcode system crucial for identifying specimens that are fragmentary or contain degraded/poor DNA. Our sequenced data and barcode systems not only enrich the genetic reference for cypresses, but also contribute to future reforestation, conservation, and forensic investigations.

Allopatric divergence and hybridization within Cupressus chengiana (Cupressaceae), a threatened conifer in the northern Hengduan Mountains of western China.

Having a comprehensive understanding of population structure, genetic differentiation and demographic history is important for the conservation and management of threatened species. High-throughput sequencing (HTS) provides exciting opportunities to address a wide range of factors for conservation genetics. Here, we generated HTS data and identified 266,884 high-quality single nucleotide polymorphisms from 82 individuals of Cupressus chengiana, to assess population genomics across the species' full range, comprising the Daduhe River (DDH), Minjiang River (MJR) and Bailongjiang River (BLJ) catchments in western China. admixture, principal components analysis and phylogenetic analyses indicated that each region contains a distinct lineage, with high levels of differentiation between them (DDH, MJR and BLJ lineages). MJR was newly distinguished compared to previous surveys, and evidence including coalescent simulations supported a hybrid origin of MJR during the Quaternary. Each of these three lineages should be recognized as an evolutionarily significant unit (ESU), due to isolation, differing genetic adaptations and different demographic history. Currently, each ESU faces distinct threats, and will require different conservation strategies. Our work shows that population genomic approaches using HTS can reconstruct the complex evolutionary history of threatened species in mountainous regions, and hence inform conservation efforts, and contribute to the understanding of high biodiversity in mountains.

Ancient introgression drives adaptation to cooler and drier mountain habitats in a cypress species complex.

Introgression may act as an important source of new genetic variation to facilitate the adaptation of organisms to new environments, yet how introgression might enable tree species to adapt to higher latitudes and elevations remains unclear. Applying whole-transcriptome sequencing and population genetic analyses, we present an example of ancient introgression from a cypress species (Cupressus gigantea) that occurs at higher latitude and elevation on the Qinghai-Tibet Plateau into a related species (C. duclouxiana), which has likely aided the latter species to extend its range by colonizing cooler and drier mountain habitats during postglacial periods. We show that 16 introgressed candidate adaptive loci could have played pivotal roles in response to diverse stresses experienced in a high-elevation environment. Our findings provide new insights into the evolutionary history of Qinghai-Tibet Plateau plants and the importance of introgression in the adaptation of species to climate change.

In-depth transcriptome characterization uncovers distinct gene family expansions for Cupressus gigantea important to this long-lived species' adaptability to environmental cues.

BACKGROUND: Cupressus gigantea, a rare and endangered tree species with remarkable medicinal value, is endemic to the Tibetan Plateau. Yet, little is known about the underlying genetics of the unique ecological adaptability of this extremely long-lived conifer with a large genome size. Here, we present its first de novo and multi-tissue transcriptome in-depth characterization. RESULTS: We performed Illumina paired-end sequencing and RNA libraries assembly derived from terminal buds, male and female strobili, biennial leaves, and cambium tissues taken from adult C. gigantea. In total, large-scale high-quality reads were assembled into 101,092 unigenes, with an average sequence length of 1029 bp, and 6848 unigenes (6.77%) were mapped against the KEGG databases to identify 292 pathways. A core set of 41,373 genes belonging to 2412 orthologous gene families shared between C. gigantea and nine other plants was revealed. In addition, we identified 2515 small to larger-size gene families containing in total 9223 genes specific to C. gigantea, and enriched for gene ontologies relating to biotic interactions. We identified an important terpene synthases gene family expansion with its 121 putative members. CONCLUSIONS: This study presents the first comprehensive transcriptome characterization of C. gigantea. Our results will facilitate functional genomic studies to support genetic improvement and conservation programs for this endangered conifer.

Phylogenomic evidence for ancient recombination between plastid genomes of the Cupressus-Juniperus-Xanthocyparis complex (Cupressaceae).

BACKGROUND: Phylogenetic relationships among Eastern Hemisphere cypresses, Western Hemisphere cypresses, junipers, and their closest relatives are controversial, and generic delimitations have been in flux for the past decade. To address relationships and attempt to produce a more robust classification, we sequenced 11 new plastid genomes (plastomes) from the five variously described genera in this complex (Callitropsis, Cupressus, Hesperocyparis, Juniperus, and Xanthocyparis) and compared them with additional plastomes from diverse members of Cupressaceae. RESULTS: Phylogenetic analysis of protein-coding genes recovered a topology in which Juniperus is sister to Cupressus, whereas a tree based on whole plastomes indicated that the Callitropsis-Hesperocyparis-Xanthocyparis (CaHX) clade is sister to Cupressus. A sliding window analysis of site-specific phylogenetic support identified a ~ 15 kb region, spanning the genes ycf1 and ycf2, which harbored an anomalous signal relative to the rest of the genome. After excluding these genes, trees based on the remainder of the genes and genome consistently recovered a topology grouping the CaHX clade and Cupressus with strong bootstrap support. In contrast, trees based on the ycf1 and ycf2 region strongly supported a sister relationship between Cupressus and Juniperus. CONCLUSIONS: These results demonstrate that standard phylogenomic analyses can result in strongly supported but conflicting trees. We suggest that the conflicting plastomic signals result from an ancient introgression event involving ycf1 and ycf2 that occurred in an ancestor of this species complex. The introgression event was facilitated by plastomic recombination in an ancestral heteroplasmic individual carrying distinct plastid haplotypes, offering further evidence that recombination occurs between plastomes. Finally, we provide strong support for previous proposals to recognize five genera in this species complex: Callitropsis, Cupressus, Hesperocyparis, Juniperus, and Xanthocyparis.

The complete chloroplast genome of Cupressus gigantea, an endemic conifer species to Qinghai-Tibetan Plateau.

The complete chloroplast genome of the wild Cupressus gigantea (Cupressaceae) is determined in this study. The circular genome is 128 244 bp in length with 115 single copy genes and two duplicated genes (trnI-CAU and trnQ-UUG). This genome contains 82 protein-coding genes, four ribosomal RNA genes and 31 transfer RNA genes. In these genes, eight genes (atpF, rpoC1, ndhA, ndhB, petB, petD, rpl16 and rpl2) harbor a single intron and two genes (rps12 and ycf3) harbor two introns. This genome does not contain canonical IRs, and the overall GC content is 34.7%. A maximum parsimony phylogenetic analysis revealed that C. gigantea and C. sempervirens are more closely related.

Determination of epigenetic inheritance, genetic inheritance, and estimation of genome DNA methylation in a full-sib family of Cupressus sempervirens L.

Genetic inheritance and epigenetic inheritance are significant determinants of plant evolution, adaptation and plasticity. We studied inheritance of restriction site polymorphisms by the f-AFLP method and epigenetic DNA cytosine methylation inheritance by the f-MSAP technique. The study involved parents and 190 progeny of a Cupressus sempervirens L. full-sib family. Results from AFLP genetic data revealed that 71.8% of the fragments studied are under Mendelian genetic control, whereas faithful Mendelian inheritance for the MSAP fragments was low (4.29%). Further, MSAP fragment analysis showed that total methylation presented a mean of 28.2%, which was higher than the midparent value, while maternal inheritance was higher (5.65%) than paternal (3.01%). Interestingly de novo methylation in the progeny was high (19.65%) compared to parental methylation. Genetic and epigenetic distances for parents and offspring were not correlated (R(2)=0.0005). Furthermore, we studied correlation of total relative methylation and CG methylation with growth (height, diameter). We found CG/CNG methylation (N: A, C, T) to be positively correlated with height and diameter, while total relative methylation and CG methylation were positively correlated with height. Results are discussed in light of further research needed and of their potential application in breeding.

Identification of low temperature stress regulated transcript sequences and gene families in Italian cypress.

Cold acclimation is a complex transcriptionally controlled process regulated by many different genes and genic-interactions in plants. The northward spreading of woody species is mainly limited by winter harshness. To increase our knowledge about the biological processes underlying cold acclimation, plants evolved in warmer climates can serve as models. In this work, a Suppression Subtractive Hybridization approach using PCR-select was used to isolate Italian cypress (Cupressus sempervirens L.) transcript sequences putatively expressed under low temperature stress. After assessing the reliability of the subtractive step, a total of 388 clones were selected and sequenced. Following sequence assembly and removal of the redundant cDNAs, 156 unique transcripts were identified and annotated in order to assign them a putative functional class. Most of the identified transcripts were functionally classified pertaining to stress in cellular and chloroplast membranes, which are previously known to be severely damaged by cold treatment. Among the identified functional gene families, the extensively represented ones were dehydrins, early light-inducible proteins, senescence-associated genes and oleosins. The last three gene families were further selected for phylogenetic analysis, with the corresponding protein sequences across the complete genomes of the model plants Populus trichocarpa, Vitis vinifera, Physcomitrella patens, and Arabidopsis thaliana. The relationship with the ortholog sequences coming from these species and their further implications are discussed.

Production of the sesquiterpenoid (+)-nootkatone by metabolic engineering of Pichia pastoris.

The sesquiterpenoid (+)-nootkatone is a highly demanded and highly valued aroma compound naturally found in grapefruit, pummelo or Nootka cypress tree. Extraction of (+)-nootkatone from plant material or its production by chemical synthesis suffers from low yields and the use of environmentally harmful methods, respectively. Lately, major attention has been paid to biotechnological approaches, using cell extracts or whole-cell systems for the production of (+)-nootkatone. In our study, the yeast Pichia pastoris initially was applied as whole-cell biocatalyst for the production of (+)-nootkatone from (+)-valencene, the abundant aroma compound of oranges. Therefore, we generated a strain co-expressing the premnaspirodiene oxygenase of Hyoscyamus muticus (HPO) and the Arabidopsis thaliana cytochrome P450 reductase (CPR) that hydroxylated extracellularly added (+)-valencene. Intracellular production of (+)-valencene by co-expression of valencene synthase from Callitropsis nootkatensis resolved the phase-transfer issues of (+)-valencene. Bi-phasic cultivations of P. pastoris resulted in the production of trans-nootkatol, which was oxidized to (+)-nootkatone by an intrinsic P. pastoris activity. Additional overexpression of a P. pastoris alcohol dehydrogenase and truncated hydroxy-methylglutaryl-CoA reductase (tHmg1p) significantly enhanced the (+)-nootkatone yield to 208mg L(-1) cell culture in bioreactor cultivations. Thus, metabolically engineered yeast P. pastoris represents a valuable, whole-cell system for high-level production of (+)-nootkatone from simple carbon sources.

Genetic diversity and conservation implications of four Cupressus species in China as revealed by microsatellite markers.

Understanding the extent and distribution of genetic diversity is crucial for the conservation and management of endangered species. Cupressus chengiana, C. duclouxiana, C. gigantea, and C. funebris are four ecologically and economically important species in China. We investigated their genetic diversity, population structure, and extant effective population size (35 populations, 484 individuals) employing six pairs of nuclear microsatellite markers (selected from 53). Their genetic diversity is moderate among conifers, and genetic differentiation among populations is much lower in C. gigantea than in the other three species; the estimated effective population size was largest for C. chengiana, at 1.70, 2.91, and 3.91 times the estimates for C. duclouxiana, C. funebris, and C. gigantea, respectively. According to Bayesian clustering analysis, the most plausible population subdivision scheme within species is two groups in C. chengiana, three groups in C. duclouxiana, and a single group for both C. funebris and C. gigantea. We propose a conservation strategy for these cypress species.

Linkage mapping of the Mediterranean cypress, Cupressus sempervirens, based on molecular and morphological markers.

Gene mapping for a Cupressus species is presented for the first time. Two linkage maps for the Mediterranean cypress (Cupressus sempervirens) varieties, C. sempervirens var. horizontalis and C. sempervirens var. pyramidalis, were constructed following the pseudo-testcross mapping strategy and employing RAPD, SCAR and morphological markers. A total of 427 loci (425 RAPDs, two SCARs) representing parents and F(1) progeny were screened for polymorphism with 32 random decamer and two SCAR primers. A morphological marker defined as "crown form" was also included. Of 274 polymorphic loci, the 188 that presented Mendelian inheritance formed the mapping dataset. Of these loci, 30% were mapped into seven linkage groups for the horizontalis (maternal) and four linkage groups for the pyramidalis (paternal) map. The putative "crown form" locus was included in a linkage group of both maps. The horizontalis and the pyramidalis maps covered 160.1 and 144.5 cM, respectively, while genome length was estimated to be 1696 cM for the former variety and 1373 cM for the latter. The four RAPD markers most tightly linked to crown form were cloned and converted to SCARs. Each of the cloned RAPD markers yielded two to three different sequences behaving as co-migrating fragments. Two SCAR markers, SC-D05(432) and SC-D09(667), produced amplified bands of the expected sizes and maintained linkage with the appropriate phenotype, but to a lesser extent compared to their original RAPD counterparts. These linkage maps represent a first step towards the localization of QTLs and genes controlling crown form and other polygenic traits in cypress.

Plastid DNA diversity is higher in the island endemic Guadalupe cypress than in the continental Tecate cypress.

BACKGROUND: Callitropsis guadalupensis (Guadalupe cypress) is endemic to Guadalupe Island, Mexico, where it is the dominant species of the only forest. The species has suffered declining numbers following the introduction of goats to the island over 150 years ago. Callitropsis guadalupensis is closely related to Callitropsis forbesii (Tecate cypress), distributed in small isolated populations in mainland Baja California and southern California. The objective of the present study was to compare the genetic diversity of the island endemic to the continental species. METHODOLOGY/PRINCIPAL FINDINGS: We measured genetic diversity in Callitropsis guadalupensis (n =54) from Guadalupe Island and in Callitropsis forbesii (n = 100) from five populations in mainland Baja California. The plastid DNA trnS-trnG spacer and the trnL-trnF region were chosen for characterization. Thirty-four haplotypes were observed, of which six were shared between both species. One of these haplotypes was also shared with three other species, Callitropsis lusitanica, Callitropsis montana, and Callitropsis stephensonii. Haplotype diversity (h) and nucleotide diversity (π) were significantly higher for Callitropsis guadalupensis (h = 0.698, π = 0.00071) than for Callitropsis forbesii (h = 0.337, π = 0.00024). CONCLUSIONS/SIGNIFICANCE: Callitropsis guadalupensis shows no evidence of a founder effect or of a genetic bottleneck, and can be added to a growing list of insular species with higher genetic diversity than their mainland relatives.

Molecular cloning and characterization of Cup a 4, a new allergen from Cupressus arizonica.

Sensitization to Cupressaceae pollen has become one of the most important causes of pollinosis in Western countries during winter and early spring. However, the characterization of the extracts, the allergens involved and the cross-reactivity with other pollen sources still remain poorly studied; in the case of Cupressus arizonica only two allergens have been described so far. A new allergen from C. arizonica pollen, Cup a 4, was cloned and expressed in Escherichia coli as an N-terminally His-tag recombinant protein that was characterized biochemically, immunologically and by circular dichroism spectroscopy. The new allergen has high sequence identity with Prickly Juniper allergen Jun o 4 and contains four EF-hand domains. The recombinant protein has structural similarities with other calcium binding allergens such as Ole e 3, Ole e 8 and Phl p 7. Cup a 4 is expressed in mature pollen grains and shares antigenic properties with the recombinant form. Sera from 9.6% C. arizonica allergic patients contain specific IgE antibodies against recombinant Cup a 4.

Phylogeography and allopatric divergence of cypress species (Cupressus L.) in the Qinghai-Tibetan Plateau and adjacent regions.

BACKGROUND: Although allopatric speciation is viewed as the most common way in which species originate, allopatric divergence among a group of closely related species has rarely been examined at the population level through phylogeographic analysis. Here we report such a case study on eight putative cypress (Cupressus) species, which each have a mainly allopatric distribution in the Qinghai-Tibetan Plateau (QTP) and adjacent regions. The analysis involved sequencing three plastid DNA fragments (trnD-trnT, trnS-trnG and trnL-trnF) in 371 individuals sampled from populations at 66 localities. RESULTS: Both phylogenetic and network analyses showed that most DNA haplotypes recovered or haplotype-clustered lineages resolved were largely species-specific. Across all species, significant phylogeographic structure (N(ST) > G(ST), P < 0.05) implied a high correlation between haplotypes/lineages and geographic distribution. Two species, C. duclouxiana and C. chengiana, which are distributed in the eastern QTP region, contained more haplotypes and higher diversity than five species with restricted distributions in the western highlands of the QTP. The remaining species, C. funebris, is widely cultivated and contained very little cpDNA diversity. CONCLUSIONS: It is concluded that the formation of high mountain barriers separating deep valleys in the QTP and adjacent regions caused by various uplifts of the plateau since the early Miocene most likely promoted allopatric divergence in Cupressus by restricting gene flow and fixing local, species-specific haplotypes in geographically isolated populations. The low levels of intraspecific diversity present in most species might stem from population bottlenecks brought about by recurrent periods of unfavorable climate and more recently by the negative impacts of human activities on species' distributions. Our findings shed new light on the importance of geographical isolation caused by the uplift of the QTP on the development of high plant species diversity in the QTP biodiversity hotspot.

High genetic variation in marginal fragmented populations at extreme climatic conditions of the Patagonian Cypress Austrocedrus chilensis.

Knowledge about current patterns of genetic structure of populations together with the evolutionary history of a species helps to understand and predict the adaptation of populations to future climate change. We assayed variation at nuclear microsatellite markers among peripheral vs. continuous populations of the temperate South American species Austrocedrus chilensis, to investigate the role of historical vs. demographical forces in shaping population genetic structure. This species occurs in continuous populations in the west and central distribution range, but becomes highly fragmented at the eastern limit, which comprised ice-free areas during Quaternary glaciations and has extreme climatic conditions at present times. Bayesian analysis methods identified two contrasting patterns of genetic structure; (I) populations from humid, mesic and peri-glacial regions formed a single deme with relatively low genetic differentiation and high admixture levels whereas (II) a highly heterogeneous genetic structure with low level of admixture was found in the steppe, towards the east and northeast limit of the distribution range. In the steppe, population fragmentation, restricted gene flow and isolation-by-distance were also inferred. In addition, several small steppe populations showed high genetic diversity and divergent gene pools, suggesting that they constitute ancient refuges from pre-Holocene glaciations with just a subgroup of them contributing significantly to post-glacial spread. These results are discussed in relation to patterns of genetic variation found for other temperate species and the contribution of the particular southern Andes topography and climate to post-glacial spread.

Molecular evidence for the natural production of homozygous Cupressus sempervirens L. lines by Cupressus dupreziana seed trees.

Paternal apomixis was recently reported in the endangered Mediterranean cypress, Cupressus dupreziana. This species acts as a surrogate mother for the development of all-paternal embryos from pollen grains. C. dupreziana production of Cupressus sempervirens haploid or diploid seedlings from C. sempervirens pollen was also demonstrated. The haploid progeny was derived from the embryogenic development of haploid gametes, but the origin of the diploid progeny remained unknown. To determine the ontogenic origin of the diploid C. sempervirens progeny, we analyzed the heterozygozity of 63 diploid all-paternal C. sempervirens seedlings using highly variable co-dominant nuclear microsatellite markers. The bi-parental inheritance of the markers was checked in C. sempervirens controlled crosses. A high level of polymorphism was observed among the diploid all-paternal trees. All but three individuals exhibited single-band profiles as expected for homozygotes, which may originate from natural diploidization of a C. sempervirens haploid embryo or from the fusion of two male gametes produced by the same C. sempervirens microgametophyte. The three heterozygous seedlings must be derived from the fusion of male gametes produced by two different C. sempervirens microgametophytes. These findings offer a unique opportunity in conifers to produce homozygous lines, highly valuable for genetic analyses or breeding.