Linking the DNA strand asymmetry to the spatio-temporal replication program: II. Accounting for neighbor-dependent substitution rates.

In paper I, we addressed the impact of the spatio-temporal program on the DNA composition evolution in the case of time homogeneous and neighbor-independent substitution rates. But substitution rates do depend on the flanking nucleotides as exemplified in vertebrates where CpG sites are hypermutable so that the substitution rate C --> T depends dramatically (ten fold) on whether the cytosine belongs to a CG dinucleotide or not. With the specific goal to account for neighbor-dependence, we revisit our minimal modeling of neutral substitution rates in the human genome. When assuming that r = CpG --> TpG and its reverse complement r(c) = CpG --> CpA are (by far) the main neighbor-dependent substitution rates, we demonstrate, using perturbative analysis, that neighbor-dependence does not affect the decomposition of the compositional asymmetry into a transcription- and a replication-associated components, the former increases in magnitude with transcription rate and changes sign with gene orientation, whereas the latter is proportional to the replication fork polarity. Indeed the neighbor dependence case differs from the neighbor-independent model by an additional source term related to the CG dinucleotide content in both the transcription and replication-associated components. We finally discuss the case of time-dependent substitution rates confirming as a very general result the fact that the skew can still be decomposed into a transcription- and a replication-associated components.

Linking the DNA strand asymmetry to the spatio-temporal replication program. I. About the role of the replication fork polarity in genome evolution.

Two key cellular processes, namely transcription and replication, require the opening of the DNA double helix and act differently on the two DNA strands, generating different mutational patterns (mutational asymmetry) that may result, after long evolutionary time, in different nucleotide compositions on the two DNA strands (compositional asymmetry). We elaborate on the simplest model of neutral substitution rates that takes into account the strand asymmetries generated by the transcription and replication processes. Using perturbation theory, we then solve the time evolution of the DNA composition under strand-asymmetric substitution rates. In our minimal model, the compositional and substitutional asymmetries are predicted to decompose into a transcription- and a replication-associated components. The transcription-associated asymmetry increases in magnitude with transcription rate and changes sign with gene orientation while the replication-associated asymmetry is proportional to the replication fork polarity. These results are confirmed experimentally in the human genome, using substitution rates obtained by aligning the human and chimpanzee genomes using macaca and orangutan as outgroups, and replication fork polarity determined in the HeLa cell line as estimated from the derivative of the mean replication timing. When further investigating the dynamics of compositional skew evolution, we show that it is not at equilibrium yet and that its evolution is an extremely slow process with characteristic time scales of several hundred Myrs.

XUTs are a class of Xrn1-sensitive antisense regulatory non-coding RNA in yeast.

Non-coding (nc)RNAs are key players in numerous biological processes such as gene regulation, chromatin domain formation and genome stability. Large ncRNAs interact with histone modifiers and are involved in cancer development, X-chromosome inactivation and autosomal gene imprinting. However, despite recent evidence showing that pervasive transcription is more widespread than previously thought, only a few examples mediating gene regulation in eukaryotes have been described. In Saccharomyces cerevisiae, the bona-fide regulatory ncRNAs are destabilized by the Xrn1 5'-3' RNA exonuclease (also known as Kem1), but the genome-wide characterization of the entire regulatory ncRNA family remains elusive. Here, using strand-specific RNA sequencing (RNA-seq), we identify a novel class of 1,658 Xrn1-sensitive unstable transcripts (XUTs) in which 66% are antisense to open reading frames. These transcripts are polyadenylated and RNA polymerase II (RNAPII)-dependent. The majority of XUTs strongly accumulate in lithium-containing media, indicating that they might have a role in adaptive responses to changes in growth conditions. Notably, RNAPII chromatin immunoprecipitation followed by DNA sequencing (ChIP-seq) analysis of Xrn1-deficient strains revealed a significant decrease of RNAPII occupancy over 273 genes with antisense XUTs. These genes show an unusual bias for H3K4me3 marks and require the Set1 histone H3 lysine 4 methyl-transferase for silencing. Furthermore, abolishing H3K4me3 triggers the silencing of other genes with antisense XUTs, supporting a model in which H3K4me3 antagonizes antisense ncRNA repressive activity. Our results demonstrate that antisense ncRNA-mediated regulation is a general regulatory pathway for gene expression in S. cerevisiae.

Bifractality of human DNA strand-asymmetry profiles results from transcription.

We use the wavelet transform modulus maxima method to investigate the multifractal properties of strand-asymmetry DNA walk profiles in the human genome. This study reveals the bifractal nature of these profiles, which involve two competing scale-invariant (up to repeat-masked distances less, or similar 40 kbp) components characterized by Hölder exponents h{1}=0.78 and h{2}=1, respectively. The former corresponds to the long-range-correlated homogeneous fluctuations previously observed in DNA walks generated with structural codings. The latter is associated with the presence of jumps in the original strand-asymmetry noisy signal S. We show that a majority of upward (downward) jumps co-locate with gene transcription start (end) sites. Here 7228 human gene transcription start sites from the refGene database are found within 2 kbp from an upward jump of amplitude DeltaS > or = 0.1 which suggests that about 36% of annotated human genes present significant transcription-induced strand asymmetry and very likely high expression rate.

DNA replication timing data corroborate in silico human replication origin predictions.

We develop a wavelet-based multiscale pattern recognition methodology to disentangle the replication- from the transcription-associated compositional strand asymmetries observed in the human genome. Comparing replication skew profiles to recent high-resolution replication timing data reveals that most of the putative replication origins that border the so-identified replication domains are replicated earlier than their surroundings whereas the central regions replicate late in the S phase. We discuss the implications of this first experimental confirmation of these replication origin predictions that are likely to be early replicating and active in most tissues.

From DNA sequence analysis to modeling replication in the human genome.

We explore the large-scale behavior of nucleotide compositional strand asymmetries along human chromosomes. As we observe for 7 of 9 origins of replication experimentally identified so far, the (TA+GC) skew displays rather sharp upward jumps, with a linear decreasing profile in between two successive jumps. We present a model of replication with well positioned replication origins and random terminations that accounts for the observed characteristic serrated skew profiles. We succeed in identifying 287 pairs of putative adjacent replication origins with an origin spacing approximately 1-2 Mbp that are likely to correspond to replication foci observed in interphase nuclei and recognized as stable structures that persist throughout subsequent cell generations.

Low frequency rhythms in human DNA sequences: a key to the organization of gene location and orientation?

We explore large-scale nucleotide compositional fluctuations of the human genome using multiresolution techniques. Analysis of the GC content and of the AT and GC skews reveals the existence of rhythms with two main periods of 110+/-20 kb and 400+/-50 kb that enlighten a remarkable cooperative gene organization. We show that the observed nonlinear oscillations are likely to display all the characteristic features of chaotic strange attractors which suggests a very attractive deterministic picture: gene orientation and location, in relation with the structure and dynamics of chromatin, might be governed by a low-dimensional nonlinear dynamical system.

Transcription-coupled TA and GC strand asymmetries in the human genome.

Analysis of the whole set of human genes reveals that most of them present TA and GC skews, that these biases are correlated to each other and are specific to gene sequences, exhibiting sharp transitions between transcribed and non-transcribed regions. The GC asymmetries cannot be explained solely by a model previously proposed for (G+T) skew based on transitions measured in a small set of human genes. We propose that the GC skew results from additional transcription-coupled mutation process that would include transversions. During evolution, both processes acting on a large majority of genes in germline cells would have produced these transcription-coupled strand asymmetries.

Long-range correlations in genomic DNA: a signature of the nucleosomal structure.

We use the "wavelet transform microscope" to carry out a comparative statistical analysis of DNA bending profiles and of the corresponding DNA texts. In the three kingdoms, one reveals on both signals a characteristic scale of 100-200 bp that separates two different regimes of power-law correlations (PLC). In the small-scale regime, PLC are observed in eukaryotic, in double-strand DNA viral, and in archaeal genomes, which contrasts with their total absence in the genomes of eubacteria and their viruses. This strongly suggests that small-scale PLC are related to the mechanisms underlying the wrapping of DNA in the nucleosomal structure. We further speculate that the large scale PLC are the signature of the higher-order structure and dynamics of chromatin.

Bigfoot. a new family of MITE elements characterized from the Medicago genus.

We have characterized from the legume plant Medicago a new family of miniature inverted-repeat transposable elements (MITE), called the Bigfoot transposable elements. Two of these insertion elements are present only in a single allele of two different M. sativa genes. Using a PCR strategy we have isolated 19 other Bigfoot elements from the M. sativa and M. truncatula genomes. They differ from the previously characterized MITEs by their sequence, a target site of 9 bp and a partially clustered genomic distribution. In addition, we show that they exhibit a significantly stable secondary structure. These elements may represent up to 0.1% of the genome of the outcrossing Medicago sativa but are present at a reduced copy number in the genome of the autogamous M. truncatula plant, revealing major differences in the genome organization of these two plants.

A multi-agent system simulating human splice site recognition.

The present paper describes a method detecting splice sites automatically on the basis of sequence data and models of site/signal recognition supported by experimental evidences. The method is designed to simulate splicing and while doing so, track prediction failures, missing information and possibly test correcting hypotheses. Correlations between nucleotides in the splice site regions and the various elements of the acceptor region are evaluated and combined to assess compensating interactions between elements of the splicing machinery. A scanning model of the acceptor region and a model of interaction between the splicing complexes (exon definition model) are also incorporated in the detection process. Subsets of sites presenting deficiencies of several splice site elements could be identified. Further examination of these sites helps to determine lacking elements and refine models.

Two co-existing mechanisms account for the large-scale deletions of mitochondrial DNA in Podospora anserina that involve the 5' border of a group-II intron.

A degenerative syndrome associated with the accumulation of site-specific deletions within mitochondrial chromosomes occurs in strains of Podospora anserina carrying the AS1-4 nuclear mutation. The site-specific deletion event has been assumed to result from the transposition of a group-II intron (intron alpha) behind an IBS motif, followed by recombination between the two intron repeats. We show here that a number of distinct deletions can accumulate in AS1-4 strains. Most of them are present in low amounts in wild-type cells where they are only detectable in PCR experiments. The deletions can be divided into two classes. In class I, intron alpha is joined to an IBS motif. In class II, the intron is not joined to an IBS site, it can be truncated or contain a few upstream exonic nucleotides; some junctions carry non-templated nucleotides. These results indicate that at least two mechanisms are involved in the generation of large-scale mitochondrial deletions in Podospora. One of them seems to be based on the transposition properties of the group-II alpha intron, the other one on illegitimate recombination. We propose that these two mechanisms use DNA double-strand breaks occurring within the 5' region of intron alpha.

Mitochondrial intronic open reading frames in Podospora: mobility and consecutive exonic sequence variations.

The mitochondrial genome of 23 wild-type strains belonging to three different species of the filamentous fungus Podospora was examined. Among the 15 optional sequences identified are two intronic reading frames, nad1-i4-orf1 and cox1-i7-orf2. We show that the presence of these sequences was strictly correlated with tightly clustered nucleotide substitutions in the adjacent exon. This correlation applies to the presence or absence of closely related open reading frames (ORFs), found at the same genetic locations, in all the Pyrenomycete genera examined. The recent gain of these optional ORFs in the evolution of the genus Podospora probably account for such sequence differences. In the homoplasmic progeny from heteroplasmons constructed between Podospora strains differing by the presence of these optional ORFs, nad1-i4-orf1 and cox1-i7-orf2 appeared highly invasive. Sequence comparisons in the nad1-i4 intron of various strains of the Pyrenomycete family led us to propose a scenario of its evolution that includes several events of loss and gain of intronic ORFs. These results strongly reinforce the idea that group 1 intronic ORFs are mobile elements and that their transfer, and concomitant modification of the adjacent exon, could participate in the modular evolution of mitochondrial genomes.

Contribution of ultra-short invasive elements to the evolution of the mitochondrial genome in the genus Podospora.

In the filamentous fungus Podospora anserina, senescence is associated with major rearrangements of the mitochondrial DNA. The undecamer GGCGCAAGCTC has been described as a preferential site for these recombination events. We show that: (i) copies of this short sequence GGCGCAAGCTC are present in unexpectedly high numbers in the mitochondrial genome of this fungus; (ii) a short cluster of this sequence, localised in a group II intronic ORF, encodes amino acids that disrupt a protein domain that is otherwise highly conserved between various species; (iii) most of the polymorphisms observed between three related species, P.anserina, P.curvicolla and P.comata, are associated with the presence/absence of this sequence; (iv) this element lies at the boundaries of major rearrangements of the mitochondrial genomes; (v) at least two other short elements in the Podospora mitochondrial genomes display similar features. We suggest that these short elements, called MUSEs (mitochondrial ultra-short elements), could be mobile and that they contribute to evolution of the mitochondrial genome in the genus Podospora. A model for mobility involving a target DNA-primed reverse transcription step is discussed.

Tissue-specific splicing of two mutually exclusive exons of the chicken beta-tropomyosin pre-mRNA: positive and negative regulations.

Alternative splicing of premessenger RNA (pre-mRNA) is a widespread process used in higher eucaryotes to regulate gene expression. A single primary transcript can generate multiple proteins with distinct functions in a tissue- and/or developmental-specific manner. A central question in alternative splicing concerns the selection of splice sites in different cell environments. In this review, we present our results on the alternative splicing of the chicken beta-tropomyosin gene which provides an interesting model for understanding mechanisms involved in splice site selection. The beta-tropomyosin gene contains in its central portion a pair of exons (6A and 6B) that are used mutually exclusively in a tissue and developmental stage-specific manner. Exon 6A is present in mRNA of non-muscle and smooth muscle tissues while exon 6B is only present in mRNA of skeletal muscle. Regulation of both exons is necessary to ensure specific expression of beta-tropomyosin gene in non-muscle cells. Several cis-acting elements involved in the repression of exon 6B and activation of exon 6A have been identified. In addition, we show that the tissue-specific choice of exon 6A is mediated through interaction with a specific class of splicing factors, the SR proteins. In the last part of this review we will focus on possible mechanisms needed to switch to exon 6B selection in skeletal muscle tissue. We propose that tissue-specific choice of exon 6B involves down regulation of exon 6A and activation of exon 6B. A G-rich enhancer sequence downstream of exon 6B has been defined that is needed for efficient recognition of the exon 6B 5' splice site. Moreover, we suggest that alteration of the ratio between proteins of the SR family contributes to tissue-specific splice site selection.

Exon prediction in eucaryotic genomes.

Two independent computer systems, NetPlantGene and AMELIE, dedicated to the identification of splice sites in plant and human genomes, respectively, are introduced here. Both methods were designed in relation to experimental work; they rely on automatically generated rules involving the nucleotide content of sequences regardless of the coding properties of exons. The specificity of plant sequences as considered in NetPlantGene is shown to enhance the quality of detection as opposed to general methods such as GRAIL. A scanning model of the acceptor site recognition is being simulated by AMELIE leading to a relatively accurate selection process of sites.

In vitro splicing deficiency induced by a C to T mutation at position -3 in the intron 10 acceptor site of the phenylalanine hydroxylase gene in a patient with phenylketonuria.

A previous study has identified a C-->U mutation at position -3 in the 3' splice site of intron 10 of the phenylalanine hydroxylase pre-mRNA in a patient with phenylketonuria. In vivo, this mutation induces the skipping of the downstream exon. This result is puzzling because both CAG and UAG have been reported to function equally as 3' splice sites. In this report, we show that the C-->U mutation affects predominantly the first step of the splicing reaction and that it blocks spliceosome assembly at an early stage. The 3' region of the phenylalanine hydroxylase intron 10 has two unusual characteristic features: multiple potential branch sites and a series of four guanosine residues, which interrupt the polypyrimidine tract at positions -8 to -11 from the 3' splice site. We show that the mutation precludes the use of the proximal branch site, while having no effect on the remote one. We also show that in the UAG transcript, the four guanosine residues inhibit the splicing of intron 10. The substitution of these purine residues by one cytosine residue, regardless of the position, increases the splicing efficiency of the mutant UAG precursor while having no effect on the wild-type CAG precursor. Substituting the four purine residues by four pyrimidines relieves the inhibition and rescues the use of the proximal branch site. These results demonstrate that according to the context, the C and U nucleotides preceding the AG are not equivalent for the splicing reaction.

enod40, a gene expressed during nodule organogenesis, codes for a non-translatable RNA involved in plant growth.

Rhizobium meliloti can interact symbiotically with Medicago plants, thereby inducing root nodules. However, certain Medicago plants can form nodules spontaneously, in the absence of rhizobia. A differential screening was performed using spontaneous nodule versus root cDNAs from Medicago sativa ssp. varia. Transcripts of a differentially expressed clone, Msenod40, were detected in all differentiating cells of nodule primordia and spontaneous nodules, but were absent in fully differentiated cells. Msenod40 showed homology to a soybean early nodulin gene, Gmenod40, although no significant open reading frame (ORF) or coding capacity was found in the Medicago sequence. Furthermore, in the sequences of cDNAs and a genomic clone (Mtenod40) isolated from Medicago truncatula, a species containing a unique copy of this gene, no ORFs were found either. In vitro translation of purified Mtenod40 transcripts did not reveal any protein product. Evaluation of the RNA secondary structure indicated that both msenod40 and Gmenod40 transcripts showed a high degree of stability, a property shared with known non-coding RNAs. The Mtenod40 RNA was localized in the cytoplasm of cells in the nodule primordium. Infection with Agrobacterium tumefaciens strains bearing antisense constructs of Mtenod40 arrested callus growth of Medicago explants, while overexpressing Mtenod40 embryos developed into teratomas. These data suggest that the enod40 genes might have a role in plant development, acting as 'riboregulators', a novel class of untranslated RNAs associated with growth control and differentiation.

Determination of an RNA structure involved in splicing inhibition of a muscle-specific exon.

We have investigated the RNA structure of the region surrounding the muscle-specific exon 6B of the chicken beta-tropomyosin gene. We have used a variety of chemical and enzymatic probes: dimethylsulfate, N-cyclohexyl-N'-(2-(N-methylmorpholino)-ethyl)-carbodiimide-p-tolu enesulfonate) , RNase T1 and RNase V1. Lead acetate was also used to obtain some information on the tertiary structure of this region. Probing the wild-type sequence suggests a model involving one-stem and three-stem-loop structures in and around this exon. Two of these, hairpin I and stem III, have previously been implicated in repression of splicing of the intron following exon 6B in a HeLa nuclear extract. Stem I includes sequences at the beginning of exon 6B and stem III results from interaction of the intron upstream from exon 6B with sequences in the middle of the intron downstream from this exon (the intron whose splicing is repressed). Neither stem I nor stem III directly involves the consensus sequences (5' splice site, branch-point, 3' splice site) of the repressed intron. Probing RNAs that are derepressed for splicing of this intron show that there are structural changes around the 5' splice site and branch-point sequence that correlate with the derepression. This is true, despite the fact that the derepressed RNAs are altered in a region far from these consensus sequences. The most striking structural correlation with splicing capacity of the intron downstream from exon 6B is seen by probing with lead acetate. Lead ions cut RNA at specific residues; these sites are very sensitive to RNA tertiary structure. Repressed and derepressed RNAs show entirely different cleavage patterns after incubation with lead acetate. Remarkably, hybridizing a derepressed RNA with an RNA comprising the ascending arm of stem III not only re-establishes repression, but also converts the pattern of susceptibility to attack by lead ions over the whole molecule. We suggest that RNA conformation plays a role in keeping exon 6B from being spliced into non-muscle cell mRNA.

RNA secondary structure repression of a muscle-specific exon in HeLa cell nuclear extracts.

The chicken beta-tropomyosin pre-messenger RNA (pre-mRNA) is spliced in a tissue-specific manner to yield messenger RNA's (mRNA's) coding for different isoforms of this protein. Exons 6A and 6B are spliced in a mutually exclusive manner; exon 6B was included in skeletal muscle, whereas exon 6A was preferred in all other tissues. The distal portion of the intron upstream of exon 6B was shown to form stable double-stranded regions with part of the intron downstream of exon 6B and with sequences in exon 6B. This structure repressed splicing of exon 6B to exon 7 in a HeLa cell extract. Derepression of splicing occurred on disruption of this structure and repression followed when the structure was re-formed, even if the structure was formed between two different RNA molecules. Repression leads to inhibition of formation of spliceosomes. Disrupting either of the two double-stranded regions could lead to derepression, whereas re-forming the helices by suppressor mutations reestablished repression. These results support a simple model of tissue-specific splicing in this region of the pre-mRNA.