Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Sex-Biased Temporal Gene Expression in Male and Female Floral Buds of Seabuckthorn (Hippophae rhamnoides)

  • Aseem Chawla,

    Affiliation Department of Biotechnology and Bioinformatics, Jaypee University of Information Technology, Waknaghat, Solan, India

  • Tsering Stobdan,

    Affiliation Defence Institute of High Altitude Research, Defence R & D Organisation, Leh, Jammu, and Kashmir, India

  • Ravi B. Srivastava,

    Affiliation Defence Institute of High Altitude Research, Defence R & D Organisation, Leh, Jammu, and Kashmir, India

  • Varun Jaiswal,

    Affiliation Department of Biotechnology and Bioinformatics, Jaypee University of Information Technology, Waknaghat, Solan, India

  • Rajinder S. Chauhan,

    Affiliation Department of Biotechnology and Bioinformatics, Jaypee University of Information Technology, Waknaghat, Solan, India

  • Anil Kant

    anil_kantv@yahoo.com

    Affiliation Department of Biotechnology and Bioinformatics, Jaypee University of Information Technology, Waknaghat, Solan, India

Abstract

Seabuckthorn is an economically important dioecious plant in which mechanism of sex determination is unknown. The study was conducted to identify seabuckthorn homologous genes involved in floral development which may have role in sex determination. Forty four putative Genes involved in sex determination (GISD) reported in model plants were shortlisted from literature survey, and twenty nine seabuckthorn homologous sequences were identified from available seabuckthorn genomic resources. Of these, 21 genes were found to differentially express in either male or female flower bud stages. HrCRY2 was significantly expressed in female flower buds only while HrCO had significant expression in male flowers only. Among the three male and female floral development stages (FDS), male stage II had significant expression of most of the GISD. Information on these sex-specific expressed genes will help in elucidating sex determination mechanism in seabuckthorn.

Introduction

Hippophae rhamnoides commonly known as seabuckthorn belongs to the family Elaeagnaceae. Seabuckthorn berries are among the most nutritious and vitamin-rich fruits found in the plant kingdom. In general, the flesh of berries contains a diverse complex of vitamins, mineral substances such as sodium salts, potassium, calcium, carbohydrates, proteins, sugars and amino acids. [1, 2]. Moreover, the oil from the sea buckthorn berry contains on average 35% of the rare and valuable palmitoleic acid (16:1n-7; omega-7 series fatty acid) [3]. Seabuckthorn has a great potential for researchers in the field of biotechnology, neutraceutical and environmental sciences [4]. Various products had been developed from the berries of seabuckthorn such as oil, juice, alcoholic beverages, candies, ice-cream, tea, jam and biscuits. [3]. Thus the demand of seabuckthorn berries has increased in past few years due to their increased use in commercial products [5]. This increase in demand warrants its intensive cultivation, instead of just collection from wild resources and requires genetic improvement to increase its productivity and quality.

For development of superior seabuckthorn, breeding projects target both females and male cultivars [3]. Moreover, the objectives for breeding male and female plants vary, since there are extra quality criteria to be met in female plants, as berries occur on female plants only [3]. The success of the breeding program in dioecious plants depends upon early identification of progeny's gender. Unfortunately, gender of seabuckthorn seedlings cannot be determined morphologically until flowering, which usually occurs after 3–4 years in the field [6]. This represents a serious problem for plant breeders who are forced retain large number of male for several years. Much of the work and money could be saved if large proportion of the males could be discarded at an early stage in evaluation process.

In dioecious plants gender determination is regulated at genetic level by X/Y chromosome system [7]. Many molecular marker based studies like RAPD, SSR, ISSR, SCAR etc. were conducted for past several years for gender identification in seabuckthorn and molecular markers can distinguish male and female plants. [812]. However, none of the marker based studies in seabuckthorn related the markers with the mechanism governing sex determination. Therefore the mechanism governing the sex determination in seabuckthorn still remains unknown [3].

Differences between male and female plants are primarily detected in reproductive organs, which occur through differential growth, repression or abortion of sex organs in unisexual flowers [13, 14]. Various category of genes like floral mersitem identity genes, floral organ identity genes and flowering time genes play a major role as Genes involved in sex determination (GISD) in development of unisexual flowers [15, 16]. In case of Thalictrum dioicum, floral organ identity genes were differentially expressed in early development stages of male and female flowers. This led to the conclusion that regulation of these homeotic genes resulted in gender determination in this species [17]. Also the role of MADS box homeotic genes was analysed in male and female flowers of Hop (Humulus lupulus). Northern hybridisation in H. lupulus showed that M1 (DEFICIENS homologue) and M2 (Petunia FLORAL BINDING PROTEIN 2 homologue) transcripts were present in the early stages of floral development of both sexes, but at later stages, expression of both genes increased in male flowers and decreased in female flowers [18, 19]. Moreover, apart from floral regulatory genes, sex determination is also dependent upon the regulatory networks which alter sex expression based on environmental cues such as photoperiod and temperature [20].

The genetic control of sex determination is well-kown in several model plant systems like Silene latifolia [2123], Cucumis sativus [2426], Salix [27, 28], etc. Moreover, molecular and genetic studies showed that the underlying mechanisms controlling flower development are largely conserved in distantly related dicotyledonous plant species [29]. Thus, genomic resources generated from these model plants could be used to identify the potential GISD in seabuckthorn. A possible strategy to identify genes essential in a development process is to screen mRNAs that are present in one sample and absent (or rare) in other ones [30]. In order to identify mRNA transcripts involved in sex determination in dioecious plants like S. latifolia, Rumex acetosa, and Actinidia chinensis,. different spatial and temporal development stages of flower were used [3133]. Numerous flowering genes like APETALA 2, CLAVATA 1 and SEPTALA 3 showed differential expression among male and female flowers of plants like Z. mays, S. latifolia, A. Officinalis [3436], which indicated their role in sex determination in the above mentioned plants. Thus for identification of potential seabuckthorn GISD, differential expression of known flowering genes was analysed using quantitative Real Time PCR (qRT-PCR) in three temporal Floral Development Stages (FDS) of both male and female seabuckthorn flowers.

Material and Methods

Plant material, RNA extraction and cDNA synthesis

The flower buds of Hippophae rhamnoides collected from Defence Institute of High Altitude Research (DIHAR), J&K, India (Geographic Coordinates—34°08’ 236” N, 77° 34’ 345” E) were used in this investigation (Permission granted by Director, DIHAR, Leh, Jammu and Kashmir, India). Three different samples of floral buds for current study were collected on the basis of phenological observations on flowering of seabuckthorn in the region of study. Flower buds start developing from the month April and flowers open in the start of May to mid-May. The flower bud samples were collected in the month of April at ten days interval, starting from dormant winter bud to when buds are about to open. This is period when female and male reproductive tissues are formed in the flower buds. Flower buds were immediately frozen in liquid nitrogen and were stored at -80 °C till further use. Male and female flower bud stages were designated as Male Stage I (MST I), Male Stage II (MST II), Male Stage III (MST III) and Female Stage I (FST I), Female Stage II (FST II) and Female Stage III (FST III) respectively as shown in Fig 1. RNA was extracted from flower buds using Bangalore Genei Plant Total RNA extraction kit as per manufacturer instructions. RNA concentration was estimated by U.V. spectrophotometry and integrity was confirmed by electrophoresing samples on a 1.2% denaturing agarose gel. First strand of cDNA was synthesised from 1 μg of total RNA using Verso cDNA Kit (Thermo Scientific). The quality of cDNA was tested by amplifying 26S gene fragment using 26S primers (S1 Table) under following amplification conditions (95°C for 4 min and then 35 cycles at 95°C for 30 s, 55.5°C for 30 s and 72°C for 50 s) and products were electrophoresed in 1.8% agarose gel.

thumbnail
Fig 1. Temporal male and female floral bud development stages in seabuckthorn.

https://doi.org/10.1371/journal.pone.0124890.g001

Identification of seabuckthorn homologues of potential GISD and phylogenetic analysis

A literature survey was undertaken to short list genes involved in flower development of Arabidopsis which could be potential candidates for sex determination in seabuckthorn (Table 1). Nucleotide sequences of floral regulatory genes well-characterized in plants like Silene latifolia, Arabidopsis thaliana, Vitis vinifera, Cucumis sativa, etc. were downloaded from NCBI Genbank database in FASTA format (S1 File). The sequence data was manually curated and redundant sequences of the same species were discarded. Quality trimmed and filtered nucleotide sequences of seabuckthorn were retrieved from seed [37], root and leaf [38] transcriptome (NCBI Accession No.SRX118240, SRX131619 and SRX131618 respectively) and ESTs from NCBI EST database. A series of BLASTN analyses with default parameters identified broadly conserved sequences of potential GISD from seabuckthorn genomic resources showing syntenic relationship with known GISD (Table 2 and S2 File). BLASTN reports were analysed manually and the sequences (showing similarity with known GISD sequences) having e-value greater than 10–4 and query coverage less than 100bp were discarded. Homologous sequences of GISD having the lowest e-value were chosen for validation through qRT-PCR (Table 2). To further confirm the identity of the seabuckthorn sequences, domains and repeats were identified within the GISD sequences. Nucleotide sequences of putative seabuckthorn GISD were translated to amino acid coding sequences (S3 File) using ExPASy translate tool (http://web.expasy.org/translate/). The sequences with longest open reading frame were used for repeats, domain and protein family identification using EBI Interpro server (http://www.ebi.ac.uk/interpro/). For Phylogenetic reconstruction of potential GISD in seabuckthorn, protein sequences of known GISD characterized in model plant species were downloaded from NCBI Genbank database (S4 File). The alignment of the sequences was done with the help of CLUSTALX [39] and the final tree was constructed using MEGA 6 (Molecular Evolutionary Genetics Analysis 6.0) software [40].

thumbnail
Table 1. List of potential genes involved in sex determination in seabuckthorn.

https://doi.org/10.1371/journal.pone.0124890.t001

thumbnail
Table 2. List of potential Seabuckthorn GISD retrieved from available seabuckthorn resources.

https://doi.org/10.1371/journal.pone.0124890.t002

Expression analysis of GISD by qRT-PCR

Primers for candidate genes were designed using the Primer3 web application (http://bioinfo.ut.ee/primer3-0.4.0/), with Tm of 55–60°C and amplicon size between 100 bp and 250 bp (S1 Table). qRT-PCR was performed with duplicate amplifications using SYBR-green-based detection system (IQ SYBR Green Supermix (Biorad) in the Biorad CFX96 Real-Time PCR Detection System). The reactions contained 100 ng cDNA template and 0.5μM of primers in total volume of 13μl. Cycle parameters of reaction were 95°C for 3 min and then 39 cycles at 95°C for 10, 60°C for 30 s and 72 for 20s. Expression data were analysed with ΔΔCT method [41]. The expression of four internal reference genes namely ubiquitin, β-actin, 26S and GAPDH was checked on four floral bud samples. 26S and GAPDH genes showed consistent expression pattern in male and female flower bud stages (Unpublished data) and were used for gene expression data normalisation. The data presented in the figures and tables are based on the average of 2 PCR samples used from 3 biological samples. Fold expression of genes was calculated between the same development stages of male and female flowers. Heat map representing the gene expression data of GISD in three developmental stages of male and female seabuckthorn flowers was generated using the GENEX Ver. 6.0 software (http://genex.gene-quantification.info).

Results

Identification of seabuckthorn homologues of potential GISD and phylogenetic analysis

The current study was focused on 44 Arabidopsis genes that were known to be involved in floral regulatory pathways (Table 1) and could be probable candidates for sex determination in H. rhamnoides. Out of 44 Arabidopsis flowering genes, 24 genes had homologous sequences in available H. rhamnoides genomic resources (Table 2). Arabidopsis genes for which homologous sequences were not present in the transcriptome data of seabuckthorn include AP3, CAL, CRC, JAG, KNU, LFY, NZZ, NUB, RBE, SPL, SVP, SUP, WUS, FLC, FLT, UFO, FIM, ER and DAD1. The identified homologous sequences of seabuckthorn GISD were compared with similar genes of other plants species deposited in NCBI genebank nucleotide database as well as EST databases of other plant species like Actnidia chinesis [90]. Results of the analysis showed that the sequences of putative seabuckthorn GISD matched with transcripts of either one plant species or the other (Table 2). Also more than one copy of homologous sequences were found for genes CO (3), CRY1 (2), FRI (2) and TFL (2) in H. rhamnoides (Table 2). Domains and repeats found in all the homologues of seabuckthorn except for EF3, GI and NEF1 were similar to those present in Arabidopsis genes (Table 2). Such an outcome signifies that the identified contigs of putative seabuckthorn GISD have similar gene structure as of A. thaliana genes and thus are likely to perform identical functions as performed by respective genes in A. thaliana. Phylogenetic reconstruction of genes (Fig 2) showed that most of the seabuckthorn GISD clustered with similar genes except for HrCRY2, HrFIL, HrAP2 and HrNEF. Clustering of putative seabuckthorn GISD along with characterized genes in model plants further confirms that putative seabuckthorn genes share high homology to well characterized genes in model plants.

thumbnail
Fig 2. Phylogenetic tree of the potential GISD in Seabuckthorn (H. rhamnoides) based on the amino acid sequence alignment.

https://doi.org/10.1371/journal.pone.0124890.g002

Expression analysis of GISD by qRT PCR

The temporal expression of the 24 potential GISD and their additional homologues (Total of 29 candidate GISD) across three developmental stages of male and female flowers of H. rhamnoides was analysed (S2 Table). 21 GISD were analysed for differential expression among temporally corresponding male and female flower developmental stages (FDS). The CT values of eight GISD which showed values greater than 35, were not considered for further investigation. Seven GISD showed elevated expression in female FDS while fourteen GISD showed higher expression in male FDS, (Fig 3 and Table 3) details of which are given below.

thumbnail
Table 3. Relative fold expression of putative GISD within temporally corresponding male and female flower development stages.

https://doi.org/10.1371/journal.pone.0124890.t003

Floral meristem identity genes.

As data presented in Table 3 and Fig 3A demonstrates, meristem identity gene HrAP1 showed female specific expression. The expression of this gene was notably 1347 fold higher in FST II as compared to MST II. On the contrary, expression of gene HrAP2 was higher in all the male flower developmental stages withthe maximum differential expression being in MST II (7.70 fold) as compared to FST II. HrLFY and HrCLV1 showed stage specific expression in male and female FDS. The expression of HrLFY was notably higher in MST I (32.16 fold) and FST II (10.15 fold) as compared to their corresponding stages. HrCLV1 was significantly expressed in MST II (9.11 fold) and FST III (4.15 fold). On the basis of this data it is concluded that expression of gene HrAP1 is female specific while that of HrAP2 is male specific. However, the expression of gene HrLFY and HrCLV1 in male and female flowers was stage dependent.

Floral organ identity genes.

Among floral organ identity genes the expression of floral organ polarity gene HrFIL was higher in all FDS of female flowers (Fig 3B). The differential expression was notably wider in FST III vs MST III (53.88 fold). On the contrary the expression of HrYAB5 and HrSEP3 was higher in all male FDS with highest differential expression of 250 fold and 1000 fold were recorded in MST II respectively. Stamen and carpel identity gene AGAMOUS (HrAG) showed stage dependent expression pattern which was higher in FST I (12.55 fold) and MST II (6.34 fold) as compared to their corresponding stages. Thus from the data recorded it can be concluded that expression of floral organ identity gene HrSEP3 and HrYAB5 was higher in male FDS and that of HrFIL was higher in female FDS. Also the relative expression of gene HrAG was flower developmental stage dependent rather than sex of flower.

Flowering time regulation genes.

The expression of Blue-Ultraviolet A receptor gene CRYPTOCHROME2 (HrCRY2) was higher across all the female FDS as compared to male FDS. The expression of this gene was 129.3 fold higher. CRYPTOCHROME1(HrCRY1) was relatively expressed higher in all male FDS with MST II and MST III showing 6.6 fold and 2.33 fold higher expression as compared to corresponding female FDS (Fig 3C). Similarly the expression of far red light receptor gene PHYTOCHROME B (HrPHYB) was higher in all male FDS notably MST II and MST III, which showed 25 fold and 7.5 fold higher expression with respect to FST II and FST III respectively (Fig 3C). The expression of CONSTANS (HrCO) responsible for flowering in long days was higher in all male FDS (Fig 3C) (9.91 fold in MST I with respect to FST I, 30 fold in MST II with respect to FST II and 113 fold in MST III with respect to FST III). The second homologue of CO (HrCOLK) showed similar pattern of expression but relative difference in a expression level was less pronounced in male and female FDS as compared to HrCO. FRIGADIA (HrFRI) and its second hoologue HrFRILK responsible for delayed flowering in absence of cold temperatures, were also found to have elevated expression in male FDS as compared to their corresponding female FDS (HrFRI 17.31 fold higher in MST II; HrFRILK 9.78 fold higher in MST II) (Fig 3C). The relative expression of genes HrGI and HrEF1 was stage dependent. Thus it is concluded that expression of most of flowering time genes including HrCRY1, HrPHYB and HrCO was higher in all male FDS while that of HrCRY2 was higher in all female FDS.

Phytohormone ethylene response pathway genes.

The expression of seabuckthorn homologues of ethylene response pathway genes ETHYLENE RESPONSE SENSOR 1 (HrERS1) and ETHYLENE RECEPTOR 1 (HrETR1) was higher in male flowers. HrERS1 (Fig 3D) showed 10.86 fold higher expression in MST II with respect to FST II while expression of HrETR1 (Fig 3D) was recorded 9.78 fold higher in MST II with respect to FST II.

Pollen exine formation genes.

HrX1 is the female specific SCAR marker which was found to show high level of similarity to Acyl CoA synthatase and other related plant ligases on the basis of BLASTn and tBLASTx analysis of sequence [8]. HrX1 (Fig 3E) expressed 43.47 fold higher in MST II with respect to FSTII. On the contrary, expression of HrNEF1 (Fig 3E) was observed to be 19.01 fold higher in FST I as compared to MST I.

Floral development stage (FDS) specific expression of GISD

In stage I of flower development the expression of genes HrAP1, HrCRY2, HrEF1, HrNEF and HrAG was higher in female flowers while expression of genes HrAP2, HrLFY, HrFRI, and HrGI was higher in male flowers (Fig 4). The expression of all putative GISD except for HrCRY2 and HrLFY was higher in 2nd developmental stage of male flowers (Fig 5). In STAGE III FDS female flowers had higher expression of HrAP1, HrCRY2, HrEF1 and HrFILF while male flowers had higher expression levels of HrCRY1, HrCO and HrPHYB (Fig 6). Moreover, the heat map of putative GISD (Fig 7) shows that male flowers have maximum GISD with higher level of expression as compared to female flowers in 2nd floral developmental stage.

thumbnail
Fig 4. Comparative expression of seabuckthorn GISD between floral Development Stages (FDS)—Male stage I (MST I) vs Female Stage I (FST I).

https://doi.org/10.1371/journal.pone.0124890.g004

thumbnail
Fig 5. Comparative expression of seabuckthorn GISD between floral Development Stages (FDS)—Male stage II (MST II) vs Female Stage II (FST II).

https://doi.org/10.1371/journal.pone.0124890.g005

thumbnail
Fig 6. Comparative expression of seabuckthorn GISD between floral Development Stages (FDS)—Male stage III (MST III) vs Female Stage III (FST III).

https://doi.org/10.1371/journal.pone.0124890.g006

thumbnail
Fig 7. Heat map of relative expression of putative GISD in seabuckthorn FDS.

https://doi.org/10.1371/journal.pone.0124890.g007

Discussion.

The identified putative GISD of seabuckthorn shared the sequence similarity with plant species like R. hybrid, V. vinefera, M. notabilis, P. trichocarpa, etc. and had similar repeats and domains as Arabidopsis floral regulatory genes. Most of these genes also clustered along with well characterised genes of model plants. Identification of homologous flowering genes in seabuckthorn reflects that flowering pathways of seabuckthorn share similarity with Arabidopsis as well as other model dioecious plants. Thus as in the case of S. latifolia, R. acetosa and A. chinensis, genes involved in these flowering pathways could be potential candidates of sex determination in seabuckthorn.

Expression pattern of MADS box genes in male and female flowers of sorrel (Rumex acetosa) suggested that these genes could play an important role in sex determination [31]. CLASS A MADS box gene HrAP1 (Fig 3A) showed female specific expression while HrAP2 (Fig 3A), expressed particularly in male flowers (Fig 3A). AP2 plays an important role for sex determination in maize [35]. It suppresses carpel in tassel of male flowers by targeting TASSELSEED 4 (TSL4). Similarly expression of another floral meristem identity gene HrCLV1 (Fig 3A) was recorded highest in MST II. In case of S. latifolia, CLV1 triggers carpel suppression in male flowers [36]. Thus HrAP2 and HrCLV1 may be involved in determining meristem identity in male flowers while HrAP1 could trigger meristem development in female flowers of seabuckthorn.

Expression of floral organ identity gene HrSEP3 (Fig 3B) was recorded highest in MST II. Higher expression of SEP3 homologue was also observed in male flowers of Asparagus officinalis [34]. On the other hand HrAG (Fig 3B) showed significant expression in FST II. Thus HrSEP3 and HrAG may have a crucial role in establishing floral organ identity in male and female flowers respectively.

The expression of cryptochrome receptor gene HrCRY2 (Fig 3C) was higher in female flowers as compared to male flowers. On the other hand, level of expression of cryptochrome gene HrCRY1 (Fig 3C), phytochrome gene HrPHYB (Fig 3C) and circardian pathway gene HrCO (Fig 3C) was higher in all male flower development stages. In dioecious plants like S. latifolia and Populus tomentosa, male and female flowers develop at different time. The photoreceptor encoding genes like CRY1, CRY2, PHYA and PHYB regulate circardian pathway genes like CO, GI and FT and could alter flowering time depending upon external cues [91]. Differential expression of CRY1, CRY2, CO and GI was observed among male and female flowers of P. tomentosa and was correlated with asynchronous development of male and female flowers [92]. Thus expression pattern of flowering time genes showed that HrCRY2 could influence time-dependent development of female flowers while HrCRY1, HrPHYB and HrCO may affect temporal development of male flowers in seabuckthorn.

Phytohormone ethylene response genes HrERS1 and HrETR1 differentially expressed in all the stages of male and female flower but without bias of expression towards particular gender. Such an outcome could be expected because, in case of dioecious plants genetic variations have a more prominent role in gender determination than internal environment and environment variation. Expression of gene containing the female specific SCAR marker HrX1 was higher in male flowers as compared to female flowers. HrX1 shares sequence similarity with known plant ligases such as acyl Coa synthatase [8]. In A. thaliana knocking out of acyl CoA synthatase led to production of unviable pollen, which in turn produced male sterile plants [93]. Another pollen exine gene HrNEF1 showed higher expression in female flowers as compared to male flowers. Disruption of NEF1 in A. thaliana affected lipid accumulation in the plastids of tapetum as well as exine formation of pollen, thus resulted in male sterility in A. thaliana [89]. Thus expression pattern of HrX1 and HrNEF1 suggested that these genes could play an important in sex determination in seabuckthorn.

The expression of genes varied throughout the development of flowers in both male and female flowers of seabuckthorn. Out of the three developmental stages, 2nd stage had the maximum number of genes with expression biased towards male flowers (Fig 4 and Fig 5 & Fig 6). Thus stage II of male and female flowers require further investigation to justify the tilt of GISD expression towards male flowers.

In conclusion, the current study showed differential expression of putative seabuckthorn GISD in all the three floral developmental stages of both male and female flowers. The expression level of HrCO gene was observed higher in the developmental stages of male flowers as compared to female flowers. Whereas HrCRY2 gene significantly showed higher expression levels in the female floral developmental stages only. Further investigation is required to understand the role of HrCO and HrCRY2 genes in development of male and female flowers respectively.

Supporting Information

S1 File. Nucleotide sequences of the known flowering pathway genes in model plants.

https://doi.org/10.1371/journal.pone.0124890.s001

(RAR)

S2 File. Nucleotide sequences of seabuckthorn putative GISD retrieved from seabuckthorn genomic resources.

https://doi.org/10.1371/journal.pone.0124890.s002

(TXT)

S3 File. Protein sequences of seabuckthorn putative GISD.

https://doi.org/10.1371/journal.pone.0124890.s003

(TXT)

S4 File. Protein sequences of the known flowering pathway genes in model plants.

https://doi.org/10.1371/journal.pone.0124890.s004

(TXT)

S1 Table. List of primers used in qRT PCR analysis of putative GISD.

https://doi.org/10.1371/journal.pone.0124890.s005

(DOCX)

S2 Table. Normalized expression values of seabuckthorn putative GISD in three temporal developmental stages of male and female flowers.

https://doi.org/10.1371/journal.pone.0124890.s006

(DOCX)

Acknowledgments

We sincerely thank Defence Research and Development Organisation (DRDO), Ministry of Defence, Govt. of India for funding the current study. We also thank Jaypee University of Information and Technology, Solan (India) for providing facilities for conducting the research.

Author Contributions

Conceived and designed the experiments: AC TS RBS RSC AK. Performed the experiments: AC. Analyzed the data: AC VJ AK. Contributed reagents/materials/analysis tools: AC TS RBS AK. Wrote the paper: AC AK. In Silico analysis of genes: AC AK. Mining of putative GISD from seabuckthorn transcriptome: AC VJ.

References

  1. 1. Fan J, Ding X, Gu W. Radical-scavenging proanthocyanidins from sea buckthorn seed. Food chemistry. 2007;102(1):168–77.
  2. 2. Li C, Xu G, Zang R, Korpelainen H, Berninger F. Sex-related differences in leaf morphological and physiological responses in Hippophae rhamnoides along an altitudinal gradient. Tree Physiology. 2007;27(3):399–406. pmid:17241981
  3. 3. Kalia RK, Singh R, Rai MK, Mishra GP, Singh SR, Dhawan A. Biotechnological interventions in sea buckthorn (Hippophae L.): current status and future prospects. Trees. 2011;25(4):559–75.
  4. 4. Stobdan T, Angchuk D, Singh SB. Seabuckthorn: An emerging storehouse for researchers in India. Bangalore, India: Indian Academy of Sciences 2008.
  5. 5. Zubarev YA. Commercial cultivation of Seabuckthorn in Western Siberia, Russia. Seabuckthorn (Hippophae L): A Multipurpose Wonder Plant. 2008;3:49–60.
  6. 6. Gupta SM, Grover A, Pandey P, Ahmed Z. Female plants of Hippophae salicifolia D. Don are more responsive to cold stress than male plants. Physiology and Molecular Biology of Plants. 2012;18(4):377–80. pmid:24082501
  7. 7. Charlesworth D. Plant sex determination and sex chromosomes. Heredity. 2002;88(2):94–101. pmid:11932767
  8. 8. Chawla A, Kant A, Stobdan T, Srivastava RB, Chauhan R. Cross-species application of sex linked markers in H. salicifolia and H. tibetana. Scientia Horticulturae. 2014;170:281–3.
  9. 9. Jain A, Ghangal R, Grover A, Raghuvanshi S, Sharma PC. Development of EST-based new SSR markers in seabuckthorn. Physiology and Molecular Biology of Plants. 2010;16(4):375–8. pmid:23572988
  10. 10. Korekar G, Sharma R, Kumar R, Meenu , Bisht N, Srivastava R, et al. Identification and validation of sex-linked SCAR markers in dioecious Hippophae rhamnoides L. (Elaeagnaceae). Biotechnol Lett. 2012;34(5):973–8. pmid:22245920
  11. 11. Persson HA, Nybom H. Genetic sex determination and RAPD marker segregation in the dioecious species sea buckthorn (Hippophae rhamnoides L.). Hereditas. 1998;129(1):45–51.
  12. 12. Sharma A, Zinta G, Rana S, Shirko P. Molecular identification of sex in Hippophae rhamnoides L. using isozyme and RAPD markers. Forestry Studies in China. 2010;12(2):62–6.
  13. 13. Chuck G. Molecular mechanisms of sex determination in monoecious and dioecious plants. Advances in botanical research. 2010;54:53–83.
  14. 14. Matsunaga S. Sex chromosome-linked genes in plants. Genes and Genetic Systems. 2006;81(4):219. pmid:17038793
  15. 15. Iovene M, Yu Q, Ming R, Jiang J. Evidence for Emergence of Sex-Determining Gene (s) in a Centromeric Region in Vasconcellea parviflora. Genetics. 2014:genetics. 114.173021.
  16. 16. Miko I. Sex chromosomes and sex determination. Nature Education. 2008;1(1):108.
  17. 17. Di Stilio VS, Kramer EM, Baum DA. Floral MADS box genes and homeotic gender dimorphism in Thalictrum dioicum (Ranunculaceae)–a new model for the study of dioecy. The Plant Journal. 2005;41(5):755–66. pmid:15703062
  18. 18. Shepard H, Parker J, Darby P, Ainsworth CC. Sex expression in hop (Humulus lupulus L. and H. japonicus Sieb. et Zucc.): floral morphology and sex chromosomes. Sex Determination in Plants. 1999;8:139–50.
  19. 19. Shephard HL, Ainsworth CC. MADS box gene expression during flower development in hop (Humulus lupulus L.) 1998 [updated 22-7-2014; cited 2014 14 October]. 70–1]. Available: http://extras.springer.com/1999/978-3-642-64248-7/ebo1998/meeting/york98/p09.pdf.
  20. 20. Zhao H, Li Y, Duan B, Korpelainen H, Li C. Sex‐related adaptive responses of Populus cathayana to photoperiod transitions. Plant, cell & environment. 2009;32(10):1401–11.
  21. 21. Delph LF, Arntz AM, Scotti‐Saintagne C, Scotti I. The genomic architecture of sexual dimorphism in the dioecious plant Silene latifolia. Evolution. 2010;64(10):2873–86. pmid:20550575
  22. 22. Lebel-Hardenack S, Hauser E, Law TF, Schmid J, Grant SR. Mapping of sex determination loci on the white campion (Silene latifolia) Y chromosome using amplified fragment length polymorphism. Genetics. 2002;160(2):717–25. pmid:11861573
  23. 23. Matsunaga S, Kawano S, Takano H, Uchida H, Sakai A, Kuroiwa T. Isolation and developmental expression of male reproductive organ‐specific genes in a dioecious campion, Melandrium album (Silene latifolia). The Plant Journal. 1996;10(4):679–89. pmid:8893544
  24. 24. Adhikari S, Bandyopadhyay TK, Ghosh P. Hormonal control of sex expression of cucumber (Cucumis sativus L.) with the identification of sex linked molecular marker. The Nucleus. 2012;55(2):115–22.
  25. 25. Foucart C, Boualem A, Lasseur B, Eleblu J, Fahraj I, Bendahmane A. Sex determination in cucurbits. Biologie aujourd'hui. 2012;206(1):57. pmid:22463996
  26. 26. Wu T, Qin Z, Zhou X, Feng Z, Du Y. Transcriptome profile analysis of floral sex determination in cucumber. Journal of plant physiology. 2010;167(11):905–13. pmid:20303197
  27. 27. Liu J, Yin T, Ye N, Chen Y, Yin T, Liu M, et al. Transcriptome Analysis of the Differentially Expressed Genes in the Male and Female Shrub Willows (Salix suchowensis). PloS one. 2013;8(4):e60181. pmid:23560075
  28. 28. Semerikov V, Lagercrantz U, Tsarouhas V, Rönnberg-Wästljung A, Alström-Rapaport C, Lascoux M. Genetic mapping of sex-linked markers in Salix viminalis L. Heredity. 2003;91(3):293–9. pmid:12939631
  29. 29. Yanofsky MF. Floral meristems to floral organs: genes controlling early events in Arabidopsis flower development. Annual review of plant biology. 1995;46(1):167–88.
  30. 30. Perl-Treves R. Male to female conversion along the cucumber shoot: approaches to studying sex genes and floral development in Cucumis sativus: Oxford, UK: Bios Scientific Publishers; 1999.
  31. 31. Ainsworth C, Crossley S, Buchanan-Wollaston V, Thangavelu M, Parker J. Male and female flowers of the dioecious plant sorrel show different patterns of MADS box gene expression. The Plant Cell Online. 1995;7(10):1583–98. pmid:7580253
  32. 32. Hardenack S, Ye D, Saedler H, Grant S. Comparison of MADS box gene expression in developing male and female flowers of the dioecious plant white campion. The Plant Cell Online. 1994;6(12):1775–87. pmid:7866023
  33. 33. Kim HB, Jun S-S, Choe S, Cho JY, Choi S-B, Kim S-C. Identification of differentially expressed genes from male and female flowers of kiwifruit. African Journal of Biotechnology. 2010;9(40):6684–94.
  34. 34. Caporali E, Spada A, Losa A, Marziani G. The MADS box gene AOM1 is expressed in reproductive meristems and flowers of the dioecious species Asparagus officinalis. Sexual Plant Reproduction. 2000;13(3):151–6.
  35. 35. Chuck G, Meeley R, Irish E, Sakai H, Hake S. The maize tasselseed4 microRNA controls sex determination and meristem cell fate by targeting Tasselseed6/indeterminate spikelet1. Nature genetics. 2007;39(12):1517–21. pmid:18026103
  36. 36. Koizumi A, Yamanaka K, Nishihara K, Kazama Y, Abe T, Kawano S. Two separate pathways including SlCLV1, SlSTM and SlCUC that control carpel development in a bisexual mutant of Silene latifolia. Plant and cell physiology. 2010;51(2):282–93. pmid:20064843
  37. 37. Fatima T, Snyder CL, Schroeder WR, Cram D, Datla R, Wishart D, et al. Fatty Acid Composition of Developing Sea Buckthorn (Hippophae rhamnoides L.) Berry and the Transcriptome of the Mature Seed. PLoS One. 2012;7(4):e34099. pmid:22558083
  38. 38. Ghangal R, Chaudhary S, Jain M, Purty RS, Sharma PC. Optimization of De Novo Short Read Assembly of Seabuckthorn (Hippophae rhamnoides L.) Transcriptome. PLoS One. 2013;8(8):e72516. pmid:23991119
  39. 39. Larkin MA, Blackshields G, Brown N, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8. pmid:17846036
  40. 40. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Molecular biology and evolution. 2013;30(12):2725–9. pmid:24132122
  41. 41. Livak KJ, Schmittgen TD. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2− ΔΔCT Method. methods. 2001;25(4):402–8. pmid:11846609
  42. 42. Huijser P, Klein J, Lönnig W, Meijer H, Saedler H, Sommer H. Bracteomania, an inflorescence anomaly, is caused by the loss of function of the MADS-box gene squamosa in Antirrhinum majus. The EMBO journal. 1992;11(4):1239. pmid:1563342
  43. 43. Irish VF, Sussex IM. Function of the apetala-1 gene during Arabidopsis floral development. The Plant Cell Online. 1990;2(8):741–53. pmid:1983792
  44. 44. Mandel MA, Gustafson-Brown C, Savidge B, Yanofsky MF. Molecular characterization of the Arabidopsis floral homeotic gene APETALA1. 1992.
  45. 45. Krizek BA, Fletcher JC. Molecular mechanisms of flower development: an armchair guide. Nature Reviews Genetics. 2005;6(9):688–98. pmid:16151374
  46. 46. Bradley D, Carpenter R, Sommer H, Hartley N, Coen E. Complementary floral homeotic phenotypes result from opposite orientations of a transposon at the plena locus of Antirrhinum. Cell. 1993;72(1):85–95. pmid:8093684
  47. 47. Yanofsky MF, Ma H, Bowman JL, Drews GN, Feldmann KA, Meyerowitz EM. The protein encoded by the Arabidopsis homeotic gene agamous resembles transcription factors. Nature. 1990;346(6279):35–9. pmid:1973265
  48. 48. Bowman JL, Alvarez J, Weigel D, Meyerowitz EM, Smyth DR. Control of flower development in Arabidopsis thaliana by APETALA 1 and interacting genes. DEVELOPMENT-CAMBRIDGE-. 1993;119:721–.
  49. 49. Alvarez J, Smyth DR. CRABS CLAW and SPATULA, two Arabidopsis genes that control carpel development in parallel with AGAMOUS. Development. 1999;126(11):2377–86. pmid:10225997
  50. 50. Clark SE, Williams RW, Meyerowitz EM. The CLAVATA1 Gene Encodes a Putative Receptor Kinase That Controls Shoot and Floral Meristem Size in Arabidopsis. Cell. 1997;89(4):575–85. pmid:9160749
  51. 51. Samach A, Onouchi H, Gold SE, Ditta GS, Schwarz-Sommer Z, Yanofsky MF, et al. Distinct roles of CONSTANS target genes in reproductive development of Arabidopsis. Science. 2000;288(5471):1613–6. pmid:10834834
  52. 52. Guo H, Yang H, Mockler TC, Lin C. Regulation of flowering time by Arabidopsis photoreceptors. Science. 1998;279(5355):1360–3. pmid:9478898
  53. 53. Erfle H, Ventzki R, Voss H, Rechmann S, Benes V, Stegemann J, et al. Sequence and analysis of chromosome 3 of the plant Arabidopsis thaliana. Nature. 2000;408(6814):820–2. pmid:11130713
  54. 54. Nole-Wilson S, Krizek BA. AINTEGUMENTA contributes to organ polarity and regulates growth of lateral organs in combination with YABBY genes. Plant Physiology. 2006;141(3):977–87. pmid:16714408
  55. 55. Dinneny JR, Yadegari R, Fischer RL, Yanofsky MF, Weigel D. The role of JAGGED in shaping lateral organs. Development. 2004;131(5):1101–10. pmid:14973282
  56. 56. Payne T, Johnson SD, Koltunow AM. KNUCKLES (KNU) encodes a C2H2 zinc-finger protein that regulates development of basal pattern elements of the Arabidopsis gynoecium. Development. 2004;131(15):3737–49. pmid:15240552
  57. 57. Coen ES, Romero J, Doyle S, Elliott R, Murphy G, Carpenter R. Floricaula: A homeotic gene required for flower development in Antirrhinum majus. Cell. 1990;63(6):1311–22. pmid:1702033
  58. 58. Reeves PH, Coupland G. Response of plant development to environment: control of flowering by daylength and temperature. Current opinion in plant biology. 2000;3(1):37–42. pmid:10679453
  59. 59. Schiefthaler U, Balasubramanian S, Sieber P, Chevalier D, Wisman E, Schneitz K. Molecular analysis of NOZZLE, a gene involved in pattern formation and early sporogenesis during sex organ development in Arabidopsis thaliana. Proceedings of the National Academy of Sciences. 1999;96(20):11664–9. pmid:10500234
  60. 60. Dinneny JR, Weigel D, Yanofsky MF. NUBBIN and JAGGED define stamen and carpel shape in Arabidopsis. Development. 2006;133(9):1645–55. pmid:16554365
  61. 61. Tröbner W, Ramirez L, Motte P, Hue I, Huijser P, Lönnig W, et al. GLOBOSA: a homeotic gene which interacts with DEFICIENS in the control of Antirrhinum floral organogenesis. The EMBO Journal. 1992;11(13):4693. pmid:1361166
  62. 62. Krizek BA, Lewis MW, Fletcher JC. RABBIT EARS is a second‐whorl repressor of AGAMOUS that maintains spatial boundaries in Arabidopsis flowers. The Plant Journal. 2006;45(3):369–83. pmid:16412084
  63. 63. Takeda S, Matsumoto N, Okada K. RABBIT EARS, encoding a SUPERMAN-like zinc finger protein, regulates petal development in Arabidopsis thaliana. Development. 2004;131(2):425–34. pmid:14681191
  64. 64. Yang W-C, Ye D, Xu J, Sundaresan V. The SPOROCYTELESS gene of Arabidopsis is required for initiation of sporogenesis and encodes a novel nuclear protein. Genes & Development. 1999;13(16):2108–17.
  65. 65. Moon J, Suh SS, Lee H, Choi KR, Hong CB, Paek NC, et al. The SOC1 MADS‐box gene integrates vernalization and gibberellin signals for flowering in Arabidopsis. The Plant Journal. 2003;35(5):613–23. pmid:12940954
  66. 66. Lee JH, Yoo SJ, Park SH, Hwang I, Lee JS, Ahn JH. Role of SVP in the control of flowering time by ambient temperature in Arabidopsis. Genes & Development. 2007;21(4):397–402.
  67. 67. Bowman JL, Sakai H, Jack T, Weigel D, Mayer U, Meyerowitz EM. SUPERMAN, a regulator of floral homeotic genes in Arabidopsis. Development. 1992;114(3):599–615. pmid:1352237
  68. 68. Gaiser JC, Robinson-Beers K, Gasser CS. The Arabidopsis SUPERMAN gene mediates asymmetric growth of the outer integument of ovules. The Plant Cell Online. 1995;7(3):333–45. pmid:12242374
  69. 69. Sakai H, Krizek BA, Jacobsen SE, Meyerowitz EM. Regulation of SUP expression identifies multiple regulators involved in Arabidopsis floral meristem development. The Plant Cell Online. 2000;12(9):1607–18. pmid:11006335
  70. 70. Sakai H, Medrano LJ, Meyerowitz EM. Role of SUPERMAN in maintaining Arabidopsis floral whorl boundaries. 1995.
  71. 71. Schultz EA, Pickett FB, Haughn GW. The FLO10 gene product regulates the expression domain of homeotic genes AP3 and PI in Arabidopsis flowers. The Plant Cell Online. 1991;3(11):1221–37. pmid:12324589
  72. 72. Alvarez J, Guli CL, Yu XH, Smyth DR. terminal flower: a gene affecting inflorescence development in Arabidopsis thaliana. The Plant Journal. 1992;2(1):103–16.
  73. 73. Esumi T, Tao R, Yonemori K. Isolation of LEAFY and TERMINAL FLOWER 1 homologues from six fruit tree species in the subfamily Maloideae of the Rosaceae. Sexual plant reproduction. 2005;17(6):277–87.
  74. 74. Laux T, Mayer K, Berger J, Jurgens G. The WUSCHEL gene is required for shoot and floral meristem integrity in Arabidopsis. Development. 1996;122(1):87–96. pmid:8565856
  75. 75. He Y, Michaels SD, Amasino RM. Regulation of flowering time by histone acetylation in Arabidopsis. Science. 2003;302(5651):1751–4. pmid:14593187
  76. 76. Kardailsky I, Shukla VK, Ahn JH, Dagenais N, Christensen SK, Nguyen JT, et al. Activation tagging of the floral inducer FT. Science. 1999;286(5446):1962–5. pmid:10583961
  77. 77. Johnson E, Bradley M, Harberd NP, Whitelam GC. Photoresponses of light-grown phyA mutants of Arabidopsis (phytochrome A is required for the perception of daylength extensions). Plant Physiology. 1994;105(1):141–9. pmid:12232194
  78. 78. Neff MM, Chory J. Genetic interactions between phytochrome A, phytochrome B, and cryptochrome 1 during Arabidopsis development. Plant Physiology. 1998;118(1):27–35. pmid:9733523
  79. 79. Reed JW, Nagatani A, Elich TD, Fagan M, Chory J. Phytochrome A and phytochrome B have overlapping but distinct functions in Arabidopsis development. Plant Physiology. 1994;104(4):1139–49. pmid:12232154
  80. 80. Sandhu KS, Hagely K, Neff MM. Genetic interactions between brassinosteroid-inactivating P450s and photomorphogenic photoreceptors in Arabidopsis thaliana. G3: Genes| Genomes| Genetics. 2012;2(12):1585–93. pmid:23275881
  81. 81. Halliday KJ, Koornneef M, Whitelam GC. Phytochrome B and at least one other phytochrome mediate the accelerated flowering response of Arabidopsis thaliana L. to low red/far-red ratio. Plant Physiology. 1994;104(4):1311–5. pmid:12232170
  82. 82. Ray A, Lang JD, Golden T, Ray S. SHORT INTEGUMENT (SIN1), a gene required for ovule development in Arabidopsis, also controls flowering time. Development. 1996;122(9):2631–8. pmid:8787738
  83. 83. Kania T, Russenberger D, Peng S, Apel K, Melzer S. FPF1 promotes flowering in Arabidopsis. The Plant Cell Online. 1997;9(8):1327–38. pmid:9286110
  84. 84. Ingram GC, Goodrich J, Wilkinson MD, Simon R, Haughn GW, Coen ES. Parallels between UNUSUAL FLORAL ORGANS and FIMBRIATA, genes controlling flower development in Arabidopsis and Antirrhinum. The Plant Cell Online. 1995;7(9):1501–10. pmid:8589630
  85. 85. Simon R, Carpenter R, Doyle S, Coen E. Fimbriata controls flower development by mediating between meristem and organ identity genes. Cell. 1994;78(1):99–107. pmid:8033217
  86. 86. Torii KU, Mitsukawa N, Oosumi T, Matsuura Y, Yokoyama R, Whittier RF, et al. The Arabidopsis ERECTA gene encodes a putative receptor protein kinase with extracellular leucine-rich repeats. The Plant Cell Online. 1996;8(4):735–46. pmid:8624444
  87. 87. Ishiguro S, Kawai-Oda A, Ueda J, Nishida I, Okada K. The DEFECTIVE IN ANTHER DEHISCENCE1 gene encodes a novel phospholipase A1 catalyzing the initial step of jasmonic acid biosynthesis, which synchronizes pollen maturation, anther dehiscence, and flower opening in Arabidopsis. The Plant Cell Online. 2001;13(10):2191–209. pmid:11595796
  88. 88. Yamasaki S, Fujii N, Takahashi H. The ethylene-regulated expression of CS-ETR2 and CS-ERS genes in cucumber plants and their possible involvement with sex expression in flowers. Plant and Cell Physiology. 2000;41(5):608–16. pmid:10929944
  89. 89. Ariizumi T, Hatakeyama K, Hinata K, Inatsugi R, Nishida I, Sato S, et al. Disruption of the novel plant protein NEF1 affects lipid accumulation in the plastids of the tapetum and exine formation of pollen, resulting in male sterility in Arabidopsis thaliana. The Plant Journal. 2004;39(2):170–81. pmid:15225283
  90. 90. Fraser LG, Tsang GK, Datson PM, De Silva HN, Harvey CF, Gill GP, et al. A gene-rich linkage map in the dioecious species Actinidia chinensis (kiwifruit) reveals putative X/Y sex-determining chromosomes. BMC genomics. 2009;10(1):102.
  91. 91. Valverde F, Mouradov A, Soppe W, Ravenscroft D, Samach A, Coupland G. Photoreceptor regulation of CONSTANS protein in photoperiodic flowering. Science. 2004;303(5660):1003–6. pmid:14963328
  92. 92. Song Y, Ma K, Ci D, Chen Q, Tian J, Zhang D. Sexual dimorphic floral development in dioecious plants revealed by transcriptome, phytohormone, and DNA methylation analysis in Populus tomentosa. Plant molecular biology. 2013;83(6):559–76. pmid:23860796
  93. 93. de Azevedo Souza C, Kim SS, Koch S, Kienow L, Schneider K, McKim SM, et al. A novel fatty Acyl-CoA Synthetase is required for pollen development and sporopollenin biosynthesis in Arabidopsis. The Plant Cell Online. 2009;21(2):507–25. pmid:19218397