Diversity and Distributions, (Diversity Distrib.) (2012) 18, 754–768
BIODIVERSITY Modelling chestnut biogeography for RESEARCH American chestnut restoration Songlin Fei1*, Liang Liang2, Frederick L. Paillet3, Kim C. Steiner4, Jingyun Fang5, Zehao Shen5, Zhiheng Wang6 and Frederick V. Hebard7
1Department of Forestry and Natural ABSTRACT Resources, Purdue University, West Lafayette, Aim Chestnuts (Castanea spp.) are ecologically and economically important IN 47907, USA, 2Department of Geography, University of Kentucky, Lexington, KY 40506, species. We studied the general biology, distribution and climatic limits of seven USA, 3Department of Geosciences, University chestnut species from around the world. We provided climatic matching of of Arkansas, Fayetteville, AR 72701, USA, Asiatic species to North America to assist the range-wide restoration of American 4School of Forest Resources, Pennsylvania chestnut [C. dentata (Marsh.) Borkh.] by incorporating blight-resistant genes State University, University Park, PA 16802, from Asiatic species. 5 USA, Department of Ecology, Peking Location North America, Europe and East Asia. University, Beijing 100871, China, 6Department of Biology, University of Methods General chestnut biology was reviewed on the basis of published DK-2100 Copenhagen, Copenhagen, literature and field observations. Chestnut distributions were established using A Journal of Conservation Biogeography Denmark, 7The American Chestnut published range maps and literature. Climatic constraints were analysed for the Foundation, 29010 Hawthorne Drive, northern and southern distribution limits and the entire range for each species Meadowview, VA 24361, USA using principal component analysis (PCA) of fourteen bioclimatic variables. Climatic envelope matching was performed for three Chinese species using Maxent modelling to predict corresponding suitable climate zones for those species in North America.
Results Chestnuts are primarily distributed in the warm-temperate and subtropical zones in the northern hemisphere. PCA results revealed that thermal gradient was the primary control of chestnut distribution. Climatic spaces of different species overlap with one another to different degrees, but strong similarities are shown especially between Chinese species and American species. Climatic envelope matching suggested that large areas in eastern North America have a favourable climate for Chinese species.
Main conclusions The general biological traits and climatic limits of the seven chestnut species are very similar. The predictions of Chinese species climatic range corresponded with most of the historical American chestnut range. Thus, a regionally adapted, blight-resistant, introgressed hybrid American chestnut appears feasible if a sufficiently diverse array of Chinese chestnut germplasm is used as a source of blight resistance. Our study provided a between-continent climate matching approach to facilitate the range-wide species restoration, which can be readily applied in planning the restoration of other threatened or endangered species. and Distributions *Correspondence: Songlin Fei, Department of Forestry and Natural Resources, Purdue Keywords University, 715 West State Street, West Lafayette, IN 47907, USA. American chestnut, climatic envelope matching, conservation biogeography, E-mail: [email protected] Fagaceae, restoration, species distribution modelling.
endangered or locally extinct species (Seddon et al., 2007; INTRODUCTION Richardson et al., 2009). Successful zoological stories include
Diversity Managed reintroduction and relocation have been recognized some charismatic species such as the gray wolf (Canis lupus L.) as viable options to assist the restoration and recovery of in the western Great Lakes region in North America
DOI:10.1111/j.1472-4642.2012.00886.x 754 http://wileyonlinelibrary.com/journal/ddi ª 2012 Blackwell Publishing Ltd Chestnut biogeography
(Mladenoff et al., 1995) and golden lion tamarins (Leontopi- To begin the process of restoring American chestnut to its thecus rosalia L.) in Brazil (Kleiman & Mallinson, 1998). While former habitats, efforts are currently underway to introduce reintroduction and relocation have also been applied to plants, blight-resistance alleles by hybridizing the species with most efforts have focused on population, community and C. mollissima and recovering the American type by back- ecosystem levels (Young, 2000). Range-wide species-level- crossing and recurrent selection (Hebard, 2006). Because the managed reintroductions of plants are rare, and literature on blight fungus is native to Asia, all Asiatic chestnut species the special problems of restoring plant species to their former carry natural levels of resistance, although the number of habitat is scant (Falk et al., 1996). As pointed out by Seddon genes involved, whether they differ among species, and how et al. (2007), successful reintroduction often involves the use of they are distributed among linkage groups is largely habitat simulation modelling and requires a well-designed unknown. It is possible that the other Asiatic species reintroduction plan. In this study, we demonstrate the use of (C. crenata, C. seguinii and C. henryi) carry resistance genes ordination and climate matching techniques to understand the that are not present in C. mollissima. If that is the case, biogeography of chestnut (Castanea, Fagaceae) to assist in the intercrosses with those species would be highly useful in range-wide ecological restoration of American chestnut backstopping and potentially strengthening the resistance [C. dentata (Marsh.) Borkh.] through the incorporation of levels achieved in the initial stages of breeding using regionally adapted and blight-resistant genetic materials from C. mollissima alone. Although the backcross and recurrent Chinese chestnut species. selection process is intended to dilute the Asiatic portion of Castanea is a genus of forest trees with rather exceptional the hybrid genome in favour of the American portion (except ecological, economic and even cultural importance, primarily for resistance alleles, of course, which are retained through because species of this genus regularly bear abundant mast strong selection pressure), some non-American traits will crops of sweet nuts. The chestnut genus includes seven remain depending upon how tightly linked they are to species with their native ranges widely distributed in the resistance loci. Furthermore, it is unknown to what extent northern hemisphere. Four species occur in Asia [C. molliss- resistance alleles may influence other fitness characteristics ima Bl., C. henryi (Skan) Rehder & Wilson, and C. seguinii (and vice versa) through pleiotropy and epistasis. Unlike Dode. in China and C. crenata Sieb. & Zucc. in Japan] (Wu other highly bred plants that are intended for propagation & Raven, 1999), and one species occurs in Europe (C. sativa under cultivation, a version of American chestnut will Mill.) (Conedera et al., 2004). At least two native chestnut ultimately be restored to the wild state to face decades or species are found in North America (C. dentata and C. pu- centuries of unmitigated competition in the natural forest mila Mill.) (Johnson, 1988). Some studies have separated environment. For these reasons, it is very important that the C. pumila into eight or more poorly defined taxa according species donating resistance alleles be reasonably well adapted to growth form, leaf morphology, bur characteristics, habitat to the environment of American chestnut lest some element and blight susceptibility (Payne et al., 1991). In this study, of unfitness be inadvertently retained through the breeding only C. pumila sensu stricto was included in the analysis. The process. A particular concern for potential blight-resistant biology and ecology of chestnuts were well studied for species backcross trees is whether they possess sufficient biogeo- such as C. mollissima and C. sativa (e.g. Zhang & He, 1999; graphical adaptability to enable restoration throughout the Conedera et al., 2004), but for other species in the genus, historical American chestnut range. Hence, a climatic match- fragmented and sometime anecdotal information was avail- ing of Asiatic species to North America will provide able. quantitative assessment of the potential climatic fitness of In North America, American chestnut was historically one of the backcross trees for restoring American chestnut forests. the most ecologically and economically important trees in the Unlike many reintroduction projects that attempt to eastern forests (Lutts, 2004; Davis, 2006), a dominant timber re-establish species within their historical ranges through the species over much of its range (Paillet, 2002) and the most release of wild or captive-bred individuals following extirpation common tree in portions of its range such as Pennsylvania or extinction in the wild (IUCN, 1998), the restoration of (Illick, 1919). Chestnut-dominated forests covered 80 million American chestnut involves the incorporation of blight- ha of the land from Maine all the way to Mississippi resistant genetic materials from Asiatic species of the same (MacDonald, 1978). During the early 20th century, the non- genus. Therefore, unlike most current reintroduction research native chestnut blight fungus [Cryphonectria parasitica (Murr.) that gains its inference mainly from post hoc interpretation of Barr, formerly Endothia parasitica] completely eliminated monitoring results (Seddon et al., 2007), a spatially explicit American chestnut trees from the overstorey (McCormick & climatic modelling approach is needed to plan for the Platt, 1980), and by the early 1970s not a hectare of the original reintroduction of blight-resistant American chestnut. forest remained blight-free (Hepting, 1974). All that remains of The geographic distribution of a species is a function of this once prominent species are root-collar sprouts that grow biotic factors (limits of genetic and physiological adaptation, into small trees between recurrent cycles of blight infection. competitive interactions, etc.), abiotic factors and movement The loss of this historically dominant and important forest factors (movement and dispersal opportunities and events) species is one of the most important events in the history of the (Peterson, 2011). Climatic variables, which constitute a major eastern North American forests. portion of abiotic factors, are first-order constraints of
Diversity and Distributions, 18, 754–768, ª 2012 Blackwell Publishing Ltd 755 S. Fei et al. species distribution, as evidenced in many range equilibrium Climatic parameters studies (e.g. Geber & Eckhart, 2005; Griffith & Watson, Studies have demonstrated that regional-level distributions of 2006). The three continental regions where chestnut species plant species are highly correlated with climatic factors (e.g. naturally occur, North America, southern and central Europe, Fang & Lechowicz, 2006; Rehfeldt et al., 2006). In this study, and East Asia, are affected by different air masses and have we adopted a suite of fourteen climatic and bioclimatic various topographic settings (Fang & Lechowicz, 2006). It is variables that likely affect the distribution of chestnut species, yet to be understood how the resultant climatic gradients in akin to an approach used in delineating climatic limits of these regions influence the distribution of chestnut species. world-wide beeches (Fagus spp.) (Fang & Lechowicz, 2006). What are the similarities and dissimilarities of climate These variables offer a comprehensive characterization of variables in forming different distribution patterns for the climatic factors that impose potential physiological constraints various chestnut species? Answering these questions will not on tree growth and regeneration. Monthly mean temperature, only fill the knowledge gap in chestnut biogeography, but monthly precipitation and geographic location data are major also provide important references for species conservation inputs for deriving these biologically related variables. For and restoration. convenience in the analysis, these variables are grouped into In this study, we investigated the general biogeography of the following three categories: thermal variables, moisture seven chestnut species, reconstructed their respective geo- variables and variability indices. graphic ranges and analysed the corresponding climatic constraints. Through bioclimatic analysis and simulation, we specifically intended to provide a comprehensive study of Thermal variables chestnut species distribution and their association with biocli- Temperature plays a fundamental role in controlling the matic factors and to further predict the climatic suitability of natural range of a tree species (Woodward, 1987). In Chinese species in North America to facilitate the efforts in particular, growing season warmth and winter coldness are American chestnut restoration. In addition, we will show that two direct limiting factors for plant growth (Fang & Lech- the bioclimatic limits of the modern genus Castanea are owicz, 2006). Here, the accumulation of growing season consistent with our knowledge of the origins and biogeography warmth is represented with the Kira’s warmth index (Kira, of the genus in the early Cenozoic and later, indicating that 1991) and Holdridge’s annual biotemperature (Holdridge, although the genus has migrated extensively, its bioclimatic 1947): envelope appears to have remained relatively unchanged for X the past tens of millions of years. WI ¼ ðT 5Þðfor months in which T>5 C in units of degree monthsÞ METHODS P Chestnut biology and distribution T ABT ¼ ðfor months in which T>0 C; We reviewed the general biology and biogeography of all 12 in units of degree monthsÞ seven chestnut species based on published literature and field observations. Spatial distributions of the two American where WI is the warmth index, ABT is the annual biotemper- species were based on Little’s (1977) maps. The distribution ature, and T is the monthly mean temperature over the year. of C. sativa was based on the natural and naturalized range Winter coldness, which mainly affects the northern boundaries composed by members of the EUFORGEN Noble Hardwoods of the distribution of tree species, is described with the mean Network (Ferna´ndez-Lo´pez & Alı´a, 2003). The distribution of temperature for the coldest month (MTCM) and Kira’s C. crenata was digitized from a map in Atlas of the Japanese coldness index (Kira, 1991): Flora (Horikawa, 1972). Distribution maps of the three X Chinese species were taken from the Database of China CI ¼ ðT 5Þðfor months in which T<5 C; Woody Plants (http://www.ecology.pku.edu.cn/plants/woody/ in units of degree monthsÞ index.asp), which contains the distribution of all 11,405 woody plants in China. In the database, these county-level where CI is the coldness index. In addition, annual mean resolution distribution maps were compiled using informa- temperature (AMT) and mean temperature for the warmest tion from all published flora records in China, including month (MTWM) were included as general thermal variables. Flora of China (Wu & Raven, 1999), additional national-level and provincial floras, as well as a great number of regional Moisture variables flora and local checklists of woody species. Distributions of the American species are putative pre-Columbian distribu- Moisture availability, another primary limiting factor for tions (Little, 1977). Modern distributions of European and species distribution, was also represented with several variables. Asian species are presumed to have been influenced by Annual precipitation (AP) is a first overall measure of moisture cultivation. availability. In addition, potential evapotranspiration (PET),
756 Diversity and Distributions, 18, 754–768, ª 2012 Blackwell Publishing Ltd Chestnut biogeography which is primarily a function of temperature, was produced to http://www.worldclim.org/). In this dataset, monthly temper- quantify potential moisture loss to the atmosphere. PET was ature and precipitation data were compiled and interpolated estimated using the Thornthwaite (1948) method: for all land regions across the world (Hijmans et al., 2005). The X resultant continuous climate data layers were available in raster = a PET ¼ 16 ð10 T IhÞ ðfor months in which T>0 CÞ format with a spatial resolution as fine as 30 arc-seconds (c. 1 km). The dataset also provides a suite of additional where Ih is the monthly heat index and a is an empirical derived variables, including AMT, MTWM, MTCM, ART, AP, exponent based on Ih as defined below X PWM and PDM. We calculated the additional bioclimatic 1:514 Ih ¼ ðÞT=5 ðfor months in which T>0 CÞ variables (WI, CI, ABT, K, PET, Im and EQ) needed for the analysis. The climatic information was extracted for the southern and 7 3 5 2 a ¼ 6:75 10 Ih 7:71 10 Ih þ 1:79 northern limits, as well as the entire distribution range for each 2 chestnut species. In geoprocessing, range polygons of chestnut 10 Ih þ 0:49 distribution were superimposed on climatic data layers. We With both precipitation and PET data available, one addi- manually placed 100 sample points along each range boundary tional variable, moisture index (Im), can be calculated by (either south or north). Climatic data were then extracted for accounting for the surplus and deficit conditions in water the sample points representing limiting climatic conditions on budget: the range boundaries of chestnut species. In addition, a X comprehensive sampling was conducted within the entire Surplus ¼ ðÞðP PET for months in which P PETÞ distribution range for each chestnut species. For sampling an entire range, a systematic grid scheme with an interval of 15 X arc-min between points along both latitudinal and longitudinal Deficit ¼ ðÞðP PET for months in which P