Live hot, die young: transmission distortion in recombination hotspots

There is strong evidence that hotspots of meiotic recombination in humans are transient features of the genome. For example, hotspot locations are not shared between human and chimpanzee. Biased gene conversion in favor of alleles that locally disrupt hotspots is a possible explanation of the short lifespan of hotspots. We investigate the implications of such a bias on human hotspots and their evolution. Our results demonstrate that gene conversion bias is a sufficiently strong force to produce the observed lack of sharing of intense hotspots between species, although sharing may be much more common for weaker hotspots. We investigate models of how hotspots arise, and find that only models in which hotspot alleles do not initially experience drive are consistent with observations of rather hot hotspots in the human genome. Mutations acting against drive cannot successfully introduce such hotspots into the population, even if there is direct selection for higher recombination rates, such as to ensure correct segregation during meiosis. We explore the impact of hotspot alleles on patterns of haplotype variation, and show that such alleles mask their presence in population genetic data, making them difficult to detect.     

Cis- and trans-acting elements regulate the mouse Psmb9 meiotic recombination hotspot

In most eukaryotes, the prophase of the first meiotic division is characterized by a high level of homologous recombination between homologous chromosomes. Recombination events are not distributed evenly within the genome, but vary both locally and at large scale. Locally, most recombination events are clustered in short intervals (a few kilobases) called hotspots, separated by large intervening regions with no or very little recombination. Despite the importance of regulating both the frequency and the distribution of recombination events, the genetic factors controlling the activity of the recombination hotspots in mammals are still poorly understood. We previously characterized a recombination hotspot located close to the Psmb9 gene in the mouse major histocompatibility complex by sperm typing, demonstrating that it is a site of recombination initiation. With the goal of uncovering some of the genetic factors controlling the activity of this initiation site, we analyzed this hotspot in both male and female germ lines and compared the level of recombination in different hybrid mice. We show that a haplotype-specific element acts at distance and in trans to activate about 2,000-fold the recombination activity at Psmb9. Another haplotype-specific element acts in cis to repress initiation of recombination, and we propose this control to be due to polymorphisms located within the initiation zone. In addition, we describe subtle variations in the frequency and distribution of recombination events related to strain and sex differences. These findings show that most regulations observed act at the level of initiation and provide the first analysis of the control of the activity of a meiotic recombination hotspot in the mouse genome that reveals the interactions of elements located both in and outside the hotspot.     

Bayesian Inference of Shared Recombination Hotspots Between Humans and Chimpanzees

Recombination generates variation and facilitates evolution. Recombination (or lack thereof) also contributes to human genetic disease. Methods for mapping genes influencing complex genetic diseases via association rely on linkage disequilibrium (LD) in human populations, which is influenced by rates of recombination across the genome. Comparative population genomic analyses of recombination using related primate species can identify factors influencing rates of recombination in humans. Such studies can indicate how variable hotspots for recombination may be both among individuals (or populations) and over evolutionary timescales. Previous studies have suggested that locations of recombination hotspots are not conserved between humans and chimpanzees. We made use of the data sets from recent resequencing projects and applied a Bayesian method for identifying hotspots and estimating recombination rates. We also reanalyzed SNP data sets for regions with known hotspots in humans using samples from the human and chimpanzee. The Bayes factors (BF) of shared recombination hotspots between human and chimpanzee across regions were obtained. Based on the analysis of the aligned regions of human chromosome 21, locations where the two species show evidence of shared recombination hotspots (with high BFs) were identified. Interestingly, previous comparative studies of human and chimpanzee that focused on the known human recombination hotspots within the β-globin and HLA regions did not find overlapping of hotspots. Our results show high BFs of shared hotspots at locations within both regions, and the estimated locations of shared hotspots overlap with the locations of human recombination hotspots obtained from sperm-typing studies.     

DNA sequence-mediated, evolutionarily rapid redistribution of meiotic recombination hotspots

Hotspots regulate the position and frequency of Spo11 (Rec12)-initiated meiotic recombination, but paradoxically they are suicidal and are somehow resurrected elsewhere in the genome. After the DNA sequence-dependent activation of hotspots was discovered in fission yeast, nearly two decades elapsed before the key realizations that (A) DNA site-dependent regulation is broadly conserved and (B) individual eukaryotes have multiple different DNA sequence motifs that activate hotspots. From our perspective, such findings provide a conceptually straightforward solution to the hotspot paradox and can explain other, seemingly complex features of meiotic recombination. We describe how a small number of single-base-pair substitutions can generate hotspots de novo and dramatically alter their distribution in the genome. This model also shows how equilibrium rate kinetics could maintain the presence of hotspots over evolutionary timescales, without strong selective pressures invoked previously, and explains why hotspots localize preferentially to intergenic regions and introns. The model is robust enough to account for all hotspots of humans and chimpanzees repositioned since their divergence from the latest common ancestor.     

Structure activity relationship modelling of milk protein-derived peptides with dipeptidyl peptidase IV (DPP-IV) inhibitory activity

Quantitative structure activity type models were developed in an attempt to predict the key features of peptide sequences having dipeptidyl peptidase IV (DPP-IV) inhibitory activity. The models were then employed to help predict the potential of peptides, which are currently reported in the literature to be present in the intestinal tract of humans following milk/dairy product ingestion, to act as inhibitors of DPP-IV. Two models (z- and v-scale) for short (2–5 amino acid residues) bovine milk peptides, behaving as competitive inhibitors of DPP-IV, were developed. The z- and the v-scale models (p < 0.05, R² of 0.829 and 0.815, respectively) were then applied to 56 milk protein-derived peptides previously reported in the literature to be found in the intestinal tract of humans which possessed a structural feature of DPP-IV inhibitory peptides (P at the N₂ position). Ten of these peptides were synthetized and tested for their in vitro DPP-IV inhibitory properties. There was no agreement between the predicted and experimentally determined DPP-IV half maximal inhibitory concentrations (IC₅₀) for the competitive peptide inhibitors. However, the ranking for DPP-IV inhibitory potency of the competitive peptide inhibitors was conserved. Furthermore, potent in vitro DPP-IV inhibitory activity was observed with two peptides, LPVPQ (IC₅₀ = 43.8 ± 8.8 μM) and IPM (IC₅₀ = 69.5 ± 8.7 μM). Peptides present within the gastrointestinal tract of human may have promise for the development of natural DPP-IV inhibitors for the management of serum glucose.     

Capsid protein evolution and comparative phylogeny of novel porcine parvoviruses

In addition to the well known “classical” porcine parvovirus (PPV1; responsible for reproductive failure of susceptible sows) several new porcine parvoviruses have been recognized (PPV2, PPV3 and PPV4) in recent years. The genetic variation, characteristics and evolutionary factors shaping these novel PPVs were studied by comparing the complete capsid (cap) genes of PPVs from domestic pigs and wild boars. Using Bayesian coalescent methods we estimated the rate of nucleotide substitution for PPV2, PPV3 and PPV4 to be of the order of 3.86 × 10^?4–8.23 × 10^?4 subs site^?1 year^?1, similar to those commonly measured for RNA viruses, although this rate in case of PPV2 is probably influenced by frequent recombination events. Given such rapid evolutionary dynamics, it is likely that novel PPVs will continue to improve their capacity to spread among Suidae hosts worldwide. The mean time to the most recent common ancestor for the sampled genetic diversity of the newly discovered porcine parvoviruses was estimated. The results indicated that novel PPVs originated within approximately the last 70 years. Incongruent phylogenetic relationships of several strains suggested recombination events supported by several recombination-detecting methods and by split-decomposition phylogenetic networks. Analyses of the selective constraints acting on each codon suggest that some regions of PPV cap genes were under positive selection. This study showed that inter- and intraspecies recombination and diversifying selection pressures are prevalent across the cap genes of novel PPVs, and beside host switching and gene flow are important driving forces of their evolution and may be significant factors in the emergence of new viral variants.     

Size-dependent palindrome-induced intrachromosomal recombination in yeast

Palindromic and quasi-palindromic sequences are important DNA motifs found in various cis-acting genetic elements, but are also known to provoke different types of genetic alterations. The instability of such motifs is clearly size-related and depends on their potential to adopt secondary structures known as hairpins and cruciforms. Here we studied the influence of palindrome size on recombination between two directly repeated copies of the yeast CYC1 gene leading to the loss of the intervening sequence (“pop-out” recombination). We show that palindromes inserted either within one copy or between the two copies of the CYC1 gene become recombinogenic only when they attain a certain critical size and we estimate this critical size to be about 70 bp. With the longest palindrome used in this study (150 bp) we observed a more than 20-fold increase in the pop-out recombination. In the sae2/com1 mutant the palindrome-stimulated recombination was completely abolished. Suppression of palindrome recombinogenicity may be crucial for the maintenance of genetic stability in organisms containing a significant number of large palindromes in their genomes, like humans.     

A coalescent model of recombination hotspots

Recent experimental findings suggest that the assumption of a homogeneous recombination rate along the human genome is too naive. These findings point to block-structured recombination rates; certain regions (called hotspots) are more prone than other regions to recombination. In this report a coalescent model incorporating hotspot or block-structured recombination is developed and investigated analytically as well as by simulation. Our main results can be summarized as follows: (1) The expected number of recombination events is much lower in a model with pure hotspot recombination than in a model with pure homogeneous recombination, (2) hotspots give rise to large variation in recombination rates along the genome as well as in the number of historical recombination events, and (3) the size of a (nonrecombining) block in the hotspot model is likely to be overestimated grossly when estimated from SNP data. The results are discussed with reference to the current debate about block-structured recombination and, in addition, the results are compared to genome-wide variation in recombination rates. A number of new analytical results about the model are derived.     

An expanded genetic linkage map of an intervarietal Agaricus bisporus var. bisporus × A. bisporus var. burnettii hybrid based on AFLP,SSR and CAPS markers sheds light on the recombination behaviour of the species

A genetic linkage map for the edible basidiomycete Agaricus bisporus was constructed from 118 haploid homokaryons derived from an intervarietal A. bisporus var. bisporus × A. bisporus var. burnettii hybrid. Two hundred and thirty-one AFLP, 21 SSR, 68 CAPS markers together with the MAT, BSN, PPC1 loci and one allozyme locus (ADH) were evenly spread over 13 linkage groups corresponding to the chromosomes of A. bisporus. The map covers 1156 cM, with an average marker spacing of 3.9 cM and encompasses nearly the whole genome. The average number of crossovers per chromosome per individual is 0.86. Normal recombination over the entire genome occurs in the heterothallic variety, burnettii, contrary to the homothallic variety, bisporus, which showed adaptive genome-wide suppressed recombination. This first comprehensive genetic linkage map for A. bisporus provides foundations for quantitative trait analyses and breeding programme monitoring, as well as genome organisation studies.     

真核生物减数分裂重组热点的研究进展

真核生物减数分裂过程中基因组中某些区域会发生较其他区域高的重组频率,这些区域被称作减数分裂重组热点。该现象首先在酵母的研究中发现,重组区域因含有启动重组的特异位点,从而使基因组中呈现出重组不均匀分布的特征。重组热点还在真菌、玉米和人类等真核生物中发现。本文列举了不同真核生物体中具有代表性鉴别重组热点的方法,总结了目前减数分裂重组热点的研究现状,探讨了引起真核生物减数分裂交换活跃的因子和机制,并就当前存在的问题和今后发展的前景进行了讨论。     

Secondary structural modeling of the second internal transcribed spacer (ITS2) from Pfiesteria-like dinoflagellates (Dinophyceae)

Pfiesteria piscicida is a harmful bloom-forming alga that has received a great deal of attention due to its potential association with large fish kills and neurological problems in humans. Since the discovery of Pfiesteria, several other Pfiesteria-like dinoflagellates (PLDs) have also been identified. Genetic identification and phylogenetic relationships among the PLDs commonly utilize sequence data from the genes and spacers of the ribosomal DNA (rDNA) operon. Of these, the internal transcribed spacers (ITSs) have been previously shown to fold into secondary structures that are critical for proper ribosomal processing. In this study, we modeled the secondary structure of the second internal transcribed spacer (ITS2) from 16 PLDs (as well as an outgroup taxon) using phylogenetic comparative methods and minimum free energy. The secondary structural models predicted for these dinoflagellates consisted of four paired helices separated by five unpaired regions, consistent with those reported from many eukaryotes. All of the structures were highly stable (ΔG = ?66.1 to ?122.3 kcal·mol at 37 °C) and several structural characters were found to be conserved either across the PLDs or were specific to monophyletic subgroups, strengthening previously inferred phylogenetic relationships among taxa. Additionally, an 18 bp motif was identified in the PLDs whose position corresponds to a ribosomal processing site described from other eukaryotes. Potential applications of these ITS2 secondary structures include utility in strain and species identification, phylogenetic inference and serving as a tool for identifying and excluding rDNA pseudogenes when assessing biodiversity within the PLDs.     

High recombination rates and hotspots in a Plasmodium falciparum genetic cross

Background

The human malaria parasite Plasmodium falciparum survives pressures from the host immune system and antimalarial drugs by modifying its genome. Genetic recombination and nucleotide substitution are the two major mechanisms that the parasite employs to generate genome diversity. A better understanding of these mechanisms may provide important information for studying parasite evolution, immune evasion and drug resistance.

Results

Here, we used a high-density tiling array to estimate the genetic recombination rate among 32 progeny of a P. falciparum genetic cross (7G8 × GB4). We detected 638 recombination events and constructed a high-resolution genetic map. Comparing genetic and physical maps, we obtained an overall recombination rate of 9.6 kb per centimorgan and identified 54 candidate recombination hotspots. Similar to centromeres in other organisms, the sequences of P. falciparum centromeres are found in chromosome regions largely devoid of recombination activity. Motifs enriched in hotspots were also identified, including a 12-bp G/C-rich motif with 3-bp periodicity that may interact with a protein containing 11 predicted zinc finger arrays.

Conclusions

These results show that the P. falciparum genome has a high recombination rate, although it also follows the overall rule of meiosis in eukaryotes with an average of approximately one crossover per chromosome per meiosis. GC-rich repetitive motifs identified in the hotspot sequences may play a role in the high recombination rate observed. The lack of recombination activity in centromeric regions is consistent with the observations of reduced recombination near the centromeres of other organisms.     

Absence of the TAP2 human recombination hotspot in chimpanzees

Recent experiments using sperm typing have demonstrated that, in several regions of the human genome, recombination does not occur uniformly but instead is concentrated in “hotspots” of 1–2 kb. Moreover, the crossover asymmetry observed in a subset of these has led to the suggestion that hotspots may be short-lived on an evolutionary time scale. To test this possibility, we focused on a region known to contain a recombination hotspot in humans, TAP2, and asked whether chimpanzees, the closest living evolutionary relatives of humans, harbor a hotspot in a similar location. Specifically, we used a new statistical approach to estimate recombination rate variation from patterns of linkage disequilibrium in a sample of 24 western chimpanzees (Pan troglodytes verus). This method has been shown to produce reliable results on simulated data and on human data from the TAP2 region. Strikingly, however, it finds very little support for recombination rate variation at TAP2 in the western chimpanzee data. Moreover, simulations suggest that there should be stronger support if there were a hotspot similar to the one characterized in humans. Thus, it appears that the human TAP2 recombination hotspot is not shared by western chimpanzees. These findings demonstrate that fine-scale recombination rates can change between very closely related species and raise the possibility that rates differ among human populations, with important implications for linkage-disequilibrium based association studies.     

Contrasted Patterns of Crossover and Non-crossover at Arabidopsis thaliana Meiotic Recombination Hotspots

The vast majority of meiotic recombination events (crossovers (COs) and non-crossovers (NCOs)) cluster in narrow hotspots surrounded by large regions devoid of recombinational activity. Here, using a new molecular approach in plants, called “pollen-typing”, we detected and characterized hundreds of CO and NCO molecules in two different hotspot regions in Arabidopsis thaliana. This analysis revealed that COs are concentrated in regions of a few kilobases where their rates reach up to 50 times the genome average. The hotspots themselves tend to cluster in regions less than 8 kilobases in size with overlapping CO distribution. Non-crossover (NCO) events also occurred in the two hotspots but at very different levels (local CO/NCO ratios of 1/1 and 30/1) and their track lengths were quite small (a few hundred base pairs). We also showed that the ZMM protein MSH4 plays a role in CO formation and somewhat unexpectedly we also found that it is involved in the generation of NCOs but with a different level of effect. Finally, factors acting in cis and in trans appear to shape the rate and distribution of COs at meiotic recombination hotspots.     

Complete nucleotide sequence and organization of the mitochondrial genome of eri-silkworm,Samia cynthia ricini (Lepidoptera: Saturniidae)

Samia cynthia ricini is a commercial silk-producing insect that is now reared year-round in Korea, with the expectation of being utilized for diverse purposes. In this report, we present the complete mitochondrial genome (mitogenome) of S. c. ricini. The 15,384-bp long S. cynthia ricini mitogenome was amplified into 26 short fragments using three long overlapping fragments using primers designed from reported lepidopteran mitogenome sequences. The genome comprises 37 genes (13 protein-coding genes, two rRNA genes, and 22 tRNA genes), and one large non-coding region termed the A + T-rich region. The A/T content of the third codon position was 91.7%, which was 18.8% and 21.6% higher than those of first and second codon positions, respectively. The high A/T content in the genome is reflected in codon usage, accounting for 39.5% of A/T-composed codons (TTA, ATT, TTT, and ATA). Unlike a previous report on the start codon for the COI gene, the S. c. ricini COI gene commences with a typical ATT codon. A total of 221 bp of non-coding sequences are dispersed in 17 regions, ranging in size from 1 to 54 bp, which comprise 1.4% of the total genome. One of the non-coding sequence located between tRNA^Gln and ND2 (54 bp) has 77% sequence homology with the 5′-sequence of the neighboring ND2 gene, suggesting partial duplication of the sequence during evolution. The 361-bp long A + T-rich region contains an 18 bp-long poly-T stretch, ATAGA motif, ATTTA element, microsatellite-like A/T sequence, poly-A stretch and one tRNA-like sequence, as typically found in Lepidoptera including Bombycoidea.     

The impact of recombination on nucleotide substitutions in the human genome

Unraveling the evolutionary forces responsible for variations of neutral substitution patterns among taxa or along genomes is a major issue for detecting selection within sequences. Mammalian genomes show large-scale regional variations of GC-content (the isochores), but the substitution processes at the origin of this structure are poorly understood. We analyzed the pattern of neutral substitutions in 1 Gb of primate non-coding regions. We show that the GC-content toward which sequences are evolving is strongly negatively correlated to the distance to telomeres and positively correlated to the rate of crossovers (R²=47%). This demonstrates that recombination has a major impact on substitution patterns in human, driving the evolution of GC-content. The evolution of GC-content correlates much more strongly with male than with female crossover rate, which rules out selectionist models for the evolution of isochores. This effect of recombination is most probably a consequence of the neutral process of biased gene conversion (BGC) occurring within recombination hotspots. We show that the predictions of this model fit very well with the observed substitution patterns in the human genome. This model notably explains the positive correlation between substitution rate and recombination rate. Theoretical calculations indicate that variations in population size or density in recombination hotspots can have a very strong impact on the evolution of base composition. Furthermore, recombination hotspots can create strong substitution hotspots. This molecular drive affects both coding and non-coding regions. We therefore conclude that along with mutation, selection and drift, BGC is one of the major factors driving genome evolution. Our results also shed light on variations in the rate of crossover relative to non-crossover events, along chromosomes and according to sex, and also on the conservation of hotspot density between human and chimp.     

Complete mitochondrial genome of Coelomactra antiquata (Mollusca: Bivalvia): The first representative from the family Mactridae with novel gene order and unusual tandem repeats

The complete mitochondrial genome plays an important role in the accurate inference of phylogenetic relationships among metazoans. Mactridae, also known as trough shells or duck clams, is an important family of marine bivalve clams in the order Veneroida. Here we present the complete mitochondrial genome sequence of the Xishishe Coelomactra antiquata (Mollusca: Bivalvia), which is the first representative from the family Mactridae. The mitochondrial genome of C. antiquata is of 17,384 bp in length, and encodes 35 genes, including 12 protein-coding, 21 transfer RNA, and 2 ribosomal RNA genes. Compared with the typical gene content of animal mitochondrial genomes, atp8 and tRNAS₂ are missing. Gene order of the mitochondrial genome of C. antiquata is unique compared with others from Veneroida. In the mitochondrial genome of the C. antiquata, a total of 2189 bp of non-coding nucleotides are scattered among 26 non-coding regions. The largest non-coding region contains one section of tandem repeats (99 bp × 11), which is the second largest tandem repeats found in the mitochondrial genomes from Veneroida. The phylogenetic trees based on mitochondrial genomes support the monophyly of Veneridae and Lucinidae, and the relationship at the family level: ((Veneridae + Mactridae) + (Cardiidae + Solecurtidae)) + Lucinidae. The phylogenetic result is consistent with the morphological classification. Meanwhile, bootstrap values are very high (BP = 94–100), suggesting that the evolutionary relationship based on mitochondrial genomes is very reliable.     

Insights into plastome of Fagonia indica Burm.f. (Zygophyllaceae): organization,annotation and phylogeny

The enhanced understanding of chloroplast genomics would facilitate various biotechnology applications; however, the chloroplast (cp) genome / plastome characteristics of plants like Fagonia indica Burm.f. (family Zygophyllaceae), which have the capability to grow in extremely hot sand desert, have been rarely understood. The de novo genome sequence of F. indica using the Illumina high-throughput sequencing technology determined 128,379 bp long cp genome, encode 115 unique coding genes. The present study added the evidence of the loss of a copy of the IR in the cp genome of the taxa capable to grow in the hot sand desert. The maximum likelihood analysis revealed two distinct sub-clades i.e. Krameriaceae and Zygophyllaceae of the order Zygophyllales, nested within fabids.