FLUXNET and modelling the global carbon cycle

Measurements of the net CO₂ flux between terrestrial ecosystems and the atmosphere using the eddy covariance technique have the potential to underpin our interpretation of regional CO₂ source–sink patterns, CO₂ flux responses to forcings, and predictions of the future terrestrial C balance. Information contained in FLUXNET eddy covariance data has multiple uses for the development and application of global carbon models, including evaluation/validation, calibration, process parameterization, and data assimilation. This paper reviews examples of these uses, compares global estimates of the dynamics of the global carbon cycle, and suggests ways of improving the utility of such data for global carbon modelling. Net ecosystem exchange of CO₂ (NEE) predicted by different terrestrial biosphere models compares favourably with FLUXNET observations at diurnal and seasonal timescales. However, complete model validation, particularly over the full annual cycle, requires information on the balance between assimilation and decomposition processes, information not readily available for most FLUXNET sites. Site history, when known, can greatly help constrain the model‐data comparison. Flux measurements made over four vegetation types were used to calibrate the land‐surface scheme of the Goddard Institute for Space Studies global climate model, significantly improving simulated climate and demonstrating the utility of diurnal FLUXNET data for climate modelling. Land‐surface temperatures in many regions cool due to higher canopy conductances and latent heat fluxes, and the spatial distribution of CO₂ uptake provides a significant additional constraint on the realism of simulated surface fluxes. FLUXNET data are used to calibrate a global production efficiency model (PEM). This model is forced by satellite‐measured absorbed radiation and suggests that global net primary production (NPP) increased 6.2% over 1982–1999. Good agreement is found between global trends in NPP estimated by the PEM and a dynamic global vegetation model (DGVM), and between the DGVM and estimates of global NEE derived from a global inversion of atmospheric CO₂ measurements. Combining the PEM, DGVM, and inversion results suggests that CO₂ fertilization is playing a major role in current increases in NPP, with lesser impacts from increasing N deposition and growing season length. Both the PEM and the inversion identify the Amazon basin as a key region for the current net terrestrial CO₂ uptake (i.e. 33% of global NEE), as well as its interannual variability. The inversion's global NEE estimate of −1.2 Pg [C] yr⁻¹ for 1982–1995 is compatible with the PEM‐ and DGVM‐predicted trends in NPP. There is, thus, a convergence in understanding derived from process‐based models, remote‐sensing‐based observations, and inversion of atmospheric data. Future advances in field measurement techniques, including eddy covariance (particularly concerning the problem of night‐time fluxes in dense canopies and of advection or flow distortion over complex terrain), will result in improved constraints on land‐atmosphere CO₂ fluxes and the rigorous attribution of mechanisms to the current terrestrial net CO₂ uptake and its spatial and temporal heterogeneity. Global ecosystem models play a fundamental role in linking information derived from FLUXNET measurements to atmospheric CO₂ variability. A number of recommendations concerning FLUXNET data are made, including a request for more comprehensive site data (particularly historical information), more measurements in undisturbed ecosystems, and the systematic provision of error estimates. The greatest value of current FLUXNET data for global carbon cycle modelling is in evaluating process representations, rather than in providing an unbiased estimate of net CO₂ exchange.     

Bootstrapping phylogenies inferred from rearrangement data

ABSTRACT: Large-scale sequencing of genomes has enabled the inference of phylogenies based on the evolution of genomic architecture, under such events as rearrangements, duplications, and losses. Many evolutionary models and associated algorithms have been designed over the last few years and have found use in comparative genomics and phylogenetic inference. However, the assessment of phylogenies built from such data has not been properly addressed to date. The standard method used in sequence-based phylogenetic inference is the bootstrap, but it relies on a large number of homologous characters that can be resampled; yet in the case of rearrangements, the entire genome is a single character. Alternatives such as the jackknife suffer from the same problem, while likelihood tests cannot be applied in the absence of well established probabilistic models. We present a new approach to the assessment of distance-based phylogenetic inference from whole-genome data; our approach combines features of the jackknife and the bootstrap and remains nonparametric. For each feature of our method, we give an equivalent feature in the sequence-based framework; we also present the results of extensive experimental testing, in both sequence-based and genome-based frameworks. Through the feature-by-feature comparison and the experimental results, we show that our bootstrapping approach is on par with the classic phylogenetic bootstrap used in sequence-based reconstruction, and we establish the clear superiority of the classic bootstrap for sequence data and of our corresponding new approach for rearrangement data over proposed variants. Finally, we test our approach on a small dataset of mammalian genomes, verifying that the support values match current thinking about the respective branches. Our method is the first to provide a standard of assessment to match that of the classic phylogenetic bootstrap for aligned sequences. Its support values follow a similar scale and its receiver-operating characteristics are nearly identical, indicating that it provides similar levels of sensitivity and specificity. Thus our assessment method makes it possible to conduct phylogenetic analyses on whole genomes with the same degree of confidence as for analyses on aligned sequences. Extensions to search-based inference methods such as maximum parsimony and maximum likelihood are possible, but remain to be thoroughly tested.     

Detection of drug binding to DNA by hydroxyl radical footprinting. Relationship of distamycin binding sites to DNA structure and positioned nucleosomes on 5S RNA genes of Xenopus.

We report the use of hydroxyl radical footprinting to analyze the interaction of distamycin and actinomycin with the 5S ribosomal RNA genes of Xenopus. There is a qualitative difference in the hydroxyl radical footprints of the two drugs. Distamycin gives a conventional (albeit high-resolution) footprint, while actinomycin does not protect DNA from hydroxyl radical attack, but instead induces discrete sites of hyperreactivity. We find concentration-dependent changes in the locations of distamycin binding sites on the somatic 5S gene of Xenopus borealis. A high-affinity site, containing a G.C base pair, is replaced at higher levels of bound drug by a periodic array of different lower affinity sites that coincide with the places where the minor groove of the DNA would face in toward a nucleosome core that is known to bind to the same sequence. These results suggest that distamycin recognizes potential binding sites more by the shape of the DNA than by the specific sequence that is contained in the site and that structures of many sequences are deformable to a shape that allows drug binding. We discuss the utility of hydroxyl radical footprinting of distamycin for investigating the underlying structure of DNA.     

Effect of naloxone on gonadotrophin secretion at various stages of development in the ewe lamb

Spring-born crossbred ewe lambs were raised in a natural photoperiod and saline (N = 6) or naloxone (1 mg/kg) in saline (N = 6) was injected (i.m.) every 2 h for 6 h at 5, 10 and 15 weeks of age and for 8 h at 20, 25 and 30 weeks of age. Blood samples were taken every 12 min during treatment periods. Naloxone had no effect on time to first oestrus (controls 235 +/- 6 days, naloxone 242 +/- 7 days). Mean serum LH concentrations and LH pulse frequency were elevated by naloxone in ewe lambs at 20, 25, and 30 weeks of age (P less than 0.05). The only FSH response to naloxone was a depression of mean serum concentrations at 30 weeks of age (P less than 0.05). LH pulse amplitude was elevated at 5 weeks of age in all ewe lambs and declined thereafter to a nadir at 30 weeks of age in control, but not in naloxone-treated animals (P less than 0.05). LH pulse frequency was elevated at 10 weeks of age in control ewe lambs and in all animals at 30 weeks of age (P less than 0.05). FSH pulse frequency declined from 5 weeks of age in control ewe lambs (P less than 0.05), with very few pulses noted in 25- and 30-week-old animals. We conclude that (1) opioidergic suppression of LH, but not FSH, secretion developed at 20 weeks of age in the growing ewe lambs used in the present study, with no obvious change in suppression before the onset of first oestrus: (2) pulsatile FSH secretion occurred in the young ewe lamb but was lost as the lamb matured: (3) attainment of sexual maturity was preceded by an elevation in LH pulse frequency.     

Cladistics and the land plants: A response to Robinson

The reconstruction of phylogenetic relationships should be based not on belief but on an explicit and logical analysis of all available characters. Hennigian phylogenetic systematics (cladistics) provides a framework for evaluating putative homologies characterizing particular hierarchical levels, determining relationships of taxa sharing congruent patterns of homologies, and constructing a classification based on this information. Fossils can and should be included in the analysis if enough of the relevant characters are preserved; this is not currently possible for early land plants because of the fragmentary fossil record. To avoid circularity, adaptive and functional considerations should be addressed only after a phylogenetic hypothesis based on patterns on patterns of shared homologies is available.     

Does decreasing extracellular Na inhibit in vitro renin secretion by affecting NaCa exchange?

It has been shown previously that in vitro renin secretion is inhibited by partial replacement of extracellular NaCl with either mannitol or choline chloride; the inhibitory effect is attributed to an increase in intracellular Ca, resulting from a decreased rate of Ca efflux via Na-Ca exchange. In the present experiments, we confirmed that partially replacing NaCl with choline chloride inhibited renin secretion from rat renal cortical slices, but we found that atropine completely blocked the effect, suggesting cholinergic mediation. Partially replacing NaCl with mannitol also inhibited renin secretion, but the effect could not be attributed specifically to a reduction in extracellular Na. Moreover, the stimulatory effect of Ca chelation on renin secretion was antagonized by either mannitol- or choline chloride -containing incubation media. These results do not support the hypothesis that lowering extracellular Na inhibits renin secretion by a mechanism involving decreased Ca efflux via Na-Ca exchange.     

THE RECONSTRUCTION OF ANCESTRAL CHARACTER STATES

The problem of error in the phylogenetic reconstruction of ancestral character states is explored by developing the model of Frumhoff and Reeve (1994). Information about the evolutionary rate of change within a character is inferred from the distribution of its character states on a known phylogeny, and this information is used to impose confidence limits on the error associated with ancestral state inference. Ancestral state inference is found to be remarkably robust under the model assumptions for a wide range of parameter values; however, the probability of error increases when the number of species within a clade is small and/or state-transition probabilities are strongly skewed in favor of the non-ancestral state. The rationale for expecting such a skew, a hypothesis of parallelism, is shown to rely on assumptions of low rates of change in at least two phylogenetically inherited characters: the tendency to occupy a particular ecological niche and the tendency to respond in a particular way to selection. A means for judging the relative likelihoods of parallelism vs. straightforward homology as explanations for a given character-state distribution is suggested. General problems with the model are discussed, as are methods for making it more realistic.