Evolutionary preservation of redundant duplicated genes

Gene duplication events produce both perfect and imperfect copies of genes. Perfect copies are said to be functionally redundant when knockout of one gene produces no 'scoreable', phenotypic effects. Preserving identical, duplicate copies of genes is problematic as all copies are prone to accumulate neutral mutations as pseudogenes, or more rarely, evolve into new genes with novel functions. We summarise theoretical treatments for the invasion and subsequent evolutionary modification of functionally redundant genes. We then consider the preservation of functionally identical copies of a gene over evolutionary time. We present several models for conserving redundancy: asymmetric mutation, asymmetric efficacy, pleiotropy, developmental buffering, allelic competition and regulatory asymmetries. In all cases, some form of symmetry breaking is required to maintain functional redundancy indefinitely.     

Accepted mutations in a gene family: Evolutionary diversification of duplicated DNA

Summary We report and compare the DNA sequences of 14 silkmoth (Antheraea polyphemus) chorion genes, derived from either cDNA or chromosomal DNA clones. Seven of these genes are members of the A multigene family, and seven are members of the B family. Where available, the previously reported (Jones and Kafatos 1980) intronic and extragenic flanking DNA sequences are also considered. Closely related sequences are compared, revealing the types of spontaneous mutations that were fixed during paralogous evolution. Segmental mutations (i.e. mutations other than substitutions) are nearly always interpretable as small duplications or deletions. related to small direct repeats. Segmental mutations are strongly constrained in the coding regions, although they do occur. Nucleotide substitutions also appear to be under selective constraints: relatively few substitutions leading to amino acid replacements are accepted, silent substitutions leading to some codons (especially purine-terminated ones) are disfavored, and different compositional biases are maintained in different parts of the sequences. Other sequence differences can be interpreted as indicative of neutral drift, including most differences in non-coding regions and most T/C transitions in third-base positions. In the non-coding regions, which are thought to be only loosely constrained by selection, transitions are observed more frequently than might be expected: they account for 52% of all substitutions, and they appear to be favored two to threefold over transversions when allowance is made for the skewed base composition of these regions.     

Evolutionary dynamics of duplicated genes in plants

Gene duplication, arising from region-specific duplication or genome-wide polyploidization, is a prominent feature in plant genome evolution. Understanding the mechanisms generating duplicate gene copies and the subsequent dynamics among gene duplicates is vital because these investigations shed light on regional and genome-wide aspects of evolutionary forces shaping intra- and interspecific genome contents, evolutionary relationships, and interactions. This review discusses recent gene duplication analyses in plants, focusing on the molecular and evolutionary dynamics occurring at three different timescales following duplication: (1). initial establishment and persistence of cytotypes, (2). interactions among duplicate gene copies, and (3). longer term differentiation between duplicated genes. These relative time points are presented in terms of their potential adaptive significance and impact on plant evolutionary genomics research.     

Evolutionary genomics of the recently duplicated amphioxus Hairy genes

Amphioxus Hairy genes have gone through a number of lineage-specific duplications, resulting in eight members, some of which are differentially expressed in the embryo. In order to gain insights into the evolution and function of this gene family we have compared their genomic structure and searched for conserved non-coding sequence elements. We have found that introns have been lost independently from these genes at least twice and after the duplication events. By carrying out phylogenetic footprinting between paralogues expressed in the embryo, we have found a differential distribution of conserved elements that could explain the limited overlap in expression patterns of Hairy genes in the amphioxus embryo. Furthermore, clustering of RBP-Jk binding sites in these conserved elements suggests that amphioxus Hairy genes are downstream targets of the Notch signaling pathway, as occurs in vertebrates. All of this evidence suggests that amphioxus Hairy genes have gone through a process of subfunctionalization shortly after their duplication, representing an extreme and rapid case of the duplication-degeneration-complementation model.     

Evolutionary dynamics of duplicated microsatellites shared by sex chromosomes

Segmental duplications on sex chromosomes constitute an important proportion of recent duplications (approximately 30%). Among those, the evolution of duplicated noncoding DNA is still poorly investigated. We focus our work on repeated DNA sequences extensively used in population genetics and evolution: microsatellites. Six duplicated (CA), microsatellite loci, located on the homologous region of human sex chromosomes, were studied at the intraspecific level in Homo sapiens and by an orthologous comparison in eight primate species. At the intraspecific level, we evaluated the congruence in paralogous divergence between the flanking sequences of the six microsatellites and the approximately 2.2-kb surrounding sequences and observed that both phylogenies are congruent. At the interspecific level (8 species of primates: 54 individuals), we analyzed the sequence polymorphism and divergence of each orthologous locus for both the flanking sequence and the microsatellite. The results showed a lower divergence of flanking sequences than expected in noncoding DNA and a relative stability of the first nucleotides close to the microsatellite. The location of each CAIII locus in a Low Copy Repeated element containing duplicated VCX/Y genes (approximately 1 kb) suggested that direct or indirect selection could explain these results. Moreover, the substitution rates in the flanking sequences and in the microsatellites were correlated. Thus, the evolutionary dynamics of microsatellites seems closely linked to the variation of spontaneous mutations in the surrounding regions.     

Functional analysis of a duplicated linked pair of ribosomal protein genes in Saccharomyces cerevisiae.

Ribosomal protein genes RP28 and S16A (RP55) are closely linked. Another set of this pair of genes exists in the genome (copy 2), genetically unlinked to copy 1. By using gene replacement techniques, we have shown that RP28 from copy 1 is required for vegetative growth and that the cells need S16A from copy 2 to achieve maximum growth rate.     

Evolutionary accessibility of mutational pathways

Functional effects of different mutations are known to combine to the total effect in highly nontrivial ways. For the trait under evolutionary selection ('fitness'), measured values over all possible combinations of a set of mutations yield a fitness landscape that determines which mutational states can be reached from a given initial genotype. Understanding the accessibility properties of fitness landscapes is conceptually important in answering questions about the predictability and repeatability of evolutionary adaptation. Here we theoretically investigate accessibility of the globally optimal state on a wide variety of model landscapes, including landscapes with tunable ruggedness as well as neutral 'holey' landscapes. We define a mutational pathway to be accessible if it contains the minimal number of mutations required to reach the target genotype, and if fitness increases in each mutational step. Under this definition accessibility is high, in the sense that at least one accessible pathway exists with a substantial probability that approaches unity as the dimensionality of the fitness landscape (set by the number of mutational loci) becomes large. At the same time the number of alternative accessible pathways grows without bounds. We test the model predictions against an empirical 8-locus fitness landscape obtained for the filamentous fungus Aspergillus niger. By analyzing subgraphs of the full landscape containing different subsets of mutations, we are able to probe the mutational distance scale in the empirical data. The predicted effect of high accessibility is supported by the empirical data and is very robust, which we argue reflects the generic topology of sequence spaces. Together with the restrictive assumptions that lie in our definition of accessibility, this implies that the globally optimal configuration should be accessible to genome wide evolution, but the repeatability of evolutionary trajectories is limited owing to the presence of a large number of alternative mutational pathways.     

Prediction of protein structure by simulating coarse-grained folding pathways: a preliminary report

A set of software tools designed to study protein structure and kinetics has been developed. The core of these tools is a program called Folding Machine (FM) which is able to generate low resolution folding pathways using modest computational resources. The FM is based on a coarse-grained kinetic ab initio Monte-Carlo sampler that can optionally use information extracted from secondary structure prediction servers or from fragment libraries of local structure. The model underpinning this algorithm contains two novel elements: (a) the conformational space is discretized using the Ramachandran basins defined in the local phi-psi energy maps; and (b) the solvent is treated implicitly by rescaling the pairwise terms of the non-bonded energy function according to the local solvent environments. The purpose of this hybrid ab initio/knowledge-based approach is threefold: to cover the long time scales of folding, to generate useful 3-dimensional models of protein structures, and to gain insight on the protein folding kinetics. Even though the algorithm is not yet fully developed, it has been used in a recent blind test of protein structure prediction (CASP5). The FM generated models within 6 A backbone rmsd for fragments of about 60-70 residues of alpha-helical proteins. For a CASP5 target that turned out to be natively unfolded, the trajectory obtained for this sequence uniquely failed to converge. Also, a new measure to evaluate structure predictions is presented and used along the standard CASP assessment methods. Finally, recent improvements in the prediction of beta-sheet structures are briefly described.     

Evolutionary genetics: The economics of mutation.

The presence of mutator genotypes in populations of bacteria may be favoured by selection because they produce rare beneficial mutations and thereby increase the rate of adaptive evolution. Recent work, however, shows that the relationship between mutation rates and adaptive evolution is more complicated.     

Evolutionary pathways to self-fertilization in a tristylous plant species

Evolutionary transitions from outcrossing to selfing occur commonly in heterostylous genera. The morphological polymorphisms that characterize heterostyly provide opportunities for different pathways for selfing to evolve. Here, we investigate the origins and pathways by which selfing has evolved in tristylous Eichhornia paniculata by providing new evidence based on morphology, DNA sequences and genetic analysis. The primary pathway from outcrossing to selfing involves the stochastic loss of the short-styled morph (S-morph) from trimorphic populations, followed by the spread of selfing variants of the mid-styled morph (M-morph). However, the discovery of selfing variants of the long-styled morph (L-morph) in Central America indicates a secondary pathway and distinct origin for selfing. Comparisons of multi-locus nucleotide sequences from 27 populations sampled from throughout the geographical range suggest multiple transitions to selfing. Genetic analysis of selfing variants of the L- and M-morphs demonstrates recessive control of the loss of herkogamy, although the number of factors appears to differ between the forms. Early stages in the establishment of selfing involve developmental instability in the formation of flowers capable of autonomous self-pollination. The relatively simple genetic control of herkogamy reduction and frequent colonizing episodes may often create demographic conditions favouring transitions to selfing in E. paniculata .     

Recovering filter-based microarray data for pathways analysis using a multipoint alignment strategy

The use of commercial microarrays is rapidly becoming the method of choice for profiling gene expression and assessing various disease states. Research Genetics has provided a series of biological and software tools to the research community for these analyses. The fidelity of data analysis using these tools is dependent on a series of well-defined reference control points in the array. During the course of our investigations, it became apparent that in some instances the reference control points that are required for analysis became lost in background noise. This effectively halted the analysis and the recovery of any information contained within that experiment. To recover this data and to increase analytical veracity, the simple strategy of superimposing a template of reference control points onto the experimental array was developed. The utility of this tool is established in this communication.     

Pysim-sv: a package for simulating structural variation data with GC-biases

Background

Structural variations (SVs) are wide-spread in human genomes and may have important implications in disease-related and evolutionary studies. High-throughput sequencing (HTS) has become a major platform for SV detection and simulation serves as a powerful and cost-effective approach for benchmarking SV detection algorithms. Accurate performance assessment by simulation requires the simulator capable of generating simulation data with all important features of real data, such GC biases in HTS data and various complexities in tumor data. However, no available package has systematically addressed all issues in data simulation for SV benchmarking.

Results

Pysim-sv is a package for simulating HTS data to evaluate performance of SV detection algorithms. Pysim-sv can introduce a wide spectrum of germline and somatic genomic variations. The package contains functionalities to simulate tumor data with aneuploidy and heterogeneous subclones, which is very useful in assessing algorithm performance in tumor studies. Furthermore, Pysim-sv can introduce GC-bias, the most important and prevalent bias in HTS data, in the simulated HTS data.

Conclusions

Pysim-sv provides an unbiased toolkit for evaluating HTS-based SV detection algorithms.

     

DNA conformational change in Gal repressor-operator complex: involvement of central G-C base pair(s) of dyad symmetry.

Gal repressor dimer binds to two gal operator sites, OE and OI, which are 16 bp long similar sequences with hyphenated dyad symmetries (11,12). Repressor occupation hinders the reactivity of the N7 atoms in the major groups of guanines, located at positions 1, 3 and 8, and the rotational 1', 3' and 8' of the symmetries. We have shown that Gal repressor binding to OE or OI DNA fragments increases the circular dichroism (CD) spectral peak in the 270 to 300 nm range. The CD change is similar to that observed for Lac repressor binding to its operator site (14). It is consistent with a DNA conformational change during complex formation between Gal repressor and OE and OI DNA. The CD spectral change was not observed when the central 8,8' G-C base pairs in the DNA-protein complex were replaced by A-T base pairs, whereas substitution of the 1,1' G-C base pairs do show the accompanying increase in the spectra during repressor binding. The absence of CD change of the Gal repressor complex with DNA mutated at the 8,8' base pairs suggest that the central G-C base pairs are required for the repressor induced conformational change.     

Mathematical analysis of endothelial sibling pair cell-cell interactions using time-lapse cinematography data

The sibling pairs from two different endothelial cell cultures were analysed by time-lapse cinematography. It was shown that wounded and regular (low density seeded) cultures differed in the behaviour patterns of their siblings. The cultures differed most significantly in the minimum interdivision time (IDT) which was 27% lower for the wounded culture. In the wounded culture there was a greater correlation of IDT values between sibling pairs. IDT values recorded both for paired and for unpaired cells were shorter for the wounded than for the regular culture. The mean IDT for unpaired cells was longer than the mean IDT for paired cells in the regular culture. Thus paired cells in the regular culture, had shorter IDTs. but not as short as in the wounded culture. It was significant that in the wounded culture the first generation of siblings were very close (less than 150 μm apart) at division. Overall the behaviour differences between the two cultures resulted in a higher rate of increase in cell numbers, and thus faster repair, of the wounded monolayer.     

Evolutionary pathways of an ancient gene recX

RecX is a regulator of RecA activity by interacting with RecA protein or RecA filaments. Genes encoding RecX were found in genomes of a wide diversity of bacteria and some plants (e.g., Arabidopsis thaliana and Oryza sativa). Our comparative genome analysis showed that although members of the RecX family are found in many bacterial species, they are not found in archaea and the only gene found in eukaryotes is likely derived from bacteria genomes. It is therefore proposed that RecX is of bacterial origin, and the gene had presented in the common ancestor of bacteria. Moreover, bacterial RecX and plant RecX domain are homologues, and RecX domain in plants may have derived from bacteria via unknown pathways. Plant RecX-like protein was formed by a gene fusion event between a unique N-terminal domain of unknown origin and RecX domain within plant cells. Finally, three possible evolutionary pathways from bacteria to plant were discussed.