Mammalian genes are transcribed with widely different bursting kinetics.

In prokaryotes and eukaryotes, most genes appear to be transcribed during short periods called transcriptional bursts, interspersed by silent intervals. We describe how such bursts generate gene-specific temporal patterns of messenger RNA (mRNA) synthesis in mammalian cells. To monitor transcription at high temporal resolution, we established various gene trap cell lines and transgenic cell lines expressing a short-lived luciferase protein from an unstable mRNA, and recorded bioluminescence in real time in single cells. Mathematical modeling identified gene-specific on- and off-switching rates in transcriptional activity and mean numbers of mRNAs produced during the bursts. Transcriptional kinetics were markedly altered by cis-regulatory DNA elements. Our analysis demonstrated that bursting kinetics are highly gene-specific, reflecting refractory periods during which genes stay inactive for a certain time before switching on again.

An autoregulatory loop controlling CYP1A1 gene expression: role of H(2)O(2) and NFI.

Cytochrome P450 1A1 (CYP1A1), like many monooxygenases, can produce reactive oxygen species during its catalytic cycle. Apart from the well-characterized xenobiotic-elicited induction, the regulatory mechanisms involved in the control of the steady-state activity of CYP1A1 have not been elucidated. We show here that reactive oxygen species generated from the activity of CYP1A1 limit the levels of induced CYP1A1 mRNAs. The mechanism involves the repression of the CYP1A1 gene promoter activity in a negative-feedback autoregulatory loop. Indeed, increasing the CYP1A1 activity by transfecting CYP1A1 expression vectors into hepatoma cells elicited an oxidative stress and led to the repression of a reporter gene driven by the CYP1A1 gene promoter. This negative autoregulation is abolished by ellipticine (an inhibitor of CYP1A1) and by catalase (which catalyzes H(2)O(2) catabolism), thus implying that H(2)O(2) is an intermediate. Down-regulation is also abolished by the mutation of the proximal nuclear factor I (NFI) site in the promoter. The transactivating domain of NFI/CTF was found to act in synergy with the arylhydrocarbon receptor pathway during the induction of CYP1A1 by 2,3,7,8-tetrachloro-p-dibenzodioxin. Using an NFI/CTF-Gal4 fusion, we show that NFI/CTF transactivating function is decreased by a high activity of CYP1A1. This regulation is also abolished by catalase or ellipticine. Consistently, the transactivating function of NFI/CTF is repressed in cells treated with H(2)O(2), a novel finding indicating that the transactivating domain of a transcription factor can be targeted by oxidative stress. In conclusion, an autoregulatory loop leads to the fine tuning of the CYP1A1 gene expression through the down-regulation of NFI activity by CYP1A1-based H(2)O(2) production. This mechanism allows a limitation of the potentially toxic CYP1A1 activity within the cell.

Perturbation expansions of multilocus fixation probabilities for frequency-dependent selection with applications to the Hill-Robertson effect and to the joint evolution of helping and punishment.

Natural populations are of finite size and organisms carry multilocus genotypes. There are, nevertheless, few results on multilocus models when both random genetic drift and natural selection affect the evolutionary dynamics. In this paper we describe a formalism to calculate systematic perturbation expansions of moments of allelic states around neutrality in populations of constant size. This allows us to evaluate multilocus fixation probabilities (long-term limits of the moments) under arbitrary strength of selection and gene action. We show that such fixation probabilities can be expressed in terms of selection coefficients weighted by mean first passages times of ancestral gene lineages within a single ancestor. These passage times extend the coalescence times that weight selection coefficients in one-locus perturbation formulas for fixation probabilities. We then apply these results to investigate the Hill-Robertson effect and the coevolution of helping and punishment. Finally, we discuss limitations and strengths of the perturbation approach. In particular, it provides accurate approximations for fixation probabilities for weak selection regimes only (Ns < or = 1), but it provides generally good prediction for the direction of selection under frequency-dependent selection.

Space-time relatedness and Hamilton's rule for long-lasting behaviors in viscous populations.

Genes affect not only the behavior and fitness of their carriers but also that of other individuals. According to Hamilton's rule, whether a mutant gene will spread in the gene pool depends on the effects of its carrier on the fitness of all individuals in the population, each weighted by its relatedness to the carrier. However, social behaviors may affect not only recipients living in the generation of the actor but also individuals living in subsequent generations. In this note, I evaluate space-time relatedness coefficients for localized dispersal. These relatedness coefficients weight the selection pressures on long-lasting behaviors, which stem from a multigenerational gap between phenotypic expression by actors and the resulting environmental feedback on the fitness of recipients. Explicit values of space-time relatedness coefficients reveal that they can be surprisingly large for typical dispersal rates, even for hundreds of generations in the future.

Ever-young sex chromosomes in European tree frogs.

Non-recombining sex chromosomes are expected to undergo evolutionary decay, ending up genetically degenerated, as has happened in birds and mammals. Why are then sex chromosomes so often homomorphic in cold-blooded vertebrates? One possible explanation is a high rate of turnover events, replacing master sex-determining genes by new ones on other chromosomes. An alternative is that X-Y similarity is maintained by occasional recombination events, occurring in sex-reversed XY females. Based on mitochondrial and nuclear gene sequences, we estimated the divergence times between European tree frogs (Hyla arborea, H. intermedia, and H. molleri) to the upper Miocene, about 5.4-7.1 million years ago. Sibship analyses of microsatellite polymorphisms revealed that all three species have the same pair of sex chromosomes, with complete absence of X-Y recombination in males. Despite this, sequences of sex-linked loci show no divergence between the X and Y chromosomes. In the phylogeny, the X and Y alleles cluster according to species, not in groups of gametologs. We conclude that sex-chromosome homomorphy in these tree frogs does not result from a recent turnover but is maintained over evolutionary timescales by occasional X-Y recombination. Seemingly young sex chromosomes may thus carry old-established sex-determining genes, a result at odds with the view that sex chromosomes necessarily decay until they are replaced. This raises intriguing perspectives regarding the evolutionary dynamics of sexually antagonistic genes and the mechanisms that control X-Y recombination.

Chromatin domain boundaries delimited by a histone-binding protein in yeast.

When located next to chromosomal elements such as telomeres, genes can be subjected to epigenetic silencing. In yeast, this is mediated by the propagation of the SIR proteins from telomeres toward more centromeric regions. Particular transcription factors can protect downstream genes from silencing when tethered between the gene and the telomere, and they may thus act as chromatin domain boundaries. Here we have studied one such transcription factor, CTF-1, that binds directly histone H3. A deletion mutagenesis localized the barrier activity to the CTF-1 histone-binding domain. A saturating point mutagenesis of this domain identified several amino acid substitutions that similarly inhibited the boundary and histone binding activities. Chromatin immunoprecipitation experiments indicated that the barrier protein efficiently prevents the spreading of SIR proteins, and that it separates domains of hypoacetylated and hyperacetylated histones. Together, these results suggest a mechanism by which proteins such as CTF-1 may interact directly with histone H3 to prevent the propagation of a silent chromatin structure, thereby defining boundaries of permissive and silent chromatin domains.

Molecular phylogenetics of shrews (Mammalia: Soricidae) reveal timing of transcontinental colonizations.

We sequenced 2167 base pairs (bp) of mitochondrial DNA cytochrome b and 16S, and 1390 bp of nuclear genes BRCA1 and ApoB in shrews taxa (Eulipotyphla, family Soricidae). The aim was to study the relationships at higher taxonomic levels within this family, and in particular the position of difficult clades such as Anourosorex and Myosorex. The data confirmed two monophyletic subfamilies, Soricinae and Crocidurinae. In the former, the tribes Anourosoricini, Blarinini, Nectogalini, Notiosoricini, and Soricini were supported. The latter was formed by the tribes Myosoricini and Crocidurini. The genus Suncus appeared to be paraphyletic and included Sylvisorex. We further suggest a biogeographical hypothesis, which shows that North America was colonized by three independent lineages of Soricinae during middle Miocene. Our hypothesis is congruent with the first fossil records for these taxa. Using molecular dating, the first exchanges between Africa and Eurasia occurred during the middle Miocene. The last one took place in the Late Miocene, with the dispersion of the genus Crocidura through the old world.

Helicase-defective RuvB(D113E) promotes RuvAB-mediated branch migration in vitro.

In Escherichia coli, the RuvA and RuvB proteins interact at Holliday junctions to promote branch migration leading to the formation of heteroduplex DNA. RuvA provides junction-binding specificity and RuvB drives ATP-dependent branch migration. Since RuvB contains sequence motifs characteristic of a DNA helicase and RuvAB exhibit helicase activity in vitro, we have analysed the role of DNA unwinding in relation to branch migration. A mutant RuvB protein, RuvB(D113E), mutated in helicase motif II (the DExx box), has been purified to homogeneity. The mutant protein forms hexameric rings on DNA similar to those formed by wild-type protein and promotes branch migration in the presence of RuvA. However, RuvB(D113E) exhibits reduced ATPase activity and is severely compromised in its DNA helicase activity. Models for RuvAB-mediated branch migration that invoke only limited DNA unwinding activity are proposed.

Adaptive divergence of ancient gene duplicates in the avian MHC class II beta.

Gene duplication and neofunctionalization are known to be important processes in the evolution of phenotypic complexity. They account for important evolutionary novelties that confer ecological adaptation, such as the major histocompatibility complex (MHC), a multigene family crucial to the vertebrate immune system. In birds, two MHC class II β (MHCIIβ) exon 3 lineages have been recently characterized, and two hypotheses for the evolutionary history of MHCIIβ lineages were proposed. These lineages could have arisen either by 1) an ancient duplication and subsequent divergence of one paralog or by 2) recent parallel duplications followed by functional convergence. Here, we compiled a data set consisting of 63 MHCIIβ exon 3 sequences from six avian orders to distinguish between these hypotheses and to understand the role of selection in the divergent evolution of the two avian MHCIIβ lineages. Based on phylogenetic reconstructions and simulations, we show that a unique duplication event preceding the major avian radiations gave rise to two ancestral MHCIIβ lineages that were each likely lost once later during avian evolution. Maximum likelihood estimation shows that following the ancestral duplication, positive selection drove a radical shift from basic to acidic amino acid composition of a protein domain facing the α-chain in the MHCII α β-heterodimer. Structural analyses of the MHCII α β-heterodimer highlight that three of these residues are potentially involved in direct interactions with the α-chain, suggesting that the shift following duplication may have been accompanied by coevolution of the interacting α- and β-chains. These results provide new insights into the long-term evolutionary relationships among avian MHC genes and open interesting perspectives for comparative and population genomic studies of avian MHC evolution.

Genometric analyses of the organization of circular chromosomes: a universal pressure determines the direction of ribosomal RNA genes transcription relative to chromosome replication.

Selective pressures related to gene function and chromosomal architecture are acting on genome sequences and can be revealed, for instance, by appropriate genometric methods. Cumulative nucleotide skew analyses, i.e., GC, TA, and ORF orientation skews, predict the location of the origin of DNA replication for 88 out of 100 completely sequenced bacterial chromosomes. These methods appear fully reliable for proteobacteria, Gram-positives, and spirochetes as well as for euryarchaeotes. Based on this genome architecture information, coorientation analyses reveal that in prokaryotes, ribosomal RNA (rRNA) genes encoding the small and large ribosomal subunits are all transcribed in the same direction as DNA replication; that is, they are located along the leading strand. This result offers a simple and reliable method for circumscribing the region containing the origin of the DNA replication and reveals a strong selective pressure acting on the orientation of rRNA genes similar to the weaker one acting on the orientation of ORFs. Rate of coorientation of transfer RNA (tRNA) genes with DNA replication appears to be taxon-specific. Analyzing nucleotide biases such as GC and TA skews of genes and plotting one against the other reveals a taxonomic clusterization of species. All ribosomal RNA genes are enriched in Gs and depleted in Cs, the only so far known exception being the rRNA genes of deuterostomian mitochondria. However, this exception can be explained by the fact that in the chromosome of the human mitochondrion, the model of the deuterostomian organelle genome, DNA replication, and rRNA transcription proceed in opposite directions. A general rule is deduced from prokaryotic and mitochondrial genomes: ribosomal RNA genes that are transcribed in the same direction as the DNA replication are enriched in Gs, and those transcribed in the opposite direction are depleted in Gs.

The role of transcription factories-mediated interchromosomal contacts in the organization of nuclear architecture.

Using numerical simulations, we investigate the underlying physical effects responsible for the overall organization of chromosomal territories in interphase nuclei. In particular, we address the following three questions: (i) why are chromosomal territories with relatively high transcriptional activity on average, closer to the centre of cell's nucleus than those with the lower activity? (ii) Why are actively transcribed genes usually located at the periphery of their chromosomal territories? (iii) Why are pair-wise contacts between active and inactive genes less frequent than those involving only active or only inactive genes? We show that transcription factories-mediated contacts between active genes belonging to different chromosomal territories are instrumental for all these features of nuclear organization to emerge spontaneously due to entropic effects arising when chromatin fibres are highly crowded.

Splicing and the evolution of proteins in mammals.

It is often supposed that a protein's rate of evolution and its amino acid content are determined by the function and anatomy of the protein. Here we examine an alternative possibility, namely that the requirement to specify in the unprocessed RNA, in the vicinity of intron-exon boundaries, information necessary for removal of introns (e.g., exonic splice enhancers) affects both amino acid usage and rates of protein evolution. We find that the majority of amino acids show skewed usage near intron-exon boundaries, and that differences in the trends for the 2-fold and 4-fold blocks of both arginine and leucine show this to be owing to effects mediated at the nucleotide level. More specifically, there is a robust relationship between the extent to which an amino acid is preferred/avoided near boundaries and its enrichment/paucity in splice enhancers. As might then be expected, the rate of evolution is lowest near intron-exon boundaries, at least in part owing to splice enhancers, such that domains flanking intron-exon junctions evolve on average at under half the rate of exon centres from the same gene. In contrast, the rate of evolution of intronless retrogenes is highest near the domains where intron-exon junctions previously resided. The proportion of sequence near intron-exon boundaries is one of the stronger predictors of a protein's rate of evolution in mammals yet described. We conclude that after intron insertion selection favours modification of amino acid content near intron-exon junctions, so as to enable efficient intron removal, these changes then being subject to strong purifying selection even if nonoptimal for protein function. Thus there exists a strong force operating on protein evolution in mammals that is not explained directly in terms of the biology of the protein.

An ontology to clarify homology-related concepts.

Although homology is a fundamental concept in biology and is one of the shared channels of communication universal to all biology, it is difficult to find a consensus definition. Indeed, the interpretations of homology have changed as biology has progressed. New terms, such as paramorphism, have been introduced into the literature with mixed success. In addition, different research fields operate with different definitions of homology, for example the mechanistic usage of evo-devo is not strictly historical and would not be acceptable in cladistics. This makes a global understanding of homology complex, whereas the integration of evolutionary concepts into bioinformatics and genomics is increasingly important. We propose an ontology organizing homology and related concepts and hope this solution will also facilitate the integration and sharing of knowledge among the community.

Heterozygote advantage and the maintenance of polymorphism for multilocus traits.

Explaining how polymorphism is maintained in the face of selection remains a puzzle since selection tends to erode genetic variation. Provided an infinitely large unsubdivided population and no frequency-dependance of selective values, heterozygote advantage is the text book explanation for the maintenance of polymorphism when selection acts at a diallelic locus. Here, we investigate whether this remains true when selection acts at multiple diallelic loci. We use five different definitions of heterozygote advantage that largely cover this concept for multiple loci. Using extensive numerical simulations, we found no clear associations between the presence of any of the five definitions of heterozygote advantage and the maintenance of polymorphism at all loci. The strength of the association decreases as the number of loci increases or as recombination decreases. We conclude that heterozygote advantage cannot be a general mechanism for the maintenance of genetic polymorphism at multiple loci. These findings suggest that a correlation between the number of heterozygote loci and fitness is not warranted on theoretical ground.

Evaluation of adding a second marker to overcome Staphylococcus aureus spa typing homoplasies.

The utility of sequencing a second highly variable locus in addition to the spa gene (e.g., double-locus sequence typing [DLST]) was investigated to overcome limitations of a Staphylococcus aureus single-locus typing method. Although adding a second locus seemed to increase discriminatory power, it was not sufficient to definitively infer evolutionary relationships within a single multilocus sequence type (ST-5).