Low-virulence Strains of Toxoplasma gondii Result in Permanent Loss of Innate Fear of Cats in Mice, Even after Parasite Clearance

Low-virulence Strains of Toxoplasma gondii Result in Permanent Loss of Innate Fear of Cats in Mice, Even after Parasite Clearance
Wendy Marie Ingram, Leeanne M Goodrich, Ellen A Robey, Michael B Eisen
(Submitted on 1 Apr 2013)

Toxoplasma gondii chronic infection in rodent secondary hosts has been reported to lead to a loss of innate, hard-wired fear toward cats, its primary host. However the generality of this response across T. gondii strains and the underlying mechanism for this pathogenmediated behavioral change remain unknown. To begin exploring these questions, we evaluated the effects of infection with isolates from the three major North American clonal lineages of T. gondii. Using an hour-long open field activity assay optimized for this purpose, we measured mouse aversion toward predator and non-predator urines. We show that loss of innate aversion of cat urine is a general trait caused by infection with all three major clonal lineages of parasite. Surprisingly, we found that infection with an attenuated Type I parasite results in sustained loss of fear at times post infection when neither parasite nor ongoing brain inflammation were detectable. This suggests that T. gondii-mediated interruption of mouse innate aversion of cats may occur during early acute infection in a permanent manner, not requiring persistence of parasite cysts or continuing brain inflammation.

A new approach to estimate directional genetic differentiation and asymmetric migration patterns

A new approach to estimate directional genetic differentiation and asymmetric migration patterns
Lisa Sundqvist, Martin Zackrisson, David Kleinhans
(Submitted on 30 Mar 2013)

In the field of population genetics measures of genetic differentiation are widely used to gather information on the structure and the amount of gene flow between populations. These indirect measures are based on a number of simplifying assumptions, for instance equal population size and symmetric migration. Structured populations with asymmetric migration patterns, frequently occur in nature and information about directional gene flow would here be of great interest. Nevertheless current measures of genetic differentiation cannot be used in such systems without violating the assumptions. To get information on asymmetric migration patterns from genetic data rather complex models using maximum likelihood or Bayesian approaches generally need to be applied. In such models a large number of parameters are estimated simultaneously and this involves complex optimization algorithms. We here introduce a new approach that intends to fill the gap between the complex approaches and the symmetric measures of genetic differentiation. Our approach makes it possible to calculate a directional component of genetic differentiation at low computational effort using any of the classical measures of genetic differentiation. The approach is based on defining a pool of migrants for any pair of populations and calculating measures for genetic differentiation between the populations and the respective pools. The directional measures of genetic differentiation can further be used to calculate asymmetric migration. The procedure is demonstrated with a simulated data set with known migration pattern. A comparison of the estimation results with the migration pattern used for simulation suggests, that our method captures relevant properties of migration patterns even at low migration frequencies and with few marker loci.

Detecting range expansions from genetic data

Detecting range expansions from genetic data

Benjamin M Peter, Montgomery Slatkin
(Submitted on 29 Mar 2013)

We propose a method that uses genetic data to test for the occurrence of a recent range expansion and to infer the location of the origin of the expansion. We introduce a statistic for pairs of populations \psi (the directionality index) that detects asymmetries in the two-dimensional allele frequency spectrum caused by the series of founder events that happen during an expansion. Such asymmetry arises because low frequency alleles tend to be lost during founder events, thus creating clines in the frequencies of surviving low-frequency alleles. Using simulations, we further show that \psi is more powerful for detecting range expansions than both F_{ST} and clines in heterozygosity. We illustrate the utility of \psi by applying it to a data set from modern humans and show how we can include more complicated scenarios such as multiple expansion origins or barriers to migration in the model.

Most viewed on Haldane’s Sieve, March 2013

The most viewed preprints on Haldane’s Sieve in March 2013 were:

The effects of transcription factor competition on gene regulation

The effects of transcription factor competition on gene regulation

Nicolae Radu Zabet, Boris Adryan
(Submitted on 27 Mar 2013)

We performed stochastic simulations of transcription factor (TF) molecules translocating by facilitated diffusion (a combination of 3D diffusion in the cytoplasm and 1D random walk on the DNA), and consider various abundances of cognate and non-cognate TFs to assess the influence of competitor molecules that also move along the DNA. We show that molecular crowding on the DNA always leads to longer times required by TF molecules to locate their target sites as well as to lower occupancy, which may confer a general mechanism to control gene activity levels globally. Finally, we show that crowding on the DNA may increase transcriptional noise through increased variability of the occupancy time of the target sites.

The influence of transcription factor competition on the relationship between occupancy and affinity

The influence of transcription factor competition on the relationship between occupancy and affinity

Nicolae Radu Zabet, Robert Foy, Boris Adryan
(Submitted on 27 Mar 2013)

Transcription factors (TFs) are proteins that bind to specific sites on the DNA and regulate gene activity. Identifying where TF molecules bind and how much time they spend on their target sites is key for understanding transcriptional regulation. It is usually assumed that the free energy of binding of a TF to the DNA (the affinity of the site) is highly correlated to the amount of time the TF remains bound (the occupancy of the site). However, knowing the binding energy is not sufficient to infer actual binding site occupancy. This mismatch between the occupancy predicted by the affinity and the observed occupancy may be caused by various factors, such as TF abundance, competition between TFs or the arrangement of the sites on the DNA. We investigated the relationship between the affinity of a TF for a set of binding sites and their occupancy. In particular, we considered the case of lac repressor (lacI) in E.coli and performed stochastic simulations of the TF dynamics on the DNA for various combinations of lacI abundance in competition with TFs that contribute to macromolecular crowding. Our results showed that for medium and high affinity sites, TF competition does not play a significant role in genomic occupancy, except in cases when the abundance of lacI is significantly increased or when a low-information content PWM was used. Nevertheless, for medium and low affinity sites, an increase in TF abundance (for both lacI or other molecules) leads to an increase in occupancy at several sites. Keywords: facilitated diffusion, Position Weight Matrix, thermodynamic equilibrium, motif information content, molecular crowding

SICLE: A high-throughput tool for extracting evolutionary relationships from phylogenetic trees

SICLE: A high-throughput tool for extracting evolutionary relationships from phylogenetic trees
Dan DeBlasio, Jennifer Wiscaver
(Submitted on 22 Mar 2013)

We present the phylogeny analysis software SICLE (Sister Clade Extractor), an easy to use, adaptable, and high-throughput tool to describe the nearest neighbors to a node of interest in a phylogenetic tree as well as the support value for the relationship. With SICLE it is possible to summarize the phylogenetic information produced by automated phylogenetic pipelines to rapidly identify and quantify the possible evolutionary relationships that merit further investigation. The program is a simple command line utility and is easy to adapt and implement in any phylogenetic pipeline. As a test case, we applied this new tool to published gene phylogenies to identify potential instances of horizontal gene transfer in Salinibacter ruber.

The Convergence of eQTL Mapping, Heritability Estimation and Polygenic Modeling: Emerging Spectrum of Risk Variation in Bipolar Disorder

The Convergence of eQTL Mapping, Heritability Estimation and Polygenic Modeling: Emerging Spectrum of Risk Variation in Bipolar Disorder
Eric R. Gamazon, Hae Kyung Im, Chunyu Liu, Members of the Bipolar Disorder Genome Study (BiGS) Consortium, Dan L. Nicolae, Nancy J. Cox
(Submitted on 25 Mar 2013)

It is widely held that a substantial genetic component underlies Bipolar Disorder (BD) and other neuropsychiatric disease traits. Recent efforts have been aimed at understanding the genetic basis of disease susceptibility, with genome-wide association studies (GWAS) unveiling some promising associations. Nevertheless, the genetic etiology of BD remains elusive with a substantial proportion of the heritability – which has been estimated to be 80% based on twin and family studies – unaccounted for by the specific genetic variants identified by large-scale GWAS. Furthermore, functional understanding of associated loci generally lags discovery. Studies we report here provide considerable support to the claim that substantially more remains to be gained from GWAS on the genetic mechanisms underlying BD susceptibility, and that a large proportion of the variation in disease risk may be uncovered through integrative functional genomic approaches. We combine recent analytic advances in heritability estimation and polygenic modeling and leverage recent technological advances in the generation of -omics data to evaluate the nature and scale of the contribution of functional classes of genetic variation to a relatively intractable disorder. We identified cis eQTLs in cerebellum and parietal cortex that capture more than half of the total heritability attributable to SNPs interrogated through GWAS and showed that eQTL-based heritability estimation is highly tissue-dependent. Our findings show that a much greater resolution may be attained than has been reported thus far on the number of common loci that capture a substantial proportion of the heritability to disease risk and that the functional nature of contributory loci may be clarified en masse.

An algebraic framework to sample the rearrangement histories of a cancer metagenome with double cut and join, duplication and deletion events

An algebraic framework to sample the rearrangement histories of a cancer metagenome with double cut and join, duplication and deletion events
Daniel R. Zerbino, Benedict Paten, Glenn Hickey, David Haussler
(Submitted on 22 Mar 2013)

Algorithms to study structural variants (SV) in whole genome sequencing (WGS) cancer datasets are currently unable to sample the entire space of rearrangements while allowing for copy number variations (CNV). In addition, rearrangement theory has up to now focused on fully assembled genomes, not on fragmentary observations on mixed genome populations. This affects the applicability of current methods to actual cancer datasets, which are produced from short read sequencing of a heterogeneous population of cells. We show how basic linear algebra can be used to describe and sample the set of possible sequences of SVs, extending the double cut and join (DCJ) model into the analysis of metagenomes. We also describe a functional pipeline which was run on simulated as well as experimental cancer datasets.

Natural selection reduced diversity on human Y chromosomes

Natural selection reduced diversity on human Y chromosomes
Melissa A. Wilson Sayres, Kirk E. Lohmueller, Rasmus Nielsen
(Submitted on 20 Mar 2013)

The human Y chromosome exhibits surprisingly low levels of genetic diversity. This could result from neutral processes if the effective population size of males is reduced relative to females due to a higher variance in the number of offspring from males than from females. Alternatively, selection acting on new mutations, and affecting linked neutral sites, could reduce variability on the Y chromosome. Here, using genome-wide analyses of X, Y, autosomal and mitochondrial DNA, in combination with extensive population genetic simulations, we show that low observed Y chromosome variability is not consistent with a purely neutral model. Instead, we show that models of purifying selection are consistent with observed Y diversity. Further, the number of sites estimated to be under purifying selection greatly exceeds the number of Y-linked coding sites, suggesting the importance of the highly repetitive ampliconic regions. Because the functional significance of the ampliconic regions is poorly understood, our findings should motivate future research in this area.