Correcting the Site Frequency Spectrum for Divergence-Based Ascertainment

Comparative genomics based on sequenced referenced genomes is essential to hypothesis generation and testing within population genetics. However, selection of candidate regions for further study on the basis of elevated or depressed divergence between species leads to a divergence-based ascertainment bias in the site frequency spectrum within selected candidate loci. Here, a method to correct this problem is developed that obtains maximum-likelihood estimates of the unascertained allele frequency distribution using numerical optimization. I show how divergence-based ascertainment may mimic the effects of natural selection and offer correction formulae for performing proper estimation into the strength of selection in candidate regions in a maximum-likelihood setting

Lineage Divergence and Historical Gene Flow in the Chinese Horseshoe Bat (Rhinolophus sinicus)

PMCID: PMC3581519This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Fluorescence correlation spectroscopy of the binding of nucleotide excision repair protein XPC-hHr23B with DNA substrates

The interaction of the nucleotide excision repair (NER) protein dimeric complex XPC-hHR23B, which is implicated in the DNA damage recognition step, with three Cy3.5 labeled 90-bp double-stranded DNA substrates (unmodified, with a central unpaired region, and cholesterol modified) and a 90-mer single-strand DNA was investigated in solution by fluorescence correlation spectroscopy. Autocorrelation functions obtained in the presence of an excess of protein show larger diffusion times (τ d) than for free DNA, indicating the presence of DNA-protein bound complexes. The fraction of DNA bound (θ), as a way to describe the percentage of protein bound to DNA, was directly estimated from FCS data. A significantly stronger binding capability for the cholesterol modified substrate (78% DNA bound) than for other double-stranded DNA substrates was observed, while the lowest affinity was found for the single-stranded DNA (27%). This is in accordance with a damage recognition role of the XPC protein. The similar affinity of XPC for undamaged and 'bubble' DNA sub

Adaptive Introgression across Species Boundaries in Heliconius Butterflies

It is widely documented that hybridisation occurs between many closely related species, but the importance of introgression in adaptive evolution remains unclear, especially in animals. Here, we have examined the role of introgressive hybridisation in transferring adaptations between mimetic Heliconius butterflies, taking advantage of the recent identification of a gene regulating red wing patterns in this genus. By sequencing regions both linked and unlinked to the red colour locus, we found a region that displays an almost perfect genotype by phenotype association across four species, H. melpomene, H. cydno, H. timareta, and H. heurippa. This particular segment is located 70 kb downstream of the red colour specification gene optix, and coalescent analysis indicates repeated introgression of adaptive alleles from H. melpomene into the H. cydno species clade. Our analytical methods complement recent genome scale data for the same region and suggest adaptive introgression has a crucial role in generating adaptive wing colour diversity in this group of butterflies

Estimation of Isolation Times of the Island Species in the Drosophila simulans Complex from Multilocus DNA Sequence Data

Background: The Drosophila simulans species complex continues to serve as an important model system for the study of new species formation. The complex is comprised of the cosmopolitan species, D. simulans, and two island endemics, D. mauritiana and D. sechellia. A substantial amount of effort has gone into reconstructing the natural history of the complex, in part to infer the context in which functional divergence among the species has arisen. In this regard, a key parameter to be estimated is the initial isolation time (t) of each island species. Loci in regions of low recombination have lower divergence within the complex than do other loci, yet divergence from D. melanogaster is similar for both classes. This might reflect gene flow of the lowrecombination loci subsequent to initial isolation, but it might also reflect differential effects of changing population size on the two recombination classes of loci when the low-recombination loci are subject to genetic hitchhiking or pseudohitchhiking Methodology/Principal Findings: New DNA sequence variation data for 17 loci corroborate the prior observation from 13 loci that DNA sequence divergence is reduced in genes of low recombination. Two models are presented to estimate t and other relevant parameters (substitution rate correction factors in lineages leading to the island species and, in the case of the 4-parameter model, the ratio of ancestral to extant effective population size) from the multilocus DNA sequence data. Conclusions/Significance: In general, it appears that both island species were isolated at about the same time, here estimated at,250,000 years ago. It also appears that the difference in divergence patterns of genes in regions of low an

Impact of Deep Coalescence on the Reliability of Species Tree Inference from Different Types of DNA Markers in Mammals

An important challenge for phylogenetic studies of closely related species is the existence of deep coalescence and gene tree heterogeneity. However, their effects can vary between species and they are often neglected in phylogenetic analyses. In addition, a practical problem in the reconstruction of shallow phylogenies is to determine the most efficient set of DNA markers for a reliable estimation. To address these questions, we conducted a multilocus simulation study using empirical values of nucleotide diversity and substitution rates obtained from a wide range of mammals and evaluated the performance of both gene tree and species tree approaches to recover the known speciation times and topological relationships. We first show that deep coalescence can be a serious problem, more than usually assumed, for the estimation of speciation times in mammals using traditional gene trees. Furthermore, we tested the performance of different sets of DNA markers in the determination of species trees using a coalescent approach. Although the best estimates of speciation times were obtained, as expected, with the use of an increasing number of nuclear loci, our results show that similar estimations can be obtained with a much lower number of genes and the incorporation of a mitochondrial marker, with its high information content. Thus, the use of the combined information of both nuclear and mitochondrial markers in a species tree framework is the most efficient option to estimate recent speciation times and, consequently, the underlying species tree

Correlated Evolution of Nearby Residues in Drosophilid Proteins

Here we investigate the correlations between coding sequence substitutions as a function of their separation along the protein sequence. We consider both substitutions between the reference genomes of several Drosophilids as well as polymorphisms in a population sample of Zimbabwean Drosophila melanogaster. We find that amino acid substitutions are “clustered” along the protein sequence, that is, the frequency of additional substitutions is strongly enhanced within ≈10 residues of a first such substitution. No such clustering is observed for synonymous substitutions, supporting a “correlation length” associated with selection on proteins as the causative mechanism. Clustering is stronger between substitutions that arose in the same lineage than it is between substitutions that arose in different lineages. We consider several possible origins of clustering, concluding that epistasis (interactions between amino acids within a protein that affect function) and positional heterogeneity in the strength of purifying selection are primarily responsible. The role of epistasis is directly supported by the tendency of nearby substitutions that arose on the same lineage to preserve the total charge of the residues within the correlation length and by the preferential cosegregation of neighboring derived alleles in our population sample. We interpret the observed length scale of clustering as a statistical reflection of the functional locality (or modularity) of proteins: amino acids that are near each other on the protein backbone are more likely to contribute to, and collaborate toward, a common subfunction

Complex aetiology of an apparently Mendelian form of Mental Retardation

<p>Abstract</p> <p>Background</p> <p>Mental Retardation is a common heterogeneous neurodevelopment condition, which causes are still largely elusive. It has been suggested that half of the phenotypic variation of intelligence is explained by genetic variation. And genetic or inherited factors indeed account for most of the cases of mental retardation with an identifiable cause. However, only a few autosomal genes have been mapped and identified to date. In this report, the genetic causes for an apparently recessive form of mental retardation, in a large nordern swedish pedigree, are investigated.</p> <p>Methods</p> <p>After extensive evaluation of the patients, which ruled out recognizable patterns of malformation and excluded known causes of MR, a comprehensive genome-wide linkage analysis, with 500 microsatellite markers, was performed in 24 members of this family. Additionally, a genome-wide copy number analysis, using an affimetrix 250 K SNP chip, was performed in this pedigree.</p> <p>Results</p> <p>No significant LOD score was found with either parametric and non-parametric linkage analysis. The highest scores are located at chromosomes 13, 15 and 17. Genome-wide copy number analysis identified no clear cause for the disorder; but rather, several variants were present in the family members, irrespective of their affected status.</p> <p>Conclusion</p> <p>These results suggest that mental retardation in this family, unlikely what was expected, has a heterogeneous aetiology; and that several lower effect genes variants might be involved. To demonstrate such effects, our family may be too small. This study also indicates that the ascertainment of the cause of MR may be challenging, and that a complex aetiology may be present even within a pedigree, constituting an additional obstacle for genetic counselling. Variants in genes involved in molecular mechanisms of cellular plasticity, in genes involved in the development of underlying neural architectures, and in genes involved in neurodevelopment and in the ongoing function of terminally differentiated neurons may underlie the phenotypic variation of intelligence and explain instances of intellectual impairment.</p

New Binding Mode to TNF-Alpha Revealed by Ubiquitin-Based Artificial Binding Protein

A variety of approaches have been employed to generate binding proteins from non-antibody scaffolds. Utilizing a beta-sheet of the human ubiquitin for paratope creation we obtained binding proteins against tumor necrosis factor (TNF)-alpha. The bioactive form of this validated pharmacological target protein is a non-covalently linked homo-trimer. This structural feature leads to the observation of a certain heterogeneity concerning the binding mode of TNF-alpha binding molecules, for instance in terms of monomer/trimer specificity. We analyzed a ubiquitin-based TNF-alpha binder, selected by ribosome display, with a particular focus on its mode of interaction. Using enzyme-linked immunosorbent assays, specific binding to TNF-alpha with nanomolar affinity was observed. In isothermal titration calorimetry we obtained comparable results regarding the affinity and detected an exothermic reaction with one ubiquitin-derived binding molecule binding one TNF-alpha trimer. Using NMR spectroscopy and other analytical methods the 1∶3 stoichiometry could be confirmed. Detailed binding analysis showed that the interaction is affected by the detergent Tween-20. Previously, this phenomenon was reported only for one other type of alternative scaffold-derived binding proteins – designed ankyrin repeat proteins – without further investigation. As demonstrated by size exclusion chromatography and NMR spectroscopy, the presence of the detergent increases the association rate significantly. Since the special architecture of TNF-alpha is known to be modulated by detergents, the access to the recognized epitope is indicated to be restricted by conformational transitions within the target protein. Our results suggest that the ubiquitin-derived binding protein targets a new epitope on TNF-alpha, which differs from the epitopes recognized by TNF-alpha neutralizing antibodies

Time in a Bottle: The Evolutionary Fate of Species Discrimination in Sibling Drosophila Species

Disadvantageous hybridization favors the evolution of prezygotic isolating behaviors, generating a geographic pattern of interspecific mate discrimination where members of different species drawn from sympatric populations exhibit stronger preference for members of their own species than do individuals drawn from allopatric populations. Geographic shifts in species' boundaries can relax local selection against hybridization; under such scenarios the fate of enhanced species preference is unknown. Lineages established from populations in the region of sympatry that have been maintained as single-species laboratory cultures represent cases where allopatry has been produced experimentally. Using such cultures dating from the 1950s, we assess how Drosophila pseudoobscura and D. persimilis mate preferences respond to relaxed natural selection against hybridization. We found that the propensity to hybridize generally declines with increasing time in experimental allopatry, suggesting that maintaining enhanced preference for conspecifics may be costly. However, our data also suggest a strong role for drift in determining mating preferences once secondary allopatry has been established. Finally, we discuss the interplay between populations in establishing the presence or absence of patterns consistent with reinforcement