Research Institute of Tropical Forestry

A DNA barcode is a DNA fragment used to identify species. For land plants, DNA fragments of plastid genome could be the primary consideration. Unfortunately, most of the plastid candidate barcodes lack species-level resolution. The identification of DNA barcodes of high resolution at species level is critical to the success of DNA barcoding in plants. We searched the available plastid genomes for the most variable regions and tested the best candidates using both a large number of tree species and seven well-sampled plant groups. Two regions of the plastid gene ycf1, ycf1a and ycf1b, were the most variable loci that were better than existing plastid candidate barcodes and can serve as a barcode of land plants. Primers were designed for the amplification of these regions, and the PCR success of these primers ranged from 82.80% to 98.17%. Of 420 tree species, 357 species could be distinguished using ycf1b, which was slightly better than the combination of matK and rbcL. For the well-sampled representative plant groups, ycf1b generally performed better than any of the matK, rbcL and trnH-psbA. We concluded that ycf1a or ycf1b is the most variable plastid genome region and can serve as a core barcode of land plants.

Global change is impacting forests worldwide, threatening biodiversity and ecosystem services including climate regulation. Understanding how forests respond is critical to forest conservation and climate protection. This review describes an international network of 59 long-term forest dynamics research sites (CTFS-ForestGEO) useful for characterizing forest responses to global change. Within very large plots (median size 25 ha), all stems ≥ 1 cm diameter are identified to species, mapped, and regularly recensused according to standardized protocols. CTFS-ForestGEO spans 25 °S-61 °N latitude, is generally representative of the range of bioclimatic, edaphic, and topographic conditions experienced by forests worldwide, and is the only forest monitoring network that applies a standardized protocol to each of the world's major forest biomes. Supplementary standardized measurements at subsets of the sites provide additional information on plants, animals, and ecosystem and environmental variables. CTFS-ForestGEO sites are experiencing multifaceted anthropogenic global change pressures including warming (average 0.61 °C), changes in precipitation (up to ± 30% change), atmospheric deposition of nitrogen and sulfur compounds (up to 3.8 g N m(-2) yr(-1) and 3.1 g S m(-2) yr(-1)), and forest fragmentation in the surrounding landscape (up to 88% reduced tree cover within 5 km). The broad suite of measurements made at CTFS-ForestGEO sites makes it possible to investigate the complex ways in which global change is impacting forest dynamics. Ongoing research across the CTFS-ForestGEO network is yielding insights into how and why the forests are changing, and continued monitoring will provide vital contributions to understanding worldwide forest diversity and dynamics in an era of global change.

Abstract Aim To examine the contribution of large‐diameter trees to biomass, stand structure, and species richness across forest biomes. Location Global. Time period Early 21st century. Major taxa studied Woody plants. Methods We examined the contribution of large trees to forest density, richness and biomass using a global network of 48 large (from 2 to 60 ha) forest plots representing 5,601,473 stems across 9,298 species and 210 plant families. This contribution was assessed using three metrics: the largest 1% of trees ≥ 1 cm diameter at breast height (DBH), all trees ≥ 60 cm DBH, and those rank‐ordered largest trees that cumulatively comprise 50% of forest biomass. Results Averaged across these 48 forest plots, the largest 1% of trees ≥ 1 cm DBH comprised 50% of aboveground live biomass, with hectare‐scale standard deviation of 26%. Trees ≥ 60 cm DBH comprised 41% of aboveground live tree biomass. The size of the largest trees correlated with total forest biomass ( r 2 = .62, p < .001). Large‐diameter trees in high biomass forests represented far fewer species relative to overall forest richness ( r 2 = .45, p < .001). Forests with more diverse large‐diameter tree communities were comprised of smaller trees ( r 2 = .33, p < .001). Lower large‐diameter richness was associated with large‐diameter trees being individuals of more common species ( r 2 = .17, p = .002). The concentration of biomass in the largest 1% of trees declined with increasing absolute latitude ( r 2 = .46, p < .001), as did forest density ( r 2 = .31, p < .001). Forest structural complexity increased with increasing absolute latitude ( r 2 = .26, p < .001). Main conclusions Because large‐diameter trees constitute roughly half of the mature forest biomass worldwide, their dynamics and sensitivities to environmental change represent potentially large controls on global forest carbon cycling. We recommend managing forests for conservation of existing large‐diameter trees or those that can soon reach large diameters as a simple way to conserve and potentially enhance ecosystem services.

BACKGROUND: The Moso bamboo (Phyllostachys edulis) is one of the most important forestry resources and plays essential ecological roles in southern China. A draft nuclear genome sequence is expected to be publicly available in the near future; an explosion of gene expression data related to the unique traits of Moso bamboo will undoubtedly follow. Reverse transcription quantitative real-time PCR ((RT-)qPCR) is a widely used method for gene expression analysis. A necessary prerequisite of exact and reliable data is the accurate choice of reference genes. RESULT: In this study, 14 candidate reference genes were chosen, and their expression levels were assessed by (RT-)qPCR in a set of six tissue samples (root, stem, mature stem, leaf, flower, and leaf sheath) and at two developmental stages (before and after flowering) in bamboo specimens obtained in three locations. The stability and suitability of the candidate reference genes were validated using the geNorm, NormFinder and BestKeeper programs. The results showed that TIP41 and NTB were suitable reference genes across all the tissues and at the different developmental stages examined in this study. While the expression of the NTB, TIP41 and UBQ were the mostly stable in different plant tissues samples, the expression of the TIP41, NTB and CAC were ranked the most stable in bamboo plants at various developmental stages. AP2-like gene was further assessed by using the reference genes TIP41 and NTB in comparison to ACT. Significant difference of the expression profile of AP2-like demonstrated the importance of choosing adequate reference genes in bamboo. CONCLUSION: TIP41 and NTB were found to be homogeneously expressed and were adequate for normalization purposes, showing equivalent transcript levels in different samples. They are therefore the recommended reference genes for measuring gene expression in P. edulis.

Extreme events often expose vulnerabilities of socioeconomic infrastructures and point to directions of much-needed policy change. Integrated impact assessment of such events can lead to finding of sustainability principles. Southern and central China has for decades been undergoing a breakneck pace of socioeconomic development. In early 2008, a massive ice storm struck this region, immobilizing millions of people. The storm was a consequence of sustained convergence between tropical maritime and continental polar air masses, caused by an anomalously stable atmospheric general circulation pattern in both low and high latitudes. Successive waves of freezing rain occurred during a month period, coating southern and central China with a layer of ice 50–160 mm in thickness. We conducted an integrated impact assessment of this event to determine whether and how the context of socioeconomic and human-disturbed natural systems may affect the transition of natural events into human disasters. We found that 1) without contingency plans, advanced technologies dependent on interrelated energy supplies can create worse problems during extreme events, 2) the weakest link in disaster response lies between science and decision making, 3) biodiversity is a form of long-term insurance for sustainable forestry against extreme events, 4) sustainable extraction of nontimber goods and services is essential to risk planning for extreme events in forest resources use, 5) extreme events can cause food shortage directly by destroying crops and indirectly by disrupting food distribution channels, 6) concentrated economic development increases societal vulnerability to extreme events, and 7) formalized institutional mechanisms are needed to ensure that unexpected opportunities to learn lessons from weather disasters are not lost in distracting circumstances.

The study of biodiversity has tended to focus primarily on relatively information‐poor measures of species diversity. Recently, many studies of local diversity (alpha diversity) have begun to use measures of functional and phylogenetic alpha diversity. Investigations into the phylogenetic and functional dissimilarity (beta diversity) of communities have been far less numerous, but these dissimilarity measures have the potential to infer the mechanisms underlying community assembly and dynamics. Here, we relate levels of phylogenetic and functional alpha diversity to levels of phylogenetic and functional beta diversity to infer the mechanism or mechanisms responsible for the assembly of tree communities in six forests located in tropical and temperate latitudes. The results show that abiotic filtering plays a role in structuring local assemblages and governing spatial turnover in community composition and that phylogenetic measures of alpha and beta diversity are not strong predictors of functional alpha and beta diversity in the forests studied.

Denitrification removes fixed nitrogen (N) from the biosphere, thereby restricting the availability of this key limiting nutrient for terrestrial plant productivity. This microbially driven process has been exceedingly difficult to measure, however, given the large background of nitrogen gas (N2) in the atmosphere and vexing scaling issues associated with heterogeneous soil systems. Here, we use natural abundance of N and oxygen isotopes in nitrate (NO3 (-)) to examine dentrification rates across six forest sites in southern China and central Japan, which span temperate to tropical climates, as well as various stand ages and N deposition regimes. Our multiple stable isotope approach across soil to watershed scales shows that traditional techniques underestimate terrestrial denitrification fluxes by up to 98%, with annual losses of 5.6-30.1 kg of N per hectare via this gaseous pathway. These N export fluxes are up to sixfold higher than NO3 (-) leaching, pointing to widespread dominance of denitrification in removing NO3 (-) from forest ecosystems across a range of conditions. Further, we report that the loss of NO3 (-) to denitrification decreased in comparison to leaching pathways in sites with the highest rates of anthropogenic N deposition.

Abstract There is increasing concern over the impact of atmospheric nitrogen (N) deposition on forest ecosystems in the tropical and subtropical areas. In this study, we quantified atmospheric N deposition and revealed current plant and soil N status in 14 forests along a 150 km urban to rural transect in southern China, with an emphasis on examining whether foliar δ 15 N can be used as an indicator of N saturation. Bulk deposition ranged from 16.2 to 38.2 kg N ha −1 yr −1 , while the throughfall covered a larger range of 11.7–65.1 kg N ha −1 yr −1 . Foliar N concentration, NO 3 − leaching to stream, and soil NO 3 − concentration were low and NO 3 − production was negligible in some rural forests, indicating that primary production in these forests may be limited by N supply. But all these N variables were enhanced in suburban and urban forests. Across the study transect, throughfall N input was correlated positively with soil nitrification and NO 3 − leaching to stream, and negatively with pH values in soil and stream water. Foliar δ 15 N was between −6.6‰ and 0.7‰, and was negatively correlated with soil NO 3 − concentration and NO 3 − leaching to stream across the entire transect, demonstrating that an increased N supply does not necessarily increase forest δ 15 N values. We proposed several potential mechanism that could contribute to the δ 15 N pattern, including (1) increased plant uptake of 15 N‐depleted soil NO 3 − , (2) foliage uptake of 15 N‐depleted NH 4 + , (3) increased utilization of soil inorganic N relative to dissolved organic N, and (4) increased fractionation during plant N uptake under higher soil N availability.

Bacteria are the highest abundant microorganisms in the soil. To investigate bacteria community structures, diversity and functions, contrasting them in four different seasons all the year round with/within two different forest type soils of China. We analyzed soil bacterial community based on 16S rRNA gene sequencing via Illumina HiSeq platform at a temperate deciduous broad-leaved forest (Baotianman, BTM) and a tropical rainforest (Jianfengling, JFL). We obtained 51,137 operational taxonomic units (OTUs) and classified them into 44 phyla and 556 known genera, 18.2% of which had a relative abundance >1%. The composition in each phylum was similar between the two forest sites. Proteobacteria and Acidobacteria were the most abundant phyla in the soil samples between the two forest sites. The Shannon index did not significantly differ among the four seasons at BTM or JFL and was higher at BTM than JFL in each season. The bacteria community at both BTM and JFL showed two significant (P<0.05) predicted functions related to carbon cycle (anoxygenic photoautotrophy sulfur oxidizing and anoxygenic photoautotrophy) and three significant (P<0.05) predicted functions related to nitrogen cycle (nitrous denitrificaton, nitrite denitrification and nitrous oxide denitrification). We provide the basis on how changes in bacterial community composition and diversity leading to differences in carbon and nitrogen cycles at the two forests.

Eucalyptus species are considered to have ectomycorrhizas (ECM), but many also have vesicular–arbuscular mycorrhizas (VAM) and their relative importance is unclear. Interactions between ECM and VAM fungi colonizing roots of Eucalyptus species were examined in a glasshouse experiment. This experiment investigated competition between these two types of fungi and compared benefits provided to Eucalyptus globulus and E. urophylla. Eucalyptus seedlings were inoculated with spores of the ECM fungus Laccaria lateritia and/or pot‐ culture soil for a VAM fungus (species of Glomus , Acaulospora or Scutellospora ). Initial inoculum levels were important, as VAM fungi became established much more rapidly than Laccaria . In plants with both types of mycorrhizas, Laccaria mycorrhizas substantially increased after 2 or 3 months and the proportion of roots with VAM declined. However, the proportion of roots with VAM also decreased significantly in plants without ECM after 2 months. Thus, the relative susceptibility of eucalypt roots to these mycorrhizal associations changed. Substantial growth responses to mycorrhizal inoculation occurred when a low concentration of phosphorus fertilization was used (5 mg kg −1 ), but not at a higher concentration (10 mg kg −1 ). Treatments where Laccaria was applied, alone or in combination with a VAM fungus, resulted in the largest growth increases relative to nonmycorrhizal plants. An Acaulospora isolate was the most effective VAM fungus for E. urophylla although other VAM fungi also increased growth relative to the control. VAM fungi alone had little effect on E. globulus growth, but plants inoculated with both ECM and VAM fungi were larger than plants only inoculated with Laccaria . ECM fungi had a major impact on root system form by reducing the proportion of fine roots (specific root length) relative to nonmycorrhizal plants or those with VAM. Changes in root colonization patterns over time must be considered in studies of ECM/VAM interactions, as there can be substantial changes in the relative importance of the two associations.

Climate is widely recognised as an important determinant of the latitudinal diversity gradient. However, most existing studies make no distinction between direct and indirect effects of climate, which substantially hinders our understanding of how climate constrains biodiversity globally. Using data from 35 large forest plots, we test hypothesised relationships amongst climate, topography, forest structural attributes (stem abundance, tree size variation and stand basal area) and tree species richness to better understand drivers of latitudinal tree diversity patterns. Climate influences tree richness both directly, with more species in warm, moist, aseasonal climates and indirectly, with more species at higher stem abundance. These results imply direct limitation of species diversity by climatic stress and more rapid (co-)evolution and narrower niche partitioning in warm climates. They also support the idea that increased numbers of individuals associated with high primary productivity are partitioned to support a greater number of species.

Phosphorus (P) is a necessary nutrient for plant growth and plays an important role in plant metabolisms; however, the majority of P in soil is in insoluble forms. Phosphate solubilizing bacteria (PSB) can convert the insoluble phosphates into plant-available forms and may have the potential for use in sustainable agricultural practices. This study examined the effects of two native PSB, namely Bacillus aryabhattai (JX285) and Pseudomonas auricularis (HN038), and a mixture of both strains (1:1) on the growth of Camellia oleifera Abel. seedlings. The results showed a significant promotion of the growth of C. oleifera plants by three inoculation treatments. All the PSB inoculation treatments could improve the leaf nitrogen (N) and P content and had positive effects on the available N, P, and potassium (K) content of rhizosphere soil. A co-inoculation of the two native PSB strains caused a synergistic effect and achieved the best benefit. In conclusion, B. aryabhattai and P. auricularis could be used as biological agents instead of chemical fertilizers for agricultural production to reduce environmental pollution and increase the yield of tea oil.

Climate models simulating continental scale deforestation suggest a warming effect of land clearing on the surface air temperature in the tropical zone and a cooling effect in the boreal zone due to different control of biogeochemical and biophysical processes. Ongoing land-use/cover changes mostly occur at local scales (hectares), and it is not clear whether the local-scale deforestation will generate temperature patterns consistent with the climate model results. Here we paired 40 and 12 flux sites with nearby weather stations in North and South America and in Eastern Asia, respectively, and quantified the temperature difference between these paired sites. Our goal was to investigate the response of the surface air temperature to local-scale (hectares) land clearing across latitudes using the surface weather stations as proxies for localized land clearing. The results show that north of 10°N, the annual mean temperature difference (open land minus forest) decreases with increasing latitude, but the temperature difference shrinks with latitude at a faster rate in the Americas [−0.079 (±0.010) °C per degree] than in Asia [−0.046 (±0.011) °C per degree]. Regression of the combined data suggests a transitional latitude of about 35.5°N that demarks deforestation warming to the south and cooling to the north. The warming in latitudes south of 35°N is associated with increase in the daily maximum temperature, with little change in the daily minimum temperature while the reverse is true in the boreal latitudes.

Diversity and community composition of soil microorganisms along the elevation climosequences have been widely studied, while the microbial metabolic potential, particularly in regard to carbon (C) cycling, remains unclear. Here, a metagenomic analysis of C related genes along five elevations ranging from 767 to 4190 m at Mount Kilimanjaro was analysed to evaluate the microbial organic C transformation capacities in various ecosystems. The highest gene abundances for decomposition of moderate mineralizable compounds, i.e. carbohydrate esters, chitin and pectin were found at the mid-elevations with hump-shaped pattern, where the genes for decompositions of recalcitrant C (i.e. lignin) and easily mineralizable C (i.e. starch) showed the opposite trend (i.e. U-shaped pattern), due to high soil pH and seasonality in both low and high elevations. Notably, the gene abundances for the decompositions of starch, carbohydrate esters, chitin and lignin had positive relationships with corresponding C compounds, indicating the consistent responses of microbial functional profiles and metabolites to elevation climosequences. Understanding of adaptation of microbial communities, potential function and metabolites to elevation climosequences and their influencing factors provided a new insight for the regulation of terrestrial C storage.

BACKGROUND: Dalbergia odorifera is an economically and culturally important species in the Fabaceae because of the high-quality lumber and traditional Chinese medicines made from this plant, however, overexploitation has increased the scarcity of D. odorifera. Given the rarity and the multiple uses of this species, it is important to expand the genomic resources for utilizing in applications such as tracking illegal logging, determining effective population size of wild stands, delineating pedigrees in marker assisted breeding programs, and resolving gene networks in functional genomics studies. Even the nuclear and chloroplast genomes have been published for D. odorifera, the complete mitochondrial genome has not been assembled or assessed for sequence transfer to other genomic compartments until now. Such work is essential in understanding structural and functional genome evolution in a lineage (Fabaceae) with frequent intergenomic sequence transfers. RESULTS: We integrated Illumina short-reads and PacBio CLR long-reads to assemble and annotate the complete mitochondrial genome of D. odorifera. The mitochondrial genome was organized as a single circular structure of 435 Kb in length containing 33 protein coding genes, 4 rRNA and 17 tRNA genes. Nearly 4.0% (17,386 bp) of the genome was annotated as repetitive DNA. From the sequence transfer analysis, it was found that 114 Kb of DNA originating from the mitochondrial genome has been transferred to the nuclear genome, with most of the transfer events having taken place relatively recently. The high frequency of sequence transfers from the mitochondria to the nuclear genome was similar to that of sequence transfer from the chloroplast to the nuclear genome. CONCLUSION: For the first-time, the complete mitochondrial genome of D. odorifera was assembled in this study, which will provide a baseline resource in understanding genomic evolution in the highly specious Fabaceae. In particular, the assessment of intergenomic sequence transfer suggests that transfers have been common and recent indicating a possible role in environmental adaptation as has been found in other lineages. The high turnover rate of genomic colinearly and large differences in mitochondrial genome size found in the comparative analyses herein providing evidence for the rapid evolution of mitochondrial genome structure compared to chloroplasts in Faboideae. While phylogenetic analyses using functional genes indicate that mitochondrial genes are very slowly evolving compared to chloroplast genes.

Pterocarpus is a genus of trees mainly distributed in tropical Asia, Africa, and South America. Some species of Pterocarpus are rosewood tree species, having important economic value for timber, and for some species, medicinal value as well. Up to now, information about this genus with regard to the genomic characteristics of the chloroplasts has been limited. Based on a combination of next-generation sequencing (Illumina Hiseq) and long-read sequencing (PacBio), the whole chloroplast genomes (cp genomes) of five species (rosewoods) in Pterocarpus (Pterocarpus macrocarpus, P. santalinus, P. indicus, P. pedatus, P. marsupium) have been assembled. The cp genomes of five species in Pterocarpus have similar structural characteristics, gene content, and sequence to other flowering plants. The cp genomes have a typical four-part structure, containing 110 unique genes (77 protein coding genes, 4 rRNAs, 29 tRNAs). Through comparative genomic analysis, abundant simple sequence repeat (SSR)loci (333–349) were detected in Pterocarpus, among which A /T single nucleotide repeats accounted for the highest proportion (72.8–76.4%). In the five cp genomes of Pterocarpus, eight hypervariable regions, including trnH-GUG_psbA, trnS-UGA_psbC, accD-psaI, ndhI-exon2_ndhI-exon1, ndhG_ndhi-exon2, rpoC2-exon2, ccsA, and trnfM-CAU, are proposed for use as DNA barcode regions. In the comparison of gene selection pressures (P. santalinus as the reference genome), purifying selection was inferred as the primary mode of selection in maintaining important biological functions. Phylogenetic analysis shows that Pterocarpus is a monophyletic group. The species P. tinctorius is resolved as early diverging in the genus. Pterocarpus was resolved as sister to the genus Tipuana.

BACKGROUND: MYB transcription factors are a kind of DNA binding protein that can specifically interact with the promoter region. Members of MYB TFs are widely involved in plant growth and development, secondary metabolism, stress response, and hormone signal transduction. However, there is no report of comprehensive bioinformatics analysis on the MYB family of Casuarina equisetifolia. RESULTS: In this study, bioinformatics methods were used to screen out 182 MYB transcription factors from the Casuarina equisetifolia genome database, including 69 1R-MYB, 107 R2R3-MYB, 4 R1R2R3-MYB, and 2 4R-MYB. The C. equisetifolia R2R3-MYB genes were divided into 29 groups based on the phylogenetic topology and the classification of the MYB superfamily in Arabidopsis thaliana, while the remaining MYB genes (1R-MYB, R1R2R3-MYB, and 4R-MYB) was divided into 19 groups. Moreover, the conserved motif and gene structure analysis shown that the members of the CeqMYBs were divided into the same subgroups with mostly similar gene structures. In addition, many conserved amino acids in the R2 and R3 domains of CeqMYBs by WebLogo analysis, especially tryptophan residues (W), with 3 conserved W in R2 repeat and 2 conserved W in R3 repeat. Combining promoter and GO annotation analysis, speculated on the various biological functions of CeqMYBs, thus 32 MYB genes were selected to further explore its response to salt stress by using qPCR analysis technique. Most CeqMYB genes were differentially regulated following multiple salt treatments. CONCLUSIONS: Seven genes (CeqMYB164, CeqMYB4, CeqMYB53, CeqMYB32, CeqMYB114, CeqMYB71 and CeqMYB177) were assigned to the "response to salt stress" by GO annotation. Among them, the expression level of CeqMYB4 was up-regulated under various salt treatments, indicating CeqMYB4 might participated in the response to salt stress. Our results provide important information for the biological function of C. equisetifolia, as well as offer candidate genes for further study of salt stress mechanism.

Forest dynamics plots, which now span longitudes, latitudes, and habitat types across the globe, offer unparalleled insights into the ecological and evolutionary processes that determine how species are assembled into communities. Understanding phylogenetic relationships among species in a community has become an important component of assessing assembly processes. However, the application of evolutionary information to questions in community ecology has been limited in large part by the lack of accurate estimates of phylogenetic relationships among individual species found within communities, and is particularly limiting in comparisons between communities. Therefore, streamlining and maximizing the information content of these community phylogenies is a priority. To test the viability and advantage of a multi-community phylogeny, we constructed a multi-plot mega-phylogeny of 1347 species of trees across 15 forest dynamics plots in the ForestGEO network using DNA barcode sequence data (rbcL, matK, and psbA-trnH) and compared community phylogenies for each individual plot with respect to support for topology and branch lengths, which affect evolutionary inference of community processes. The levels of taxonomic differentiation across the phylogeny were examined by quantifying the frequency of resolved nodes throughout. In addition, three phylogenetic distance (PD) metrics that are commonly used to infer assembly processes were estimated for each plot [PD, Mean Phylogenetic Distance (MPD), and Mean Nearest Taxon Distance (MNTD)]. Lastly, we examine the partitioning of phylogenetic diversity among community plots through quantification of inter-community MPD and MNTD. Overall, evolutionary relationships were highly resolved across the DNA barcode-based mega-phylogeny, and phylogenetic resolution for each community plot was improved when estimated within the context of the mega-phylogeny. Likewise, when compared with phylogenies for individual plots, estimates of phylogenetic diversity in the mega-phylogeny were more consistent, thereby removing a potential source of bias at the plot-level, and demonstrating the value of assessing phylogenetic relationships simultaneously within a mega-phylogeny. An unexpected result of the comparisons among plots based on the mega-phylogeny was that the communities in the ForestGEO plots in general appear to be assemblages of more closely related species than expected by chance, and that differentiation among communities is very low, suggesting deep floristic connections among communities and new avenues for future analyses in community ecology.

Climate change could alter plant aboveground and belowground resource allocation. Compared with shoots, we know much less about how roots, especially root system architecture (RSA) and their interactions, may respond to temperature changes. Such responses could have great influence on species’ acquisition of resources and their competition with neighbors. We used a gel-based transparent growth system to in situ observe the responses of RSA and root interactions of three common subtropical plant species seedlings in Asia differing in growth forms (herb, shrub and tree) under a wide growth temperature range of 18-34°C, including low and supraoptimal temperatures. Results showed that the RSA, especially root depth and root width, of the three species varied significantly in response to increased temperature although the response of their aboveground shoot traits was very similar. Increased temperature was also observed to have little impact on shoot/root resource allocation pattern. The variations in RSA responses among species could lead to both the intensity and direction change of root interactions. Under high temperature, negative root interactions could be intensified and species with larger root size and fast early root expansion had competitive advantages. In summary, our findings indicate that greater root resilience play a key role in plant adapting to high temperature. The varied intensity and direction of root interactions suggest changed temperatures could alter plant competition. Seedlings with larger root size and fast early root expansion may better adapt to warmer climates.

, belongs to Pyronemataceae (Pezizomycetes, Pezizales).