Comparative Transcriptome Analysis of Male and Female Conelets and Development of Microsatellite Markers in Pinus bungeana, an Endemic Conifer in China

Sex-specific responses of Taxus mairei to UV-B radiation involved altering the interactions between the microbiota assembly and host secondary metabolism

BACKGROUND

To adapt to constantly changing environments, ancient gymnosperms have coevolved with diverse endophytic fungi that are essential for the fitness and adaptability of the plant host. However, the effect of sex on plant-endophyte interactions in response to environmental stressors remains unknown. RNA-seq integrated with ITS analysis was applied to reveal the potential mechanisms underlying the sex-specific responses of Taxus mairei to ultraviolet (UV)-B radiation.

RESULTS

Enrichment analysis suggested that sex influenced the expression of several genes related to the oxidation-reduction system, which might play potential roles in sex-mediated responses to UV-B radiations. ITS-seq analysis clarified the effects of UV-B radiation and sex on the composition of endophytic fungal communities. Sex influenced various secondary metabolic pathways, thereby providing chemicals for T. mairei host to produce attractants and/or inhibitors to filter microbial taxa. Analysis of fungal biomarkers suggested that UV-B radiation reduced the effect of sex on fungal communities. Moreover, Guignardia isolate #1 was purified to investigate the role of endophytic fungi in sex-mediated responses to UV-B radiation. Inoculation with spores produced by isolate #1 significantly altered various oxidation-reduction systems of the host by regulating the expression of APX2, GST7 NCED1, ZE1, CS1, and CM1.

CONCLUSION

These results revealed the roles of endophytic fungi in sex-mediated responses to UV-B radiation and provided novel insights into the sex-specific responses of Taxus trees to environmental stressors. Video Abstract.

Genetic diversity, population structure, and DNA fingerprinting of Ailanthus altissima var. erythrocarpa based on EST-SSR markers

Ailanthus altissima var. erythrocarpa is an A. altissima variety with high economic, ecological and ornamental value, but there have been no reports on the development of SSR primers for it. According to the SSR primer information provided by the transcriptome of A. altissima var. erythrocarpa, 120 individuals with different redness levels were used to screen polymorphic primers. Transcriptomic analysis revealed 10,681 SSR loci, of which mononucleotide repeats were dominant (58.3%), followed by dinucleotide and trinucleotide repeats (16.6%, 15.1%) and pentanucleotide repeats (0.2%). Among 140 pairs of randomly selected primers, nineteen pairs of core primers with high polymorphism were obtained. The average number of alleles (Na), average number of effective alleles (Ne), average Shannon's diversity index (I), average observed heterozygosity (Ho), average expected heterozygosity (He), fixation index (F) and polymorphic information content (PIC) were 11.623, 4.098, 1.626, 0.516, 0.696, 0.232 and 0.671, respectively. Nineteen EST-SSR markers were used to study the genetic diversity and population structure of A. altissima var. erythrocarpa. The phylogenetic tree, PCoA, and structure analysis all divided the tested resources into two categories, clearly showing the genetic variation between individuals. The population showed high genetic diversity, mainly derived from intraspecific variation. Among nineteen pairs of primers, 4 pairs (p33, p15, p46, p92) could effectively distinguish and be used for fingerprinting of the tested materials. This study is of great significance for genetic diversity analysis and molecular-assisted breeding of A. altissima var. erythrocarpa.

Geographic isolation and environmental heterogeneity contribute to genetic differentiation in Cephalotaxus oliveri

Evaluating the contributions of geographic distance and environmental heterogeneity to the genetic divergence can inform the demographic history and responses to environmental change of natural populations. The isolation-by-distance (IBD) reveals that genetic differentiation among populations increases with geographic distance, while the isolation-by-environment (IBE) assumes a linear relationship between genetic variation and environmental differences among populations. Here, we sampled and genotyped 330 individuals from 18 natural populations of Cephalotaxus oliveri throughout the species' distribution. Twenty-eight EST-SSR markers were applied to analyze population genetics, for the investigation of the driving factors that shaped spatial structure. In addition, we identified the outlier loci under positive selection and tested their association with environmental factors. The results showed a moderate genetic diversity in C. oliveri and high genetic differentiation among populations. Population structure analyses indicated that 18 populations were clustered into two major groups. We observed that the genetic diversity of central populations decreased and the genetic differentiation increased towards the marginal populations. Additionally, the signatures of IBD and IBE were detected in C. oliveri, and IBE provided a better contribution to genetic differentiation. Six outlier loci under positive selection were demonstrated to be closely correlated with environmental variables, among which bio8 was associated with the greatest number of loci. Genetic evidence suggests the consistency of the central-marginal hypothesis (CMH) for C. oliveri. Furthermore, our results suggest that temperature-related variables played an important role in shaping genetic differentiation.

Alternaria alternata, the Causal Agent of a New Needle Blight Disease on Pinus bungeana

Pinus bungeana, an endangered and native coniferous tree species in China, has considerable timber and horticulture value. However, little is known about needle diseases in P. bungeana. A needle blight of P. bungeana has been observed in Hebei Province, China. P. bungeana inoculated with mycelial plugs of fungal isolates presented symptoms similar to those observed under field conditions. Ten virulent fungal isolates were identified as a small-spored Alternaria species based on morphological observations. Maximum likelihood and Bayesian phylogenetic analyses carried out with multilocus sequence typing of eight regions (SSU, LSU, ITS, gapdh, tef1, Alt a 1, endoPG, OPA10-2) assigned the pathogen to Alternaria alternata. This is the first report of A. alternata causing needle blight on P. bungeana in China.

Full-length transcriptome analysis of Zanthoxylum nitidum (Roxb.) DC

Zanthoxylum nitidum (Roxb.) DC. (Z. nitidum) is a type of Chinese Dao-di herb, also called Liangmianzhen, which is widely used to treat arthralgia, rheumatic arthralgia, and stomach pain. However, genomic resources for Z. nitidum are still scarce. This study provides transcriptomic resources for Z. nitidum by applying single-molecule real-time (SMRT) sequencing technology. In total, 456,109 circular consensus sequencing (CCS) reads were generated with a mean length of 2,216 bp from Z. nitidum roots, old stems, young branches, leaves, flowers, and fruits. Of these total reads, 353,932 were full-length nonchimeric (FLNC) reads with an average length of 1,996 bp. A total of 16,163 transcripts with a mean length of 1,171 bp were acquired. Of these transcripts, 14,231 (88%) were successfully annotated using public databases. Across all the 16,163 transcripts, we identified 6,255 long non-coding RNAs (lncRNAs) and 22,780 simple sequence repeats (SSRs). Furthermore, 3,482 transcription factors were identified. Among the SSR loci, 1-3 nucleotide repeats were dominant, occupying 99.36% of the total SSR loci, with mono-, di-, and tri-nucleotide repeats accounting for 61.80%, 19.89%, and 5.02% of the total SSR loci, respectively. A total of 36 out of 100 randomly selected primer pairs were verified to be positive, 20 of which showed polymorphism. These findings enrich the genetic resources available for facilitating future studies and research on relevant topics such as population genetics in Z. nitidum.

Regulation of lignin biosynthesis by an atypical bHLH protein CmHLB in Chrysanthemum

Stem mechanical strength is one of the most important agronomic traits that affects the resistance of plants against insects and lodging, and plays an essential role in the quality and yield of plants. Several transcription factors regulate mechanical strength in crops. However, mechanisms of stem strength formation and regulation remain largely unexplored, especially in ornamental plants. In this study, we identified an atypical bHLH transcription factor CmHLB (HLH PROTEIN INVOLVED IN LIGNIN BIOSYNTHESIS) in chrysanthemum, belonging to a small bHLH sub-family - the PACLOBUTRAZOL RESISTANCE (PRE) family. Overexpression of CmHLB in chrysanthemum significantly increased mechanical strength of the stem, cell wall thickness, and lignin content, compared with the wild type. In contrast, CmHLB RNA interference lines exhibited the opposite phenotypes. RNA-seq analysis indicated that CmHLB promoted the expression of genes involved in lignin biosynthesis. Furthermore, we demonstrated that CmHLB interacted with Chrysanthemum KNOTTED ARABIDOPSIS THALIANA7 (CmKNAT7) through the KNOX2 domain, which has a conserved function, i.e. it negatively regulates secondary cell wall formation of fibres and lignin biosynthesis. Collectively, our results reveal a novel role for CmHLB in regulating lignin biosynthesis by interacting with CmKNAT7 and affecting stem mechanical strength in Chrysanthemum.

Transcriptional Regulation of Pine Male and Female Cone Initiation and Development: Key Players Identified Through Comparative Transcriptomics

With long reproductive timescales, large complex genomes, and a lack of reliable reference genomes, understanding gene function in conifers is extremely challenging. Consequently, our understanding of which genetic factors influence the development of reproductive structures (cones) in monoecious conifers remains limited. Genes with inferred roles in conifer reproduction have mostly been identified through homology and phylogenetic reconstruction with their angiosperm counterparts. We used RNA-sequencing to generate transcriptomes of the early morphological stages of cone development in the conifer species Pinus densiflora and used these to gain a deeper insight into the transcriptional changes during male and female cone development. Paired-end Illumina sequencing was used to generate transcriptomes from non-reproductive tissue and male and female cones at four time points with a total of 382.82 Gbp of data generated. After assembly and stringent filtering, a total of 37,164 transcripts were retrieved, of which a third were functionally annotated using the Mercator plant pipeline. Differentially expressed gene (DEG) analysis resulted in the identification of 172,092 DEGs in the nine tissue types. This, alongside GO gene enrichment analyses, pinpointed transcripts putatively involved in conifer reproductive structure development, including co-orthologs of several angiosperm flowering genes and several that have not been previously reported in conifers. This study provides a comprehensive transcriptome resource for male and early female cone development in the gymnosperm species Pinus densiflora. Characterisation of this resource has allowed the identification of potential key players and thus provides valuable insights into the molecular regulation of reproductive structure development in monoecious conifers.

Development and Application of EST-SSR Markers in Cephalotaxus oliveri From Transcriptome Sequences

Cephalotaxus oliveri is an endemic conifer of China, which has medicinal and ornamental value. However, the limited molecular markers and genetic information are insufficient for further genetic studies of this species. In this study, we characterized and developed the EST-SSRs from transcriptome sequences for the first time. The results showed that a total of 5089 SSRs were identified from 36446 unigenes with a density of one SSR per 11.1 kb. The most common type was trinucleotide repeats, excluding mononucleotide repeats, followed by dinucleotide repeats. AAG/CTT and AT/AT exhibited the highest frequency in the trinucleotide and dinucleotide repeats, respectively. Of the identified SSRs, 671, 1125, and 1958 SSRs were located in CDS, 3′UTR, and 5′UTR, respectively. Functional annotation showed that the SSR-containing unigenes were involved in growth and development with various biological functions. Among successfully designed primer pairs, 238 primer pairs were randomly selected for amplification and validation of EST-SSR markers and 47 primer pairs were identified as polymorphic. Finally, 28 high-polymorphic primers were used for genetic analysis and revealed a moderate level of genetic diversity. Seven natural C. oliveri sampling sites were divided into two genetic groups. Furthermore, the 28 EST-SSRs had 96.43, 71.43, and 78.57% of transferability rate in Cephalotaxus fortune, Ametotaxus argotaenia, and Pseudotaxus chienii, respectively. These markers developed in this study lay the foundation for further genetic and adaptive evolution studies in C. oliveri and related species.

EST-SSR-based landscape genetics of Pseudotaxus chienii, a tertiary relict conifer endemic to China

Pseudotaxus chienii, belonging to the monotypic genus Pseudotaxus (Taxaceae), is a relict conifer endemic to China. Its populations are usually small and patchily distributed, having a low capacity of natural regeneration. To gain a clearer understanding of how landscape variables affect the local adaptation of P. chienii, we applied EST-SSR markers in conjunction with landscape genetics methods: (a) to examine the population genetic pattern and spatial genetic structure; (b) to perform genome scan and selection scan to identify outlier loci and the associated landscape variables; and (c) to model the ecological niche under climate change. As a result, P. chienii was found to have a moderate level of genetic variation and a high level of genetic differentiation. Its populations displayed a significant positive relationship between the genetic and geographical distance (i.e., "isolation by distance" pattern) and a strong fine-scale spatial genetic structure within 2 km. A putatively adaptive locus EMS6 (functionally annotated to cellulose synthase A catalytic subunit 7) was identified, which was found significantly associated with soil Cu, K, and Pb content and the combined effects of temperature and precipitation. Moreover, P. chienii was predicted to experience significant range contractions in future climate change scenarios. Our results highlight the potential of specific soil metal content and climate variables as the driving force of adaptive genetic differentiation in P. chienii. The data would also be useful to develop a conservation action plan for P. chienii.

Genetic Diversity and Population Differentiation of Pinus koraiensis in China

Pinus koraiensis is a well-known precious tree species in East Asia with high economic, ornamental and ecological value. More than fifty percent of the P. koraiensis forests in the world are distributed in northeast China, a region with abundant germplasm resources. However, these natural P. koraiensis sources are in danger of genetic erosion caused by continuous climate changes, natural disturbances such as wildfire and frequent human activity. Little work has been conducted on the population genetic structure and genetic differentiation of P. koraiensis in China because of the lack of genetic information. In this study, 480 P. koraiensis individuals from 16 natural populations were sampled and genotyped. Fifteen polymorphic expressed sequence tag-simple sequence repeat (EST-SSR) markers were used to evaluate genetic diversity, population structure and differentiation in P. koraiensis. Analysis of molecular variance (AMOVA) of the EST-SSR marker data showed that 33% of the total genetic variation was among populations and 67% was within populations. A high level of genetic diversity was found across the P. koraiensis populations, and the highest levels of genetic diversity were found in HH, ZH, LS and TL populations. Moreover, pairwise Fst values revealed significant genetic differentiation among populations (mean Fst = 0.177). According to the results of the STRUCTURE and Neighbor-joining (NJ) tree analyses and principal component analysis (PCA), the studied geographical populations cluster into two genetic clusters: cluster 1 from Xiaoxinganling Mountains and cluster 2 from Changbaishan Mountains. These results are consistent with the geographical distributions of the populations. The results provide new genetic information for future genome-wide association studies (GWAS), marker-assisted selection (MAS) and genomic selection (GS) in natural P. koraiensis breeding programs and can aid the development of conservation and management strategies for this valuable conifer species.

Development and characterization of Novel EST-based single-copy genic microsatellite DNA markers in white spruce and black spruce

Due mainly to large genome size and prevalence of repetitive sequences in the nuclear genome of spruce (Picea Mill.), it is very difficult to develop single-copy genomic microsatellite markers. We have developed and characterized 25 polymorphic, single-copy genic microsatellites from white spruce (Picea glauca (Moench) Voss) EST sequences and determined their informativeness in white spruce and black spruce (Picea mariana (Mill.) B.S.P.) and inheritance in black spruce. White spruce EST sequences from NCBI dbEST were searched for the presence of microsatellite repeats. Forty-seven sequences containing dinucleotide, trinucleotide, tetranucleotide and compound repeats were selected to develop primers. Twenty-five of the designed primer pairs yielded scorable amplicons, with single-locus patterns, and were characterized in 20 individuals each of white spruce and black spruce. All 25 microsatellites were polymorphic in white spruce and 24 in black spruce. The number of alleles at a locus ranged from two to 18, with a mean of 8.8 in white spruce, and from one to 17, with a mean of 7.6 in black spruce. The expected heterozygosity/polymorphic information content ranged from 0.10 to 0.92, with a mean of 0.67 in white spruce, and from 0 to 0.93, with a mean of 0.59 in black spruce. Microsatellites with dinucleotide and compound repeats were more informative than those with trinucleotide and tetranucleotide repeats. Eighteen microsatellite markers polymorphic between the parents of a black spruce controlled cross inherited in a single-locus Mendelian fashion. The microsatellite markers developed can be applied for various genetics, genomics, breeding, and conservation studies and applications.

New insight into the rapid growth of the Mikania micrantha stem based on DIA proteomic and RNA-Seq analysis

Mikania micrantha is one of the world's most invasive plants, which causes severe damage to natural ecosystems and agroforestry systems due to its rapid stem growth. This work investigated the proteomic and transcriptomic profiles of M. micrantha in different stem tissues (pre-internode, post-internode, and internode), as well as in adventitious roots and primary roots with the final goal of elucidating differentially expressed genes and proteins responsible for the rapid growth of stem. The objective was approached by using DIA-based proteomic and RNA-Seq technologies. More than seven giga-transcriptome clean reads were sequenced, and 5196 protein species were identified. Differentially expressed genes identified in all stem tissues were significantly enriched in photosynthesis and carbon fixation, suggesting that the stem possesses a strong photosynthetic capacity in order to maintain the energy supply for this species. Analysis of differentially expressed proteins showed that proteins related to photosystem I/II and the cytochrome b6/f complex, such as D1, D2, and cp43, were also highly accumulated in the adventitious roots, corroborating the transcriptome analysis results. These results provided basic proteomic and transcriptional expression information about the M. micrantha stem and adventitious root, thereby improving our understanding of the molecular mechanism underlying rapid growth in this species. SIGNIFICANCE: This is the first study to investigate the proteomic and transcriptomic profiles of Mikania micrantha, a highly invasive plant, in different stem tissues (pre-internode, post-internode, and internode), as well as in adventitious and primary roots, using the latest DIA-based (data-independent acquisition mode) proteomic and RNA-Seq technologies. A comprehensive study was carried out, and differentially expressed genes and differentially expressed proteins identified in the pre-internode, post-internode, and internode tissues were significantly enriched during photosynthesis and carbon fixation, suggesting that the M. micrantha stem possesses a strong photosynthetic capacity that allows the plant to maintain a high energy supply. Enriched plant hormone signal transduction pathway analysis revealed an interaction between auxin and other phytohormones involved in adventitious root development. The study provided basic data on the molecular mechanism of M. micrantha vegetative propagation and the rapid growth of its stem. The novel scientific content of this study successfully builds upon the limited information currently available on the subject, therefore warranting publication.

Transcriptome Analysis in Male Strobilus Induction by Gibberellin Treatment in Cryptomeria japonica D. Don

The plant hormone gibberellin (GA) is known to regulate elongating growth, seed germination, and the initiation of flower bud formation, and it has been postulated that GAs originally had functions in reproductive processes. Studies on the mechanism of induction of flowering by GA have been performed in Arabidopsis and other model plants. In coniferous trees, reproductive organ induction by GAs is known to occur, but there are few reports on the molecular mechanism in this system. To clarify the gene expression dynamics of the GA induction of the male strobilus in Cryptomeria japonica, we performed comprehensive gene expression analysis using a microarray. A GA-treated group and a nontreated group were allowed to set, and individual trees were sampled over a 6-week time course. A total of 881 genes exhibiting changed expression was identified. In the GA-treated group, genes related to ‘stress response’ and to ‘cell wall’ were initially enriched, and genes related to ‘transcription’ and ‘transcription factor activity’ were enriched at later stages. This analysis also clarified the dynamics of the expression of genes related to GA signaling transduction following GA treatment, permitting us to compare and contrast with the expression dynamics of genes implicated in signal transduction responses to other plant hormones. These results suggested that various plant hormones have complex influences on the male strobilus induction. Additionally, principal component analysis (PCA) using expression patterns of the genes that exhibited sequence similarity with flower bud or floral organ formation-related genes of Arabidopsis was performed. PCA suggested that gene expression leading to male strobilus formation in C. japonica became conspicuous within one week of GA treatment. Together, these findings help to clarify the evolution of the mechanism of induction of reproductive organs by GA.

Development and Transferability of EST-SSR Markers for Pinus koraiensis from Cold-Stressed Transcriptome through Illumina Sequencing

Pinus koraiensis has significant economic and ecological value in Northeast China. However, due to the lack of suitable molecular markers, only a few available microsatellite markers were developed for further population genetics studies. In this study, for the first time we developed expressed sequence tag-simple sequence repeat (EST-SSR) markers from the cold-stressed transcriptome of P. koraiensis using Illumina Sequencing. We identified a total of 7,235 EST-SSRs from 97,376 sequences, and we tested their transferability among seven related Pinus species. The results showed that trinucleotides were the most abundant type of repeat (1287, 18.74%) excluding mononucleotides, followed by dinucleotides (1284, 18.7%) and tetranucleotides (72, 1.05%). The most dominant dinucleotides and trinucleotide repeat motifs were AT/AT (535, 7.79%) and AAT/ATT (103, 1.5%). The observed heterozygosity (Ho) and expected heterozygosity (He) ranged from 0.002 to 0.986 and 0.017 to 0.743, respectively, and the polymorphism information content (PIC) values and number of alleles (Na) varied from 0.029 to 0.794 and 2 to 23, respectively. A total of 8 natural P. koraiensis populations were divided into two main genetic clusters. Furthermore, nine of twenty polymorphic primer pairs were successfully amplified in seven Pinus species, and at least 80% of the successful P. koraiensis EST-SSR primers could be amplified in more than four species (16, 80%). Combined results for the development of EST-SSR markers in P. koraiensis and transferability among related species would contribute to improved studies on the genetic diversity and population structure in P. koraiensis and phylogenetic relationships among Pinus species. They would also provide a significant source for quantitative trait locus analysis.

De novo assembly of the Platycladus orientalis (L.) Franco transcriptome provides insight into the development and pollination mechanism of female cone based on RNA-Seq data

For seed-bearing plants, the basis of seed and fruit formation is pollination. The normal progression of pollination is through advances in continuous signal exchange and material transfer, which occur mainly in female reproductive organs; thus, the molecular mechanism of development in female reproductive organs is vital for understanding the principle of pollination. However, molecular biology studies on the development of female cones related to pollination are rare and unclear in gymnosperms, especially in Cupressaceae. In this study, Platycladus orientalis, a monotypic genus within Cupressaceae, was chosen to examine female cone transcriptomes at pre-pollination and pollination stages by Illumina paired-end sequencing technology to de novo sequence six libraries with 3 biological replicates. These libraries were used to construct a P. orientalis transcriptome database containing 71,669 unigenes (4,963 upregulated unigenes and 11,747 downregulated unigenes at the pollination stage) for subsequent analysis. Based on the annotations and expression levels, the functions of differentially expressed unigenes and enriched pathways between the developmental processes of female cones were analysed to detail the preliminary development and pollination mechanism of the female cone. Targeted investigations were specifically performed to determine the elementary mechanism of secretion and functioning of the pollination drop, a vital ovule secretion at the pollination stage. Ultimately, the expression of 15 unigenes selected between two stages were further assessed and confirmed using qRT-PCR, which demonstrated reliable data and significant differences in the expression profiles of key genes. As one of the largest available transcriptomic resources of this species, the database is constructed to prospectively adapt to the physiological and genomic data of woody plants. This work provided the first transcriptome profile of P. orientalis female cones at different developmental stages, and will promote the illumination of the pollination mechanism of P. orientalis, and will serve as the basis for in-depth genomic study in the Cupressaceae family. This initiative will arouse the interest and attention of scholars and pave the way for future studies.

An Integrated Transcriptome and Proteome Analysis Reveals Putative Regulators of Adventitious Root Formation in Taxodium 'Zhongshanshan'

Adventitious root (AR) formation from cuttings is the primary manner for the commercial vegetative propagation of trees. Cuttings is also the main method for the vegetative reproduction of Taxodium ‘Zhongshanshan’, while knowledge of the molecular mechanisms regulating the processes is limited. Here, we used mRNA sequencing and an isobaric tag for relative and absolute quantitation-based quantitative proteomic (iTRAQ) analysis to measure changes in gene and protein expression levels during AR formation in Taxodium ‘Zhongshanshan’. Three comparison groups were established to represent the three developmental stages in the AR formation process. At the transcript level, 4743 genes showed an expression difference in the comparison groups as detected by RNA sequencing. At the protein level, 4005 proteins differed in their relative abundance levels, as indicated by the quantitative proteomic analysis. A comparison of the transcriptome and proteome data revealed regulatory aspects of metabolism during AR formation and development. In summary, hormonal signal transduction is different at different developmental stages during AR formation. Other factors related to carbohydrate and energy metabolism and protein degradation and some transcription factor activity levels, were also correlated with AR formation. Studying the identified genes and proteins will provide further insights into the molecular mechanisms controlling AR formation.

Pinus massoniana Introgression Hybrids Display Differential Expression of Reproductive Genes

Pinus massoniana and P. hwangshanensis are two conifer species located in southern China, which are of both economic and ornamental value. Around the middle and lower reaches of the Yangtze River, P. massoniana occurs mainly at altitudes below 700 m, while P. hwangshanensis can be found above 900 m. At altitudes where the distribution of both pines overlaps, a natural introgression hybrid exists, which we will further refer to as the Z pine. This pine has a morphological character that shares attributes of both P. massoniana and P. hwangshanensis. However, compared to the other two pines, its reproductive structure, the pinecone, has an ultra-low ripening rate with seeds that germinate poorly. In this study, we aimed to find the reason for the impaired cone maturation by comparing transcriptome libraries of P. massoniana and Z pine cones at seven successive growth stages. After sequencing and assembly, we obtained unigenes and then annotated them against NCBI’s non-redundant nucleotide and protein sequences, Swiss-Prot, Clusters of Orthologous Groups, Gene Ontology and KEGG Orthology databases. Gene expression levels were estimated and differentially expressed genes (DEGs) of the two pines were mined and analyzed. We found that several of them indeed relate to reproductive process. At every growth stage, these genes are expressed at a higher level in P. massoniana than in the Z pine. These data provide insight into understanding which molecular mechanisms are altered between P. massoniana and the Z pine that might cause changes in the reproductive process.