Cryptic MHC Polymorphism Revealed but Not Explained by Selection on the Class IIB Peptide-Binding Region

The immune genes of the major histocompatibility complex (MHC) are characterized by extraordinarily high levels of nucleotide and haplotype diversity. This variation is maintained by pathogen-mediated balancing selection that is operating on the peptide-binding region (PBR). Several recent studies have found, however, that some populations possess large clusters of alleles that are translated into virtually identical proteins. Here, we address the question of how this nucleotide polymorphism is maintained with little or no functional variation for selection to operate on. We investigate circa 750–850 bp of MHC class II DAB genes in four wild populations of the guppy Poecilia reticulata. By sequencing an extended region, we uncovered 40.9% more sequences (alleles), which would have been missed if we had amplified the exon 2 alone. We found evidence of several gene conversion events that may have homogenized sequence variation. This reduces the visible copy number variation (CNV) and can result in a systematic underestimation of the CNV in studies of the MHC and perhaps other multigene families. We then focus on a single cluster, which comprises 27 (of a total of 66) sequences. These sequences are virtually identical and show no signal of selection. We use microsatellites to reconstruct the populations' demography and employ simulations to examine whether so many similar nucleotide sequences can be maintained in the populations. Simulations show that this variation does not behave neutrally. We propose that selection operates outside the PBR, for example, on linked immune genes or on the “sheltered load” that is thought to be associated to the MHC. Future studies on the MHC would benefit from extending the amplicon size to include polymorphisms outside the exon with the PBR. This may capture otherwise cryptic haplotype variation and CNV, and it may help detect other regions in the MHC that are under selection

Bayesian co-estimation of selfing rate and locus-specific mutation rates for a partially selfing population

We present a Bayesian method for characterizing the mating system of populations reproducing through a mixture of self-fertilization and random outcrossing. Our method uses patterns of genetic variation across the genome as a basis for inference about pure hermaphroditism, androdioecy, and gynodioecy. We extend the standard coalescence model to accommodate these mating systems, accounting explicitly for multilocus identity disequilibrium, inbreeding depression, and variation in fertility among mating types. We incorporate the Ewens Sampling Formula (ESF) under the infinite-alleles model of mutation to obtain a novel expression for the likelihood of mating system parameters. Our Markov chain Monte Carlo (MCMC) algorithm assigns locus-specific mutation rates, drawn from a common mutation rate distribution that is itself estimated from the data using a Dirichlet Process Prior (DPP) model. Among the parameters jointly inferred are the population-wide rate of self-fertilization, locus-specific mutation rates, and the number of generations since the most recent outcrossing event for each sampled individual

On the Evolutionary Modification of Self-Incompatibility: Implications of Partial Clonality for Allelic Diversity and Genealogical Structure

Experimental investigations of homomorphic self-incompatibility (SI) have revealed an unanticipated level of complexity in its expression, permitting fine regulation over the course of a lifetime or a range of environmental conditions. Many flowering plants express some level of clonal reproduction, and phylogenetic analyses suggest that clonality evolves in a correlated fashion with SI in Solanum (Solanaceae). Here, we use a diffusion approximation to explore the effects on the evolutionary dynamics of SI of vegetative propagation with SI restricted to reproduction through seed. While clonality reduces the strength of frequency-dependent selection maintaining S-allele diversity, much of the great depth typical of S-allele genealogies is preserved. Our results suggest that clonality can play an important role in the evolution of SI systems, and may afford insight into unexplained features of allele genealogies in the Solanaceae

Hierarchical selection theory and sex ratios I. General solutions for structured populations

Models of sex-ratio evolution in structured populations are derived withG. R. Price's covariance form for the hierarchical analysis of natural selection (1970, Nature 227, 520-521). Previous work on competition among related males for mates (local mate competition), competition among related females for a limiting resource (local resource competition), inbreeding, group selection, and asymmetry of genetic inheritance between males and females, are subsumed under a general formulation for sex-ratio biases in structured populations. I found that the evolutionarily stable strategy sex ratio (males:females) for diploids is 1 - [varrho]m:1 - [varrho]f, where [varrho]m is the regression coefficient of relatedness of the controlling genotypes on males competing for mates, [varrho]f is the regression of controlling genotypes on females that compete for a fixed, limiting resource, and there is no inbreeding. For inbreeding and no competition among females, the evolutionarily stable strategy is 1 - [varrho]m:1 + [varrho]mf, where [varrho]mf is the regression of controlling genotypes on females' mates.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/26154/1/0000231.pd

A Bayesian Approach to Inferring Rates of Selfing and Locus-Specific Mutation.

We present a Bayesian method for characterizing the mating system of populations reproducing through a mixture of self-fertilization and random outcrossing. Our method uses patterns of genetic variation across the genome as a basis for inference about reproduction under pure hermaphroditism, gynodioecy, and a model developed to describe the self-fertilizing killifish Kryptolebias marmoratus. We extend the standard coalescence model to accommodate these mating systems, accounting explicitly for multilocus identity disequilibrium, inbreeding depression, and variation in fertility among mating types. We incorporate the Ewens sampling formula (ESF) under the infinite-alleles model of mutation to obtain a novel expression for the likelihood of mating system parameters. Our Markov chain Monte Carlo (MCMC) algorithm assigns locus-specific mutation rates, drawn from a common mutation rate distribution that is itself estimated from the data using a Dirichlet process prior model. Our sampler is designed to accommodate additional information, including observations pertaining to the sex ratio, the intensity of inbreeding depression, and other aspects of reproduction. It can provide joint posterior distributions for the population-wide proportion of uniparental individuals, locus-specific mutation rates, and the number of generations since the most recent outcrossing event for each sampled individual. Further, estimation of all basic parameters of a given model permits estimation of functions of those parameters, including the proportion of the gene pool contributed by each sex and relative effective numbers

Sex-Specific Incompatibility Generates Locus-Specific Rates of Introgression Between Species

Disruption of interactions among ensembles of epistatic loci has been shown to contribute to reproductive isolation among various animal and plant species. Under the Bateson–Dobzhansky–Muller model, such interspecific incompatibility arises as a by-product of genetic divergence in each species, and the Orr–Turelli model indicates that the number of loci involved in incompatible interactions may “snowball” over time. We address the combined effect of multiple incompatibility loci on the rate of introgression at neutral marker loci across the genome. Our analysis extends previous work by accommodating sex specificity: differences between the sexes in the expression of incompatibility, in rates of crossing over between neutral markers and incompatibility loci, and in transmission of markers or incompatibility factors. We show that the evolutionary process at neutral markers in a genome subject to incompatibility selection is well approximated by a purely neutral process with migration rates appropriately scaled to reflect the influence of selection targeted to incompatibility factors. We confirm that in the absence of sex specificity and functional epistasis among incompatibility factors, the barrier to introgression induced by multiple incompatibility factors corresponds to the product of the barriers induced by the factors individually. A new finding is that barriers to introgression due to sex-specific incompatibility depart in general from multiplicativity. Our partitioning of variation in relative reproductive rate suggests that such departures derive from associations between sex and incompatibility and between sex and neutral markers. Concordant sex-specific incompatibility (for example, greater impairment of male hybrids or longer map lengths in females) induces lower barriers (higher rates of introgression) than expected under multiplicativity, and discordant sex-specific incompatibility induces higher barriers

The Dawn of Open Access to Phylogenetic Data

The scientific enterprise depends critically on the preservation of and open access to published data. This basic tenet applies acutely to phylogenies (estimates of evolutionary relationships among species). Increasingly, phylogenies are estimated from increasingly large, genome-scale datasets using increasingly complex statistical methods that require increasing levels of expertise and computational investment. Moreover, the resulting phylogenetic data provide an explicit historical perspective that critically informs research in a vast and growing number of scientific disciplines. One such use is the study of changes in rates of lineage diversification (speciation - extinction) through time. As part of a meta-analysis in this area, we sought to collect phylogenetic data (comprising nucleotide sequence alignment and tree files) from 217 studies published in 46 journals over a 13-year period. We document our attempts to procure those data (from online archives and by direct request to corresponding authors), and report results of analyses (using Bayesian logistic regression) to assess the impact of various factors on the success of our efforts. Overall, complete phylogenetic data for ~60% of these studies are effectively lost to science. Our study indicates that phylogenetic data are more likely to be deposited in online archives and/or shared upon request when: (1) the publishing journal has a strong data-sharing policy; (2) the publishing journal has a higher impact factor, and; (3) the data are requested from faculty rather than students. Although the situation appears dire, our analyses suggest that it is far from hopeless: recent initiatives by the scientific community -- including policy changes by journals and funding agencies -- are improving the state of affairs

Individual and population sex allocation patterns

A variety of sex allocation models is considered in which (i) the reproductive returns on investment in males differ from the returns on investment in females, (ii) the amounts of resources available for reproduction vary in the population, (iii) the costs of making male and female reproductive structures differ, and (iv) the conception sex ratio may be fixed and there may be an initial minimum investment per offspring. Results of these models include quantatitive predictions for both individual- and population-level sex allocation, an opportunity to study the magnitude of changes in predicted patterns as key variables change, and therefore an analysis of the robustness of Fisher's equal investment theory. One example is that Fisher's argument is extremely robust for high fecundity organisms, but, in low fecundity organisms, is sensitive to differences between the sexes in reproductive returns on investment per offspring, a situation that occurs in many vertebrates to which Fisher's theory is often applied. A second example is that individual- and population-level patterns often depend strongly on the distribution of resources available for reproduction among individuals in the population.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/26816/1/0000374.pd