Genetic analysis and QTL mapping for pericarp thickness in maize (Zea mays L.).

Proper pericarp thickness protects the maize kernel against pests and diseases, moreover, thinner pericarp improves the eating quality in fresh corn. In this study, we aimed to investigate the dynamic changes in maize pericarp during kernel development and identified the major quantitative trait loci (QTLs) for maize pericarp thickness. It was observed that maize pericarp thickness first increased and then decreased. During the growth and formation stages, the pericarp thickness gradually increased and reached the maximum, after which it gradually decreased and reached the minimum during maturity. To identify the QTLs for pericarp thickness, a BC4F4 population was constructed using maize inbred lines B73 (recurrent parent with thick pericarp) and Baimaya (donor parent with thin pericarp). In addition, a high-density genetic map was constructed using maize 10 K SNP microarray. A total of 17 QTLs related to pericarp thickness were identified in combination with the phenotypic data. The results revealed that the heritability of the thickness of upper germinal side of pericarp (UG) was 0.63. The major QTL controlling UG was qPT1-1, which was located on chromosome 1 (212,215,145-212,948,882). The heritability of the thickness of upper abgerminal side of pericarp (UA) was 0.70. The major QTL controlling UA was qPT2-1, which was located on chromosome 2 (2,550,197-14,732,993). In addition, a combination of functional annotation, DNA sequencing analysis and quantitative real-time PCR (qPCR) screened two candidate genes, Zm00001d001964 and Zm00001d002283, that could potentially control maize pericarp thickness. This study provides valuable insights into the improvement of maize pericarp thickness during breeding.

Coordinated control for the auricle asymmetric development by ZmIDD14 and ZmIDD15 fine-tune the high-density planting adaption in maize.

Multiple distinct specialized regions shape the architecture of maize leaves. Among them, the fringe-like and wedge-shaped auricles alter the angle between the leaf and stalk, which is a key trait in crop plant architecture. As planting density increased, a small leaf angle (LA) was typically selected to promote crop light capture efficiency and yield. In the present study, we characterized two paralogous INDETERMINATE DOMAIN (IDD) genes, ZmIDD14 and ZmIDD15, which contain the Cys2-His2 zinc finger domain and function redundantly to regulate auricle development and LA in maize. Loss-of-function mutants showed decreased LA by reducing adaxial sclerenchyma thickness and increasing the colourless cell layers. In addition, the idd14;idd15 double mutant exhibited asymmetrically smaller auricles, which might cause by a failed maintenance of symmetric expression of the key auricle size controlling gene, LIGULELESS(LG1). The transcripts of ZmIDD14 and ZmIDD15 enriched in the ligular region, where LG1 was highly expressed, and both proteins physically interacted with ZmILI1 to promote LG1 transcription. Notably, the idd14;idd15 enhanced the grain yield of hybrids under high planting densities by shaping the plant architecture with a smaller LA. These findings demonstrate the functions of ZmIDD14 and ZmIDD15 in controlling the abaxial/adaxial development of sclerenchyma in the midrib and polar development along the medial-lateral axes of auricles and provide an available tool for high-density and high-yield breeding in maize.

A critical thermal transition driving spring phenology of Northern Hemisphere conifers.

Despite growing interest in predicting plant phenological shifts, advanced spring phenology by global climate change remains debated. Evidence documenting either small or large advancement of spring phenology to rising temperature over the spatio-temporal scales implies a potential existence of a thermal threshold in the responses of forests to global warming. We collected a unique data set of xylem cell-wall-thickening onset dates in 20 coniferous species covering a broad mean annual temperature (MAT) gradient (-3.05 to 22.9°C) across the Northern Hemisphere (latitudes 23°-66° N). Along the MAT gradient, we identified a threshold temperature (using segmented regression) of 4.9 ± 1.1°C, above which the response of xylem phenology to rising temperatures significantly decline. This threshold separates the Northern Hemisphere conifers into cold and warm thermal niches, with MAT and spring forcing being the primary drivers for the onset dates (estimated by linear and Bayesian mixed-effect models), respectively. The identified thermal threshold should be integrated into the Earth-System-Models for a better understanding of spring phenology in response to global warming and an improved prediction of global climate-carbon feedbacks.

Variations in the TAS2R38 gene among college students in Hubei.

BACKGROUND: The bitter taste receptor gene TAS2R38 is a member of the human TAS2R gene family. Polymorphisms in TAS2R38 affect the ability to taste the bitterness of phenylthiourea (PTC) compounds, thus affecting an individual's food preference and health status. METHODS: We investigated polymorphisms in the TAS2R38 gene and the sensitivity to PTC bitterness among healthy Chinese college students in Hubei province. The association of TAS2R38 polymorphisms and PTC sensitivity with body mass index (BMI), food preference, and health status was also analyzed. A total of 320 healthy college students were enrolled (male: 133, female: 187; aged 18-23 years). The threshold value method was used to measure the perception of PTC bitterness, and a questionnaire was used to analyze dietary preferences and health status. Polymerase chain reaction (PCR) was used to analyze polymorphisms at three common TAS2R38 loci (rs713598, rs1726866, and rs10246939). RESULTS: In our study population, 65.00% of individuals had medium sensitivity to the bitterness of PTC; in contrast, 20.94% were highly sensitive to PTC bitterness, and 14.06% were not sensitive. For the TAS2R38 gene, the PAV/PAV and PAV/AAI diplotypes were the most common (42.19% and 40.63%, respectively), followed by the homozygous AVI/AVI (8.75%) and PAV/AVI (5.00%) diplotypes. CONCLUSION: There was a significant correlation between the sensitivity to PTC bitterness and sex, but there was no correlation between the common diplotypes of TAS2R38 and gender. Polymorphisms in the TAS2R38 gene were associated with the preference for tea, but not with one's native place, BMI, health status, or other dietary preferences. There was no significant correlation between the perception of PTC bitterness and one's native place, BMI, dietary preference, or health status. We hope to find out the relationship between PTC sensitivity and TAS2R38 gene polymorphisms and dietary preference and health status of Chinese population through this study, providing relevant guidance and suggestions for dietary guidance and prevention of some chronic diseases in Chinese population.

Genome-Wide Identification and Characterization of Actin-Depolymerizing Factor (ADF) Family Genes and Expression Analysis of Responses to Various Stresses in Zea Mays L.

Actin-depolymerizing factor (ADF) is a small class of actin-binding proteins that regulates the dynamics of actin in cells. Moreover, it is well known that the plant ADF family plays key roles in growth, development and defense-related functions. Results: Thirteen maize (Zea mays L., ZmADFs) ADF genes were identified using Hidden Markov Model. Phylogenetic analysis indicated that the 36 identified ADF genes in Physcomitrella patens, Arabidopsis thaliana, Oryza sativa japonica, and Zea mays were clustered into five groups. Four pairs of segmental genes were found in the maize ADF gene family. The tissue-specific expression of ZmADFs and OsADFs was analyzed using microarray data obtained from the Maize and Rice eFP Browsers. Five ZmADFs (ZmADF1/2/7/12/13) from group V exhibited specifically high expression in tassel, pollen, and anther. The expression patterns of 13 ZmADFs in seedlings under five abiotic stresses were analyzed using qRT-PCR, and we found that the ADFs mainly responded to heat, salt, drought, and ABA. Conclusions: In our study, we identified ADF genes in maize and analyzed the gene structure and phylogenetic relationships. The results of expression analysis demonstrated that the expression level of ADF genes was diverse in various tissues and different stimuli, including abiotic and phytohormone stresses, indicating their different roles in plant growth, development, and response to external stimulus. This report extends our knowledge to understand the function of ADF genes in maize.

Genome-wide identification of AGO18b-bound miRNAs and phasiRNAs in maize by cRIP-seq.

BACKGROUND: Argonaute proteins (AGOs) are important players in the regulation of plant development by directing sRNAs to target mRNAs. In maize (Zea mays), AGO18b is a tassel-enriched and grass-specific AGO. Previous studies have shown that AGO18b is highly expressed in tassels during meiosis and negatively regulates determinacy of spikelet meristems. However, binding profile on RNAs and acting mechanisms of AGO18b remain unknown. RESULTS: In this study, we explored the binding profile of AGO18b in maize tassel by UV cross-linking RNA immunoprecipitation, followed by deep sequencing of these cDNA libraries (cRIP-seq), and systematically studied AGO18b-associated small RNAs and mRNAs by bioinformatics analysis. By globally analyzing the phased small-interfering RNA (phasiRNA) and miRNA abundance bound by AGO18b, we found AGO18b primarily binds to 21-nt phasiRNAs/miRNAs with a 5'-uridine and binds less strongly to 24-nt phasiRNAs with a 5'-adenosine in the premeiotic tassels. The abundance profile of AGO18b-associated miRNAs was different from their expression profile. Moreover, AGO18b strongly binds to miR166a-3p. We then obtained the AGO18b-bound mRNA targets of miR166a-3p by cRIP-seq, and confirmed the molecular function of AGO18b in regulating spikelet meristems. CONCLUSIONS: Our results indicated that AGO18b binds to phasiRNAs with obvious 5 prime end bias under different sRNA length. MiRNAs and their target mRNAs associated with AGO18b indicated the molecular mechanisms of AGO18b as a negative regulator of inflorescence meristem and tassel development through integrating both phasiRNAs and miRNA pathways, which extended our view of sRNA regulation in flower development and provided potential methods to control pollination in the future.

Genome-wide identification, phylogenetic and expression analysis of the maize HECT E3 ubiquitin ligase genes.

HECT (homologous to the E6AP carboxyl terminus) ubiquitin ligase genes (E3s) are enzymes with diverse functions influencing plant growth, development, and responses to abiotic stresses. However, there is relatively little information available regarding the maize HECT E3 gene family. In the present study, 12 maize HECT E3 genes (ZmUPL1 to ZmUPL12) were identified at the whole-genome level. The phylogenetic relationships, structures, and expression levels of the maize HECT E3 genes were then analyzed. On the basis of the constructed maximum likelihood phylogenetic tree, the HECT E3 genes were divided into six groups. The quantitative real-time polymerase chain reaction assay results revealed that all of the maize ZmUPL genes were expressed in most of the examined tissues and were responsive to three abiotic stresses. Considered together, the study results may provide a useful foundation for future investigations of maize stress-tolerance genes as well as functional analyses of the E3 enzymes in diverse agriculturally important crop species.

Expression analysis of Argonaute genes in maize (Zea mays L.) in response to abiotic stress.

BACKGROUND: Argonaute (AGO) protein is a kind of RNA binding protein that plays an integral role in the gene-silencing pathways guided by small RNAs. But there are few studies about the regulation of AGO genes responded to diverse abiotic stress in maize. RESULTS: In this study, we analyzed the expression of seventeen ZmAGO genes under heat, cold, salinity, drought and ABA treatments using quantitative PCR (qPCR). All ZmAGOs showed differential expression modes under various abiotic stress treatments. Two ZmAGOs (ZmAGO1a and ZmAGO5d) and other fifteen ZmAGOs exhibited specific up-regulation in response to heat separately. Several ZmAGO genes are very sensitive to cold stress, but many ZmAGO genes are slow to respond to NaCl treatment. Nine ZmAGO genes (ZmAGO1f, ZmAGO2b, ZmAGO4, ZmAGO5a/b/c, ZmAGO7, ZmAGO9 and ZmAGO18a/b) presented definite up-regulation in response to drought, which were similar to the pattern of gene regulation under abscisic acid (ABA) treatment. CONCLUSIONS: Various ZmAGO genes respond to different abiotic stress treatments. These results provide fundamental information and insights for the further study on the role of abiotic stress resistance genes in maize and provide basis for further study on the function of AGO genes in response to abiotic stress in maize.

Contributions of insects and droughts to growth decline of trembling aspen mixed boreal forest of western Canada.

Insects, diseases, fire and drought and other disturbances associated with global climate change contribute to forest decline and mortality in many parts of the world. Forest decline and mortality related to drought or insect outbreaks have been observed in North American aspen forests. However, little research has been done to partition and estimate their relative contributions to growth declines. In this study, we combined tree-ring width and basal area increment series from 40 trembling aspen (Populus tremuloides Michx.) sites along a latitudinal gradient (from 52° to 58°N) in western Canada and attempted to investigate the effect of drought and insect outbreaks on growth decline, and simultaneously partition and quantify their relative contributions. Results indicated that the influence of drought on forest decline was stronger than insect outbreaks, although both had significant effects. Furthermore, the influence of drought and insect outbreaks showed spatiotemporal variability. In addition, our data suggest that insect outbreaks could be triggered by warmer early spring temperature instead of drought, implicating that potentially increased insect outbreaks are expected with continued warming springs, which may further exacerbate growth decline and death in North America aspen mixed forests.

AGO18b negatively regulates determinacy of spikelet meristems on the tassel central spike in maize.

The maize tassel represents an indeterminate male inflorescence. The number of primordia that a given inflorescence meristem produces is related to its determinacy, i.e., capacity for continued meristem activity. Transcription factors (TFs) controlling determinacy in tassel axillary meristems are well studied in maize, and small RNAs are known to influence tassel development by repressing targets, including tassel-related TFs. As core components of the RNA-inducible silence complex (RISC), Argonaute (AGO) proteins are required for small RNA-mediated repression. Here, we characterized the biological function of AGO18b, a tassel-enriched AGO. The abundance of AGO18b transcripts gradually increased during tassel development from inception to gametogenesis and were enriched in the inflorescence meristem and axillary meristems of the tassel. Repressing AGO18b expression resulted in more spikelets, which contributed to a longer central spike of the tassel. Additionally, the transcripts of several HD-ZIP III TFs that were canonical targets of microRNA166 (miR166) accumulated in the AGO18b-repressed lines. We propose that AGO18b is a negative regulator of the determinacy of inflorescence and axillary meristems, and that it acts by interacting with the miR166-HD-ZIP III TF regulatory pathway.

Drought causes reduced growth of trembling aspen in western Canada.

Adequate and advance knowledge of the response of forest ecosystems to temperature-induced drought is critical for a comprehensive understanding of the impacts of global climate change on forest ecosystem structure and function. Recent massive decline in aspen-dominated forests and an increased aspen mortality in boreal forests have been associated with global warming, but it is still uncertain whether the decline and mortality are driven by drought. We used a series of ring-width chronologies from 40 trembling aspen (Populus tremuloides Michx.) sites along a latitudinal gradient (from 52° to 58°N) in western Canada, in an attempt to clarify the impacts of drought on aspen growth by using Standardized Precipitation Index (SPI) and Standardized Precipitation Evapotranspiration Index (SPEI). Results indicated that prolonged and large-scale droughts had a strong negative impact on trembling aspen growth. Furthermore, the spatiotemporal variability of drought indices is useful for explaining the spatial heterogeneity in the radial growth of trembling aspen. Due to ongoing global warming and rising temperatures, it is likely that severer droughts with a higher frequency will occur in western Canada. As trembling aspen is sensitive to drought, we suggest that drought indices could be applied to monitor the potential effects of increased drought stress on aspen trees growth, achieve classification of eco-regions and develop effective mitigation strategies to maintain western Canadian boreal forests.

Argonaute and Argonaute-Bound Small RNAs in Stem Cells.

Small RNAs are essential for a variety of cellular functions. Argonaute (AGO) proteins are associated with all of the different classes of small RNAs, and are indispensable in small RNA-mediated regulatory pathways. AGO proteins have been identified in various types of stem cells in diverse species from plants and animals. This review article highlights recent progress on how AGO proteins and AGO-bound small RNAs regulate the self-renewal and differentiation of distinct stem cell types, including pluripotent, germline, somatic, and cancer stem cells.

ZmGA3ox2, a candidate gene for a major QTL, qPH3.1, for plant height in maize.

Maize plant height is closely associated with biomass, lodging resistance and grain yield. Determining the genetic basis of plant height by characterizing and cloning plant height genes will guide the genetic improvement of crops. In this study, a quantitative trait locus (QTL) for plant height, qPH3.1, was identified on chromosome 3 using populations derived from a cross between Zong3 and its chromosome segment substitution line, SL15. The plant height of the two lines was obviously different, and application of exogenous gibberellin A(3) removed this difference. QTL mapping placed qPH3.1 within a 4.0 cM interval, explaining 32.3% of the phenotypic variance. Furthermore, eight homozygous segmental isolines (SILs) developed from two larger F(2) populations further narrowed down qPH3.1 to within a 12.6 kb interval. ZmGA3ox2, an ortholog of OsGA3ox2, which encodes a GA3 ß-hydroxylase, was positionally cloned. Association mapping identified two polymorphisms in ZmGA3ox2 that were significantly associated with plant height across two experiments. Quantitative RT-PCR showed that SL15 had higher ZmGA3ox2 expression relative to Zong3. The resultant higher GA(1) accumulation led to longer internodes in SL15 because of increased cell lengths. Moreover, a large deletion in the coding region of ZmGA3ox2 is responsible for the dwarf mutant d1-6016. The successfully isolated qPH3.1 enriches our knowledge on the genetic basis of plant height in maize, and provides an opportunity for improvement of plant architecture in maize breeding.

Identification and characterization of Argonaute gene family and meiosis-enriched Argonaute during sporogenesis in maize.

Argonaute (AGO) proteins play a key role in regulation of gene expression through small RNA-directed RNA cleavage and translational repression, and are essential for multiple developmental processes. In the present study, 17 AGO genes of maize (Zea mays L., ZmAGOs) were identified using a Hidden Markov Model and validated by rapid amplification of cDNA ends assay. Subsequently, quantitative PCR revealed that expressions of these genes were higher in reproductive than in vegetative tissues. AGOs presented five temporal and spatial expression patterns, which were likely modulated by DNA methylation, 5'-untranslated exons and microRNA-mediated feedback loops. Intriguingly, ZmAGO18b was highly expressed in tassels during meiosis. Furthermore, in situ hybridization and immunofluorescence showed that ZmAGO18b was enriched in the tapetum and germ cells in meiotic anthers. We hypothesized that ZmAGOs are highly expressed in reproductive tissues, and that ZmAGO18b is a tapetum and germ cell-specific member of the AGO family in maize.

High preseason temperature variability drives convergence of xylem phenology in the Northern Hemisphere conifers.

Wood growth is key to understanding the feedback of forest ecosystems to the ongoing climate warming. An increase in spatial synchrony (i.e., coincident changes in distant populations) of spring phenology is one of the most prominent climate responses of forest trees. However, whether temperature variability contributes to an increase in the spatial synchrony of spring phenology and its underlying mechanisms remains largely unknown. Here, we analyzed an extensive dataset of xylem phenology observations of 20 conifer species from 75 sites over the Northern Hemisphere. Along the gradient of increase in temperature variability in the 75 sites, we observed a convergence in the onset of cell enlargement roughly toward the 5th of June, with a convergence in the onset of cell wall thickening toward the summer solstice. The increase in rainfall since the 5th of June is favorable for cell division and expansion, and as the most hours of sunlight are received around the summer solstice, it allows the optimization of carbon assimilation for cell wall thickening. Hence, the convergences can be considered as the result of matching xylem phenological activities to favorable conditions in regions with high temperature variability. Yet, forest trees relying on such consistent seasonal cues for xylem growth could constrain their ability to respond to climate warming, with consequences for the potential growing season length and, ultimately, forest productivity and survival in the future.

Genome-wide identification and analysis of microRNA responding to long-term waterlogging in crown roots of maize seedlings.

MicroRNAs (miRNAs) are critical post-transcriptional modulators of gene expression involving in plant responses to abiotic stress. However, the regulation of miRNA in the morphological response to waterlogging is poorly understood in maize. In this study, we detected miRNAs and their targets that expressed in waterlogged crown roots of maize seedlings in two inbred lines (Hz32 and Mo17) by RNA sequencing. A total of 61 mature miRNAs were found including 36 known maize (zma) miRNAs and 25 potential novel miRNA candidates. Comparison of miRNA expression in both waterlogged and control crown roots revealed 32 waterlogging-responsive miRNAs, most were consistently downregulated under waterlogging in the two inbred lines. We identified the miRNA targets through degradome sequencing. Many known miRNA targets involving in transcription regulation and reactive oxygen species elimination were found in the degradome libraries, and 17 targets of 10 newly detected miRNAs were identified as well. Moreover, the miRNA-mediated pathways that respond to waterlogging and regulate the induction of crown roots were discussed. This study is a comprehensive survey of responsive miRNAs in waterlogged maize crown roots. The results will help to understand the miRNA expression in response to waterlogging and miRNA-mediated regulation of morphological adaptation to waterlogging in maize.

Kirenol inhibits inflammation challenged by lipopolysaccharide through the AMPK-mTOR-ULK1 autophagy pathway.

Kirenol is a bioactive substance isolated from Herba Siegesbeckiae. Although the anti-inflammatory activity of kirenol has been well documented, its role in autophagy remains unknown. The present study aimed to investigate the protective role of kirenol on inflammation challenged by lipopolysaccharide (LPS) in acute lung injury (ALI) cell and mouse models and unravel the underlying mechanisms, with a particular focus on autophagy. For this purpose, an ALI cell and mouse models were established, and the effects of kirenol on the expression of molecules related to inflammation and autophagy were examined. The present results revealed that kirenol could significantly inhibit inflammatory cytokines secretion in cells and in the mice injured by LPS; this effect may be attributed to enhanced autophagy as evidenced by the up-regulation of LC3-II and the down-regulation of p62 both in vitro and in vivo. Phosphorylated AMPK and ULK1 increased, while phosphorylated mTOR decreased in the kirenol-treated ALI cell model. Moreover, inhibition of autophagy using AMPK inhibitor or 3-MA or chloroquine (CQ) reversed the anti-inflammatory and autophagy-enhancement effects of kirenol exposure in vitro, indicating that kirenol could enhance autophagy by activating the AMPK-mTOR-ULK1 pathway. The results of RNA sequencing suggested that kirenol was strongly related to the biological functions of acute inflammatory response and the AMPK signaling pathway. Further in vivo ALI mouse model studies demonstrated the protective role of kirenol against lung inflammation, such as improved histopathology, decreased lung edema, and leukocyte infiltration were abolished by 3-MA. These findings implicate that kirenol can inhibit LPS-induced inflammation via the AMPK-mTOR-ULK1 autophagy pathway.

Pilots' Spatial Visualization Ability Assessment Based on Virtual Reality.

BACKGROUND: The aim of this study was to investigate the effectiveness of the mental rotation test (MRT) based on virtual reality (VR) in predicting pilots' spatial visualization ability (SVA).METHODS: Based on VR, 118 healthy pilots' SVA were evaluated by MRT. The pilot flight ability evaluation scale was used as the criterion of test validity. According to the scale score, pilots were divided into high, middle, or low spatial ability groups pursuant to the 27% allocation principle. Differences in reaction time (RT), correct rate (CR), and correct number per second (CNPS) of MRT between groups were compared. Correlations between scale scores and MRT scores were analyzed. RT, CR, and CNPS of MRT among different age groups and between genders were also compared.RESULTS: The RT of the high spatial ability group was remarkably slower than that of the low spatial ability group (363.4 ± 140.2 s, 458.1 ± 151.7 s). The CNPS of the high spatial ability group was dramatically higher than that of the low spatial ability group (0.111 ± 0.045 s, 0.086 ± 0.001 s). There were no significant differences in RT, CR, and CNPS between different genders. Pilots in the 29-35 yr old age group had considerably slower RT than those in the 22-28 yr old age group (330.8 ± 140.3 s, 417.2 ± 132.7 s). Pilots in the 29-35 yr old age group had conspicuously higher CNPS than pilots in the 22-28 yr old age group (0.119 ± 0.040 s, 0.096 ± 0.036 s). All pilots' scale scores were positively correlated with CNPS (r = 0.254) and negatively correlated with RT (r = -0.234).DISCUSSION: MRT based on VR has a good discrimination efficacy for SVA of pilots and is a good indicator for the SVA component measurement.Zhang M, Wang M, Feng H, Liu X, Zhai L, Xu X, Jin Z. Pilots' spatial visualization ability assessment based on virtual reality. Aerosp Med Hum Perform. 2023; 94(6):422-428.

Variation in intra-annual wood formation, and foliage and shoot development of three major Canadian boreal tree species.

PREMISE OF THE STUDY: In a warming climate, boreal trees may have adjusted their growth strategy (e.g., onset and coordination of growth among different organs such as stem, shoot, and foliage, within and among species) to cope with the extended growing seasons. A detailed investigation into growth of different organs during a growing season may help us assess the potential effects of climate change on tree growth in the boreal forest. METHODS: The intra-annual growth of stem xylem, shoot tips, and foliage area of Pinus banksiana, Populus tremuloides, and Betula papyrifera was monitored in a boreal forest in Quebec, Canada during the growing season of 2007. Xylem formation was measured at weekly intervals, and shoot elongation and foliage expansion were measured three times per week from May to September. Growth indices for stem, shoot, and foliage were calculated and used to identify any climate-growth dependence. KEY RESULTS: The time periods required for stem growth, branch extension, and foliage expansion differed among species. Of the three species, P. banksiana had the earliest budburst (20 May) yet the latest completion date of the foliage growth (2 August); P. tremuloides had the latest budburst (27 May) yet the earliest completion date of the foliage growth (10 July). Air temperature positively affected shoot extension growth of all three species. Precipitation positively influenced stem growth of the two broadleaf species, whereas growing season temperature positively impacted stem growth of P. banksiana. CONCLUSION: The results show that both the timing of growth processes and environmental dependences differ among co-occurring species, thereby leading to different adaptive capability of these boreal tree species to climate change.