Isolation and functional analysis of the Larix olgensis LoNAC3 transcription factor gene.

BACKGROUND: Larch is an important timber tree species. The traditional methods of tree genetic breeding have been progressing slowly. It is necessary to carry out gene function analysis and genetically modified breeding research. The NAC transcription factor family is a plant-specific transcription factor family with various biological functions, as shown in recent research. However, there are few studies on the NAC gene among gymnosperm coniferous species. RESULTS: LoNAC3 with complete cds was identified and isolated from the cDNA of Larix olgensis based on transcriptome data. The cDNA length of LoNAC3 is 1185 bp, encoding 394 amino acids, with a conserved NAM domain located at the N-terminus, and subcellular localization in the nucleus. The results of real-time quantitative PCR analysis showed that at different growth stages and in different tissues of L. olgensis, the relative expression level of LoNAC3 was highest in the needles. After drought, salt, alkali stress and hormone treatment, expression was induced to different degrees. The expression level of LoNAC3 was significantly increased under drought and salt conditions. The relative expression level changed under methyl jasmonate (MeJA) and abscisic acid (ABA) treatment. By observing the phenotype of overexpressed LoNAC3 tobacco, it was found that overexpressed tobacco is shorter and blooms earlier than wild-type tobacco. Under abiotic stress, LoNAC3 overexpressed tobacco has lower germination rates and poorer growth status. Transgenic tobacco under stress treatment has a higher malondialdehyde (MDA) content than wild-type tobacco, while peroxidase (POD) activity is lower than wild-type tobacco. CONCLUSIONS: Through the analysis of LoNAC3 sequence and promoter expression, it can be concluded that LoNAC3 is involved in the drought and salt stress response processes of L. olgensis, and is induced by ABA and MeJA expression. Overexpression of LoNAC3 leads to stunted tobacco growth and negatively regulates its tolerance to drought and salt stress through the reactive oxygen species pathway. The preliminary analysis of the expression pattern and function of the LoNAC3 can provide a theoretical basis and high-quality materials for genetic improvement of larch in later stages.

Differential phenotypic plasticity of subalpine trees predicts trait integration under climate warming.

Understanding limiting factors of phenotypic plasticity is essential given its critical role in shaping biological adaptation and evolution in changing environments. It has been proposed that the pattern of phenotypic correlation could constrain trait plasticity. However, the interplay between phenotypic plasticity and integration has remained contentious. We experimentally simulated climate warming in juveniles of three subalpine tree species by exposing them to three-year in situ open-top chambers (OTCs), and then measured functional plasticity of 72 eco-physiological traits to evaluate whether phenotypic integration constituted an intrinsic constraint to plasticity. We also tested the relationship between the differences in plasticity and maintenance in trait integration. Phenotypic plasticity was positively associated with integration in deciduous tree species under warming. The difference in the plasticity of two paired traits could predict their integration in different environments, where traits displaying more similar plasticity were more likely to be correlated. Our study showed no indication that phenotypic integration constrained plasticity. More importantly, we demonstrated that differential plasticity between traits might result in a notable reorganization of the trait associations, and that warming commonly induced a tighter phenotype. Our study provides new insights into the interplay between phenotypic plasticity and integration in subalpine trees under climate warming.

Anatomical and morphological changes in Pinus sylvestris and Larix sibirica needles under impact of emissions from a large aluminum enterprise.

Species-specific anatomical and morphological characteristics of Pinus sylvestris and Larix sibirica needles were studied at different levels of tree stand pollution by aluminum smelter emissions. The anatomical characteristics of the needle were studied using light microscopy. The level of tree stand pollution was determined using the cluster analysis outcomes of the pollutant elements content (fluorine, sulfur, and heavy metals) in the needles. Four levels of tree stand pollution were separated: low, moderate, high, and critical, as well as background tree stand in unpolluted areas. It was found that the state of tree phytomass deteriorated with increasing levels of pollution (from low to critical): pine crown defoliation increased to 85%, and larch defoliation increased to 65%. The life span of pine needles was reduced to 2-3 years, with a background value of 6-7 years. The change of morphological parameters was more pronounced in P. sylvestris: the weight and length of the 2-year-old shoot decreased by 2.7-3.1 times compared to the background values; the weight of needles on the shoot and the number of needle pairs on the shoot-by 1.9-2.1 times. The length of the needle and shoot and the number of L. sibirica brachyblasts decreased by 1.8-1.9 times. The anatomical parameters of the needle also changed to a greater extent in P. sylvestris. Up to the high level of tree pollution, we observed a decrease in the cross-sectional area of the needle, central cylinder, vascular bundle, area and thickness of mesophyll, number and diameter of resin ducts by 18-66% compared to background values. At the critical pollution level, when the content of pollutant elements in pine needles reached maximum values, the anatomical parameters of the remaining few green needles were close to background values. In our opinion, this may be due to the activation of mechanisms aimed at maintaining the viability of trees. A reduction in thickness and area of assimilation tissue in the L. sibirica needle was detected only at the critical pollution level. An upward trend in these parameters was found at low, medium, and high pollution levels of tree stand, which may indicate an adaptive nature. The results suggested that at a similar pollution level of trees, the greatest amount of negative anatomical and morphological changes were recorded in pine needles, which indicates a greater sensitivity of this species to technogenic emissions.

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.

Needle and Branch Trait Variation Analysis and Associated SNP Loci Mining in Larix olgensis.

Needles play key roles in photosynthesis and branch growth in Larix olgensis. However, genetic variation and SNP marker mining associated with needle and branch-related traits have not been reported yet. In this study, we examined 131 samples of unrelated genotypes from L. olgensis provenance trails. We investigated phenotypic data for seven needle and one branch-related traits before whole genome resequencing (WGRS) was employed to perform a genome-wide association study (GWAS). Subsequently, the results were used to screen single nucleotide polymorphism (SNP) loci that were significantly correlated with the studied traits. We identified a total of 243,090,868 SNP loci, and among them, we discovered a total of 161 SNP loci that were significantly associated with these traits using a general linear model (GLM). Based on the GWAS results, Kompetitive Allele-Specific PCR (KASP), designed based on the DNA of population samples, were used to validate the loci associated with L. olgensis phenotypes. In total, 20 KASP markers were selected from the 161 SNPs loci, and BSBM01000635.1_4693780, BSBM01000114.1_5114757, and BSBM01000114.1_5128586 were successfully amplified, were polymorphic, and were associated with the phenotypic variation. These developed KASP markers could be used for the genetic improvement of needle and branch-related traits in L. olgensis.

Large-Scale Quantitative Proteomic Analysis during Different Stages of Somatic Embryogenesis in Larix olgensis.

Larix olgensis is an economically important tree species native to northeastern China. The use of somatic embryogenesis (SE) is efficient and enables the rapid production of varieties with desirable qualities. Here, isobaric labeling via tandem mass tags was used to conduct a large-scale quantitative proteomic analysis of proteins in three critically important stages of SE in L. olgensis: the primary embryogenic callus, the single embryo, and the cotyledon embryo. We identified 6269 proteins, including 176 shared differentially expressed proteins across the three groups. Many of these proteins are involved in glycolipid metabolism, hormone response/signal transduction, cell synthesis and differentiation, and water transport; proteins involved in stress resistance and secondary metabolism, as well as transcription factors, play key regulatory roles in SE. The results of this study provide new insights into the key pathways and proteins involved in SE in Larix. Our findings have implications for the expression of totipotency, the preparation of synthetic seeds, and genetic transformation.

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.

Predictive accuracy of post-fire conifer death declines over time in models based on crown and bole injury.

A key uncertainty of empirical models of post-fire tree mortality is understanding the drivers of elevated post-fire mortality several years following fire, known as delayed mortality. Delayed mortality can represent a substantial fraction of mortality, particularly for large trees that are a conservation focus in western US coniferous forests. Current post-fire tree mortality models have undergone limited evaluation of how injury level and time since fire interact to influence model accuracy and predictor variable importance. Less severe injuries potentially serve as an indicator for vulnerability to additional stressors such as bark beetle attack or moisture stress. We used a collection of 164,293 individual tree records to examine post-fire tree mortality in eight western USA conifers: Abies concolor, Abies grandis, Calocedrus decurrens, Larix occidentalis, Pinus contorta, Pinus lambertiana, Pinus ponderosa, and Pseudotsuga menziesii. We evaluated the importance of fire injury predictors on discriminating between surviving trees versus immediate and delayed post-fire mortality. We fit balanced random forest models for each species using cumulative tree mortality from 1 to 5-years post-fire. We compared these results to multi-class random forest models using first-year mortality, 2-5-year mortality, and survival 5-years post-fire as a response variable. Crown volume scorched, diameter at breast height, and relative bark char height, were used as predictor variables. The cumulative mortality models all predicted trees that died within 1-year of fire with high accuracy but failed to predict 2-5-year mortality. The multi-class models were an improvement but had lower accuracy for predicting 2-5-year mortality. Multi-class model accuracies ranged from 85% to 95% across all species for predicting 1-year post-fire mortality, 42%-71% for predicting 2-5-year mortality, and 64%-85% for predicting trees that lived past 5-years. Our study highlights the differences in tree species tolerance to fire injury and suggests that including second-order predictors such as beetle attack or climatic water stress before and after fire will be critical to improve accuracy and better understand the mechanisms and patterns of fire-caused tree death. Random forest models have potential for management applications such as post-fire harvesting and simulating future stand dynamics.

Turpentine Ointment for the Treatment of Folliculitis: An Open, Prospective, Randomized, Placebo- and Comparator-Controlled Multicenter Trial.

INTRODUCTION: Folliculitis is a painful infection and inflammation of the hair follicles, mostly caused by bacterial, fungal, or, more rarely, viral infections. Turpentine derivatives have been used traditionally to treat various skin infections and could thus also be effective in treating folliculitis. We carried out an open, prospective, randomized, placebo- and comparator-controlled multicenter trial to evaluate the efficacy and safety of an ointment containing pine turpentine oil, larch turpentine, and eucalyptus oil in the treatment of acute folliculitis. METHODS: Seventy outpatients with acute folliculitis were treated with the turpentine ointment, a comparator (povidone iodine solution), or a placebo (Vaseline) for 7 days. Photographs of the affected skin areas were taken by the physicians at four visits and by the patients on a daily basis. Photographs were evaluated by blinded observers. Primary efficacy endpoint was the change in total hair follicle lesion counts. Secondary endpoints included the evolution of the lesion counts in the course of the study, responder rate (improvement of follicle lesions by at least one count), and the patient's global assessment. Safety endpoints were the tolerability of the treatments and adverse event recording. RESULTS: A decrease of follicle lesions counts was detected for both active treatments but not for placebo, but the differences among groups were not statistically significant. As for the secondary endpoints, the ointment showed statistically significant superiority over placebo for the evolution of the lesions during the course of the study (p = 0.017), the responder rate (p = 0.032), and the subjective efficacy assessment by patients (p = 0.029). All treatments were equally well tolerated, with a similar number of treatment-emergent adverse events. CONCLUSION: The turpentine ointment is an effective and safe option for the treatment of folliculitis.

Evolutionary adaptation of trees and modelled future larch forest extent in Siberia.

With changing climate, the boreal forest could potentially migrate north and become threatened by droughts in the south. However, whether larches, the dominant tree species in eastern Siberia, can adapt to novel situations is largely unknown but is crucial for predicting future population dynamics. Exploring variable traits and trait adaptation through inheritance in an individual-based model can improve our understanding and help future projections. We updated the individual-based spatially explicit vegetation model LAVESI (Larix Vegetation Simulator), used for forest predictions in Eastern Siberia, with trait value variation and incorporated inheritance of parental values to their offspring. Forcing the model with both past and future climate projections, we simulated two areas - the expanding northern treeline and a southerly area experiencing drought. While the specific trait of 'seed weight' regulates migration, the abstract 'drought resistance' protects stands. We show that trait variation with inheritance leads to an increase in migration rate (∼ 3% area increase until 2100). The drought resistance simulations show that, under increasing stress, including adaptive traits leads to larger surviving populations (17% of threatened under RCP 4.5 (Representative Concentration Pathway)). We show that with the increase expected under the RCP 8.5 scenario vast areas (80% of the extrapolated area) of larch forest are threatened and could disappear due to drought as adaptation plays only a minor role under strong warming. We conclude that variable traits facilitate the availability of variants under environmental changes. Inheritance allows populations to adapt to environments and promote successful traits, which leads to populations that can spread faster and be more resilient, provided the changes are not too drastic in both time and magnitude. We show that trait variation and inheritance contribute to more accurate models that can improve our understanding of responses of boreal forests to global change.

Larches of Kuzhanovo Have a Unique Mutation in the atpF-atpH Intergenic Spacer.

The larches of Kuzhanovo (Larix sibirica Ledeb.) are protected trees with a round crown growing in the Southern Urals. In 2020 vandals sawed the sapwood of these trees, which exposed the problem of insufficient conservation measures. Their origin and genetic characteristics have been of particular interest to breeders and scientists. The larches of Kuzhanovo were screened for polymorphisms using SSR and ISSR analyses and the sequencing of genetic markers and genes GIGANTEA and mTERF, associated with wider crown shape. A unique mutation was discovered in the atpF-atpH intergenic spacer of all protected trees, but it was absent in some of their descendants and larches with similar crown shape. Mutations were discovered in the rpoC1 and mTERF genes of all samples. Flow cytometry did not reveal any changes in genome size. Our results suggest that the unique phenotype arose from point mutations in L. sibirica, but they are yet to be found in the nuclear genome. The concurrent mutations in the rpoC1 and mTERF genes may indicate that the round crown shape originates from the Southern Urals. The atpF-atpH and rpoC1 genetic markers are not common in studies of Larix sp., but their wider use could help to establish the origin of these endangered plants. The discovery of the unique atpF-atpH mutation also allows for stronger conservation and crime detection efforts.

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.

Genetic and Epigenetic Mechanisms of Longevity in Forest Trees.

Trees are unique in terms of development, sustainability and longevity. Some species have a record lifespan in the living world, reaching several millennia. The aim of this review is to summarize the available data on the genetic and epigenetic mechanisms of longevity in forest trees. In this review, we have focused on the genetic aspects of longevity of a few well-studied forest tree species, such as Quercus robur, Ginkgo biloba, Ficus benghalensis and F. religiosa, Populus, Welwitschia and Dracaena, as well as on interspecific genetic traits associated with plant longevity. A key trait associated with plant longevity is the enhanced immune defense, with the increase in gene families such as RLK, RLP and NLR in Quercus robur, the expansion of the CC-NBS-LRR disease resistance families in Ficus species and the steady expression of R-genes in Ginkgo biloba. A high copy number ratio of the PARP1 family genes involved in DNA repair and defense response was found in Pseudotsuga menziesii, Pinus sylvestris and Malus domestica. An increase in the number of copies of the epigenetic regulators BRU1/TSK/MGO3 (maintenance of meristems and genome integrity) and SDE3 (antiviral protection) was also found in long-lived trees. CHG methylation gradually declines in the DAL 1 gene in Pinus tabuliformis, a conservative age biomarker in conifers, as the age increases. It was shown in Larix kaempferi that grafting, cutting and pruning change the expression of age-related genes and rejuvenate plants. Thus, the main genetic and epigenetic mechanisms of longevity in forest trees were considered, among which there are both general and individual processes.

Genetic Adaptation of Siberian Larch (Larix sibirica Ledeb.) to High Altitudes.

Forest trees growing in high altitude conditions offer a convenient model for studying adaptation processes. They are subject to a whole range of adverse factors that are likely to cause local adaptation and related genetic changes. Siberian larch (Larix sibirica Ledeb.), whose distribution covers different altitudes, makes it possible to directly compare lowland with highland populations. This paper presents for the first time the results of studying the genetic differentiation of Siberian larch populations, presumably associated with adaptation to the altitudinal gradient of climatic conditions, based on a joint analysis of altitude and six other bioclimatic variables, together with a large number of genetic markers, single nucleotide polymorphisms (SNPs), obtained from double digest restriction-site-associated DNA sequencing (ddRADseq). In total, 25,143 SNPs were genotyped in 231 trees. In addition, a dataset of 761 supposedly selectively neutral SNPs was assembled by selecting SNPs located outside coding regions in the Siberian larch genome and mapped to different contigs. The analysis using four different methods (PCAdapt, LFMM, BayeScEnv and RDA) revealed 550 outlier SNPs, including 207 SNPs whose variation was significantly correlated with the variation of some of environmental factors and presumably associated with local adaptation, including 67 SNPs that correlated with altitude based on either LFMM or BayeScEnv and 23 SNPs based on both of them. Twenty SNPs were found in the coding regions of genes, and 16 of them represented non-synonymous nucleotide substitutions. They are located in genes involved in the processes of macromolecular cell metabolism and organic biosynthesis associated with reproduction and development, as well as organismal response to stress. Among these 20 SNPs, nine were possibly associated with altitude, but only one of them was identified as associated with altitude by all four methods used in the study, a nonsynonymous SNP in scaffold_31130 in position 28092, a gene encoding a cell membrane protein with uncertain function. Among the studied populations, at least two main groups (clusters), the Altai populations and all others, were significantly genetically different according to the admixture analysis based on any of the three SNP datasets as follows: 761 supposedly selectively neutral SNPs, all 25,143 SNPs and 550 adaptive SNPs. In general, according to the AMOVA results, genetic differentiation between transects or regions or between population samples was relatively low, although statistically significant, based on 761 neutral SNPs (FST = 0.036) and all 25,143 SNPs (FST = 0.017). Meanwhile, the differentiation based on 550 adaptive SNPs was much higher (FST = 0.218). The data showed a relatively weak but highly significant linear correlation between genetic and geographic distances (r = 0.206, p = 0.001).

Comparative transcriptomic responses of European and Japanese larches to infection by Phytophthora ramorum.

BACKGROUND AND OBJECTIVES: Phytophthora ramorum severely affects both European larch (EL) and Japanese larch (JL) trees as indicated by high levels of mortality particularly in the UK. Field observations suggested that EL is less severely affected and so may be less susceptible to P. ramorum than JL; however, controlled inoculations have produced inconsistent or non-statistically significant differences. The present study aimed to compare RNA transcript accumulation profiles in EL and JL in response to inoculation with P. ramorum to improve our understanding of their defence responses. METHODOLOGY: RNA-sequencing was carried out on bark tissues following the inoculation with P. ramorum of potted saplings in both EL and JL carried out under controlled environment conditions, with sampling at 1, 3, 10, and 25 days post inoculation in infected and control plants. RESULTS: All of the inoculated trees rapidly developed lesions but no statistically significant differences were found in lesion lengths between EL and JL. RNA-Sequencing comparing control and inoculate saplings identified key differences in differentially expressed genes (DEGs) between the two larch species. European larch had rapid induction of defence genes within 24 hours of infection followed by sustained expression until 25 days after inoculation. Results in JL were more varied; upregulation was stronger but more transient and represented fewer defence pathways. Gene enrichment analyses highlighted differences in jasmonate signalling and regulation including NPR1 upregulation in EL only, and specific aspects of secondary metabolism. Some DEGs were represented by multiple responsive copies including lipoxygenase, chalcone synthase and nucleotide-binding, leucine-rich-repeat genes. CONCLUSION: The variations between EL and JL in responsive DEGs of interest as potentially related to differences seen in the field and should be considered in the selection of trees for planting and future breeding.

Pedigree reconstruction and spatial analysis for genetic testing and selection in a Larix kaempferi (Lamb.) Carrière plantation.

BACKGROUND: Larix kaempferi is one of the major timber species in Northeast Asia. Demand for the reforestation of the species is rising in South Korea due to an increase in large timber production and utilization. However, progeny trials for the species have not been explored, making it challenging to foster advanced generations of tree improvement. In the present study, genetic testing and selection for diameter growth were conducted using pedigree reconstruction and phenotypic spatial distribution analysis in a plantation of L. kaempferi. The aim of the present study was to select the superior larch individuals using the pedigree reconstruction and phenotypic spatial distribution to substitute progeny trials. The plantation of seed orchard crops was established in 1990 and one-hundred and eighty-eight trees were selected as the study material. Genetic variation was investigated first to validate its adequacy as breeding material. Genetic testing was carried out using a model considering pedigree information and spatial autoregression of the phenotypes. RESULTS: The expected heterozygosity of the mother trees and offspring were 0.672 and 0.681 presenting the corresponding level of genetic variation between two groups. The pedigree reconstruction using maternity analysis assigned one to six progenies to ninety-two candidate mothers. The accuracy of genetic testing was exceedingly increased with the animal model considering AR1 ⊗ AR1 structure compared to the animal model only. The estimated genetic variance of the former was 9.086 whereas that of the latter was 4.9E-5 for DBH. The predicted breeding values of the offspring for DBH were ranged from -5.937 cm to 5.655 cm and the estimated heritability of diameter growth was 0.344. CONCLUSIONS: The genetic testing approach based on pedigree reconstruction and phenotypic spatial distribution analysis was considered a useful analytical scheme that could replace or supplement progeny trials.

Comprehensive collection of genes and comparative analysis of full-length transcriptome sequences from Japanese larch (Larix kaempferi) and Kuril larch (Larix gmelinii var. japonica).

BACKGROUND: Japanese larch (Larix kaempferi) is an economically important deciduous conifer species that grows in cool-temperate forests and is endemic to Japan. Kuril larch (L. gmelinii var. japonica) is a variety of Dahurian larch that is naturally distributed in the Kuril Islands and Sakhalin. The hybrid larch (L. gmelinii var. japonica × L. kaempferi) exhibits heterosis, which manifests as rapid juvenile growth and high resistance to vole grazing. Since these superior characteristics have been valued by forestry managers, the hybrid larch is one of the most important plantation species in Hokkaido. To accelerate molecular breeding in these species, we collected and compared full-length cDNA isoforms (Iso-Seq) and RNA-Seq short-read, and merged them to construct candidate gene as reference for both Larix species. To validate the results, candidate protein-coding genes (ORFs) related to some flowering signal-related genes ​were screened from the reference sequences, and the phylogenetic relationship with closely related species was elucidated. RESULTS: Using the isoform sequencing of PacBio RS ll and the de novo assembly of RNA-Seq short-read sequences, we identified 50,690 and 38,684 ORFs in Japanese larch and Kuril larch, respectively. BUSCO completeness values were 90.5% and 92.1% in the Japanese and Kuril larches, respectively. After comparing the collected ORFs from the two larch species, a total of 19,813 clusters, comprising 22,571 Japanese larch ORFs and 22,667 Kuril larch ORFs, were contained in the intersection of the Venn diagram. In addition, we screened several ORFs related to flowering signals (SUPPRESSER OF OVEREXPRESSION OF CO1: SOC1, LEAFY: LFY, FLOWERING Locus T: FT, CONSTANCE: CO) from both reference sequences, and very similar found in other species. CONCLUSIONS: The collected ORFs will be useful as reference sequences for molecular breeding of Japanese and Kuril larches, and also for clarifying the evolution of the conifer genome and investigating functional genomics.

Temporal dynamics of fine root production, mortality and turnover deviate across branch orders in a larch stand.

Fine roots play a key role in carbon, nutrient, and water biogeochemical cycles in forest ecosystems. However, inter-annual dynamics of fine root production, mortality, and turnover on the basis of long-term measurement have been less studied. Here, field scanning rhizotrons were employed for tracking fine root by branch order over a 6 years period in a larch plantation. For total fine roots, from the first- to the fifth-order roots, annual root length production, length mortality, standing crops, and turnover rate varied up to 3.4, 2.3, 1.5, and 2.3-folds during the study period, respectively. The inter-annual variability of those roots indices in the first-order and the second-order roots were greater than that of the higher order (third- to fifth-order) roots. The turnover rate was markedly larger for the first-order roots than for the higher order roots, showing the greatest variability up to 20 times. Seasonal dynamics of root length production followed a general concentrated pattern with peak typically occurring in June or July, whereas root length mortality followed a general bimodal mortality pattern with the dominant peak in May and the secondary peak in August or October. Furthermore, the seasonal patterns of root length production and mortality were similar across years, especially for the first-order and the second-order roots. These results from long-term observation were beneficial for reducing uncertainty of characterizing fine root demography in consideration of large variation among years. Our findings highlight it is important for better understanding of fine root dynamics and determining root demography through distinguishing observation years and root branch orders.

Activity of low-molecular weight components of Larix sibirica antioxidant system under exposure to technogenic pollution.

Changes in the antioxidant protection system of Larix sibirica Ledeb at different pollution levels caused by emissions from a large aluminum smelter (BrAS) have been studied. We revealed that the content of peroxide (H2O2) in the needles is a reliable marker of oxidative stress in the trees under pollution. The crucial role of non-enzymatic components, in particular, proline, phenolic compounds, ascorbic acid, glutathione, in reducing the level of free radicals in the needles cells was found. Proline concentration in the needles significantly rises with the increase in pollution levels from low to high. Under critical level pollution, it decreases by 40% compared to the background. The total content of ascorbic acid (ASC) in the needles of polluted trees varies slightly; however, there are significant changes in its various forms. With an increase in pollution to a high level, the content of the reduced form of ASC in the needles increases by 1.5-2.9 times compared to the background content. At a critical level of pollution, the total level of ascorbic acid and its reduced form falls, the content of the oxidized form reaches minimum values. The total content of phenolic compounds in the needles increased by 50-55%, concentration of flavonoids by 1.5-1.8 times, catechins by 1.9-2.5 times, proanthocyanidins by 45% compared to the background level under low, moderate, high pollution, whereas under critical pollution their content decreased. The absolute concentration of the reduced form glutathione in the needles falls by 1.9-3.0 times, the oxidized form increases by 1.5-2.0 times compared to the background. The ratio of reduced glutathione to oxidized glutathione decreased, especially during critical pollution. The data obtained show significant activation of Siberian larch biochemical protection at low, moderate and high levels of pollution by the aluminum smelter emissions. At a critical levels of contamination, a significant depletion of the pool of low-molecular antioxidants was observed.

Comparative Genomics of Seasonal Senescence in Forest Trees.

In the course of evolution, both flowering plants and some gymnosperms have developed such an adaptation to winter and unfavorable living conditions as deciduousness. Of particular interest is Siberian larch (Larix sibirica Ledeb.), which is the only species in the pine family (Pinaceae) with a seasonal deciduousness. New generation sequencing technologies make it possible to study this phenomenon at the genomic level and to reveal the genetic mechanisms of leaf and needle aging in angiosperms and gymnosperms. Using a comparative analysis of the genomes of evergreen and deciduous trees, it was found that the genes that control EXORDIUM LIKE 2 (EXL2) and DORMANCY-ASSOCIATED PROTEIN 1 (DRM1) proteins are most represented in Siberian larch, while an excess of genes that control proteins acting as immune receptors were found in evergreens. Orthologs from the family of genes that control leucine-rich repeat receptor-like kinases (LRR-RLK) contributed mostly to the distinction between evergreens and deciduous plants.