Relative roles of mutation and recombination in generating allelic polymorphism at an MHC class II locus in Peromyscus maniculatus

The MHC class II loci encoding cell surface antigens exhibit extremely high allelic polymorphism. There is considerable uncertainty in the literature over the relative roles of recombination and de novo mutation in generating this diversity. We studied class II sequence diversity and allelic polymorphism in two populations of Peromyscus maniculatus, which are among the most widespread and abundant mammals of North America. We find that intragenic recombination (or gene conversion) has been the predominant mode for the generation of allelic polymorphism in this species, with the amount of population recombination per base pair exceeding mutation by at least an order of magnitude during the history of the sample. Despite this, patchwork motifs of sites with high linkage disequilibrium are observed. This does not appear to be consistent with the much larger amount of recombination versus mutation in the history of the sample, unless the recombination rate is highly non-uniform over the sequence or selection maintains certain sites in linkage disequilibrium. We conclude that selection is most likely to be responsible for preserving sequence motifs in the presence of abundant recombination.     

Length Variation in 18S rRNA Expansion Segment 43/e4 of <Emphasis Type="Italic">Daphnia obtusa</Emphasis>: Ancient or Recurring Polymorphism?

Expansion segments in ribosomal DNA (rDNA) can show length variation at the level of the individual, yet our understanding of the evolutionary forces shaping this variation is incomplete. Previous studies of expansion segment 43/e4 of the 18S rRNA gene in Daphnia obtusa have examined this variation in six individuals; however, it is not known if the variation documented at this locus is representative of variation across the species’ geographic range. Furthermore, it is unclear whether length variants found in multiple individuals share common ancestry, or were generated de novo through recombination. We quantified expansion segment length variant frequencies in 134 individual D. obtusa from 33 populations at 15 sites across the species range in the US, and used a phylogeographic approach to determine whether recombination continues to add to the standing crop of variation at this locus. We identified seven length variants across the sampling range, which spans almost 3000 km. Based on the phylogeographic distribution of length variants in the expansion segment, we conclude that they are shared ancient polymorphisms that have persisted despite the operation of molecular mechanisms that cause the concerted evolution of multigene families such as rDNA. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Direct Measure of the De Novo Mutation Rate in Autism and Schizophrenia Cohorts

The role of de novo mutations (DNMs) in common diseases remains largely unknown. Nonetheless, the rate of de novo deleterious mutations and the strength of selection against de novo mutations are critical to understanding the genetic architecture of a disease. Discovery of high-impact DNMs requires substantial high-resolution interrogation of partial or complete genomes of families via resequencing. We hypothesized that deleterious DNMs may play a role in cases of autism spectrum disorders (ASD) and schizophrenia (SCZ), two etiologically heterogeneous disorders with significantly reduced reproductive fitness. We present a direct measure of the de novo mutation rate (μ) and selective constraints from DNMs estimated from a deep resequencing data set generated from a large cohort of ASD and SCZ cases (n = 285) and population control individuals (n = 285) with available parental DNA. A survey of ∼430 Mb of DNA from 401 synapse-expressed genes across all cases and 25 Mb of DNA in controls found 28 candidate DNMs, 13 of which were cell line artifacts. Our calculated direct neutral mutation rate (1.36 × 10⁻⁸) is similar to previous indirect estimates, but we observed a significant excess of potentially deleterious DNMs in ASD and SCZ individuals. Our results emphasize the importance of DNMs as genetic mechanisms in ASD and SCZ and the limitations of using DNA from archived cell lines to identify functional variants.     

Mechanism,Prevalence, and More Severe Neuropathy Phenotype of the Charcot-Marie-Tooth Type 1A Triplication

Copy-number variations cause genomic disorders. Triplications, unlike deletions and duplications, are poorly understood because of challenges in molecular identification, the choice of a proper model system for study, and awareness of their phenotypic consequences. We investigated the genomic disorder Charcot-Marie-Tooth disease type 1A (CMT1A), a dominant peripheral neuropathy caused by a 1.4 Mb recurrent duplication occurring by nonallelic homologous recombination. We identified CMT1A triplications in families in which the duplication segregates. The triplications arose de novo from maternally transmitted duplications and caused a more severe distal symmetric polyneuropathy phenotype. The recombination that generated the triplication occurred between sister chromatids on the duplication-bearing chromosome and could accompany gene conversions with the homologous chromosome. Diagnostic testing for CMT1A (n = 20,661 individuals) identified 13% (n = 2,752 individuals) with duplication and 0.024% (n = 5 individuals) with segmental tetrasomy, suggesting that triplications emerge from duplications at a rate as high as ∼1:550, which is more frequent than the rate of de novo duplication. We propose that individuals with duplications are predisposed to acquiring triplications and that the population prevalence of triplication is underascertained.     

Enhanced Maternal Origin of the 22q11.2 Deletion in Velocardiofacial and DiGeorge Syndromes

Velocardiofacial and DiGeorge syndromes, also known as 22q11.2 deletion syndrome (22q11DS), are congenital-anomaly disorders caused by a de novo hemizygous 22q11.2 deletion mediated by meiotic nonallelic homologous recombination events between low-copy repeats, also known as segmental duplications. Although previous studies exist, each was of small size, and it remains to be determined whether there are parent-of-origin biases for the de novo 22q11.2 deletion. To address this question, we genotyped a total of 389 DNA samples from 22q11DS-affected families. A total of 219 (56%) individuals with 22q11DS had maternal origin and 170 (44%) had paternal origin of the de novo deletion, which represents a statistically significant bias for maternal origin (p = 0.0151). Combined with many smaller, previous studies, 465 (57%) individuals had maternal origin and 345 (43%) had paternal origin, amounting to a ratio of 1.35 or a 35% increase in maternal compared to paternal origin (p = 0.000028). Among 1,892 probands with the de novo 22q11.2 deletion, the average maternal age at time of conception was 29.5, and this is similar to data for the general population in individual countries. Of interest, the female recombination rate in the 22q11.2 region was about 1.6–1.7 times greater than that for males, suggesting that for this region in the genome, enhanced meiotic recombination rates, as well as other as-of-yet undefined 22q11.2-specific features, could be responsible for the observed excess in maternal origin.     

Detection of somaclonal variation during cocoa somatic embryogenesis characterised using cleaved amplified polymorphic sequence and the new freeware Artbio

The scarcity and stochastic nature of genetic mutations presents a significant challenge for scientists seeking to characterise de novo mutation frequency at specific loci. Such mutations can be particularly numerous during regeneration of plants from in vitro culture and can undermine the value of germplasm conservation efforts. We used cleaved amplified polymorphic sequence (CAPS) analysis to characterise new mutations amongst a clonal population of cocoa plants regenerated via a somatic embryogenesis protocol used previously for cocoa cryopreservation. Efficacy of the CAPS system for mutation detection was greatly improved after an ‘a priori’ in silico screen of reference target sequences for actual and potential restriction enzyme recognition sites using a new freely available software called Artbio. Artbio surveys known sequences for existing restriction enzyme recognition sites but also identifies all single nucleotide polymorphism (SNP) deviations from such motifs. Using this software, we performed an in silico screen of seven loci for restriction sites and their potential mutant SNP variants that were possible from 21 restriction enzymes. The four most informative locus-enzyme combinations were then used to survey the regenerant populations for de novo mutants. We characterised the pattern of point mutations and, using the outputs of Artbio, calculated the ratio of base substitution in 114 somatic embryo-derived cocoa regenerants originating from two explant genotypes. We found 49 polymorphisms, comprising 26.3% of the samples screened, with an inferred rate of 2.8 × 10⁻³ substitutions/screened base. This elevated rate is of a similar order of magnitude to previous reports of de novo microsatellite length mutations arising in the crop and suggests caution should be exercised when applying somatic embryogenesis for the conservation of plant germplasm.     

Structural variation in two human genomes mapped at single-nucleotide resolution by whole genome de novo assembly

Here we use whole-genome de novo assembly of second-generation sequencing reads to map structural variation (SV) in an Asian genome and an African genome. Our approach identifies small- and intermediate-size homozygous variants (1-50 kb) including insertions, deletions, inversions and their precise breakpoints, and in contrast to other methods, can resolve complex rearrangements. In total, we identified 277,243 SVs ranging in length from 1-23 kb. Validation using computational and experimental methods suggests that we achieve overall <6% false-positive rate and <10% false-negative rate in genomic regions that can be assembled, which outperforms other methods. Analysis of the SVs in the genomes of 106 individuals sequenced as part of the 1000 Genomes Project suggests that SVs account for a greater fraction of the diversity between individuals than do single-nucleotide polymorphisms (SNPs). These findings demonstrate that whole-genome de novo assembly is a feasible approach to deriving more comprehensive maps of genetic variation.     

Gene conversion rapidly generates major histocompatibility complex diversity in recently founded bird populations

Population bottlenecks can restrict variation at functional genes, reducing the ability of populations to adapt to new and changing environments. Understanding how populations generate adaptive genetic variation following bottlenecks is therefore central to evolutionary biology. Genes of the major histocompatibility complex (MHC) are ideal models for studying adaptive genetic variation due to their central role in pathogen recognition. While de novo MHC sequence variation is generated by point mutation, gene conversion can generate new haplotypes by transferring sections of DNA within and across duplicated MHC loci. However, the extent to which gene conversion generates new MHC haplotypes in wild populations is poorly understood. We developed a 454 sequencing protocol to screen MHC class I exon 3 variation across all 13 island populations of Berthelot's pipit (Anthus berthelotii). We reveal that just 11-15 MHC haplotypes were retained when the Berthelot's pipit dispersed across its island range in the North Atlantic ca. 75,000 years ago. Since then, at least 26 new haplotypes have been generated in situ across populations. We show that most of these haplotypes were generated by gene conversion across divergent lineages, and that the rate of gene conversion exceeded that of point mutation by an order of magnitude. Gene conversion resulted in significantly more changes at nucleotide sites directly involved with pathogen recognition, indicating selection for functional variants. We suggest that the creation of new variants by gene conversion is the predominant mechanism generating MHC variation in genetically depauperate populations, thus allowing them to respond to pathogenic challenges.     

Detection of new MHC mutations in mice by skin grafting, tumor transplantation and monoclonal antibodies: a comparison

Two mechanisms of major histocompatibility complex (MHC) mutations have been described in mice; gene conversion and homologous but unequal recombination. However, our knowledge of mutations in MHC is incomplete because studies have been limited almost exclusively to two haplotypes, H-2b and H-2d, while hundreds of haplotypes exist in nature; it has been biased by the use of only one procedure of screening for mutation, skin grafting. We used three procedures to screen for MHC mutations: (1) conventional techniques of skin grafting, (2) syngeneic tumor transplantation and (3) typing with monoclonal anti-MHC antibodies (mAbs) and complement. The faster technique of tumor transplantation detected mutants similar to those discovered by skin grafting technique. Screening with mAbs allowed us to detect both mutants that are capable of rejecting standard skin grafts and those that are silent in skin grafting tests, and which therefore resulted in a higher apparent mutation frequency. Two mutants of the H-2a haplotype were found that carry concomitant class I and class II antigenic alterations. Both MHC mutants silent in skin grafting tests and mutants carrying concomitant class I and class II alterations have never been studied before and are expected to reveal new mechanisms of generating MHC mutations. 1-Ethyl-1-nitrosourea (ENU) failed to induce de novo MHC mutations in male mice in our skin grafting series.     

De Novo Nonsense Mutations in KAT6A,a Lysine Acetyl-Transferase Gene,Cause a Syndrome Including Microcephaly and Global Developmental Delay

Chromatin remodeling through histone acetyltransferase (HAT) and histone deactylase (HDAC) enzymes affects fundamental cellular processes including the cell-cycle, cell differentiation, metabolism, and apoptosis. Nonsense mutations in genes that are involved in histone acetylation and deacetylation result in multiple congenital anomalies with most individuals displaying significant developmental delay, microcephaly and dysmorphism. Here, we report a syndrome caused by de novo heterozygous nonsense mutations in KAT6A (a.k.a., MOZ, MYST3) identified by clinical exome sequencing (CES) in four independent families. The same de novo nonsense mutation (c.3385C>T [p.Arg1129^∗]) was observed in three individuals, and the fourth individual had a nearby de novo nonsense mutation (c.3070C>T [p.Arg1024^∗]). Neither of these variants was present in 1,815 in-house exomes or in public databases. Common features among all four probands include primary microcephaly, global developmental delay including profound speech delay, and craniofacial dysmorphism, as well as more varied features such as feeding difficulties, cardiac defects, and ocular anomalies. We further demonstrate that KAT6A mutations result in dysregulation of H3K9 and H3K18 acetylation and altered P53 signaling. Through histone and non-histone acetylation, KAT6A affects multiple cellular processes and illustrates the complex role of acetylation in regulating development and disease.     

Mutational History of a Human Cell Lineage from Somatic to Induced Pluripotent Stem Cells

The accuracy of replicating the genetic code is fundamental. DNA repair mechanisms protect the fidelity of the genome ensuring a low error rate between generations. This sustains the similarity of individuals whilst providing a repertoire of variants for evolution. The mutation rate in the human genome has recently been measured to be 50–70 de novo single nucleotide variants (SNVs) between generations. During development mutations accumulate in somatic cells so that an organism is a mosaic. However, variation within a tissue and between tissues has not been analysed. By reprogramming somatic cells into induced pluripotent stem cells (iPSCs), their genomes and the associated mutational history are captured. By sequencing the genomes of polyclonal and monoclonal somatic cells and derived iPSCs we have determined the mutation rates and show how the patterns change from a somatic lineage in vivo through to iPSCs. Somatic cells have a mutation rate of 14 SNVs per cell per generation while iPSCs exhibited a ten-fold lower rate. Analyses of mutational signatures suggested that deamination of methylated cytosine may be the major mutagenic source in vivo, whilst oxidative DNA damage becomes dominant in vitro. Our results provide insights for better understanding of mutational processes and lineage relationships between human somatic cells. Furthermore it provides a foundation for interpretation of elevated mutation rates and patterns in cancer.     

Human SNP variability and mutation rate are higher in regions of high recombination

Understanding the co-variation of nucleotide diversity and local recombination rates is important both for the mapping of disease-associated loci and in understanding the causes of sequence evolution. It is known that single nucleotide polymorphisms (SNPs) around protein coding genes show higher diversity in regions of high recombination. Here, we find that this correlation holds for SNPs across the entire human genome, the great majority of which are not near exons or control elements. Contrasting with results from coding regions, we provide evidence that the higher nucleotide diversity in regions of high recombination is most likely due, at least in part, to a higher mutation rate. One possible explanation for this is that recombination is mutagenic.     

Evidence for effective suppression of recombination in the chromosome 17q21 segment spanning RNU2-BRCA1

Characterization of associations between polymorphic sites located throughout the approximately 200-400-kb variable-length region spanning RNU2-BRCA1 reveals nearly complete linkage disequilibrium. This segment spans the RNU2 array, which includes 6-30 tandem copies of the U2 snRNA gene, and an adjacent region containing NBR1, the LBRCA1 pseudogene, NBR2, and BRCA1 in a tandemly duplicated structure. A series of biallelic polymorphisms define two common haplotypes that do not vary significantly, in structure or frequency, between populations of primarily European (n=275) or Asian (n=34) ancestry. Lower-frequency variants occurring at distantly located sites within this region also show very strong associations. The rarer haplotype classes appear to be distinguished by mutational alteration and are not recombination products of the two major classes. The two major haplotypes also exhibit significantly different allele-length distributions for local simple tandem-repeat markers. The conservation of extensive distinct chromosomal haplotypes during a long period of human population expansion and divergence indicates that selective forces or specific chromosomal mechanisms result in effective recombination suppression. The extreme degree of long-range linkage disequilibrium at this locus may be exceeded only by that reported for the human MHC locus, where allele-specific functional interactions are believed to be significant. These findings have implications for the estimation of the time of origin of BRCA1 mutations having a founder effect, the interpretation of the significance of rare allelic variants, and the study of the origins of modern populations.     

De Novo Mutations in Synaptic Transmission Genes Including DNM1 Cause Epileptic Encephalopathies

Emerging evidence indicates that epileptic encephalopathies are genetically highly heterogeneous, underscoring the need for large cohorts of well-characterized individuals to further define the genetic landscape. Through a collaboration between two consortia (EuroEPINOMICS and Epi4K/EPGP), we analyzed exome-sequencing data of 356 trios with the “classical” epileptic encephalopathies, infantile spasms and Lennox Gastaut syndrome, including 264 trios previously analyzed by the Epi4K/EPGP consortium. In this expanded cohort, we find 429 de novo mutations, including de novo mutations in DNM1 in five individuals and de novo mutations in GABBR2, FASN, and RYR3 in two individuals each. Unlike previous studies, this cohort is sufficiently large to show a significant excess of de novo mutations in epileptic encephalopathy probands compared to the general population using a likelihood analysis (p = 8.2 × 10⁻⁴), supporting a prominent role for de novo mutations in epileptic encephalopathies. We bring statistical evidence that mutations in DNM1 cause epileptic encephalopathy, find suggestive evidence for a role of three additional genes, and show that at least 12% of analyzed individuals have an identifiable causal de novo mutation. Strikingly, 75% of mutations in these probands are predicted to disrupt a protein involved in regulating synaptic transmission, and there is a significant enrichment of de novo mutations in genes in this pathway in the entire cohort as well. These findings emphasize an important role for synaptic dysregulation in epileptic encephalopathies, above and beyond that caused by ion channel dysfunction.     

Evaluation of PRDM9 variation as a risk factor for recurrent genomic disorders and chromosomal non-disjunction

Recent studies have identified PRDM9, a zinc finger (ZF) protein, as a key regulator of meiotic recombination. As both recurrent genomic disorders and chromosomal non-disjunction are known to be associated with specific unusual patterns of recombination, we hypothesized a possible link between PRDM9 ZF variation and susceptibility to microdeletion syndromes and/or trisomy. We sequenced the PRDM9 ZF domain in 271 parents of patients with de novo microdeletions of known parental origin (velocardiofacial syndrome, the 17q21.31 microdeletion syndrome, Prader-Willi/Angelman syndrome and Williams-Beuren syndrome), and in 61 parents of individuals with a supernumerary X chromosome. We compared PRDM9 ZF genotype frequencies between parents in whose germ line the de novo rearrangement occurred and their spouses. We observed a significantly increased frequency (p?=?0.006) of PRDM9 variants in parents who transmitted de novo 7q11.23 deletions to their offspring. These data suggest that certain PRDM9 alleles may be associated with an increased susceptibility to recurrent 7q11.23 microdeletions that cause Williams-Beuren syndrome. However, as the majority of parents who transmitted a de novo microdeletion/supernumerary X chromosome to their offspring have the common AA genotype, we conclude that none of the rearrangements we have studied are dependent on specific non-A PRDM9 alleles.     

DVL3 Alleles Resulting in a ?1 Frameshift of the Last Exon Mediate Autosomal-Dominant Robinow Syndrome

Robinow syndrome is a rare congenital disorder characterized by mesomelic limb shortening, genital hypoplasia, and distinctive facial features. Recent reports have identified, in individuals with dominant Robinow syndrome, a specific type of variant characterized by being uniformly located in the penultimate exon of DVL1 and resulting in a −1 frameshift allele with a premature termination codon that escapes nonsense-mediated decay. Here, we studied a cohort of individuals who had been clinically diagnosed with Robinow syndrome but who had not received a molecular diagnosis from variant studies of DVL1, WNT5A, and ROR2. Because of the uniform location of frameshift variants in DVL1-mediated Robinow syndrome and the functional redundancy of DVL1, DVL2, and DVL3, we elected to pursue direct Sanger sequencing of the penultimate exon of DVL1 and its paralogs DVL2 and DVL3 to search for potential disease-associated variants. Remarkably, targeted sequencing identified five unrelated individuals harboring heterozygous, de novo frameshift variants in DVL3, including two splice acceptor mutations and three 1 bp deletions. Similar to the variants observed in DVL1-mediated Robinow syndrome, all variants in DVL3 result in a −1 frameshift, indicating that these highly specific alterations might be a common cause of dominant Robinow syndrome. Here, we review the current knowledge of these peculiar variant alleles in DVL1- and DVL3-mediated Robinow syndrome and further elucidate the phenotypic features present in subjects with DVL1 and DVL3 frameshift mutations.     

The B chromosome polymorphism of the grasshopper Eyprepocnemis plorans in North Africa: III. mutation rate of B chromosomes

B chromosome variation in nine Moroccan populations of the grasshopper Eyprepocnemis plorans was analysed for 3 consecutive years. In addition to B1, which was the predominant B chromosome in all nine populations, we found 15 other B variants, albeit at very low frequency. Eight variants were found in adults caught in the wild, four appeared in adults reared in the laboratory and seven were found in embryo progeny of controlled crosses between a 0B male and a B-carrying female. Some variants were found in more than one kind of material. At least the seven B variants that appeared in embryo progeny of females carrying a different B type arose de novo through mutation of the maternal B chromosome. The mutation rate of B chromosomes was 0.73%, on average, which explains the high variety of morphs and banding patterns found. The most frequent de novo mutations observed in these chromosomes were centromere misdivision with or without chromatid nondisjunction, which generates iso-B-chromosomes or telocentric Bs, respectively, as well as translocations with A and B chromosomes and deletions. But the whole variation observed, including that found in adult individuals, suggests that other mutations such as duplications, inversions and centric fusions do usually affect B chromosomes. Finally, B chromosome mutation rate was remarkably similar in both Moroccan and Spanish populations, which suggests that it might be dependent on B chromosome intrinsic factors.     

The influence of mutation, recombination, population history, and selection on patterns of genetic diversity in Neisseria meningitidis

Patterns of genetic diversity within populations of human pathogens, shaped by the ecology of host-microbe interactions, contain important information about the epidemiological history of infectious disease. Exploiting this information, however, requires a systematic approach that distinguishes the genetic signal generated by epidemiological processes from the effects of other forces, such as recombination, mutation, and population history. Here, a variety of quantitative techniques were employed to investigate multilocus sequence information from isolate collections of Neisseria meningitidis, a major cause of meningitis and septicemia world wide. This allowed quantitative evaluation of alternative explanations for the observed population structure. A coalescent-based approach was employed to estimate the rate of mutation, the rate of recombination, and the size distribution of recombination fragments from samples from disease-associated and carried meningococci obtained in the Czech Republic in 1993 and a global collection of disease-associated isolates collected globally from 1937 to 1996. The parameter estimates were used to reject a model in which genetic structure arose by chance in small populations, and analysis of molecular variation showed that geographically restricted gene flow was unlikely to be the cause of the genetic structure. The genetic differentiation between disease and carriage isolate collections indicated that, whereas certain genotypes were overrepresented among the disease-isolate collections (the "hyperinvasive" lineages), disease-associated and carried meningococci exhibited remarkably little differentiation at the level of individual nucleotide polymorphisms. In combination, these results indicated the repeated action of natural selection on meningococcal populations, possibly arising from the coevolutionary dynamic of host-pathogen interactions.     

TMEM43 mutations associated with arrhythmogenic right ventricular cardiomyopathy in non-Newfoundland populations

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a myocardial disease characterized by fibro-fatty replacement of right ventricular free wall myocardium and life-threatening ventricular arrhythmias. A missense mutation, c.1073C>T (p.S358L) in the transmembrane protein 43 (TMEM43) gene, has been genetically identified to cause ARVC type 5 in a founder population from Newfoundland. It is unclear whether this mutation occurs in other populations outside of this founder population or if other variants of TMEM43 are associated with ARVC disease. We sought to identify non-Newfoundland individuals with TMEM43 variants among patient samples sent for genetic assessment for possible ARVC. Of 195 unrelated individuals with suspected ARVC, mutation of desmosomal proteins was seen in 28 and the p.S358L TMEM43 mutation in six. We identified a de novo p.S358L mutation in a non-Newfoundland patient and five separate rare TMEM43 (four novel) sequence variants in non-Newfoundland patients, each occurring in an evolutionarily conserved amino acid. TMEM43 mutations occur outside of the founder population of the island of Newfoundland where it was originally described. TMEM43 sequencing should be incorporated into clinical genetic testing for ARVC patients.