Matrix dimensions bias demographic inferences: implications for comparative plant demography

While the wealth of projection matrices in plant demography permits comparative studies, variation in matrix dimensions complicates interspecific comparisons. Collapsing matrices to a common dimension may facilitate such comparisons but may also bias the inferred demographic parameters. Here we examine how matrix dimension affects inferred demographic elasticities and how different collapsing criteria perform. We analyzed 13 x 13 matrices representing nine plant species, collapsing these matrices (i) into even 7 x 7, 5 x 5, 4 x 4, and 3 x 3 matrices and (ii) into 5 x 5 matrices using different criteria. Stasis and fecundity elasticities increased when matrix dimension was reduced, whereas those of progression and retrogression decreased. We suggest a collapsing criterion that minimizes dissimilarities between the original- and collapsed-matrix elasticities and apply it to 66 plant species to study how life span and growth form influence the relationship between matrix dimension and elasticities. Our analysis demonstrates that (i) projection matrix dimension has significant effects on inferred demographic parameters, (ii) there are better-performing methods than previously suggested for standardizing matrix dimension, and (iii) herbaceous perennial projection matrices are particularly sensitive to changes in matrix dimensionality. For comparative demographic studies, we recommend normalizing matrices to a common dimension by collapsing higher classes and leaving the first few classes unaltered.     

Inferences about human demography based on multilocus analyses of noncoding sequences

Data from 10 unlinked autosomal noncoding regions, resequenced in 15 individuals from each of three populations, were used in a multilocus analysis to test models of human demography. Each of the 10 regions consisted of approximately 2500 bp. The multilocus analysis, based on summary statistics (average and variance of Tajima's D and Fu and Li's D*), was used to test a family of models with recent population expansion. The African sample (Hausa of Cameroon) is compatible with a constant population size model and a range of models with recent expansion. For this population sample, we estimated confidence sets that showed the limited range of parameter values compatible with growth. For an exponential growth rate as low as 1 x 10(-3)/generation, population growth is unlikely to have started prior to 50,000 years ago. For higher growth rates, the onset of growth must be more recent. On the basis of the average value of Tajima's D, our sample from an Italian population was found to be incompatible with a constant population size model or any simple expansion model. In the Chinese sample, the variance of Tajima's D was too large to be compatible with the constant population size model or any simple expansion model.     

Relatedness and demography of African forest elephants: inferences from noninvasive fecal DNA analyses

African forest elephants (Loxodonta cyclotis) are genetically and morphologically distinct from their savannah counterparts, but their biology remains poorly understood. In this study, I use noninvasive fecal DNA analyses to examine the relatedness structure and historical demography of forest elephants at 2 sites in SW Gabon, central Africa. Pairwise relatedness values calculated between 162 elephant individuals genotyped at 8 microsatellite loci were significantly higher within spatially associated dung piles than between random pairings for one site. First- and second-order relatives were most commonly detected among dung piles from adult female pairs and adult females and juveniles. Pairwise relatedness estimates suggested that, like savannah elephants, forest groups are largely composed of adult females, their sisters, and juvenile offspring. Associations between males, and groups containing juveniles from multiple related females, were detected but at much lower frequency. Analysis of mitochondrial d-loop sequences from 70 elephant individuals identified 2 haplogroups in SW Gabon.     

Drawing inferences about the coancestry coefficient

The coancestry coefficient, also known as the population structure parameter, is of great interest in population genetics. It can be thought of as the intraclass correlation of pairs of alleles within populations and it can serve as a measure of genetic distance between populations. For a general class of evolutionary models it determines the distribution of allele frequencies among populations. Under more restrictive models it can be regarded as the probability of identity by descent of any pair of alleles at a locus within a random mating population. In this paper we review estimation procedures that use the method of moments or are maximum likelihood under the assumption of normally distributed allele frequencies. We then consider the problem of testing hypotheses about this parameter. In addition to parametric and non-parametric bootstrap tests we present an asymptotically-distributed chi-square test. This test reduces to the contingency-table test for equal sample sizes across populations. Our new test appears to be more powerful than previous tests, especially for loci with multiple alleles. We apply our methods to HapMap SNP data to confirm that the coancestry coefficient for humans is strictly positive.     

Structural inferences for Cholera toxin mutations in Vibrio cholerae

Cholera is a global disease that has persisted for millennia. The cholera toxin (CT) from Vibrio cholerae is responsible for the clinical symptoms of cholera. This toxin is a hetero-hexamer (AB(5)) complex consisting of a subunit A (CTA) with a pentamer (B(5)) of subunit B (CTB). The importance of the AB(5) complex for pathogenesis is established for the wild type O1 serogroup using known structural and functional data. However, its role is not yet documented in other known serogroups harboring sequence level residue mutations. The sequences for the toxin from different serogroups are available in GenBank (release 177). Sequence analysis reveals mutations at several sequence positions in the toxin across serogroups. Therefore, it is of interest to locate the position of these mutations in the AB(5) structure to infer complex assembly for its functional role in different serogroups. We show that mutations in the CTA are at the solvent exposed regions of the AB(5) complex, whereas those in the CTB are at the CTB/CTB interface of the homo-pentamer complex. Thus, the role of mutations at the CTB/CTB interface for B(5) complex assembly is implied. It is observed that these mutations are often non-synonymous (e.g. polar to non-polar or vice versa). The formation of the AB(5) complex involves inter-subunit residue-residue interactions at the protein-protein interfaces. Hence, these mutations, at the structurally relevant positions, are of importance for the understanding of pathogenesis by several serogroups. This is also of significance in the improvement of recombinant CT protein complex analogs for vaccine design and their use against multiple serogroups.     

Properties of ethylmethane sulfonate-induced mutations affecting life-history traits in Caenorhabditis elegans and inferences about bivariate distributions of mutation effects

The homozygous effects of ethylmethane sulfonate (EMS)-induced mutations in Caenorhabditis elegans are compared across life-history traits. Mutagenesis has a greater effect on early than late reproductive output, since EMS-induced mutations tend to cause delayed reproduction. Mutagenesis changes the mean and variance of longevity much less than reproductive output traits. Mutations that increase total or early productivity are not detected, but the net effect of mutations is to increase and decrease late productivity to approximately equal extents. Although most mutations decrease longevity, a mutant line with increased longevity was found. A flattening of mortality curves with age is noted, particularly in EMS lines. We infer that less than one-tenth of mutations that have fitness effects in natural conditions are detected in the laboratory, and such mutations have moderately large effects ( approximately 20% of the mean). Mutational correlations for life-history traits are strong and positive. Correlations between early or late productivity and longevity are of similar magnitude. We develop a maximum-likelihood procedure to infer bivariate distributions of mutation effects. We show that strong mutation-induced genetic correlations do not necessarily imply strong directional correlations between mutational effects, since correlation is also generated by lines carrying different numbers of mutations.     

Rational inferences about departures from Hardy-Weinberg equilibrium

Previous studies have explored the use of departure from Hardy-Weinberg equilibrium (DHW) for fine mapping Mendelian disorders and for general fine mapping. Other studies have used Hardy-Weinberg tests for genotyping quality control. To enable investigators to make rational decisions about whether DHW is due to genotyping error or to underlying biology, we developed an analytic framework and software to determine the parameter values for which DHW might be expected for common diseases. We show analytically that, for a general disease model, the difference between population and Hardy-Weinberg expected genotypic frequencies (delta) at the susceptibility locus is a function of the susceptibility-allele frequency (q), heterozygote relative risk (beta), and homozygote relative risk (gamma). For unaffected control samples, is a function of risk in nonsusceptible homozygotes (alpha), the population prevalence of disease (KP), q, beta, and gamma. We used these analytic functions to calculate and the number of cases or controls needed to detect DHW for a range of genetic models consistent with common diseases (1.1 < or = gamma < or = 10 and 0.005 < or = KP < or = 0.2). Results suggest that significant DHW can be expected in relatively small samples of patients over a range of genetic models. We also propose a goodness-of-fit test to aid investigators in determining whether a DHW observed in the context of a case-control study is consistent with a genetic disease model. We illustrate how the analytic framework and software can be used to help investigators interpret DHW in the context of association studies of common diseases.     

A comprehensive model of mutations affecting fitness and inferences for Arabidopsis thaliana

As the ultimate source of genetic variation, spontaneous mutation is essential to evolutionary change. Theoretical studies over several decades have revealed the dependence of evolutionary consequences of mutation on specific mutational properties, including genomic mutation rates, U, and the effects of newly arising mutations on individual fitness, s. The recent resurgence of empirical effort to infer these properties for diverse organisms has not achieved consensus. Estimates, which have been obtained by methods that assume mutations are unidirectional in their effects on fitness, are imprecise. Both because a general approach must allow for occurrence of fitness-enhancing mutations, even if these are rare, and because recent evidence demands it, we present a new method for inferring mutational parameters. For the distribution of mutational effects, we retain Keightley's assumption of the gamma distribution, to take advantage of the flexibility of its shape. Because the conventional gamma is one sided, restricting it to unidirectional effects, we include an additional parameter, rho, as an amount it is displaced from zero. Estimation is accomplished by Markov chain Monte Carlo maximum likelihood. Through a limited set of simulations, we verify the accuracy of this approach. We apply it to analyze data on two reproductive fitness components from a 17-generation mutation-accumulation study of a Columbia accession of Arabidopsis thaliana in which 40 lines sampled in three generations were assayed simultaneously. For these traits, U approximately/= 0.1-0.2, with distributions of mutational effects broadly spanning zero, such that roughly half the mutations reduce reproductive fitness. One evolutionary consequence of these results is lower extinction risks of small populations of A. thaliana than expected from the process of mutational meltdown. A comprehensive view of the evolutionary consequences of mutation will depend on quantitatively accounting for fitness-enhancing, as well as fitness-reducing, mutations.     

Positive selection versus demography: evolutionary inferences based on an unusual haplotype structure in Drosophila simulans

Coalescent simulations were used to investigate the possible role of population subdivision and history in shaping nucleotide variation in a recombining 88-kb genomic fragment of Drosophila simulans displaying an unusual large-scale haplotype structure. The multilocus analysis, based on summary statistics using specific demographic null models under recombination, indicates that the observed levels of linkage disequilibrium differed significantly from the values expected under different bottleneck and population admixture scenarios. These results indicate that demography alone may not account for the observed pattern of variation and support the previous claim that the data are better described by a model in which an adaptive mutation has not yet gone to fixation.     

1 Taxon sampling and inferences about diatom phylogeny

Proper taxon sampling is one of the greatest challenges to understanding phylogenetic relationships, perhaps as important as choice of optimality criterion or data type. This has been demonstrated in diatoms where centric diatoms may either be strongly supported as monophyletic or paraphyletic when analyzing SSU rDNA sequences using the same optimality criterion. The effect of ingroup and outgroup taxon sampling on relationships of diatoms is explored for diatoms as a whole and for the order Thalassiosirales. In the latter case, SSU rDNA and rbcL sequence data result in phylogenetic relationships that appear to be strongly incongruent with morphology and broadly incongruent with the fossil record. For example, Cyclotella stelligera Cleve & Grunow behaves like a rogue taxon, jumping from place to place throughout the tree. Morphological data place C. stelligera near the base of the freshwater group as sister to the extinct genus Mesodictyon Theriot and Bradbury, suggesting that it is an old, long branch that might be expected to “misbehave” in poorly sampled trees. Cyclotella stelligera and C. bodanica Grunow delimit the diameter of morphological diversity in Cyclotella, so increased sampling of intermediate taxa will be critical to resolving this part of the tree. Morphology is sampled for a much greater number of taxa and many transitional states of putative synapomorphies seem to suggest a robust morphological hypothesis. The Thalassiosirales are unstable with regards to taxon sampling in the genetic data, suggesting that perhaps the morphological hypothesis is (for now) preferable.     

The Leinster and Cobbold indices improve inferences about microbial diversity

An increasing number of ecological studies compare the diversity of microbial taxa along environmental gradients or between imposed treatments. Estimates are often based on analysis-of-variance of taxon-richness inferred from pyrosequencing data. We conducted a reanalysis of three 454-pyrosequencing studies on arbuscular-mycorrhizal-fungal diversity to evaluate the suitability of using the Leinster and Cobbold diversity-indices (LCdis) to assess diversity. We expected that the potential of LCdis to consider phylogenic relationships could resolve problems arising from ambiguous species-delineation in microbial-systems. Our reanalysis showed that comparisons between studies differing considerably in sequencing depth may be risky. Moreover, we show that LCdis not only reproduce the results of analyses of variance but can also resolve issues connected to variation in sequence read number, while additionally representing a less conservative metric of diversity than analysis-of-variance of taxa-richness. Based on these results we advocate the use of inclusive diversity indices in ecological studies targeting microbial communities.     

Microsatellite mutations in litters of the Australian lizard Egernia stokesii

Seven tetranucleotide (AAAG) loci were analysed in a population (including 19 litters) of the Australian lizard Egernia stokesii. In an examination of 76 offspring we observed 13 mutations involving five loci. Two of our loci were highly mutable, with observed mutation rates of 2.7% and 4.2%, representing some of the highest mutations rates reported so far. A high proportion of mutations (46.2%) could not be assigned to changes involving only a single repeat, suggesting that mutations in at least two of the loci follow a multistep model. There was no significant bias of mutations leading to an increase or decrease in allele size; however, all multistep mutations involved a loss of repeats. These results add to increasing evidence casting doubt on microsatellite mutations being primarily single step mutations.     

Microsatellite mutations in the germline: implications for evolutionary inference

Microsatellite DNA sequences mutate at rates several orders of magnitude higher than that of the bulk of DNA. Such high rates mean that spontaneous mutations that form new-length variants can realistically be seen in pedigree analysis. Data on observed mutation events from various organisms are now accumulating, allowing inferences on DNA sequence evolution to be made through an unusually direct approach. Here I discuss and integrate microsatellite mutation data in an evolutionary context. A striking feature of the mutation process is that it seems highly heterogeneous, with distinct differences between species, repeat types, loci and alleles. Age and sex also affect the mutation rate. Within genomes at equilibrium, the microsatellite-length distribution is a delicate balance between biased mutation processes and point mutations acting towards the decay of repetitive DNA. Indeed, simple repeats do not evolve simply.     

Ecological inferences about marine mammals from passive acoustic data

Monitoring on the basis of sound recordings, or passive acoustic monitoring, can complement or serve as an alternative to real-time visual or aural monitoring of marine mammals and other animals by human observers. Passive acoustic data can support the estimation of common, individual-level ecological metrics, such as presence, detection-weighted occupancy, abundance and density, population viability and structure, and behaviour. Passive acoustic data also can support estimation of some community-level metrics, such as species richness and composition. The feasibility of estimation and certainty of estimates is highly context dependent, and understanding the factors that affect the reliability of measurements is useful for those considering whether to use passive acoustic data. Here, we review basic concepts and methods of passive acoustic sampling in marine systems that often are applicable to marine mammal research and conservation. Our ultimate aim is to facilitate collaboration among ecologists, bioacousticians, and data analysts. Ecological applications of passive acoustics require one to make decisions about sampling design, which in turn requires consideration of sound propagation, sampling of signals, and data storage. One also must make decisions about signal detection and classification and evaluation of the performance of algorithms for these tasks. Investment in the research and development of systems that automate detection and classification, including machine learning, are increasing. Passive acoustic monitoring is more reliable for detection of species presence than for estimation of other species-level metrics. Use of passive acoustic monitoring to distinguish among individual animals remains difficult. However, information about detection probability, vocalisation or cue rate, and relations between vocalisations and the number and behaviour of animals increases the feasibility of estimating abundance or density. Most sensor deployments are fixed in space or are sporadic, making temporal turnover in species composition more tractable to estimate than spatial turnover. Collaborations between acousticians and ecologists are most likely to be successful and rewarding when all partners critically examine and share a fundamental understanding of the target variables, sampling process, and analytical methods.     

Integron diversity in heavy-metal-contaminated mine tailings and inferences about integron evolution

Integrons are horizontal gene transfer (HGT) systems containing elements necessary for site-specific recombination and expression of foreign DNA. The overall phylogenetic distribution of integrons and range of genes that can be transferred by integrons are unknown. This report contains an exploration of integrons in an environmental microbial community and an investigation of integron evolution. First, using culture-independent techniques, we explored the diversity of integrons and integron-transferred genes in heavy-metal-contaminated mine tailings. Using degenerate primers, we amplified integron integrase genes from the tailings. We discovered 14 previously undescribed integrase genes, including six novel gene lineages. In addition, we found 11 novel gene cassettes in this sample. One of the gene cassettes that we sequenced is similar to a gene that codes for a step in a pathway for nitroaromatic catabolism, a group of compounds associated with mining activity. This suggests that integrons may be important for gene transfer in response to selective pressures other than the presence of antibiotics. We also investigated the evolution of integrons by statistically comparing the phylogenies of 16S rRNA and integrase genes from the same organisms, using sequences from GenBank and various sequencing projects. We found significant differences between the organismal (16S rRNA) and integrase trees, and we suggest that these differences may be due to HGT.     

Divergence,demography and gene loss along the human lineage

Genomic DNA sequences are an irreplaceable source for reconstructing the vanished past of living organisms. Based on updated sequence data, this paper summarizes our studies on species divergence time, ancient population size and functional loss of genes in the primate lineage leading to modern humans (Homo sapiens sapiens). The inter- and intraspecific comparisons of DNA sequences suggest that the human lineage experienced a rather severe bottleneck in the Middle Pleistocene, throughout which period the subdivided African population played a predominant role in shaping the genetic architecture of modern humans. Also, published and newly identified human-specific pseudogenes (HSPs) are enumerated in order to infer their significance for human evolution. Of the 121 candidate genes obtained, authentic HSPs turn out to comprise only 25 olfactory receptor genes, four T cell receptor genes and nine other genes. The fixation of HSPs has been too rare over the past 6–7 Myr to account for species differences between humans and chimpanzees.     

Disentangling the effects of demography and selection in human history

Demographic events affect all genes in a genome, whereas natural selection has only local effects. Using publicly available data from 151 loci sequenced in both European-American and African-American populations, we attempt to distinguish the effects of demography and selection. To analyze large sets of population genetic data such as this one, we introduce "Perlymorphism," a Unix-based suite of analysis tools. Our analyses show that the demographic histories of human populations can account for a large proportion of effects on the level and frequency of variation across the genome. The African-American population shows both a higher level of nucleotide diversity and more negative values of Tajima's D statistic than does a European-American population. Using coalescent simulations, we show that the significantly negative values of the D statistic in African-Americans and the positive values in European-Americans are well explained by relatively simple models of population admixture and bottleneck, respectively. Working within these nonequilibrium frameworks, we are still able to show deviations from neutral expectations at a number of loci, including ABO and TRPV6. In addition, we show that the frequency spectrum of mutations--corrected for levels of polymorphism--is correlated with recombination rate only in European-Americans. These results are consistent with repeated selective sweeps in non-African populations, in agreement with recent reports using microsatellite data.