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.

Genetic transformation of LoHDZ2 and analysis of its function to enhance stress resistance in Larix olgensis.

To study the function of LoHDZ2 in larch, we first constructed a VB191103-LoHDZ2::GUS overexpression vector. Through Agrobacterium-mediated infection, the expression vector was transferred into a larch embryogenic cell line. A stable resistant cell line was subsequently screened, and mature embryos were induced to grow until they developed into seedlings. Antagonistic cell lines were identified at both the DNA and RNA levels. The transgenic cell lines were then subjected to GUS staining, and transgenic cell lines were ultimately identified and obtained. These transgenic cell lines were sequenced to identify differentially expressed genes, and a cluster analysis was performed. The resistant cell lines were cultured under stress conditions involving 20% PEG6000 and 200 mM NaCl proliferation media (1/10-BM). After the stress treatment, the contents of peroxidase (POD), malondialdehyde (MDA) and superoxide dismutase (SOD) in both wild-type and transgenic cell lines were measured. The results are summarized below: (1) When the specific fragment of the target gene in the genome of the resistant cell line was amplified. At the RNA level, the expression of the fragment in four resistant lines increased. In addition, GUS staining showed a blue reaction, indicating that LoHDZ2 was successfully integrated into the larch embryonic cell lines. (2) To verify the accuracy and reliability of the transcriptome data, 10 differentially expressed genes (5 upregulated and 5 down regulated genes) were subjected to qRT-PCR verification. The results showed that the expression trend of the 10 differentially expressed genes was the same as that revealed by RNA-Seq, indicating that the transcriptome data were reliable. (3) The transcriptome sequencing showed that 176 genes were upregulated and that 140 genes were down regulated. Through GO enrichment analysis and KEGG metabolic pathway analysis, the screened differentially expressed genes were related to biological processes such as larch metabolism and response to stimuli, indicating that these genes may be closely involved in the regulation of the larch response to external stimuli, including heat stress, drought stress, metal ion stress and bacterial infection, and may participate in the growth process. (4) After 20% PEG6000 treatment, the POD enzyme activity of the transgenic cell line was greater than that of the wild-type; this activity could effectively remove the amount of peroxide produced. The MDA content of the transgenic cell lines was lower than that of the wild-type cell lines, and the accumulation degree of harmful substances was low, indicating that the degree of oxidative damage of the transgenic cell lines was lower than that of the wild-type cell lines. The SOD content of the transgenic cell lines was lower than that of the wild-type cell lines, indicating that the drought resistance of the transgenic cell lines was enhanced. After 200 mM NaCl treatment, although the increase in SOD content was not obvious, the same trend was detected, indicating that the resistance of the transgenic cell lines was indeed stronger than that of the wild-type cell lines. According to the results of previous experiments, after this gene was overexpressed in tobacco, the transformed plants showed obvious dwarfing, which may indicate that the stress resistance of the plant was enhanced. In conclusion, a transgenic larch cell line was successfully obtained, and transgenic larch seedlings were successfully induced. LoHDZ2 may participate in the response of plants to the external environment, and may participate in the growth and development of Larix olgensis by affecting plant metabolic pathways.

Genetic transformation and growth index determination of the Larix olgensis LoHDZ2 transcription factor gene in tobacco.

Homeodomain-leucine zippers (HD-Zip) are plant-specific transcription factors that participate in different plant development processes and differentially regulate metabolic processes. LoHDZ2 is an HD-ZipII subfamily transcription factor gene that we identified from a transcriptomic analysis of Larix olgensis. To understand its function, we built a LoHDZ2 expression vector and then inserted it into tobacco by genetic transformation. Transgenic plants were identified at the DNA and RNA levels. Phenotypic index analysis of transgenic tobacco showed dwarfed growth with larger leaves and earlier flowering than the wild type. LoHDZ2 was expressed differently after hormone treatment with IAA, MeJA and 2,4-D. The results suggested that LoHDZ2 may respond to hormones and be involved in regulating growth and metabolism. These results helped us better understand the function of LoHDZ2 and provided a candidate gene for Larix olgensis molecular breeding.

Embryogenic callus induction from immature zygotic embryos and genetic transformation of Larix kaempferi 3x Larix gmelinii 9.

To date, there are few reports of the successful genetic transformation of larch and other conifers, mainly because it is difficult to transform and integrate exogenous genes. In this study, hybrid larch Larix kaempferi 3x Larix gmelinii 9 cones were collected on June 27, July 1, July 4, July 7 and July 16, 2017. Embryogenic callus induction was studied using a combination of different plant growth regulators and concentrations. The results showed that July 1 was the best stage; the highest induction rate was 10.83%, which cultured in BM medium (Button medium, which formula was listed in S1 Table) with 1.0 mg/L 2,4-D (2,4-dichlorophenoxyacetic acid) and 0.2 mg/L KT(kinetin). When cultured on a proliferation medium for 12 days, proliferation was the fastest, reaching 323.08%, which could also maintain the freshness and vitality. The suitable pre-culture medium for somatic embryogenesis was 1/4 BM medium containing 10 g/L inositol and 60 g/L sucrose. The combination of 45 mg/L ABA (abscisic acid) and 75 g/L PEG4000 (Polyethyene glycol 4000) could promote the number of somatic embryos, and reached the maximum, 210 140 per 1 g FW. The genetic transformation was carried out by the Agrobacterium-mediated transformation method with embryogenic callus cultured for 12 days. The results showed the optimal OD600 of the infection solution(suspension of A. tumefaciens) was 0.5, co-culture time was 2 days, and screening concentration of Hyg (hygromycin B) was 4 mg/L. In this study, the transformation rate of resistance callus was 32.1%. It provides a reference for low genetic transformation efficiency of larch at present. This study could be beneficial for the innovation and breeding of larch by genetic engineering and provides a certain basis for rapid propagation of excellent larch germplasm resources and genetic engineering breeding of larch and other conifers.