TY - JOUR AU - Dawson, Todd E. AB - Abstract Water transport in conifers occurs through single-celled tracheids that are connected to one another via intertracheid pit membranes. These membranes have two components: the porous margo, which allows water to pass through the membrane, and the impermeable torus, which functions to isolate gas-filled tracheids. During drought, tracheids can become air filled and thus hydraulically dysfunctional, a result of air entering through the pit membrane and nucleating cavitation in the water column. What are the hydraulic tradeoffs associated with cavitation resistance at the pit level, and how do they vary within the structural components of the intertracheid pit? To address these questions, we examined pit structure in 15 species of Cupressaceae exhibiting a broad range of cavitation resistances. Across species, cavitation resistance was most closely correlated to the ratio of the torus to pit aperture diameter but did not vary systematically with margo porosity. Furthermore, our data indicate that constraints on pit hydraulic efficiency are shared: the pit aperture limits pit conductivity in more drought-resistant taxa, while increased margo resistance is more likely to control pit conductivity in species that are more vulnerable to cavitation. These results are coupled with additional data concerning pit membrane structure and function and are discussed in the context of the evolutionary biogeography of the Cupressaceae. Water transport in conifers occurs through narrow, single-celled conduits (tracheids) that are organized in overlapping, longitudinal files. This simple, homoxylous arrangement represents an ancestral vascular design that has remained remarkably consistent since its first appearance in the progymnosperms of the Mid-Devonian (Taylor et al., 2009). However, the small size of tracheids can impose a high resistance to water transport as compared with the large, hydraulically efficient vessels present in many angiosperms (Hacke et al., 2004; Sperry et al., 2006). Despite this handicap, conifer tracheids can be just as hydraulically efficient as angiosperm xylem for a given conduit diameter, a result that can be wholly attributed to the distinctive structure of the conifer intertracheid pit membrane (Pittermann et al., 2005; Sperry et al., 2006). Because pit membranes also function to limit the spread of air from one conduit to another (cavitation), the physiological consequences of the transport efficiency versus cavitation safety tradeoffs in conifer and angiosperm pit membranes have received considerable attention at the pit and xylem levels, whereby cavitation resistance in north temperate woody plants appears to come at the cost of hydraulic efficiency (Pittermann et al., 2006a, 2006b; Sperry et al., 2006; Choat et al., 2008; Domec et al., 2008; Jansen et al., 2009; Schoonmaker et al., 2010). Previous work has shown that the integrated vascular performance of plants is key to understanding species distributions (Sperry et al., 1994; Brodribb and Hill, 1999; Pockman and Sperry, 2000; Choat et al., 2007), and within this framework, pit membranes have the potential to act as the nexus of the cavitation safety versus transport efficiency compromise. Yet, despite our progress, we are just starting to learn how these tradeoffs play out at the level of the pit membrane, particularly in one as complex as that of conifers. Hence, the goals of this study were to determine whether selection has acted to optimize conifer pit membrane performance in a manner that reflects species cavitation resistance and habitat distribution as well as to examine the role, if any, of evolutionary lineage. Unlike the homogenous pit membrane of angiosperm vessels, the conifer pit membrane is composed of two distinct regions: a thickened, centrally located torus and a porous margo region that surrounds it (Fig. 1; for study species, see Table I; Hacke et al., 2004; Choat et al., 2008; Choat and Pittermann, 2009). When tracheids are water filled, the pit membrane is centrally located in the pit chamber and water moves from tracheid to tracheid through the margo. Should an air-seeding event (cavitation) occur, causing a tracheid to become air filled, (i.e. embolized), the negative xylem pressure in the water-filled tracheid will act on the air-water interface in the margo pores by deflecting the pit membrane in the direction of the functional tracheid, thereby appressing the torus against the pit aperture border (Bailey, 1913; Liese, 1965; Liese and Bauch, 1967; Petty, 1972). This valve action of the membrane can create an effective seal that prevents further spread of air in the xylem. Cavitation is thought to occur when the water potential of the water-filled tracheid becomes negative enough to dislodge the torus from its sealing position, allowing air to enter the conduit. Overall, the structure of the torus-margo pit membrane must optimize what at first glance appear to be conflicting functional requirements: on the one hand, cavitation resistance selects for a combination of large tori and small apertures, but on the other hand, hydraulic efficiency favors porous margos, large apertures, and small tori. Figure 1. Open in new tabDownload slide SEM images of intertracheid pit membranes belonging to nine Cupressaceae species (of 15) that represent the broad range of observed cavitation pressures. The opaque torus region of the membrane (T) is held centrally by the microfibrils of the margo (M). Visually, increased cavitation resistance appears to be associated with increased margo porosity, but quantitative estimates of margo resistance made on the most intact regions of the pit membranes (Fig. 8) revealed no differences among the species surveyed. Figure 1. Open in new tabDownload slide SEM images of intertracheid pit membranes belonging to nine Cupressaceae species (of 15) that represent the broad range of observed cavitation pressures. The opaque torus region of the membrane (T) is held centrally by the microfibrils of the margo (M). Visually, increased cavitation resistance appears to be associated with increased margo porosity, but quantitative estimates of margo resistance made on the most intact regions of the pit membranes (Fig. 8) revealed no differences among the species surveyed. Study species, figure abbreviations (Fig. Abbrevs.), locations (SFBG, San Francisco Botanical Garden, San Francisco; UCBG, University of California Botanical Garden, Berkeley, CA; UCSC, University of California, Santa Cruz, Arboretum, Santa Cruz, CA), and species natural history (Farjon, 2005) Table I. Study species, figure abbreviations (Fig. Abbrevs.), locations (SFBG, San Francisco Botanical Garden, San Francisco; UCBG, University of California Botanical Garden, Berkeley, CA; UCSC, University of California, Santa Cruz, Arboretum, Santa Cruz, CA), and species natural history (Farjon, 2005) Species Fig. Abbrevs. Location and Accession Phenology, Mature Tree Height, Native Elevational Range, and Habitat Athrotaxis laxifolia AL SFBG not cataloged Evergreen, 10–15 m, 1,000–1,200 m, montane forests, Tasmania Callitris rhomboidea CR SFBG #1999-0290 Evergreen, 10–15 m, 0–1,250 m, open woodland, Southeast Australia Calocedrus decurrens CD SFBG #XY-2004 Evergreen, 60 m, 50–2,960 m, mixed conifer forests, Oregon to Baja California Norte Cryptomeria japonica CJ SFBG not cataloged Evergreen, 50–60 m, 1–2,050 m, mixed evergreen forests, Japan Cupressus forbesii CF SFBG #1980-0055 Evergreen, 10 m, 210–1,400 m, chaparral, Baja California Norte Fitzroya cupressoides FC UCBG #2007.0165 Evergreen, 50–60 m, 4–1,000 m, emergent tree is evergreen rainforest, Chile Glyptostrobus pensilis GP UCBG #70.0169 Deciduous, 15–25 m, 1–730 m, river floodplains, deltas, Southern China Juniperus californica JC UCBG #83.0567 Evergreen, 7–10 m, 500–1,400 m, desert scrubland, Southern California to Baja California Norte Libocedrus plumosa LP UCSC #81.1172 Evergreen, 30–35 m, 1–600 m, lowland mixed angiosperm and conifer rainforests, New Zealand Metasequoia glyptostroboides MGS = stem MGR = root UCBG #49.0500 Deciduous, 35–50 m, 750 m, ravines and moist temperate forests, Central China Sequoiadendron giganteum SG UCBG #2002.1062 Evergreen, over 100 m, 1,400–2,150 m, western slopes of Sierra Nevada, California Sequoia sempervirens SS Campus, University of California, Berkeley, CA Evergreen, over 100 m, 1–750 m, moist and foggy climates, Central to Northern California coasts Taxodium distichum TD UCBG #60.1174 Evergreen, 90–95 m, 1,400–2,150 m, mixed conifer montane forests, Eastern California Taiwania cryptomerioides TC SFBG #1984.93, #1990.616 Evergreen, 60–65 m, 1,750–2,900 m, cool temperate forests, Asia Widdringtonia cedarbergensis WC SFBG #2004-0570 Evergreen, 20–22 m, 1,000–1,500 m, fynbos vegetation, South Africa Species Fig. Abbrevs. Location and Accession Phenology, Mature Tree Height, Native Elevational Range, and Habitat Athrotaxis laxifolia AL SFBG not cataloged Evergreen, 10–15 m, 1,000–1,200 m, montane forests, Tasmania Callitris rhomboidea CR SFBG #1999-0290 Evergreen, 10–15 m, 0–1,250 m, open woodland, Southeast Australia Calocedrus decurrens CD SFBG #XY-2004 Evergreen, 60 m, 50–2,960 m, mixed conifer forests, Oregon to Baja California Norte Cryptomeria japonica CJ SFBG not cataloged Evergreen, 50–60 m, 1–2,050 m, mixed evergreen forests, Japan Cupressus forbesii CF SFBG #1980-0055 Evergreen, 10 m, 210–1,400 m, chaparral, Baja California Norte Fitzroya cupressoides FC UCBG #2007.0165 Evergreen, 50–60 m, 4–1,000 m, emergent tree is evergreen rainforest, Chile Glyptostrobus pensilis GP UCBG #70.0169 Deciduous, 15–25 m, 1–730 m, river floodplains, deltas, Southern China Juniperus californica JC UCBG #83.0567 Evergreen, 7–10 m, 500–1,400 m, desert scrubland, Southern California to Baja California Norte Libocedrus plumosa LP UCSC #81.1172 Evergreen, 30–35 m, 1–600 m, lowland mixed angiosperm and conifer rainforests, New Zealand Metasequoia glyptostroboides MGS = stem MGR = root UCBG #49.0500 Deciduous, 35–50 m, 750 m, ravines and moist temperate forests, Central China Sequoiadendron giganteum SG UCBG #2002.1062 Evergreen, over 100 m, 1,400–2,150 m, western slopes of Sierra Nevada, California Sequoia sempervirens SS Campus, University of California, Berkeley, CA Evergreen, over 100 m, 1–750 m, moist and foggy climates, Central to Northern California coasts Taxodium distichum TD UCBG #60.1174 Evergreen, 90–95 m, 1,400–2,150 m, mixed conifer montane forests, Eastern California Taiwania cryptomerioides TC SFBG #1984.93, #1990.616 Evergreen, 60–65 m, 1,750–2,900 m, cool temperate forests, Asia Widdringtonia cedarbergensis WC SFBG #2004-0570 Evergreen, 20–22 m, 1,000–1,500 m, fynbos vegetation, South Africa Open in new tab Table I. Study species, figure abbreviations (Fig. Abbrevs.), locations (SFBG, San Francisco Botanical Garden, San Francisco; UCBG, University of California Botanical Garden, Berkeley, CA; UCSC, University of California, Santa Cruz, Arboretum, Santa Cruz, CA), and species natural history (Farjon, 2005) Species Fig. Abbrevs. Location and Accession Phenology, Mature Tree Height, Native Elevational Range, and Habitat Athrotaxis laxifolia AL SFBG not cataloged Evergreen, 10–15 m, 1,000–1,200 m, montane forests, Tasmania Callitris rhomboidea CR SFBG #1999-0290 Evergreen, 10–15 m, 0–1,250 m, open woodland, Southeast Australia Calocedrus decurrens CD SFBG #XY-2004 Evergreen, 60 m, 50–2,960 m, mixed conifer forests, Oregon to Baja California Norte Cryptomeria japonica CJ SFBG not cataloged Evergreen, 50–60 m, 1–2,050 m, mixed evergreen forests, Japan Cupressus forbesii CF SFBG #1980-0055 Evergreen, 10 m, 210–1,400 m, chaparral, Baja California Norte Fitzroya cupressoides FC UCBG #2007.0165 Evergreen, 50–60 m, 4–1,000 m, emergent tree is evergreen rainforest, Chile Glyptostrobus pensilis GP UCBG #70.0169 Deciduous, 15–25 m, 1–730 m, river floodplains, deltas, Southern China Juniperus californica JC UCBG #83.0567 Evergreen, 7–10 m, 500–1,400 m, desert scrubland, Southern California to Baja California Norte Libocedrus plumosa LP UCSC #81.1172 Evergreen, 30–35 m, 1–600 m, lowland mixed angiosperm and conifer rainforests, New Zealand Metasequoia glyptostroboides MGS = stem MGR = root UCBG #49.0500 Deciduous, 35–50 m, 750 m, ravines and moist temperate forests, Central China Sequoiadendron giganteum SG UCBG #2002.1062 Evergreen, over 100 m, 1,400–2,150 m, western slopes of Sierra Nevada, California Sequoia sempervirens SS Campus, University of California, Berkeley, CA Evergreen, over 100 m, 1–750 m, moist and foggy climates, Central to Northern California coasts Taxodium distichum TD UCBG #60.1174 Evergreen, 90–95 m, 1,400–2,150 m, mixed conifer montane forests, Eastern California Taiwania cryptomerioides TC SFBG #1984.93, #1990.616 Evergreen, 60–65 m, 1,750–2,900 m, cool temperate forests, Asia Widdringtonia cedarbergensis WC SFBG #2004-0570 Evergreen, 20–22 m, 1,000–1,500 m, fynbos vegetation, South Africa Species Fig. Abbrevs. Location and Accession Phenology, Mature Tree Height, Native Elevational Range, and Habitat Athrotaxis laxifolia AL SFBG not cataloged Evergreen, 10–15 m, 1,000–1,200 m, montane forests, Tasmania Callitris rhomboidea CR SFBG #1999-0290 Evergreen, 10–15 m, 0–1,250 m, open woodland, Southeast Australia Calocedrus decurrens CD SFBG #XY-2004 Evergreen, 60 m, 50–2,960 m, mixed conifer forests, Oregon to Baja California Norte Cryptomeria japonica CJ SFBG not cataloged Evergreen, 50–60 m, 1–2,050 m, mixed evergreen forests, Japan Cupressus forbesii CF SFBG #1980-0055 Evergreen, 10 m, 210–1,400 m, chaparral, Baja California Norte Fitzroya cupressoides FC UCBG #2007.0165 Evergreen, 50–60 m, 4–1,000 m, emergent tree is evergreen rainforest, Chile Glyptostrobus pensilis GP UCBG #70.0169 Deciduous, 15–25 m, 1–730 m, river floodplains, deltas, Southern China Juniperus californica JC UCBG #83.0567 Evergreen, 7–10 m, 500–1,400 m, desert scrubland, Southern California to Baja California Norte Libocedrus plumosa LP UCSC #81.1172 Evergreen, 30–35 m, 1–600 m, lowland mixed angiosperm and conifer rainforests, New Zealand Metasequoia glyptostroboides MGS = stem MGR = root UCBG #49.0500 Deciduous, 35–50 m, 750 m, ravines and moist temperate forests, Central China Sequoiadendron giganteum SG UCBG #2002.1062 Evergreen, over 100 m, 1,400–2,150 m, western slopes of Sierra Nevada, California Sequoia sempervirens SS Campus, University of California, Berkeley, CA Evergreen, over 100 m, 1–750 m, moist and foggy climates, Central to Northern California coasts Taxodium distichum TD UCBG #60.1174 Evergreen, 90–95 m, 1,400–2,150 m, mixed conifer montane forests, Eastern California Taiwania cryptomerioides TC SFBG #1984.93, #1990.616 Evergreen, 60–65 m, 1,750–2,900 m, cool temperate forests, Asia Widdringtonia cedarbergensis WC SFBG #2004-0570 Evergreen, 20–22 m, 1,000–1,500 m, fynbos vegetation, South Africa Open in new tab One of the first studies to examine the hydraulic resistance of the conifer pit used a physical model to show that 28% and 44% of pit resistance is explained by the torus and pit border (aperture), respectively, with the remaining 28% of pit resistance residing in the margo (Lancashire and Ennos, 2002). By contrast, computational fluid dynamics suggested that the pit aperture explains only 25% of pit resistance, with 25% to 38% resulting from the margo (Valli et al., 2002). The balance of remaining resistances was attributed to the internal architecture of the pit chamber. Both studies relied on physical or computational models that treated the margo as a homogeneously porous mesh, an approach that may have overestimated or underestimated the margo's contribution to pit resistance. This is not unexpected, because the margo is an intricate, irregularly porous structure that is difficult to replicate in a model. Compounding this complexity is an additional problem: despite one qualitative survey of pits from 120 gymnosperms (Bauch et al., 1972), very little is actually known about the structural variation of the margo, and even less about how this variation could relate to cavitation resistance. This uncertainty was broadly quantified by Hacke et al. (2004), who combined empirical data of cavitation resistance from a wide sampling of conifers with a model that treated the margo as a heterogeneous, but organized, mesh composed of pores of varying diameters. The assumption was that the xylem pressure at which membrane aspiration occurred (Pasp) was directly related to the porosity of the margo. Hence, an increase in the number of margo microfibril “spokes” reduced margo porosity, which increased Pasp, stabilized the torus, and thus conferred a higher resistance to cavitation. Consequently, reduced margo porosity was associated with greater cavitation resistance. Given that conifer tracheids are, on the whole, significantly shorter and narrower than angiosperm vessels, Hacke et al. (2004) recognized that it is essential for the structure of the pit membrane to be optimized for hydraulic efficiency, in addition to the basic requirement of cavitation safety (Hacke et al., 2004). To this end, the models of Hacke et al. (2004) suggested tight scaling between the pit aperture and torus diameter, whereby the torus-aperture overlap was sufficient to achieve a required resistance to cavitation without compromising pit hydraulic efficiency. Specifically, insufficient torus-aperture overlap required a dense margo to achieve a given air-seed pressure, while excessive overlap (due to increased torus diameter and smaller aperture diameter) reduced both the margo area available for water transport and the aperture conductance (Hacke et al., 2004). Interestingly, increased rupture of the margo microfibrils during membrane aspiration was another consequence of excessive overlap, because shorter microfibrils were subject to stretching beyond their inherent tensile strength. Since conifers can experience repeated cycles of cavitation and embolism (Sperry et al., 1994; Mayr et al., 2002), suggesting that the pit membrane can rebound from an aspirated position (Sperry and Tyree, 1990), it seems reasonable to assume that the key constituents of the torus-margo pit membrane have evolved to scale in a manner that optimizes the safety/efficiency tradeoff in light of fixed, biomechanical limitations imposed by the properties of cellulose. Although the degree of margo variation was unknown to Hacke et al. (2004), they presented important ideas about the functional morphology of conifer pit membranes that have since been confirmed. First, the notion that it is the torus-aperture overlap that determines cavitation safety was recently validated in three species of Pinaceae as well as in stems of Douglas fir (Pseudotsuga menziesii) at different heights (Domec et al., 2008; Hacke and Jansen, 2009). Second, recent studies have shown that the pit aperture controls pit hydraulic conductivity and transport efficiency in the distally located xylem of tall Douglas fir trees as well as across a range of cavitation pressures in different Douglas fir organs (Domec et al., 2006, 2008). By contrast, estimates of the margo's contribution to pit resistance have come about by indirect calculations rather than empirical observation. In a survey study that evaluated pit resistance across 19 species of conifers, Pittermann et al. (2006b) concluded that, on average, the pit membrane (torus and margo) probably explains a large fraction of total pit resistance, with less than 10% attributed to the aperture, a major deviation from the results of Hacke et al. (2004) and the other studies cited above. Just how does margo porosity contribute to pit membrane resistance? In this study, we partitioned the margo and aperture contribution to pit hydraulic resistance by first measuring cavitation resistance in the distal stems and one root belonging to 15 species of Cupressaceae and then combining these data with anatomical measurements obtained via scanning electron microscopy (SEM) and transmission electron microscopy (TEM) of earlywood intertracheid pits. Specifically, we tested the hypothesis that an efficiency-versus-safety tradeoff exists at the pit level and that it is related to the variation in margo structure. We also examined additional features of the pit membrane in order to identify the anatomical characters that combine to affect cavitation resistance and pit hydraulic performance. The Cupressaceae are an ideal system for exploring the tradeoffs between safety and efficiency in pit membrane structure. They span a wide range of environments, from riparian habitats with ample moisture to desert habitats where precipitation is variable and droughts occur frequently. They are an excellent model for evolutionary studies because their phylogeny is well understood (Gadek et al., 2000) and they have a rich fossil record dating to the Jurassic (Stockey et al., 2005). Conclusions drawn from their anatomy may be all the more relevant as they are the only family of conifers with a worldwide distribution (Farjon, 2005). By sampling broadly across the Cupressaceae phylogeny, we present the functional adaptations of pit membranes across a range of cavitation pressures and can comment on the evolutionary trends of pit morphology in this family. RESULTS Species vulnerability to cavitation was determined by measuring the segment percentage loss of conductivity in response to progressively more negative xylem pressures induced by centrifugation (Alder et al., 1997; see “Materials and Methods”). The xylem pressures at which stems showed 50% loss of conductivity (cavitation pressure [P50]) spanned a broad range from −2.8 ± 0.62 MPa (mean ± sd) in the semiriparian Glyptostrobus pensilis to a low of −11.3 ± 3.52 MPa recorded in xeric-adapted Widdringtonia cedarbergensis. The P50 of Metasequoia glyptostroboides roots was −0.16 ± 0.06 MPa. SEM images of earlywood intertracheid pit membranes revealed a high degree of variation in membrane structure (Fig. 1). Qualitatively, increasing cavitation resistance was associated with more distinct torus-margo boundaries, thicker tori, torus extensions (Juniperus, Widdringtonia), and the presence of protrusions on the walls of the pit chamber. We also observed that membranes belonging to vulnerable taxa such as Taxodium, Metasequoia (roots especially), and Glyptostrobus were much more fragile and had a tendency to tear and deform when subjected to the electron beam, in contrast to the more cavitation-resistant species, whose membranes remained mostly intact. Although cavitation resistance visually appears to be related to greater margo porosity (compare Athrotaxis with Widdringtonia in Fig. 1), the actual margo structure proved to be so variable that any trends related to cavitation resistance were obscured (see below). Our measurements of intertracheid pit features on these SEM images revealed no correlation between pit diameter and P50 and only a weak relationship between torus diameter and P50 (Fig. 2, A and B). However, we found that aperture diameter decreased nearly 2-fold with more negative cavitation pressures, ranging from 4.5 ± 0.54 μm for M. glyptostroboides roots to 2.3 ± 0.35 μm from stems of Cupressus forbesii (P50 = −11.2 ± 1.78 MPa; Fig. 2C). Figure 2. Open in new tabDownload slide The relationship between species cavitation pressure and pit diameter (A), torus diameter (B), and pit aperture diameter (C). Data represent average values per species ± 1 sd (n = 15–20 pits per species). Figure 2. Open in new tabDownload slide The relationship between species cavitation pressure and pit diameter (A), torus diameter (B), and pit aperture diameter (C). Data represent average values per species ± 1 sd (n = 15–20 pits per species). The ratio of the torus to aperture diameter showed the strongest and most statistically significant relationship with P50 across the 15 species surveyed, suggesting that the sealing function of the torus is essential to controlling air seeding in Cupressaceous conifers (Fig. 3). Specifically, it is the reduction in aperture diameter that determined the extent of torus-aperture overlap, which was calculated according to Hacke et al. (2004) and ranged from 0.01 in Athrotaxis laxifolia to 0.45 in W. cedarbergensis (data not shown). We calculated an average overlap of 0.19 ± 0.13, not dissimilar from the 0.24 to 0.3 range predicted to optimize pit conductivity while minimizing air seeding via margo strand rupture (Hacke et al., 2004). Figure 3. Open in new tabDownload slide The torus diameter-to-pit aperture diameter ratio as a function of cavitation pressure calculated using the anatomical data presented in Figure 2. Overlap between the torus and aperture increases with more negative cavitation pressure. Figure 3. Open in new tabDownload slide The torus diameter-to-pit aperture diameter ratio as a function of cavitation pressure calculated using the anatomical data presented in Figure 2. Overlap between the torus and aperture increases with more negative cavitation pressure. Following standard statistical analysis, the physiological and anatomical data were also analyzed in a phylogenetic context in order to examine whether any of the traits under consideration might have evolved in a correlated manner. In addition, because this study looks at a group of closely related plants, the Cupressaceae, we wanted to test whether phylogenetic relationships within our data either obscured any relationships between traits or, alternatively, led to a false impression of correlation (Felsenstein, 1985; Garland et al., 1992). In order to test for both these effects, we used phylogenetically independent contrasts (PICs; Felsenstein, 1985). This approach has previously been used with hydraulic traits by several authors (Preston and Ackerly, 2003, Maherali et al., 2004; Pratt et al., 2007; Willson et al., 2008). Results from the PIC analysis were largely similar to those from the nonphylogenetic (Tip) analysis, and all relationships between P50 and anatomical data exhibited the same direction whether analyzed as Tips or PICs (Supplemental Figs. S1 and S2; Supplemental Table S1). Surprisingly, there was little evidence of correlated evolution in any of the pit anatomical characters. However, we found that the key torus-aperture ratio showed the strongest evidence of correlated evolution with P50 (Fig. 4). As in Figure 2, there was an indication that this ratio was controlled primarily by changes in aperture (P = 0.014, r 2 = 0.38) rather than by torus diameter (P = 0.047, r 2 = 0.27), but other than this, none of the other independent anatomical characters had any correlation coefficient greater than 0.3. Figure 4. Open in new tabDownload slide The phylogeny (left) and phylogenetic independent contrasts (right) of the Cupressaceae species considered in this study. The number on each node corresponds to the independent contrasts data (right) calculated from the P50 and torus-aperture ratio values presented in Figure 3. Figure 4. Open in new tabDownload slide The phylogeny (left) and phylogenetic independent contrasts (right) of the Cupressaceae species considered in this study. The number on each node corresponds to the independent contrasts data (right) calculated from the P50 and torus-aperture ratio values presented in Figure 3. TEM measurements of seven select species that spanned the range of P50 values revealed no relationship between margo thickness and cavitation resistance (Figs. 5 and 6), with thickness averaging 194.8 ± 27.1 nm. By contrast, the thickness of the torus doubled in response to cavitation resistance, ranging from 288 ± 215 nm in A. laxifolia (P50 = −3.45 MPa) to 666 ± 118 nm in W. cedarbergensis (P50 = −11.3 MPa). Interestingly, we observed some variation in the electron density of the pit membranes, whereby some tori showed regions of less lignified wall material, such as in A. laxifolia and Taxodium distichum and the middle lamella of W. cedarbergensis (Fig. 5). These qualitative differences in membrane composition, however, could not be functionally linked to a species’ cavitation resistance. Lastly, our TEM analysis revealed that, overall, pit border thickness and aperture depth are weakly correlated with cavitation resistance (Fig. 6), mostly because P50, membrane features, and pit border thickness are decoupled in A. laxifolia. Figure 5. Open in new tabDownload slide Transmission electron micrographs of pit membrane cross-sections belonging to six (of seven) selected species of Cupressaceae exhibiting a range of cavitation pressures. A high degree of variation exists in aperture (A) diameter and depth and the thickness of the torus (T), but aside from the torus extension (TE) present in W. cedarbergensis, margo (M) thickness remained invariable. Figure 5. Open in new tabDownload slide Transmission electron micrographs of pit membrane cross-sections belonging to six (of seven) selected species of Cupressaceae exhibiting a range of cavitation pressures. A high degree of variation exists in aperture (A) diameter and depth and the thickness of the torus (T), but aside from the torus extension (TE) present in W. cedarbergensis, margo (M) thickness remained invariable. Figure 6. Open in new tabDownload slide The relationship between cavitation pressure and pit border thickness (A), aperture depth (B), margo thickness (C), and torus thickness (D) in seven species of Cupressaceae (n = 7–10 pits per species). Figure 6. Open in new tabDownload slide The relationship between cavitation pressure and pit border thickness (A), aperture depth (B), margo thickness (C), and torus thickness (D) in seven species of Cupressaceae (n = 7–10 pits per species). Our approach of combined SEM and TEM image analysis with hydraulic methods allowed us to directly test the specific predictions of Hacke et al. (2004) about the structure-function relationships of torus-margo pit membranes. By directly measuring features such as aperture diameter, aperture depth, margo strand length, thickness, and margo porosity, we could identify the hydraulic and biomechanical tradeoffs associated with cavitation resistance across the Cupressaceae. The first of these analyses revealed that pit aperture hydraulic resistance is tightly and significantly correlated with cavitation resistance, decreasing nearly 8-fold over the 10 MPa range of a species’ P50 values (Fig. 7A). This is consistent with the observed adjustments in aperture diameter described in Figure 2 and in other studies (Domec et al., 2008). Second, margo hydraulic resistance was shown to be invariable across the spectrum of cavitation pressures, indicating that margo structure is generally unrelated to cavitation resistance in these plants (Fig. 7B). Although surprising, we are confident that these results are an accurate and methodologically unbiased reflection of margo resistance, because in order to minimize variation due to artifacts, margo porosity was consistently measured in the most intact portions (30%–50% of total margo area) of the membrane across all of the species we surveyed (although the best photographs are shown in Fig. 1). Thus, we conclude that the reduction in total pit resistance with less negative P50 values is most strongly related to increasing pit aperture diameter (Fig. 7C). It is important to note, however, that the proportional contribution of the margo to overall pit resistance increases in pits of more vulnerable taxa, presumably because the relative pit aperture resistance is so low (Fig. 8). Figure 7. Open in new tabDownload slide Pit aperture resistance (A), margo resistance (B), and total pit resistance (C) plotted as a function of cavitation pressure. The aperture resistance is calculated on the basis of two apertures per pit. Figure 7. Open in new tabDownload slide Pit aperture resistance (A), margo resistance (B), and total pit resistance (C) plotted as a function of cavitation pressure. The aperture resistance is calculated on the basis of two apertures per pit. Figure 8. Open in new tabDownload slide The percentage contribution of the margo to the total resistance of the torus-margo pit membrane increases in taxa that are more vulnerable to cavitation (see text for details). Figure 8. Open in new tabDownload slide The percentage contribution of the margo to the total resistance of the torus-margo pit membrane increases in taxa that are more vulnerable to cavitation (see text for details). Our final analysis examined the membrane-level biomechanical costs associated with varying vulnerability to cavitation. We predicted that the margo fibrils of aspirated membranes might be subject to increasing tensions in species with progressively greater resistance to cavitation. However, we observed no relationship between maximum margo tension at deflection and P50 across the seven species surveyed (Fig. 9). This result is probably due to the scaling of margo strand length (l) with the maximum torus displacement (y), whereby displacement is roughly one-third the length of the average margo strand (Fig. 9). Consequently, the margo microfibrils experience constant tension across a range of P50 values. To look at this from a different perspective, we also calculated the margo strain at aspiration (ea) using pit membrane and torus diameter (Hacke et al., 2004) to show that ea is, in fact, invariable across P50 values (Fig. 9). Figure 9. Open in new tabDownload slide A, The maximum tension calculated from margo strand length and torus displacement indicates that individual margo strand is invariable across a range of cavitation pressures. B, Presumably, this is because the average length of the margo strands scales with the mean maximum displacement of the pit membrane during aspiration. C, Calculations of margo strain at aspiration based on measurements of membrane and torus diameter suggest that deformation of the margo is decoupled from the cavitation pressure. Figure 9. Open in new tabDownload slide A, The maximum tension calculated from margo strand length and torus displacement indicates that individual margo strand is invariable across a range of cavitation pressures. B, Presumably, this is because the average length of the margo strands scales with the mean maximum displacement of the pit membrane during aspiration. C, Calculations of margo strain at aspiration based on measurements of membrane and torus diameter suggest that deformation of the margo is decoupled from the cavitation pressure. DISCUSSION Our primary objective was to evaluate the pit-level hydraulic tradeoffs associated with cavitation resistance across the Cupressaceae and to specifically address the respective contributions of the margo and pit aperture to total pit hydraulic resistance. Surprisingly, we found no decrease in margo hydraulic efficiency with increasing resistance to cavitation, a result that stems from the high degree of variation apparent in margo strand arrangement in the surveyed taxa (Fig. 7). It seems reasonable to predict that when the membrane is aspirated, a margo composed of thicker strands and/or a denser strand meshwork would stabilize the torus against the aperture and thus prevent air seeding from occurring. Indeed, a denser margo is an important assumption in the modeling work by Hacke et al. (2004) and was, in fact, observed in open-grown, cavitation-resistant individuals of Pinus contorta and Picea glauca relative to those grown in the shade (Schoonmaker et al., 2010). In the same study, shade-grown plants exhibited greater vulnerability to cavitation, a response thought to be a function of the sparse and thin-stranded margo that may have rendered the membrane susceptible to slipping or tearing under negative xylem pressures (Schoonmaker et al., 2010). Given that margo structure in the Cupressaceae appears to be functionally decoupled from P50, pit function in the Cupressaceae requires an alternative explanation in addition to the current model, which appears to apply to the Pinaceae. Recent data suggest that aside from the margo, Pinaceae and Cupressaceae also appear to differ in the structure of the torus, despite showing similar torus-aperture ratios at intermediate P50 values (Fig. 3; Hacke and Jansen, 2009). Hacke and Jansen (2009) reported a decrease in torus thickness associated with greater cavitation resistance in three species of Pinaceae, suggesting that thinner, and presumably more flexible, tori may form a tighter seal against the pit aperture during aspiration. By contrast, our results indicate a substantial thickening of the torus in response to more negative P50 values (Figs. 5 and 6). In cavitation-resistant Cupressaceae, once the torus is aspirated and resting against the pit border, the margo could be, up to a point, unnecessary for keeping the torus in place, because the tori of these plants tend to be thick and potentially less susceptible to distortion (Fig. 3). Second, the presence of knobby wall protrusions on the internal side of the pit border may improve the wettability and adhesive properties of the chamber wall (Fig. 1; Kohonen, 2006; Heady et al., 2008; Kohonen and Helland, 2009), so the aspirated torus may actually be quite stable under negative xylem pressures. However, the pit membranes of highly cavitation-resistant species such as Widdringtonia and Juniperus consistently possess torus extensions, which presumably keep the torus from slipping under very negative xylem tensions. Interestingly, Schoonmaker et al. (2010) interpreted the torus extensions in shade-grown Pinaceae differently, regarding them as lesions that increase vulnerability to air seeding by creating a weaker seal between the torus and the aperture. The presence of thicker and extended tori in cavitation-resistant Cupressoid and Callitroid clades suggests that in these groups, selection favors pit membranes that minimize air seeding through the seal between the torus and the pit border. In less cavitation-resistant species such as T. distichum and G. pensilis, the role of the torus in protecting against air seeding is much less clear, because not only is the torus-aperture ratio much lower but the torus itself is thin, poorly differentiated from the margo, vulnerable to deformation as evidenced by its tendency to tear under the SEM electron beam, and often porous, not unlike the homogenous pit membrane of angiosperms (Jansen et al., 2009). It is not unreasonable to suggest that in semiriparian Cupressaceae, air seeding may occur through the minute pores of an aspirated torus or by rupture of the membrane under negative pressure, as suggested by Hacke et al. (2004) and Cochard et al. (2009). In these deciduous conifers, where consistent water availability selects for large pit apertures and high pit transport efficiency (Fig. 7), the margo probably lends little mechanical stability to the aspirated torus. We were surprised at the high degree of variation in margo resistance (R m), but our measurements of margo porosity were consistently made on the most intact regions of the margo across all species, so our estimates of R m are conservative. The variation in R m arises from the presence of one or several large pores, which can have a large impact on R m, because the hydraulic conductance of a pore (D) is a function of D 3 (see Eq. 3 below). Yet, despite this variation, it is not unreasonable to conclude that in the Cupressaceae, the margo's contribution to overall pit hydraulic resistance progressively increases in less cavitation-resistant taxa (Fig. 8). Although it is difficult to identify a precise P50 at which margo resistance exceeds aperture resistance, our data suggest that the margo may constrain pit transport efficiency in more mesic-riparian Cupressaceae, or in parts of the tree such as the trunk, roots, or latewood that are typically less resistant to cavitation (Domec et al., 2006). Generally speaking, margo resistance may be higher than aperture resistance in Cupressaceae with P50 values that are less negative than about −5.5 MPa or that have apertures greater than 4 μm in diameter. This increased role of the margo in less cavitation-resistant Cupressaceae is a direct consequence of a concurrent decrease in pit aperture resistance rather than a margo-level trend associated with P50. The idea that the pit aperture controls most of the pit hydraulic efficiency was first proposed by Lancashire and Ennos (2002) using a physical model and later developed by Hacke et al. (2004), who showed that aperture resistance was consistently higher than the modeled margo resistance in pits of north temperate conifers. Their work suggested that a reduction in aperture diameter was responsible for a 1.7-fold increase in pit resistance across a range of cavitation pressures (Hacke et al., 2004). The similar conclusion of Domec et al. (2008) that the aperture controls pit hydraulic efficiency across a range of P50 values in Douglas fir echoed what seems an obvious and fundamental pit-level hydraulic tradeoff associated with cavitation resistance in conifers. At first glance, our data also support these findings, but closer inspection with SEM suggests that the underlying assumptions of what controls pit hydraulic resistance in the Cupressaceae are more complex, probably because of developmental limitations on the control of margo porosity. Manipulative experiments suggest that margo porosity may be more tightly regulated in the Pinaceae (Schoonmaker et al., 2010), but the degree of pit membrane plasticity in the Cupressaceae is unknown. Consistent with previous work, our data show that vulnerability to cavitation is strongly controlled by the ratio of torus to aperture diameter, that is, the degree of torus-aperture overlap (Fig. 3; Domec et al., 2008; Hacke and Jansen, 2009). While torus diameter changes modestly across P50, it is the variation in aperture diameter that sets the torus-aperture overlap and ultimately determines the air-seeding resistance of torus-margo pit membranes. The observed increase in torus thickness is consistent with the need for the membrane to structurally withstand more negative xylem pressures in xeric-adapted taxa such as Widdringtonia and Juniperus. A similar correlation between pit membrane thickness and vulnerability to air seeding was recently found for angiosperms (Jansen et al., 2009). By contrast, Hacke and Jansen (2009) observed that thinner tori were associated with an increase in cavitation resistance in three species of Pinaceae, so the phylogenetic influence on the structure and function of conifer pit membranes needs to be broadly examined. Although our survey of pit membranes in the Cupressaceae is by no means exhaustive, we can comment on the general evolutionary trends in pit structure and function across members of this family. Broadly speaking, the higher cavitation resistance of the Cupressoid and Callitroid grades is a derived character, reflected in part by the trend toward higher torus-aperture ratios, thicker tori, and a progressively more targeted digestion and distinct separation of the torus and the margo (Fig. 1). More quantitatively, the PIC analyses of the cavitation pressure and torus-aperture ratio give us some capacity to trace the evolution of drought resistance in this family. In particular, drought tolerance would most likely have been selected for as the Callitroid and Cupressoid grades of the Cupressaceae radiated into drier habitats of the Paleogene from the semiriparian habitats of their Cretaceous-era Sequoioid and Taxodioid ancestors, as suggested by the fossil record (Florin, 1963; Farjon, 2005; Stockey et al., 2005; see also Willis and McElwain, 2002). Indeed, one of the largest P50 divergences between sister clades is found at the divergence between the Taxodioid clade and the Cupressoid and Callitroid clades at node 21. The PIC results also suggest that drought tolerance has been separately derived in the Cupressoid and Callitroid clades, with the largest differences in P50 and the torus-aperture ratio occurring at the split between arid and mesic taxa, such as in nodes 26 and 27 (two of the largest contrasts) of the Callitroids and node 28 of the Cupressoids (Fig. 4). By comparison, the differences between the riparian taxa (Taxodioid clade) and their mesic relatives (Sequoioid clade) are smaller, as indicated by smaller P50 contrasts at nodes 22 and 19, respectively. It seems likely that the greater vulnerability to cavitation found in the basal Cupressaceae is associated with the high water availability in the mesic-riparian habitats in which these taxa are often found and that this is the ancestral condition for the entire family. Consistent with this, the deciduous Glyptostrobus, Taxodium, and Metasequoia, which grow in semiriparian regions, tend to exhibit thin, fragile, and seemingly diaphanous pit membranes that appear poorly equipped to withstand low water potentials. Thus, both anatomical and independent contrast results support the conclusion that low cavitation resistance is the ancestral condition for the Cupressaceae and that drought resistance is a derived character that evolved separately in the Callitroid and Cupressoid clades. Future studies linking additional xylem-level traits with ancestral state reconstructions may add to our understanding of the evolutionary trajectory of the Cupressaceae. Lastly, we observed no biomechanical tradeoffs associated with cavitation resistance and the tension imposed on the margo strands of the pit membrane (Fig. 9). The forces imposed upon the margo during aspiration remain invariable because the length of the strands (l) scales with the membrane displacement distance (y) and torus diameter, such that the ratio of l to y remains constant. Our data suggest that, structurally, cavitation-resistant Cupressaceae pit membranes can be classified as type 1 pits, which are for the most part strong, flexible, and resistant to rupture upon aspiration (Hacke et al., 2004). Whether or not “stretch seeding” consistently occurs (i.e. air seeding via stretching of the margo and displacement of the torus) is debatable, because more cavitation-vulnerable taxa such as Taxodium and Glyptostrobus possess pit membranes that are rather fragile, so capillary failure may not be limited to the torus-pit border interface (Cochard et al., 2009). The scaling between pit features is also consistent with the observed proportions between the P50 and the conduit diameter-to-wall thickness ratio (Hacke et al., 2001; Pittermann et al., 2006a, 2006b), aperture and torus diameter, as well as aperture depth and pit border thickness (Fig. 5). This lends support to the idea that, in gymnosperms, where transport is limited to unicellular tracheids, selection has acted to optimize pit function to fulfill both safety and transport requirements (Hacke et al., 2004; Domec et al., 2006, 2008; Sperry et al., 2006). Where are we with regard to linking pit structure with cavitation resistance across woody plants? The structural heterogeneity of the torus-margo pit makes it an attractive feature to study because theoretical predictions can be reasonably tested using SEM and hydraulic methods, so key traits across all torus-margo pit membranes as well as important differences at the level of the family or genus can be identified among extant conifers. Second, if we are to understand the acclimation response of cavitation resistance, it is important to also consider phenotypic plasticity via strategically manipulative experiments such as that of Schoonmaker et al. (2010). For example, our study was conducted on common garden plants grown in a mild climate, but we expect that some species grown in their native environment may have dramatically different P50 values (i.e. Callitris rhomboidea; Brodribb and Cochard, 2009), and this may be reflected in pit structure. Lastly, it is becoming apparent that the process of air seeding may be, in some sense, convergent in both angiosperms and conifers and dependent in part on the thickness of the pit membrane. Jansen et al. (2009) show a tight relationship between air seeding, membrane thickness, and porosity in the homogenous pit membranes of angiosperms that conforms to the predictions made by Sperry and Hacke (2004). Our observations of thin membranes with poorly distinguished tori in deciduous, riparian Cupressaceae suggest that similar patterns may hold in these plants. In this regard, advancing our understanding of homogenous pit membranes, especially the safety/efficiency tradeoffs with respect to cavitation, will fill in an important gap in our understanding of the adaptive nature of interconduit pit membranes in general. This project focused on the structure/function tradeoffs in pit-level traits, but it is clear that cavitation resistance encompasses a suite of well-studied developmentally coordinated functional traits such as conduit wall thickness, conduit allometry, wood density, and overall physiology (Hacke et al., 2001, 2004; Santiago et al., 2004; Pittermann et al., 2006a, 2006b; Sperry et al., 2006; Jacobsen et al., 2008; Schoonmaker et al., 2010), all of which combine to determine the safety/efficiency tradeoff. Functionally, this means that cavitation-resistant pits cannot be associated with large, thin-walled tracheids that are hydraulically efficient and vulnerable to implosion under tension. A manifestation of this scenario would require that species with efficient xylem such as T. distichum possess pit membranes with narrow apertures, resulting in a combination of traits that would significantly diminish any hydraulic gains associated with large tracheids (Hacke et al., 2004; Sperry et al., 2006). Future work incorporating both the fine structure of pit membranes and tissue-level traits will be required to fully unravel the complex relationships that underlie adaptive variation in cavitation resistance and hydraulic efficiency at the whole plant level. MATERIALS AND METHODS Plant Material All Cupressaceae species were collected from arboreta located on the central California coast, so all plants were subject to similar climatic conditions (Supplemental Table S1). Stem samples from a total of 15 species were collected from Strybing Arboretum in San Francisco's Golden Gate Park (37.78°N 122.43°W) and from the nearby University of California's Botanical Garden in Berkeley. Libocedrus plumosa was sampled at the University of California's arboretum in Santa Cruz (36.97°N 122.03°W). In order to examine the intertracheid pits of purely riparian conifer xylem, we collected submerged root material from Metasequoia glyptostroboides trees growing creekside at the University of California's Botanical Garden. In general, the climate of these three sites is described as Mediterranean, with mild temperatures, wet winters, dry summers, abundant summer fog input, and very rare freezing events (mean annual precipitation is 500–580 mm, and mean annual temperature is 15°C–17°C). However, species’ watering regimes reflected their natural history, such that individuals of Taxodium distichum were typically found near ponds, while xeric-adapted taxa such as Widdringtonia cedarbergensis grew in gardens that were less frequently watered. All species were clearly labeled, associated with accession numbers, and their age and planting history were well catalogued. Juvenile stems ranging from 5 to 8 mm in diameter were collected from sun-exposed, vigorous regions of mature tree canopies, from heights of 1.5 to 4 m. Five to eight branches were clipped from two to six individuals (depending on permission) and were immediately wrapped in a plastic sac with a wet paper towel until further processing for hydraulic measurements in the laboratory. With the exception of the deciduous conifers (M. glyptostroboides, T. distichum, Glyptostrobus pensilis), which were sampled in June 2006, random species sampling occurred from April to August in 2006 and 2007 in order to account for seasonal variation in cavitation resistance (Kolb and Sperry, 1999). Hydraulic Measurements Sampled stems were recut under water to a length of 142 mm, and the distal ends were shaved smooth with a razor blade. Although native embolism is typically less than 10% in well-watered conifer stems, air in the xylem was removed by submersing the segments in distilled and filtered water (0.22 μm; E-Pure filtration system; Barnstead International) and degassing them overnight under “house” vacuum. Stems were degassed rather than flushed because the latter method often reduces hydraulic conductivity (Mayr et al., 2002; Pittermann and Sperry, 2003; Willson et al., 2008). Hydraulic conductivity (k) was measured according to the method of Sperry (1993) and calculated as the flow rate for a given pressure gradient per unit of stem length. The segments were mounted on a tubing apparatus where k was measured gravimetrically under a pressure of 5 to 6 kPa using filtered water. The flow rate through the segments was determined without a pressure head before and after each gravimetric flow measurement. These background flows were averaged and subtracted from the pressure-induced flow in order to improve accuracy. We used the centrifuge method to determine species vulnerability to cavitation in response to a range of xylem pressures (Pockman et al., 1995; Alder et al., 1997). Stems were secured in a custom rotor designed to fit a Sorvall RC-5C centrifuge and spun for 3 min at speeds that induce a known xylem pressure (Px). The percentage loss of conductivity (PLC) caused by centrifugation at each Px was calculated from the k measured after spinning, relative to the maximum conductivity (k max) at Px = 0 MPa, such that (1) where k max was determined at Px = 0 MPa following degassing. The segments were spun to progressively more negative Px until the PLC exceeded 90% or, alternatively, until Px = −10 MPa, which is the most negative Px that can be achieved using the centrifuge. Five to six segments were used to generate the same number of vulnerability curves, each of which was subject to the curve-fitting routine of Pammenter and Vander Willigen (1998). The xylem pressure at which segments exhibited a 50% loss of conductivity (P50) was computed as average ± 1 sd per species. SEM Wood samples for SEM were collected from the same individuals and similar canopy locations as described above. Four stem segments that were 5 to 8 mm in diameter and 10 to 15 cm in length were collected from each species, wrapped in a wet paper towel, and transported to the laboratory, where they were recut under water to a 5-cm length and placed in a −10°C freezer for 2 d. These samples were then freeze dried overnight. Freeze-dried twig samples were split in half and mounted on aluminum stubs using nail polish to increase conductivity between the sample and stub. Samples were coated with gold-palladium for 4 min at 20 mA using a sputter coater (Emitech). All samples were observed with a field emission scanning electron microscope (Zeiss UltraPlus Analytical) with an accelerating voltage of 3 kV. A minimum of 15 to 20 earlywood intertracheid pit membranes were photographed for analysis. TEM Wood samples for TEM were collected from similar canopy locations from the following species: Athrotaxis laxifolia, Callitris rhomboidea, Cryptomeria japonica, Cupressus forbesii, Sequoia sempervirens, T. distichum, and W. cedarbergensis. Stem samples measuring approximately 5 to 8 mm in diameter and 5 cm in length were clipped (two stems per species), wrapped in moist paper towels, and immediately mailed by overnight courier to the Jodrell Laboratory at the Royal Botanic Gardens in Kew. Due to customs inspection, the samples were in transit for 3 d. At Kew, one sample per species was prepared for TEM. Wood from the last two growth rings was cut into thin longitudinal slivers, cut into 1-mm3 blocks, and fixed overnight in Karnovsky's fixative at room temperature (Karnovsky, 1965). After washing in 0.05 m phosphate buffer, the specimens were postfixed in 1% buffered osmium tetroxide for 4 h at room temperature, washed again, and dehydrated through a graded ethanol series (30%, 50%, 70%, 90%, 100%). The ethanol was gradually replaced with LR White resin (London Resin Co.) over several days, with the resin being changed approximately every 12 h. The resin was polymerized in a Gallenkamp vacuum oven at 60°C and 1,000 mbar for 24 h. Embedded samples were trimmed with a Leica EM Specimen Trimmer (Leica Microsystems) and sectioned on an ultramicrotome (Ultracut; Reichert-Jung). Transverse sections of earlywood conduits about 1 and 2 μm thick were cut with a glass knife, heat fixed to glass slides, stained with 0.5% toluidine blue-O in 0.1 m phosphate buffer, and mounted in DPX (containing distyrene, dibutyl phtalate, and xylene; Agar Scientific). Ultrathin sections were cut between 60 and 90 nm using a diamond knife. These sections were attached to Formvar (Agar Scientific) and copper grids (100 mesh) and stained with uranyl acetate and lead citrate using a Leica EM Stain Ultrostainer (Leica Microsystems). Observations were carried out using a JEOL JEM-1210 transmission electron microscope at 80 kV accelerating voltage, and digital images were taken using a MegaView III camera (Soft Imaging System). Measurements were made on at least 15 pit membranes per species. Phylogenetic Reconstruction Sequences for two markers (rbcL and MatK) from Gadek et al. (2000) were downloaded from GenBank. Sequences were aligned using ClustalW (Thompson et al., 1994; EMBL-EBI, 2009). Mesquite (Maddison and Maddison, 2009) was used to touch up alignments and concatenate both genes into a matrix, which was then analyzed using maximum parsimony in PAUP* beta 11 (Swofford, 1998), with a random addition sequence and tree bisection and reconnection branch swapping. We were unable to recover the topology found by Gadek et al. (2000) using the smaller number of taxa included in this physiological study. The data were reanalyzed using a constraint tree based on that presented by Gadek et al. (2000), resulting in a single phylogenetic topology with molecular branch lengths. This topology was used for the subsequent comparative analysis. Comparative Phylogenetic Analysis The phylogeny was imported into the R statistical environment (R Development Core Team, 2008) using the APE package (Paradis et al., 2004) and used to analyze the physiological data presented in this paper using PICs (Felsenstein, 1985). A single mean measurement for each genus and each trait was used in the PIC analysis, with the exception that only the branch data and not the root data from M. glyptostroboides was used. PIC analysis assumes a Brownian model of branch length evolution. Diagnostic tests (from the Picante package; Kembel et al., 2010) showed that the molecular branch lengths resulting from the combined PAUP* analysis did not obey this assumption. Nonparametric rate smoothing (Sanderson, 2002; Paradis et al., 2004) did not address this problem, so unit branch lengths were used, which, when significance levels for correlation coefficients are adjusted, do not bias estimates of correlated evolution (Ackerly, 2000). Further statistical analysis was conducted on the resulting contrasts using APE and Picante. R scripts used are available from S.A. Stuart on request. Estimates of Pit Hydraulic Resistance and Margo Strand Tension In order to quantify the hydraulic tradeoffs associated with cavitation resistance at the pit level, we calculated the aperture, margo, and pit resistances using methods described by Sperry and Hacke (2004) and Hacke et al. (2004). The hydraulic resistance of one pit aperture (R a) was calculated according to: (2) where t a is the pit aperture depth, ν is the viscosity of water (0.001 Pa × s at 20°C), and D a is the aperture diameter. For taxa where t a could not be obtained from TEM images, t a was assigned the average value of 0.7 μm. The margo resistance (R m) was computed according to Sperry and Hacke (2004), whereby (3) where n is the number of pores in the membrane, D p is the hydraulic diameter for all the margo pores measured (assuming an equivalent circle diameter of each of the pores), and f(h) is the proportion of the margo occupied by pores. SEM images were used to estimate margo porosity, and pore areas were measured within an intact, nonappressed section of the margo from an average of four pits using image-analysis software. An average of 350 pores was measured on each photograph. Total pit area-specific resistance was calculated as the sum of R m and 2 × R a, where 2 × R a accounts for the two apertures present in one pit. We applied the equation of Petty (1972) to calculate the maximum tension experienced by a single margo strand when the membrane is fully appressed against the pit border. For one radial strand, the tension (T) is estimated according to (4) where E is the elastic modulus of the margo strand (estimated at 5 GPa from Hacke et al. [2004]), A is the average strand cross-sectional area (2.987 × 10−18 m2) based on the average margo strand thickness (195 nm), l is the length of the margo strand, and y is the displacement of the torus against walls of the pit border. Displacement was measured according to Petty (1972) as the distance between the unaspirated, neutrally placed pit membrane and the pit border location where the torus presses against the pit border during aspiration (approximately the center of the pit border). Margo spoke strain at aspiration (ea) was calculated as (5) where D m and D t represent membrane and torus diameter, respectively (Hacke et al., 2004). Sequence data from this article can be found in the GenBank/EMBL data libraries under accession numbers AF127426 to AF127431, L12537, L12538, L12569 to L12573, L12575 to L12578, L12580, L13172, L25750, L25751, L25754 to L25757, S75127, and X63660. Supplemental Data The following materials are available in the online version of this article. Supplemental Figure S1. A matrix of hydraulic and SEM-based anatomical trait relationships analyzed using PICs. Supplemental Figure S2. A matrix of hydraulic and TEM-based anatomical trait relationships analyzed using PICs. Supplemental Table S1. Comparison of standard statistical analysis (Tips) and PICs. Supplemental Table S2. GenBank accessions used in the phylogenetic reconstructions. ACKNOWLEDGMENTS We are grateful to the curators and staff of the San Francisco Botanical Garden, the Japanese Tea Garden (San Francisco), and the University of California Botanical Garden (Berkeley, CA) for kindly allowing access to the collections. Holly Forbes, Barbara Keller, Mona Bourell, and Tony Morrosco provided much help with species sampling. We thank Cheng Huang and Frank Brink of the Center for Advanced Microscopy (Australian National University) for assistance with SEM. LITERATURE CITED Ackerly DD ( 2000 ) Taxon sampling, correlated evolution, and independent contrasts . Evolution 54 : 1480 – 1492 Google Scholar Crossref Search ADS PubMed WorldCat Alder NN Pockman WT Sperry JS Nuismer S ( 1997 ) Use of centrifugal force in the study of xylem cavitation . J Exp Bot 48 : 665 – 674 Google Scholar Crossref Search ADS WorldCat Bailey IW ( 1913 ) The preservative treatment of wood. II. The structure of the pit membranes in the tracheids of conifers and their relation to the penetration of gases, liquids and finely divided solids into green and seasoned wood . Forestry Quarterly 11 : 12 – 20 Google Scholar OpenURL Placeholder Text WorldCat Bauch JW Liese W Schultze R ( 1972 ) The morphological variability of the bordered pit membranes in gymnosperms . Wood Sci Technol 6 : 165 – 184 Google Scholar Crossref Search ADS WorldCat Brodribb T Cochard H ( 2009 ) Hydraulic failure defines the recovery and point of death in water-stressed conifers . Plant Physiol 149 : 575 – 584 Google Scholar Crossref Search ADS PubMed WorldCat Brodribb T Hill RS ( 1999 ) The importance of xylem constraints in the distribution of conifer species . New Phytol 143 : 364 – 372 Google Scholar Crossref Search ADS WorldCat Choat B Cobb AR Jansen S ( 2008 ) Structure and function of bordered pits: new discoveries and impacts on whole-plant hydraulic function . New Phytol 177 : 608 – 626 Google Scholar Crossref Search ADS PubMed WorldCat Choat B Pittermann J ( 2009 ) New insights into bordered pit structure and cavitation resistance in angiosperms and conifers . New Phytol 182 : 557 – 560 Google Scholar Crossref Search ADS PubMed WorldCat Choat B Sack L Holbrook NM ( 2007 ) Diversity of hydraulic traits in nine Cordia species growing in tropical forests with contrasting precipitation . New Phytol 175 : 686 – 698 Google Scholar Crossref Search ADS PubMed WorldCat Cochard H Holtta T Herbette S Delzon S Mencuccini M ( 2009 ) New insights into the mechanisms of water-stress-induced cavitation in conifers . Plant Physiol 151 : 949 – 954 Google Scholar Crossref Search ADS PubMed WorldCat Domec J Lachenbruch B Meinzer F ( 2006 ) Bordered pit structure and function determine spatial patterns of air-seeding thresholds in xylem of Douglas fir (Pseudotsuga menziesii; Pinaceae) trees . Am J Bot 93 : 1588 – 1600 Google Scholar Crossref Search ADS PubMed WorldCat Domec JC Lachenbruch B Meinzer F Woodruff D Warren JM McCulloh K ( 2008 ) Maximum height in a conifer is associated with conflicting requirements for xylem design . Proc Natl Acad Sci USA 105 : 12069 – 12074 Google Scholar Crossref Search ADS PubMed WorldCat EMBL-EBI ( 2009 ) Align sequences using ClustalW2 . http://www.ebi.ac.uk/Tools/clustalw2/index.html (January 28, 2010) Farjon A, editor ( 2005 ) A Monograph of Cupressaceae and Sciadopitys . Royal Botanic Gardens , Kew, UK Google Scholar Felsenstein J ( 1985 ) Phylogenies and the comparative method . Am Nat 125 : 1 – 15 Google Scholar Crossref Search ADS WorldCat Florin R ( 1963 ) The distribution of conifer and taxad genera in time and space . Acta Hort Bergiani 20 : 121 – 312 Google Scholar OpenURL Placeholder Text WorldCat Gadek PA Alpers DL Heslewood MM Quinn CJ ( 2000 ) Relationships between the Cupressaceae sensu lato: a combined morphological and molecular approach . Am J Bot 87 : 1044 – 1057 Google Scholar Crossref Search ADS PubMed WorldCat Garland T Harvey PH Ives AR ( 1992 ) Procedures for the analysis of comparative data using phylogenetically independent contrasts . Syst Biol 41 : 18 – 32 Google Scholar Crossref Search ADS WorldCat Hacke U Jansen S ( 2009 ) Embolism resistance of three boreal conifer species varies with pit structure . New Phytol 182 : 675 – 686 Google Scholar Crossref Search ADS PubMed WorldCat Hacke UG Sperry JS Pittermann J ( 2004 ) Analysis of circular bordered pit function. II. Gymnosperm tracheids with torus-margo pit membranes . Am J Bot 91 : 386 – 400 Google Scholar Crossref Search ADS PubMed WorldCat Hacke UG Stiller V Sperry JS Pittermann J McCulloh KA ( 2001 ) Cavitation fatigue: embolism and refilling cycles can weaken the cavitation resistance of xylem . Plant Physiol 125 : 779 – 786 Google Scholar Crossref Search ADS PubMed WorldCat Heady R Cunningham RB Evans PD ( 2008 ) Differences in the ultrastructure of their large warts allow white cypress pine (Callitris glaucophylla) to be distinguished from black cypress pine (C. endlicheri) . Wood Sci Technol 42 : 313 – 323 Google Scholar Crossref Search ADS WorldCat Karnovsky MJ ( 1965 ) A formaldehyde-glutaraldehyde fixative of high osmolality for use in electron microscopy . J Cell Biol 27 : 137 – 138 Google Scholar OpenURL Placeholder Text WorldCat Kembel SW Cowan PD Helmus MR Cornwell WK Morlon H Ackerly DD Blomberg SP Webb CO ( 2010 ) Picante: R tools for integrating phylogenies and ecology . Bioinformatics 26 : 1463 – 1464 Google Scholar Crossref Search ADS PubMed WorldCat Jacobsen A Pratt R Davis S Ewers F ( 2008 ) Comparative community physiology: nonconvergence in water relations among three semi-arid shrub communities . New Phytol 180 : 100 – 113 Google Scholar Crossref Search ADS PubMed WorldCat Jansen S Choat B Pletsers A ( 2009 ) Morphological variation of intervessel pit membranes and implications to xylem function in angiosperms . Am J Bot 96 : 409 – 419 Google Scholar Crossref Search ADS PubMed WorldCat Kohonen MM ( 2006 ) Engineered wettability in tree capillaries . Langmuir 22 : 3148 – 3153 Google Scholar Crossref Search ADS PubMed WorldCat Kohonen MM Helland A ( 2009 ) On the function of wall sculpturing in xylem conduits . J Bionics Eng 6 : 324 – 329 Google Scholar Crossref Search ADS WorldCat Kolb KJ Sperry JS ( 1999 ) Transport constraints on water use by the Great Basin shrub, Artemisia tridentata . Plant Cell Environ 22 : 925 – 935 Google Scholar Crossref Search ADS WorldCat Lancashire JR Ennos AR ( 2002 ) Modelling the hydrodynamic resistance of bordered pits . J Exp Bot 53 : 1485 – 1493 Google Scholar Crossref Search ADS PubMed WorldCat Liese W ( 1965 ) The fine structure of bordered pits in softwoods . In Cote WA , ed , Cellular Ultrastructure of Woody Plants . Syracuse University Press , Syracuse, NY , pp 270 – 289 Google Scholar Liese W Bauch J ( 1967 ) On the closure of bordered pits in conifers . Wood Sci Technol 1 : 1 – 13 Google Scholar Crossref Search ADS WorldCat Maddison WP Maddison DR ( 2009 ) Mesquite: a modular system for evolutionary analysis, version 2.71 . http://mesquiteproject.org (January 28, 2010) Maherali H Pockman WT Jackson RB ( 2004 ) Adaptive variation in the vulnerability of woody plants to xylem cavitation . Ecology 85 : 2184 – 2199 Google Scholar Crossref Search ADS WorldCat Mayr S Wolfschwenger M Bauer H ( 2002 ) Winter-drought induced embolism in Norway spruce (Picea abies) at the alpine timberline . Physiol Plant 115 : 74 – 80 Google Scholar Crossref Search ADS PubMed WorldCat Pammenter NW Vander Willigen C ( 1998 ) A mathematical and statistical analysis of the curves illustrating vulnerability of xylem to cavitation . Tree Physiol 18 : 589 – 593 Google Scholar Crossref Search ADS PubMed WorldCat Paradis E Claude J Strimmer K ( 2004 ) APE: analyses of phylogenetics and evolution in R language . Bioinformatics 20 : 289 – 290 Google Scholar Crossref Search ADS PubMed WorldCat Petty JA ( 1972 ) The aspiration of bordered pits in conifer wood . Proc R Soc Lond B Biol Sci 181 : 395 – 406 Google Scholar OpenURL Placeholder Text WorldCat Pittermann J Sperry JS ( 2003 ) Tracheid diameter is the key trait determining the extent of freezing-induced embolism in conifers . Tree Physiol 23 : 907 – 914 Google Scholar Crossref Search ADS PubMed WorldCat Pittermann J Sperry JS Hacke UG Wheeler JK Sikkema EH ( 2005 ) The torus-margo pit valve makes conifers hydraulically competitive with angiosperms . Science 310 : 1924 Google Scholar Crossref Search ADS PubMed WorldCat Pittermann J Sperry JS Hacke UG Wheeler JK Sikkema EH ( 2006a ) Mechanical reinforcement against tracheid implosion compromises the hydraulic efficiency of conifer xylem . Plant Cell Environ 29 : 1618 – 1628 Google Scholar Crossref Search ADS WorldCat Pittermann J Sperry JS Wheeler JK Hacke UG Sikkema EH ( 2006b ) Inter-tracheid pitting and the hydraulic efficiency of conifer wood: the role of tracheid allometry and cavitation protection . Am J Bot 93 : 1265 – 1273 Google Scholar Crossref Search ADS WorldCat Pockman WT Sperry JS ( 2000 ) Vulnerability to xylem cavitation and the distribution of Sonoran Desert vegetation . Am J Bot 87 : 1287 – 1299 Google Scholar Crossref Search ADS PubMed WorldCat Pockman WT Sperry JS O'Leary JW ( 1995 ) Sustained and significant negative water pressure in xylem . Nature 378 : 715 – 716 Google Scholar Crossref Search ADS WorldCat Pratt RB Jacobsen AL Ewers FW Davis SD ( 2007 ) Relationships among xylem transport, biomechanics and storage in stems and roots of nine Rhamnaceae species of the California chaparral . New Phytol 174 : 787 – 798 Google Scholar Crossref Search ADS PubMed WorldCat Preston KA Ackerly DD ( 2003 ) Hydraulic architecture and the evolution of shoot allometry in contrasting climates . Am J Bot 90 : 1502 – 1511 Google Scholar Crossref Search ADS PubMed WorldCat R Development Core Team ( 2008 ) R: A Language and Environment for Statistical Computing . R Foundation for Statistical Computing , Vienna Google Scholar Sanderson MJ ( 2002 ) Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach . Mol Biol Evol 19 : 101 – 109 Google Scholar Crossref Search ADS PubMed WorldCat Santiago LS Goldstein G Meinzer FC Fisher JB Machado K Woodruff D Jones T ( 2004 ) Leaf photosynthetic traits scale with hydraulic conductivity and wood density in Panamanian forest canopy trees . Oecologia 140 : 543 – 550 Google Scholar Crossref Search ADS PubMed WorldCat Schoonmaker AL Hacke UG Landhausser SM Lieffers VJ Tyree MT ( 2010 ) Hydraulic acclimation to shading in boreal conifers of varying shade tolerance . Plant Cell Environ 33 : 382 – 393 Google Scholar Crossref Search ADS PubMed WorldCat Sperry JS ( 1993 ) Winter xylem embolism and spring recovery in Betula cordifolia, Fagus grandifolia, Abies balsamea, and Picea rubens . In Borghetti M Grace J Raschi A , eds , Water Transport in Plants under Climatic Stress . Cambridge University Press , Cambridge, UK , pp 86 – 98 Google Scholar Crossref Search ADS Sperry JS Hacke UG ( 2004 ) Analysis of circular bordered pit function. I. Angiosperm vessels with homogenous pit membranes . Am J Bot 91 : 369 – 385 Google Scholar Crossref Search ADS PubMed WorldCat Sperry JS Hacke UG Pittermann J ( 2006 ) Size and function in conifer tracheids and angiosperm vessels . Am J Bot 93 : 1490 – 1500 Google Scholar Crossref Search ADS PubMed WorldCat Sperry JS Nichols KL Sullivan JEM Eastlack SE ( 1994 ) Xylem embolism in ring-porous, diffuse-porous, and coniferous trees of northern Utah and interior Alaska . Ecology 75 : 1736 – 1752 Google Scholar Crossref Search ADS WorldCat Sperry JS Tyree MT ( 1990 ) Water-stress-induced xylem embolism in three species of conifers . Plant Cell Environ 13 : 427 – 436 Google Scholar Crossref Search ADS WorldCat Stockey RA Kvacek J Hill R Rothwell GH Kvacek Z ( 2005 ) The fossil record of the Cupressaceae . In Farjon A , ed , A Monograph of the Cupressaceae and Sciadopitys . Royal Botanic Gardens , Kew, UK , pp 54 – 68 Google Scholar Swofford DL ( 1998 ) PAUP*: Phylogenetic Analysis Using Parsimony (*and Other Methods), 4.0 Beta 5 . Sinauer Associates , Sunderland, MA Google Scholar Taylor TN Taylor EL Krings M ( 2009 ) Paleobotany: The Biology and Evolution of Fossil Plants . Academic Press , Oxford Google Scholar Thompson JD Higgins DG Gibson TJ ( 1994 ) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice . Nucleic Acids Res 22 : 4673 – 4680 Google Scholar Crossref Search ADS PubMed WorldCat Valli A Koponen A Vesala T Timonen S ( 2002 ) Simulations of water flow through bordered pits of conifer xylem . J Stat Phys 107 : 121 – 142 Google Scholar Crossref Search ADS WorldCat Willis KJ McElwain JC ( 2002 ) The Evolution of Plants . Oxford University Press , Oxford Google Scholar Willson CJ Manos PS Jackson RB ( 2008 ) Hydraulic traits are influenced by phylogenetic history in the drought-resistant, invasive genus Juniperus (Cupressaceae) . Am J Bot 95 : 299 – 314 Google Scholar Crossref Search ADS PubMed WorldCat Author notes 1 This work was supported by the Miller Institute for Basic Research at the University of California, Berkeley (to J.P. and T.D.), and the Ministry of Science, Research, and the Arts from Baden-Wurttemberg (to S.J.). 2 These authors contributed equally to the article. * Corresponding author; e-mail pittermann@biology.ucsc.edu. The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantphysiol.org) is: Jarmila Pittermann (pittermann@biology.ucsc.edu). [W] The online version of this article contains Web-only data. [OA] Open Access articles can be viewed online without a subscription. www.plantphysiol.org/cgi/doi/10.1104/pp.110.158824 © 2010 American Society of Plant Biologists © The Author(s) 2010. Published by Oxford University Press on behalf of American Society of Plant Biologists. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. TI - The Relationships between Xylem Safety and Hydraulic Efficiency in the Cupressaceae: The Evolution of Pit Membrane Form and Function     JF - Plant Physiology DO - 10.1104/pp.110.158824 DA - 2010-08-03 UR - https://www.deepdyve.com/lp/oxford-university-press/the-relationships-between-xylem-safety-and-hydraulic-efficiency-in-the-dPSqvMnoHy SP - 1919 EP - 1931 VL - 153 IS - 4 DP - DeepDyve ER -