Sketching algorithms for genomic data analysis and querying in a secure enclave.

Genome-wide association studies (GWAS), especially on rare diseases, may necessitate exchange of sensitive genomic data between multiple institutions. Since genomic data sharing is often infeasible due to privacy concerns, cryptographic methods, such as secure multiparty computation (SMC) protocols, have been developed with the aim of offering privacy-preserving collaborative GWAS. Unfortunately, the computational overhead of these methods remain prohibitive for human-genome-scale data. Here we introduce SkSES (https://github.com/ndokmai/sgx-genome-variants-search), a hardware-software hybrid approach for privacy-preserving collaborative GWAS, which improves the running time of the most advanced cryptographic protocols by two orders of magnitude. The SkSES approach is based on trusted execution environments (TEEs) offered by current-generation microprocessors-in particular, Intel's SGX. To overcome the severe memory limitation of the TEEs, SkSES employs novel 'sketching' algorithms that maintain essential statistical information on genomic variants in input VCF files. By additionally incorporating efficient data compression and population stratification reduction methods, SkSES identifies the top k genomic variants in a cohort quickly, accurately and in a privacy-preserving manner.

The complexity of biological disturbance agents, fuels heterogeneity, and fire in coniferous forests of the western United States.

Forest biological disturbance agents (BDAs) are insects, pathogens, and parasitic plants that affect tree decline, mortality, and forest ecosystems processes. BDAs are commonly thought to increase the likelihood and severity of fire by converting live standing trees to more flammable, dead and downed fuel. However, recent research indicates that BDAs do not necessarily increase, and can reduce, the likelihood or severity of fire. This has led to confusion regarding the role of BDAs in influencing fuels and fire in fire-prone western United States forests. Here, we review the existing literature on BDAs and their effects on fuels and fire in the western US and develop a conceptual framework to better understand the complex relationships between BDAs, fuels and fire. We ask: 1) What are the major BDA groups in western US forests that affect fuels? and 2) How do BDA-affected fuels influence fire risk and outcomes? The conceptual framework is rooted in the spatiotemporal aspects of BDA life histories, which drive forest impacts, fuel characteristics and if ignited, fire outcomes. Life histories vary among BDAs from episodic, landscape-scale outbreaks (bark beetles, defoliators), to chronic, localized disturbance effects (dwarf mistletoes, root rots). Generally, BDAs convert aboveground live biomass to dead biomass, decreasing canopy fuels and increasing surface fuels. However, the rate of conversion varies with time-since-event and among BDAs and forest types, resulting in a wide range of effects on the amount of dead fuels at any given time and place, which interacts with the structure and composition of the stand before and subsequent to BDA events. A major influence on fuels may be that BDAs have emerged as dominant agents of forest heterogeneity creation. Because BDAs play complex roles in fuels and fire heterogeneity across the western US which are further complicated by interactions with climate change, drought, and forest management (fire suppression), their impacts on fuels, fire and ecological consequences cannot be categorized simply as positive or negative but need to be evaluated within the context of BDA life histories and ecosystem dynamics.

How to reduce dimension with PCA and random projections?

In our "big data" age, the size and complexity of data is steadily increasing. Methods for dimension reduction are ever more popular and useful. Two distinct types of dimension reduction are "data-oblivious" methods such as random projections and sketching, and "data-aware" methods such as principal component analysis (PCA). Both have their strengths, such as speed for random projections, and data-adaptivity for PCA. In this work, we study how to combine them to get the best of both. We study "sketch and solve" methods that take a random projection (or sketch) first, and compute PCA after. We compute the performance of several popular sketching methods (random iid projections, random sampling, subsampled Hadamard transform, CountSketch, etc) in a general "signal-plus-noise" (or spiked) data model. Compared to well-known works, our results (1) give asymptotically exact results, and (2) apply when the signal components are only slightly above the noise, but the projection dimension is non-negligible. We also study stronger signals allowing more general covariance structures. We find that (a) signal strength decreases under projection in a delicate way depending on the structure of the data and the sketching method, (b) orthogonal projections are slightly more accurate, (c) randomization does not hurt too much, due to concentration of measure, (d) CountSketch can be somewhat improved by a normalization method. Our results have implications for statistical learning and data analysis. We also illustrate that the results are highly accurate in simulations and in analyzing empirical data.

Water potential control of turgor-driven tracheid enlargement in Scots pine at its xeric distribution edge.

The extent to which water availability can be used to predict the enlargement and final dimensions of xylem conduits remains an open issue. We reconstructed the time course of tracheid enlargement in Pinus sylvestris trees in central Spain by repeated measurements of tracheid diameter on microcores sampled weekly during a 2 yr period. We analyzed the role of water availability in these dynamics empirically through time-series correlation analysis and mechanistically by building a model that simulates daily tracheid enlargement rate and duration based on Lockhart's equation and water potential as the sole input. Tracheid enlargement followed a sigmoid-like time course, which varied intra- and interannually. Our empirical analysis showed that final tracheid diameter was strongly related to water availability during tracheid enlargement. The mechanistic model was calibrated and successfully validated (R2  = 0.92) against the observed tracheid enlargement time course. The model was also able to reproduce the seasonal variations of tracheid enlargement rate, duration and final diameter (R2  = 0.84-0.99). Our results support the hypothesis that tracheid enlargement and final dimensions can be modeled based on the direct effect of water potential on turgor-driven cell expansion. We argue that such a mechanism is consistent with other reported patterns of tracheid dimension variation.

Coordination and trade-offs between leaf and stem hydraulic traits and stomatal regulation along a spectrum of isohydry to anisohydry.

The degree of plant iso/anisohydry, a widely used framework for classifying species-specific hydraulic strategies, integrates multiple components of the whole-plant hydraulic pathway. However, little is known about how it associates with coordination of functional and structural traits within and across different organs. We examined stem and leaf hydraulic capacitance and conductivity/conductance, stem xylem anatomical features, stomatal regulation of daily minimum leaf and stem water potential (Ψ), and the kinetics of stomatal responses to vapour pressure deficit (VPD) in six diverse woody species differing markedly in their degree of iso/anisohydry. At the stem level, more anisohydric species had higher wood density and lower native capacitance and conductivity. Like stems, leaves of more anisohydric species had lower hydraulic conductance; however, unlike stems, their leaves had higher native capacitance at their daily minimum values of leaf Ψ. Moreover, rates of VPD-induced stomatal closure were related to intrinsic rather than native leaf capacitance and were not associated with species' degree of iso/anisohydry. Our results suggest a trade-off between hydraulic storage and efficiency in the leaf, but a coordination between hydraulic storage and efficiency in the stem along a spectrum of plant iso/anisohydry.

Co-occurring woody species have diverse hydraulic strategies and mortality rates during an extreme drought.

From 2011 to 2013, Texas experienced its worst drought in recorded history. This event provided a unique natural experiment to assess species-specific responses to extreme drought and mortality of four co-occurring woody species: Quercus fusiformis, Diospyros texana, Prosopis glandulosa, and Juniperus ashei. We examined hypothesized mechanisms that could promote these species' diverse mortality patterns using postdrought measurements on surviving trees coupled to retrospective process modelling. The species exhibited a wide range of gas exchange responses, hydraulic strategies, and mortality rates. Multiple proposed indices of mortality mechanisms were inconsistent with the observed mortality patterns across species, including measures of the degree of iso/anisohydry, photosynthesis, carbohydrate depletion, and hydraulic safety margins. Large losses of spring and summer whole-tree conductance (driven by belowground losses of conductance) and shallower rooting depths were associated with species that exhibited greater mortality. Based on this retrospective analysis, we suggest that species more vulnerable to drought were more likely to have succumbed to hydraulic failure belowground.

Interactions of predominant insects and diseases with climate change in Douglas-fir forests of western Oregon and Washington, U.S.A.

Forest disturbance regimes are beginning to show evidence of climate-mediated changes, such as increasing severity of droughts and insect outbreaks. We review the major insects and pathogens affecting the disturbance regime for coastal Douglas-fir forests in western Oregon and Washington State, USA, and ask how future climate changes may influence their role in disturbance ecology. Although the physiological constraints of light, temperature, and moisture largely control tree growth, episodic and chronic disturbances interacting with biological factors have substantial impacts on the structure and functioning of forest ecosystems in this region. Understanding insect and disease interactions is critical to predicting forest response to climate change and the consequences for ecosystem services, such as timber, clean water, fish and wildlife. We focused on future predictions for warmer wetter winters, hotter drier summers, and elevated atmospheric CO2 to hypothesize the response of Douglas-fir forests to the major insects and diseases influencing this forest type: Douglas-fir beetle, Swiss needle cast, black stain root disease, and laminated root rot. We hypothesize that 1) Douglas-fir beetle and black stain root disease could become more prevalent with increasing, fire, temperature stress, and moisture stress, 2) future impacts of Swiss needle cast are difficult to predict due to uncertainties in May-July leaf wetness, but warmer winters could contribute to intensification at higher elevations, and 3) laminated root rot will be influenced primarily by forest management, rather than climatic change. Furthermore, these biotic disturbance agents interact in complex ways that are poorly understood. Consequently, to inform management decisions, insect and disease influences on disturbance regimes must be characterized specifically by forest type and region in order to accurately capture these interactions in light of future climate-mediated changes.

Stomatal kinetics and photosynthetic gas exchange along a continuum of isohydric to anisohydric regulation of plant water status.

Species' differences in the stringency of stomatal control of plant water potential represent a continuum of isohydric to anisohydric behaviours. However, little is known about how quasi-steady-state stomatal regulation of water potential may relate to dynamic behaviour of stomata and photosynthetic gas exchange in species operating at different positions along this continuum. Here, we evaluated kinetics of light-induced stomatal opening, activation of photosynthesis and features of quasi-steady-state photosynthetic gas exchange in 10 woody species selected to represent different degrees of anisohydry. Based on a previously developed proxy for the degree of anisohydry, species' leaf water potentials at turgor loss, we found consistent trends in photosynthetic gas exchange traits across a spectrum of isohydry to anisohydry. More anisohydric species had faster kinetics of stomatal opening and activation of photosynthesis, and these kinetics were closely coordinated within species. Quasi-steady-state stomatal conductance and measures of photosynthetic capacity and performance were also greater in more anisohydric species. Intrinsic water-use efficiency estimated from leaf gas exchange and stable carbon isotope ratios was lowest in the most anisohydric species. In comparisons between gas exchange traits, species rankings were highly consistent, leading to species-independent scaling relationships over the range of isohydry to anisohydry observed.

Mapping 'hydroscapes' along the iso- to anisohydric continuum of stomatal regulation of plant water status.

The concept of iso- vs. anisohydry has been used to describe the stringency of stomatal regulation of plant water potential (ψ). However, metrics that accurately and consistently quantify species' operating ranges along a continuum of iso- to anisohydry have been elusive. Additionally, most approaches to quantifying iso/anisohydry require labour-intensive measurements during prolonged drought. We evaluated new and previously developed metrics of stringency of stomatal regulation of ψ during soil drying in eight woody species and determined whether easily-determined leaf pressure-volume traits could serve as proxies for their degree of iso- vs. anisohydry. Two metrics of stringency of stomatal control of ψ, (1) a 'hydroscape' incorporating the landscape of ψ over which stomata control ψ, and (2) the slope of the daily range of ψ as pre-dawn ψ declined, were strongly correlated with each other and with the leaf osmotic potential at full and zero turgor derived from pressure-volume curves.

Linking nonstructural carbohydrate dynamics to gas exchange and leaf hydraulic behavior in Pinus edulis and Juniperus monosperma.

Leaf hydraulics, gas exchange and carbon storage in Pinus edulis and Juniperus monosperma, two tree species on opposite ends of the isohydry-anisohydry spectrum, were analyzed to examine relationships between hydraulic function and carbohydrate dynamics. Leaf hydraulic vulnerability, leaf water potential (Ψl ), leaf hydraulic conductance (Kleaf ), photosynthesis (A), stomatal conductance (gs) and nonstructural carbohydrate (NSC) content were analyzed throughout the growing season. Leaf hydraulic vulnerability was significantly lower in the relatively anisohydric J. monosperma than in the more isohydric P. edulis. In P. edulis, Ψl dropped and stayed below 50% loss of leaf hydraulic conductance (P50) early in the day during May, August and around midday in September, leading to sustained reductions in Kleaf . In J. monosperma, Ψl dropped below P50 only during August, resulting in the maintenance of Kleaf during much of the growing season. Mean A and gs during September were significantly lower in P. edulis than in J. monosperma. Foliar total NSC was two to three times greater in J. monosperma than in P. edulis in June, August and September. Consistently lower levels of total NSC in P. edulis suggest that its isohydric strategy pushes it towards the exhaustion of carbon reserves during much of the growing season.

The dynamic pipeline: hydraulic capacitance and xylem hydraulic safety in four tall conifer species.

Recent work has suggested that plants differ in their relative reliance on structural avoidance of embolism versus maintenance of the xylem water column through dynamic traits such as capacitance, but we still know little about how and why species differ along this continuum. It is even less clear how or if different parts of a plant vary along this spectrum. Here we examined how traits such as hydraulic conductivity or conductance, xylem vulnerability curves, and capacitance differ in trunks, large- and small-diameter branches, and foliated shoots of four species of co-occurring conifers. We found striking similarities among species in most traits, but large differences among plant parts. Vulnerability to embolism was high in shoots, low in small- and large-diameter branches, and high again in the trunks. Safety margins, defined as the pressure causing 50% loss of hydraulic conductivity or conductance minus the midday water potential, were large in small-diameter branches, small in trunks and negative in shoots. Sapwood capacitance increased with stem diameter, and was correlated with stem vulnerability, wood density and latewood proportion. Capacitive release of water is a dynamic aspect of plant hydraulics that is integral to maintenance of long-distance water transport.

Dynamics of leaf water relations components in co-occurring iso- and anisohydric conifer species.

Because iso- and anisohydric species differ in stomatal regulation of the rate and magnitude of fluctuations in shoot water potential, they may be expected to show differences in the plasticity of their shoot water relations components, but explicit comparisons of this nature have rarely been made. We subjected excised shoots of co-occurring anisohydric Juniperus monosperma and isohydric Pinus edulis to pressure-volume analysis with and without prior artificial rehydration. In J. monosperma, the shoot water potential at turgor loss (Ψ(TLP)) ranged from -3.4 MPa in artificially rehydrated shoots to -6.6 MPa in shoots with an initial Ψ of -5.5 MPa, whereas in P. edulis mean Ψ(TLP) remained at ∼ -3.0 MPa over a range of initial Ψ from -0.1 to -2.3 MPa. The shoot osmotic potential at full turgor and the bulk modulus of elasticity also declined sharply with shoot Ψ in J. monosperma, but not in P. edulis. The contrasting behaviour of J. monosperma and P. edulis reflects differences in their capacity for homeostatic regulation of turgor that may be representative of aniso- and isohydric species in general, and may also be associated with the greater capacity of J. monosperma to withstand severe drought.

Learning-Augmented Sketching Offers Improved Performance for Privacy Preserving and Secure GWAS.

The introduction of trusted execution environments (TEEs), such as secure enclaves provided by the Intel SGX technology has enabled secure and privacy-preserving computation on the cloud. The stringent resource limitations, such as memory constraints, required by some TEEs necessitates the development of computational approaches with reduced memory usage, such as sketching. One example is the SkSES method for GWAS on a cohort of case and control samples from multiple institutions, which identifies the most significant SNPs in a privacy-preserving manner without disclosing sensitive genotype information to other institutions or the cloud service provider. Here we show how to improve the performance of SkSES on large datasets by augmenting it with a learning-augmented approach. Specifically, we show how individual institutions can perform smaller scale GWAS on their own datasets and identify two sets of variants according to certain criteria, which are then used to guide the sketching process to more accurately identify significant variants over the collective dataset. The new method achieves up to 40% accuracy gain compared to the original SkSES method under the same memory constraints on datasets we tested on. The code is available at https://github.com/alreadydone/sgx-genome-variants-search.

Contrasting hydraulic strategies in two tropical lianas and their host trees.

PREMISE OF THE STUDY: Tropical liana abundance has been increasing over the past 40 yr, which has been associated with reduced rainfall. The proposed mechanism allowing lianas to thrive in dry conditions is deeper root systems than co-occurring trees, although we know very little about the fundamental hydraulic physiology of lianas. METHODS: To test the hypothesis that two abundant liana species would physiologically outperform their host tree under reduced water availability, we measured rooting depth, hydraulic properties, plant water status, and leaf gas exchange during the dry season in a seasonally dry tropical forest. We also used a model to compare water use by one of the liana species and the host tree during drought. KEY RESULTS: All species measured were shallowly rooted. The liana species were more vulnerable to embolism than host trees and experienced water potentials that were predicted to result in substantial hydraulic losses in both leaves and stems. Water potentials measured in host trees were not negative enough to result in significant hydraulic losses. Model results predicted the liana to have greater gas exchange than its host tree during drought and nondrought conditions. CONCLUSIONS: The host tree species had a more conservative strategy for maintenance of the soil-to-leaf hydraulic pathway than the lianas it supported. The two liana species experienced embolism in stems and leaves, based on vulnerability curves and water potentials. These emboli were presumably repaired before the next morning. However, in the host tree species, reduced stomatal conductance prevented leaf or stem embolism.

Evidence for xylem embolism as a primary factor in dehydration-induced declines in leaf hydraulic conductance.

Hydraulic conductance of leaves (K(leaf)) typically decreases with increasing water stress and recent studies have proposed different mechanisms responsible for decreasing K(leaf) . We measured K(leaf) concurrently with ultrasonic acoustic emissions (UAEs) in dehydrating leaves of several species to determine whether declining K(leaf) was associated with xylem embolism. In addition, we performed experiments in which the surface tension of water in the leaf xylem was reduced by using a surfactant solution. Finally, we compared the hydraulic vulnerability of entire leaves with the leaf lamina in three species. Leaf hydraulic vulnerability based on rehydration kinetics and UAE was very similar, except in Quercus garryana. However, water potentials corresponding to the initial decline in K(leaf) and the onset of UAE in Q. garryana were similar. In all species tested, reducing the surface tension of water caused K(leaf) to decline at less negative water potentials compared with leaves supplied with water. Microscopy revealed that as the fraction of embolized xylem increased, K(leaf) declined sharply in Q. garryana. Measurements on leaf discs revealed that reductions in lamina hydraulic conductance with dehydration were not as great as those observed in intact leaves, suggesting that embolism was the primary mechanism for reductions in K(leaf) during dehydration.

A Fast, Provably Accurate Approximation Algorithm for Sparse Principal Component Analysis Reveals Human Genetic Variation Across the World.

Principal component analysis (PCA) is a widely used dimensionality reduction technique in machine learning and multivariate statistics. To improve the interpretability of PCA, various approaches to obtain sparse principal direction loadings have been proposed, which are termed Sparse Principal Component Analysis (SPCA). In this paper, we present ThreSPCA, a provably accurate algorithm based on thresholding the Singular Value Decomposition for the SPCA problem, without imposing any restrictive assumptions on the input covariance matrix. Our thresholding algorithm is conceptually simple; much faster than current state-of-the-art; and performs well in practice. When applied to genotype data from the 1000 Genomes Project, ThreSPCA is faster than previous benchmarks, at least as accurate, and leads to a set of interpretable biomarkers, revealing genetic diversity across the world.

Thermodynamic Driving Forces for Substrate Atom Extraction by Adsorption of Strong Electron Acceptor Molecules.

A quantitative structural investigation is reported, aimed at resolving the issue of whether substrate adatoms are incorporated into the monolayers formed by strong molecular electron acceptors deposited onto metallic electrodes. A combination of normal-incidence X-ray standing waves, low-energy electron diffraction, scanning tunnelling microscopy, and X-ray photoelectron spectroscopy measurements demonstrate that the systems TCNQ and F4TCNQ on Ag(100) lie at the boundary between these two possibilities and thus represent ideal model systems with which to study this effect. A room-temperature commensurate phase of adsorbed TCNQ is found not to involve Ag adatoms, but to adopt an inverted bowl configuration, long predicted but not previously identified experimentally. By contrast, a similar phase of adsorbed F4TCNQ does lead to Ag adatom incorporation in the overlayer, the cyano end groups of the molecule being twisted relative to the planar quinoid ring. Density functional theory (DFT) calculations show that this behavior is consistent with the adsorption energetics. Annealing of the commensurate TCNQ overlayer phase leads to an incommensurate phase that does appear to incorporate Ag adatoms. Our results indicate that the inclusion (or exclusion) of metal atoms into the organic monolayers is the result of both thermodynamic and kinetic factors.

Cross-biome synthesis of source versus sink limits to tree growth.

Uncertainties surrounding tree carbon allocation to growth are a major limitation to projections of forest carbon sequestration and response to climate change. The prevalence and extent to which carbon assimilation (source) or cambial activity (sink) mediate wood production are fundamentally important and remain elusive. We quantified source-sink relations across biomes by combining eddy-covariance gross primary production with extensive on-site and regional tree ring observations. We found widespread temporal decoupling between carbon assimilation and tree growth, underpinned by contrasting climatic sensitivities of these two processes. Substantial differences in assimilation-growth decoupling between angiosperms and gymnosperms were determined, as well as stronger decoupling with canopy closure, aridity, and decreasing temperatures. Our results reveal pervasive sink control over tree growth that is likely to be increasingly prominent under global climate change.

Water stress, shoot growth and storage of non-structural carbohydrates along a tree height gradient in a tall conifer.

We analysed concentrations of starch, sucrose, glucose and fructose in upper branch wood, foliage and trunk sapwood of Douglas-fir trees in height classes ranging from ~2 to ~57 m. Mean concentrations of non-structural carbohydrates (NSC) for all tissues were highest in the tallest height class and lowest in the lowest height class, and height-related trends in NSC were most pronounced in branches. Throughout a 17-month sampling period, mean values of branch NSC from the 57 m trees ranged between 30 and 377% greater than the 2 m trees. Branch NSC was inversely correlated with midday shoot water potential (Ψ(l)), shoot osmotic potential at full turgor (Ψ) and shoot extension. Temporal fluctuation in branch NSC was inversely correlated with height, and positively correlated with midday Ψ(l) , Ψ and shoot extension. The positive correlation between height and storage of NSC, and the negative correlation between NSC storage and shoot extension provide evidence that size-related growth decline in trees is not strongly associated with constraints on photosynthesis. The negative correlation between height and fluctuation in NSC suggests that mobilization of photosynthate in taller trees is constrained by some factor such as reductions in turgor-driven cell expansion or constraints on phloem transport.

Climate-related trends in sapwood biophysical properties in two conifers: avoidance of hydraulic dysfunction through coordinated adjustments in xylem efficiency, safety and capacitance.

In the Pacific north-west, the Cascade Mountain Range blocks much of the precipitation and maritime influence of the Pacific Ocean, resulting in distinct climates east and west of the mountains. The current study aimed to investigate relationships between water storage and transport properties in populations of Douglas-fir (Pseudotsuga menziesii) and ponderosa pine (Pinus ponderosa) adapted to both climates. Sapwood thickness, capacitance, vulnerability to embolism, and axial and radial conductivity were measured on samples collected from trunks of mature trees. The sapwood of ponderosa pine was three to four times thicker than Douglas-fir. Radial conductivity was higher in west-side populations of both species, but axial conductivity was higher in the east-side populations and in Douglas-fir. Eastern populations of both species had sapwood that was more vulnerable to embolism than west-side populations. Sapwood capacitance was similar between species, but was about twice as great in east-side populations (580 kg m⁻³ MPa⁻¹) as in west-side populations (274 kg m⁻³ MPa⁻¹). Capacitance was positively correlated with both mean embolism pressure and axial conductivity across species and populations, suggesting that coordinated adjustments in xylem efficiency, safety and water storage capacity may serve to avoid embolism along a gradient of increasing aridity.