Evidence for Depressed Reproduction of Golden Eagles in Washington

Beginning in 1977 the Washington Department of Fish and Wildlife conducted annual surveys to determine statewide golden eagle (Aquila chrysaetos) occupancy and productivity. Current interest in the regional and national status of the species prompted our investigation to determine utility of historical data in assessing trends in reproduction, and to test efficacy of a sampling protocol that surveyed randomly selected territories and also accounted for detection probability. We found evidence indicating poor reproduction from 38 annual surveys conducted at 301 known territories statewide between 1977 and 2014. At 256 territories in eastern Washington, USA, apparent occupancy was low ( = 50.9%) and nesting success declined by 22%. All reproductive parameters were higher than at 45 territories in western Washington. We tested efficacy of a sampling protocol in 2013 and 2014 by surveying 108 randomly selected eastern territories. Probability of detecting eagles for these years from ground (= 89%) was greater than from air (= 66%). Our estimate of territory occupancy, corrected by probability of detection, was lower in 2013 (= 56.7%, 95% CI = 46.3–66.7%) than in 2014 (= 73.7%, 95% CI = 64.8–81.7%), as was the estimated number of breeding pairs (2013: = 158, 95% CI = 151–164; 2014: = 187, 95% CI = 182–192). Higher productivity (young/occupied territory) in 2013 (= 0.59, 95% CI = 0.40–0.82) than in 2014 (= 0.41, 95% CI = 0.27–0.59) and lower proportions of ≥1 immature eagle among nesting pairs in 2013 (16%) than in 2014 (31%), suggested higher immature pairing among sampled pairs contributed to inter-year differences in these reproductive parameters. Current and historical evidence for depressed golden eagle nesting in Washington is consistent with documented effects from habitat conversion, prey declines, lead contamination, and wind power development. We recommend future surveys in eastern Washington adhere to the random sampling protocol and conduct surveys at regular intervals to allow for trend analysis of reproductive parameters to better monitor golden eagle status. Surveys in western Washington, conducted exclusively from ground at all nests, will improve detection and cost efficiency. © 2020 The Wildlife Society.     

MATERNAL INHERITANCE AND ITS EFFECT ON ADAPTIVE EVOLUTION: A QUANTITATIVE GENETIC ANALYSIS OF MATERNAL EFFECTS IN A NATURAL PLANT POPULATION

A mother can influence a trait in her offspring both by the genes she transmits (Mendelian inheritance) and by maternal attributes that directly affect that trait in her offspring (maternal inheritance). Maternal inheritance can alter the direction, rate, and duration of adaptive evolution from standard Mendelian models and its impact on adaptive evolution is virtually unexplored in natural populations. In a hierarchical quantitative genetic analysis to determine the magnitude and structure of maternal inheritance in the winter annual plant, Collinsia verna, I consider three potential models of inheritance. These range from a standard Mendelian model estimating only direct (i.e., Mendelian) additive and environmental variance components to a maternal inheritance model estimating six additive and environmental variance components: direct additive and environmental variances; maternal additive and environmental variances; and the direct-maternal additive () and environmental covariances. The structure of maternal inheritance differs among the 10 traits considered at four stages in the life cycle. Early in the life cycle, seed weight and embryo weight display substantial , a negative , and a positive . Subsequently, cotyledon diameter displays and of roughly the same magnitude and negative . For fall rosettes, leaf number and length are best described by a Mendelian model. In the spring, leaf length displays maternal inheritance with significant and and a negative . All maternally inherited traits show significant negative . Predicted response to selection under maternal inheritance depends on and as well as . Negative results in predicted responses in the opposite direction to selection for seed weight and embryo weight and predicted responses near zero for all subsequent maternally inherited traits. Maternal inheritance persists through the life cycle of this annual plant for a number of size-related traits and will alter the direction and rate of evolutionary response in this population.     

RISK OF POPULATION EXTINCTION FROM FIXATION OF NEW DELETERIOUS MUTATIONS

The fixation of new deleterious mutations is analyzed for a randomly mating population of constant size with no environmental or demographic stochasticity. Mildly deleterious mutations are far more important in causing loss of fitness and eventual extinction than are lethal and semilethal mutations in populations with effective sizes, N_e, larger than a few individuals. If all mildly deleterious mutations have the same selection coefficient, s against heterozygotes and 2s against homozygotes, the mean time to extinction, , is asymptotically proportional to for 4N_es > 1. Nearly neutral mutations pose the greatest risk of extinction for stable populations, because the magnitude of selection coefficient that minimizes is about ? = 0.4/N_e. The influence of variance in selection coefficients among mutations is analyzed assuming a gamma distribution of s, with mean and variance . The mean time to extinction increases with variance in selection coefficients if is near ?, but can decrease greatly if is much larger than ?. For a given coefficient of variation of , the mean time to extinction is asymptotically proportional to for . When s is exponentially distributed, (c = 1) is asymptotically proportional to . These results in conjunction with data on the rate and magnitude of mildly deleterious mutations in Drosophila melanogaster indicate that even moderately large populations, with effective sizes on the order of N_e = 10³, may incur a substantial risk of extinction from the fixation of new mutations.     

Activation of Terminal Components of Human Complement by a Trypsin-Activated Complex of Human Factor B and Cobra Venom Factor

Cleavage of C3 by CVF-B? was demonstrated by hemolytic, immunoelectrophoretic and immune adherence reactions. No cleavage of C5 was detected by immunoelectrophoresis, but C5 hemolytic activity, assayed with decreased although less than C3 hemolytic activity. The co-existence of C3 with limiting amounts of C5 did not reduce the final degree of hemolysis of guinea pig erythrocytes (GPE) induced by late-acting components C6 through C9 and CVF-B?. Thus, a CVF-B? hemolytic system composed of GPE, C5 through C9 and CVF-B? provided a method for titration of terminal components of human complement. CVF-B? was able to generate hemolytically active sites of on GPE by activation of C5, C6 and C7. The complex in the fluid-phase decayed within 1 min but on GPE was quite stable. Originally insensitive sheep erythrocytes became sensitive to the CVF-B? hemolytic system if C3b sites were present, suggesting that cell-bound C3b played a role in orienting the positions of to be fixed. CVF-B? could be recovered quantitatively from the supernatant of the reaction mixture in which the hemolytically active intermediate GPEC- had been formed through the interaction between C5 to C8 and CVF-B?.     

Application of the\mathop m\limits^* - m method to the analysis of spatial patterns by changing the quadrat sizemethod to the analysis of spatial patterns by changing the quadrat size

A method for the analysis of spatial pattern using quadrats of different sizes is developed on the basis of the relationship of mean crowding () to mean density (m). The -on-m regression obtained by successive changes in quadrat size in a single population (unit-size relation) shows a characteristic pattern according to the type of distribution. By aid of the ρ-index proposed here, we can distinguish the random, aggregated and uniform distributions of the basic components (individual or group of individuals). The ρ serves as an index of spatial correlation between neighbouring quadrats, and it also provides information on the approximate area occupied by clump (colony), distribution pattern of individuals within clumps, and possibly the distribution pattern of clumps themselves. Even in a specified type of distribution, the unit-size relation is not necessarily identical with the relation for a series of populations at a particular quadrat size (series relation). The changes in the series relationship with successive changes of quadrat sizes are also considered for some basic patterns of distributions. The combined use of the unit-size and the series relations for a set of populations of the species under study may provide a satisfactory picture of the spatial pattern characteristic of the species. Application of the method is illustrated by using distribution data of several species of animals and plants. The advantage of the present method over other methods are discussed, and the formulae for determining the optimum quadrat unit in sampling surveys are given.     

THE MAINTENANCE OF POLYGENIC VARIATION IN FINITE POPULATIONS

Models of the maintenance of genetic variance in a polygenic trait have usually assumed that population size is infinite and that selection is weak. Consequently, they will overestimate the amount of variation maintained in finite populations. I derive approximations for the equilibrium genetic variance, in finite populations under weak stabilizing selection for triallelic loci and for an infinite “rare alleles” model. These are compared to results for neutral characters, to the “Gaussian allelic” model, and to Wright's approximation for a biallelic locus under arbitrary selection pressures. For a variety of parameter values, the three-allele, Gaussian, and Wrightian approximations all converge on the neutral model when population size is small. As expected, far less equilibrium genetic variance can be maintained if effective population size, N, is on the order of a few hundred than if N is infinite. All of the models predict that comparisons among populations with N less than about 10⁴ should show substantial differences in . While it is easier to maintain absolute when alleles interact to yield dominance or overdominance for fitness, less additivity also makes more susceptible to differences in N. I argue that experimental data do not seem to reflect the predicted degree of relationship between N and . This calls into question the ability of mutation-selection balance or simple balancing selection to explain observed . The dependence of on N could be used to test the adequacy of mutation-selection balance models.     

Mortality process in relation to aggregation in the southern green stink bug

A simple experiment of simulation was done to analyze the natural mortality process of young larval colonies and egg masses of the southern green stink bug. In this experiment, a degree of contagiousness was allowed in regard to the action of a mortality factor, and was defined as the mean number killed per a colony or an egg mass by the mortality factor within a unit time and the number killed per a colony was assumed to follow the Poisson series with the mean . Thus each component of the Poisson series was opposed to each colony or egg mass which was taken at random from 162 egg masses, 135 and 117 colonies of the first and the second instar larvae, respectively. It was revealed that mortality factors in the field did not act with a small degree of contagiosness, e. g., on all colonies or egg masses, but acted with a large degree of contagiousness, e. g., on some of the colonies or egg masses. Thus differential survival somewhat in all or none way occurred among the insect colonies irrespective of their initial sizes. These results were well explained by taking actual mortality factors into account.     

A new proposal for measurement of the adjusted mean crowding through consideration of size variability in habitat units

The mean crowding has previously been measured under the assumption that all quadrats or habitat units have the same size, even though the actual habitat units such as seeds or leaves are generally variable in size. A new index, ‘adjusted mean crowding’, which is adjusted for this variability can be given as where Q is the total number of habitat units in the whole area, x_j the number of individuals in the jth habitat unit, and a_j is defined as the ‘relative size’ of the jth habitat unit, i.e. a_y=y_y/(∑y_j/Q) where y_j is the actually measured size of the jth habitat unit. It is expected that and for the uniform distribution and the random distribution ‘per unit size’, respectively. The comparison between and regressions ( analysis) for the egg distribution pattern of Callosobruchus chinensis or C. maculatus proved that the regression is biased by a positive correlation between the egg number per seed and seed size rather than by a density-dependent change in the ovipositional behavior.     

Variance component tests of multivariate mediation effects under composite null hypotheses

Mediation effects of multiple mediators are determined by two associations: one between an exposure and mediators (‐) and the other between the mediators and an outcome conditional on the exposure (‐). The test for mediation effects is conducted under a composite null hypothesis, that is, either one of the ‐ and ‐ associations is zero or both are zeros. Without accounting for the composite null, the type 1 error rate within a study containing a large number of multimediator tests may be much less than the expected. We propose a novel test to address the issue. For each mediation test , , we examine the ‐ and ‐ associations using two separate variance component tests. Assuming a zero‐mean working distribution with a common variance for the element‐wise ‐ (and ‐) associations, score tests for the variance components are constructed. We transform the test statistics into two normally distributed statistics under the null. Using a recently developed result, we conduct hypothesis tests accounting for the composite null hypothesis by adjusting for the variances of the normally distributed statistics for the ‐ and ‐ associations. Advantages of the proposed test over other methods are illustrated in simulation studies and a data application where we analyze lung cancer data from The Cancer Genome Atlas to investigate the smoking effect on gene expression through DNA methylation in 15 114 genes.     

Anthropogenic Disturbance and Population Viability of Woodland Caribou in Ontario

One of the most challenging tasks in wildlife conservation and management is to clarify how spatial variation in land cover due to anthropogenic disturbance influences wildlife demography and long-term viability. To evaluate this, we compared rates of survival and population growth by woodland caribou (Rangifer tarandus caribou) from 2 study sites in northern Ontario, Canada that differed in the degree of anthropogenic disturbance because of commercial logging and road development, resulting in differences in predation risk due to gray wolves (Canis lupus). We used an individual-based model for population viability analysis (PVA) that incorporated adaptive patterns of caribou movement in relation to predation risk and food availability to predict stochastic variation in rates of caribou survival. Field estimates of annual survival rates for adult female caribou in the unlogged ( 0.90) and logged ( 0.76) study sites recorded during 2010–2014 did not differ significantly (P > 0.05) from values predicted by the individual-based PVA model (unlogged: = 0.87; logged: 0.79). Outcomes from the individual-based PVA model and a simpler stage-structured matrix model suggest that substantial differences in adult survival largely due to wolf predation are likely to lead to long-term decline of woodland caribou in the commercially logged landscape, whereas the unlogged landscape should be considerably more capable of sustaining caribou. Estimates of population growth rates (λ) for the 2010–2014 period differed little between the matrix model and the individual-based PVA model for the unlogged (matrix model = 1.01; individual-based model = 0.98) and logged landscape (matrix model = 0.88; individual-based model = 0.89). We applied the spatially explicit PVA model to assess the viability of woodland caribou across 14 woodland caribou ranges in Ontario. Outcomes of these simulations suggest that woodland caribou ranges that have experienced significant levels of commercial forestry activities in the past had annual growth rates <0.89, whereas caribou ranges that had not experienced commercial forestry operations had population growth rates >0.96. These differences were strongly related to regional variation in wolf densities. Our results suggest that increased wolf predation risk due to anthropogenic disturbance is of sufficient magnitude to cause appreciable risk of population decline in woodland caribou in Ontario. © 2020 The Authors. The Journal of Wildlife Management published by Wiley Periodicals, Inc. on behalf of The Wildlife Society.     

Different fates of deposited and in a temperate forest in northeast China: a 15N tracer study

Increasing atmospheric reactive nitrogen (N) deposition due to human activities could change N cycling in terrestrial ecosystems. However, the differences between the fates of deposited and are still not fully understood. Here, we investigated the fates of deposited and , respectively, via the application of ¹⁵NH₄NO₃ and NH₄¹⁵NO₃ in a temperate forest ecosystem. Results showed that at 410 days after tracer application, most was immobilized in litter layer (50 ± 2%), while a considerable amount of penetrated into 0–5 cm mineral soil (42 ± 2%), indicating that litter layer and 0–5 cm mineral soil were the major N sinks of and , respectively. Broad‐leaved trees assimilated more ¹⁵N under NH₄¹⁵NO₃ treatment compared to under ¹⁵NH₄NO₃ treatment, indicating their preference for –N. At 410 days after tracer application, 16 ± 4% added ¹⁵N was found in aboveground biomass under treatment, which was twice more than that under treatment (6 ± 1%). At the same time, approximately 80% added ¹⁵N was recovered in soil and plants under both treatments, which suggested that this forest had high potential for retention of deposited N. These results provided evidence that there were great differences between the fates of deposited and , which could help us better understand the mechanisms and capability of forest ecosystems as a sink of reactive nitrogen.     

Biophysical analysis of the dynamics of calmodulin interactions with neurogranin and Ca2+/calmodulin‐dependent kinase II

Calmodulin (CaM) functions depend on interactions with CaM‐binding proteins, regulated by . Induced structural changes influence the affinity, kinetics, and specificities of the interactions. The dynamics of CaM interactions with neurogranin (Ng) and the CaM‐binding region of /calmodulin‐dependent kinase II (CaMKII_290−309) have been studied using biophysical methods. These proteins have opposite dependencies for CaM binding. Surface plasmon resonance biosensor analysis confirmed that and CaM interact very rapidly, and with moderate affinity ( ). Calmodulin‐CaMKII_290−309 interactions were only detected in the presence of , exhibiting fast kinetics and nanomolar affinity ( ). The CaM–Ng interaction had higher affinity under ‐depleted ( and k ₋₁ = 1.6 × 10⁻¹s⁻¹) than ‐saturated conditions ( ). The IQ motif of Ng (Ng₂₇₋₅₀) had similar affinity for CaM as Ng under ‐saturated conditions ( ), but no interaction was seen under ‐depleted conditions. Microscale thermophoresis using fluorescently labeled CaM confirmed the surface plasmon resonance results qualitatively, but estimated lower affinities for the Ng ( ) and CaMKII_290−309( ) interactions. Although CaMKII_290−309 showed expected interaction characteristics, they may be different for full‐length CaMKII. The data for full‐length Ng, but not Ng₂₇₋₅₀, agree with the current model on Ng regulation of /CaM signaling.     

Optimal Strategies for Managing Wildlife Harvest Under System Change

Wildlife populations are experiencing shifting dynamics due to climate and landscape change. Management policies that fail to account for non-stationary dynamics may fail to achieve management objectives. We establish a framework for understanding optimal strategies for managing a theoretical harvested population under non-stationarity. Building from harvest theory, we develop scenarios representing changes in population growth rate () or carrying capacity () and derive time-dependent optimal harvest policies using stochastic dynamic programming. We then evaluate the cost of falsely assuming stationarity by comparing the outcomes of forward projections in which either the optimal policy or a stationary policy is applied. When declines over time, the stationary policy leads to an underharvest of the population, resulting in less harvest over the short term but leaving the population in a higher-value state. When declines over time, the stationary policy leads to overharvest, resulting in greater harvest returns in the short term but leaving the population in a lower and potentially more vulnerable state. This work demonstrates the basic properties of time-dependent harvest management and provides a framework for evaluating the many outstanding questions about optimal management strategies under climate change. Published 2021. This article is a U.S. Government work and is in the public domain in the USA.     

Analysis of spatial association between two species based on the interspecies mean crowding

The measurement of spatial association between two species is considered on the basis of interspecies mean crowding. Two indices of overlapping, γ and C_p, are derived as geometric and weighted arithmetic means of the same component ratios related to inter-and intraspecies mean crowdings. Both indices behave in a similar way, ranging from 1 when the distributions of two species are completely overlapped to 0 when they are completely exclusive with each other. The former is essentially identical with indices proposed byKuno (1968) andPianka (1973), and the latter is a modified form ofMorisita's (1959)C_σ index. Indices to measure the degree of spatial correlation between species, Ω and R_μ, are then derived for both kinds of overlapping indices, which vary from 1 in complete overlapping, through 0 in independent occurrence, to −1 in complete exclusion. Various kinds of interspecies association are analyzed using these indices and an extended form of the regression graph which provides a convenient way of indicating the spatial interrelation between two species as well as distribution patterns of respective species. The method presented in this paper may also be applicable to compare temporal distribution patterns between species, similarity between communities, etc. For such a wider application which includes continuous as well as discrete distributions, the interpretation of intra-and interspecies mean crowdings is not necessarily appropriate, and hence the concept of mean concentration with the symbols and for intraspecies relation and and for interspecies relation is suggested.     

Nitrate or ammonium: Influences of nitrogen source on the physiology of a green alga

In freshwaters, algal species are exposed to different inorganic nitrogen (N_i) sources whose incorporation varies in biochemical energy demand. We hypothesized that due to the lesser energy requirement of ammonium ()‐use, in contrast to nitrate ()‐use, more energy remains for other metabolic processes, especially under CO₂‐ and phosphorus (P_i) limiting conditions. Therefore, we tested differences in cell characteristics of the green alga Chlamydomonas acidophila grown on or under covariation of CO₂ and P_i‐supply in order to determine limitations, in a full‐factorial design. As expected, results revealed higher carbon fixation rates for ‐grown cells compared to growth with under low CO₂ conditions. ‐grown cells accumulated more of the nine analyzed amino acids, especially under P_i‐limited conditions, compared to cells provided with . This is probably due to a slower protein synthesis in cells provided with . In contrast to our expectations, compared to ‐grown cells ‐grown cells had higher photosynthetic efficiency under P_i‐limitation. In conclusion, growth on the N_i‐source did not result in a clearly enhanced C_i‐assimilation, as it was highly dependent on P_i and CO₂ conditions (replete or limited). Results are potentially connected to the fact that C. acidophila is able to use only CO₂ as its inorganic carbon (C_i) source.     

Dispersal of juvenile bobcats in a recolonizing population

Natural recolonization of bobcat (Lynx rufus) populations in the midwestern United States presents challenges for managers with limited knowledge of the species’ population dynamics in a highly fragmented agricultural landscape. Dispersal is a component of population dynamics of a recolonizing population, which is likely influenced by landscape features. To better understand population expansion, we examined dispersal characteristics of juvenile (<2 yr) bobcats in recently recolonized south-central Iowa, USA, from 2003–2009. We radio-collared and tracked dispersal fates of 61 individuals (34 males, 27 females), with 29 (22 males, 7 females) dispersing by approximately 2 years of age and the rest remaining philopatric. Most (65%) juvenile males dispersed, whereas only 26% of females dispersed. Initiation of dispersal varied, but none occurred in July–August. Dispersal duration ranged from 4–240 days. Average age at dispersal was 16.9 ± 1.1 (SE) months. Dispersal was most prevalent in an east-west direction and straight-line distances ranged from 6.6–203.2 km ( = 57.9 km). On average, males dispersed 44 km farther than females. Land cover composition was similar in natal and post-dispersal core ranges and consisted predominantly (70%) of forest and grassland. Lower abundance of forest and grassland in areas north of the study area may ultimately limit the ability of juvenile bobcats to successfully disperse and expand throughout the Corn Belt of Iowa. © 2019 The Wildlife Society.     

Tiller population dynamics of reciprocally transplanted Eriophorum vaginatum L. ecotypes in a changing climate

Moist tussock tundra, dominated by the sedge Eriophorum vaginatum L., covers approximately 3.36 × 10⁸ km² of arctic surface area along with large amounts of subarctic land area. Eriophorum vaginatum exhibits ecotypic differentiation along latitudinal gradients in Alaska. While ecotypic differentiation may be beneficial during periods of climate stability, it may be detrimental as climate changes, causing adaptational lag. Following harvest of a 30-year reciprocal transplant experiment, age-specific demographic data on E. vaginatum tillers were collected to parameterize a Leslie matrix. Yellow Taxi analysis, based on Tukey's Jackknife, was used to determine mean pseudovalues of tiller population growth rate () for four source populations of E. vaginatum tussocks that were transplanted to each of three gardens along a latitudinal gradient. Source populations responded differentially along the latitudinal gradient. Survival and daughter tiller production influenced differences seen at the mid-latitude garden, and the overall tiller population performance was generally improved by northward transplanting relative to southward transplanting. A comparison of home-source and away-source within the same transplant garden indicates no home-site advantage. Although populations were still growing when transplanted to home-sites ( = 1.056), tiller population growth rate increased as ΔGDD became more negative relative to the home site (i.e., as tussocks were transplanted north). These results imply that populations are affected by climate gradients in a manner consistent with adaptational lag. This study documenting the response of high-latitude ecotypes to climate gradients may be an indication of the possible future effects of climate shift in more southern latitudes.     

The Arabidopsis root stele transporter NPF2.3 contributes to nitrate translocation to shoots under salt stress

In most plants, constitutes the major source of nitrogen, and its assimilation into amino acids is mainly achieved in shoots. Furthermore, recent reports have revealed that reduction of translocation from roots to shoots is involved in plant acclimation to abiotic stress. NPF2.3, a member of the NAXT (nitrate excretion transporter) sub‐group of the NRT1/PTR family (NPF) from Arabidopsis, is expressed in root pericycle cells, where it is targeted to the plasma membrane. Transport assays using NPF2.3‐enriched Lactococcus lactis membranes showed that this protein is endowed with transport activity, displaying a strong selectivity for against Cl^?. In response to salt stress, translocation to shoots is reduced, at least partly because expression of the root stele transporter gene NPF7.3 is decreased. In contrast, NPF2.3 expression was maintained under these conditions. A loss‐of‐function mutation in NPF2.3 resulted in decreased root‐to‐shoot translocation and reduced shoot content in plants grown under salt stress. Also, the mutant displayed impaired shoot biomass production when plants were grown under mild salt stress. These mutant phenotypes were dependent on the presence of Na⁺ in the external medium. Our data indicate that NPF2.3 is a constitutively expressed transporter whose contribution to translocation to the shoots is quantitatively and physiologically significant under salinity.     

Minireview: demystifying microbial reaction energetics

The biology literature is rife with misleading information on how to quantify catabolic reaction energetics. The principal misconception is that the sign and value of the standard Gibbs energy () define the direction and energy yield of a reaction; they do not. is one part of the actual Gibbs energy of a reaction (ΔG_r ), with a second part accounting for deviations from the standard composition. It is also frequently assumed that applies only to 25 °C and 1 bar; it does not. is a function of temperature and pressure. Here, we review how to determine ΔG_r as a function of temperature, pressure and chemical composition for microbial catabolic reactions, including a discussion of the effects of ionic strength on ΔG_r and highlighting the large effects when multi-valent ions are part of the reaction. We also calculate ΔG_r for five example catabolisms at specific environmental conditions: aerobic respiration of glucose in freshwater, anaerobic respiration of acetate in marine sediment, hydrogenotrophic methanogenesis in a laboratory batch reactor, anaerobic ammonia oxidation in a wastewater reactor and aerobic pyrite oxidation in acid mine drainage. These examples serve as templates to determine the energy yields of other catabolic reactions at environmentally relevant conditions.