Universal Temporal Profile of Replication Origin Activation in Eukaryotes

Although replication proteins are conserved among eukaryotes, the sequence requirements for replication initiation differ between species. In all species, however, replication origins fire asynchronously throughout S phase. The temporal program of origin firing is reproducible in cell populations but largely probabilistic at the single-cell level. The mechanisms and the significance of this program are unclear. Replication timing has been correlated with gene activity in metazoans but not in yeast. One potential role for a temporal regulation of origin firing is to minimize fluctuations in replication end time and avoid persistence of unreplicated DNA in mitosis. Here, we have extracted the population-averaged temporal profiles of replication initiation rates for S. cerevisiae, S. pombe, D. melanogaster, X. laevis and H. sapiens from genome-wide replication timing and DNA combing data. All the profiles have a strikingly similar shape, increasing during the first half of S phase then decreasing before its end. A previously proposed minimal model of stochastic initiation modulated by accumulation of a recyclable, limiting replication-fork factor and fork-promoted initiation of new origins, quantitatively described the observed profiles without requiring new implementations

A Simple Model for the Influence of Meiotic Conversion Tracts on GC Content

A strong correlation between GC content and recombination rate is observed in many eukaryotes, which is thought to be due to conversion events linked to the repair of meiotic double-strand breaks. In several organisms, the length of conversion tracts has been shown to decrease exponentially with increasing distance from the sites of meiotic double-strand breaks. I show here that this behavior leads to a simple analytical model for the evolution and the equilibrium state of the GC content of sequences devoid of meiotic double-strand break sites. In the yeast Saccharomyces cerevisiae, meiotic double-strand breaks are practically excluded from protein-coding sequences. A good fit was observed between the predictions of the model and the variations of the average GC content of the third codon position (GC3) of S. cerevisiae genes. Moreover, recombination parameters that can be extracted by fitting the data to the model coincide with experimentally determined values. These results thus indicate that meiotic recombination plays an important part in determining the fluctuations of GC content in yeast coding sequences. The model also accounted for the different patterns of GC variations observed in the genes of Candida species that exhibit a variety of sexual lifestyles, and hence a wide range of meiotic recombination rates. Finally, the variations of the average GC3 content of human and chicken coding sequences could also be fitted by the model. These results suggest the existence of a widespread pattern of GC variation in eukaryotic genes due to meiotic recombination, which would imply the generality of two features of meiotic recombination: its association with GC-biased gene conversion and the quasi-exclusion of meiotic double-strand breaks from coding sequences. Moreover, the model points out to specific constraints on protein fragments encoded by exon terminal sequences, which are the most affected by the GC bias

The checkpoint Saccharomyces cerevisiae Rad9 protein contains a tandem tudor domain that recognizes DNA.

International audienceDNA damage checkpoints are signal transduction pathways that are activated after genotoxic insults to protect genomic integrity. At the site of DNA damage, 'mediator' proteins are in charge of recruiting 'signal transducers' to molecules 'sensing' the damage. Budding yeast Rad9, fission yeast Crb2 and metazoan 53BP1 are presented as mediators involved in the activation of checkpoint kinases. Here we show that, despite low sequence conservation, Rad9 exhibits a tandem tudor domain structurally close to those found in human/mouse 53BP1 and fission yeast Crb2. Moreover, this region is important for the resistance of Saccharomyces cerevisiae to different genotoxic stresses. It does not mediate direct binding to a histone H3 peptide dimethylated on K79, nor to a histone H4 peptide dimethylated on lysine 20, as was demonstrated for 53BP1. However, the tandem tudor region of Rad9 directly interacts with single-stranded DNA and double-stranded DNAs of various lengths and sequences through a positively charged region absent from 53BP1 and Crb2 but present in several yeast Rad9 homologs. Our results argue that the tandem tudor domains of Rad9, Crb2 and 53BP1 mediate chromatin binding next to double-strand breaks. However, their modes of chromatin recognition are different, suggesting that the corresponding interactions are differently regulated

Genome-Wide Crossover Distribution in Arabidopsis thaliana Meiosis Reveals Sex-Specific Patterns along Chromosomes

In most species, crossovers (COs) are essential for the accurate segregation of homologous chromosomes at the first meiotic division. Their number and location are tightly regulated. Here, we report a detailed, genome-wide characterization of the rate and localization of COs in Arabidopsis thaliana, in male and female meiosis. We observed dramatic differences between male and female meiosis which included: (i) genetic map length; 575 cM versus 332 cM respectively; (ii) CO distribution patterns: male CO rates were very high at both ends of each chromosome, whereas female CO rates were very low; (iii) correlations between CO rates and various chromosome features: female CO rates correlated strongly and negatively with GC content and gene density but positively with transposable elements (TEs) density, whereas male CO rates correlated positively with the CpG ratio. However, except for CpG, the correlations could be explained by the unequal repartition of these sequences along the Arabidopsis chromosome. For both male and female meiosis, the number of COs per chromosome correlates with chromosome size expressed either in base pairs or as synaptonemal complex length. Finally, we show that interference modulates the CO distribution both in male and female meiosis

The Sex-Specific Impact of Meiotic Recombination on Nucleotide Composition

Meiotic recombination is an important evolutionary force shaping the nucleotide landscape of genomes. For most vertebrates, the frequency of recombination varies slightly or considerably between the sexes (heterochiasmy). In humans, male, rather than female, recombination rate has been found to be more highly correlated with the guanine and cytosine (GC) content across the genome. In the present study, we review the results in human and extend the examination of the evolutionary impact of heterochiasmy beyond primates to include four additional eutherian mammals (mouse, dog, pig, and sheep), a metatherian mammal (opossum), and a bird (chicken). Specifically, we compared sex-specific recombination rates (RRs) with nucleotide substitution patterns evaluated in transposable elements. Our results, based on a comparative approach, reveal a great diversity in the relationship between heterochiasmy and nucleotide composition. We find that the stronger male impact on this relationship is a conserved feature of human, mouse, dog, and sheep. In contrast, variation in genomic GC content in pig and opossum is more strongly correlated with female, rather than male, RR. Moreover, we show that the sex-differential impact of recombination is mainly driven by the chromosomal localization of recombination events. Independent of sex, the higher the RR in a genomic region and the longer this recombination activity is conserved in time, the stronger the bias in nucleotide substitution pattern, through such mechanisms as biased gene conversion. Over time, this bias will increase the local GC content of the region

Late Replicating Domains Are Highly Recombining in Females but Have Low Male Recombination Rates: Implications for Isochore Evolution

In mammals sequences that are either late replicating or highly recombining have high rates of evolution at putatively neutral sites. As early replicating domains and highly recombining domains both tend to be GC rich we a priori expect these two variables to covary. If so, the relative contribution of either of these variables to the local neutral substitution rate might have been wrongly estimated owing to covariance with the other. Against our expectations, we find that sex-averaged recombination rates show little or no correlation with replication timing, suggesting that they are independent determinants of substitution rates. However, this result masks significant sex-specific complexity: late replicating domains tend to have high recombination rates in females but low recombination rates in males. That these trends are antagonistic explains why sex-averaged recombination is not correlated with replication timing. This unexpected result has several important implications. First, although both male and female recombination rates covary significantly with intronic substitution rates, the magnitude of this correlation is moderately underestimated for male recombination and slightly overestimated for female recombination, owing to covariance with replicating timing. Second, the result could explain why male recombination is strongly correlated with GC content but female recombination is not. If to explain the correlation between GC content and replication timing we suppose that late replication forces reduced GC content, then GC promotion by biased gene conversion during female recombination is partly countered by the antagonistic effect of later replicating sequence tending increase AT content. Indeed, the strength of the correlation between female recombination rate and local GC content is more than doubled by control for replication timing. Our results underpin the need to consider sex-specific recombination rates and potential covariates in analysis of GC content and rates of evolution

Quantitative Trait Locus (QTL) Mapping Reveals a Role for Unstudied Genes in Aspergillus Virulence

Infections caused by the fungus Aspergillus are a major cause of morbidity and mortality in immunocompromised populations. To identify genes required for virulence that could be used as targets for novel treatments, we mapped quantitative trait loci (QTL) affecting virulence in the progeny of a cross between two strains of A. nidulans (FGSC strains A4 and A91). We genotyped 61 progeny at 739 single nucleotide polymorphisms (SNP) spread throughout the genome, and constructed a linkage map that was largely consistent with the genomic sequence, with the exception of one potential inversion of ∼527 kb on Chromosome V. The estimated genome size was 3705 cM and the average intermarker spacing was 5.0 cM. The average ratio of physical distance to genetic distance was 8.1 kb/cM, which is similar to previous estimates, and variation in recombination rate was significantly positively correlated with GC content, a pattern seen in other taxa. To map QTL affecting virulence, we measured the ability of each progeny strain to kill model hosts, larvae of the wax moth Galleria mellonella. We detected three QTL affecting in vivo virulence that were distinct from QTL affecting in vitro growth, and mapped the virulence QTL to regions containing 7–24 genes, excluding genes with no sequence variation between the parental strains and genes with only synonymous SNPs. None of the genes in our QTL target regions have been previously associated with virulence in Aspergillus, and almost half of these genes are currently annotated as “hypothetical”. This study is the first to map QTL affecting the virulence of a fungal pathogen in an animal host, and our results illustrate the power of this approach to identify a short list of unknown genes for further investigation

Concerted cutting by Spo11 illuminates meiotic DNA break mechanics

Genetic recombination arises during meiosis through the repair of DNA double-strand breaks (DSBs) that are created by Spo11, a topoisomerase-like protein1,2. Spo11 DSBs form preferentially in nucleosome-depleted regions termed hotspots3,4, yet how Spo11 engages with its DNA substrate to catalyse DNA cleavage is poorly understood. Although most recombination events are initiated by a single Spo11 cut, here we show in Saccharomyces cerevisiae that hyperlocalized, concerted Spo11 DSBs separated by 33 to more than 100 base pairs also form, which we term ‘double cuts’. Notably, the lengths of double cuts vary with a periodicity of 10.5 base pairs, which is conserved in yeast and mice. This finding suggests a model in which the orientation of adjacent Spo11 molecules is fixed relative to the DNA helix—a proposal supported by the in vitro DNA-binding properties of the Spo11 core complex. Deep sequencing of meiotic progeny identifies recombination scars that are consistent with repair initiated from gaps generated by adjacent Spo11 DSBs. Collectively, these results revise our present understanding of the mechanics of Spo11-DSB formation and expand on the original concepts of gap repair during meiosis to include DNA gaps that are generated by Spo11 itself

Schematics illustrating the computation of asymmetry indices.

<p>(A) Definition of the four strand segments. (B) Four simple cases are shown for different positions <i>x</i>. Segments synthesized as leading and lagging strands are shown as blue and red lines, respectively. The small, black rectangles represent coding sequences with their transcriptional orientation given by the associated arrows. The large, black arrowheads indicate the orientation of the replication forks.</p

Variations of GC3 content in <i>S. cerevisiae</i> protein-coding sequences.

<p>The mean GC3 contents (A) or (B) of 66-codon segments are plotted as a function of the positions or relative to the ATG or to the stop codon, respectively, on which the segments are centered. For a given position or , and are averaged over classes of genes binned by their lengths. The genes of set 1 (, blue dotted line) contain 3 to 5 66-codon segments, the genes of set 2 (, blue dashed line) contain 6 to 8 segments, the genes of set 3 (, blue dot-dash line) contain 9 to 11 segments, the genes of set 4 (, blue long-dash line) contain 12 to 14 segments, and the genes of set 5 (, blue solid line) contain at least 15 segments. Only the average values for the segments common to all the genes of a given set are plotted (<i>i.e.</i> only the segments corresponding to the shortest genes of the set). The thick, red, solid lines represent the mean or averaged over all genes containing at least one 66-codon segment (). The thin, red, solid lines represent the mean theoretical values or averaged over all genes containing at least one 66-codon segment and whose values of and could be determined (). and were calculated as functions of , , and using Equations 7 and 8 (and the estimates of , , and given in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0016109#pone-0016109-t001" target="_blank">Table 1</a>), and averaged over segments with the same position or , respectively.</p