Gene conversion: mechanisms, evolution and human disease

High mutational rates of large-scale duplication and deletion in Daphnia pulex

Knowledge of the genome-wide rate and spectrum of mutations is necessary to understand the origin of disease and the genetic variation driving all evolutionary processes. Here, we provide a genome-wide analysis of the rate and spectrum of mutations obtained in two Daphnia pulex genotypes via separate mutation-accumulation (MA) experiments. Unlike most MA studies that utilize haploid, homozygous, or self-fertilizing lines, D. pulex can be propagated ameiotically while maintaining a naturally heterozygous, diploid genome, allowing the capture of the full spectrum of genomic changes that arise in a heterozygous state. While base-substitution mutation rates are similar to those in other multicellular eukaryotes (about 4 × 10(-9) per site per generation), we find that the rates of large-scale (>100 kb) de novo copy-number variants (CNVs) are significantly elevated relative to those seen in previous MA studies. The heterozygosity maintained in this experiment allowed for estimates of gene-conversion processes. While most of the conversion tract lengths we report are similar to those generated by meiotic processes, we also find larger tract lengths that are indicative of mitotic processes. Comparison of MA lines to natural isolates reveals that a majority of large-scale CNVs in natural populations are removed by purifying selection. The mutations observed here share similarities with disease-causing, complex, large-scale CNVs, thereby demonstrating that MA studies in D. pulex serve as a system for studying the processes leading to such alterations.

Comparative Genomic Analyses of the Human NPHP1 Locus Reveal Complex Genomic Architecture and Its Regional Evolution in Primates

Many loci in the human genome harbor complex genomic structures that can result in susceptibility to genomic rearrangements leading to various genomic disorders. Nephronophthisis 1 (NPHP1, MIM# 256100) is an autosomal recessive disorder that can be caused by defects of NPHP1; the gene maps within the human 2q13 region where low copy repeats (LCRs) are abundant. Loss of function of NPHP1 is responsible for approximately 85% of the NPHP1 cases—about 80% of such individuals carry a large recurrent homozygous NPHP1 deletion that occurs via nonallelic homologous recombination (NAHR) between two flanking directly oriented ~45 kb LCRs. Published data revealed a non-pathogenic inversion polymorphism involving the NPHP1 gene flanked by two inverted ~358 kb LCRs. Using optical mapping and array-comparative genomic hybridization, we identified three potential novel structural variant (SV) haplotypes at the NPHP1 locus that may protect a haploid genome from the NPHP1 deletion. Inter-species comparative genomic analyses among primate genomes revealed massive genomic changes during evolution. The aggregated data suggest that dynamic genomic rearrangements occurred historically within the NPHP1 locus and generated SV haplotypes observed in the human population today, which may confer differential susceptibility to genomic instability and the NPHP1 deletion within a personal genome. Our study documents diverse SV haplotypes at a complex LCR-laden human genomic region. Comparative analyses provide a model for how this complex region arose during primate evolution, and studies among humans suggest that intra-species polymorphism may potentially modulate an individual’s susceptibility to acquiring disease-associated alleles.

Genomic instability due to the intrinsic sequence architecture of the genome, such as low copy repeats (LCRs), is a major contributor to de novo mutations that can occur in the process of human genome evolution. LCRs can mediate genomic rearrangements associated with genomic disorders by acting as substrates for nonallelic homologous recombination. Juvenile-onset nephronophthisis 1 is the most frequent genetic cause of renal failure in children. An LCR-mediated, homozygous common recurrent deletion encompassing NPHP1 is found in the majority of affected subjects, while heterozygous deletion representing the nephronophthisis 1 recessive carrier state is frequently observed amongst world populations. Interestingly, the human NPHP1 locus is located proximal to the head-to-head fusion site of two ancestral chromosomes that occurred in the great apes, which resulted in a reduction of chromosome number from 48 in nonhuman primates to the current 46 in humans. In this study, we characterized and provided evidence for the diverse genomic architecture at the NPHP1 locus and potential structural variant haplotypes in the human population. Furthermore, our analyses of primate genomes shed light on the massive changes of genomic architecture at the human NPHP1 locus and delineated a model for the emergence of the LCRs during primate evolution.

Understanding the contrasting spatial haplotype patterns of malaria-protective β-globin polymorphisms

The malaria-protective β-globin polymorphisms, sickle-cell (β(S)) and β(0)-thalassaemia, are canonical examples of human adaptation to infectious disease. Occurring on distinct genetic backgrounds, they vary markedly in their patterns of linked genetic variation at the population level, suggesting different evolutionary histories. β(S) is associated with five classical restriction fragment length polymorphism haplotypes that exhibit remarkable specificity in their geographical distributions; by contrast, β(0)-thalassaemia mutations are found on haplotypes whose distributions overlap considerably. Here, we explore why these two polymorphisms display contrasting spatial haplotypic distributions, despite having malaria as a common selective pressure. We present a meta-population genetic model, incorporating individual-based processes, which tracks the evolution of β-globin polymorphisms on different haplotypic backgrounds. Our simulations reveal that, depending on the rate of mutation, a large population size and/or high population growth rate are required for both the β(S)- and the β(0)-thalassaemia-like patterns. However, whilst the β(S)-like pattern is more likely when population subdivision is high, migration low and long-distance migration absent, the opposite is true for β(0)-thalassaemia. Including gene conversion has little effect on the overall probability of each pattern; however, when inter-haplotype fitness variation exists, gene conversion is more likely to have contributed to the diversity of haplotypes actually present in the population. Our findings highlight how the contrasting spatial haplotype patterns exhibited by β(S) and β(0)-thalassaemia may provide important indications as to the evolution of these adaptive alleles and the demographic history of the populations in which they have evolved.

Gene conversion between cationic trypsinogen (PRSS1) and the pseudogene trypsinogen 6 (PRSS3P2) in patients with chronic pancreatitis

Mutations of the human cationic trypsinogen gene (PRSS1) are frequently found in association with hereditary pancreatitis. The most frequent variants p.N29I and p.R122H are recognized as disease-causing mutations. Three pseudogene paralogs in the human trypsinogen family, including trypsinogen 6 (PRSS3P2), carry sequence variations in exon 3 that mimic the p.R122H mutation. In routine genetic testing of patients with chronic pancreatitis, we identified in two unrelated individuals similar gene conversion events of 24-71 nucleotides length between exon 3 of the PRSS1 (acceptor) and PRSS3P2 (donor) genes. The converted allele resulted in three nonsynonymous alterations c.343T>A (p.S115T), c.347G>C (p.R116P), and c.365_366delinsAT (p.R122H). Functional analysis of the conversion triple mutant revealed markedly increased autoactivation resulting in high and sustained trypsin activity in the presence of chymotrypsin C. This activation phenotype was identical to that of the p.R122H mutant. In addition, cellular secretion of the triple mutant from transfected HEK 293T cells was increased about twofold and this effect was attributable to mutation p.R116P. Our observations confirm and extend the notion that recombination events between members of the trypsinogen family can generate high-risk PRSS1 alleles. The pathogenic phenotype of the novel conversion is explained by a unique combination of increased trypsinogen activation and secretion.

Concurrent nucleotide substitution mutations in the human genome are characterized by a significantly decreased transition/transversion ratio

There is accumulating evidence that the number of multiple-nucleotide substitutions (MNS) occurring in closely spaced sites in eukaryotic genomes is significantly higher than would be predicted from the random accumulation of independently generated single-nucleotide substitutions (SNS). Although this excess can in principle be accounted for by the concept of transient hypermutability, a general mutational signature of concurrent MNS mutations has not so far been evident. Employing a dataset (N = 449) of "concurrent" double MNS mutations causing human inherited disease, we have identified just such a mutational signature: concurrently generated double MNS mutations exhibit a >twofold lower transition/transversion ratio (termed RTs/Tv ) than independently generated de novo SNS mutations (<0.80 vs. 2.10; P = 2.69 × 10(-14) ). We replicated this novel finding through a similar analysis employing two double MNS variant datasets with differing abundances of concurrent events (150,521 variants with both substitutions on the same haplotypic lineage vs. 94,875 variants whose component substitutions were on different haplotypic lineages) plus 5,430,874 SNS variants, all being derived from the whole-genome sequencing of seven Chinese individuals. Evaluation of the newly observed mutational signature in diverse contexts provides solid support for the postulated role of translesion synthesis DNA polymerases in transient hypermutability.

Genotypic Diversity and Population Structure of Vibrio vulnificus Strains Isolated in Taiwan and Korea as Determined by Multilocus Sequence Typing

The genetic diversity and population structure of Vibrio vulnificus isolates from Korea and Taiwan were investigated using PCR-based assays targeting putative virulence-related genes and multilocus sequence typing (MLST). BOX-PCR genomic fingerprinting identified 52 unique genotypes in 84 environmental and clinical V. vulnificus isolates. The majority (> 50%) of strains had pathogenic genotypes for all loci tested; moreover, many environmental strains had pathogenic genotypes. Although significant (p < 0.05) inter-relationships among the genotypes were observed, the association between genotype and strain source (environmental or clinical) was not significant, indicating that genotypic characteristics alone are not sufficient to predict the isolation source or the virulence of a given V. vulnificus strain and vice versa. MLST revealed 23–35 allelic types per locus analyzed, resulting in a total of 44 unique sequence types (STs). Two major monophyletic groups (lineages A and B) corresponding to the two known lineages of V. vulnificus were observed; lineage A had six STs that were exclusively environmental, whereas lineage B had STs from both environmental and clinical sources. Pathogenic and nonpathogenic genotypes predominated in MLST lineages B and A, respectively. In addition, V. vulnificus was shown to be in linkage disequilibrium (p < 0.05), although two different recombination tests (PHI and Sawyer’s tests) detected significant evidence of recombination. Tajima’s D test also indicated that V. vulnificus might be comprised of recently sub-divided lineages. These results suggested that the two lineages revealed by MLST correspond to two distinct ecotypes of V. vulnificus.

Patterns of Gene Conversion in Duplicated Yeast Histones Suggest Strong Selection on a Coadapted Macromolecular Complex

We find evidence for interlocus gene conversion in five duplicated histone genes from six yeast species. The sequences of these duplicated genes, surviving from the ancient genome duplication, show phylogenetic patterns inconsistent with the well-resolved orthology relationships inferred from a likelihood model of gene loss after the genome duplication. Instead, these paralogous genes are more closely related to each other than any is to its nearest ortholog. In addition to simulations supporting gene conversion, we also present evidence for elevated rates of radical amino acid substitutions along the branches implicated in the conversion events. As these patterns are similar to those seen in ribosomal proteins that have undergone gene conversion, we speculate that in cases where duplicated genes code for proteins that are a part of tightly interacting complexes, selection may favor the fixation of gene conversion events in order to maintain high protein identities between duplicated copies.

A de novo mutation in CYP21A2 gene in a case of in vitro fertilization

Congenital adrenal hyperplasia, one of the most frequent autosome recessive disorders, is caused by defects in steroidogenic enzymes involved in the cortisol biosynthesis. Approximately 95% of the cases are caused by abnormal function of the 21-hydroxylase enzyme. This deficiency leads to androgen excess, consequently, to virilization and rapid somatic growth with accelerated skeletal maturation. Mutations in CYP21A2 are responsible for different forms of 21-hydroxylase deficiency. Mild impairment in the enzymatic activity causes the non-classic or late-onset congenital adrenal hyperplasia that is observed with a prevalence of 1 in 1000 subjects in different populations. The present paper describes a de novo mutation that occurred in the paternal meiosis. The child, who was conceived by in vitro fertilization, presented with precocious puberty and diagnosed with non-classical 21-hydroxylase deficiency. DNA sequencing showed the compound heterozygosis for a de novo CYP21A1P/A2 chimeric gene and the p.Val281Leu mutation inherited from her mother, who was heterozygous for the mutation. The chimeric gene showed pseudogene-derived sequence from 5′-end to intron 3 and CYP21A2 sequences from intron 3 to 3′-end of the gene. Sequencing analysis of the father did not show any mutation. The multiplex ligation-dependent probe amplification (MLPA) assay did not indicate loss of DNA discarding gene deletion but confirmed the chimeric gene. In addition, supernumerary copies of CYP21A1P were observed for both parents and for the affect child. Since paternity has been confirmed, those results suggest that a de novo large gene conversion in the paternal meiosis could have occurred by misalignment of alleles bearing different copy numbers of genes in CYP21 locus.

The type I interferons: Basic concepts and clinical relevance in immune-mediated inflammatory diseases

There is increasing scientific and clinical interest in elucidating the biology of type I Interferons, which began approximately 60 years ago with the concept of "viral interference", a property that reduces the ability of a virus to infect cells. Although our understanding of the multiple cellular and molecular functions of interferons has advanced significantly, much remains to be learned and type I Interferons remain an active and fascinating area of inquiry. In this review, we cover some general aspects of type I interferon genes, with emphasis on interferon-alpha, and various aspects of molecular mechanisms triggered by type I interferons and toll-like receptor signaling by the Janus activated kinase/signal transducer activation of transcription (JAK-STAT) pathway and interferon regulatory factor pathway. We will also describe the role of type I interferons in autoimmune and inflammatory diseases, and its potential use as therapeutic agent.

Co-option and evolution of non-olfactory proteinaceous pheromones in a terrestrial lungless salamander

Gene co-option is a major force in the evolution of novel biological functions. In plethodontid salamanders, males deliver proteinaceous courtship pheromones to the female olfactory system or transdermally to the bloodstream. Molecular studies identified three families of highly duplicated, rapidly evolving pheromones (PRF, PMF, and SPF). Analyses for Plethodon salamanders revealed pheromone mixtures of primarily PRF and PMF. The current study demonstrates that in Desmognathus ocoee--a plesiomorphic species with transdermal delivery--SPF is the major pheromone component representing >30% of total protein. Chromatographic profiles of D. ocoee pheromones were consistent from May through October. LC/MS-MS analysis suggested uniform SPF isoform expression between individual male D. ocoee. A gene ancestry for SPF with the Three-Finger Protein superfamily was supported by intron-exon boundaries, but not by the disulfide bonding pattern. Further analysis of the pheromone mixture revealed paralogs to peptide hormones that contained mutations in receptor binding regions, such that these novel molecules may alter female physiology by acting as hormone agonists/antagonists. Cumulatively, gene co-option, duplication, and neofunctionalization have permitted recruitment of additional gene families for pheromone activity. Such independent co-option events may be playing a key role in salamander speciation by altering male traits that influence reproductive success.

Suppression of Meiotic Recombination by CENP-B Homologs in Schizosaccharomyces pombe

Meiotic homologous recombination (HR) is not uniform across eukaryotic genomes, creating regions of HR hot- and coldspots. Previous study reveals that the Spo11 homolog Rec12 responsible for initiation of meiotic double-strand breaks in the fission yeast Schizosaccharomyces pombe is not targeted to Tf2 retrotransposons. However, whether Tf2s are HR coldspots is not known. Here, we show that the rates of HR across Tf2s are similar to a genome average but substantially increase in mutants deficient for the CENP-B homologs. Abp1, which is the most prominent of the CENP-B family members and acts as the primary determinant of HR suppression at Tf2s, is required to prevent gene conversion and maintain proper recombination exchange of homologous alleles flanking Tf2s. In addition, Abp1-mediated suppression of HR at Tf2s requires all three of its domains with distinct functions in transcriptional repression and higher-order genome organization. We demonstrate that HR suppression of Tf2s can be robustly maintained despite disruption to chromatin factors essential for transcriptional repression and nuclear organization of Tf2s. Intriguingly, we uncover a surprising cooperation between the histone methyltransferase Set1 responsible for histone H3 lysine 4 methylation and the nonhomologous end joining pathway in ensuring the suppression of HR at Tf2s. Our study identifies a molecular pathway involving functional cooperation between a transcription factor with epigenetic regulators and a DNA repair pathway to regulate meiotic recombination at interspersed repeats.

The role of fibrinogen-related proteins in the gastropod immune response

Fibrinogen-related proteins or FREPs constitute a large family of molecules, defined by the presence of a fibrinogen-related domain (FReD). These molecules are found in all animals and are diverse in both form and function. Here, we review the current understanding of gastropod FREPs, which are characterized by the presence of a fibrinogen domain connected to one or two immunoglobulin superfamily domains by way of a short interceding region. We present a historical perspective on the discovery of FREPs in gastropods followed by a summary of advances made in the nearly two decades of research focused on the characterization of FREPs in Biomphalaria glabrata (BgFREPs). Topics covered include BgFREP genomic architecture, predicted structure and known functions, structural comparisons between BgFREPs, and evidence of somatic diversification. Also examined are the expression patterns of BgFREPs during snail development and immunological challenges. Recent functional characterization of the role BgFREPs play in the defence response against digenean trematodes is also presented, as well as new data investigating the nucleotide-level genomic conservation of FREPs among Pulmonate gastropods. Finally, we identify areas in need of further research. These include confirming and identifying the specific binding targets of BgFREPs and elucidating how they later engage snail haemocytes to elicit an immunological response, precise mechanisms and importance of BgFREP diversification, characterizing the tissue expression patterns of BgFREPs, as well as addressing whether gastropod FREPs retain immunological importance in alternative snail-trematode associations or more broadly in snail-pathogen interactions.

Complex Multiple-Nucleotide Substitution Mutations Causing Human Inherited Disease Reveal Novel Insights into the Action of Translesion Synthesis DNA Polymerases

Translesion synthesis (TLS) DNA polymerases allow the bypass of unrepaired lesions during DNA replication. Based upon mutational signatures of a subtype of multiple-nucleotide substitution (MNS) mutations causing human inherited disease, we have recently postulated two properties of TLS DNA polymerases in DNA repair, namely, the generation of neo-microhomologies potentiating strand-misalignment, and additional microlesions within the templated inserts when recruited to stalled replication forks. To provide further support for this postulate, we analyzed the mutational signatures of a new and complex subtype of pathogenic MNS mutation. Several mutations containing long templated inserts (8-19 bp) that are highly informative with regard to their underlying mutational mechanisms, harbor imprints of TLS DNA polymerase action. Dissecting the mechanism underlying the generation of the 19-bp insert implicated repeated participation of TLS DNA polymerases in the conversion of a damaged base into a complex MNS lesion through a process of successive template switching and bypass repair.

The molecular evolution of spiggin nesting glue in sticklebacks

Gene duplication and subsequent divergence can lead to the evolution of new functions and lineage-specific traits. In sticklebacks, the successive duplication of a mucin gene (MUC19) into a tandemly arrayed, multigene family has enabled the production of copious amounts of 'spiggin', a secreted adhesive protein essential for nest construction. Here, we examine divergence between spiggin genes among three-spined sticklebacks (Gasterosteus aculeatus) from ancestral marine and derived freshwater populations, and propose underpinning gene duplication mechanisms. Sanger sequencing revealed substantial diversity among spiggin transcripts, including alternatively spliced variants and interchromosomal spiggin chimeric genes. Comparative analysis of the sequenced transcripts and all other spiggin genes in the public domain support the presence of three main spiggin lineages (spiggin A, spiggin B and spiggin C) with further subdivisions within spiggin B (B1, B2) and spiggin C (C1, C2). Spiggin A had diverged least from the ancestral MUC19, while the spiggin C duplicates had diversified most substantially. In silico translations of the spiggin gene open reading frames predicted that spiggins A and B are secreted as long mucin-like polymers, while spiggins C1 and C2 are secreted as short monomers, with putative antimicrobial properties. We propose that diversification of duplicated spiggin genes has facilitated local adaptation of spiggin to a range of aquatic habitats.

RAD mapping reveals an evolving, polymorphic and fuzzy boundary of a plant pseudoautosomal region

How loss of genetic exchanges (recombination) evolves between sex chromosomes is a long-standing question. Suppressed recombination may evolve when a sexually antagonistic (SA) polymorphism occurs in a partially sex-linked 'pseudoautosomal' region (or 'PAR'), maintaining allele frequency differences between the two sexes, and creating selection for closer linkage with the fully sex-linked region of the Y chromosome in XY systems, or the W in ZW sex chromosome systems. Most evidence consistent with the SA polymorphism hypothesis is currently indirect, and more studies of the genetics and population genetics of PAR genes are clearly needed. The sex chromosomes of the plant Silene latifolia are suitable for such studies, as they evolved recently and the loss of recombination could still be ongoing. Here, we used RAD sequencing to genetically map sequences in this plant, which has a large genome (c. 3 gigabases) and no available whole-genome sequence. We mapped 83 genes on the sex chromosomes, and comparative mapping in the related species S. vulgaris supports previous evidence for additions to an ancestral PAR and identified at least 12 PAR genes. We describe evidence that recombination rates have been reduced in meiosis of both sexes, and differences in recombination between S. latifolia families suggest ongoing recombination suppression. Large allele frequency differences between the sexes were found at several loci closely linked to the PAR boundary, and genes in different regions of the PAR showed striking sequence diversity patterns that help illuminate the evolution of the PAR.

A Novel Terminal-Repeat Retrotransposon in Miniature (TRIM) Is Massively Expressed in Echinococcus multilocularis Stem Cells

Taeniid cestodes (including the human parasites Echinococcus spp. and Taenia solium) have very few mobile genetic elements (MGEs) in their genome, despite lacking a canonical PIWI pathway. The MGEs of these parasites are virtually unexplored, and nothing is known about their expression and silencing. In this work, we report the discovery of a novel family of small nonautonomous long terminal repeat retrotransposons (also known as terminal-repeat retrotransposons in miniature, TRIMs) which we have named ta-TRIM (taeniid TRIM). ta-TRIMs are only the second family of TRIM elements discovered in animals, and are likely the result of convergent reductive evolution in different taxonomic groups. These elements originated at the base of the taeniid tree and have expanded during taeniid diversification, including after the divergence of closely related species such as Echinococcus multilocularis and Echinococcus granulosus. They are massively expressed in larval stages, from a small proportion of full-length copies and from isolated terminal repeats that show transcriptional read-through into downstream regions, generating novel noncoding RNAs and transcriptional fusions to coding genes. In E. multilocularis, ta-TRIMs are specifically expressed in the germinative cells (the somatic stem cells) during asexual reproduction of metacestode larvae. This would provide a developmental mechanism for insertion of ta-TRIMs into cells that will eventually generate the adult germ line. Future studies of active and inactive ta-TRIM elements could give the first clues on MGE silencing mechanisms in cestodes.

Interlocus gene conversion explains at least 2.7% of single nucleotide variants in human segmental duplications

Background

Interlocus gene conversion (IGC) is a recombination-based mechanism that results in the unidirectional transfer of short stretches of sequence between paralogous loci. Although IGC is a well-established mechanism of human disease, the extent to which this mutagenic process has shaped overall patterns of segregating variation in multi-copy regions of the human genome remains unknown. One expected manifestation of IGC in population genomic data is the presence of one-to-one paralogous SNPs that segregate identical alleles.

Results

Here, I use SNP genotype calls from the low-coverage phase 3 release of the 1000 Genomes Project to identify 15,790 parallel, shared SNPs in duplicated regions of the human genome. My approach for identifying these sites accounts for the potential redundancy of short read mapping in multi-copy genomic regions, thereby effectively eliminating false positive SNP calls arising from paralogous sequence variation. I demonstrate that independent mutation events to identical nucleotides at paralogous sites are not a significant source of shared polymorphisms in the human genome, consistent with the interpretation that these sites are the outcome of historical IGC events. These putative signals of IGC are enriched in genomic contexts previously associated with non-allelic homologous recombination, including clear signals in gene families that form tandem intra-chromosomal clusters.

Conclusions

Taken together, my analyses implicate IGC, not point mutation, as the mechanism generating at least 2.7 % of single nucleotide variants in duplicated regions of the human genome.

Electronic supplementary material

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

High-resolution skim genotyping by sequencing reveals the distribution of crossovers and gene conversions in Cicer arietinum and Brassica napus

We characterise the distribution of crossover and non-crossover recombination in Brassica napus and Cicer arietinum using a low-coverage genotyping by sequencing pipeline SkimGBS. The growth of next-generation DNA sequencing technologies has led to a rapid increase in sequence-based genotyping for applications including diversity assessment, genome structure validation and gene-trait association. We have established a skim-based genotyping by sequencing method for crop plants and applied this approach to genotype-segregating populations of Brassica napus and Cicer arietinum. Comparison of progeny genotypes with those of the parental individuals allowed the identification of crossover and non-crossover (gene conversion) events. Our results identify the positions of recombination events with high resolution, permitting the mapping and frequency assessment of recombination in segregating populations.

Trans-species polymorphism at antimicrobial innate immunity cathelicidin genes of Atlantic cod and related species

Natural selection, the most important force in evolution, comes in three forms. Negative purifying selection removes deleterious variation and maintains adaptations. Positive directional selection fixes beneficial variants, producing new adaptations. Balancing selection maintains variation in a population. Important mechanisms of balancing selection include heterozygote advantage, frequency-dependent advantage of rarity, and local and fluctuating episodic selection. A rare pathogen gains an advantage because host defenses are predominantly effective against prevalent types. Similarly, a rare immune variant gives its host an advantage because the prevalent pathogens cannot escape the host's apostatic defense. Due to the stochastic nature of evolution, neutral variation may accumulate on genealogical branches, but trans-species polymorphisms are rare under neutrality and are strong evidence for balancing selection. Balanced polymorphism maintains diversity at the major histocompatibility complex (MHC) in vertebrates. The Atlantic cod is missing genes for both MHC-II and CD4, vital parts of the adaptive immune system. Nevertheless, cod are healthy in their ecological niche, maintaining large populations that support major commercial fisheries. Innate immunity is of interest from an evolutionary perspective, particularly in taxa lacking adaptive immunity. Here, we analyze extensive amino acid and nucleotide polymorphisms of the cathelicidin gene family in Atlantic cod and closely related taxa. There are three major clusters, Cath1, Cath2, and Cath3, that we consider to be paralogous genes. There is extensive nucleotide and amino acid allelic variation between and within clusters. The major feature of the results is that the variation clusters by alleles and not by species in phylogenetic trees and discriminant analysis of principal components. Variation within the three groups shows trans-species polymorphism that is older than speciation and that is suggestive of balancing selection maintaining the variation. Using Bayesian and likelihood methods positive and negative selection is evident at sites in the conserved part of the genes and, to a larger extent, in the active part which also shows episodic diversifying selection, further supporting the argument for balancing selection.

Sequence variation within the KIV-2 copy number polymorphism of the human LPA gene in African, Asian, and European populations

Amazingly little sequence variation is reported for the kringle IV 2 copy number variation (KIV 2 CNV) in the human LPA gene. Apart from whole genome sequencing projects, this region has only been analyzed in some detail in samples of European populations. We have performed a systematic resequencing study of the exonic and flanking intron regions within the KIV 2 CNV in 90 alleles from Asian, European, and four different African populations. Alleles have been separated according to their CNV length by pulsed field gel electrophoresis prior to unbiased specific PCR amplification of the target regions. These amplicons covered all KIV 2 copies of an individual allele simultaneously. In addition, cloned amplicons from genomic DNA of an African individual were sequenced. Our data suggest that sequence variation in this genomic region may be higher than previously appreciated. Detection probability of variants appeared to depend on the KIV 2 copy number of the analyzed DNA and on the proportion of copies carrying the variant. Asians had a high frequency of so-called KIV 2 type B and type C (together 70% of alleles), which differ by three or two synonymous substitutions respectively from the reference type A. This is most likely explained by the strong bottleneck suggested to have occurred when modern humans migrated to East Asia. A higher frequency of variable sites was detected in the Africans. In particular, two previously unreported splice site variants were found. One was associated with non-detectable Lp(a). The other was observed at high population frequencies (10% to 40%). Like the KIV 2 type B and C variants, this latter variant was also found in a high proportion of KIV 2 repeats in the affected alleles and in alleles differing in copy numbers. Our findings may have implications for the interpretation of SNP analyses in other repetitive loci of the human genome.

Dissecting meiotic recombination based on tetrad analysis by single-microspore sequencing in maize

Meiotic recombination drives eukaryotic sexual reproduction and the generation of genome diversity. Tetrad analysis, which examines the four chromatids resulting from a single meiosis, is an ideal method to study the mechanisms of homologous recombination. Here we develop a method to isolate the four microspores from a single tetrad in maize for the purpose of whole-genome sequencing. A high-resolution recombination map reveals that crossovers are unevenly distributed across the genome and are more likely to occur in the genic than intergenic regions, especially common in the 5′- and 3′-end regions of annotated genes. The direct detection of genomic exchanges suggests that conversions likely occur in most crossover tracts. Negative crossover interference and weak chromatid interference are observed at the population level. Overall, our findings further our understanding of meiotic recombination with implications for both basic and applied research.

The crossovers and gene conversions that occur during meiotic recombination contribute to genome diversity in eukaryotes. Here Li et al. describe a method of isolating individual microspores for whole-genome sequencing, providing new insights into the generation of genome diversity through sexual reproduction.

Non-allelic gene conversion enables rapid evolutionary change at multiple regulatory sites encoded by transposable elements

Transposable elements (TEs) allow rewiring of regulatory networks, and the recent amplification of the ISX element dispersed 77 functional but suboptimal binding sites for the dosage compensation complex to a newly formed X chromosome in Drosophila. Here we identify two linked refining mutations within ISX that interact epistatically to increase binding affinity to the dosage compensation complex. Selection has increased the frequency of this derived haplotype in the population, which is fixed at 30% of ISX insertions and polymorphic among another 41%. Sharing of this haplotype indicates that high levels of gene conversion among ISX elements allow them to ‘crowd-source’ refining mutations, and a refining mutation that occurs at any single ISX element can spread in two dimensions: horizontally across insertion sites by non-allelic gene conversion, and vertically through the population by natural selection. These results describe a novel route by which fully functional regulatory elements can arise rapidly from TEs and implicate non-allelic gene conversion as having an important role in accelerating the evolutionary fine-tuning of regulatory networks.

DOI:http://dx.doi.org/10.7554/eLife.05899.001

Mutations change genes and provide the raw material for evolution. Genes are sections of DNA that contain the instructions for making proteins or other molecules, and so determine the physical characteristics of each organism. Genetic mutations that increase an organism's number of offspring and chances of survival are more likely to be passed on to future generations. Changes to when or where a gene is switched on (so-called regulatory mutations) can also provide fitness benefits and can therefore be selected for during evolution.

Transposable elements are sequences of DNA that are also called ‘jumping genes’ because they can make copies of themselves and these copies of the transposable element can move to other locations in the genome. Some transposable elements contain sequences that switch on nearby genes. If different copies of a transposable element that contains such a regulatory sequence insert themselves in more than one place, it can result in a network of genes that can all be controlled in the same way. The regulatory sequences contained within transposable elements are not always optimal, but they can be fine-tuned through evolution.

A fruit fly called Drosophila miranda has a transposable element called ISX that has, over time, placed up to 77 regulatory sequences around one of this species' sex chromosomes. Just as in humans, female flies are XX and males are XY; but having only one copy of the X chromosome means that male flies need to increase the expression of certain genes to produce a full-dose of the molecules made by the genes. This process is called dosage compensation and in 2013 the 77 ISX regulatory sequences on the fruit fly's X chromosome were shown to help recruit the molecular machinery that carries out dosage compensation to nearby genes, albeit inefficiently. Now Ellison and Bachtrog—who also conducted the 2013 study—report how these transposable elements have been fine-tuned to make them more effective for dosage compensation.

Ellison and Bachtrog uncovered two mutations that make the ISX transposable element better at recruiting the dosage compensation molecular machinery. ISX spread around different locations along the fly's X chromosome before these mutations arose; this means that initially none of the 77 insertions carried the two mutations, but now 30% of the 77 elements have the mutations in all flies, and 41% have them in only some flies.

The same mutations have spread between the different ISX elements because transposable elements with the mutations have been used to directly convert other ISX elements without them. These mutations have also become more common in the fruit fly population by being passed on to offspring and increasing their survival. These two routes have accelerated the fine-tuning of these transposable elements for use in gene regulation. This implies that regulatory sequences derived from transposable elements evolve in a way that is fundamentally different from those that arise by other means, as the direct conversion between these insertions allows fine-tuning mutations to spread more rapidly.

DOI:http://dx.doi.org/10.7554/eLife.05899.002

Causes and consequences of crossing-over evidenced via a high-resolution recombinational landscape of the honey bee

Background

Social hymenoptera, the honey bee (Apis mellifera) in particular, have ultra-high crossover rates and a large degree of intra-genomic variation in crossover rates. Aligned with haploid genomics of males, this makes them a potential model for examining the causes and consequences of crossing over. To address why social insects have such high crossing-over rates and the consequences of this, we constructed a high-resolution recombination atlas by sequencing 55 individuals from three colonies with an average marker density of 314 bp/marker.

Results

We find crossing over to be especially high in proximity to genes upregulated in worker brains, but see no evidence for a coupling with immune-related functioning. We detect only a low rate of non-crossover gene conversion, contrary to current evidence. This is in striking contrast to the ultrahigh crossing-over rate, almost double that previously estimated from lower resolution data. We robustly recover the predicted intragenomic correlations between crossing over and both population level diversity and GC content, which could be best explained as indirect and direct consequences of crossing over, respectively.

Conclusions

Our data are consistent with the view that diversification of worker behavior, but not immune function, is a driver of the high crossing-over rate in bees. While we see both high diversity and high GC content associated with high crossing-over rates, our estimate of the low non-crossover rate demonstrates that high non-crossover rates are not a necessary consequence of high recombination rates.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0566-0) contains supplementary material, which is available to authorized users.

A Host-Pathogen Interaction Reduced to First Principles: Antigenic Variation in T. brucei

Antigenic variation is a common microbial survival strategy, powered by diversity in expressed surface antigens across the pathogen population over the course of infection. Even so, among pathogens, African trypanosomes have the most comprehensive system of antigenic variation described. African trypanosomes (Trypanosoma brucei spp.) are unicellular parasites native to sub-Saharan Africa, and the causative agents of sleeping sickness in humans and of n'agana in livestock. They cycle between two habitats: a specific species of fly (Glossina spp. or, colloquially, the tsetse) and the bloodstream of their mammalian hosts, by assuming a succession of proliferative and quiescent developmental forms, which vary widely in cell architecture and function. Key to each of the developmental forms that arise during these transitions is the composition of the surface coat that covers the plasma membrane. The trypanosome surface coat is extremely dense, covered by millions of repeats of developmentally specified proteins: procyclin gene products cover the organism while it resides in the tsetse and metacyclic gene products cover it while in the fly salivary glands, ready to make the transition to the mammalian bloodstream. But by far the most interesting coat is the Variant Surface Glycoprotein (VSG) coat that covers the organism in its infectious form (during which it must survive free living in the mammalian bloodstream). This coat is highly antigenic and elicits robust VSG-specific antibodies that mediate efficient opsonization and complement mediated lysis of the parasites carrying the coat against which the response was made. Meanwhile, a small proportion of the parasite population switches coats, which stimulates a new antibody response to the prevalent (new) VSG species and this process repeats until immune system failure. The disease is fatal unless treated, and treatment at the later stages is extremely toxic. Because the organism is free living in the blood, the VSG:antibody surface represents the interface between pathogen and host, and defines the interaction of the parasite with the immune response. This interaction (cycles of VSG switching, antibody generation, and parasite deletion) results in stereotypical peaks and troughs of parasitemia that were first recognized more than 100 years ago. Essentially, the mechanism of antigenic variation in T. brucei results from a need, at the population level, to maintain an extensive repertoire, to evade the antibody response. In this chapter, we will examine what is currently known about the VSG repertoire, its depth, and the mechanisms that diversify it both at the molecular (DNA) and at the phenotypic (surface displayed) level, as well as how it could interact with antibodies raised specifically against it in the host.

Dosage Sensitivity of RPL9 and Concerted Evolution of Ribosomal Protein Genes in Plants

The ribosome in higher eukaryotes is a large macromolecular complex composed of four rRNAs and eighty different ribosomal proteins. In plants, each ribosomal protein is encoded by multiple genes. Duplicate genes within a family are often necessary to provide a threshold dose of a ribosomal protein but in some instances appear to have non-redundant functions. Here, we addressed whether divergent members of the RPL9 gene family are dosage sensitive or whether these genes have non-overlapping functions. The RPL9 family in Arabidopsis thaliana comprises two nearly identical members, RPL9B and RPL9C, and a more divergent member, RPL9D. Mutations in RPL9C and RPL9D genes lead to delayed growth early in development, and loss of both genes is embryo lethal, indicating that these are dosage-sensitive and redundant genes. Phylogenetic analysis of RPL9 as well as RPL4, RPL5, RPL27a, RPL36a, and RPS6 family genes in the Brassicaceae indicated that multicopy ribosomal protein genes have been largely retained following whole genome duplication. However, these gene families also show instances of tandem duplication, small scale deletion, and evidence of gene conversion. Furthermore, phylogenetic analysis of RPL9 genes in angiosperm species showed that genes within a species are more closely related to each other than to RPL9 genes in other species, suggesting ribosomal protein genes undergo convergent evolution. Our analysis indicates that ribosomal protein gene retention following whole genome duplication contributes to the number of genes in a family. However, small scale rearrangements influence copy number and likely drive concerted evolution of these dosage-sensitive genes.

Fibrinogen-Related Proteins (FREPs) in Mollusks

Anti-parasite responses of the snail Biomphalaria glabrata involve antigen-reactive plasma lectins termed fibrinogen-related proteins (FREPs) comprising a C-terminal fibrinogen (FBG) domain and one or two upstream immunoglobulin domains. FREPs are highly polymorphic; they derive from several gene families with multiple loci and alleles that are diversified by exon loss, alternative splicing, and random somatic mutation (gene conversion and point mutations). Individual B. glabrata snails have dynamically distinct FREP sequence repertoires. The immune relevance of B. glabrata FREPs is indicated by FREP binding to polymorphic antigens of (snail-specific) digenean parasites and altered resistance of B. glabrata to digeneans following RNAi knockdown of FREPs. The compatibility polymorphism hypothesis proposes that FREP mutation increases the range of germline-encoded immune recognition in B. glabrata to counter antigenically-varied parasites. Somatic mutation may result from sequence exchange among tandemly arranged FREP genes in the genome, and analysis of sequence variants also suggests involvement of cytidine deaminase-like activity or epigenetic regulation. Without current indications of selection or retention of effective sequence variants toward immunological memory, FREP diversification is thought to afford B. glabrata immunity that is anticipatory but not adaptive. More remains to be learned about this system; other mollusks elaborate diversified lectins consisting of single FBG domains, and bona fide FREPs were reported from additional gastropod species, but these may not be diversified. Future comparative immunological studies and gene discovery driven by next-generation sequencing will further clarify taxonomic distribution of FREP diversification and the underlying mutator mechanisms as a component of immune function in mollusks.

NextGen sequencing reveals short double crossovers contribute disproportionately to genetic diversity in Toxoplasma gondii

Background

Toxoplasma gondii is a widespread protozoan parasite of animals that causes zoonotic disease in humans. Three clonal variants predominate in North America and Europe, while South American strains are genetically diverse, and undergo more frequent recombination. All three northern clonal variants share a monomorphic version of chromosome Ia (ChrIa), which is also found in unrelated, but successful southern lineages. Although this pattern could reflect a selective advantage, it might also arise from non-Mendelian segregation during meiosis. To understand the inheritance of ChrIa, we performed a genetic cross between the northern clonal type 2 ME49 strain and a divergent southern type 10 strain called VAND, which harbors a divergent ChrIa.

Results

NextGen sequencing of haploid F1 progeny was used to generate a genetic map revealing a low level of conventional recombination, with an unexpectedly high frequency of short, double crossovers. Notably, both the monomorphic and divergent versions of ChrIa were isolated with equal frequency. As well, ChrIa showed no evidence of being a sex chromosome, of harboring an inversion, or distorting patterns of segregation. Although VAND was unable to self fertilize in the cat, it underwent successful out-crossing with ME49 and hybrid survival was strongly associated with inheritance of ChrIII from ME49 and ChrIb from VAND.

Conclusions

Our findings suggest that the successful spread of the monomorphic ChrIa in the wild has not been driven by meiotic drive or related processes, but rather is due to a fitness advantage. As well, the high frequency of short double crossovers is expected to greatly increase genetic diversity among progeny from genetic crosses, thereby providing an unexpected and likely important source of diversity.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1168) contains supplementary material, which is available to authorized users.

Generation of antigenic diversity in Plasmodium falciparum by structured rearrangement of Var genes during mitosis

The most polymorphic gene family in P. falciparum is the ∼60 var genes distributed across parasite chromosomes, both in the subtelomeres and in internal regions. They encode hypervariable surface proteins known as P. falciparum erythrocyte membrane protein 1 (PfEMP1) that are critical for pathogenesis and immune evasion in Plasmodium falciparum. How var gene sequence diversity is generated is not currently completely understood. To address this, we constructed large clone trees and performed whole genome sequence analysis to study the generation of novel var gene sequences in asexually replicating parasites. While single nucleotide polymorphisms (SNPs) were scattered across the genome, structural variants (deletions, duplications, translocations) were focused in and around var genes, with considerable variation in frequency between strains. Analysis of more than 100 recombination events involving var exon 1 revealed that the average nucleotide sequence identity of two recombining exons was only 63% (range: 52.7–72.4%) yet the crossovers were error-free and occurred in such a way that the resulting sequence was in frame and domain architecture was preserved. Var exon 1, which encodes the immunologically exposed part of the protein, recombined in up to 0.2% of infected erythrocytes in vitro per life cycle. The high rate of var exon 1 recombination indicates that millions of new antigenic structures could potentially be generated each day in a single infected individual. We propose a model whereby var gene sequence polymorphism is mainly generated during the asexual part of the life cycle.

Malaria kills >600,000 people each year, with most deaths caused by Plasmodium falciparum. A family of proteins known as P. falciparum erythrocyte membrane protein 1, PfEMP1, is expressed on the surface of infected erythrocytes and plays an important role in pathogenesis. Each P. falciparum genome contains approximately 60 highly polymorphic var genes encoding the PfEMP1 proteins, and monoallelic expression with periodic switching results in immune evasion. Var gene polymorphism is thus critical to this survival strategy. We investigated how var gene diversity is generated by performing an in vitro evolution experiment, tracking var gene mutation in ‘real-time’ with whole genome sequencing. We found that genome structural variation is focused in and around var genes. These genetic rearrangements created new ‘chimeric’ var gene sequences during the mitotic part of the life cycle, and were consistent with processes of mitotic non-allelic homologous recombination. The recombinant var genes were always in frame and with conserved overall var gene architecture, and the recombination rate implies that many millions of rearranged var gene sequences are produced every 48-hour life cycle within infected individuals. In conclusion, we provide a detailed description of how new var gene sequences are continuously generated in the parasite genome, helping to explain long-term parasite survival within infected human hosts.

Seven novel HLA alleles reflect different mechanisms involved in the evolution of HLA diversity: description of the new alleles and review of the literature

The human leukocyte antigen (HLA) loci are among the most polymorphic genes in the human genome. The diversity of these genes is thought to be generated by different mechanisms including point mutation, gene conversion and crossing-over. During routine HLA typing, we discovered seven novel HLA alleles which were probably generated by different evolutionary mechanisms. HLA-B*41:21, HLA-DQB1*02:10 and HLA-DQA1*01:12 likely emerged from the common alleles of their groups by point mutations, all of which caused non-synonymous amino acid substitutions. In contrast, a deletion of one nucleotide leading to a frame shift with subsequent generation of a stop codon is responsible for the appearance of a null allele, HLA-A*01:123N. Whereas HLA-B*35:231 and HLA-B*53:31 were probably products of intralocus gene conversion between HLA-B alleles, HLA-C*07:294 presumably evolved by interlocus gene conversion between an HLA-C and an HLA-B allele. Our analysis of these novel alleles illustrates the different mechanisms which may have contributed to the evolution of HLA polymorphism.

A new hemoglobin variant: Hb Meylan [β73(E17)Asp → Phe; HBB: c.220G>T; c.221A>T] with a double base mutation at the same codon

We report a new β-globin chain variant: Hb Meylan [β73(E17)Asp → Phe; HBB: c.220G>T; c.221A>T]. The new variant results from a double nucleotide mutation at the same codon. The possible molecular mechanisms are discussed.

Biased gene conversion skews allele frequencies in human populations, increasing the disease burden of recessive alleles

Gene conversion results in the nonreciprocal transfer of genetic information between two recombining sequences, and there is evidence that this process is biased toward G and C alleles. However, the strength of GC-biased gene conversion (gBGC) in human populations and its effects on hereditary disease have yet to be assessed on a genomic scale. Using high-coverage whole-genome sequences of African hunter-gatherers, agricultural populations, and primate outgroups, we quantified the effects of GC-biased gene conversion on population genomic data sets. We find that genetic distances (FST and population branch statistics) are modified by gBGC. In addition, the site frequency spectrum is left-shifted when ancestral alleles are favored by gBGC and right-shifted when derived alleles are favored by gBGC. Allele frequency shifts due to gBGC mimic the effects of natural selection. As expected, these effects are strongest in high-recombination regions of the human genome. By comparing the relative rates of fixation of unbiased and biased sites, the strength of gene conversion was estimated to be on the order of Nb ≈ 0.05 to 0.09. We also find that derived alleles favored by gBGC are much more likely to be homozygous than derived alleles at unbiased SNPs (+42.2% to 62.8%). This results in a curse of the converted, whereby gBGC causes substantial increases in hereditary disease risks. Taken together, our findings reveal that GC-biased gene conversion has important population genetic and public health implications.

Origins and functional diversification of salinity-responsive Na(+) , K(+) ATPase α1 paralogs in salmonids

The Salmoniform whole-genome duplication is hypothesized to have facilitated the evolution of anadromy, but little is known about the contribution of paralogs from this event to the physiological performance traits required for anadromy, such as salinity tolerance. Here, we determined when two candidate, salinity-responsive paralogs of the Na(+) , K(+) ATPase α subunit (α1a and α1b) evolved and studied their evolutionary trajectories and tissue-specific expression patterns. We found that these paralogs arose during a small-scale duplication event prior to the Salmoniform, but after the teleost, whole-genome duplication. The 'freshwater paralog' (α1a) is primarily expressed in the gills of Salmoniformes and an unduplicated freshwater sister species (Esox lucius) and experienced positive selection in the freshwater ancestor of Salmoniformes and Esociformes. Contrary to our predictions, the 'saltwater paralog' (α1b), which is more widely expressed than α1a, did not experience positive selection during the evolution of anadromy in the Coregoninae and Salmonine. To determine whether parallel mutations in Na(+) , K(+) ATPase α1 may contribute to salinity tolerance in other fishes, we studied independently evolved salinity-responsive Na(+) , K(+) ATPase α1 paralogs in Anabas testudineus and Oreochromis mossambicus. We found that a quarter of the mutations occurring between salmonid α1a and α1b in functionally important sites also evolved in parallel in at least one of these species. Together, these data argue that paralogs contributing to salinity tolerance evolved prior to the Salmoniform whole-genome duplication and that strong selection and/or functional constraints have led to parallel evolution in salinity-responsive Na(+) , K(+) ATPase α1 paralogs in fishes.

Gene conversions are frequent but not under positive selection in the Siglec gene families of primates

Siglecs are cell surface proteins that belong to the immunoglobulin superfamily and which bind sialic acids. They are composed of two groups, the conserved Siglecs and the CD33-related Siglecs. Previous studies have reported the occurrence of gene conversions between human CD33-related Siglecs and suggested that these conversions are adaptive because they increase the diversity of these immunoglobulin-related genes. Here, we analyze the Siglec genes of five primate species and show that gene conversions are not observed between conserved Siglec genes but that they are frequent between primate CD33-related Siglecs. The gene conversions between CD33-related Siglec genes only occur between similar genes and equally frequently in sialic acid binding and nonbinding domains. Furthermore, dN/dS ratio tests show that most of the Ig-like V-type 1 and the Ig-like C2-type 1 domains of Siglec genes evolve either neutrally or under purifying selection and that gene conversions were not responsible for the positively selected regions detected in the Ig-like V-type1 domain of the human SIGLEC7 and SIGLEC9 genes. Our results suggest that the frequent gene conversions between CD33-related Siglec genes are simply a consequence of the high degree of sequence similarity of these genes and that they are not adaptive.

Plant genetics. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome

Oilseed rape (Brassica napus L.) was formed ~7500 years ago by hybridization between B. rapa and B. oleracea, followed by chromosome doubling, a process known as allopolyploidy. Together with more ancient polyploidizations, this conferred an aggregate 72× genome multiplication since the origin of angiosperms and high gene content. We examined the B. napus genome and the consequences of its recent duplication. The constituent An and Cn subgenomes are engaged in subtle structural, functional, and epigenetic cross-talk, with abundant homeologous exchanges. Incipient gene loss and expression divergence have begun. Selection in B. napus oilseed types has accelerated the loss of glucosinolate genes, while preserving expansion of oil biosynthesis genes. These processes provide insights into allopolyploid evolution and its relationship with crop domestication and improvement.

Crossing-over in a hypervariable species preferentially occurs in regions of high local similarity

Recombination between double-stranded DNA molecules is a key genetic process which occurs in a wide variety of organisms. Usually, crossing-over (CO) occurs during meiosis between genotypes with 98.0–99.9% sequence identity, because within-population nucleotide diversity only rarely exceeds 2%. However, some species are hypervariable and it is unclear how CO can occur between genotypes with less than 90% sequence identity. Here, we study CO in Schizophyllum commune, a hypervariable cosmopolitan basidiomycete mushroom, a frequently encountered decayer of woody substrates. We crossed two haploid individuals, from the United States and from Russia, and obtained genome sequences for their 17 offspring. The average genetic distance between the parents was 14%, making it possible to study CO at very high resolution. We found reduced levels of linkage disequilibrium between loci flanking the CO sites indicating that they are mostly confined to hotspots of recombination. Furthermore, CO events preferentially occurred in regions under stronger negative selection, in particular within exons that showed reduced levels of nucleotide diversity. Apparently, in hypervariable species CO must avoid regions of higher divergence between the recombining genomes due to limitations imposed by the mismatch repair system, with regions under strong negative selection providing the opportunity for recombination. These patterns are opposite to those observed in a number of less variable species indicating that population genomics of hypervariable species may reveal novel biological phenomena.

Whole-genome duplication in teleost fishes and its evolutionary consequences

Whole-genome duplication (WGD) events have shaped the history of many evolutionary lineages. One such duplication has been implicated in the evolution of teleost fishes, by far the most species-rich vertebrate clade. After initial controversy, there is now solid evidence that such event took place in the common ancestor of all extant teleosts. It is termed teleost-specific (TS) WGD. After WGD, duplicate genes have different fates. The most likely outcome is non-functionalization of one duplicate gene due to the lack of selective constraint on preserving both. Mechanisms that act on preservation of duplicates are subfunctionalization (partitioning of ancestral gene functions on the duplicates), neofunctionalization (assigning a novel function to one of the duplicates) and dosage selection (preserving genes to maintain dosage balance between interconnected components). Since the frequency of these mechanisms is influenced by the genes' properties, there are over-retained classes of genes, such as highly expressed ones and genes involved in neural function. The consequences of the TS-WGD, especially its impact on the massive radiation of teleosts, have been matter of controversial debate. It is evident that gene duplications are crucial for generating complexity and that WGDs provide large amounts of raw material for evolutionary adaptation and innovation. However, it is less clear whether the TS-WGD is directly linked to the evolutionary success of teleosts and their radiation. Recent studies let us conclude that TS-WGD has been important in generating teleost complexity, but that more recent ecological adaptations only marginally related to TS-WGD might have even contributed more to diversification. It is likely, however, that TS-WGD provided teleosts with diversification potential that can become effective much later, such as during phases of environmental change.

[High gene conversion frequency between genes encoding 2-deoxyglucose-6-phosphate phosphatase in 3 Saccharomyces species]

Gene conversions are nonreciprocal sequence exchanges between genes. They are relatively common in Saccharomyces cerevisiae, but few studies have investigated the evolutionary fate of gene conversions or their functional impacts. Here, we analyze the evolution and impact of gene conversions between the two genes encoding 2-deoxyglucose-6-phosphate phosphatase in S. cerevisiae, Saccharomyces paradoxus and Saccharomyces mikatae. Our results demonstrate that the last half of these genes are subject to gene conversions among these three species. The greater similarity and the greater percentage of GC nucleotides in the converted regions, as well as the absence of long regions of adjacent common converted sites, suggest that these gene conversions are frequent and occur independently in all three species. The high frequency of these conversions probably result from the fact that they have little impact on the protein sequences encoded by these genes.

Coherent somatic mutation in autoimmune disease

BACKGROUND

Many aspects of autoimmune disease are not well understood, including the specificities of autoimmune targets, and patterns of co-morbidity and cross-heritability across diseases. Prior work has provided evidence that somatic mutation caused by gene conversion and deletion at segmentally duplicated loci is relevant to several diseases. Simple tandem repeat (STR) sequence is highly mutable, both somatically and in the germ-line, and somatic STR mutations are observed under inflammation.

RESULTS

Protein-coding genes spanning STRs having markers of mutability, including germ-line variability, high total length, repeat count and/or repeat similarity, are evaluated in the context of autoimmunity. For the initiation of autoimmune disease, antigens whose autoantibodies are the first observed in a disease, termed primary autoantigens, are informative. Three primary autoantigens, thyroid peroxidase (TPO), phogrin (PTPRN2) and filaggrin (FLG), include STRs that are among the eleven longest STRs spanned by protein-coding genes. This association of primary autoantigens with long STR sequence is highly significant (p<3.0x10(-7)). Long STRs occur within twenty genes that are associated with sixteen common autoimmune diseases and atherosclerosis. The repeat within the TTC34 gene is an outlier in terms of length and a link with systemic lupus erythematosus is proposed.

CONCLUSIONS

The results support the hypothesis that many autoimmune diseases are triggered by immune responses to proteins whose DNA sequence mutates somatically in a coherent, consistent fashion. Other autoimmune diseases may be caused by coherent somatic mutations in immune cells. The coherent somatic mutation hypothesis has the potential to be a comprehensive explanation for the initiation of many autoimmune diseases.

Recombination in the human Pseudoautosomal region PAR1

The pseudoautosomal region (PAR) is a short region of homology between the mammalian X and Y chromosomes, which has undergone rapid evolution. A crossover in the PAR is essential for the proper disjunction of X and Y chromosomes in male meiosis, and PAR deletion results in male sterility. This leads the human PAR with the obligatory crossover, PAR1, to having an exceptionally high male crossover rate, which is 17-fold higher than the genome-wide average. However, the mechanism by which this obligatory crossover occurs remains unknown, as does the fine-scale positioning of crossovers across this region. Recent research in mice has suggested that crossovers in PAR may be mediated independently of the protein PRDM9, which localises virtually all crossovers in the autosomes. To investigate recombination in this region, we construct the most fine-scale genetic map containing directly observed crossovers to date using African-American pedigrees. We leverage recombination rates inferred from the breakdown of linkage disequilibrium in human populations and investigate the signatures of DNA evolution due to recombination. Further, we identify direct PRDM9 binding sites using ChIP-seq in human cells. Using these independent lines of evidence, we show that, in contrast with mouse, PRDM9 does localise peaks of recombination in the human PAR1. We find that recombination is a far more rapid and intense driver of sequence evolution in PAR1 than it is on the autosomes. We also show that PAR1 hotspot activities differ significantly among human populations. Finally, we find evidence that PAR1 hotspot positions have changed between human and chimpanzee, with no evidence of sharing among the hottest hotspots. We anticipate that the genetic maps built and validated in this work will aid research on this vital and fascinating region of the genome.

Genome-wide linkage-disequilibrium profiles from single individuals

Although the analysis of linkage disequilibrium (LD) plays a central role in many areas of population genetics, the sampling variance of LD is known to be very large with high sensitivity to numbers of nucleotide sites and individuals sampled. Here we show that a genome-wide analysis of the distribution of heterozygous sites within a single diploid genome can yield highly informative patterns of LD as a function of physical distance. The proposed statistic, the correlation of zygosity, is closely related to the conventional population-level measure of LD, but is agnostic with respect to allele frequencies and hence likely less prone to outlier artifacts. Application of the method to several vertebrate species leads to the conclusion that >80% of recombination events are typically resolved by gene-conversion-like processes unaccompanied by crossovers, with the average lengths of conversion patches being on the order of one to several kilobases in length. Thus, contrary to common assumptions, the recombination rate between sites does not scale linearly with distance, often even up to distances of 100 kb. In addition, the amount of LD between sites separated by <200 bp is uniformly much greater than can be explained by the conventional neutral model, possibly because of the nonindependent origin of mutations within this spatial scale. These results raise questions about the application of conventional population-genetic interpretations to LD on short spatial scales and also about the use of spatial patterns of LD to infer demographic histories.

Is homologous recombination really an error-free process?

Homologous recombination (HR) is an evolutionarily conserved process that plays a pivotal role in the equilibrium between genetic stability and diversity. HR is commonly considered to be error-free, but several studies have shown that HR can be error-prone. Here, we discuss the actual accuracy of HR. First, we present the product of genetic exchanges (gene conversion, GC, and crossing over, CO) and the mechanisms of HR during double strand break repair and replication restart. We discuss the intrinsic capacities of HR to generate genome rearrangements by GC or CO, either during DSB repair or replication restart. During this process, abortive HR intermediates generate genetic instability and cell toxicity. In addition to genome rearrangements, HR also primes error-prone DNA synthesis and favors mutagenesis on single stranded DNA, a key DNA intermediate during the HR process. The fact that cells have developed several mechanisms protecting against HR excess emphasize its potential risks. Consistent with this duality, several pro-oncogenic situations have been consistently associated with either decreased or increased HR levels. Nevertheless, this versatility also has advantages that we outline here. We conclude that HR is a double-edged sword, which on one hand controls the equilibrium between genome stability and diversity but, on the other hand, can jeopardize the maintenance of genomic integrity. Therefore, whether non-homologous end joining (which, in contrast with HR, is not intrinsically mutagenic) or HR is the more mutagenic process is a question that should be re-evaluated. Both processes can be "Dr. Jekyll" in maintaining genome stability/variability and "Mr. Hyde" in jeopardizing genome integrity.

Extensive and biased intergenomic nonreciprocal DNA exchanges shaped a nascent polyploid genome, Gossypium (cotton)

Genome duplication is thought to be central to the evolution of morphological complexity, and some polyploids enjoy a variety of capabilities that transgress those of their diploid progenitors. Comparison of genomic sequences from several tetraploid (AtDt) Gossypium species and genotypes with putative diploid A- and D-genome progenitor species revealed that unidirectional DNA exchanges between homeologous chromosomes were the predominant mechanism responsible for allelic differences between the Gossypium tetraploids and their diploid progenitors. Homeologous gene conversion events (HeGCEs) gradually subsided, declining to rates similar to random mutation during radiation of the polyploid into multiple clades and species. Despite occurring in a common nucleus, preservation of HeGCE is asymmetric in the two tetraploid subgenomes. At-to-Dt conversion is far more abundant than the reciprocal, is enriched in heterochromatin, is highly correlated with GC content and transposon distribution, and may silence abundant A-genome-derived retrotransposons. Dt-to-At conversion is abundant in euchromatin and genes, frequently reversing losses of gene function. The long-standing observation that the nonspinnable-fibered D-genome contributes to the superior yield and quality of tetraploid cotton fibers may be explained by accelerated Dt to At conversion during cotton domestication and improvement, increasing dosage of alleles from the spinnable-fibered A-genome. HeGCE may provide an alternative to (rare) reciprocal DNA exchanges between chromosomes in heterochromatin, where genes have approximately five times greater abundance of Dt-to-At conversion than does adjacent intergenic DNA. Spanning exon-to-gene-sized regions, HeGCE is a natural noninvasive means of gene transfer with the precision of transformation, potentially important in genetic improvement of many crop plants.

How population growth affects linkage disequilibrium

The "LD curve" relates the linkage disequilibrium (LD) between pairs of nucleotide sites to the distance that separates them along the chromosome. The shape of this curve reflects natural selection, admixture between populations, and the history of population size. This article derives new results about the last of these effects. When a population expands in size, the LD curve grows steeper, and this effect is especially pronounced following a bottleneck in population size. When a population shrinks, the LD curve rises but remains relatively flat. As LD converges toward a new equilibrium, its time path may not be monotonic. Following an episode of growth, for example, it declines to a low value before rising toward the new equilibrium. These changes happen at different rates for different LD statistics. They are especially slow for estimates of [Formula: see text], which therefore allow inferences about ancient population history. For the human population of Europe, these results suggest a history of population growth.

Interplay of interlocus gene conversion and crossover in segmental duplications under a neutral scenario

Interlocus gene conversion is a major evolutionary force that drives the concerted evolution of duplicated genomic regions. Theoretical models successfully have addressed the effects of interlocus gene conversion and the importance of crossover in the evolutionary fate of gene families and duplications but have not considered complex recombination scenarios, such as the presence of hotspots. To study the interplay between interlocus gene conversion and crossover, we have developed a forward-time simulator that allows the exploration of a wide range of interlocus gene conversion rates under different crossover models. Using it, we have analyzed patterns of nucleotide variation and linkage disequilibrium within and between duplicate regions, focusing on a neutral scenario with constant population size and validating our results with the existing theoretical models. We show that the interaction of gene conversion and crossover is nontrivial and that the location of crossover junctions is a fundamental determinant of levels of variation and linkage disequilibrium in duplicated regions. We also show that if crossover activity between duplications is strong enough, recurrent interlocus gene conversion events can break linkage disequilibrium within duplicates. Given the complex nature of interlocus gene conversion and crossover, we provide a framework to explore their interplay to help increase knowledge on molecular evolution within segmental duplications under more complex scenarios, such as demographic changes or natural selection.

Repeated evolution of chimeric fusion genes in the β-globin gene family of laurasiatherian mammals

The evolutionary fate of chimeric fusion genes may be strongly influenced by their recombinational mode of origin and the nature of functional divergence between the parental genes. In the β-globin gene family of placental mammals, the two postnatally expressed δ- and β-globin genes (HBD and HBB, respectively) have a propensity for recombinational exchange via gene conversion and unequal crossing-over. In the latter case, there are good reasons to expect differences in retention rates for the reciprocal HBB/HBD and HBD/HBB fusion genes due to thalassemia pathologies associated with the HBD/HBB "Lepore" deletion mutant in humans. Here, we report a comparative genomic analysis of the mammalian β-globin gene cluster, which revealed that chimeric HBB/HBD fusion genes originated independently in four separate lineages of laurasiatherian mammals: Eulipotyphlans (shrews, moles, and hedgehogs), carnivores, microchiropteran bats, and cetaceans. In cases where an independently derived "anti-Lepore" duplication mutant has become fixed, the parental HBD and/or HBB genes have typically been inactivated or deleted, so that the newly created HBB/HBD fusion gene is primarily responsible for synthesizing the β-type subunits of adult and fetal hemoglobin (Hb). Contrary to conventional wisdom that the HBD gene is a vestigial relict that is typically inactivated or expressed at negligible levels, we show that HBD-like genes often encode a substantial fraction (20-100%) of β-chain Hbs in laurasiatherian taxa. Our results indicate that the ascendancy or resuscitation of genes with HBD-like coding sequence requires the secondary acquisition of HBB-like promoter sequence via unequal crossing-over or interparalog gene conversion.