Impact of Mt. St. Helens ashfall on fruit yield of Mountain Huckleberry,Vaccinium membranaceum,important native American food

Huckleberry plants evaluated at 13 sites in the southern Cascade Mountains of Washington State during August 1980 showed significantly lower fruit yields where subjected to heavy ash deposition following the 18 May 1980 eruption of Mt. St. Helens. High insect pollinator mortality is suspected as the major causal factor. The impact of such effects of volcanic activity on Native American subsistence is discussed.     

Antagonism of the ATP-gated P2X7 receptor: a potential therapeutic strategy for cancer

Purinergic Signalling - The P2X receptor 7 (P2X7R) is a plasma membrane receptor sensing extracellular ATP associated with a wide variety of cellular functions. It is most commonly expressed on...     

Lanthanide titanates and LnxBi4−xTi3O12 solid solutions; their raman spectra

The formation of Ln_xBi_4?xTi₃O₁₂ solid solutions, where x_max ranges from 0.08 for Lu to 2.92 for La, has been examined by Raman spectroscopy. A review of the literature concerning Ln₂O₃Ln₂Ti₂O₇, solid solutions and Ln₂TiO₅ is presented and the Raman spectra of many of these phases are reported.     

The IRG protein-based resistance mechanism in mice and its relation to virulence in Toxoplasma gondii

IRG proteins (immunity-related GTPases) provide an early defense mechanism in mice against the protozoal pathogen, Toxoplasma gondii. This is a particularly suitable time to provide a brief review of this host-pathogen interaction because the nature of the IRG resistance system, and to some extent its mode of action, have become known in the past few years. Likewise, forward genetic screens have recently drawn attention to a number of loci contributing to the differential virulence of T. gondii strains in mice. It is now clear that at least some important virulence mechanisms exert their action against components of the IRG resistance system. Thus these two mechanisms form the two poles of a dynamic host-pathogen virulence-resistance relationship with interesting and accessible properties.     

Gene Genealogies Strongly Distorted by Weakly Interfering Mutations in Constant Environments

Neutral nucleotide diversity does not scale with population size as expected, and this “paradox of variation” is especially severe for animal mitochondria. Adaptive selective sweeps are often proposed as a major cause, but a plausible alternative is selection against large numbers of weakly deleterious mutations subject to Hill–Robertson interference. The mitochondrial genealogies of several species of whale lice (Amphipoda: Cyamus) are consistently too short relative to neutral-theory expectations, and they are also distorted in shape (branch-length proportions) and topology (relative sister-clade sizes). This pattern is not easily explained by adaptive sweeps or demographic history, but it can be reproduced in models of interference among forward and back mutations at large numbers of sites on a nonrecombining chromosome. A coalescent simulation algorithm was used to study this model over a wide range of parameter values. The genealogical distortions are all maximized when the selection coefficients are of critical intermediate sizes, such that Muller''s ratchet begins to turn. In this regime, linked neutral nucleotide diversity becomes nearly insensitive to N. Mutations of this size dominate the dynamics even if there are also large numbers of more strongly and more weakly selected sites in the genome. A genealogical perspective on Hill–Robertson interference leads directly to a generalized background-selection model in which the effective population size is progressively reduced going back in time from the present.OBSERVED levels of apparently neutral nucleotide diversity (π_n) are typically lower than expected under the assumptions of standard equilibrium theories, and they vary much less among species than do estimates of long-term effective population sizes (Nei and Grauer 1984; Bazin et al. 2006; Nabholz et al. 2008). Many explanations have been proposed for the apparent shortfalls and the lack of proportionality with population size, including (1) complex demographic histories (e.g., recurring population bottlenecks), (2) adaptive selective sweeps (Maynard Smith and Haigh 1974; Gillespie 1999), and (3) selection against deleterious mutations (Charlesworth et al. 1993, 1995; McVean and Charlesworth 2000; Comeron et al. 2008). Of these three possibilities, bottlenecks and sweeps are by far the most frequently mentioned, even though deleterious mutations occur at high rates in all species, regardless of ecological circumstances (Eyre-Walker and Keightley 2007). Here we show that weakly deleterious mutations can distort genealogies in three different ways and dramatically reduce nucleotide diversities in large populations of nonrecombining chromosomes. The mitochondrial genealogies of several species of whale lice (Kaliszewska et al. 2005) are distorted in exactly these ways, and several lines of evidence suggest that bottlenecks and adaptive sweeps are not likely to be the primary causes.Mitochondria have been proposed to be especially sensitive to selective sweeps. Animal mitochondrial genomes contain more than three dozen essential protein and structural RNA genes, so they are large targets for both mutation and selection (Ballard and Whitlock 2004). They do not undergo sexual recombination, so every advantageous mutation that fixes will reduce variation throughout the genome. Mitochondrial nucleotide diversity therefore could depend strongly on rates of environmental change, which could be similar for species with very different population sizes. Indeed, if rates of mitochondrial adaptation were mutation limited, then larger populations might actually experience higher rates of adaptive substitution and as a result show lower average levels of neutral diversity than smaller populations (Gillespie 2000, 2001). This idea was recently invoked to explain the remarkable similarity of average levels of mitochondrial nucleotide diversity among the major animal classes which appear to have very different average population sizes and substantially different average levels of nuclear nucleotide and amino acid diversity (Bazin et al. 2006).Unconditionally deleterious mutations can also depress linked neutral diversity by reducing the effective population size either through (1) background selection against relatively strongly selected mutations (Charlesworth et al. 1993, 1995) or (2) Hill–Robertson interference (Hill and Robertson 1966) among large numbers of relatively weakly selected mutations (reviewed by Comeron et al. 2008). The second of these processes, called “weak-selection Hill–Robertson interference” (wsHRi) by McVean and Charlesworth (2000) and “interference selection” (IS) by Comeron and Kreitman (2002), can shorten genealogies, give them strongly nonneutral branch-length proportions, and skew their topologies (Higgs and Woodcock 1995; Maia et al. 2004).To date, weak interference has mainly been studied by forward simulation, with the aim of assessing its possible effects on patterns of optimal synonymous codon use within eukaryotic nuclear genes and genomes, in the presence of recombination (Comeron and Guthrie 2005; Loewe and Charlesworth 2007; Comeron et al. 2008). In an attempt to understand the striking genealogical distortions seen in whale-louse mitochondria (Kaliszewska et al. 2005), we have developed a structured-coalescent algorithm that accurately models selection of arbitrary strength on a nonrecombining chromosome of finite length. All of the distortions seen in the whale-louse mitochondria are replicated under parameters that might plausibly apply to whale lice and many other animal species, and these distortions scale only weakly with population size.Whale lice are permanent, obligate ectoparasites of cetaceans. They feed on the dead outer surface of their host''s skin, and they appear to be harmless. They are amphipod Crustacea comprising a monophyletic family, Cyamidae, with ∼50 described species in several genera. Three of these species (Cyamus ovalis, C. gracilis, and C. erraticus) occur on right whales (Eubalaena spp.) but not regularly on any other hosts. Most adult right whales carry large populations of all three species.Right whales in the North Pacific, the North Atlantic, and the southern hemisphere have been separated for ∼5 million years, and so have their cyamids (Rosenbaum et al. 2000; Gaines et al. 2005; Kaliszewska et al. 2005). For this reason the right whales in different ocean systems are now considered distinct species (Eubalaena japonica, E. glacialis, and E. australis), and we refer to their cyamids as North Pacific C. ovalis, North Atlantic C. ovalis, southern C. ovalis, and so on, in anticipation that a future revision of the genus Cyamus will recognize them as “triplet” sibling species (3 × 3 = 9 species in all). We studied their mitochondrial population genetics with the initial aim of quantifying patterns of genetic differentiation among the cyamid populations on individual whales within local populations (Kaliszewska et al. 2005).We had reasoned (incorrectly) that the pattern of differentiation among whales might say something about their social interactions, since cyamids can transfer only between whales that are in direct physical contact with each other. We found very low levels of differentiation among whales and to our surprise literally no differentiation among the major southern hemisphere breeding aggregations that calve off the coasts of South America, South Africa, Australia, and New Zealand (Kaliszewska et al. 2005). This absence of population structure seems remarkable by terrestrial standards but is easily explained by modest rates of cyamid exchange among whales within local populations and between the major breeding aggregations, given the enormous sizes of cyamid populations. Right whales are highly gregarious (spending hours per day in social interactions), mobile (traveling thousands of kilometers per year on annual foraging migrations), and mortal (carrying their cyamid populations to the sea floor when they die). Thus cyamids have many opportunities to transfer between whales, and they might be expected to have evolved an inclination to do so when the opportunity presents itself (Hamilton and May 1977).The well-defined ecology of right-whale cyamids allows their population sizes to be estimated directly. The number of adult cyamids per whale (∼500–10,000, varying by species) times the number of whales per ocean (∼50,000–200,000, prior to human exploitation) equals the number of cyamids per species (Kaliszewska et al. 2005). Thus for all three nominal species of right-whale cyamids, long-term census population sizes are expected to have been in the range 2.5 × 10⁷–2 × 10⁹. Given conservative estimates of the per-generation mitochondrial mutation rate, even the lower end of this range predicts levels of synonymous nucleotide diversity at least an order of magnitude larger than those actually seen in the cyamids, which are consistently modest and similar to those seen in typical terrestrial arthropods (Kaliszewska et al. 2005). The three North Atlantic and southern hemisphere sibling-species pairs are strongly reciprocally monophyletic, as illustrated for C. ovalis in Figure 1. This is not expected at mutation–drift equilibrium, given their very large population sizes.Open in a separate windowFigure 1.—Mitochondrial gene genealogies for North Atlantic and southern hemisphere Cyamus ovalis, estimated by UPGMA from partial COI sequences. Left: The intraspecific genealogies coalesce globally at ∼0.5 and 1 MY and share an ancestor at ∼5 MY (Kaliszewska et al. 2005). Center and right: Each intraspecific genealogy is aligned with a generalized skyline plot (Strimmer and Pybus 2001) showing estimates of θ at different times in the past under a piecewise constant model of population size change fit by GENIE 3.0 (Pybus and Rambaut 2002). Three different point estimates of present-day θ are also indicated on the plots (synonymous-site nucleotide diversity π, Watterson''s θ estimated from synonymous sites, and θ estimated jointly with the apparent exponential growth rate by the MCMC coalescent algorithm in LAMARC). Values of Tajima''s D are lower (−1.5 to −1.6) when estimated from the sequences than when estimated from the branch lengths of the trees (−2.1 to −2.4), as expected because multiple substitutions occur at some sites. Similar values of D_T (−2.3 for both species) were obtained from the highest-likelihood trees found by BEAST with a fully parameterized GTR substitution model and a coalescent prior. Those trees have standardized imbalance statistics (−I_S) of −3.7 (southern hemisphere C. ovalis) and −2.2 (North Atlantic C. ovalis). The trees shown here have more extreme values of −I_S (−5.5 and −4.0, respectively), probably as a consequence of artificially pectinate branching orders induced by UPGMA among sets of identical sequences. Slightly different sets of sequences were used to make the two-species genealogy on the left and the single-species genealogies in the center. The long interspecific branches in the two-species tree are based on a multispecies maximum-likelihood analysis involving smaller numbers of much longer (4.1 kb) sequences (Kaliszewskaet al. 2005, Figure 3).These dramatic deficits of variation are not easily explained by population bottlenecks or adaptive selective sweeps. The bottleneck hypothesis is especially problematic because it requires the long-term near extinction of right whales in all three ocean systems. (Short bottlenecks such as those caused by human exploitation of right whales are not expected to have a noticeable effect on cyamid genetic diversity because cyamid populations remain large, and rates of genetic drift low, even when there are few whales.) That all three right whales have survived for millions of years suggests that they have maintained reasonably large population sizes, and the mitochondrial nucleotide diversity of southern right whales is consistent with this assumption (Kaliszewska et al. 2005), as is the nucleotide diversity of a cyamid nuclear gene (described below). Right whales eat copepods and krill, which are relatively close to the base of marine food webs, and right-whale populations are thought to be food limited. Thus a long-term, severe depression of their numbers would also seem to imply a collapse of marine ecosystems worldwide, for which there is no evidence.Several features of cyamid mitochondrial nucleotide diversity are also inconsistent with the bottleneck model and with adaptive sweeps as well. The most obvious of these features is the uniformity of cyamid mitochondrial diversity among species (π = 0.007–0.015 for COI sequences in the seven species surveyed by Kaliszewska et al. 2005). Gene genealogies estimated from these sequences also seem remarkably uniform in total depth, with last common ancestors differing in age by only a factor of 3 and in six of the seven species by less than a factor of 2 (see Kaliszewska et al. 2005, Figure 4). Adaptive sweeps might be expected to occur at roughly random intervals and not to be well coordinated in time among seven species in three different ocean systems. Interspecific coordination of such sweeps (over the whole globe) would seem to be required if they were to be a plausible primary cause of the genealogical shortening.Open in a separate windowFigure 4.—Apparent average effective sizes of ancestral populations, for models with different values of s. Values of N_e are given on the vertical axis (logarithmically scaled). They are estimated from the variance of expected contributions to the present, for the adults of any given generation. The parameters are those of Figure 2 (N = 65,536 = “64k,” μ = 1.5 × 10⁻⁶, L_s = 2048, U = 0.0031). In theory the curve for s = 0 should be perfectly horizontal. The discrepancy appears to be caused by subtle flaws in the shape of the very broad “idealized” simulated distribution that was used in this calculation. The curves for strong-selection cases are perfectly horizontal at times beyond a few thousand generations because the mutation-number distributions are compact, with little stochastic variation. The two lowest curves are those for the selection coefficients (s = 2⁻¹⁰, U/s = 3.2, and s = 2⁻¹¹, U/s = 6.3) that produce the most extreme values of the polymorphism and tree-shape statistics, given the other parameters (Figure 2).In addition to showing too little nucleotide variation, the cyamid mitochondrial genomes show strong and consistent excesses of rare nucleotide states, reflecting the “comb-like” or “star-like” shapes of the genealogies, in which deeper branches tend to be much too short relative to terminal branches (as if the trees had been “squished” from behind). This kind of distortion causes negative values of Tajima''s (1989) D and related statistics. It can be caused by population expansion from a bottleneck or by lineage expansion under positive selection (Kaplan et al. 1989; Slatkin and Hudson 1991; Rogers and Harpending 1992; Bamshad and Wooding 2003). However, the form of branch-length distortion seen in the cyamid genealogies suggests a slow, steady, roughly exponential form of population or lineage growth, not the relatively sudden increases suggested by the bottleneck and selective-sweep hypotheses. Generalized skyline plots (Strimmer and Pybus 2001) describing the histories of population size implied by the shapes of the northern and southern C. ovalis gene genealogies are shown in Figure 1. They are remarkably similar, as are the growth rates and estimates of present-day θ (= 2N_fμ_n) obtained by fitting exponential growth models using the coalescent algorithms in LAMARC (Kuhner et al. 1998, 2004) or BEAST (Drummond and Rambaut 2007).Cyamid populations cannot have grown in numbers as seemingly implied by these analyses. The number of cyamids on each whale appears to be set mainly by microhabitat limitations (e.g., by the area of rough callosity tissue on the head, where C. ovalis and C. gracilis live), and these features of their environment have hardly changed for millions of years, as demonstrated by the strong similarities of northern and southern right whales and their cyamids. Likewise, the numbers of right whales cannot have increased gradually from vanishingly small numbers over several hundred thousand years, for the reasons discussed above.The genealogical signals of “growth” therefore seem likely to be caused by selection. Environmental change is the most obvious potential cause of selection, but the apparent rate of growth seen here is strangely slow—in fact, slower than glacial. The orbitally forced Plio-Pleistocene glacial climate cycles have a major period of ∼100,000 years (Lambert et al. 2008 and references therein), but all seven of the right-whale cyamids for which we have mitochondrial population samples appear to have been “expanding,” more or less continuously, through at least several such cycles. The seemingly fairly consistent rate of branch-length foreshortening seen in the genealogies therefore suggests the action of a process that is relatively homogeneous in time, in addition to being very slow overall.The cyamid mitochondrial genealogies also appear to be topologically skewed, with sister clades too unequal in size, on average. In random bifurcating trees, the distribution of sister-clade sizes is uniform (Yule 1924; Heard 1992; Rogers 1994). Deviations from this null expectation can be quantified by statistics such as Colless''s (1982) index of tree imbalance (Shao and Sokal 1990; Rogers 1996). Our estimates of the cyamid genealogies tend to be excessively imbalanced (Figure 1). Strong topological imbalance is not caused by classic adaptive sweeps or by population growth following a bottleneck, but previous theoretical work has indicated that it can be caused by selection (Higgs and Woodcock 1995; Maia et al. 2004).To summarize, the cyamid mitochondrial genealogies are consistently much too short, too squished, and too skewed, relative to neutral-theory expectations. Owing to several special features of cyamid and right-whale biology, selection seems to be the only plausible explanation for this set of distortions, but conventional adaptive sweeps do not seem likely to be the primary cause. We therefore asked whether weakly deleterious mutations might be sufficient to generate the observed combination of patterns, in the absence of environmental change. Previous work (mentioned above) showed that interference among weak mutations at many sites can strongly affect linked neutral variation, but this work did not fully explore the parameter space relevant to our system or connect all the patterns in a genealogical setting.To address this question we first carried out forward simulations of populations of nonrecombining chromosomes with large numbers of nucleotide positions subject to forward and back mutations with unconditional fitness effects of size s. Large numbers of linked neutral sites were used to estimate genealogies and to calculate population statistics of interest. We found that for a range of intermediate values of s, considerable fitness variation was maintained and all three of the genealogical distortions (and the signal of apparent exponential growth) seen in the cyamid genealogies reached impressively large maxima. However, the computational burden of full forward simulation prevented us from considering realistic parameter values (i.e., large N and small μ), and it was not obvious that extrapolations based on compound parameters (e.g., Nμ and Ns) would work as hoped in all respects (see Comeron et al. 2008). We developed an equivalent coalescent algorithm that accurately reproduces all results of the forward simulations and allows for realistic parameter values. Under parameters relevant to cyamid mitochondria, the distortions of genealogical depth, proportions, and topology can be even more extreme than those seen in the cyamids, and the mean pairwise coalescence times (and resulting neutral nucleotide diversities) associated with maximally distorting intermediate values of s (U_s/s ∼ 5, where U_s is the total genomic mutation rate at sites with selection coefficents of size s) depend only weakly on N.All parameters of this model (including those of the environment) remain constant over time, yet in some respects it displays apparently nonequilibirum behavior. Under weak to intermediate selection (U_s/s > 10), the effective population size appears to become progressively smaller as time recedes into the past, giving rise to the illusion of growth. And in the maximally distorting range of intermediate selection coefficients, the distributions of deleterious mutation numbers and the shapes of genealogies show conspicuous dynamical instability of a form that could be taken to suggest “adaptive evolution” in response to episodes of environmental change. Adaptive mutations contribute importantly to this process, but they are reversions at some of the many sites previously mutated to mildly deleterious states. Subtle patterns of environmental change that converted previously optimal nucleotide states to slightly suboptimal states could give rise to a category of “virtual reversions” that would augment (or even outnumber) simple reversions, and the effects of such a process might well be consistent with the distortions seen in the cyamid mitochondrial genealogies. However, models with no environmental change of any kind appear to explain the observations surprisingly well.     

Immunity-related GTPase M (IRGM) proteins influence the localization of guanylate-binding protein 2 (GBP2) by modulating macroautophagy

The immunity-related GTPases (IRGs) are a family of proteins induced by interferon-γ that play a crucial role in innate resistance to intracellular pathogens. The M subfamily of IRG proteins (IRGM) plays a profound role in this context, in part because of the ability of its members to regulate the localization and expression of other IRG proteins. We present here evidence that IRGM proteins affect the localization of the guanylate-binding proteins (GBPs), a second family of interferon-induced GTP-binding proteins that also function in innate immunity. Absence of Irgm1 or Irgm3 led to accumulation of Gbp2 in intracellular compartments that were positive for both the macroautophagy (hereafter referred to as autophagy) marker LC3 and the autophagic adapter molecule p62/Sqstm1. Gbp2 was similarly relocalized in cells in which autophagy was impaired because of the absence of Atg5. Both in Atg5- and IRGM-deficient cells, the IRG protein Irga6 relocalized to the same compartments as Gbp2, raising the possibility of a common regulatory mechanism. However, other data indicated that Irga6, but not Gbp2, was ubiquitinated in IRGM-deficient cells. Similarly, coimmunoprecipitation studies indicated that although Irgm3 did interact directly with Irgb6, it did not interact with Gbp2. Collectively, these data suggest that IRGM proteins indirectly modulate the localization of GBPs through a distinct mechanism from that through which they regulate IRG protein localization. Further, these results suggest that a core function of IRGM proteins is to regulate autophagic flux, which influences the localization of GBPs and possibly other factors that instruct cell-autonomous immune resistance.