The Effect of Single Recombination Events on Coalescent Tree Height and Shape

The coalescent with recombination is a fundamental model to describe the genealogical history of DNA sequence samples from recombining organisms. Considering recombination as a process which acts along genomes and which creates sequence segments with shared ancestry, we study the influence of single recombination events upon tree characteristics of the coalescent. We focus on properties such as tree height and tree balance and quantify analytically the changes in these quantities incurred by recombination in terms of probability distributions. We find that changes in tree topology are often relatively mild under conditions of neutral evolution, while changes in tree height are on average quite large. Our results add to a quantitative understanding of the spatial coalescent and provide the neutral reference to which the impact by other evolutionary scenarios, for instance tree distortion by selective sweeps, can be compared.     

Identification of selective sweeps using a dynamically adjusted number of linked microsatellites

There is currently large interest in distinguishing the signatures of genetic variation produced by demographic events from those produced by natural selection. We propose a simple multilocus statistical test to identify candidate sites of selective sweeps with high power. The test is based on the variability profile measured in an array of linked microsatellites. We also show that the analysis of flanking markers drastically reduces the number of false positives among the candidates that are identified in a genomewide survey of unlinked loci and find that this property is maintained in many population-bottleneck scenarios. However, for a certain range of intermediately severe population bottlenecks we find genomic signatures that are very similar to those produced by a selective sweep. While in these worst-case scenarios the power of the proposed test remains high, the false-positive rate reaches values close to 50%. Hence, selective sweeps may be hard to identify even if multiple linked loci are analyzed. Nevertheless, the integration of information from multiple linked loci always leads to a considerable reduction of the false-positive rate compared to a genome scan of unlinked loci. We discuss the application of this test to experimental data from Drosophila melanogaster.     

CaGE: cardiac gene expression knowledgebase

CaGE is a Cardiac Gene Expression knowledgebase we have developed to facilitate the analysis of genes important to human cardiac function. CaGE integrates the functionality of the LocusLink database with data from several human cardiac expression libraries, phenotypic data from OMIM and data from large-scale microarray gene expression studies to create a knowledgebase of gene expression in human cardiac tissue. The knowledgebase is fully searchable via the web using several intuitive query interfaces. Results can be displayed in several concise easy to navigate formats. AVAILABILITY: CaGE is located at http://www.cage.wbmei.jhu.edu     

Applying genetic programming to the prediction of alternative mRNA splice variants

Genetic programming (GP) can be used to classify a given gene sequence as either constitutively or alternatively spliced. We describe the principles of GP and apply it to a well-defined data set of alternatively spliced genes. A feature matrix of sequence properties, such as nucleotide composition or exon length, was passed to the GP system "Discipulus." To test its performance we concentrated on cassette exons (SCE) and retained introns (SIR). We analyzed 27,519 constitutively spliced and 9641 cassette exons including their neighboring introns; in addition we analyzed 33,316 constitutively spliced introns compared to 2712 retained introns. We find that the classifier yields highly accurate predictions on the SIR data with a sensitivity of 92.1% and a specificity of 79.2%. Prediction accuracies on the SCE data are lower, 47.3% (sensitivity) and 70.9% (specificity), indicating that alternative splicing of introns can be better captured by sequence properties than that of exons.     

A machine learning approach for the prediction of protein surface loop flexibility

Proteins often undergo conformational changes when binding to each other. A major fraction of backbone conformational changes involves motion on the protein surface, particularly in loops. Accounting for the motion of protein surface loops represents a challenge for protein-protein docking algorithms. A first step in addressing this challenge is to distinguish protein surface loops that are likely to undergo backbone conformational changes upon protein-protein binding (mobile loops) from those that are not (stationary loops). In this study, we developed a machine learning strategy based on support vector machines (SVMs). Our SVM uses three features of loop residues in the unbound protein structures-Ramachandran angles, crystallographic B-factors, and relative accessible surface area-to distinguish mobile loops from stationary ones. This method yields an average prediction accuracy of 75.3% compared with a random prediction accuracy of 50%, and an average of 0.79 area under the receiver operating characteristic (ROC) curve using cross-validation. Testing the method on an independent dataset, we obtained a prediction accuracy of 70.5%. Finally, we applied the method to 11 complexes that involve members from the Ras superfamily and achieved prediction accuracy of 92.8% for the Ras superfamily proteins and 74.4% for their binding partners.     

How repetitive are genomes?

Background  

Genome sequences vary strongly in their repetitiveness and the causes for this are still debated. Here we propose a novel measure of genome repetitiveness, the index of repetitiveness, I _r, which can be computed in time proportional to the length of the sequences analyzed. We apply it to 336 genomes from all three domains of life.     

Functional Relevance of Improbable Antibody Mutations for HIV Broadly Neutralizing Antibody Development

1. Download high-res image (257KB)
2. Download full-size image

     

Anthropometric Indices and the Incidence of Hypertension: A Comparative Analysis

Objective: To investigate the association between several anthropometric measurements of obesity with the incidence of hypertension. Research Methods and Procedures: Participants were 592 individuals free of hypertension, selected at random from the community. In the baseline evaluation, they were submitted to completed measures of demographics, anthropometrics, blood pressure, and other risk factors for hypertension. Incident hypertension was defined by blood pressure equal or higher than 140/90 mm Hg or use of blood pressure‐lowering drugs. Results: During a mean follow‐up time of 5.6 ± 1.1 years, 127 developed hypertension. The hazard ratios for the development of hypertension, adjusted for age, baseline blood pressure, gender, and alcohol consumption, were 1.042 (p = 0.091) for BMI, 1.023 (p = 0.028) for waist circumference, 1.042 (p = 0.013) for waist‐to‐height ratio, 1.061 (p = 0.014) for waist‐to‐height² index, 1.079 (p = 0.022) for waist‐to‐height³ index, and 1.033 (p = 0.006) for the waist‐to‐hip ratio. Discussion: The correction of the circumference of waist for stature or hip circumference improves its performance in the prediction of the incidence of hypertension.     

Reversion and T Cell Escape Mutations Compensate the Fitness Loss of a CD8+ T Cell Escape Mutant in Their Cognate Transmitted/Founder Virus

Immune escape mutations that revert back to the consensus sequence frequently occur in newly HIV-1-infected individuals and have been thought to render the viruses more fit. However, their impact on viral fitness and their interaction with other immune escape mutations have not been evaluated in the background of their cognate transmitted/founder (T/F) viral genomes. To precisely determine the role of reversion mutations, we introduced reversion mutations alone or together with CD8⁺ T cell escape mutations in their unmodified cognate T/F viral genome and determined their impact on viral fitness in primary CD4⁺ T cells. Two reversion mutations, V247I and I64T, were identified in Gag and Tat, respectively, but neither had measurable effect on the fitness of their cognate T/F virus. The V247I and G248A mutations that were detected before and concurrently with the potent T cell escape mutation T242N, respectively, were selected by early T cell responses. The V247I or the G248A mutation alone partially restored the fitness loss caused by the T242N mutation. Together they could fully restore the fitness of the T242N mutant to the T/F level. These results demonstrate that the fitness loss caused by a T cell escape mutation could be compensated by preexisting or concurrent reversion and other T cell escape mutations. Our findings indicate that the overall viral fitness is modulated by the complex interplay among T cell escape, compensatory and reversion mutations to maintain the balance between immune escape and viral replication capacity.     

Bifurcations in haploid and diploid sequence space models

 Deterministic models of mutation and selection in the space of (binary) nucleotide-type sequences have been investigated for haploid populations during the past 25 years, and, recently, for diploid populations as well. These models, in particular their ‘error thresholds’, have mainly been analyzed by numerical methods and perturbation techniques. We consider them here by means of bifurcation theory, which improves our understanding of both equilibrium and dynamical properties. In a caricature obtained from the original model by neglecting back mutation to the favourable allele, the familiar error threshold of the haploid two-class model turns out to be a simple transcritical bifurcation, whereas its diploid counterpart exhibits an additional saddle node. This corresponds to a second error threshold. Three-class models with neutral spaces of unequal size introduce further features. Such are a global bifurcation in haploid populations, and simple examples of Hopf bifurcations (as predicted by Akin’s theorem) in the diploid case. Received 13 June 1995; received in revised form 26 July 1996