Joint evolution of dispersal and inbreeding load.

Inbreeding avoidance is often invoked to explain observed patterns of dispersal, and theoretical models indeed point to a possibly important role. However, while inbreeding load is usually assumed constant in these models, it is actually bound to vary dynamically under the combined influences of mutation, drift, and selection and thus to evolve jointly with dispersal. Here we report the results of individual-based stochastic simulations allowing such a joint evolution. We show that strongly deleterious mutations should play no significant role, owing to the low genomic mutation rate for such mutations. Mildly deleterious mutations, by contrast, may create enough heterosis to affect the evolution of dispersal as an inbreeding-avoidance mechanism, but only provided that they are also strongly recessive. If slightly recessive, they will spread among demes and accumulate at the metapopulation level, thus contributing to mutational load, but not to heterosis. The resulting loss of viability may then combine with demographic stochasticity to promote population fluctuations, which foster indirect incentives for dispersal. Our simulations suggest that, under biologically realistic parameter values, deleterious mutations have a limited impact on the evolution of dispersal, which on average exceeds by only one-third the values expected from kin-competition avoidance.

Numerical simulation of gel electrophoresis of DNA knots in weak and strong electric fields.

Gel electrophoresis allows one to separate knotted DNA (nicked circular) of equal length according to the knot type. At low electric fields, complex knots, being more compact, drift faster than simpler knots. Recent experiments have shown that the drift velocity dependence on the knot type is inverted when changing from low to high electric fields. We present a computer simulation on a lattice of a closed, knotted, charged DNA chain drifting in an external electric field in a topologically restricted medium. Using a Monte Carlo algorithm, the dependence of the electrophoretic migration of the DNA molecules on the knot type and on the electric field intensity is investigated. The results are in qualitative and quantitative agreement with electrophoretic experiments done under conditions of low and high electric fields.

One billion years of bZIP transcription factor evolution: conservation and change in dimerization and DNA-binding site specificity.

The genomic era has revealed that the large repertoire of observed animal phenotypes is dependent on changes in the expression patterns of a finite number of genes, which are mediated by a plethora of transcription factors (TFs) with distinct specificities. The dimerization of TFs can also increase the complexity of a genetic regulatory network manifold, by combining a small number of monomers into dimers with distinct functions. Therefore, studying the evolution of these dimerizing TFs is vital for understanding how complexity increased during animal evolution. We focus on the second largest family of dimerizing TFs, the basic-region leucine zipper (bZIP), and infer when it expanded and how bZIP DNA-binding and dimerization functions evolved during the major phases of animal evolution. Specifically, we classify the metazoan bZIPs into 19 families and confirm the ancient nature of at least 13 of these families, predating the split of the cnidaria. We observe fixation of a core dimerization network in the last common ancestor of protostomes-deuterostomes. This was followed by an expansion of the number of proteins in the network, but no major dimerization changes in interaction partners, during the emergence of vertebrates. In conclusion, the bZIPs are an excellent model with which to understand how DNA binding and protein interactions of TFs evolved during animal evolution.

Computational problems in perfect phylogeny haplotyping: typing without calling the allele.

A haplotype is an m-long binary vector. The XOR-genotype of two haplotypes is the m-vector of their coordinate-wise XOR. We study the following problem: Given a set of XOR-genotypes, reconstruct their haplotypes so that the set of resulting haplotypes can be mapped onto a perfect phylogeny (PP) tree. The question is motivated by studying population evolution in human genetics, and is a variant of the perfect phylogeny haplotyping problem that has received intensive attention recently. Unlike the latter problem, in which the input is "full" genotypes, here we assume less informative input, and so may be more economical to obtain experimentally. Building on ideas of Gusfield, we show how to solve the problem in polynomial time, by a reduction to the graph realization problem. The actual haplotypes are not uniquely determined by that tree they map onto, and the tree itself may or may not be unique. We show that tree uniqueness implies uniquely determined haplotypes, up to inherent degrees of freedom, and give a sufficient condition for the uniqueness. To actually determine the haplotypes given the tree, additional information is necessary. We show that two or three full genotypes suffice to reconstruct all the haplotypes, and present a linear algorithm for identifying those genotypes.

Binding of double-strand breaks in DNA by human Rad52 protein.

Double-strand breaks (DSBs) in DNA are caused by ionizing radiation. These chromosomal breaks can kill the cell unless repaired efficiently, and inefficient or inappropriate repair can lead to mutation, gene translocation and cancer. Two proteins that participate in the repair of DSBs are Rad52 and Ku: in lower eukaryotes such as yeast, DSBs are repaired by Rad52-dependent homologous recombination, whereas vertebrates repair DSBs primarily by Ku-dependent non-homologous end-joining. The contribution of homologous recombination to vertebrate DSB repair, however, is important. Biochemical studies indicate that Ku binds to DNA ends and facilitates end-joining. Here we show that human Rad52, like Ku, binds directly to DSBs, protects them from exonuclease attack and facilitates end-to-end interactions. A model for repair is proposed in which either Ku or Rad52 binds the DSB. Ku directs DSBs into the non-homologous end-joining repair pathway, whereas Rad52 initiates repair by homologous recombination. Ku and Rad52, therefore, direct entry into alternative pathways for the repair of DNA breaks.

Quantitative antibiogram as a typing method for the prospective epidemiological surveillance and control of MRSA: comparison with molecular typing.

OBJECTIVE: Evaluation of the quantitative antibiogram as an epidemiological tool for the prospective typing of methicillin-resistant Staphylococcus aureus (MRSA), and comparison with ribotyping. METHODS: The method is based on the multivariate analysis of inhibition zone diameters of antibiotics in disk diffusion tests. Five antibiotics were used (erythromycin, clindamycin, cotrimoxazole, gentamicin, and ciprofloxacin). Ribotyping was performed using seven restriction enzymes (EcoRV, HindIII, KpnI, PstI, EcoRI, SfuI, and BamHI). SETTING: 1,000-bed tertiary university medical center. RESULTS: During a 1-year period, 31 patients were found to be infected or colonized with MRSA. Cluster analysis of antibiogram data showed nine distinct antibiotypes. Four antibiotypes were isolated from multiple patients (2, 4, 7, and 13, respectively). Five additional antibiotypes were isolated from the remaining five patients. When analyzed with respect to the epidemiological data, the method was found to be equivalent to ribotyping. Among 206 staff members who were screened, six were carriers of MRSA. Both typing methods identified concordant of MRSA types in staff members and in the patients under their care. CONCLUSIONS: The quantitative antibiogram was found to be equivalent to ribotyping as an epidemiological tool for typing of MRSA in our setting. Thus, this simple, rapid, and readily available method appears to be suitable for the prospective surveillance and control of MRSA for hospitals that do not have molecular typing facilities and in which MRSA isolates are not uniformly resistant or susceptible to the antibiotics tested.

Similarities and differences of polyadenylation signals in human and fly.

BACKGROUND: Cleavage of messenger RNA (mRNA) precursors is an essential step in mRNA maturation. The signal recognized by the cleavage enzyme complex has been characterized as an A rich region upstream of the cleavage site containing a motif with consensus AAUAAA, followed by a U or UG rich region downstream of the cleavage site. RESULTS: We studied these signals using exhaustive databases of cleavage sites obtained from aligning raw expressed sequence tags (EST) sequences to genomic sequences in Homo sapiens and Drosophila melanogaster. These data show that the polyadenylation signal is highly conserved in human and fly. In addition, de novo motif searches generated a refined description of the U-rich downstream sequence (DSE) element, which shows more divergence between the two species. These refined motifs are applied, within a Hidden Markov Model (HMM) framework, to predict mRNA cleavage sites. CONCLUSION: We demonstrate that the DSE is a specific motif in both human and Drosophila. These findings shed light on the sequence correlates of a highly conserved biological process, and improve in silico prediction of 3' mRNA cleavage and polyadenylation sites.

The expansion of amino-acid repeats is not associated to adaptive evolution in mammalian genes.

BACKGROUND: The expansion of amino acid repeats is determined by a high mutation rate and can be increased or limited by selection. It has been suggested that recent expansions could be associated with the potential of adaptation to new environments. In this work, we quantify the strength of this association, as well as the contribution of potential confounding factors. RESULTS: Mammalian positively selected genes have accumulated more recent amino acid repeats than other mammalian genes. However, we found little support for an accelerated evolutionary rate as the main driver for the expansion of amino acid repeats. The most significant predictors of amino acid repeats are gene function and GC content. There is no correlation with expression level. CONCLUSIONS: Our analyses show that amino acid repeat expansions are causally independent from protein adaptive evolution in mammalian genomes. Relaxed purifying selection or positive selection do not associate with more or more recent amino acid repeats. Their occurrence is slightly favoured by the sequence context but mainly determined by the molecular function of the gene.

A multiplicity of factors contributes to selective RNA polymerase III occupancy of a subset of RNA polymerase III genes in mouse liver.

The genomic loci occupied by RNA polymerase (RNAP) III have been characterized in human culture cells by genome-wide chromatin immunoprecipitations, followed by deep sequencing (ChIP-seq). These studies have shown that only ∼40% of the annotated 622 human tRNA genes and pseudogenes are occupied by RNAP-III, and that these genes are often in open chromatin regions rich in active RNAP-II transcription units. We have used ChIP-seq to characterize RNAP-III-occupied loci in a differentiated tissue, the mouse liver. Our studies define the mouse liver RNAP-III-occupied loci including a conserved mammalian interspersed repeat (MIR) as a potential regulator of an RNAP-III subunit-encoding gene. They reveal that synteny relationships can be established between a number of human and mouse RNAP-III genes, and that the expression levels of these genes are significantly linked. They establish that variations within the A and B promoter boxes, as well as the strength of the terminator sequence, can strongly affect RNAP-III occupancy of tRNA genes. They reveal correlations with various genomic features that explain the observed variation of 81% of tRNA scores. In mouse liver, loci represented in the NCBI37/mm9 genome assembly that are clearly occupied by RNAP-III comprise 50 Rn5s (5S RNA) genes, 14 known non-tRNA RNAP-III genes, nine Rn4.5s (4.5S RNA) genes, and 29 SINEs. Moreover, out of the 433 annotated tRNA genes, half are occupied by RNAP-III. Transfer RNA gene expression levels reflect both an underlying genomic organization conserved in dividing human culture cells and resting mouse liver cells, and the particular promoter and terminator strengths of individual genes.

Going the distance: human population genetics in a clinal world.

Global human genetic variation is greatly influenced by geography, with genetic differentiation between populations increasing with geographic distance and within-population diversity decreasing with distance from Africa. In fact, these 'clines' can explain most of the variation in human populations. Despite this, population genetics inferences often rely on models that do not take geography into account, which could result in misleading conclusions when working at global geographic scales. Geographically explicit approaches have great potential for the study of human population genetics. Here, we discuss the most promising avenues of research in the context of human settlement history and the detection of genomic elements under natural selection. We also review recent technical advances and address the challenges of integrating geography and genetics.

How demography, life history, and kinship shape the evolution of genomic imprinting.

How phenomena like helping, dispersal, or the sex ratio evolve depends critically on demographic and life-history factors. One phenotype that is of particular interest to biologists is genomic imprinting, which results in parent-of-origin-specific gene expression and thus deviates from the predictions of Mendel's rules. The most prominent explanation for the evolution of genomic imprinting, the kinship theory, originally specified that multiple paternity can cause the evolution of imprinting when offspring affect maternal resource provisioning. Most models of the kinship theory do not detail how population subdivision, demography, and life history affect the evolution of imprinting. In this work, we embed the classic kinship theory within an island model of population structure and allow for diverse demographic and life-history features to affect the direction of selection on imprinting. We find that population structure does not change how multiple paternity affects the evolution of imprinting under the classic kinship theory. However, if the degree of multiple paternity is not too large, we find that sex-specific migration and survival and generation overlap are the primary factors determining which allele is silenced. This indicates that imprinting can evolve purely as a result of sex-related asymmetries in the demographic structure or life history of a species.

Sustained transgene expression using MAR elements.

Matrix attachment regions (MARs) are DNA sequences that may be involved in anchoring DNA/chromatin to the nuclear matrix and they have been described in both mammalian and plant species. MARs possess a number of features that facilitate the opening and maintenance of euchromatin. When incorporated into viral or non-viral vectors MARs can increase transgene expression and limit position-effects. They have been used extensively to improve transgene expression and recombinant protein production and promising studies on the potential use of MAR elements for mammalian gene therapy have appeared. These illustrate how MARs may be used to mediate sustained or higher levels of expression of therapeutic genes and/or to reduce the viral vector multiplicity of infection required to achieve consistent expression. More recently, the discovery of potent MAR elements and the development of improved vectors for transgene delivery, notably non-viral episomal vectors, has strengthened interest in their use to mediate expression of therapeutic transgenes. This article will describe the progress made in this field, and it will discuss future directions and issues to be addressed.

The fasting-induced adipose factor/angiopoietin-like protein 4 is physically associated with lipoproteins and governs plasma lipid levels and adiposity.

Proteins secreted from adipose tissue are increasingly recognized to play an important role in the regulation of glucose metabolism. However, much less is known about their effect on lipid metabolism. The fasting-induced adipose factor (FIAF/angiopoietin-like protein 4/peroxisome proliferator-activated receptor gamma angiopoietin-related protein) was previously identified as a target of hypolipidemic fibrate drugs and insulin-sensitizing thiazolidinediones. Using transgenic mice that mildly overexpress FIAF in peripheral tissues we show that FIAF is an extremely powerful regulator of lipid metabolism and adiposity. FIAF overexpression caused a 50% reduction in adipose tissue weight, partly by stimulating fatty acid oxidation and uncoupling in fat. In addition, FIAF overexpression increased plasma levels of triglycerides, free fatty acids, glycerol, total cholesterol, and high density lipoprotein (HDL)-cholesterol. Functional tests indicated that FIAF overexpression severely impaired plasma triglyceride clearance but had no effect on very low density lipoprotein production. The effects of FIAF overexpression were amplified by a high fat diet, resulting in markedly elevated plasma and liver triglycerides, plasma free fatty acids, and plasma glycerol levels, and impaired glucose tolerance in FIAF transgenic mice fed a high fat diet. Remarkably, in mice the full-length form of FIAF was physically associated with HDL, whereas truncated FIAF was associated with low density lipoprotein. In human both full-length and truncated FIAF were associated with HDL. The composite data suggest that via physical association with plasma lipoproteins, FIAF acts as a powerful signal from fat and other tissues to prevent fat storage and stimulate fat mobilization. Our data indicate that disturbances in FIAF signaling might be involved in dyslipidemia.

MYC and metastasis.

Aggressive primary tumors express transcriptional signatures that correlate with their metastatic propensity. A number of these signatures have been deployed in the clinic as risk stratification tools. However, the molecular basis of these clinically useful prognostic signatures has remained a largely unresolved area of controversy. We recently found that many prognostic signatures reflect the activity of the MYC oncogene, which in turn regulates tumor metastasis through specific effects on cancer cell invasion and migration. These findings offer a general framework for understanding the molecular basis of clinically prognostic transcriptional signatures and suggest potentially new avenues for studying metastasis.