Chromosome-level genome assembly of Hippophae rhamnoides variety.

Fructus hippophae (Hippophae rhamnoides spp. mongolica×Hippophae rhamnoides sinensis), a hybrid variety of sea buckthorn that Hippophae rhamnoides spp. mongolica serves as the female parent and Hippophae rhamnoides sinensis serves as the male parent, is a traditional plant with great potentials of economic and medical values. Herein, we gained a chromosome-level genome of Fructus hippophae about 918.59 Mb, with the scaffolds N50 reaching 83.65 Mb. Then, we anchored 440 contigs with 97.17% of the total genome sequences onto 12 pseudochromosomes. Next, de-novo, homology and transcriptome assembly strategies were adopted for gene structure prediction. This predicted 36475 protein-coding genes, of which 36226 genes could be functionally annotated. Simultaneously, various strategies were used for quality assessment, both the complete BUSCO value (98.80%) and the mapping rate indicated the high assembly quality. Repetitive elements, which occupied 63.68% of the genome, and 1483600 bp of non-coding RNA were annotated. Here, we provide genomic information on female plants of a popular variety, which can provide data for pan-genomic construction of sea buckthorn and for the resolution of the mechanism of sex differentiation.

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.

Distribution of simple sequence repeats, transcription factors, and differentially expressed genes in the NGS-based transcriptome of male and female seabuckthorn (Hippophae salicifolia).

Seabuckthorn (Hippophae salicifolia) is a perennial, multipurpose wonder plant, popular for its immense medicinal, nutritional, and therapeutic properties. However, due to the lack of whole-genome-based studies, the molecular mechanism governing distinct sexual phenotypes is still not clear. We employed the high-throughput NGS Illumina NovaSeq paired-end technology to generate whole transcriptome profiles of male and female plants of H. salicifolia. In total, 3.2 million raw short reads were generated with an average length of 150 bp, including 50911358 reads from the male leaf tissue samples and 45850364 reads from the female leaf tissue samples. Clustering of the high-quality reads yielded de novo short read assembly of 50259 transcripts of >100 bp length. The final transcripts were assigned Gene Ontology (GO) terms. The digital expression of genes was studied using the DESeq2 of R package that identified 7180 differentially expressed genes (DEGs) between the male and female plant samples. Further, 10,850 simple sequence repeats, and 8,351 transcription factors, distributed in more than 85 transcription families, were also mined from the final assembled transcriptome. Next, COG and KEGG pathway analyses were performed to assign biological functional terms to the DEGs. The findings of the present study will provide a valuable resource for gene expression discovery and other functional genomics studies aiming towards the selection of candidate genes for the development of sex-specific markers in seabuckthorn and other closely related species.Communicated by Ramaswamy H. Sarma.

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.

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.

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.

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.

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-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.

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.

Tandem Mass Tag Based Quantitative Proteomics of Developing Sea Buckthorn Berries Reveals Candidate Proteins Related to Lipid Metabolism.

Sea buckthorn ( Hippophae L.) is an economically important shrub or small tree distributed in Eurasia. Most of its well-recognized medicinal and nutraceutical products are derived from its berry oil, which is rich in monounsaturated omega-7 (C16:1) fatty acid and polyunsaturated omega-6 (C18:2) and omega-3 (C18:3) fatty acids. In this study, tandem mass tags (TMT)-based quantitative analysis was used to investigate protein profiles of lipid metabolism in sea buckthorn berries harvested 30, 50, and 70 days after flowering. In total, 8626 proteins were identified, 6170 of which were quantified. Deep analysis results for the proteins identified and related pathways revealed initial fatty acid accumulation during whole-berry development. The abundance of most key enzymes involved in fatty acid and triacylglycerol (TAG) biosynthesis peaked at 50 days after flowering, but TAG synthesis through the PDAT (phospholipid: diacylglycerol acyltransferase) pathway mostly occurred early in berry development. In addition, the patterns of proteins involved in lipid metabolism were confirmed by combined quantitative real-time polymerase chain reaction, enzyme-linked immunosorbent assay, and parallel reaction monitoring analyses. Our data on the proteomic spectrum of sea buckthorn berries provide a scientific basic for understanding lipid metabolism and related pathways in the developing berries.

Identification and characterization of circular RNAs during the sea buckthorn fruit development.

As a rising star of noncoding RNA, circular RNAs (circRNAs) have a covalently closed loop structure, which formed by 3'-5' ligation during splicing. A few circRNAs were identified and thought to be transcriptional noise due to cognitive defect over the past 40 years. Recently, with the development of high-throughput RNA sequencing techniques and specific algorithms for circRNA detection and quantification, plenty of potential circRNAs were identified in many species which play important roles in various biological processes. However, researches on circRNAs in fruit ripening process were lacking. Here, we totally identified 2616 circRNAs in sea buckthorn fruit development process, which uniformly distributed in sea buckthorn chromosome. Among them, 1721 (65.8%) circRNAs were arising from the exons of their host genes, 252 circRNAs were identified as the differentially expressed circRNAs (DEcircRNAs) between three different development stages, and 181 (71.8%) DEcircRNAs had sequence similarity with 235 identified circRNAs from five know plant species. Functional annotation revealed that host genes of DEcircRNAs were predicted to be involved in carotenoid biosynthesis, lipid synthesis and plant hormone signal transduction. Additionally, 53 DEcircRNAs were predicted as the corresponding nine miRNAs sponges in sea buckthorn. Divergent reverse-transcription PCR and RT-qPCR were used for validate the differential expression and back-splicing sites of six DEcircRNAs. These results revealed the role of circRNAs in sea buckthorn fruit ripening process and promoted the noncoding RNA researches in plants.

Occurrence of subdioecy and scarcity of gender-specific markers reveal an ongoing transition to dioecy in Himalayan seabuckthorn (Hippophae rhamnoides ssp. turkestanica).

Dioecy and the dynamics of its evolution are intensely investigated aspects of plant reproduction. Seabuckthorn (Hippophae rhamnoides ssp. turkestanica) is an alpine shrub growing wild in certain parts of western Himalaya. The previous studies have reported heteromorphic sex chromosomes in the species and yet marker-based studies indicate high similarity between the male and female genomes. Lack of information on sexual system in the species has further complicated the situation. A systematic study was thus undertaken to understand the sexual system in seabuckthorn and to discern the extent of similarity/dissimilarity between the male and female genomes by generating a large number of markers using amplified fragment length polymorphism and representational difference analysis. Floral biology and regular monitoring of species revealed the presence of polygamomonoecious (PGM) plants in most populations at a low percentage (~2-4%). PGM plants showed low pollen production and overall low fertility, suggesting a monoecy-paradioecy pathway at function. The results of the marker study demonstrated that there are limited differences between male and female genomes and these differences were not uniform across the populations in the Leh-Ladakh region, especially when the geographical distance increases. Results also suggest that a dynamic partitioning of genomes is operational between the two genders of seabuckthorn and differences are not homogenized across the populations. Both reproductive biology-based and DNA marker-based studies indicate that genders have separated recently. The present study proposes seabuckthorn as a promising model system to study evolution of dioecy and sex determination.

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.

Transcriptomic and functional analyses unveil the role of long non-coding RNAs in anthocyanin biosynthesis during sea buckthorn fruit ripening.

Fruit ripening is a developmental process regulated by a complex network of endogenous and exogenous cues. Sea buckthorn is an excellent material for fruit ripening studies due to its dramatic ripening process and high contents of nutritional and anti-oxidant compounds in berries. Here, the whole transcriptome of sea buckthorn fruit at three development stages were analysed using multiple high-throughput sequencings. We assembled and annotated 9,008 long non-coding RNAs (lncRNAs) in sea buckthorn fruits, and identified 118 differentially expressed lncRNAs (DE-lncRNAs) and 32 differentially expressed microRNAs in fruit developmental process. In addition, we predicted 1,061 cis-regulated and 782 trans-regulated targets of DE-lncRNAs, and these DE-lncRNAs are specifically enriched in the biosynthesis of ascorbic acid, carotenoids and flavonoids. Moreover, the silencing of two lncRNAs (LNC1 and LNC2) in vivo and expression analysis revealed that LNC1 and LNC2 can act as endogenous target mimics of miR156a and miR828a to reduce SPL9 and induce MYB114 expression, respectively, which lead to increased and decreased anthocyanin content as revealed by high-performance liquid chromatography analysis. Our results present the first global functional analysis of lncRNA in sea buckthorn and provide two essential regulators of anthocyanin biosynthesis, which provides new insights into the regulation of fruit quality.