REDHORSE-REcombination and Double crossover detection in Haploid Organisms using next-geneRation SEquencing data

Background

Next-generation sequencing technology provides a means to study genetic exchange at a higher resolution than was possible using earlier technologies. However, this improvement presents challenges as the alignments of next generation sequence data to a reference genome cannot be directly used as input to existing detection algorithms, which instead typically use multiple sequence alignments as input. We therefore designed a software suite called REDHORSE that uses genomic alignments, extracts genetic markers, and generates multiple sequence alignments that can be used as input to existing recombination detection algorithms. In addition, REDHORSE implements a custom recombination detection algorithm that makes use of sequence information and genomic positions to accurately detect crossovers. REDHORSE is a portable and platform independent suite that provides efficient analysis of genetic crosses based on Next-generation sequencing data.

Results

We demonstrated the utility of REDHORSE using simulated data and real Next-generation sequencing data. The simulated dataset mimicked recombination between two known haploid parental strains and allowed comparison of detected break points against known true break points to assess performance of recombination detection algorithms. A newly generated NGS dataset from a genetic cross of Toxoplasma gondii allowed us to demonstrate our pipeline. REDHORSE successfully extracted the relevant genetic markers and was able to transform the read alignments from NGS to the genome to generate multiple sequence alignments. Recombination detection algorithm in REDHORSE was able to detect conventional crossovers and double crossovers typically associated with gene conversions whilst filtering out artifacts that might have been introduced during sequencing or alignment. REDHORSE outperformed other commonly used recombination detection algorithms in finding conventional crossovers. In addition, REDHORSE was the only algorithm that was able to detect double crossovers.

Conclusion

REDHORSE is an efficient analytical pipeline that serves as a bridge between genomic alignments and existing recombination detection algorithms. Moreover, REDHORSE is equipped with a recombination detection algorithm specifically designed for Next-generation sequencing data. REDHORSE is portable, platform independent Java based utility that provides efficient analysis of genetic crosses based on Next-generation sequencing data. REDHORSE is available at http://redhorse.sourceforge.net/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1309-7) contains supplementary material, which is available to authorized users.     

The Theory and Applications of Measuring Broad-Range and Chromosome-Wide Recombination Rate from Allele Frequency Decay around a Selected Locus

Recombination is the exchange of genetic material between homologous chromosomes via physical crossovers. High-throughput sequencing approaches detect crossovers genome wide to produce recombination rate maps but are difficult to scale as they require large numbers of recombinants individually sequenced. We present a simple and scalable pooled-sequencing approach to experimentally infer near chromosome-wide recombination rates by taking advantage of non-Mendelian allele frequency generated from a fitness differential at a locus under selection. As more crossovers decouple the selected locus from distal loci, the distorted allele frequency attenuates distally toward Mendelian and can be used to estimate the genetic distance. Here, we use marker selection to generate distorted allele frequency and theoretically derive the mathematical relationships between allele frequency attenuation, genetic distance, and recombination rate in marker-selected pools. We implemented nonlinear curve-fitting methods that robustly estimate the allele frequency decay from batch sequencing of pooled individuals and derive chromosome-wide genetic distance and recombination rates. Empirically, we show that marker-selected pools closely recapitulate genetic distances inferred from scoring recombinants. Using this method, we generated novel recombination rate maps of three wild-derived strains of Drosophila melanogaster, which strongly correlate with previous measurements. Moreover, we show that this approach can be extended to estimate chromosome-wide crossover interference with reciprocal marker selection and discuss how it can be applied in the absence of visible markers. Altogether, we find that our method is a simple and cost-effective approach to generate chromosome-wide recombination rate maps requiring only one or two libraries.     

Propagation of Biochirality: Crossovers and Nonclassical Crystallization Kinetics of Aspartic Acid in Water

All experimental procedures discussed could be treated as a screening tool for probing the existence of molecular association among the chiral molecules and the solvent system. The molecular association phases of a racemic conglomerate solution (CS) and a racemic compound solution (RCS), and the templating effect of aspartic acid solid surface were observed to minimize the chance of redissolving racemic conglomerate and racemic compound aspartic acid in water and reforming an RCS in crossovers experiments. Only 1 %wt% of l‐aspartic acid was adequate enough to induce a transformation from a racemic compound aspartic acid to a racemic conglomerate aspartic acid. This would make the propagation of biochirality more feasible and sound. However, tetrapeptide, (l‐aspartic acid)₄, failed to induce enantioseparation as templates purely by crystallization. Nonclassical crystallization theory was needed to take into account the existence of a CS. Fundamental parameters of the crystallization kinetics such as the induction time, interfacial energy, Gibbs energetic barrier, nucleation rate, and critical size of stable nuclei of: (i) racemic compound aspartic acid, (ii) racemic compound aspartic acid seeded with 1 %wt% l‐aspartic acid, (iii) racemic conglomerate aspartic acid, and (iv) l‐aspartic acid were evaluated and compared with different initial supersaturation ratios. Morphological studies of crystals grown from the crystallization kinetics were also carried out.Chirality 25:768–779, 2013. © 2013 Wiley Periodicals, Inc.     

Engineering meiotic recombination pathways in rice

In the last 15 years, outstanding progress has been made in understanding the function of meiotic genes in the model dicot and monocot plants Arabidopsis and rice (Oryza sativa L.), respectively. This knowledge allowed to modulate meiotic recombination in Arabidopsis and, more recently, in rice. For instance, the overall frequency of crossovers (COs) has been stimulated 2.3‐ and 3.2‐fold through the inactivation of the rice FANCM and RECQ4 DNA helicases, respectively, two genes involved in the repair of DNA double‐strand breaks (DSBs) as noncrossovers (NCOs) of the Class II crossover pathway. Differently, the programmed induction of DSBs and COs at desired sites is currently explored by guiding the SPO11‐1 topoisomerase‐like transesterase, initiating meiotic recombination in all eukaryotes, to specific target regions of the rice genome. Furthermore, the inactivation of 3 meiosis‐specific genes, namely PAIR1, OsREC8 and OsOSD1, in the Mitosis instead of Meiosis (MiMe) mutant turned rice meiosis into mitosis, thereby abolishing recombination and achieving the first component of apomixis, apomeiosis. The successful translation of Arabidopsis results into a crop further allowed the implementation of two breakthrough strategies that triggered parthenogenesis from the MiMe unreduced clonal egg cell and completed the second component of diplosporous apomixis. Here, we review the most recent advances in and future prospects of the manipulation of meiotic recombination in rice and potentially other major crops, all essential for global food security.     

Where are the frame-shifted gene ‘duplications’?

Evolution by gene duplication has been well documented, mostly by the discovery of two or more similar amino-acid sequences in proteins translated from one haploid genome. The sequences differ by interpolation and deletion of amino-acids as well as substitution. This is evidence for triplet (codon) addition or subtraction; do deletions/interpolations of one or two nucleotides (frame-shifting mutations) occur? If they are not seen (as in haemoglobins) either this is because they are not generated (the genetic code constrains DNA events); they are generated but edited (gene conversion?); gametes containing them are selected against; or they are actively and rapidly selected against evolutionarily. If they are seen (as in immunoglobulins, possibly) they may generate somatic diversity (as in antibodies) and there may be a limbo of more-or-less changed structural sequences, returning to translation (and selection) only by another frame-shifting event. This suggestion combines neutral mutation theory with pan-selectionism.     

Sources and Structures of Mitotic Crossovers That Arise When BLM Helicase Is Absent in Drosophila

The Bloom syndrome helicase, BLM, has numerous functions that prevent mitotic crossovers. We used unique features of Drosophila melanogaster to investigate origins and properties of mitotic crossovers that occur when BLM is absent. Induction of lesions that block replication forks increased crossover frequencies, consistent with functions for BLM in responding to fork blockage. In contrast, treatment with hydroxyurea, which stalls forks, did not elevate crossovers, even though mutants lacking BLM are sensitive to killing by this agent. To learn about sources of spontaneous recombination, we mapped mitotic crossovers in mutants lacking BLM. In the male germline, irradiation-induced crossovers were distributed randomly across the euchromatin, but spontaneous crossovers were nonrandom. We suggest that regions of the genome with a high frequency of mitotic crossovers may be analogous to common fragile sites in the human genome. Interestingly, in the male germline there is a paucity of crossovers in the interval that spans the pericentric heterochromatin, but in the female germline this interval is more prone to crossing over. Finally, our system allowed us to recover pairs of reciprocal crossover chromosomes. Sequencing of these revealed the existence of gene conversion tracts and did not provide any evidence for mutations associated with crossovers. These findings provide important new insights into sources and structures of mitotic crossovers and functions of BLM helicase.     

DNA disentangling by type-2 topoisomerases

A type-2 topoisomerase cleaves a DNA strand, passes another through the break, and then rejoins the severed ends. Because it appears that this action is as likely to increase as to decrease entanglements, the question is: how are entanglements removed? We argue that type-2 topoisomerases have evolved to act at "hooked" juxtapositions of strands (where the strands are curved toward each other). This type of juxtaposition is a natural consequence of entangled long strands. Our model accounts for the observed preference for unlinking and unknotting of short DNA plasmids by type-2 topoisomerases and well explains experimental observations.     

Pericentromeric chromatin reorganisation follows the initiation of recombination and coincides with early events of synapsis in cereals

The reciprocal exchange of genetic information between homologous chromosomes during meiotic recombination is essential to secure balanced chromosome segregation and to promote genetic diversity. The chromosomal position and frequency of reciprocal genetic exchange shapes the efficiency of breeding programmes and influences crop improvement under a changing climate. In large genome cereals, such as wheat and barley, crossovers are consistently restricted to subtelomeric chromosomal regions, thus preventing favourable allele combinations being formed within a considerable proportion of the genome, including interstitial and pericentromeric chromatin. Understanding the key elements driving crossover designation is therefore essential to broaden the regions available for crossovers. Here, we followed early meiotic chromatin dynamism in cereals through the visualisation of a homologous barley chromosome arm pair stably transferred into the wheat genetic background. By capturing the dynamics of a single chromosome arm at the same time as detecting the undergoing events of meiotic recombination and synapsis, we showed that subtelomeric chromatin of homologues synchronously transitions to an open chromatin structure during recombination initiation. By contrast, pericentromeric and interstitial regions preserved their closed chromatin organisation and become unpackaged only later, concomitant with initiation of recombinatorial repair and the initial assembly of the synaptonemal complex. Our results raise the possibility that the closed pericentromeric chromatin structure in cereals may influence the fate decision during recombination initiation, as well as the spatial development of synapsis, and may also explain the suppression of crossover events in the proximity of the centromeres.     

The genetic architecture of genome‐wide recombination rate variation in allopolyploid wheat revealed by nested association mapping

Recombination affects the fate of alleles in populations by imposing constraints on the reshuffling of genetic information. Understanding the genetic basis of these constraints is critical for manipulating the recombination process to improve the resolution of genetic mapping, and reducing the negative effects of linkage drag and deleterious genetic load in breeding. Using sequence‐based genotyping of a wheat nested association mapping (NAM) population of 2,100 recombinant inbred lines created by crossing 29 diverse lines, we mapped QTL affecting the distribution and frequency of 102 000 crossovers (CO). Genome‐wide recombination rate variation was mostly defined by rare alleles with small effects together explaining up to 48.6% of variation. Most QTL were additive and showed predominantly trans‐acting effects. The QTL affecting the proximal COs also acted additively without increasing the frequency of distal COs. We showed that the regions with decreased recombination carry more single nucleotide polymorphisms (SNPs) with possible deleterious effects than the regions with a high recombination rate. Therefore, our study offers insights into the genetic basis of recombination rate variation in wheat and its effect on the distribution of deleterious SNPs across the genome. The identified trans‐acting additive QTL can be utilized to manipulate CO frequency and distribution in the large polyploid wheat genome opening the possibility to improve the efficiency of gene pyramiding and reducing the deleterious genetic load in the low‐recombining pericentromeric regions of chromosomes.