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

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

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

Theileria parasites secrete a prolyl isomerase to maintain host leukocyte transformation

Infectious agents develop intricate mechanisms to interact with host cell pathways and hijack the genetic and epigenetic machinery to change phenotypic states. Amongst the Apicomplexa phylum of obligate intracellular parasites which cause veterinary and human diseases, Theileria is the only genus which transforms its mammalian host cells1. Theileria infection of bovine leukocytes induces proliferative and invasive phenotypes associated with activated signalling pathways, notably JNK and AP-12. The transformed phenotypes are reversed by treatment with the theilericidal drug Buparvaquone3. We used comparative genomics to identify a homologue of the Peptidyl Prolyl Isomerase Pin1 (designated TaPin1) in T. annulata which is secreted into the host cell and modulates oncogenic signalling pathways. Here we show that TaPin1 is a bona fide prolyl isomerase and that it interacts with the host ubiquitin ligase FBW7 leading to its degradation and subsequent stabilization of c-Jun which promotes transformation. We performed in vitro analysis and in vivo zebrafish xenograft experiments to demonstrate that TaPin1 is directly inhibited by the anti-parasite drug Buparvaquone (and other known Pin1 inhibitors) and is mutated in a drug-resistant strain. Prolyl isomerisation is thus a conserved mechanism which is important in cancer and is used by Theileria parasites to manipulate host oncogenic signaling