Comprehensive Proteome and Acetylome Analysis of Needle Senescence in Larix gmelinii.

Leaf senescence is essential for the growth and development of deciduous trees in the next season. Larix gmelinii, a deciduous coniferous tree, exhibits its most distinctive feature by turning yellow in the autumn and eventually shedding its leaves, resulting in significant changes in its appearance during the fall. Lysine acetylation plays an important role in diverse cellular processes; however, limited knowledge is available regarding acetylations in the needle senescence of L. gmelinii. In this study, the proteomics and acetylated modification omics of two phenotypic leaves, yellow and green (senescent and non-senescent) needles, were analyzed before autumn defoliation. In total, 5022 proteins and 4469 unique acetylation sites in 2414 lysine acylated proteins were identified, and this resulted in the discovery of 1335 differentially expressed proteins (DEPs) and 605 differentially expressed acetylated proteins (DAPs) in yellow versus green needles. There are significant differences between the proteome and acetylome; only 269 proteins were found to be DEP and DAP, of which 136 proteins were consistently expressed in both the DEP and DAP, 91 proteins were upregulated, and 45 proteins were down-regulated. The DEPs participate in the metabolism of starch and sucrose, while the DAPs are involved in glycolysis and the tricarboxylic acid cycle. Among them, DEPs underwent significant changes in glycolysis and citric acid cycling. Most of the enzymes involved in glycolysis and the citrate cycle were acetylated. DAPs were down-regulated in glycolysis and up-regulated in the citrate cycle. In all, the results of this study reveal the important role of lysine acetylation in the senescence of L. gmelinii needles and provide a new perspective for understanding the molecular mechanism of leaf senescence and tree seasonal growth.

Effects of bacteria on early-stage litter decomposition in Wudalianchi volcanic forest.

To understand the role of microorganisms in litter decomposition and nutrient cycling in volcanic forest ecosystem, we conducted in-situ litterbag decomposition experiment and used Illumina MiSeq high-throughput sequencing to analyze the response of bacterial community structure and diversity during the decomposition of litters from Larix gmelinii, Betula platyphylla and Populus davidiana, the dominant tree species in volcanic lava plateau of Wudalianchi. The results showed that mass remaining percentage of litters of three species after 18-month decomposition was 63.9%-68.1%. Litter of B. platyphylla decomposed the fastest, with significant difference in N, C:N, and N:P before and after decomposition. The richness of bacterial species and diversity index differed significantly among the three litters. Proteobacteria, Actinomycetes, and Bacteroidetes were the dominant bacterial groups at the phylum level, while Rhizobium, Sphingomonas, and Pseudomonas were the dominant groups at the genus level, with significant difference among the three litters. After 18 months, the dominant bacterial groups in litter tended to be consistent with those in volcanic lava platform soil. In the volcanic forest ecosystem, bacterial diversity and community structure were mainly affected by P, C:N, and N:P in the litter.

Boreal conifers maintain carbon uptake with warming despite failure to track optimal temperatures.

Warming shifts the thermal optimum of net photosynthesis (ToptA) to higher temperatures. However, our knowledge of this shift is mainly derived from seedlings grown in greenhouses under ambient atmospheric carbon dioxide (CO2) conditions. It is unclear whether shifts in ToptA of field-grown trees will keep pace with the temperatures predicted for the 21st century under elevated atmospheric CO2 concentrations. Here, using a whole-ecosystem warming controlled experiment under either ambient or elevated CO2 levels, we show that ToptA of mature boreal conifers increased with warming. However, shifts in ToptA did not keep pace with warming as ToptA only increased by 0.26-0.35 °C per 1 °C of warming. Net photosynthetic rates estimated at the mean growth temperature increased with warming in elevated CO2 spruce, while remaining constant in ambient CO2 spruce and in both ambient CO2 and elevated CO2 tamarack with warming. Although shifts in ToptA of these two species are insufficient to keep pace with warming, these boreal conifers can thermally acclimate photosynthesis to maintain carbon uptake in future air temperatures.

Identification of C2H2 zinc finger genes through genome-wide association study and functional analyses of LkZFPs in response to stresses in Larix kaempferi.

BACKGROUND: C2H2 zinc finger proteins (C2H2-ZFPs), one of the largest transcription factors, play a variety of roles in plant development and growth as well as stress response. While, the evolutionary history and expression profile of the C2H2-ZFP genes in Larix kaempferi (LkZFPs) have not been reported so far. RESULTS: In this study, the whole genome of the LkZFPs was identified and characterized, including physicochemical properties, phylogenetic relationships, conservative motifs, the promoter cis-elements and Gene Ontology (GO) annotation. We identified 47 LkZFPs and divided them into four subfamilies based on phylogenetic analysis and conserved motifs. Subcellular localization prediction showed that most of the LkZFPs were located in the nucleus. Promoter cis-element analysis suggested that the LkZFPs may be involved in the regulation of stress responses. Moreover, Real-time quantitative PCR (RT-qPCR) results showed that Q-type LkZFP genes were involved in the response to abiotic stress, such as salt, drought and hormone stresses. Subcellular localization results showed that LkZFP7 and LkZFP37 were located in the nucleus, LkZFP32 was located in both cytoplasm and nucleus. CONCLUSION: The identification and functional analysis of LkZFPs suggested that some LkZFP genes might play important roles in coping with both biological and abiotic stresses. These results could further increase understanding of the function of the LkZFPs, and provide some research direction and theoretical support.

Characterization of hybrid larch somatic embryo maturation by biochemical analyses and by a novel, fast mid-infrared approach.

The morphology of somatic embryos (SE) is not a sufficient criterion to determine the level of maturation and the optimal stage to transfer embryos for germination, unlike the biochemical components. This composition characterization in the laboratory is too restrictive to be considered at each maturation cycle, as would be necessary. It is, therefore, essential to consider alternative methods. The objectives of this work were to achieve a complete biochemical characterization of the embryos during their development, to serve as a reference and develop a characterization based on infrared spectrometry and chemometrics. During the precotyledonary stage (0-3 weeks of maturation), water content and glucose and fructose levels were high, which is consistent with SE development. After 4 weeks, the cotyledonary SE had a metabolism oriented towards the storage accumulation of lipids, proteins and starch, whereas raffinose only appeared from 8 weeks. Mid-infrared calibration models were developed for water, proteins, lipids, carbohydrates, glucose, fructose, inositols, raffinose, stachyose and starch contents with an r2 average of 0.84. A model was also developed to discriminate the weeks of SE maturation. Different classes of age were discriminated with at least 72% of accuracy. Infrared analysis of the SE based on their full biochemical spectral fingerprint revealed a very slight variation in composition between 7 and 9 weeks, information that is very difficult to obtain by conventional analysis methods. These results provide novel insights into the maturation of conifer SE and indicate that mid-infrared spectrometry could be an easy and effective method for SE characterization.

Comprehensive Analysis of the NF-YB Gene Family and Expression under Abiotic Stress and Hormone Treatment in Larix kaempferi.

NF-YB, a subfamily of Nuclear Factor Y (NF-Y) transcription factor, play crucial role in many biological processes of plant growth and development and abiotic stress responses, and they can therefore be good candidate factors for breeding stress-resistant plants. However, the NF-YB proteins have not yet been explored in Larix kaempferi, a tree species with high economic and ecological values in northeast China and other regions, limiting the breeding of anti-stress L. kaempferi. In order to explore the roles of NF-YB transcription factors in L. kaempferi, we identified 20 LkNF-YB family genes from L. kaempferi full-length transcriptome data and carried out preliminary characterization of them through series of analyses on their phylogenetic relationships, conserved motif structure, subcellular localization prediction, GO annotation, promoter cis-acting elements as well as expression profiles under treatment of phytohormones (ABA, SA, MeJA) and abiotic stresses (salt and drought). The LkNF-YB genes were classified into three clades through phylogenetic analysis and belong to non-LEC1 type NF-YB transcription factors. They have 10 conserved motifs; all genes contain a common motif, and their promoters have various phytohormones and abiotic stress related cis-acting elements. Quantitative real time reverse transcription PCR (RT-qPCR) analysis showed that the sensitivity of the LkNF-YB genes to drought and salt stresses was higher in leaves than roots. The sensitivity of LKNF-YB genes to ABA, MeJA, SA stresses was much lower than that to abiotic stress. Among the LkNF-YBs, LkNF-YB3 showed the strongest responses to drought and ABA treatments. Further protein interaction prediction analysis for LkNF-YB3 revealed that LkNF-YB3 interacts with various factors associated with stress responses and epigenetic regulation as well as NF-YA/NF-YC factors. Taken together, these results unveiled novel L. kaempferi NF-YB family genes and their characteristics, providing the basic knowledge for further in-depth studies on their roles in abiotic stress responses of L. kaempferi.

LkARF7 and LkARF19 overexpression promote adventitious root formation in a heterologous poplar model by positively regulating LkBBM1.

Cuttage propagation involves adventitious root formation induced by auxin. In our previous study, Larix kaempferi BABY BOOM 1 (LkBBM1), which is known to regulate adventitious root formation, was affected by auxin. However, the relationship between LkBBM1 and auxin remains unclear. Auxin response factors (ARFs) are a class of important transcription factors in the auxin signaling pathway and modulate the expression of early auxin-responsive genes by binding to auxin response elements. In the present study, we identified 14 L. kaempferi ARFs (LkARFs), and found LkARF7 and LkARF19 bound to LkBBM1 promoter and enhanced its transcription using yeast one-hybrid, ChIP-qPCR, and dual-luciferase assays. In addition, the treatment with naphthalene acetic acid promoted the expression of LkARF7 and LkARF19. We also found that overexpression of these two genes in poplar promoted adventitious root formation. Furthermore, LkARF19 interacted with the DEAD-box ATP-dependent RNA helicase 53-like protein to form a heterodimer to regulate adventitious root formation. Altogether, our results reveal an additional regulatory mechanism underlying the control of adventitious root formation by auxin.

Antioxidant Properties of Larch Tannins with Different Mean Polymerization Degrees: Controlled Degradation Based on Hydroxyl Radical Degradation.

Hydroxyl radical produced by hydrogen peroxide decomposition under UV radiation was used to degrade larch tannins in an environmentally friendly manner. The formaldehyde reactivity of the degraded products was used as an index to control the mean degree of polymerization (mDP) of the degraded products, and the effects of different mDP on the antioxidant activity of tannins were studied. Results showed that hydroxyl radical could significantly reduce the degree of polymerization (DP) and molecular weight (Mw) of larch tannins, and the mDP and Mw of degraded products could be controlled by considering the formaldehyde reactivity as the index. The antioxidant activity of larch tannins increased with the decrease in mDP. When the degradation time was 6 h, the formaldehyde reactivity was the highest at 0.823. The antioxidant activity of the degraded product was excellent, and the free radical scavenging rate was more than 98%.

The Larix kaempferi genome reveals new insights into wood properties.

Here, through single-molecule real-time sequencing, we present a high-quality genome sequence of the Japanese larch (Larix kaempferi), a conifer species with great value for wood production and ecological afforestation. The assembled genome is 10.97 Gb in size, harboring 45,828 protein-coding genes. Of the genome, 66.8% consists of repeat sequences, of which long terminal repeat retrotransposons are dominant and make up 69.86%. We find that tandem duplications have been responsible for the expansion of genes involved in transcriptional regulation and stress responses, unveiling their crucial roles in adaptive evolution. Population transcriptome analysis reveals that lignin content in L. kaempferi is mainly determined by the process of monolignol polymerization. The expression values of six genes (LkCOMT7, LkCOMT8, LkLAC23, LkLAC102, LkPRX148, and LkPRX166) have significantly positive correlations with lignin content. These results indicated that the increased expression of these six genes might be responsible for the high lignin content of the larches' wood. Overall, this study provides new genome resources for investigating the evolution and biological function of conifer trees, and also offers new insights into wood properties of larches.

Functional nuclear retention of pre-mRNA involving Cajal bodies during meiotic prophase in European larch (Larix decidua).

Gene regulation ensures that the appropriate genes are expressed at the proper time. Nuclear retention of incompletely spliced or mature mRNAs is emerging as a novel, previously underappreciated layer of posttranscriptional regulation. Studies on this phenomenon indicated that it exerts a significant influence on the regulation of gene expression by regulating export and translation delay, which allows the synthesis of specific proteins in response to a stimulus or at strictly controlled time points, for example, during cell differentiation or development. Here, we show that transcription in microsporocytes of European larch (Larix decidua) occurs in a pulsatile manner during prophase of the first meiotic division. Transcriptional activity was then silenced after each pulse. However, the transcripts synthesized were not exported immediately to the cytoplasm but were retained in the nucleoplasm and Cajal bodies (CBs). In contrast to the nucleoplasm, we did not detect mature transcripts in CBs, which only stored nonfully spliced transcripts with retained introns. Notably, the retained introns were spliced at precisely defined times, and fully mature mRNAs were released into the cytoplasm for translation. As similar processes have been observed during spermatogenesis in animals, our results illustrate an evolutionarily conserved mechanism of gene expression regulation during generative cells development in Eukaryota.

LaDAL1 Coordinates Age and Environmental Signals in the Life Cycle of Larix kaempferi.

Perennial woody plants are long-lived, and their life-cycle events occur in order in each generation, but what drives the occurrence and restart of these events in their offspring is unknown. Based on its age-dependent expression pattern and function, Larix kaempferi DEFICIENS-AGAMOUS-LIKE 1 (LaDAL1), a MADS transcription factor has been suggested to be a time recorder and life-cycle event coordinator. Here, we studied the dynamic spatiotemporal expression pattern of LaDAL1 in the life cycle of L. kaempferi to analyze the molecular mechanism of life-cycle progression. In full view of the life cycle, LaDAL1 transcription was related with life-cycle progression, and its transcript level increased sharply from age 3 to 5 years, which might be the molecular characteristic of the vegetative phase change, and then stayed at a high level. During sexual reproduction, LaDAL1 transcript level decreased sequentially during meiosis and embryogenesis, suggesting that meiosis rapidly lowers the age signal, and after fertilization, the age signal was reset to "0" with the embryogenesis. When a seed germinates, the next generation restarts, and the age is re-counted. Altogether, these results not only provide important and novel insights into the life-cycle progression and transgeneration in perennial woody plants, but also advance our understanding of age recording.

Glutathione, carbohydrate and other metabolites of Larix olgensis A. Henry reponse to polyethylene glycol-simulated drought stress.

Drought stress in trees limits their growth, survival, and productivity and it negatively affects the afforestation survival rate. Our study focused on the molecular responses to drought stress in a coniferous species Larix olgensis A. Henry. Drought stress was simulated in one-year-old seedlings using 25% polyethylene glycol 6000. The drought stress response in these seedlings was assessed by analyzing select biochemical parameters, along with gene expression and metabolite profiles. The soluble protein content, peroxidase activity, and malondialdehyde content of L. olgensis were significantly changed during drought stress. Quantitative gene expression analysis identified a total of 8172 differentially expressed genes in seedlings processed after 24 h, 48 h, and 96 h of drought stress treatment. Compared with the gene expression profile of the untreated control, the number of up-regulated genes was higher than that of down-regulated genes, indicating that L. olgensis mainly responded to drought stress through positive regulation. Metabolite analysis of the control and stress-treated samples showed that under drought stress, the increased abundance of linoleic acid was the highest among up-regulated metabolites, which also included some saccharides. A combined analysis of the transcriptome and metabolome revealed that genes dominating the differential expression profile were involved in glutathione metabolism, galactose metabolism, and starch and sucrose metabolism. Moreover, the relative abundance of specific metabolites of these pathways was also altered. Thus, our results indicated that L. olgensis prevented free radical-induced damage through glutathione metabolism and responded to drought through sugar accumulation.

Root foraging of birch and larch in heterogeneous soil nutrient patches under water deficit.

Water and nutrient are two critical factors that limit plant growth to spatial-temporal extents. Tree root foraging behavior has not received adequate attention in heterogeneous soil environments in temperate forest under drought pressure. In this study, birch (Betula platyphylla) and larch (Larix olgensis) seedlings were raised in pots in a split-root system with artificially heterogeneous soil environments to study the root foraging response to drought. Potted space was split into two halves where substrates were mixed with fertilizers in 67.5 mg nitrogen (N) plant-1 (N-P2O5-K2O, 14-13-13) to both halves as to create a homogeneous condition. Otherwise, a rate of 135 mg N plant-1 of fertilizers was delivered to a random half to create a heterogeneous condition. Half of seedlings were fully sub-irrigated every three days with the other half received the drought treatment by being watered every six days. Both birch and larch seedlings showed greater net shoot growth and biomass increment in well-watered condition, while root morphology was promoted by drought. Both species placed more fine roots with higher root N concentration in nutrient-enriched patches. In the heterogeneous pattern, birch showed a higher foraging precision assessed by biomass and greater foraging plasticity assessed in morphology and physiology. In contrast, larch seedlings had higher root N concentration in the well-watered condition. Neither species showed a significant response of N utilization to the heterogeneous pattern, but both used more N when water supply was improved. Overall, birch is better at acclimating to heterogeneous soil conditions, but its ability to seize N was lower than larch when drought was alleviated.

Different Patterns of mRNA Nuclear Retention during Meiotic Prophase in Larch Microsporocytes.

Recent studies show a crucial role of post-transcriptional processes in the regulation of gene expression. Our research has shown that mRNA retention in the nucleus plays a significant role in such regulation. We studied larch microsporocytes during meiotic prophase, characterized by pulsatile transcriptional activity. After each pulse, the transcriptional activity is silenced, but the transcripts synthesized at this time are not exported immediately to the cytoplasm but are retained in the cell nucleus and especially in Cajal bodies, where non-fully-spliced transcripts with retained introns are accumulated. Analysis of the transcriptome of these cells and detailed analysis of the nuclear retention and transport dynamics of several mRNAs revealed two main patterns of nuclear accumulation and transport. The majority of studied transcripts followed the first one, consisting of a more extended retention period and slow release to the cytoplasm. We have shown this in detail for the pre-mRNA and mRNA encoding RNA pol II subunit 10. In this pre-mRNA, a second (retained) intron is posttranscriptionally spliced at a precisely defined time. Fully mature mRNA is then released into the cytoplasm, where the RNA pol II complexes are produced. These proteins are necessary for transcription in the next pulse to occur.mRNAs encoding translation factors and SERRATE followed the second pattern, in which the retention period was shorter and transcripts were rapidly transferred to the cytoplasm. The presence of such a mechanism in various cell types from a diverse range of organisms suggests that it is an evolutionarily conserved mechanism of gene regulation.

Concerted control of the LaRAV1-LaCDKB1;3 module by temperature during dormancy release and reactivation of larch.

Dormancy release and reactivation of temperate-zone trees involve the temperature-modulated expression of cell-cycle genes. However, information on the detailed regulatory mechanism is limited. Here, we compared the transcriptomes of the stems of active and dormant larch trees, emphasizing the expression patterns of cell-cycle genes and transcription factors and assessed their relationships and responses to temperatures. Twelve cell-cycle genes and 31 transcription factors were strongly expressed in the active stage. Promoter analysis suggested that these 12 genes might be regulated by transcription factors from 10 families. Altogether, 73 cases of regulation between 16 transcription factors and 12 cell-cycle genes were predicted, while the regulatory interactions between LaMYB20 and LaCYCB1;1, and LaRAV1 and LaCDKB1;3 were confirmed by yeast one-hybrid and dual-luciferase assays. Last, we found that LaRAV1 and LaCDKB1;3 had almost the same expression patterns during dormancy release and reactivation induced naturally or artificially by temperature, indicating that the LaRAV1-LaCDKB1;3 module functions in the temperature-modulated dormancy release and reactivation of larch trees. These results provide new insights into the link between temperature and cell-cycle gene expression, helping to understand the temperature control of tree growth and development in the context of climate change.

Spatial and Temporal Distribution of Arabinogalactan Proteins during Larix decidua Mill. Male Gametophyte and Ovule Interaction.

The role of ArabinoGalactan Proteins (AGPs) in the sexual reproduction of gymnosperms is not as well documented as that of angiosperms. In earlier studies, we demonstrated that AGPs play important roles during ovule differentiation in Larix decidua Mill. The presented results encouraged us to carry out further studies focused on the functions of these unique glycoproteins during pollen/pollen tube and ovule interactions in Larix. We identified and analyzed the localization of AGPs epitopes by JIM4, JIM8, JIM13 and LM2 antibodies (Abs) in male gametophytes and ovule tissue during pollination, the progamic phase, and after fertilization and in vitro growing pollen tubes. Our results indicated that (1) AGPs recognized by JIM4 Abs play an essential role in the interaction of male gametophytes and ovules because their appearance in ovule cells is induced by physical contact between reproductive partners; (2) after pollination, AGPs are secreted from the pollen cytoplasm into the pollen wall and contact the extracellular matrix of stigmatic tip cells followed by micropylar canal cells; (3) AGPs synthesized in nucellus cells before pollen grain germination are secreted during pollen tube growth into the extracellular matrix, where they can directly interact with male gametophytes; (4) in vitro cultured pollen tube AGPs labeled with LM2 Abs participate in the germination of pollen grain, while AGPs recognized by JIM8 Abs are essential for pollen tube tip growth.

Microbial community composition in the rhizosphere of Larix decidua under different light regimes with additional focus on methane cycling microorganisms.

Microbial community and diversity in the rhizosphere is strongly influenced by biotic and/or abiotic factors, like root exudates, nutrient availability, edaphon and climate. Here we report on the microbial diversity within the rhizosphere of Larix decidua, a dominant tree species in the Alps, as compared with the microbiome within the surrounding soil. We describe how increased light intensity influenced the rhizobiome and put emphasize on methane cycling microorganisms. Microbial taxa were classified into 26 bacterial, 4 archaeal and 6 fungal phyla revealing significant differences between bulk and rhizosphere soils. The dominant prokaryotic phyla were Proteobacteria, Acidobacteria, Actinobacteria (both, rhizosphere and bulk soil) and Bacteroidetes (rhizosphere soil only) and dominant fungal phyla in both fractions included Ascomycota and Basidiomycota. The rhizosphere community was indicated by Suillus sp., plant growth-promoting bacteria and Candidatus Saccharibacteria. Predicted genes in membrane transport and carbohydrate metabolism were significantly more abundant in rhizosphere soils while genes connected with energy metabolisms and cell motility increased in bulk soils. Dominant methanotrophic microorganisms were Upland Soil Cluster (USC) α methanotrophs, Methylogaea spp. and Methylosinus spp., while most methanogens belonged to Methanomassiliicoccales. The overall abundance of methanotrophs distinctly increased in the rhizosphere but to a very different species-specific extent. The increased light intensity only led to minor changes in the rhizobiome, nevertheless a couple of indicator species (e.g. Pseudomonas sp.) for intensified light conditions were established.

Core spliceosomal Sm proteins as constituents of cytoplasmic mRNPs in plants.

In recent years, research has increasingly focused on the key role of post-transcriptional regulation of messenger ribonucleoprotein (mRNP) function and turnover. As a result of the complexity and dynamic nature of mRNPs, the full composition of a single mRNP complex remains unrevealed and mRNPs are poorly described in plants. Here we identify canonical Sm proteins as part of the cytoplasmic mRNP complex, indicating their function in the post-transcriptional regulation of gene expression in plants. Sm proteins comprise an evolutionarily ancient family of small RNA-binding proteins involved in pre-mRNA splicing. The latest research indicates that Sm could also impact on mRNA at subsequent stages of its life cycle. In this work we show that in the microsporocyte cytoplasm of Larix decidua, the European larch, Sm proteins accumulate within distinct cytoplasmic bodies, also containing polyadenylated RNA. To date, several types of cytoplasmic bodies involved in the post-transcriptional regulation of gene expression have been described, mainly in animal cells. Their role and molecular composition in plants remain less well established, however. A total of 222 mRNA transcripts have been identified as cytoplasmic partners for Sm proteins. The specific colocalization of these mRNAs with Sm proteins within cytoplasmic bodies has been confirmed via microscopic analysis. The results from this work support the hypothesis, that evolutionarily conserved Sm proteins have been adapted to perform a whole repertoire of functions related to the post-transcriptional regulation of gene expression in Eukaryota. This adaptation presumably enabled them to coordinate the interdependent processes of splicing element assembly, mRNA maturation and processing, and mRNA translation regulation, and its degradation.

Identification of miRNAs and their target genes in Larix olgensis and verified of differential expression miRNAs.

BACKGROUND: MiRNAs (microRNA) are 18-24 nt endogenous noncoding RNAs that regulate gene expression at the post-transcriptional level, including tissue-specific, developmental timing and evolutionary conservation gene expression. RESULTS: This study used high-throughput sequencing technology for the first time in Larix olgensis, predicted 78 miRNAs, including 12,229,003 reads sRNA, screened differentially expressed miRNAs. Predicting target genes was helpful for understanding the miRNA regulation function and obtained 333 corresponding target genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional annotation were analysed, mostly including nucleic acid binding, plant hormone signal transduction, pantothenate and CoA biosynthesis, and cellulose synthase. This study will lay the foundation for clarifying the complex miRNA-mediated regulatory network for growth and development. In view of this, spatio-temporal expression of miR396, miR950, miR164, miR166 and miR160 were analysed in Larix olgensis during the growth stages of not lignified, beginning of lignification, and completely lignified in different tissues (root, stem, and leaf) by quantitative real-time PCR (qRT-PCR). There were differences in the expression of miRNAs in roots, stems and leaves in the same growth period. At 60 days, miR160, miR166 and miR396-2 exhibited the highest expression in leaves. At 120 days, most miRNAs in roots and stems decreased significantly. At 180 days, miRNAs were abundantly expressed in roots and stems. Meanwhile, analysis of the expression of miRNAs in leaves revealed that miR396-2 was reduced as time went on, whereas other miRNAs increased initially and then decreased. On the other hand, in the stems, miR166-1 was increase, whereas other miRNAs, especially miR160, miR164, miR396 and miR950-1, first decreased and then increased. Similarly, in the roots, miR950-2 first decreased and then increased, whereas other miRNAs exhibited a trend of continuous increase. CONCLUSIONS: The present investigation included rapid isolation and identification of miRNAs in Larix olgensis through construction of a sRNA library using Solexa and predicted 78 novel miRNAs, which showed differential expression levels in different tissues and stages. These results provided a theoretical basis for further revealing the genetic regulation mechanism of miRNA in the growth and development of conifers and the verification of function in target genes.

Long-term effects of post-fire restoration types on nitrogen mineralisation in a Dahurian larch (Larix gmelinii) forest in boreal China.

Fire is an important disturbance agent in Chinese boreal forests but the long-term effects of wildfires on soil nitrogen (N) net mineralisation rates (Rmin) in natural versus human-assisted restorations are not well understood. In this study, we analysed upper (0-10cm) and lower layer (10-20cm) soil samples from natural restoration and afforestation plots in a Dahurian larch (Larix gmelinii) forest in north-eastern China 29years after a mega fire disturbance. Our results showed that the soil inorganic N (NH4+-N and NO3--N) pool of the upper and lower layers of the regenerated plots remained significantly lower than in unburned control plots. This suggests that the effects of a high burn severity fire on soil N availability were still significant almost 30years after the event. Restoration type (natural restoration versus afforestation) also had significant effects on upper layer soil N availability; compared with afforestation, natural restoration was more beneficial for the accumulation of soil inorganic N and the recovery of Rmin after fire disturbance. Specifically, the concentration of inorganic N and the mean Rmin in upper layer soils in the natural restoration plots were approximately 41% greater and 3.6 times greater, respectively, than in the afforestation plots. The differences in soil N availability between the two restoration types were attributed to differences in soil water content (SWC), soil microbial biomass nitrogen (MBN), and the recovery of vegetation after the fire disturbance. Our study demonstrates that natural restoration has proved more successful than afforestation in countering soil N losses in boreal forests in China resulting from a high burn severity fire disturbance.