On the regulation of populations of mammals, birds, fish, and insects

A key unresolved question in population ecology concerns the relationship between a population's size and its growth rate. We estimated this relationship for 1780 time series of mammals, birds, fish, and insects. We found that rates of population growth are high at low population densities but, contrary to previous predictions, decline rapidly with increasing population size and then flatten out, for all four taxa. This produces a strongly concave relationship between a population's growth rate and its size. These findings have fundamental implications for our understanding of animals' lives, suggesting in particular that many animals in these taxa will be found living at densities above the carrying capacity of their environments.

Identifying adaptive genetic divergence among populations from genome scans

The identification of signatures of natural selection in genomic surveys has become an area of intense research, stimulated by the increasing ease with which genetic markers can be typed. Loci identified as subject to selection may be functionally important, and hence (weak) candidates for involvement in disease causation. They can also be useful in determining the adaptive differentiation of populations, and exploring hypotheses about speciation. Adaptive differentiation has traditionally been identified from differences in allele frequencies among different populations, summarised by an estimate of F-ST. Low outliers relative to an appropriate neutral population-genetics model indicate loci subject to balancing selection, whereas high outliers suggest adaptive (directional) selection. However, the problem of identifying statistically significant departures from neutrality is complicated by confounding effects on the distribution of F-ST estimates, and current methods have not yet been tested in large-scale simulation experiments. Here, we simulate data from a structured population at many unlinked, diallelic loci that are predominantly neutral but with some loci subject to adaptive or balancing selection. We develop a hierarchical-Bayesian method, implemented via Markov chain Monte Carlo (MCMC), and assess its performance in distinguishing the loci simulated under selection from the neutral loci. We also compare this performance with that of a frequentist method, based on moment-based estimates of F-ST. We find that both methods can identify loci subject to adaptive selection when the selection coefficient is at least five times the migration rate. Neither method could reliably distinguish loci under balancing selection in our simulations, even when the selection coefficient is twenty times the migration rate.

Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

Conserving marginal populations of the foodplant (Impatiens noli-tangere) of an endangered moth (Eustroma reticulatum) in a changing climate

Impatiens noli-tangere is scarce in the UK and probably only native to the Lake District and Wales. It is the sole food plant for the endangered moth Eustroma reticulattum. Significant annual fluctuations in the size of I. noli-tangere populations endanger the continued presence of E. reticulatum in the UK. In this study, variation in population size was monitored across native populations of L noli-tangere in the English Lake District and Wales. In 1998, there was a crash in the population size of all metapopulations in the Lake District but not of those found in Wales. A molecular survey of the genetic affinities of samples in 1999 from both regions and a reference population from Switzerland was performed using AFLP and ISSR analyses. The consensus UPGMA dendrogram and a PCO scatter plot revealed clear differentiation between the populations of L noli-tangere in Wales and those in the Lake District. Most of the genetic variation in the UK (H-T= 0.064) was partitioned between (G(ST) = 0.455) rather than within (H-S = 0.034) regions, inferring little gene flow occurs between regions. There was similar bias towards differentiation between metapopulations in Wales, again consistent with low levels of interpopulation gene flow. This contrasts with far lower levels of differentiation in the Lake District which suggests modest rates of gene flow may occur between populations. It is concluded that in the event of local extinction of sites or populations, reintroductions should be restricted to samples collected from the same region. We then surveyed climatic variables to identify those most likely to cause local extinctions. Climatic correlates of population size were sought from two Lake District metapopulations situated close to a meteorological station. A combination of three climatic variables common to both sites explained 81-84% of the variation in plant number between 1990 and 2001. Projected trends for these climatic variables were used in a Monte Carlo simulation which suggested an increased risk of I. noli-tangere population crashes by 2050 at Coniston Water. but not at Derwentwater. Implications of these findings for practical conservation strategies are explored. (C) 2003 Elsevier Ltd. All rights reserved.

Combining competition from Lolium perenne and an insect-fungus combination to control Rumex obtusifolius seedlings

Although adult Rumex obtusifolius are problematic weeds, their seedlings are poor competitors against Lolium perenne, particularly in established swards. We investigated the possibility of using this weakness to augment control of R. obtusifolius seedlings with combinations of Gastrophysa viridula (Coleoptera: Chrysomelidae) and the rust fungus Uromyces rumicis. Rumex obtusifolius seedlings were grown in competition with L. perenne sown at different rates and times after R. obtusifolius: they competed successfully with L. perenne when sown 21 days before the grass. Sowing both species at the same time resulted in a dominant grass sward, with R. obtusifolius becoming dominant when sown 42 days prior to L. perenne. Grass sowing rate had no effect on R. obtusifolius growth or biomass. A second experiment investigated how competition from L. perenne sown 21 days after R. obtusifolius combined with damage from G. viridula and/or U. rumicis (applied at either the 3-4- or 10-13-leaf stage, or at both stages) affected the growth and final biomass of R. obtusifolius. Beetle grazing at the latter leaf stage was the only treatment that reduced R. obtusifolius biomass, although rust infection at the earlier application led to an increase in shoot and root weight. The results are discussed in terms of the potential for use of these agents in the field.

Climate change and coastal waterbird populations - past declines and future impacts

Considerable attention has been given to the impact of climate change on avian populations over the last decade. In this paper we examine two issues with respect to coastal bird populations in the UK: (1) is there any evidence that current populations are declining due to climate change, and (2) how might we predict the response of populations in the future? We review the cause of population decline in two species associated with saltmarsh habitats. The abundance of Common Redshank Tringa totanus breeding on saltmarsh declined by about 23% between the mid-1980s and mid-1990s, but the decline appears to have been caused by an increase in grazing pressure. The number of Twite Carduelis flavirostris wintering on the coast of East Anglia has declined dramatically over recent decades; there is evidence linking this decline with habitat loss but a causal role for climate change is unclear. These examples illustrate that climate change could be having population-level impacts now, but also show that it is dangerous to become too narrowly focused on single issues affecting coastal birds. Making predictions about how populations might respond to future climate change depends on an adequate understanding of important ecological processes at an appropriate spatial scale. We illustrate this with recent work conducted on the Icelandic population of Black-tailed Godwits Limosa limosa islandica that shows large-scale regulatory processes. Most predictive models to date have focused on local populations (single estuary or a group of neighbouring estuaries). We discuss the role such models might play in risk assessment, and the need for them to be linked to larger-scale ecological processes. We argue that future work needs to focus on spatial scale issues and on linking physical models of coastal environments with important ecological processes.

Rapid parallel expression in E-coli and insect cells: Analysis of five lef gene products of the Autographa californica multiple nuclear polyhedrosis virus (AcMNPV)

A number of strategies are emerging for the high throughput (HTP) expression of recombinant proteins to enable structural and functional study. Here we describe a workable HTP strategy based on parallel protein expression in E. coli and insect cells. Using this system we provide comparative expression data for five proteins derived from the Autographa californica polyhedrosis virus genome that vary in amino acid composition and in molecular weight. Although the proteins are part of a set of factors known to be required for viral late gene expression, the precise function of three of the five, late expression factors (lefs) 6, 7 and 10, is unknown. Rapid expression and characterisation has allowed the determination of their ability to bind DNA and shown a cellular location consistent with their properties. Our data point to the utility of a parallel expression strategy to rapidly obtain workable protein expression levels from many open reading frames (ORFs).

A block in virus-like particle maturation following assembly of murine leukaemia virus in insect cells

Expression of the murine leukaemia virus (MLV) major Gag antigen p65(Gag) using the baculovirus expression system leads to efficient assembly and release of virus-like particles (VLP) representative of immature MLV. Expression of P180(Gag-Pol), facilitated normally in mammalian cells by readthrough of the p65(Gag) termination codon, also occurs efficiently in insect cells to provide a source of the MLV protease and a pattern of p65(Gag) processing similar to that observed in mammalian cells. VLP release from P180(Gag-Pol) expressing cells however remains essentially immature with disproportionate levels of the uncleaved p65(Gag) precursor when compared to the intracellular Gag profile. Changing the p65(Gag) termination codon altered the level of p65(Gag) and p180(Gag-Pol) within expressing cells but did not alter the pattern of released VLP, which remained immature. Coexpression of p65(Gag) with a fixed readthrough p180(Gag-Pol) also led to only immature VLP release despite high intracellular protease levels. Our data suggest a mechanism that preferentially selects uncleaved p65(Gag) for the assembly of MLV in this heterologous expression system and implies that, in addition to their relative levels, active sorting of the correct p65(Gag) and p180(Gag-Pol) ratios may occur in producer cells. (C) 2003 Elsevier Inc. All rights reserved.

Control of human immunodeficiency virus type-1 protease activity in insect cells expressing Gag-Pol rescues assembly of immature but not mature virus-like particles

Expression of human immunodeficiency virus type 1 (HIV-1) Gag protein in insect cells using baculovirus vectors leads to the abundant production of virus-like particles (VLPs) that represent the immature form of the virus. When Gag-Pol is included, however, VLP production is abolished, a result attributed to premature protease activation degrading the intracellular pool of Gag precursor before particle assembly can occur. As large-scale synthesis of mature noninfectious VLPs would be useful, we have sought to control HIV protease activity in insect cells to give a balance of Gag and Gag-Pol that is compatible with mature particle formation. We show here that intermediate levels of protease activity in insect cells can be attained through site-directed mutagenesis of the protease and through antiprotease drug treatment. However, despite Gag cleavage patterns that mimicked those seen in mammalian cells, VLP synthesis exhibited an essentially all-or-none response in which VLP synthesis occurred but was immature or failed completely. Our data are consistent with a requirement for specific cellular factors in addition to the correct ratio of Gag and Gag-Pol for assembly of mature retrovirus particles in heterologous cell types. (C) 2003 Elsevier Science (USA). All rights reserved.

Estimation of population growth or decline in genetically monitored populations

This article introduces a new general method for genealogical inference that samples independent genealogical histories using importance sampling (IS) and then samples other parameters with Markov chain Monte Carlo (MCMC). It is then possible to more easily utilize the advantages of importance sampling in a fully Bayesian framework. The method is applied to the problem of estimating recent changes in effective population size from temporally spaced gene frequency data. The method gives the posterior distribution of effective population size at the time of the oldest sample and at the time of the most recent sample, assuming a model of exponential growth or decline during the interval. The effect of changes in number of alleles, number of loci, and sample size on the accuracy of the method is described using test simulations, and it is concluded that these have an approximately equivalent effect. The method is used on three example data sets and problems in interpreting the posterior densities are highlighted and discussed.

Analysis of the sulfate-reducing bacterial and methanogenic archaeal populations in contrasting Antarctic sediments

The distribution and activity of communities of sulfate-reducing bacteria (SRB) and methanogenic archaea in two contrasting Antarctic sediments were investigated. Methanogenesis dominated in freshwater Lake Heywood, while sulfate reduction dominated in marine Shallow Bay. Slurry experiments indicated that 90% of the methanogenesis in Lake Heywood was acetoclastic. This finding was supported by the limited diversity of clones detected in a Lake Heywood archaeal clone library, in which most clones were closely related to the obligate acetate-utilizing Methanosaeta concilii. The Shallow Bay archaeal clone library contained clones related to the C-1-utilizing Methanolobus and Methanococcoides and the H-2-utilizing Methanogenium. Oligonucleotide probing of RNA extracted directly from sediment indicated that archaea represented 34% of the total prokaryotic signal in Lake Heywood and that Methanosaeta was a major component (13.2%) of this signal. Archaea represented only 0.2% of the total prokaryotic signal in RNA extracted from Shallow Bay sediments. In the Shallow Bay bacterial clone library, 10.3% of the clones were SRB-like, related to Desulfotalea/Desulforhopalus, Desulfofaba, Desulfosarcina, and Desulfobacter as well as to the sulfur and metal oxidizers comprising the Desulfuromonas cluster. Oligonucleotide probes for specific SRB clusters indicated that SRB represented 14.7% of the total prokaryotic signal, with Desulfotalea/Desulforhopalus being the dominant SRB group (10.7% of the total prokaryotic signal) in the Shallow Bay sediments; these results support previous results obtained for Arctic sediments. Methanosaeta and Desulfotalea/Desulforhopalus appear to be important in Lake Heywood and Shallow Bay, respectively, and may be globally important in permanently low-temperature sediments.

Development of a fermentation system to model sessile bacterial populations in the human colon

A fermentation system was designed to model the human colonic microflora in vitro. The system provided a framework of mucin beads to encourage the adhesion of bacteria, which was encased within a dialysis membrane. The void between the beads was inoculated with faeces from human donors. Water and metabolites were removed from the fermentation by osmosis using a solution of polyethylene glycol (PEG). The system was concomitantly inoculated alongside a conventional single-stage chemostat. Three fermentations were carried out using inocula from three healthy human donors. Bacterial populations from the chemostat and biofilm system were enumerated using fluorescence in situ hybridization. The culture fluid was also analysed for its short-chain fatty acid (SCFA) content. A higher cell density was achieved in the biofilm fermentation system (taking into account the contribution made by the bead-associated bacteria) as compared with the chemostat, owing to the removal of water and metabolites. Evaluation of the bacterial populations revealed that the biofilm system was able to support two distinct groups of bacteria: bacteria growing in association with the mucin beads and planktonic bacteria in the culture fluid. Furthermore, distinct differences were observed between populations in the biofilm fermenter system and the chemostat, with the former supporting higher populations of clostridia and Escherichia coli. SCFA levels were lower in the biofilm system than in the chemostat, as in the former they were removed via the osmotic effect of the PEG. These experiments demonstrated the potential usefulness of the biofilm system for investigating the complexity of the human colonic microflora and the contribution made by sessile bacterial populations.

Influence of environmental stress on distributions of times to first division in Escherichia coli populations, as determined by digital-image analysis of individual cells

The distributions of times to first cell division were determined for populations of Escherichia coli stationary-phase cells inoculated onto agar media. This was accomplished by using automated analysis of digital images of individual cells growing on agar and calculation of the "box area ratio." Using approximately 300 cells per experiment, the mean time to first division and standard deviation for cells grown in liquid medium at 37 degrees C and inoculated on agar and incubated at 20 degrees C were determined as 3.0 h and 0.7 h, respectively. Distributions were observed to tail toward the higher values, but no definitive model distribution was identified. Both preinoculation stress by heating cultures at 50 degrees C and postinoculation stress by growth in the presence of higher concentrations of NaCl increased mean times to first division. Both stresses also resulted in an increase in the spread of the distributions that was proportional to the mean division time, the coefficient of variation being constant at approximately 0.2 in all cases. The "relative division time," which is the time to first division for individual cells expressed in terms of the cell size doubling time, was used as measure of the "work to be done" to prepare for cell division. Relative division times were greater for heat-stressed cells than for those growing under osmotic stress.