A novel selection system for recombinational and mutational events within an intron of a eucaryotic gene.

In order to identify a poison sequence that might be useful in studying illegitimate recombination of mammalian cell chromosomes, several DNA segments were tested for their ability to interfere with gene expression when placed in an intron. A tRNA gene and its flanking sequences (267 bp) were shown to inhibit SV40 plaque formation 100-fold, when inserted into the intron in the T-antigen gene. Similarly, when the same DNA segment was placed in the second intron of the adenosine phosphoribosyl transferase (APRT) gene from CHO cells, it inhibited transformation of APRT-CHO cells 500-fold. These two tests indicated that the 267-bp DNA segment contained a poison sequence. The poison sequence did not affect replication since the replication of poisoned SV40 genomes was complemented by viable SV40 genomes and poisoned APRT genes were stably integrated into cell chromosomes. Cleavage of the poison sequence in the SV40 T-antigen intron by restriction enzymes indicated that the tRNA structural sequences and the 5' flanking sequences were not required for inhibition of SV40 plaque formation. Sequence analysis of viable mutant SV40, which arose after transfection of poisoned genomes, localized the poison sequence to a 35 bp segment immediately 3' of the tRNA structural sequences.     

A simple algorithm to infer gene duplication and speciation events on a gene tree

MOTIVATION: When analyzing protein sequences using sequence similarity searches, orthologous sequences (that diverged by speciation) are more reliable predictors of a new protein's function than paralogous sequences (that diverged by gene duplication), because duplication enables functional diversification. The utility of phylogenetic information in high-throughput genome annotation ('phylogenomics') is widely recognized, but existing approaches are either manual or indirect (e.g. not based on phylogenetic trees). Our goal is to automate phylogenomics using explicit phylogenetic inference. A necessary component is an algorithm to infer speciation and duplication events in a given gene tree. RESULTS: We give an algorithm to infer speciation and duplication events on a gene tree by comparison to a trusted species tree. This algorithm has a worst-case running time of O(n(2)) which is inferior to two previous algorithms that are approximately O(n) for a gene tree of sequences. However, our algorithm is extremely simple, and its asymptotic worst case behavior is only realized on pathological data sets. We show empirically, using 1750 gene trees constructed from the Pfam protein family database, that it appears to be a practical (and often superior) algorithm for analyzing real gene trees. AVAILABILITY: http://www.genetics.wustl.edu/eddy/forester.     

Leader of the pack: gene mapping in dogs and other model organisms

The domestic dog offers a unique opportunity to explore the genetic basis of disease, morphology and behaviour. We share many diseases with our canine companions, including cancer, diabetes and epilepsy, making the dog an ideal model organism for comparative disease genetics. Using newly developed resources, whole-genome association in dog breeds is proving to be exceptionally powerful. Here, we review the different trait-mapping strategies, some key biological findings emerging from recent studies and the implications for human health. We also discuss the development of similar resources for other vertebrate organisms.     

Stochastic search gene suggestion: a Bayesian hierarchical model for gene mapping

Mapping the genes for a complex disease, such as diabetes or rheumatoid arthritis (RA), involves finding multiple genetic loci that may contribute to the onset of the disease. Pairwise testing of the loci leads to the problem of multiple testing. Looking at haplotypes, or linear sets of loci, avoids multiple tests but results in a contingency table with sparse counts, especially when using marker loci with multiple alleles. We propose a hierarchical Bayesian model for case-parent triad data that uses a conditional logistic regression likelihood to model the probability of transmission to a diseased child. We define hierarchical prior distributions on the allele main effects to model the genetic dependencies present in the human leukocyte antigen (HLA) region of chromosome 6. First, we add a hierarchical level for model selection that accounts for both locus and allele selection. This allows us to cast the problem of identifying genetic loci relevant to the disease into a problem of Bayesian variable selection. Second, we attempt to include linkage disequilibrium as a covariance structure in the prior for model coefficients. We evaluate the performance of the procedure with some simulated examples and then apply our procedure to identifying genetic markers in the HLA region that influence risk for RA. Our software is available on the website http://www.epigenetic.org/Linkage/ssgs-public/.     

The best of both worlds: Phylogenetic eigenvector regression and mapping

Eigenfunction analyses have been widely used to model patterns of autocorrelation in time, space and phylogeny. In a phylogenetic context, Diniz-Filho et al. (1998) proposed what they called Phylogenetic Eigenvector Regression (PVR), in which pairwise phylogenetic distances among species are submitted to a Principal Coordinate Analysis, and eigenvectors are then used as explanatory variables in regression, correlation or ANOVAs. More recently, a new approach called Phylogenetic Eigenvector Mapping (PEM) was proposed, with the main advantage of explicitly incorporating a model-based warping in phylogenetic distance in which an Ornstein-Uhlenbeck (O-U) process is fitted to data before eigenvector extraction. Here we compared PVR and PEM in respect to estimated phylogenetic signal, correlated evolution under alternative evolutionary models and phylogenetic imputation, using simulated data. Despite similarity between the two approaches, PEM has a slightly higher prediction ability and is more general than the original PVR. Even so, in a conceptual sense, PEM may provide a technique in the best of both worlds, combining the flexibility of data-driven and empirical eigenfunction analyses and the sounding insights provided by evolutionary models well known in comparative analyses.     

BranchClust: a phylogenetic algorithm for selecting gene families

Background  

Automated methods for assembling families of orthologous genes include those based on sequence similarity scores and those based on phylogenetic approaches. The first are easy to automate but usually they do not distinguish between paralogs and orthologs or have restriction on the number of taxa. Phylogenetic methods often are based on reconciliation of a gene tree with a known rooted species tree; a limitation of this approach, especially in case of prokaryotes, is that the species tree is often unknown, and that from the analyses of single gene families the branching order between related organisms frequently is unresolved.     

Comparative mapping of a major aluminum tolerance gene in sorghum and other species in the poaceae

In several crop species within the Triticeae tribe of the grass family Poaceae, single major aluminum (Al) tolerance genes have been identified that effectively mitigate Al toxicity, a major abiotic constraint to crop production on acidic soils. However, the trait is quantitatively inherited in species within other tribes, and the possible ancestral relationships between major Al tolerance genes and QTL in the grasses remain unresolved. To help establish these relationships, we conducted a molecular genetic analysis of Al tolerance in sorghum and integrated our findings with those from previous studies performed in crop species belonging to different grass tribes. A single locus, AltSB, was found to control Al tolerance in two highly Al tolerant sorghum cultivars. Significant macrosynteny between sorghum and the Triticeae was observed for molecular markers closely linked to putatively orthologous Al tolerance loci present in the group 4 chromosomes of wheat, barley, and rye. However, AltSB was not located within the homeologous region of sorghum but rather mapped near the end of sorghum chromosome 3. Thus, AltSB not only is the first major Al tolerance gene mapped in a grass species that does not belong to the Triticeae, but also appears to be different from the major Al tolerance locus in the Triticeae. Intertribe map comparisons suggest that a major Al tolerance QTL on rice chromosome 1 is likely to be orthologous to AltSB, whereas another rice QTL on chromosome 3 is likely to correspond to the Triticeae group 4 Al tolerance locus. Therefore, this study demonstrates a clear evolutionary link between genes and QTL encoding the same trait in distantly related species within a single plant family.     

Phylogenetic and evolutionary relationships among torovirus field variants: evidence for multiple intertypic recombination events

Toroviruses (family Coronaviridae, order Nidovirales) are enveloped, positive-stranded RNA viruses that have been implicated in enteric disease in cattle and possibly in humans. Despite their potential veterinary and clinical relevance, little is known about torovirus epidemiology and molecular genetics. Here, we present the first study into the diversity among toroviruses currently present in European swine and cattle herds. Comparative sequence analysis was performed focusing on the genes for the structural proteins S, M, HE, and N, with fecal specimens serving as sources of viral RNA. Sequence data published for animal and human torovirus variants were included. Four genotypes, displaying 30 to 40% divergence, were readily distinguished, exemplified by bovine torovirus (BToV) Breda, porcine torovirus (PToV) Markelo, equine torovirus Berne, and the putative human torovirus. The ungulate toroviruses apparently display host species preference. In phylogenetic analyses, all PToV variants clustered, while the recent European BToVs mostly resembled the New World BToV variant Breda, identified 19 years ago. However, we found ample evidence for recurring intertypic recombination. All newly characterized BToV variants seem to have arisen from a genetic exchange, during which the 3' end of the HE gene, the N gene, and the 3' nontranslated region of a Breda virus-like parent had been swapped for those of PToV. Moreover, some PToV and BToV variants carried chimeric HE genes, which apparently resulted from recombination events involving hitherto unknown toroviruses. From these observations, the existence of two additional torovirus genotypes can be inferred. Toroviruses may be even more promiscuous than their closest relatives, the coronaviruses and arteriviruses.     

Epitope mapping using combinatorial phage-display libraries: a graph-based algorithm

A phage-display library of random peptides is a combinatorial experimental technique that can be harnessed for studying antibody–antigen interactions. In this technique, a phage peptide library is scanned against an antibody molecule to obtain a set of peptides that are bound by the antibody with high affinity. This set of peptides is regarded as mimicking the genuine epitope of the antibody's interacting antigen and can be used to define it. Here we present PepSurf, an algorithm for mapping a set of affinity-selected peptides onto the solved structure of the antigen. The problem of epitope mapping is converted into the task of aligning a set of query peptides to a graph representing the surface of the antigen. The best match of each peptide is found by aligning it against virtually all possible paths in the graph. Following a clustering step, which combines the most significant matches, a predicted epitope is inferred. We show that PepSurf accurately predicts the epitope in four cases for which the epitope is known from a solved antibody–antigen co-crystal complex. We further examine the capabilities of PepSurf for predicting other types of protein–protein interfaces. The performance of PepSurf is compared to other available epitope mapping programs.     

Comparative gene mapping: a valuable new tool for mammalian developmental studies

Technological advances in the 1970s encouraged the mapping of homologous gene loci in different mammalian species, including mouse and man. One hundred eighty-five homologous loci have now been mapped in these two species. Conservation of linkage is sufficient to identify substantial segments of the two genomes that have been left intact since their divergence from a common ancestor. The recognition of these conserved segments allows experimental manipulation of mouse chromosomes or chromosomal regions to produce models of human chromosomal anomalies of medical importance. Comparative gene mapping has been extended beyond mouse and man and the genomes of some species, including domestic cattle, appear to be more highly conserved relative to humans than the mouse. Such species may be particularly useful in providing models of human chromosomal anomalies that cannot be duplicated in laboratory mice.     

Dideoxy polymorphism scanning,a gene-based method for marker development for genetic linkage mapping

One of the fastest growing areas of biotechnology research today is marker-assisted breeding of crops. As a prerequisite to marker assisted breeding, genetic linkage maps are currently being developed for many species. For many purposes gene-based markers are the marker type of choice. The biggest problem in genetic linkage mapping using gene-based markers is the identification of polymorphisms between the parents of the population. To improve the efficiency of marker generation, we have developed a simple, and reasonable-cost method of polymorphism detection termed dideoxy polymorphism scanning. Since most of the time required to develop a gene-based linkage map is spent in identification of useful polymorphisms, this method will significantly shorten the time required for map generation and therefore reduce the overall cost.     

ALES: cell lineage analysis and mapping of developmental events

MOTIVATION: Animals build their bodies by altering the fates of cells. The way in which they do so is reflected in the topology of cell lineages and the fates of terminal cells. Cell lineages should, therefore, contain information about the molecular events that determined them. Here we introduce new tools for visualizing, manipulating, and extracting the information contained in cell lineages. Our tools enable us to analyze very large cell lineages, where previously analyses have only been carried out on cell lineages no larger than a few dozen cells. RESULTS: Ales (A Lineage Evaluation System) allows the display, evaluation and comparison of cell lineages with the aim of identifying molecular and cellular events underlying development. Ales introduces a series of algorithms that locate putative developmental events. The distribution of these predicted events can then be compared to gene expression patterns or other cellular characteristics. In addition, artificial lineages can be generated, or existing lineages modified, according to a range of models, in order to test hypotheses about lineage evolution. AVAILABILITY: The program can run on any operating system with a compliant Java 2 environment. Ales is free for academic use and can be downloaded from http://mbi.dkfz-heidelberg.de/mbi/research/cellsim/ales.