Microsatellite mutation models: insights from a comparison of humans and chimpanzees

Using genomic data from homologous microsatellite loci of pure AC repeats in humans and chimpanzees, several models of microsatellite evolution are tested and compared using likelihood-ratio tests and the Akaike information criterion. A proportional-rate, linear-biased, one-phase model emerges as the best model. A focal length toward which the mutational and/or substitutional process is linearly biased is a crucial feature of microsatellite evolution. We find that two-phase models do not lead to a significantly better fit than their one-phase counterparts. The performance of models based on the fit of their stationary distributions to the empirical distribution of microsatellite lengths in the human genome is consistent with that based on the human-chimp comparison. Microsatellites interrupted by even a single point mutation exhibit a twofold decrease in their mutation rate when compared to pure AC repeats. In general, models that allow chimps to have a larger per-repeat unit slippage rate and/or a shorter focal length compared to humans give a better fit to the human-chimp data as well as the human genomic data.     

Interaction of the anaphase-promoting complex/cyclosome and proteasome protein complexes with multiubiquitin chain-binding proteins

Fission yeast Rhp23 and Pus1 represent two families of multiubiquitin chain-binding proteins that associate with the proteasome. We show that both proteins bind to different regions of the proteasome subunit Mts4. The binding site for Pus1 was mapped to a cluster of repetitive sequences also found in the proteasome subunit SpRpn2 and the anaphase-promoting complex/cyclosome (APC/C) subunit Cut4. The putative role of Pus1 as a factor involved in allocation of ubiquitinylated substrates for the proteasome is discussed.     

The role of individual spikes and spike patterns in population coding of stimulus location in rat somatosensory cortex

This report addresses the nature of population coding in sensory cortex by applying information theoretic analysis to data recorded simultaneously from neuron pairs located in primary somatosensory cortex of anaesthetised rats. We studied how cortical spike trains code for the location of a whisker stimulus on the rat's snout. We found that substantially more information was conveyed by 10 ms precision spike timing compared with that conveyed by the number of spikes counted over a 40 ms response interval. Most of this information was accounted for by the timing of individual spikes. In particular, it was the first post-stimulus spikes that were crucial. Spike patterns within individual cells played a smaller role; spike patterns across cells were negligible. This pattern of results was robust both to the exact nature of the stimulus set and to the precision at which spikes were binned.     

Mapping mutations on phylogenies

Mapping of mutations on a phylogeny has been a commonly used analytical tool in phylogenetics and molecular evolution. However, the common approaches for mapping mutations based on parsimony have lacked a solid statistical foundation. Here, I present a Bayesian method for mapping mutations on a phylogeny. I illustrate some of the common problems associated with using parsimony and suggest instead that inferences in molecular evolution can be made on the basis of the posterior distribution of the mappings of mutations. A method for simulating a mapping from the posterior distribution of mappings is also presented, and the utility of the method is illustrated on two previously published data sets. Applications include a method for testing for variation in the substitution rate along the sequence and a method for testing whether the d(N)/d(S) ratio varies among lineages in the phylogeny.     

MAXIMUM-LIKELIHOOD ESTIMATION OF POPULATION DIVERGENCE TIMES AND POPULATION PHYLOGENY IN MODELS WITHOUT MUTATION

In this paper we present a method for estimating population divergence times by maximum likelihood in models without mutation. The maximum-likelihood estimator is compared to a commonly applied estimator based on Wright's F_ST statistic. Simulations suggest that the maximum-likelihood estimator is less biased and has a lower variance than the F_ST-based estimator. The maximum-likelihood estimator provides a statistical framework for the analysis of population history given genetic data. We demonstrate how maximum-likelihood estimates of the branching pattern of divergence of multiple populations may be obtained. We also describe how the method may be applied to test hypotheses such as whether populations have maintained equal population sizes. We illustrate the method by applying it to two previously published sets of human restriction fragment length polymorphism (RFLP) data.     

Contrasting ecosystem vegetation response in global drylands under drying and wetting conditions

Increasing aridity is one major consequence of ongoing global climate change and is expected to cause widespread changes in key ecosystem attributes, functions, and dynamics. This is especially the case in naturally vulnerable ecosystems, such as drylands. While we have an overall understanding of past aridity trends, the linkage between temporal dynamics in aridity and dryland ecosystem responses remain largely unknown. Here, we examined recent trends in aridity over the past two decades within global drylands as a basis for exploring the response of ecosystem state variables associated with land and atmosphere processes (e.g., vegetation cover, vegetation functioning, soil water availability, land cover, burned area, and vapor-pressure deficit) to these trends. We identified five clusters, characterizing spatiotemporal patterns in aridity between 2000 and 2020. Overall, we observe that 44.5% of all areas are getting dryer, 31.6% getting wetter, and 23.8% have no trends in aridity. Our results show strongest correlations between trends in ecosystem state variables and aridity in clusters with increasing aridity, which matches expectations of systemic acclimatization of the ecosystem to a reduction in water availability/water stress. Trends in vegetation (expressed by leaf area index [LAI]) are affected differently by potential driving factors (e.g., environmental, and climatic factors, soil properties, and population density) in areas experiencing water-related stress as compared to areas not exposed to water-related stress. Canopy height for example, has a positive impact on trends in LAI when the system is stressed but does not impact the trends in non-stressed systems. Conversely, opposite relationships were found for soil parameters such as root-zone water storage capacity and organic carbon density. How potential driving factors impact dryland vegetation differently depending on water-related stress (or no stress) is important, for example within management strategies to maintain and restore dryland vegetation.     

Estimating admixture pedigrees of recent hybrids without a contiguous reference genome

The genome of recently admixed individuals or hybrids has characteristic genetic patterns that can be used to learn about their recent admixture history. One of these are patterns of interancestry heterozygosity, which can be inferred from SNP data from either called genotypes or genotype likelihoods, without the need for information on genomic location. This makes them applicable to a wide range of data that are often used in evolutionary and conservation genomic studies, such as low-depth sequencing mapped to scaffolds and reduced representation sequencing. Here we implement maximum likelihood estimation of interancestry heterozygosity patterns using two complementary models. We furthermore develop apoh (Admixture Pedigrees of Hybrids), a software that uses estimates of paired ancestry proportions to detect recently admixed individuals or hybrids, and to suggest possible admixture pedigrees. It furthermore calculates several hybrid indices that make it easier to identify and rank possible admixture pedigrees that could give rise to the estimated patterns. We implemented apoh both as a command line tool and as a Graphical User Interface that allows the user to automatically and interactively explore, rank and visualize compatible recent admixture pedigrees, and calculate the different summary indices. We validate the performance of the method using admixed family trios from the 1000 Genomes Project. In addition, we show its applicability on identifying recent hybrids from RAD-seq data of Grant's gazelle (Nanger granti and Nanger petersii) and whole genome low-depth data of waterbuck (Kobus ellipsiprymnus) which shows complex admixture of up to four populations.     

Variation in the Pattern of Nucleotide Substitution Across Sites

A model of nucleotide substitution that allows the transition/transversion rate bias to vary across sites was constructed. We examined the fit of this model using likelihood-ratio tests by analyzing 13 protein coding genes and 1 pseudogene. Likelihood-ratio testing indicated that a model that allows variation in the transition/transversion rate bias across sites provided a significant improvement in fit for most protein coding genes but not for the pseudogene. When the analysis was repeated with parameters estimated separately for first, second, and third codon positions, strong heterogeneity was uncovered for the first and second codon positions; the variation in the transition/transversion rate was generally weaker at the third codon position. The transition rate bias and branch lengths are underestimated when variation in the transition/transversion rate was not accommodated, suggesting that it may be important to accommodate variation in the pattern of nucleotide substitution for accurate estimation of evolutionary parameters. Received: 4 November 1997 / Accepted: 19 May 1998     

Human-Mediated Admixture and Selection Shape the Diversity on the Modern Swine (Sus scrofa) Y Chromosomes

Throughout its distribution across Eurasia, domestic pig (Sus scrofa) populations have acquired differences through natural and artificial selection, and have often interbred. We resequenced 80 Eurasian pigs from nine different Asian and European breeds; we identify 42,288 reliable SNPs on the Y chromosome in a panel of 103 males, among which 96.1% are newly detected. Based on these new data, we elucidate the evolutionary history of pigs through the lens of the Y chromosome. We identify two highly divergent haplogroups: one present only in Asia and one fixed in Europe but present in some Asian populations. Analyzing the European haplotypes present in Asian populations, we find evidence of three independent waves of introgression from Europe to Asia in last 200 years, agreeing well with the literature and historical records. The diverse European lineages were brought in China by humans and left significant imprints not only on the autosomes but also on the Y chromosome of geographically and genetically distinct Chinese pig breeds. We also find a general excess of European ancestry on Y chromosomes relative to autosomes in Chinese pigs, an observation that cannot be explained solely by sex-biased migration and genetic drift. The European Y haplotype is associated with leaner meat production, and we hypothesize that the European Y chromosome increased in frequency in Chinese populations due to artificial selection. We find evidence of Y chromosomal gene flow between Sumatran wild boar and Chinese pigs. Our results demonstrate how human-mediated admixture and selection shaped the distribution of modern swine Y chromosomes.