Interpreting the genomic landscape of speciation: a road map for finding barriers to gene flow

Speciation, the evolution of reproductive isolation among populations, is continuous, complex, and involves multiple, interacting barriers. Until it is complete, the effects of this process vary along the genome and can lead to a heterogeneous genomic landscape with peaks and troughs of differentiation and divergence. When gene flow occurs during speciation, barriers restricting gene flow locally in the genome lead to patterns of heterogeneity. However, genomic heterogeneity can also be produced or modified by variation in factors such as background selection and selective sweeps, recombination and mutation rate variation, and heterogeneous gene density. Extracting the effects of gene flow, divergent selection and reproductive isolation from such modifying factors presents a major challenge to speciation genomics. We argue one of the principal aims of the field is to identify the barrier loci involved in limiting gene flow. We first summarize the expected signatures of selection at barrier loci, at the genomic regions linked to them and across the entire genome. We then discuss the modifying factors that complicate the interpretation of the observed genomic landscape. Finally, we end with a road map for future speciation research: a proposal for how to account for these modifying factors and to progress towards understanding the nature of barrier loci. Despite the difficulties of interpreting empirical data, we argue that the availability of promising technical and analytical methods will shed further light on the important roles that gene flow and divergent selection have in shaping the genomic landscape of speciation.     

Three aldo-keto reductases of the yeast Saccharomyces cerevisiae

Saccharomyces cerevisiae is an industrially important yeast, which is also used extensively as a model eukaryote. The S. cerevisiae genome has been sequenced in its entirety and therefore represents an ideal organism in which to carry out functional analysis of genes. We have identified several open reading frames in the S. cerevisiae genome which show significant similarity to members of the aldo-keto reductase superfamily. The physiological roles of these gene products have not been previously determined, but their similarity to other enzymes suggests they may perform roles in carbohydrate metabolism and detoxification pathways. Cloning and expression of three of these enzymes has allowed their substrate specificities to be determined. Expression profiling and gene disruption analysis will allow potential roles for these enzymes within the cell to be examined.     

Mapping the Human Reference Genome’s Missing Sequence by Three-Way Admixture in Latino Genomes

A principal obstacle to completing maps and analyses of the human genome involves the genome’s “inaccessible” regions: sequences (often euchromatic and containing genes) that are isolated from the rest of the euchromatic genome by heterochromatin and other repeat-rich sequence. We describe a way to localize these sequences by using ancestry linkage disequilibrium in populations that derive ancestry from at least three continents, as is the case for Latinos. We used this approach to map the genomic locations of almost 20 megabases of sequence unlocalized or missing from the current human genome reference (NCBI Genome GRCh37)—a substantial fraction of the human genome’s remaining unmapped sequence. We show that the genomic locations of most sequences that originated from fosmids and larger clones can be admixture mapped in this way, by using publicly available whole-genome sequence data. Genome assembly efforts and future builds of the human genome reference will be strongly informed by this localization of genes and other euchromatic sequences that are embedded within highly repetitive pericentromeric regions.     

A computational study of Shewanella oneidensis MR-1: structural prediction and functional inference of hypothetical proteins

The genomes of many organisms have been sequenced in the last 5 years. Typically about 30% of predicted genes from a newly sequenced genome cannot be given functional assignments using sequence comparison methods. In these situations three-dimensional structural predictions combined with a suite of computational tools can suggest possible functions for these hypothetical proteins. Suggesting functions may allow better interpretation of experimental data (e.g., microarray data and mass spectroscopy data) and help experimentalists design new experiments. In this paper, we focus on three hypothetical proteins of Shewanella oneidensis MR-1 that are potentially related to iron transport/metabolism based on microarray experiments. The threading program PROSPECT was used for protein structural predictions and functional annotation, in conjunction with literature search and other computational tools. Computational tools were used to perform transmembrane domain predictions, coiled coil predictions, signal peptide predictions, sub-cellular localization predictions, motif prediction, and operon structure evaluations. Combined computational results from all tools were used to predict roles for the hypothetical proteins. This method, which uses a suite of computational tools that are freely available to academic users, can be used to annotate hypothetical proteins in general.     

Improved analytical methods for microarray-based genome-composition analysis

Background  

Whereas genome sequencing has given us high-resolution pictures of many different species of bacteria, microarrays provide a means of obtaining information on genome composition for many strains of a given species. Genome-composition analysis using microarrays, or 'genomotyping', can be used to categorize genes into 'present' and 'divergent' categories based on the level of hybridization signal. This typically involves selecting a signal value that is used as a cutoff to discriminate present (high signal) and divergent (low signal) genes. Current methodology uses empirical determination of cutoffs for classification into these categories, but this methodology is subject to several problems that can result in the misclassification of many genes.     

An integrated web interface for large-scale characterization of sequence data

Large-scale genome projects require the analysis of large amounts of raw data. This analysis often involves the application of a chain of biology-based programs. Many of these programs are difficult to operate because they are non-integrated, command-line driven, and platform-dependent. The problem is compounded when the number of data files involved is large, making navigation and status-tracking difficult. To demonstrate how this problem can be addressed, we have created a platform-independent Web front end that integrates a set of programs used in a genomic project analyzing gene function by transposon mutagenesis in Saccharomyces cerevisiae. In particular, these programs help define a large number of transposon insertion events within the yeast genome, identifying both the precise site of transposon insertion as well as potential open reading frames disrupted by this insertion event. Our Web interface facilitates this analysis by performing the following tasks. Firstly, it allows each of the analysis programs to be launched against multiple directories of data files. Secondly, it allows the user to view, download, and upload files generated by the programs. Thirdly, it indicates which sets of data directories have been processed by each program. Although designed specifically to aid in this project, our interface exemplifies a general approach by which independent software programs may be integrated into an efficient protocol for large-scale genomic data processing. Electronic Publication     

Preparing for the flood: evolutionary biology in the age of genomics

Genome sequencing projects will soon provide a large number of gene sequences from a wide variety of species. Although the impetus of these genome projects involves biotechnology rather than evolutionary biology, several substantial evolutionary problems may potentially be addressed using genomic data.     

Protein-protein interaction through beta-strand addition

Protein-protein interactions have essential roles at almost every level of organization and communication in living cells. During complex formation, proteins can interact via covalent, surface-surface or peptide-surface contacts. Many protein complexes are now known to involve the binding of linear motifs in one of the binding partners. An emerging mechanism of such non-covalent peptide-surface interaction involves the donation or addition of a beta strand in the ligand to a beta sheet or a beta strand in the receptor. Such 'beta-strand addition' contacts can dictate or modulate binding specificity and affinity, or can be used in more promiscuous protein-protein contacts. Three main classes of beta-strand addition can be distinguished: beta-sheet augmentation; beta-strand insertion and fold complementation; and beta-strand zippering. A survey of protein-protein complexes in the protein data bank identifies beta-strand additions in many important metabolic pathways. Targeting these interactions might, thus, provide novel routes for rational drug design.     

Mouse models of male infertility

Spermatogenesis is a complex process that involves stem-cell renewal, genome reorganization and genome repackaging, and that culminates in the production of motile gametes. Problems at all stages of spermatogenesis contribute to human infertility, but few of them can be modelled in vitro or in cell culture. Targeted mutagenesis in the mouse provides a powerful method to analyse these steps and has provided new insights into the origins of male infertility.     

Three aldo–keto reductases of the yeast Saccharomyces cerevisiae

Saccharomyces cerevisiae is an industrially important yeast, which is also used extensively as a model eukaryote. The S. cerevisiae genome has been sequenced in its entirety and therefore represents an ideal organism in which to carry out functional analysis of genes. We have identified several open reading frames in the S. cerevisiae genome which show significant similarity to members of the aldo–keto reductase superfamily. The physiological roles of these gene products have not been previously determined, but their similarity to other enzymes suggests they may perform roles in carbohydrate metabolism and detoxification pathways. Cloning and expression of three of these enzymes has allowed their substrate specificities to be determined. Expression profiling and gene disruption analysis will allow potential roles for these enzymes within the cell to be examined.     

Chromatin remodelling in mammalian differentiation: lessons from ATP-dependent remodellers

The initiation of cellular differentiation involves alterations in gene expression that depend on chromatin changes, at the level of both higher-order structures and individual genes. Consistent with this, chromatin-remodelling enzymes have key roles in differentiation and development. The functions of ATP-dependent chromatin-remodelling enzymes have been studied in several mammalian differentiation pathways, revealing cell-type-specific and gene-specific roles for these proteins that add another layer of precision to the regulation of differentiation. Recent studies have also revealed a role for ATP-dependent remodelling in regulating the balance between proliferation and differentiation, and have uncovered intriguing links between chromatin remodelling and other cellular processes during differentiation, including recombination, genome organization and the cell cycle.     

Determination of protease cleavage site motifs using mixture-based oriented peptide libraries.

The number of known proteases is increasing at a tremendous rate as a consequence of genome sequencing projects. Although one can guess at the functions of these novel enzymes by considering sequence homology to known proteases, there is a need for new tools to rapidly provide functional information on large numbers of proteins. We describe a method for determining the cleavage site specificity of proteolytic enzymes that involves pooled sequencing of peptide library mixtures. The method was used to determine cleavage site motifs for six enzymes in the matrix metalloprotease (MMP) family. The results were validated by comparison with previous literature and by analyzing the cleavage of individually synthesized peptide substrates. The library data led us to identify the proteoglycan neurocan as a novel MMP-2 substrate. Our results indicate that a small set of libraries can be used to quickly profile an expanding protease family, providing information applicable to the design of inhibitors and to the identification of protein substrates.     

Ensuring the fidelity of recombination in mammalian chromosomes

Mammalian cells frequently depend on homologous recombination (HR) to repair DNA damage accurately and to help rescue stalled or collapsed replication forks. The essence of HR is an exchange of nucleotides between identical or nearly identical sequences. Although HR fulfills important biological roles, recombination between inappropriate sequence partners can lead to translocations or other deleterious rearrangements and such events must be avoided. For example, the recombination machinery must follow stringent rules to preclude recombination between the many repetitive elements in a mammalian genome that share significant but imperfect homology. This paper takes a conceptual approach in addressing the homology requirements for recombination in mammalian genomes as well as the general strategy used by cells to reject recombination between similar but imperfectly matched sequences. A mechanism of heteroduplex rejection that involves the unwinding of recombination intermediates that may form between mismatched sequences is discussed.