Chromosome-level genome assembly of Hippophae gyantsensis

Hippophae gyantsensis, which is a native tree species in China, is ideal for windbreak and sand-fixing forests. It is an economically and ecologically valuable tree species distributed exclusively in the Qinghai-Tibet Plateau in China. In our study, we assembled a chromosome-level genome of H. gyantsensis using Illumina sequencing, Nanopore sequencing and chromosome structure capture technique. The genome was 716.32 Mb in size with scaffold N50 length of 64.84 Mb. A total of 716.25 Mb genome data was anchored and orientated onto 12 chromosomes with a mounting rate of up to 99.99%. Additionally, the genome was found to comprise approximately 56.84% repeat sequences, of which long terminal repeats(LTRs) that accounted for 33.19% of the entire genome. Meanwhile, a total of 32,316 protein-coding genes were predicted, and 91.07% of these genes were functionally annotated. We also completed a series of comparative genomic analyses to provide researchers with useful reference material for future studies on seabuckthorn.

The Systematics, Reproductive Biology, Biochemistry, and Breeding of Sea Buckthorn-A Review

Both the fruit flesh and seeds of sea buckthorn have multiple uses for medicinal and culinary purposes, including the valuable market for supplementary health foods. Bioactive compounds, such as essential amino acids, vitamins B, C, and E, carotenoids, polyphenols, ursolic acid, unsaturated fatty acids, and other active substances, are now being analyzed in detail for their medicinal properties. Domestication with commercial orchards and processing plants is undertaken in many countries, but there is a large need for improved plant material with high yield, tolerance to environmental stress, diseases, and pests, suitability for efficient harvesting methods, and high contents of compounds that have medicinal and/or culinary values. Applied breeding is based mainly on directed crosses between different subspecies of Hippophae rhamnoides. DNA markers have been applied to analyses of systematics and population genetics as well as for the discrimination of cultivars, but very few DNA markers have as yet been developed for use in selection and breeding. Several key genes in important metabolic pathways have, however, been identified, and four genomes have recently been sequenced.

Genome of Hippophae rhamnoides provides insights into a conserved molecular mechanism in actinorhizal and rhizobial symbioses

Sea buckthorn (Hippophae rhamnoides), a horticulturally multipurpose species in the family Elaeagnaceae, can build associations with Frankia actinomycetes to enable symbiotic nitrogen-fixing. Currently, no high-quality reference genome is available for an actinorhizal plant, which greatly hinders the study of actinorhizal symbiotic nodulation. Here, by combining short-read, long-read and Hi-C sequencing technologies, we generated a chromosome-level reference genome of H. rhamnoides (scaffold N50: 65 Mb, and genome size: 730 Mb) and predicted 30 812 protein-coding genes mainly on 12 pseudochromosomes. Hippophae rhamnoides was found to share a high proportion of symbiotic nodulation genes with Medicago truncatula, implying a shared molecular mechanism between actinorhizal and rhizobial symbioses. Phylogenetic analysis clustered the three paralogous NODULE INCEPTION (NIN) genes of H. rhamnoides with those of other nodulating species, forming the NIN group that most likely evolved from the ancestral NLP group. The genome of H. rhamnoides will help us to decipher the underlying genetic programming of actinorhizal symbiosis, and our high-quality genome and transcriptomic resources will make H. rhamnoides a new excellent model plant for actinorhizal symbiosis research.

Genome sequence and population genomics provide insights into chromosomal evolution and phytochemical innovation of Hippophae rhamnoides

Plants of the Elaeagnaceae family are widely used to treat various health disorders owing to their natural phytochemicals. Seabuckthorn (Hippophae rhamnoides L.) is an economically and ecologically important species within the family with richness of biologically and pharmacologically active substances. Here, we present a chromosome-level genome assembly of seabuckthorn (http://hipp.shengxin.ren/), the first genome sequence of Elaeagnaceae, which has a total length of 849.04 Mb with scaffold N50 of 69.52 Mb and 30 864 annotated genes. Two sequential tetraploidizations with one occurring ~36-41 million years ago (Mya) and the last ~24-27 Mya were inferred, resulting in expansion of genes related to ascorbate and aldarate metabolism, lipid biosynthesis, and fatty acid elongation. Comparative genomic analysis reconstructed the evolutionary trajectories of the seabuckthorn genome with the predicted ancestral genome of 14 proto-chromosomes. Comparative transcriptomic and metabonomic analyses identified some key genes contributing to high content of polyunsaturated fatty acids and ascorbic acid (AsA). Additionally, we generated and analysed 55 whole-genome sequences of diverse accessions, and identified 9.80 million genetic variants in the seabuckthorn germplasms. Intriguingly, genes in selective sweep regions identified through population genomic analysis appeared to contribute to the richness of AsA and fatty acid in seabuckthorn fruits, among which GalLDH, GMPase and ACC, TER were the potentially major-effect causative genes controlling AsA and fatty acid content of the fruit, respectively. Our research offers novel insights into the molecular basis underlying phytochemical innovation of seabuckthorn, and provides valuable resources for exploring the evolution of the Elaeagnaceae family and molecular breeding.

Genome-wide DNA N6-adenine methylation in sea buckthorn (Hippophae rhamnoides L.) fruit development

As a new epigenetic mark, DNA N6-adenine (6mA) methylation plays an important role in various biological processes and has been reported in many prokaryotic organisms in recent years. However, the distribution patterns and functions of DNA 6mA modification have been poorly studied in non-model crops. In this study, we observed that the methylation ratio of 6mA was about 0.016% in the sea buckthorn (Hippophae rhamnoides L.) genome using mass spectrometry. We first constructed a comprehensive 6mA landscape in sea buckthorn genome using nanopore sequencing at single-base resolution. Distribution analysis suggested that 6mA methylated sites were widely distributed in the sea buckthorn chromosomes, which were similar to those in Arabidopsis and rice. Furthermore, reduced 6mA DNA methylation is associated with different expression of genes related to the fruit-ripening process in sea buckthorn. Our results revealed that 6mA DNA modification could be considered an important epigenomic mark and contributes to the fruit ripening process in plants.

Transcriptome and DNA methylome provide insights into the molecular regulation of drought stress in sea buckthorn

Sea buckthorn is a typical drought-resistant tree species. However, there is a general lack of understanding of the pattern of DNA methylation linked with sea buckthorn responses to drought, and its relationship with drought tolerance mechanisms. In this study, we performed whole-transcriptome RNA sequencing and methylome sequencing in response to drought stress to explore differentially expressed mRNAs, miRNAs, lncRNAs and circRNAs in sea buckthorn leaves. Based on predicted DE pairs, we constructed a competitive endogenous RNA network, which revealed potential transcriptional regulatory roles in response to drought stress. The results of methylome sequencing revealed that the DNA methylation level was increased in sea buckthorn leaves under drought stress. We identified 13,405 differentially methylated regions between CK and TR. We found one DMR-associated DEG (Vacuolar-sorting receptor 6) involved in the ABA accumulation pathway. In addition, two DNA methyltransferases (HrMET1 and HrDRM1) were closely associated with drought-induced hypermethylation in sea buckthorn. Together, we firstly conducted a comprehensive transcriptomic and epigenetic analysis of sea buckthorn under drought stress, providing a resource for further study of the potential functions of genes, miRNAs, lncRNAs, circRNAs and DNA methyltransferases.

HrTCP20 dramatically enhance drought tolerance of sea buckthorn (Hippophae rhamnoides L). by mediating the JA signaling pathway

Sea buckthorn, an important ecological and economical tree species, have remarkable drought and salt resistance. The plant-specific transcription factor TCPs play important roles in plant growth, development, and stress responses. However, in sea buckthorn, the molecular mechanism of TCP proteins and their involvement in drought stress are unknown. Here, we found that the expression of HrTCP20 was significantly up-regulated in sea buckthorn under drought stress. Overexpression of HrTCP20 in Arabidopsis thaliana showed that the superoxide dismutase (SOD), polyphenol oxidase (POD), and chlorophyll (SPAD) content was significantly increased by 1.37 and 1.35 times. However, the malondialdehyde (MDA) content decreased by 0.51 times. Our studies further confirmed that silencing HrTCP20 by virus-induced gene silencing (VIGS) led to a decrease in the content of defense enzymes, relative water content (RWC), and an increase of relative electrical conductivity (REC). Silencing HrTCP20 also caused the jasmonic acid (JA) content to decrease in the VIGS-treated tree. Interestingly, we found that JA accumulation content and the expression of HrLOX2, an essential enzyme for JA synthesis, was significantly inhibited in HrTCP20-silenced sea buckthorn under drought stress. Yeast two-hybrid analysis also showed that HrTCP20 is directly bound to HrLOX2. Taken together, the HrTCP20 transcription factor was a positive regulator in drought resistance of sea buckthorn. Further, our findings will provide comprehensive insights into the forest tree defence system of drought stress.

Oligo-FISH Can Identify Chromosomes and Distinguish Hippophaë rhamnoides L. Taxa

Oligo-fluorescence in situ hybridization (FISH) facilitates precise chromosome identification and comparative cytogenetic analysis. Detection of autosomal chromosomes of Hippophaë rhamnoides has not been achieved using oligonucleotide sequences. Here, the chromosomes of five H. rhamnoides taxa in the mitotic metaphase and mitotic metaphase to anaphase were detected using the oligo-FISH probes (AG3T3)3, 5S rDNA, and (TTG)6. In total, 24 small chromosomes were clearly observed in the mitotic metaphase (0.89–3.03 μm), whereas 24–48 small chromosomes were observed in the mitotic metaphase to anaphase (0.94–3.10 μm). The signal number and intensity of (AG3T3)3, 5S rDNA, and (TTG)6 in the mitotic metaphase to anaphase chromosomes were nearly consistent with those in the mitotic metaphase chromosomes when the two split chromosomes were integrated as one unit. Of note, 14 chromosomes (there is a high chance that sex chromosomes are included) were exclusively identified by (AG3T3)3, 5S rDNA, and (TTG)6. The other 10 also showed a terminal signal with (AG3T3)3. Moreover, these oligo-probes were able to distinguish one wild H. rhamnoides taxon from four H. rhamnoides taxa. These chromosome identification and taxa differentiation data will help in elucidating visual and elaborate physical mapping and guide breeders’ utilization of wild resources of H. rhamnoides.

Explicating genetic diversity based on ITS characterization and determination of antioxidant potential in sea buckthorn (Hippophae spp.)

BACKGROUND

Sea buckthorn (Hippophae) is in the focus of interest mainly for its positive effects on health of both human and animal organisms. Due to the similarities in vegetative morphology, Hippophae species are often misidentified. Therefore, current study was focused on ITS based sequence characterization of sea buckthorn species and comparative biochemical evaluation for its antioxidant properties.

METHODS AND RESULTS

DNA was extracted from leaf samples. Primer pairs K-Lab-SeaBukRhm-ITS1F1- K-Lab-SeaBukRhm-ITS1R1 and K-LabSeaBukTib- ITSF1- K-LabSeaBukTib-ITSR1 were used for PCR amplification. The purified PCR products were outsourced for sequencing. Phylogenetic tree was constructed based on neighbor-joining (NJ) method. Moreover, comparison of antioxidant potential of leaves of two sea buckthorn species (Hippophae rhamnoides and Hippophae tibetana) collected from different regions of Ladakh viz., Stakna, Nubra, DRDO Leh and Zanskar was determined by 2,2-diphenyl-1-picrylhydrazyl (DPPH), 2,2-azino-bis (3- ethylbenzothiazoline-6-sulphonic acid) diammonium salt (ABTS), and Total antioxidant capacity (TAC) by phosphomolybdenum assays. The present investigation led to the differentiation of two sea buckthorn species viz., H. rhamnoides and H. tibetana based on Internal Transcribed Spacer (ITS) region. Moreover, significant variation was observed in antioxidant potential of leaf extracts collected from different regions.

CONCLUSIONS

Primary ITS sequence analysis was found to be powerful tool for identification and genetic diversity studies in sea buckthorn. Leaves of sea buckthorn have pronounced antioxidant properties and can be used in food, neutraceuticals and pharmaceutical industries etc. The current study will pave the way to discover small bioactive molecules responsible for antioxidant and anticancer properties in sea buckthorn.

An ABA-flavonoid relationship contributes to the differences in drought resistance between different sea buckthorn subspecies

Drought is the most severe abiotic stress and hinders the normal growth and development of plants. Sea buckthorn (Hippophae rhamnoides Linn.) is a typical drought-resistant tree species. In this study, the leaves of the H. rhamnoides ssp. sinensis ('FN') and H. rhamnoides ssp. mongolica ('XY') were selected during drought-recovery cycles for RNA sequencing, and physiological and biochemical analyses. The results revealed that drought stress significantly decreased leaf water potential, net photosynthetic rate and stomatal conductance in both sea buckthorn subspecies. Similarly, the contents of flavone, flavonol, isoflavone and flavanone significantly decreased under drought stress in 'XY'. Conversely, in 'FN', the flavone and abscisic acid (ABA) contents were significantly higher under drought stress and recovered after rehydration. Meanwhile, 4618 and 6100 differentially expressed genes (DEGs) were identified under drought stress in 'FN' and 'XY', respectively. In total, 5164 DEGs were observed in the comparison between 'FN' and 'XY' under drought stress. This was more than the 3821 and 3387 DEGs found when comparing the subspecies under control and rehydration conditions, respectively. These DEGs were mainly associated with carotenoid biosynthesis, flavonoid biosynthesis, photosynthesis and plant hormone signal transduction. Six hub DEGs (ABCG5, ABCG22, ABCG32, ABCG36, ABF2 and PYL4) were identified to respond to drought stress based on weighted gene co-expression network analysis and Basic Local Alignment Search Tool (BLAST) analysis using DroughtDB. These six DEGs were annotated to play roles in the ABA-dependent signaling pathway. Sixteen RNA sequencing results involving eight genes and similar expression patterns (12/16) were validated using quantitative real-time Polymerase Chain Reaction (PCR). The biochemical and molecular mechanisms underlying the regulation of drought responses by ABA and flavonoids in sea buckthorn were clarified. In this study, gene co-expression networks were constructed, and the results suggested that the mutual regulation of ABA and flavonoid signaling contributed to the difference in drought resistance between the different sea buckthorn subspecies.

Diversity in sea buckthorn (Hippophae rhamnoides L.) accessions with different origins based on morphological characteristics, oil traits, and microsatellite markers

Sea buckthorn (Hippophae rhamnoides) is an ecologically and economically important species. Here, we assessed the diversity of 78 accessions cultivated in northern China using 8 agronomic characteristics, oil traits (including oil content and fatty acid composition) in seeds and fruit pulp, and SSR markers at 23 loci. The 78 accessions included 52 from ssp. mongolica, 6 from ssp. sinensis, and 20 hybrids. To assess the phenotypic diversity of these accessions, 8 agronomic fruit traits were recorded and analyzed using principal component analysis (PCA). The first two PCs accounted for approximately 78% of the variation among accessions. The oil contents were higher in pulp (3.46–38.56%) than in seeds (3.88–8.82%), especially in ssp. mongolica accessions. The polyunsaturated fatty acid (PUFA) ratio was slightly lower in the seed oil of hybrids (76.06%) than that of in ssp. mongolica (77.66%) and higher than that of in ssp. sinensis (72.22%). The monounsaturated fatty acid (MUFA) ratio in the pulp oil of ssp. sinensis (57.00%) was highest, and that in ssp. mongolica (51.00%) was equal to the ratio in the hybrids (51.20%). Using canonical correspondence analysis (CCA), we examined the correlation between agronomic traits and oil characteristics in pulp and seeds. Oil traits in pulp from different origins were correlated with morphological groupings (r = 0.8725, p = 0.0000). To assess the genotypic diversity, 23 SSR markers (including 17 loci previously reported) were used among the 78 accessions with 59 polymorphic amplified fragments obtained and an average PIC value of 0.2845. All accessions were classified into two groups based on the UPGMA method. The accessions of ssp. sinensis and ssp. mongolica were genetically distant. The hybrid accessions were close to ssp. mongolica accessions. The 8 agronomic traits, oil characteristics in seed and pulp oils, and 23 SSR markers successfully distinguished the 78 accessions. These results will be valuable for cultivar identification and genetic diversity analysis in cultivated sea buckthorn.

Characterization of the Hippophae rhamnoides WRKY gene family and functional analysis of the role of the HrWRKY21 gene in resistance to abiotic stresses

Sea buckthorn (Hippophae rhamnoides L.) is a plant with economic and ecological value. It is uniquely capable of growing well under salt and drought stress. WRKY transcription factors play important roles in the ability of plants to resist stress. In this study, 48 HrWRKY genes were identified based on RNA sequencing of H. rhamnoides. Evaluation of expression pattern of HrWRKY1, HrWRKY17, HrWRKY18, HrWRKY21, HrWRKY33-2, HrWRKY40-2, HrWRKY41, and HrWRKY71 suggested that they were involved in abiotic stress. Interestingly, HrWRKY21, one of eight HrWRKY genes, was a positive regulator of abiotic stress tolerance in H. rhamnoides. In addition, most morphological attributes of roots in transgenic Nicotiana tabacum lines (overexpressing HrWRKY21) were also markedly increased compared with the wild-type (WT), including total lengths, specific root lengths, and surface areas. Stress tolerance of transgenic lines was also correlated with higher antioxidant activity (SOD and POD), lower percentage of relative conductivity (REC), and lower activity of malondialdehyde (MDA) under stress conditions. These findings represent a foundation of knowledge about the molecular mechanisms driving resistance to adverse conditions in plants; they are a promising step towards development of tree cultivars with improved tolerance to abiotic stress.

RNA-seq data reveals a coordinated regulation mechanism of multigenes involved in the high accumulation of palmitoleic acid and oil in sea buckthorn berry pulp

BACKGROUND

Sea buckthorn is a woody oil crop in which palmitoleic acid (C16:1n7, an omega-7 fatty acid (FA)) contributes approximately 40% of the total FA content in berry pulp (non-seed tissue). However, the molecular mechanisms contributing to the high accumulation of C16:1n7 in developing sea buckthorn berry pulp (SBP) remain poorly understood.

RESULTS

We identified 1737 unigenes associated with lipid metabolism through RNA-sequencing analysis of the four developmental stages of berry pulp in two sea buckthorn lines, 'Za56' and 'TF2-36'; 139 differentially expressed genes were detected between the different berry pulp developmental stages in the two lines. Analyses of the FA composition showed that the C16:1n7 contents were significantly higher in line 'Za56' than in line 'TF2-36' in the mid-late developmental stages of SBP. Additionally, qRT-PCR analyses of 15 genes involved in FA and triacylglycerol (TAG) biosynthesis in both lines revealed that delta9-ACP-desaturase (ACP-Δ9D) competed with 3-ketoacyl-ACP-synthase II (KASII) for the substrate C16:0-ACP and that ACP-Δ9D and delta9-CoA-desaturase (CoA-Δ9D) gene expression positively correlated with C16:1n7 content; KASII and fatty acid elongation 1 (FAE1) gene expression positively correlated with C18:0 content in developing SBP. Specifically, the abundance of ACP-Δ9D and CoA-Δ9D transcripts in line 'Za56', which had a higher C16:1n7 content than line 'TF2-36', suggests that these two genes play an important role in C16:1n7 biosynthesis. Furthermore, the high expressions of the glycerol-3-phosphate dehydrogenase (GPD1) gene and the WRINKLED1 (WRI1) transcription factor contributed to increased biosynthesis of TAG precursor and FAs, respectively, in the early developmental stages of SBP, and the high expression of the diacylglycerol O-acyltransferase 1 (DGAT1) gene increased TAG assembly in the later developmental stages of SBP. Overall, we concluded that increased ACP-Δ9D and CoA-Δ9D levels coupled with decreased KASII and FAE1 activity is a critical event for high C16:1n7 accumulation and that the coordinated high expression of WRI1, GPD1, and DGAT1 genes resulted in high oil accumulation in SBP.

CONCLUSION

Our results provide a scientific basis for understanding the mechanism of high C16:1n7 accumulation in berry pulp (non-seed tissue) and are valuable to the genetic breeding programme for achieving a high quality and yield of SBP oil.

Transcriptomic analysis of drought stress responses of sea buckthorn (Hippophae rhamnoidessubsp. sinensis) by RNA-Seq

Sea buckthorn is one of the most important eco-economic tree species in China due to its ability to grow and produce acceptable yields under limited water and fertilizer availability. In this study, the differentially expressed genes under drought stress (DS) of sea buckthorn were identified and compared with control (CK) by RNA-Seq. A total of 122,803 unigenes were identified in sea buckthorn, and 70,025 unigenes significantly matched a sequence in at least one of the seven databases. A total of 24,060 (19.59%) unigenes can be assigned to 19 KEGG pathways, and 1,644 unigenes were differentially expressed between DS and CK, of which 519 unigenes were up-regulated and 1,125 unigenes down-regulated. Of the 47 significantly enriched GO terms, 14, 7 and 26 items were related to BP, CC and MF, respectively. KEGG enrichment analysis showed 398 DEGs involved in 97 different pathways, of which 119 DEGs were up-regulated and 279 DEGs were down-regulated under drought stress. In addition, we found 4438 transcriptor factors (TFs) in sea buckthorn, of which 100 were differentially expressed between DS and CK. These results lay a first foundation for further investigations of the very specific functions of these unigenes in sea buckthorn in response to drought stress.

Molecular Sex Identification in Dioecious Hippophae rhamnoides L. via RAPD and SCAR Markers

The dioecious property of the sea buckthorn (Hippophae rhamnoides L.) prevents sex recognition via traditional observation at the juvenile stage, thus impeding breeding and economic cropping; A random amplified polymorphic DNA (RAPD) and a sequence characterized amplified region (SCAR) markers were used to identify the sexes. A total of 45 random decamer primers were used to screen genomic DNA pools of staminate and pistillate genotypes for genetic polymorphisms. One female sex-linked marker was identified. D15 (5′-CATCCGTGCT-3′) amplified a particular band of 885 bp, which showed polymorphism among staminate and pistillate genotype plants. The SCAR marker Hrcx-15 was obtained by sequencing the fragment. The alleles of 140 pistillate genotypes were examined but not of the 140 staminate genotypes discerned via taxonomy. Staminate and pistillate genotypes of sea buckthorn plants can be distinguished, using Hrcx-15 as a genetic marker for sex identification and for expediting cultivation for commercial applications.

Testing the effect of the Himalayan mountains as a physical barrier to gene flow in Hippophae tibetana Schlect. (Elaeagnaceae)

Hippophae tibetana is a small, dioecious wind-pollinated shrub endemic to the Tibetan-Qinghai Plateau. It is one of the shrubs that occur at very high elevations (5250 m a.s.l.). The Himalayan mountains provides a significant geographical barrier to the Qinghai-Tibetan Plateau, dividing the Himalayan area into two regions with Nepal to the south and Tibet to the north. There is no information on how the Himalayan mountains influence gene flow and population differentiation of alpine plants. In this study, we analyzed eight nuclear microsatellite markers and cpDNA trnT-trnF regions to test the role of the Himalayan mountains as a barrier to gene flow between populations of H. tibetana. We also examined the fine-scale genetic structure within a population of H. tibetana on the north slope of Mount (Mt.) Everest. For microsatellite analyses, a total of 241 individuals were sampled from seven populations in our study area (4 from Nepal, 3 from Tibet), including 121 individuals that were spatially mapped within a 100 m × 100 m plot. To test for seed flow, the cpDNA trnT-trnF regions of 100 individuals from 6 populations (4 from Nepal, 2 from Tibet) were also sequenced. Significant genetic differentiation was detected between the two regions by both microsatellite and cpDNA data analyses. These two datasets agree about southern and northern population differentiation, indicating that the Himalayan mountains represent a barrier to H. tibetana limiting gene flow between these two areas. At a fine scale, spatial autocorrelation analysis suggests significant genetic structure within a distance of less than 45 m, which may be attributed mainly to vegetative reproduction and habitat fragmentation, as well as limited gene flow.

ISSR markers for gender identification and genetic diagnosis of Hippophae rhamnoides ssp. turkestanica growing at high altitudes in Ladakh region (Jammu and Kashmir)

Hippophae rhamnoides L. ssp. turkestanica (Elaeagnaceae) is a predominantly dioecious and wind-pollinated medicinal plant species. The mature fruits of the species possess antioxidative, anti-inflammatory, antimicrobial, anticancerous, and antistimulatory properties that are believed to improve the immune system. The identification of male and female plants in H. rhamnoides ssp. turkestanica is quite difficult until flowering which usually takes 3-4 years or more. A sex-linked marker can be helpful in establishing the orchards through identification of genders at an early stage of development. Therefore, we studied the genetic diversity of populations in Ladakh with the aim to identify a gender-specific marker using ISSR markers. Fifty-eight ISSR primers were used to characterize the genome of H. rhamnoides ssp. turkestanica, of which eight primers generated 12 sex-specific fragments specific to one or more populations. The ISSR primer (P-45) produced a fragment which faithfully segregates all the males from the female plants across all the three valleys surveyed. This male-specific locus was converted into a SCAR. Forward and reverse primers designed from this fragment amplified a 750-bp sequence in males only, thus specifying it as an informative male-specific sex-linked marker. This SCAR marker was further validated for its capability to differentiate gender on an additional collection of plants, representing three geographically isolated valleys (Nubra, Suru, and Indus) from Ladakh region of India. The results confirmed sex-linked specificity of the marker suggesting that this conserved sequence at the Y chromosome is well preserved through the populations in Ladakh region. At present, there are no reliable markers which can differentiate male from female plants across all the three valleys of Ladakh region at an early stage of plant development. It is therefore envisaged that the developed SCAR marker shall provide a reliable molecular tool for early identification of the sex in this commercial crop. The genetic diversity of populations as surveyed by ISSR primers revealed 85.71 % polymorphism at the population level. The dendrogram obtained divided the genotypes into three different clusters, and the distribution of male and female genotypes in all the clusters was random. The Nei's genetic similarity index was in the range of 0.63-0.96.

Satellite DNA and Transposable Elements in Seabuckthorn (Hippophae rhamnoides), a Dioecious Plant with Small Y and Large X Chromosomes

Seabuckthorn (Hippophae rhamnoides) is a dioecious shrub commonly used in the pharmaceutical, cosmetic, and environmental industry as a source of oil, minerals and vitamins. In this study, we analyzed the transposable elements and satellites in its genome. We carried out Illumina DNA sequencing and reconstructed the main repetitive DNA sequences. For data analysis, we developed a new bioinformatics approach for advanced satellite DNA analysis and showed that about 25% of the genome consists of satellite DNA and about 24% is formed of transposable elements, dominated by Ty3/Gypsy and Ty1/Copia LTR retrotransposons. FISH mapping revealed X chromosome-accumulated, Y chromosome-specific or both sex chromosomes-accumulated satellites but most satellites were found on autosomes. Transposable elements were located mostly in the subtelomeres of all chromosomes. The 5S rDNA and 45S rDNA were localized on one autosomal locus each. Although we demonstrated the small size of the Y chromosome of the seabuckthorn and accumulated satellite DNA there, we were unable to estimate the age and extent of the Y chromosome degeneration. Analysis of dioecious relatives such as Shepherdia would shed more light on the evolution of these sex chromosomes.

Proanthocyanidins in Sea Buckthorn (Hippophaë rhamnoides L.) Berries of Different Origins with Special Reference to the Influence of Genetic Background and Growth Location

Wild sea buckthorn berries from Finland (Hippophaë rhamnoides ssp. rhamnoides) and China (ssp. sinensis) as well as berries of two varieties of ssp. rhamnoides cultivated in Finland and five of ssp. mongolica cultivated in Canada were compared on the basis of the content and composition of proanthocyanidins (PAs). Among all of the samples, only B-type PAs were found. The contents of dimeric, trimeric, tetrameric, and total PAs were in the range of 1.4-8.9, 1.3-9.5, 1.0-7.1, and 390-1940 mg/100 g of dry weight, respectively. The three subspecies were separated by three validated factors (R(2), 0.724; Q(2), 0.677) in the partial least squares discriminant analysis model. Significant differences in total PAs were found between the ssp. rhamnoides and mongolica samples (p < 0.05). In ssp. rhamnoides, samples grown in northern Finland were characterized by a high amount of total PAs, typically 2-3 times higher than that in the level found in southern Finland. In ssp. sinensis, altitude did not have a systematic effect on the PA composition, suggesting the significance of the interaction between genetic background and growth location.

[Genetic and chemical discrimination of traditional Tibetan medicine seabuckthorn based on DNA barcode and ¹H-NMR metabolic method]

To differentiate three medicinal Hippopahe species of seabuckthorn, a combined genetic and chemical identification method was established in this study. ITS2 and psbA-trnH were tested for identification of 3 species of seabuckthorn. Detection of the kimura 2-parameter (K2P) distance, the neighbor-joining (NJ) tree and the barcoding gap were used to assess the identification efficiency. ¹H-NMR based metabolic method was applied to acquire the profile of metabolites. PCA was used to analysis the metabolite data. The results indicated that DNA barcode combined ¹H-NMR based metabolic method is a powerful tool for the identification of 3 medicinal Hippopahe species of seabuckthorn. The finding demonstrated that different genetic variation and chemical constituents existed among 3 medicinal Hippopahe species of seabuckthorn. The combined identification method will improve the reliability of species discrimination and could be applicable to much other ethnic medicine which has various origins in China.

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

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.

DeepSAGE based differential gene expression analysis under cold and freeze stress in seabuckthorn (Hippophae rhamnoides L.)

Seabuckthorn (Hippophae rhamnoides L.), an important plant species of Indian Himalayas, is well known for its immense medicinal and nutritional value. The plant has the ability to sustain growth in harsh environments of extreme temperatures, drought and salinity. We employed DeepSAGE, a tag based approach, to identify differentially expressed genes under cold and freeze stress in seabuckthorn. In total 36.2 million raw tags including 13.9 million distinct tags were generated using Illumina sequencing platform for three leaf tissue libraries including control (CON), cold stress (CS) and freeze stress (FS). After discarding low quality tags, 35.5 million clean tags including 7 million distinct clean tags were obtained. In all, 11922 differentially expressed genes (DEGs) including 6539 up regulated and 5383 down regulated genes were identified in three comparative setups i.e. CON vs CS, CON vs FS and CS vs FS. Gene ontology and KEGG pathway analysis were performed to assign gene ontology term to DEGs and ascertain their biological functions. DEGs were mapped back to our existing seabuckthorn transcriptome assembly comprising of 88,297 putative unigenes leading to the identification of 428 cold and freeze stress responsive genes. Expression of randomly selected 22 DEGs was validated using qRT-PCR that further supported our DeepSAGE results. The present study provided a comprehensive view of global gene expression profile of seabuckthorn under cold and freeze stresses. The DeepSAGE data could also serve as a valuable resource for further functional genomics studies aiming selection of candidate genes for development of abiotic stress tolerant transgenic plants.

Microrefugia and Shifts of Hippophae tibetana (Elaeagnaceae) on the north side of Mt. Qomolangma (Mt. Everest) during the last 25000 years

Microrefugia at high altitudes or high latitudes are thought to play an important role in the post-glacial colonization of species. However, how populations in such microrefugia have responded to climate changes in alternating cold glacial and warm interglacial stages remain unclear. Here we present evidence to indicate the Rongbuk Valley of the Mt. Qomolangma (Mt. Everest) area, the highest region on earth, had microrefugia for Hippophae tibetana and discuss how this low shrub was adapted to the extreme climate fluctuations of the last 25,000 years by shifts. By integrating geological, glaciological, meteorological, and genetic information, we found that the Rongbuk Valley was not only a glacial microrefugium but also an interglacial microrefugium for H. tibetana: the former was located on the riverbank below 4800 m above sea level (asl) or lower area and the latter at ∼5000 m asl. Our results show that after the Last Glacial Maximum (LGM), H. tibetana in the valley has undergone upward and downward migrations around ∼5000 m driven by climate fluctuations and the population in the glacial microrefugium has suffered extinction or extreme contraction. Moreover, with the rise of temperature in the last four decades, the upper limit of H. tibetana has shifted at least 30 m upward. Combining population history and recent range shift of this species is important in predicting the fate of this endemic species to future climate changes.

Optimization of de novo short read assembly of seabuckthorn (Hippophae rhamnoides L.) transcriptome

Seabuckthorn (Hippophaerhamnoides L.) is known for its medicinal, nutritional and environmental importance since ancient times. However, very limited efforts have been made to characterize the genome and transcriptome of this wonder plant. Here, we report the use of next generation massive parallel sequencing technology (Illumina platform) and de novo assembly to gain a comprehensive view of the seabuckthorn transcriptome. We assembled 86,253,874 high quality short reads using six assembly tools. At our hand, assembly of non-redundant short reads following a two-step procedure was found to be the best considering various assembly quality parameters. Initially, ABySS tool was used following an additive k-mer approach. The assembled transcripts were subsequently subjected to TGICL suite. Finally, de novo short read assembly yielded 88,297 transcripts (> 100 bp), representing about 53 Mb of seabuckthorn transcriptome. The average length of transcripts was 610 bp, N50 length 1198 BP and 91% of the short reads uniquely mapped back to seabuckthorn transcriptome. A total of 41,340 (46.8%) transcripts showed significant similarity with sequences present in nr protein databases of NCBI (E-value < 1E-06). We also screened the assembled transcripts for the presence of transcription factors and simple sequence repeats. Our strategy involving the use of short read assembler (ABySS) followed by TGICL will be useful for the researchers working with a non-model organism’s transcriptome in terms of saving time and reducing complexity in data management. The seabuckthorn transcriptome data generated here provide a valuable resource for gene discovery and development of functional molecular markers.

Effects of genotype, latitude, and weather conditions on the composition of sugars, sugar alcohols, fruit acids, and ascorbic acid in sea buckthorn (Hippophaë rhamnoides ssp. mongolica) berry juice

Sea buckthorn berries (Hippophaë rhamnoides ssp. mongolica) of nine varieties were collected from three growth locations in five inconsecutive years (n = 152) to study the compositional differences of sugars, sugar alcohols, fruit acids, and ascorbic acid in berries of different genotypes. Fructose and glucose (major sugars) were highest in Chuiskaya and Vitaminaya among the varieties studied, respectively. Malic acid and quinic acid (major acids) were highest in Pertsik and Vitaminaya, respectively. Ascorbic acid was highest in Oranzhevaya and lowest in Vitaminaya. Berry samples of nine varieties collected from two growth locations in five years (n = 124) were combined to study the effects of latitude and weather conditions on the composition of H. rhamnoides ssp. mongolica. Sea buckthorn berries grown at lower latitude had higher levels of total sugar and sugar/acid ratio and a lower level of total acid and were supposed to have better sensory properties than those grown at higher latitude. Glucose, quinic acid, and ascorbic acid were hardly influenced by weather conditions. The other components showed various correlations with temperature, radiation, precipitation, and humidity variables. In addition, fructose, sucrose, and myo-inositol correlated positively with each other and showed negative correlation with malic acid on the basis of all the samples studied (n = 152).

Effect of chitinase on resistance to fungal pathogens in sea buckthorn, Hippophae rhamnoides, and cloning of Class I and III chitinase genes

Sea buckthorn (Hippophae rhamnoides L.) is naturally distributed from Asia to Europe. It has been widely planted as an ornamental shrub and is rich in nutritional and medicinal compounds. Fungal pathogens that cause diseases such as dried-shrink disease are threats to the production of this plant. In this study, we isolated the dried-shrink disease pathogen from bark and total chitinase protein from leaves of infected plants. The results of the Oxford Cup experiment suggested that chitinase protein inhibited the growth of this pathogen. To improve pathogen resistance, we cloned chitinase Class I and III genes in H. rhamnoides, designated Hrchi1 and Hrchi3. The full-length cDNA of the open reading frame region of Hrchi1 contained 903 bp encoding 300 amino acids and Hrchi3 contained 894 bp encoding 297 amino acids. Active domain analysis, protein types, and secondary and 3D structures were predicted using online software.

Expressed sequence tag based identification and expression analysis of some cold inducible elements in seabuckthorn (Hippophae rhamnoides L.)

A cDNA library was constructed from the mature leaves of seabuckthorn (Hippophae rhamnoides). Expressed Sequence Tags (ESTs) were generated by single pass sequencing of 4500 cDNA clones. We submitted 3412 ESTs to dbEST of NCBI. Clustering of these ESTs yielded 1665 unigenes comprising of 345 contigs and 1320 singletons. Out of 1665 unigenes, 1278 unigenes were annotated by similarity search while the remaining 387 unannotated unigenes were considered as organism specific. Gene Ontology (GO) analysis of the unigene dataset showed 691 unigenes related to biological processes, 727 to molecular functions and 588 to cellular component category. On the basis of similarity search and GO annotation, 43 unigenes were found responsive to biotic and abiotic stresses. To validate this observation, 13 genes that are known to be associated with cold stress tolerance from previous studies in Arabidopsis and 3 novel transcripts were examined by Real time RT-PCR to understand the change in expression pattern under cold/freeze stress. In silico study of occurrence of microsatellites in these ESTs revealed the presence of 62 Simple Sequence Repeats (SSRs), some of which are being explored to assess genetic diversity among seabuckthorn collections. This is the first report of generation of transcriptome data providing information about genes involved in managing plant abiotic stress in seabuckthorn, a plant known for its enormous medicinal and ecological value.

Isolation of good quality RNA from a medicinal plant seabuckthorn, rich in secondary metabolites

Medicinal plants are being widely investigated owing to their ability to produce molecules of therapeutic significance. Isolation of good quality RNA is a tedious but primary step towards undertaking molecular biology experiments. However, medicinal plants are rich in secondary metabolites and not amenable to standard RNA isolation protocols involving Guanidine isothiocyanate (GITC). So an RNA isolation protocol from difficult samples (richer in secondary metabolites) is of highest desiderata. Here we propose a new protocol suitable for isolating RNA from plant tissues rich in secondary metabolites. To standard CTAB (Cetyl Trimethyl Ammonium Bromide) buffer, addition of 2% PVPP (polyvinyl polypyrrolidone) and 350 mM beta-mercaptoethanol was found useful. Use of glacial acetic acid (1M) along with ethanol for precipitation after phenolization and chloroform extraction enhanced the RNA yield. This is the first report of using glacial acetic acid in a CTAB based protocol for the precipitation of RNA. This protocol has been validated in medicinal plant Hippophae rhamnoides vern. seabuckthorn, where standard RNA isolation methods involving GITC and TRIZol extraction buffers failed. The RNA isolated by this method was of good quality as gauged by spectrophotometric readings and denaturing agarose gel electrophoresis. To the best of our knowledge, this RNA isolation protocol has never been published before. The RNA thus obtained could be suitably used for the downstream molecular procedures like Reverse Transcription Polymerase Chain Reaction (RT-PCR), Real Time-PCR, cDNA library construction, etc.

[Genetic diversity of Hippophae rhamnoides L. sub-populations on loess plateau under effects of varying meteorologic conditions]

The genetic diversity and genetic differentiation of eight Hippophae rhamnoides L. populations in eastern Gansu Province, as well as the effects of varying meteorologic conditions on them were assessed by using inter-simple sequence repeat (ISSR) molecular marker method. A total of 240 individuals in the populations were sampled. Using eleven primers, 165 bands were generated, ranging in size from 300 to 1500 bp, and 157 (95.76%) were found to be polymorphic. The analysis of molecular variance (AMOVA) demonstrated that there was a relatively high level (76.5%) of genetic variation within the populations, with the gene differentiation coefficient (Gst) and gene flow being 0.2418 and 1.5675, respectively. Therefore, to protect the gene resources of H. rhamnoides, the individuals within the populations should be first considered. Mantel test showed that genetic distance was significantly positively correlated with geographical distance (r = 0.65, p = 0.002), and regression modeling between genetic diversity and meteorologic factors suggested that there was a significant positive correlation between wind speed during blooming and genetic diversity of H. rhamnoides, illustrating that wind speed in blooming period and geographic distance were the vital factors affecting the genetic diversity of H. rhamnoides population.

High-performance liquid chromatographic fingerprint analysis for different origins of sea buckthorn berries

Using high-performance liquid chromatography (HPLC), a chemical fingerprint method was developed for investigating and demonstrating the variance of flavonoids among different origins of sea buckthorn berries. Thirty-four samples were analyzed including 15 RS (Hippophae rhamnoides ssp. sinensis) samples, 7 RY (H. rhamnoindes ssp. yunnanensis) samples, 5 RW (H. rhamnoides ssp. wolongensis) samples, 4 NS (H. neurocarpa ssp. stellatopilosa) samples and 3 TI (H. tibetana) samples. In the HPLC chromatograms, 12 compounds were identified as flavonoids, including quercetin 3-O-sophoroside-7-rhamnoside, kaempferol 3-O-sophoroside-7-O-rhamnoside, isorhamnetin 3-O-sophoroside-7-O-rhamnoside, isorhamnetin 3-O-glucoside-7-O-rhamnoside, quercetin 3-O-rutinoside, quercetin 3-O-glucoside, isorhamnetin 3-O-rutinoside, isorhamnetin 3-O-glucoside, quercetin, kaempferol 7-O-rhamnoside, kaempferol and isorhamnetin. Both correlation coefficient of similarity in chromatograms and relative peak areas of characteristic compounds were calculated for quantitative expression of the HPLC fingerprints. Our results revealed that the chromatographic fingerprint combining similarity evaluation could efficiently identify and distinguish sea buckthorn berries from different species. However, no obvious difference between RS and RY suggested that the two subspecies might have very close relationship in terms of chemotaxonomy. The established method was considered to be suitable for fingerprint analysis to check the genuine origin and control the quality of sea buckthorn berries and extracts.

Late Quaternary history of Hippophaë rhamnoides L. (Elaeagnaceae) inferred from chalcone synthase intron (Chsi) sequences and chloroplast DNA variation

Fossil pollen records indicate that Hippophaë rhamnoides (Sea Buckthorn) was widespread on late- and early postglacial raw soils throughout much of central and northern Europe, but that Early Holocene reforestation restricted populations to northern coastal habitats, or along mountain streams in the Alps, Pyrenees, and Carpathians. We used sequence variation at the nuclear chalcone synthase intron (Chsi), in conjunction with chloroplast DNA-restriction fragment length polymorphism data, to investigate the intraspecific phylogeny, phylogeographic structure, and expansion demographic history of this dioecious and wind-pollinated shrub at its range-wide scale in Europe and Asia Minor. Four major Chsi phylogroups of unresolved relationships were identified with estimated divergences approximately 172,000 years ago. Large-scale phylogeographic structures of nuclear and cytoplasmic markers were congruent in identifying (i) southeastern Europe as the most likely source of colonization into central Europe and Scandinavia, and (ii) the area just north of the Alps as a contact zone between populations from the Alps and the east/central European-Scandinavian lineage. Coalescence-based analyses (i.e. nested clade analysis and mismatch distributions) of Chsi variation were able to detect at least four major episodes of population growth, all within about the last 40,000 years. In particular, these analyses identified a nearly synchronized timing of population expansions in various parts of the species' range in central-eastern Europe/Asia Minor, most likely correlating with the Younger Dryas Stadial ( approximately 13,000-11,600 years ago). It remains to be established whether the phylogeographic history of H. rhamnoides, and particularly its rapid response to the rapid environmental changes of the Younger Dryas cold snap, is unique to the species, or whether it is shared with other cold-tolerant shrub (or grassland) species known from late-glacial raw soils in Europe.

Genetic Variation in Hippophae rhamnoides ssp. sinensis (Elaeagnaceae) Revealed by RAPD Markers

Hippophae rhamnoides ssp. sinensis is endemic to China, and it is a dioecious, outcrossing plant. Although many studies have been undertaken mainly on its agricultural, nutritional, medical, and ornamental value, little is known about its population genetics. This study uses random amplified polymorphic DNA to investigate the genetic diversity and population genetic structure of 13 natural populations of the subspecies sinensis. Fifteen primers amplified 107 reproducible bands, with 95 (88.79%) being polymorphic. The gene diversity within population was 0.168, considerably lower than that of tree species and most perennial, outcrossing species, but higher than that of annual or short-lived, selfing species. The Gst value showed that 18.3% of the total genetic variation resided among populations, a little lower than that of outcrossing species. The present results are quite similar to those previously reported in another subspecies, H ssp. . rhamnoides rhamnoides. The low genetic differentiation among populations in ssp. sinensis may be attributed to the long-distance dispersal of seeds facilitated by birds, in addition to its characteristics of outcrossing, wind pollination, and widespread distribution. No association between genetic distance and geographical distribution was found. The population relationships revealed by the UPGMA dendrogram parallel this result, in that genetic distance did not increase with geographic separation. This pattern of population differentiation may imply the adaptation of ssp. s populations to the local environment, given that its habitats vary greatly across its distribution.

Molecular genetic variation in Chinese populations of three subspecies of Hippophae rhamnoides

To provide a population-level genetic profile for investigation and conservation of genetic diversity of the sea buckthorn (Hippophae rhamnoides L.), 300 individuals of fifteen natural populations of sea buckthorn in China were analyzed by using ISSR (inter-simple sequence repeats) markers. On the basis of Shannon's index and Nei's genetic diversity, the mean genetic diversity detected in the natural populations of Hippophae rhamnoides ssp. yunnanensis, ssp. sinensis, and ssp. gyantsensis was 0.1944, 0.2169, and 0.1372, respectively. The coefficient of gene differentiation (Gst) of seven ssp. yunnanensis populations is 0.2790, and that of seven ssp. sinensis populations is 0.4184. This means that 28% of the total molecular variance of seven ssp. yunnanensis populations existed among populations, and 42% for seven ssp. sinensis populations, suggesting that the subspecies have different genetic structures. No significant correlation between genetic and geographic distances of the populations was found using ISSR markers.

Sex ratio of some long-lived dioecious plants in a sand dune area

In dioecious plants the fraction of males among flowering plants in the field (the secondary sex ratio) is the result of the fraction of males in the seeds (the primary sex ratio) and the subsequent survival and age at first reproduction of the two genders. It has been assumed that survival and age at first reproduction are the main determinants of biased secondary sex ratio but, especially for long-lived perennials, few data are available. We address this issue for natural populations of four long-lived perennials in a dune area. In Asparagus officinale and Bryonia dioica, the secondary sex ratio was unbiased. In Salix repens the secondary sex ratio was female-biased (0.337). Hippophae rhamnoides populations were male-biased; the average sex ratio of flowering plants was 0.658, while the fraction of males varied between 0.39 near the sea to 0.84 at the inland side of the dunes. The primary sex ratio was estimated by germinating seeds and growing plants under favourable conditions with minimal mortality. In S. repens the primary sex ratio in seeds was variable among mother plants and was, on average, female-biased (0.289). This is close to the secondary sex ratio, suggesting that the female bias already originates in the seed stage. In Hippophae rhamnoides the primary sex ratio was slightly male-biased (0.564). We argue that in this species, apart from the primary sex ratio, higher mortality and a later age at first reproduction for females contribute to the strong male bias among flowering plants in the field.