native woody species often demonstrate inherited latitudinal variation in the timing of growth cessation and cold hardiness in response to the latitudinal gradient in temperature ( pauley and perry 1954 ; howe et al .
can an introduced species develop such a pattern in a century or two ( weber and schmid 1998 ) ?
development of clinal variation could be accelerated by multiple introductions of populations from different latitudes ( novak and mack 2001 ) , by hybridization between closely related species from different climates ( hurka et al .
2003 ) , or by epigenetic inheritance of cold hardiness ( saxe et al . 2001 )
. this question can be addressed by studies comparing clinal variation in native and introduced species .
the diploid genus tamarix was introduced to north america in the mid-1800s to control erosion and to serve as a drought - tolerant ornamental flowering shrub or tree ( robinson 1965 ) .
although there are now several tamarix species in the united states ( usa ) , tamarix chinensis lour . and tamarix ramosissima ledeb .
tamarix chinensis is native to china , korea , and japan at latitude 2342n and longitude 79110e , whereas t. ramosissima occurs more widely across temperate asia at latitude 3053n and longitude 42127e ( baum 1978 ) .
although these two species are morphologically distinct in asia ( baum 1978 ) , north american specimens can not always be readily distinguished ( crins 1989 ) .
hereafter , we refer to the complex of t. ramosissima , t. chinensis , and their hybrids as saltcedar . in a study of a single nuclear locus in saltcedar , gaskin and schaal ( 2002 )
found that the most common genotype in north america is a heterozygote containing one allele that in asia is found only in t. ramosissima and another allele that in asia is found only in t. chinensis .
in other words , these north american heterozygotes have a hybrid genotype that has not been found in asia .
although homozygotes for either of the two alleles are also common in north america , it is not possible using a single locus to determine whether they are true t. chinensis and t. ramosissima or cryptic hybrids .
for this reason , there is a strong need for analysis of additional loci to characterize the taxonomic status of north american saltcedar .
saltcedar is now the second most abundant riparian woody plant in the interior western usa ( friedman et al .
2005 ) , occupying 470 000650 000 ha ( zavaleta 2000 ) . because replacement of native cottonwood - willow ( populus - salix ) communities by saltcedar can be associated with degradation of habitat for vertebrates and increased water loss from evapotranspiration ( shafroth et al .
2005 ) , there is an intense interest in controlling its spread ( hart et al .
2005 ) . to inform such efforts , it is important to know what factors limit saltcedar distribution and to understand how that distribution would be altered by changing temperatures .
although the cold hardiness of saltcedar is unknown , there is evidence that cold temperatures may limit its expansion in north america . in the us north of about 39n latitude
, saltcedar is relatively scarce , individuals are small , and dieback is frequent ( lesica and miles 2001 ; friedman et al .
on the other hand , some populations of saltcedar in montana are growing rapidly ( pearce and smith 2003 ; sexton et al .
furthermore , growth chamber experiments comparing northern and southern saltcedar have demonstrated inherited temperature - dependent differences in root 
genotypic differences among t. chinensis , t. ramosissima , and their hybrids could cause variation in cold hardiness across north america .
if a brief residence in north america has limited development of clinal variation in saltcedar , we might expect a stronger cline in a comparable native species . the native riparian tree plains cottonwood [ populus deltoides marshall subsp .
monilifera ( aiton ) eckenwalder ] is an important competitor of saltcedar ranging from northern texas to southern manitoba , saskatchewan , and alberta at latitude 3055n and longitude 96114w ( van haverbeke 1990 ; friedman et al .
populus deltoides is known to be cold hardy ; buds and stems collected in winter can survive chilling to 80c ( sakai and weiser 1973 ) .
common - garden studies of cottonwood including p. deltoides have shown inherited latitudinal variation in the timing of fall growth cessation ( howe et al .
when grown in massachusetts , p. deltoides from 43 to 46n ceased growing an average of 1 month earlier than p. deltoides from 30 to 32n ( pauley and perry 1954 ) .
northern cottonwood , therefore , may develop cold resistance earlier in the fall than southern cottonwood .
we used a common - garden approach to compare cold hardiness in plains cottonwood and saltcedar collected along a latitudinal gradient in the central usa .
we hypothesized that plains cottonwood is more cold tolerant than saltcedar , and that the physiological limits imposed by cold temperatures are consistent with the northern limit to saltcedar distribution in the usa .
we expected to find strong latitudinal variation in cold hardiness in plains cottonwood , but not in saltcedar , because there has been little time for evolution of such variation in this introduced species .
finally , we used microsatellite analysis in plains cottonwood and saltcedar to investigate genetic variation along the latitudinal gradient . in saltcedar , it was necessary to develop the microsatellite primers to perform the analysis ( gaskin et al .
this work is especially important in saltcedar , a widespread dominant tree in north america whose taxonomic status is still uncertain ( gaskin and schaal 2002 ) .
we hoped to determine to what extent saltcedar in north america behaves as a single population as opposed to two or more isolated populations .
we collected plains cottonwood and saltcedar from 15 sites along a latitudinal gradient from 29 to 48n in the central usa ( table 1 ) . in february and march of 2005
, we collected samples from 25 ( or more ) distinct genetic individuals of each species at each site , except for the three southernmost sites where plains cottonwood was scarce or absent ( table 1 ) .
the sampled populations occurred as bands 2200 m wide along the shore of a river or reservoir .
we walked along the band and sampled the first 25 plants that had young , rapidly growing stems suitable for vegetative propagation .
we planted the cuttings in pots ( top diameter = 7 cm , depth = 25 cm ) filled with fafard f2-sv growing medium with 25% perlite and no added fertilizer or rooting hormone .
we rooted the cuttings under mist in a greenhouse in fort collins , co. on may 31 , we moved the plants to a shadehouse in fort collins at latitude 40.58n and longitude 105.14w .
the shadehouse is an outdoor space with no walls and a roof of widely spaced boards to reduce direct sun .
also shown are the number of individuals collected at each location and the number tested for cold hardiness on april 17 , 2006 we determined the change in cold hardiness over time by periodically testing plants from four species  latitude combinations : northern and southern plains cottonwood and saltcedar . to conduct a test ,
we selected a single individual from each of the four species  latitude combinations ; different individuals were selected for each test .
we performed the test 17 times from september 13 , 2005 to june 5 , 2006 .
for the first eight tests , ending with january 23 , 2006 , the selected species 
latitude combinations were plains cottonwood from latitudes 33.06 and 47.60n and saltcedar from latitudes 29.18 and 47.60n ; for the final nine tests , beginning with february 6 , 2006 , we used plains cottonwood and saltcedar from latitudes 34.90 and 45.43n .
the change in latitudes was necessary to provide sufficient stem material to perform the tests .
we measured cold hardiness of the selected individuals using freeze - induced electrolyte leakage ( burr et al .
this method is based on the fact that freezing of stem tissues causes leakage of electrolytes , which can be quantified by measurement of electrical conductivity . to perform a test , we collected a 36-cm piece of current - year stem from each of the four selected individuals .
we cut each piece into 1-cm samples and placed each sample in a vial ( or test tube ) with 1 ml of deionized water and a piece of lead shot to promote nucleation of ice crystals .
the 36 samples from each selected individual were divided into six sets of six tubes ( or vials ) .
we placed one set ( controls ) in a refrigerated water bath ( fall ) or an ice water bath ( spring ) .
we placed the remaining five sets in a programmable temperature chamber ( testequity model 115 , thousand oaks , ca ) , held them at 2c ( fall ) or 4c ( spring ) overnight to ensure that all were frozen , and then exposed them to gradually declining temperatures ( 5c per hour ) . at five different temperatures
, we removed a set of the samples and temporarily stored them at 4c until all were removed from the chamber .
the five temperatures varied seasonally and were chosen to be close to the warmest lethal temperatures for all of the selected individuals . after adding 10 ml ( vials ) or 4 ml ( tubes ) of deionized water , we agitated samples for 20 h on a shaker table , measured the conductivity of each vial or tube using an amber science , conductivity meter model 1056 ( amber science inc . ,
eugene , or ) , boiled samples for 20 min , agitated again for 20 h and remeasured the conductivity .
1967 ) for a frozen 1-cm sample was where it is index of injury resulting from exposure to temperature t , lt is specific conductance of leachate from sample frozen at temperature t , lk is specific conductance of leachate from sample frozen at temperature t and then heat - killed , l0 is specific conductance of leachate from unfrozen control sample , and ld is specific conductance of leachate from unfrozen sample , heat killed .
the index of injury is designed to correct for error introduced by variation between samples in the quantity of electrolytes available for leaching following injury ( flint et al .
for an individual plant on a given day , this parameter displayed a sigmoid ( s - shaped ) relationship with temperature .
it approached a lower asymptote of 0 at warm temperatures that did not harm the sample , and approached an upper asymptote between 0 and 1 at cold temperatures that killed the sample .
comparison to frozen stems held for 2 weeks at room temperature for direct observation of cold - induced mortality indicated that stem samples died when it reached a point midway between the low and high asymptotes .
for each of the four individuals selected in a test , we determined the killing temperature using a nonlinear sigmoid regression of the form where t is temperature ; a , b , and c are estimated parameters , such that a is the upper asymptote and ( ln(b))/c , the inflection point of the curve , is the killing temperature .
we estimated these intervals for the nonlinear model following the format for linear models where an inverse prediction interval for a new value of the independent variable ( here temperature , t ) takes the form of a prediction interval for a new observation of the dependent variable ( here index of injury , it ) divided by the slope estimate ( neter et al .
1996:167169 ) . because the nonlinear functional form prohibited this simple algebraic manipulation , we estimated 95% prediction intervals for a new observation of it in the nonlinear model for an interval of values of t and determined the smallest and largest values of t that had 95% prediction intervals for it that included the estimated killing temperature .
we measured the variation in cold hardiness over latitude by comparing all individuals with 35 cm of stem ( cottonwood , n = 187 ; saltcedar , n = 192 ) on april 17 , 2006 .
we cut a 35-cm stem piece from each individual , cut each piece into seven 5-cm samples and distributed the samples into zip - lock plastic bags . for each latitude , there were seven plastic bags , each containing one 5-cm sample from each of the selected individuals from that latitude .
we placed one bag from each latitude ( control ) in an ice water bath and the remaining bags in the temperature chamber at 20c .
we reduced the temperature to 3.3c in 15 min and held the chamber at that temperature for 14.5 h to ensure that all twigs were frozen .
we then lowered the temperature 5c per hour to 71.1c and held this temperature for 11.5 h. we removed one bag from each latitude at 8.5 , 20.3 , 29.4 , 41.9 , 55.4 , and 71.1c ( at the end of the 11.5-h hold ) .
we then placed samples upright in moist sand at room temperature for 2 weeks before determining survival by visual observation ( calkins and swanson 1990 ) .
live samples had moist bright green inner bark , and often grew new roots and leaves .
dead samples had brown  green inner bark , which was often dry , and grew no roots or leaves . for a given species at a given latitude
, there was a sigmoid relationship between temperature and percent survival similar to the sigmoid relationship observed in the conductivity trials .
this allowed us to use an analogous approach designed for categorical ( live or dead ) data , in which we defined the killing temperature as the temperature at which half of stem samples were killed .
for each species  latitude combination , we used proc probit in sas ( sas institute , inc . , cary , nc ) to calculate the killing temperature ( n = 323 individuals , depending on species and latitude , table 1 ) .
finally , for each species , we performed a multivariate probit analysis on all individuals combined using temperature and latitude as independent variables .
we assessed genetic similarity among plants using neighbor - joining analysis of microsatellite loci for each species .
extraction of dnas , pcr amplications , and fragment analysis were performed as in gaskin et al .
( 2006 ) , and for plains cottonwood we used the nine microsatellites : pmgc 61 , pmgc 223 , pmgc 667 , pmgc 2060 , pmgc 2088 , pmgc 2105 and pmgc 2573 ( from international populus genome consortium ssr resource 2007 , http://www.ornl.gov/sci/ipgc/ssr_resource.htm ) ; wpms14 ( smulders et al .
2001 ) ; and ptr2 ( dayanandan et al . 1998 ) . for each species
, we calculated chord distances ( cavalli - sforza and edwards 1967 ) among latitudes and used these to prepare an unrooted dendrogram .
the saltcedar dendrogram included 15 populations ( 369 individuals ) from north america and two populations ( 12 individuals ) from asia .
the 12 individuals collected in asia ( six t. ramosissima and six t. chinensis ) were never introduced to the shadehouse or used for cold - hardiness testing , but were included in this analysis to explore genetic relatedness between north american and asian saltcedar .
the plains cottonwood dendrogram included 13 populations ( all 304 cottonwood individuals introduced to the shadehouse ) .
all calculations were carried out using the programs seqboot , gendist , neighbor , and consense from the phylip ver .
felsenstein 2004 ) . analyses of molecular variance ( amova ) were carried out for both species using the distance method ( weir and cockerham 1984 ; excoffier et al . 1992 ; weir 1996 ) . as a dioecious species
saltcedar has perfect flowers and has been categorized as self - compatible ( brotherson and winkel 1986 ) , but stevens ( 1989 ) in arizona , found no seeds produced by flowers enclosed in mesh bags to exclude pollinators .
using some of our shadehouse individuals planted in a nearby garden , in fort collins , colorado , we also observed no seeds produced by bagged flowers even when we attempted to promote self - pollination using bees on a stick ( a. blair , colorado state university , personal communication , 2007 ) .
both plains cottonwood and saltcedar were more resistant to cold temperatures in the winter than in the summer , but this seasonal cycle was more extreme in cottonwood ( fig .
1 ) . although saltcedar was slightly more cold hardy in the early fall , cottonwood hardened off more rapidly and deeply ( fig .
in early september , northern and southern individuals of both species were killed by temperatures ranging from 4 to 19c . by mid october ,
both northern and southern cottonwood survived cooling to 70c , the coldest temperature attainable in our temperature chamber .
even when we extended the coldest treatment by holding the chamber at 70c for 13 h , we observed no cold damage to cottonwood during the winter ( november through february , data not shown ) .
in contrast , even in mid - winter saltcedar was killed at temperatures ranging from 33 to 47c .
fall ( a ) and spring ( b ) temporal variation in killing temperature for twigs of northern and southern plains cottonwood ( populus deltoides subsp .
monilifera ) and saltcedar ( tamarix ramosissima , tamarix chinensis , and hybrids ) grown in fort collins , colorado .
saltcedar ( sc ) data are solid lines , and cottonwood ( cw ) data are dashed lines .
different individuals were used on each sampling date . note that different latitudes were used in the spring and fall .
daily minimum temperature data are from the christman field weather station ( colorado climate center ) about 1 km from the shadehouse .
error bars are 95% confidence intervals for an individual sample . from september through january , there was latitudinal variation in cold hardiness of saltcedar ( fig .
during this period , the killing temperature for saltcedar from 29.18n latitude was 521c higher than that for saltcedar from 47.60n .
latitudinal variation in saltcedar was not apparent in the time sequence of cold hardiness measurements from february to may ( fig .
this may be because of weaker latitudinal variation in the spring or because the latitudes tested in the spring were not as extreme as those tested in the fall ( fig .
1 ) . although latitudinal variation was not apparent in the spring at this coarse scale ( fig .
1 ) , the detailed data from april 17 showed latitudinal variation for saltcedar ( fig .
2 ) , and multivariate probit analysis demonstrated that this latitudinal effect was significant ( p < 0.0001 , table 2 ) . killing temperature for twigs of plains cottonwood ( populus deltoides subsp .
monilifera ) and saltcedar ( tamarix ramosissima , tamarix chinensis , and hybrids ) as a function of latitude of origin on april 17 , 2006 .
for saltcedar at latitude 29.18n , there was insufficient data to calculate the killing temperature .
analysis of parameter estimates for multivariate probit analysis of the effects of temperature and latitude on survival of stems of plains cottonwood ( 187 individuals ) and saltcedar ( 192 individuals ) on april 17 , 2006 numbers of individuals by species and latitude are shown in table 1 .
there was strong latitudinal variation in cold hardiness of cottonwood during the fall and spring ( fig .
on september 26 and october 11 , killing temperatures for cottonwood from 33.06n were > 53 and > 43c higher than those for cottonwood from 47.60n . from november 11 to march 13 , however , cottonwood from all latitudes survived temperatures down to 70c , the coldest temperature attainable in our temperature chamber , which prevented observations of latitudinal variation during this period ( fig .
multivariate probit analysis of the detailed data set from plants collected on april 17 showed a significant effect ( p < 0.0001 ) of latitude of origin on the killing temperature for cottonwood ( fig .
overwinter survival of whole saltcedar plants in the shadehouse was mostly below 0.4 and strongly correlated with latitude of origin ( r = 0.58 , p = 0.0009 , fig .
3 ) , which is consistent with the hypotheses that frost damage is an important source of mortality in saltcedar , and that southern plants are more susceptible to cold than northern plants .
in contrast , overwinter survival of whole cottonwood plants in the shadehouse was mostly above 0.4 and uncorrelated to latitude ( p = 0.51 , fig .
3 ) , suggesting that frost damage was relatively unimportant in plains cottonwood in the winter of 20052006 .
heavy winter mortality of southern saltcedar is inconsistent with the fact that temperatures from september 2005 through may 2006 were always above the killing temperatures we measured for twigs of these species ( fig .
1 ) , suggesting either that saltcedar roots are more sensitive to cold than stems , that prolonged exposure to cold is more damaging than the brief exposures in our freezing trials , or that our application of freeze - induced electrolyte leakage underestimated killing temperatures for saltcedar .
monilifera ) and saltcedar ( tamarix ramosissima , tamarix chinensis , and hybrids ) in the shadehouse as a function of latitude of origin . survival is the proportion of plants alive on june 1 , 2006 .
survival was correlated to latitude in saltcedar ( p = 0.0009 , r = 0.58 ) , but not in plains cottonwood ( p = 0.51 , r = 0.04 ) .
microsatellite analysis demonstrated that neither saltcedar nor cottonwood have strong barriers to gene flow in the central usa .
fixation index ( fst ) was low for both species ( 0.049 for plains cottonwood and 0.075 for saltcedar ) , indicating that genetic variation within populations is large compared to variation between populations ( table 3 ) .
the dendrogram for saltcedar shows a gradual north  south genetic gradient with the eurasian parent populations as approximate end members ( fig .
northern saltcedar is most similar to t. ramosissima and southern saltcedar is most similar to t. chinensis .
4b ) . analysis of molecular variance ( amova ) calculated using the distance method for plains cottonwood ( populus deltoides subsp .
monilifera , 304 individuals ) and saltcedar ( tamarix ramosissima , tamarix chinensis , and hybrids , 369 individuals ) collected along a latitudinal gradient in the central usa unrooted dendrogram of chord distances ( cavalli - sforza and edwards 1967 ) inferred from neighbor - joining analysis of microsatellite loci data .
bootstrap values above 50% ( vertically oriented ) are shown to the left of the respective nodes , and were derived from a consensus of 1000 trees .
ultimate dendrogram branches are populations identified by latitude ( table 1 ) . ( a ) dendrogram of 15 us populations ( 369 plants ) of saltcedar ( tamarix chinensis , tamarix ramosissima , and hybrids ) based on nine microsatellite loci ( gaskin et al .
populations of six t. chinensis from china and six t. ramosissima from asia were also included in the analysis for comparison .
( b ) dendrogram of 13 us populations ( 304 plants ) of plains cottonwood ( populus deltoides subsp .
monilifera ) based on nine microsatellite loci ( international populus genome consortium ssr resource 2007 , http://www.ornl.gov/sci/ipgc/ssr_resource.htm , smulders et al .
all of our plains cottonwoods survived colder temperatures than are known to occur in the range of this species , suggesting that short - term exposure to extreme cold in mid - winter is not an important mortality factor .
other studies have reported extreme cold hardiness in populus ( sakai and weiser 1973 ) .
cold - related mortality is still possible in cottonwood when plants that are not completely hardened off encounter sudden extreme drops in temperature .
in fact , plains cottonwood may be more susceptible to cold damage than saltcedar in early fall and late spring ( fig .
1 ) . saltcedar , even when completely hardened off , did not survive below 33 to 47c , which is within the temperature range of the northern great plains .
in addition , there was heavy mortality in the shadehouse of southern saltcedar ( but not cottonwood ) during the winter of 20052006 ( fig .
whereas saltcedar grows to large size in the southern usa , northern stems are small in size and often older below ground than above ( lesica and miles 2001 ) , suggesting that winter die - back may be common in the north . finally , saltcedar is much more abundant in the south than in the north , and its occurrence is strongly correlated to mild winter temperatures ( friedman et al . 2005 ) .
we conclude that winter cold is a significant factor influencing distribution of saltcedar in the usa .
increases in winter low temperatures could promote northward spread of saltcedar in the future . we observed latitudinal variation in cold hardiness in both plains cottonwood and saltcedar .
this variation was expected for the native cottonwood , which is known to exhibit inherited latitudinal variation in the timing of growth cessation ( pauley and perry 1954 ; kaszkurewicz and fogg 1967 ; howe et al .
1995 , 1999 ; dunlap and stettler 1996 ) , but not for saltcedar , a woody species that has been present in north america for only about 150 years ( robinson 1965 ) .
such latitudinal variation in an introduced species may be a result of multiple introductions of genotypically and phenotypically distinct individuals or populations , or of rapid evolution after introduction ( maron et al .
the dominance in the north american invasion of hybrids ( t. ramosissima  t. chinensis ) that have not been found in eurasia ( gaskin and schaal 2002 ) , the absence of genetic isolation in the north american population ( fig .
4 ) , and the fact that north american tamarix are more similar to each other than to eurasian tamarix ( fig .
4 ) support the hypothesis that the observed latitudinal gradient in north america evolved after introduction .
hybridization involving t. ramosissima , t. chinensis and possibly other species may have accelerated this process by producing a population that was more heterogeneous with respect to cold hardiness than any of the originally introduced populations .
an alternative explanation for the observed latitudinal gradient in saltcedar cold hardiness is a carryover effect .
cuttings might have a chemical memory of the climate of the source plant that influences the phenology and cold hardiness of the clones grown in our colorado shadehouse ( van zandt and mopper 2002 ; heide 2003 ) .
if a carryover effect exists , it could be expected to fade over time as the plants are influenced by the colorado climate .
there is a need for multi - year common - garden studies of cold hardiness to investigate whether the observed latitudinal variation fades over time .
development of clinal variation in saltcedar cold hardiness could be accelerated by maternal effects ( galloway and etterson 2007 ) .
for example , in picea abies ( norway spruce ) temperature and photoperiod experienced by a tree during seed production can affect the cold hardiness of progeny ( saxe et al .
spread of introduced species across long environmental gradients is commonly attributed to phenotypic plasticity and multiple introductions ( neuffer and hurka 1999 ; novak and mack 2001 ) .
the results of this and other studies , however , suggest that invasion can also be facilitated by rapid evolution of clinal variation ( weber and schmid 1998 ; sexton et al . 2002 ; maron et al .
our microsatellite data show a lack of strong genetic isolation among tamarix populations in the central usa , and a gradual trend from individuals that resemble asian t. chinensis in the south to individuals that resemble asian t. ramosissima in the north .
this result argues against the possibility that the observed latitudinal variation in cold hardiness reflects multiple introductions of reproductively isolated taxa , and against the possibility that hybrid tamarix in north america are reproductively isolated from sympatric populations of either parent .
our results are consistent with the possibility of adaptive radiation from a single , mostly hybrid population , and with the possibility of hybridization and backcrossing following introduction of t. ramosissima in the north and t. chinensis in the south .
therefore , some of the latitudinal gradient in saltcedar cold hardiness may be a result of patterns of introduction as opposed to natural selection .
future genetic work should focus first on analyzing additional asian samples of t. chinensis and t. ramosissima and other tamarix species to pinpoint the locations of origin and genotypes of material brought to north america , second on developing additional genetic loci to further quantify the genetic structure of the saltcedar population in north america , and third on characterization of specific genes related to cold hardiness ( mckay and latta 2002 ) .
 to investigate the evolution of clinal variation in an invasive plant , we compared cold hardiness in the introduced saltcedar ( tamarix ramosissima , tamarix chinensis , and hybrids ) and the native plains cottonwood ( populus deltoides subsp . 
 monilifera ) . in a shadehouse in colorado ( 41n ) 
 , we grew plants collected along a latitudinal gradient in the central united states ( 2948n ) . on 17 occasions between september 2005 and june 2006 , we determined killing temperatures using freeze - induced electrolyte leakage and direct observation . in midwinter , cottonwood survived cooling to 70c , while saltcedar was killed at 33 to 47c . 
 frost sensitivity , therefore , may limit northward expansion of saltcedar in north america . 
 both species demonstrated inherited latitudinal variation in cold hardiness . 
 for example , from september through january killing temperatures for saltcedar from 29.18n were 521c higher than those for saltcedar from 47.60n , and on september 26 and october 11 , killing temperatures for cottonwood from 33.06n were > 43c higher than those for cottonwood from 47.60n . 
 analysis of nine microsatellite loci showed that southern saltcedars are more closely related to t. chinensis while northern plants are more closely related to t. ramosissima . 
 hybridization may have introduced the genetic variability necessary for rapid evolution of the cline in saltcedar cold hardiness . 
Introduction
Methods
Results
Discussion