Genome-wide patterns and properties of de novo mutations in humans

Applicability of the Mutation-Selection Balance Model to Population Genetics of Heterozygous Protein-Truncating Variants in Humans

The fate of alleles in the human population is believed to be highly affected by the stochastic force of genetic drift. Estimation of the strength of natural selection in humans generally necessitates a careful modeling of drift including complex effects of the population history and structure. Protein-truncating variants (PTVs) are expected to evolve under strong purifying selection and to have a relatively high per-gene mutation rate. Thus, it is appealing to model the population genetics of PTVs under a simple deterministic mutation-selection balance, as has been proposed earlier (Cassa et al. 2017). Here, we investigated the limits of this approximation using both computer simulations and data-driven approaches. Our simulations rely on a model of demographic history estimated from 33,370 individual exomes of the Non-Finnish European subset of the ExAC data set (Lek et al. 2016). Additionally, we compared the African and European subset of the ExAC study and analyzed de novo PTVs. We show that the mutation-selection balance model is applicable to the majority of human genes, but not to genes under the weakest selection.

Analysis of De Novo Mutations in Sporadic Cardiomyopathies Emphasizes Their Clinical Relevance and Points to Novel Candidate Genes

The vast majority of cardiomyopathies have an autosomal dominant inheritance; hence, genetic testing is typically offered to patients with a positive family history. A de novo mutation is a new germline mutation not inherited from either parent. The purpose of our study was to search for de novo mutations in patients with cardiomyopathy and no evidence of the disease in the family. Using next-generation sequencing, we analyzed cardiomyopathy genes in 12 probands. In 8 (66.7%), we found de novo variants in known cardiomyopathy genes (TTN, DSP, SCN5A, TNNC1, TPM1, CRYAB, MYH7). In the remaining probands, the analysis was extended to whole exome sequencing in a trio (proband and parents). We found de novo variants in genes that, so far, were not associated with any disease (TRIB3, SLC2A6), a possible disease-causing biallelic genotype (APOBEC gene family), and a de novo mosaic variant without strong evidence of pathogenicity (UNC45A). The high prevalence of de novo mutations emphasizes that genetic screening is also indicated in cases of sporadic cardiomyopathy. Moreover, we have identified novel cardiomyopathy candidate genes that are likely to affect immunological function and/or reaction to stress that could be especially relevant in patients with disease onset associated with infection/infestation.

Elucidation of de novo small insertion/deletion biology with parent-of-origin phasing

The mechanisms underlying de novo insertion/deletion (indel) genesis, such as polymerase slippage, have been hypothesized but not well characterized in the human genome. We implemented two methodological improvements, which were leveraged to dissect indel mutagenesis. We assigned de novo variants to parent-of-origin (i.e., phasing) with low-coverage long-read whole-genome sequencing, achieving better phasing compared to short-read sequencing (medians of 84% and 23%, respectively). We then wrote an application programming interface to classify indels into three subtypes according to sequence context. Across three cohorts with different phasing methods (N_trios  = 540, all cohorts), we observed that one de novo indel subtype, change in copy count (CCC), was significantly correlated with father's (p = 7.1 × 10^-4 ) but not mother's (p = .45) age at conception. We replicated this effect in three cohorts without de novo phasing (p_paternal  = 1.9 × 10^-9 , p_maternal  = .61; N_trios  = 3,391, all cohorts). Although this is consistent with polymerase slippage during spermatogenesis, the percentage of variance explained by paternal age was low, and we did not observe an association with replication timing. These results suggest that spermatogenesis-specific events have a minor role in CCC indel mutagenesis, one not observed for other indel subtypes nor for maternal age in general. These results have implications for indel modeling in evolution and disease.

Investigation of de novo variation in pediatric cardiomyopathy

Pediatric cardiomyopathies can be caused by variants in genes encoding the sarcomere and cytoskeleton in cardiomyocytes. Variants are typically inherited in an autosomal dominant manner with variable expressivity. De novo variants have been reported, however their overall frequency is largely unknown. We sought to determine the rate of de novo, pathogenic and likely pathogenic (P/LP) variants in children with a diagnosis of hypertrophic, dilated, or restrictive cardiomyopathy (HCM, DCM, or RCM), and to compare disease outcomes between individuals with and without a de novo variant. A retrospective record review identified 126 individuals with HCM (55%), DCM (37%), or RCM (8%) ≤18 years of age who had genetic testing. Overall, 50 (40%) had positive genetic testing and 18% of P/LP variants occurred de novo. The rate of de novo variation in those with RCM (80%) was higher than in those with HCM (9%) or DCM (20%). There was evidence of germline mosaicism in one family with RCM. Individuals with de novo variants were more likely than those without to have a history of arrhythmia (p = .049), sudden cardiac arrest (p = .024), hospitalization (p = .041), and cardiac transplantation (p = .030). The likelihood of de novo variation and impact on family risk and screening should be integrated into genetic counseling.

Impact of Mutation Rate and Selection at Linked Sites on DNA Variation across the Genomes of Humans and Other Homininae

DNA diversity varies across the genome of many species. Variation in diversity across a genome might arise from regional variation in the mutation rate, variation in the intensity and mode of natural selection, and regional variation in the recombination rate. We show that both noncoding and nonsynonymous diversity are positively correlated to a measure of the mutation rate and the recombination rate and negatively correlated to the density of conserved sequences in 50 kb windows across the genomes of humans and nonhuman homininae. Interestingly, we find that although noncoding diversity is equally affected by these three genomic variables, nonsynonymous diversity is mostly dominated by the density of conserved sequences. The positive correlation between diversity and our measure of the mutation rate seems to be largely a direct consequence of regions with higher mutation rates having more diversity. However, the positive correlation with recombination rate and the negative correlation with the density of conserved sequences suggest that selection at linked sites also affect levels of diversity. This is supported by the observation that the ratio of the number of nonsynonymous to noncoding polymorphisms is negatively correlated to a measure of the effective population size across the genome. We show these patterns persist even when we restrict our analysis to GC-conservative mutations, demonstrating that the patterns are not driven by GC biased gene conversion. In conclusion, our comparative analyses describe how recombination rate, gene density, and mutation rate interact to produce the patterns of DNA diversity that we observe along the hominine genomes.

A Return to the Origin of the EMGS: Rejuvenating the Quest for Human Germ Cell Mutagens and Determining the Risk to Future Generations

Fifty years ago, the Environmental Mutagen Society (now Environmental Mutagenesis and Genomics Society) was founded with a laser-focus on germ cell mutagenesis and the protection of "our most vital assets"-the sperm and egg genomes. Yet, five decades on, despite the fact that many agents have been demonstrated to induce inherited changes in the offspring of exposed laboratory rodents, there is no consensus on whether human germ cell mutagens exist. We argue that it is time to reevaluate the available data and conclude that we already have evidence for the existence of environmental exposures that impact human germ cells. What is missing are definite data to demonstrate a significant increase in de novo mutations in the offspring of exposed parents. We believe that with over two decades of research advancing knowledge and technologies in genomics, we are at the cusp of generating data to conclusively show that environmental exposures cause heritable de novo changes in the human offspring. We call on the research community to harness our technologies, synergize our efforts, and return to our Founders' original focus. The next 50 years must involve collaborative work between clinicians, epidemiologists, genetic toxicologists, genomics experts and bioinformaticians to precisely define how environmental exposures impact germ cell genomes. It is time for the research and regulatory communities to prepare to interpret the coming outpouring of data and develop a framework for managing, communicating and mitigating the risk of exposure to human germ cell mutagens. Environ. Mol. Mutagen. 61:42-54, 2020. © 2019 Her Majesty the Queen in Right of Canada.

Detection of Clinically Relevant Genetic Variants in Chinese Patients With Nanophthalmos by Trio-Based Whole-Genome Sequencing Study

Purpose

Nanophthalmos is a rare genetic disorder commonly characterized by a short axial length (AL) and severe hyperopia. Mutations that have been identified through Mendelian genetic analysis can only explain a fraction of nanophthalmic cases. We investigate the clinically relevant genetic variants in nanophthalmos by whole-genome sequencing (WGS), including de novo mutations (DNMs) and inherited mutations.

Methods

Clinically relevant genetic variants of 11 trios (11 nanophthalmic probands and their unaffected parents) from the Zhongshan Ophthalmic Center, China, were analyzed by WGS. We further screened three trios and 10 sporadic cases to identify the MYRF mutations.

Results

In two of 11 trios, without evidence of the presence of deleterious inherited autosomal variants, two DNMs of MYRF (c.789delC, p.S264fs and c.789dupC, p.S264fs) were identified in the probands. These loss-of-function DNMs were predicted to result in premature stop codons and protein structure damage in both probands. In addition, deleterious inherited genetic variants in PRSS56 and MFRP were found in eight probands of the other nine trios. Expanded screening found an additional MYRF DNM (c.1433G>C, p.R478P) in one trio and a stop-gain MYRF mutation (c.2956C>T, p.R986X) in one sporadic case, suggesting the recurrence of MYRF mutations in nanophthalmic patients.

Conclusions

This is the first trio-based WGS study for nanophthalmos, revealing the potential role of DNMs in MYRF and rare inherited genetic variants in PRSS56 and MFRP. The underlying mechanism of MYRF in the development of nanophthalmos needs to be further investigated.

De Novo Mutations Reflect Development and Aging of the Human Germline

Human germline de novo mutations (DNMs) are both a driver of evolution and an important cause of genetic diseases. In the past few years, whole-genome sequencing (WGS) of parent-offspring trios has facilitated the large-scale detection and study of human DNMs, which has led to exciting discoveries. The overarching theme of all of these studies is that the DNMs of an individual are a complex mixture of mutations that arise through different biological processes acting at different times during human development and life.

Inference of recombination maps from a single pair of genomes and its application to ancient samples

Understanding the causes and consequences of recombination landscape evolution is a fundamental goal in genetics that requires recombination maps from across the tree of life. Such maps can be obtained from population genomic datasets, but require large sample sizes. Alternative methods are therefore necessary to research organisms where such datasets cannot be generated easily, such as non-model or ancient species. Here we extend the sequentially Markovian coalescent model to jointly infer demography and the spatial variation in recombination rate. Using extensive simulations and sequence data from humans, fruit-flies and a fungal pathogen, we demonstrate that iSMC accurately infers recombination maps under a wide range of scenarios-remarkably, even from a single pair of unphased genomes. We exploit this possibility and reconstruct the recombination maps of ancient hominins. We report that the ancient and modern maps are correlated in a manner that reflects the established phylogeny of Neanderthals, Denisovans, and modern human populations.

Brain Development and Stochastic Processes During Prenatal and Early Life: You Can't Lose It if You've Never Had It; But It's Better to Have It and Lose It, Than Never to Have Had It at All

Brain development, although largely driven by genetic processes, also is influenced by environmental factors. However, there has been little discussion in the psychiatric literature on the role of stochastic, or chance, events that take place during neurodevelopment. Studies suggest that the brain capitalizes on and regulates the extent of stochastic processes during development. Furthermore, because neurodevelopment is influenced by environmental factors, there is emerging evidence that fostering those positive environmental factors during prenatal and early life could optimize neurodevelopment and provide greater resilience, including those potentially resulting from stochastic processes. Evidence for the role of environmental factors in optimizing early brain development is supported by work in large population-based studies of child development, randomized control trials in high-risk populations, and early-life adoption studies. The public health message is that creating an environment that fosters optimal brain development during prenatal and early life could prevent psychopathology and provide the developing brain the best chance against negative stochastic processes and potential stressors that are inevitable later in life.

Large, three-generation human families reveal post-zygotic mosaicism and variability in germline mutation accumulation

The number of de novo mutations (DNMs) found in an offspring's genome increases with both paternal and maternal age. But does the rate of mutation accumulation in human gametes differ across families? Using sequencing data from 33 large, three-generation CEPH families, we observed significant variability in parental age effects on DNM counts across families, ranging from 0.19 to 3.24 DNMs per year. Additionally, we found that ~3% of DNMs originated following primordial germ cell specification in a parent, and differed from non-mosaic germline DNMs in their mutational spectra. We also discovered that nearly 10% of candidate DNMs in the second generation were post-zygotic, and present in both somatic and germ cells; these gonosomal mutations occurred at equivalent frequencies on both parental haplotypes. Our results demonstrate that rates of germline mutation accumulation vary among families with similar ancestry, and confirm that post-zygotic mosaicism is a substantial source of human DNM.

Flanking heterozygosity influences the relative probability of different base substitutions in humans

Understanding when, where and which mutations are mostly likely to occur impacts many areas of evolutionary biology, from genetic diseases to phylogenetic reconstruction. Africans and non-African humans differ in the mutability of different triplet base combinations. Africans and non-Africans also differ in mutation rate, possibly because heterozygosity is mutagenic, such that diversity lost when humans expanded out of Africa also lowered the mutation rate. I show that these phenomena are linked: as flanking heterozygosity increases, some triplets become progressively more mutable while others become less so. Africans and non-African show near-identical patterns of dependence on heterozygosity. Thus, the striking differences in triplet mutation frequency between Africans and non-Africans, at least in part, seem to be an emergent property, driven by the way changes in heterozygosity 'out of Africa' have differentially impacted the mutability of different triplets. As heterozygosity decreased, the mutation spectrum outside Africa became enriched for triplet mutations that are favoured by low heterozygosity while those favoured by high heterozygosity became relatively rarer.

Signatures of replication timing, recombination, and sex in the spectrum of rare variants on the human X chromosome and autosomes

The sources of human germline mutations are poorly understood. Part of the difficulty is that mutations occur very rarely, and so direct pedigree-based approaches remain limited in the numbers that they can examine. To address this problem, we consider the spectrum of low-frequency variants in a dataset (Genome Aggregation Database, gnomAD) of 13,860 human X chromosomes and autosomes. X-autosome differences are reflective of germline sex differences and have been used extensively to learn about male versus female mutational processes; what is less appreciated is that they also reflect chromosome-level biochemical features that differ between the X and autosomes. We tease these components apart by comparing the mutation spectrum in multiple genomic compartments on the autosomes and between the X and autosomes. In so doing, we are able to ascribe specific mutation patterns to replication timing and recombination and to identify differences in the types of mutations that accrue in males and females. In particular, we identify C > G as a mutagenic signature of male meiotic double-strand breaks on the X, which may result from late repair. Our results show how biochemical processes of damage and repair in the germline interact with sex-specific life history traits to shape mutation patterns on both the X chromosome and autosomes.

Inherited and De Novo Genetic Risk for Autism Impacts Shared Networks

We performed a comprehensive assessment of rare inherited variation in autism spectrum disorder (ASD) by analyzing whole-genome sequences of 2,308 individuals from families with multiple affected children. We implicate 69 genes in ASD risk, including 24 passing genome-wide Bonferroni correction and 16 new ASD risk genes, most supported by rare inherited variants, a substantial extension of previous findings. Biological pathways enriched for genes harboring inherited variants represent cytoskeletal organization and ion transport, which are distinct from pathways implicated in previous studies. Nevertheless, the de novo and inherited genes contribute to a common protein-protein interaction network. We also identified structural variants (SVs) affecting non-coding regions, implicating recurrent deletions in the promoters of DLG2 and NR3C2. Loss of nr3c2 function in zebrafish disrupts sleep and social function, overlapping with human ASD-related phenotypes. These data support the utility of studying multiplex families in ASD and are available through the Hartwell Autism Research and Technology portal.

Mapping and Making Sense of Noncoding Mutations in the Genome

Whole-genome sequencing efforts of tumors and normal tissues have identified numerous genetic mutations, both somatic and germline, that do not overlap with coding genomic sequences. Attributing a functional role to these noncoding mutations and characterizing them using experimental methods has been more challenging compared with coding mutations. In this review, we provide a brief introduction to the world of noncoding mutations. We discuss recent progress in identifying noncoding mutations and the analytic and experimental approaches utilized to interpret their functional roles. We also highlight the potential mechanisms by which a noncoding mutation may exert its effect and discuss future challenges and opportunities.

Parental Age and Offspring Childhood Mental Health: A Multi-Cohort, Population-Based Investigation

To examine the contributions of maternal and paternal age on offspring externalizing and internalizing problems, this study analyzed problem behaviors at age 10-12 years from four Dutch population-based cohorts (N = 32,892) by a multiple informant design. Bayesian evidence synthesis was used to combine results across cohorts with 50% of the data analyzed for discovery and 50% for confirmation. There was evidence of a robust negative linear relation between parental age and externalizing problems as reported by parents. In teacher-reports, this relation was largely explained by parental socio-economic status. Parental age had limited to no association with internalizing problems. Thus, in this large population-based study, either a beneficial or no effect of advanced parenthood on child problem behavior was observed.

Somatic and Germline Mutation Periodicity Follow the Orientation of the DNA Minor Groove around Nucleosomes

Mutation rates along the genome are highly variable and influenced by several chromatin features. Here, we addressed how nucleosomes, the most pervasive chromatin structure in eukaryotes, affect the generation of mutations. We discovered that within nucleosomes, the somatic mutation rate across several tumor cohorts exhibits a strong 10 base pair (bp) periodicity. This periodic pattern tracks the alternation of the DNA minor groove facing toward and away from the histones. The strength and phase of the mutation rate periodicity are determined by the mutational processes active in tumors. We uncovered similar periodic patterns in the genetic variation among human and Arabidopsis populations, also detectable in their divergence from close species, indicating that the same principles underlie germline and somatic mutation rates. We propose that differential DNA damage and repair processes dependent on the minor groove orientation in nucleosome-bound DNA contribute to the 10-bp periodicity in AT/CG content in eukaryotic genomes.

Scales and mechanisms of somatic mutation rate variation across the human genome

Cancer genome sequencing has revealed that somatic mutation rates vary substantially across the human genome and at scales from megabase-sized domains to individual nucleotides. Here we review recent work that has both revealed the major mutation biases that operate across the genome and the molecular mechanisms that cause them. The default mutation rate landscape in mammalian genomes results in active genes having low mutation rates because of a combination of factors that increase DNA repair: early DNA replication, transcription, active chromatin modifications and accessible chromatin. Therefore, either an increase in the global mutation rate or a redistribution of mutations from inactive to active DNA can increase the rate at which consequential mutations are acquired in active genes. Several environmental carcinogens and intrinsic mechanisms operating in tumor cells likely cause cancer by this second mechanism: by specifically increasing the mutation rate in active regions of the genome.

Exploring background mutational processes to decipher cancer genetic heterogeneity

Much remains unknown about the progression and heterogeneity of mutational processes in different cancers and their diagnostic and clinical potential. A growing body of evidence supports mutation rate dependence on the local DNA sequence context for various types of mutations. We propose several tools for the analysis of cancer context-dependent mutations, which are implemented in an online computational framework MutaGene. The framework explores DNA context-dependent mutational patterns and underlying somatic cancer mutagenesis, analyzes mutational profiles of cancer samples, identifies the combinations of underlying mutagenic processes including those related to infidelity of DNA replication and repair machinery, and various other endogenous and exogenous mutagenic factors. As a result, the combination of mutagenic processes can be identified in any query sample with subsequent comparison to mutational profiles derived from malignant and benign samples. In addition, mutagen or cancer-specific mutational background models are applied to calculate expected DNA and protein site mutability to decouple relative contributions of mutagenesis and selection in carcinogenesis, thus elucidating the site-specific driving events in cancer. MutaGene is freely available at https://www.ncbi.nlm.nih.gov/projects/mutagene/.

Mutational Landscape of Spontaneous Base Substitutions and Small Indels in Experimental Caenorhabditis elegans Populations of Differing Size

Experimental investigations into the rates and fitness effects of spontaneous mutations are fundamental to our understanding of the evolutionary process. To gain insights into the molecular and fitness consequences of spontaneous mutations, we conducted a mutation accumulation (MA) experiment at varying population sizes in the nematode Caenorhabditis elegans, evolving 35 lines in parallel for 409 generations at three population sizes (N = 1, 10, and 100 individuals). Here, we focus on nuclear SNPs and small insertion/deletions (indels) under minimal influence of selection, as well as their accrual rates in larger populations under greater selection efficacy. The spontaneous rates of base substitutions and small indels are 1.84 (95% C.I. ± 0.14) × 10-9 substitutions and 6.84 (95% C.I. ± 0.97) × 10-10 changes/site/generation, respectively. Small indels exhibit a deletion bias with deletions exceeding insertions by threefold. Notably, there was no correlation between the frequency of base substitutions, nonsynonymous substitutions, or small indels with population size. These results contrast with our previous analysis of mitochondrial DNA mutations and nuclear copy-number changes in these MA lines, and suggest that nuclear base substitutions and small indels are under less stringent purifying selection compared to the former mutational classes. A transition bias was observed in exons as was a near universal base substitution bias toward A/T. Strongly context-dependent base substitutions, where 5'-Ts and 3'-As increase the frequency of A/T → T/A transversions, especially at the boundaries of A or T homopolymeric runs, manifest as higher mutation rates in (i) introns and intergenic regions relative to exons, (ii) chromosomal cores vs. arms and tips, and (iii) germline-expressed genes.

Estimation of allele-specific fitness effects across human protein-coding sequences and implications for disease

A central challenge in human genomics is to understand the cellular, evolutionary, and clinical significance of genetic variants. Here, we introduce a unified population-genetic and machine-learning model, called Linear Allele-Specific Selection InferencE (LASSIE), for estimating the fitness effects of all observed and potential single-nucleotide variants, based on polymorphism data and predictive genomic features. We applied LASSIE to 51 high-coverage genome sequences annotated with 33 genomic features and constructed a map of allele-specific selection coefficients across all protein-coding sequences in the human genome. This map is generally consistent with previous inferences of the bulk distribution of fitness effects but reveals pervasive weak negative selection against synonymous mutations. In addition, the estimated selection coefficients are highly predictive of inherited pathogenic variants and cancer driver mutations, outperforming state-of-the-art variant prioritization methods. By contrasting our estimated model with ultrahigh coverage ExAC exome-sequencing data, we identified 1118 genes under unusually strong negative selection, which tend to be exclusively expressed in the central nervous system or associated with autism spectrum disorder, as well as 773 genes under unusually weak selection, which tend to be associated with metabolism. This combination of classical population genetic theory with modern machine-learning and large-scale genomic data is a powerful paradigm for the study of both human evolution and disease.

What can be learned by scanning the genome for molecular convergence in wild populations?

Convergent evolution, where independent lineages evolve similar phenotypes in response to similar challenges, can provide valuable insight into how selection operates and the limitations it encounters. However, it has only recently become possible to explore how convergent evolution is reflected at the genomic level. The overlapping outlier approach (OOA), where genome scans of multiple independent lineages are used to find outliers that overlap and therefore identify convergently evolving loci, is becoming popular. Here, we present a quantitative analysis of 34 studies that used this approach across many sampling designs, taxa, and sampling intensities. We found that OOA studies with increased biological sampling power within replicates have increased likelihood of finding overlapping, "convergent" signals of adaptation between them. When identifying convergent loci as overlapping outliers, it is tempting to assume that any false-positive outliers derived from individual scans will fail to overlap across replicates, but this cannot be guaranteed. We highlight how population demographics and genomic context can contribute toward both true convergence and false positives in OOA studies. We finish with an exploration of emerging methods that couple genome scans with phenotype and environmental measures, leveraging added information from genome data to more directly test hypotheses of the likelihood of convergent evolution.

A pilot study of the implementation of pharmacogenomic pharmacist initiated pre-emptive testing in primary care

Despite the nationwide availability of pharmacogenomic (PGx) guidelines in electronic medication surveillance systems in The Netherlands, PGx guided prescribing is still uncommon in primary care. We set out to investigate the adoption of pharmacist initiated PGx testing in primary care. Community pharmacists were offered a free PGx test covering 40 variants in 8 genes to test patients receiving an incident prescription (IRx) of a selection of 10 drugs. Results of the PGx test along with predicted phenotypes and a therapeutic recommendation based on the Dutch Pharmacogenetics Working Group (DPWG) guidelines were transferred to the pharmacist and physician. Adoption was defined as the percentage of eligible patients that received genotyping. From November 2014-July 2016, 200 patients were included with an adoption of 18.0%. Of the included patients 57.5% received an IRx for atorvastatin, 14.5% started with simvastatin and 28.0% received an IRx for amitriptyline, (es)citalopram, nortriptyline, or venlafaxine. 90% of the patients carried at least one actionable PGx test result in the selected PGx-panel. In 31.0% of the incident prescriptions a combination between a drug with a known gene-drug interaction and an actionable genotype was present and a therapeutic recommendation was provided. The provided recommendations were accepted by the clinicians in 88.7% of the patients. Pharmacist initiated implementation of PGx in primary care is feasible, and the frequency of actionable gene-drug interactions for the selected drugs is high.

Paternal exposure to benzo(a)pyrene induces genome-wide mutations in mouse offspring

Understanding the effects of environmental exposures on germline mutation rates has been a decades-long pursuit in genetics. We used next-generation sequencing and comparative genomic hybridization arrays to investigate genome-wide mutations in the offspring of male mice exposed to benzo(a)pyrene (BaP), a common environmental pollutant. We demonstrate that offspring developing from sperm exposed during the mitotic or post-mitotic phases of spermatogenesis have significantly more de novo single nucleotide variants (1.8-fold; P < 0.01) than controls. Both phases of spermatogenesis are susceptible to the induction of heritable mutations, although mutations arising from post-fertilization events are more common after post-mitotic exposure. In addition, the mutation spectra in sperm and offspring of BaP-exposed males are consistent. Finally, we report a significant increase in transmitted copy number duplications (P = 0.001) in BaP-exposed sires. Our study demonstrates that germ cell mutagen exposures induce genome-wide mutations in the offspring that may be associated with adverse health outcomes.

Conditional Approximate Bayesian Computation: A New Approach for Across-Site Dependency in High-Dimensional Mutation-Selection Models

A key question in molecular evolutionary biology concerns the relative roles of mutation and selection in shaping genomic data. Moreover, features of mutation and selection are heterogeneous along the genome and over time. Mechanistic codon substitution models based on the mutation-selection framework are promising approaches to separating these effects. In practice, however, several complications arise, since accounting for such heterogeneities often implies handling models of high dimensionality (e.g., amino acid preferences), or leads to across-site dependence (e.g., CpG hypermutability), making the likelihood function intractable. Approximate Bayesian Computation (ABC) could address this latter issue. Here, we propose a new approach, named Conditional ABC (CABC), which combines the sampling efficiency of MCMC and the flexibility of ABC. To illustrate the potential of the CABC approach, we apply it to the study of mammalian CpG hypermutability based on a new mutation-level parameter implying dependence across adjacent sites, combined with site-specific purifying selection on amino-acids captured by a Dirichlet process. Our proof-of-concept of the CABC methodology opens new modeling perspectives. Our application of the method reveals a high level of heterogeneity of CpG hypermutability across loci and mild heterogeneity across taxonomic groups; and finally, we show that CpG hypermutability is an important evolutionary factor in rendering relative synonymous codon usage. All source code is available as a GitHub repository (https://github.com/Simonll/LikelihoodFreePhylogenetics.git).

Evolution of the Mutational Process under Relaxed Selection in Caenorhabditis elegans

The mutational process varies at many levels, from within genomes to among taxa. Many mechanisms have been linked to variation in mutation, but understanding of the evolution of the mutational process is rudimentary. Physiological condition is often implicated as a source of variation in microbial mutation rate and may contribute to mutation rate variation in multicellular organisms.Deleterious mutations are an ubiquitous source of variation in condition. We test the hypothesis that the mutational process depends on the underlying mutation load in two groups of Caenorhabditis elegans mutation accumulation (MA) lines that differ in their starting mutation loads. "First-order MA" (O1MA) lines maintained under minimal selection for ∼250 generations were divided into high-fitness and low-fitness groups and sets of "second-order MA" (O2MA) lines derived from each O1MA line were maintained for ∼150 additional generations. Genomes of 48 O2MA lines and their progenitors were sequenced. There is significant variation among O2MA lines in base-substitution rate (µbs), but no effect of initial fitness; the indel rate is greater in high-fitness O2MA lines. Overall, µbs is positively correlated with recombination and proximity to short tandem repeats and negatively correlated with 10 bp and 1 kb GC content. However, probability of mutation is sufficiently predicted by the three-nucleotide motif alone. Approximately 90% of the variance in standing nucleotide variation is explained by mutability. Total mutation rate increased in the O2MA lines, as predicted by the "drift barrier" model of mutation rate evolution. These data, combined with experimental estimates of fitness, suggest that epistasis is synergistic.

Germline Structural Variations Are Preferential Sites of DNA Replication Timing Plasticity during Development

The DNA replication timing program is modulated throughout development and is also one of the main factors influencing the distribution of mutation rates across the genome. However, the relationship between the mutagenic influence of replication timing and its developmental plasticity remains unexplored. Here, we studied the distribution of copy number variations (CNVs) and single nucleotide polymorphisms across the zebrafish genome in relation to changes in DNA replication timing during embryonic development in this model vertebrate species. We show that CNV sites exhibit strong replication timing plasticity during development, replicating significantly early during early development but significantly late during more advanced developmental stages. Reciprocally, genomic regions that changed their replication timing during development contained a higher proportion of CNVs than developmentally constant regions. Developmentally plastic CNV sites, in particular those that become delayed in their replication timing, were enriched for the clustered protocadherins, a set of genes important for neuronal development that have undergone extensive genetic and epigenetic diversification during zebrafish evolution. In contrast, single nucleotide polymorphism sites replicated consistently early throughout embryonic development, highlighting a unique aspect of the zebrafish genome. Our results uncover a hitherto unrecognized interface between development and evolution.

Overlooked roles of DNA damage and maternal age in generating human germline mutations

The textbook view that most germline mutations in mammals arise from replication errors is indirectly supported by the fact that there are both more mutations and more cell divisions in the male than in the female germline. When analyzing large de novo mutation datasets in humans, we find multiple lines of evidence that call that view into question. Notably, despite the drastic increase in the ratio of male to female germ cell divisions after the onset of spermatogenesis, even young fathers contribute three times more mutations than young mothers, and this ratio barely increases with parental age. This surprising finding points to a substantial contribution of damage-induced mutations. Indeed, C-to-G transversions and CpG transitions, which together constitute over one-fourth of all base substitution mutations, show genomic distributions and sex-specific age dependencies indicative of double-strand break repair and methylation-associated damage, respectively. Moreover, we find evidence that maternal age at conception influences the mutation rate both because of the accumulation of damage in oocytes and potentially through an influence on the number of postzygotic mutations in the embryo. These findings reveal underappreciated roles of DNA damage and maternal age in the genesis of human germline mutations.

Significant Strain Variation in the Mutation Spectra of Inbred Laboratory Mice

Mutation provides the ultimate source of all new alleles in populations, including variants that cause disease and fuel adaptation. Recent whole genome sequencing studies have uncovered variation in the mutation rate among individuals and differences in the relative frequency of specific nucleotide changes (the mutation spectrum) between populations. Although parental age is a major driver of differences in overall mutation rate among individuals, the causes of variation in the mutation spectrum remain less well understood. Here, I use high-quality whole genome sequences from 29 inbred laboratory mouse strains to explore the root causes of strain variation in the mutation spectrum. My analysis leverages the unique, mosaic patterns of genetic relatedness among inbred mouse strains to identify strain private variants residing on haplotypes shared between multiple strains due to their recent descent from a common ancestor. I show that these strain-private alleles are strongly enriched for recent de novo mutations and lack signals of widespread purifying selection, suggesting their faithful recapitulation of the spontaneous mutation landscape in single strains. The spectrum of strain-private variants varies significantly among inbred mouse strains reared under standardized laboratory conditions. This variation is not solely explained by strain differences in age at reproduction, raising the possibility that segregating genetic differences affect the constellation of new mutations that arise in a given strain. Collectively, these findings imply the action of remarkably precise nucleotide-specific genetic mechanisms for tuning the de novo mutation landscape in mammals and underscore the genetic complexity of mutation rate control.

SNP and indel frequencies at transcription start sites and at canonical and alternative translation initiation sites in the human genome

Single-nucleotide polymorphisms (SNPs) are the most common form of genetic variation in humans and drive phenotypic variation. Due to evolutionary conservation, SNPs and indels (insertion and deletions) are depleted in functionally important sequence elements. Recently, population-scale sequencing efforts such as the 1000 Genomes Project and the Genome of the Netherlands Project have catalogued large numbers of sequence variants. Here, we present a systematic analysis of the polymorphisms reported by these two projects in different coding and non-coding genomic elements of the human genome (intergenic regions, CpG islands, promoters, 5’ UTRs, coding exons, 3’ UTRs, introns, and intragenic regions). Furthermore, we were especially interested in the distribution of SNPs and indels in direct vicinity to the transcription start site (TSS) and translation start site (CSS). Thereby, we discovered an enrichment of dinucleotides CpG and CpA and an accumulation of SNPs at base position −1 relative to the TSS that involved primarily CpG and CpA dinucleotides. Genes having a CpG dinucleotide at TSS position -1 were enriched in the functional GO terms “Phosphoprotein”, “Alternative splicing”, and “Protein binding”. Focusing on the CSS, we compared SNP patterns in the flanking regions of canonical and alternative AUG and near-cognate start sites where we considered alternative starts previously identified by experimental ribosome profiling. We observed similar conservation patterns of canonical and alternative translation start sites, which underlines the importance of alternative translation mechanisms for cellular function.

De novo single-nucleotide and copy number variation in discordant monozygotic twins reveals disease-related genes

Recent studies have demonstrated genetic differences between monozygotic (MZ) twins. To test the hypothesis that early post-twinning mutational events associate with phenotypic discordance, we investigated a cohort of 13 twin pairs (n = 26) discordant for various clinical phenotypes using whole-exome sequencing and screened for copy number variation (CNV). We identified a de novo variant in PLCB1, a gene involved in the hydrolysis of lipid phosphorus in milk from dairy cows, associated with lactase non-persistence, and a variant in the mitochondrial complex I gene MT-ND5 associated with amyotrophic lateral sclerosis (ALS). We also found somatic variants in multiple genes (TMEM225B, KBTBD3, TUBGCP4, TFIP11) in another MZ twin pair discordant for ALS. Based on the assumption that discordance between twins could be explained by a common variant with variable penetrance or expressivity, we screened the twin samples for known pathogenic variants that are shared and identified a rare deletion overlapping ARHGAP11B, in the twin pair manifesting with either schizotypal personality disorder or schizophrenia. Parent-offspring trio analysis was implemented for two twin pairs to assess potential association of variants of parental origin with susceptibility to disease. We identified a de novo variant in RASD2 shared by 8-year-old male twins with a suspected diagnosis of autism spectrum disorder (ASD) manifesting as different traits. A de novo CNV duplication was also identified in these twins overlapping CD38, a gene previously implicated in ASD. In twins discordant for Tourette's syndrome, a paternally inherited stop loss variant was detected in AADAC, a known candidate gene for the disorder.

Excess of Rare Missense Variants in Hearing Loss Genes in Sporadic Meniere Disease

Meniere's disease (MD) is a clinical spectrum of rare disorders characterized by vertigo attacks, associated with sensorineural hearing loss (SNHL) and tinnitus involving low to medium frequencies. Although it shows familial aggregation with incomplete phenotypic forms and variable expressivity, most cases are considered sporadic. The aim of this study was to investigate the burden for rare variation in SNHL genes in patients with sporadic MD. We conducted a targeted-sequencing study including SNHL and familial MD genes in 890 MD patients to compare the frequency of rare variants in cases using three independent public datasets as controls. Patients with sporadic MD showed a significant enrichment of missense variants in SNHL genes that was not found in the controls. The list of genes includes GJB2, USH1G, SLC26A4, ESRRB, and CLDN14. A rare synonymous variant with unknown significance was found in the MARVELD2 gene in several unrelated patients with MD. There is a burden of rare variation in certain SNHL genes in sporadic MD. Furthermore, the interaction of common and rare variants in SNHL genes may have an additive effect on MD phenotype. This study will contribute to design a gene panel for the genetic diagnosis of MD.

Prenatal diagnosis for de novo mutations: Experience from a tertiary center over a 10-year period

Background

This study summarizes the results of prenatal diagnosis due to a history of de novo mutation in a previous pregnancy, in a tertiary center in Israel, over a 10‐year period.

Methods

We sorted all cases of de novo mutations from a pool of 2,260 pregnancies for which prenatal molecular diagnosis was applied, between the years 2008 and 2017. We identified 122 molecular prenatal diagnosis performed for de novo mutations, in 90 women.

Results

While the total number of yearly prenatal diagnoses stayed stable, a linear increase was detected in the number of cases for which the procedure was done due to a previous de novo mutation: from 3 cases in 2008 to 24 cases in 2017. The most common diseases were Rett syndrome (19), neurofibromatosis Type‐1 (12) and Tuberous sclerosis (5). Recurrence occurred in 3 of the 90 women (3.3%) and hotspot mutations were identified in two genes accounting for 11 cases. We did not find a difference in paternal age at first occurrence of the de novo mutation between the study group and the control group.

Conclusion

The large increase in the annual number of prenatal diagnoses performed due to a previous pregnancy with a de novo mutation reflects the growing understanding regarding the role of these mutations in the pathogenesis of genetic diseases.

Characterizing the Major Structural Variant Alleles of the Human Genome

In order to provide a comprehensive resource for human structural variants (SVs), we generated long-read sequence data and analyzed SVs for fifteen human genomes. We sequence resolved 99,604 insertions, deletions, and inversions including 2,238 (1.6 Mbp) that are shared among all discovery genomes with an additional 13,053 (6.9 Mbp) present in the majority, indicating minor alleles or errors in the reference. Genotyping in 440 additional genomes confirms the most common SVs in unique euchromatin are now sequence resolved. We report a ninefold SV bias toward the last 5 Mbp of human chromosomes with nearly 55% of all VNTRs (variable number of tandem repeats) mapping to this portion of the genome. We identify SVs affecting coding and noncoding regulatory loci improving annotation and interpretation of functional variation. These data provide the framework to construct a canonical human reference and a resource for developing advanced representations capable of capturing allelic diversity.

Characterizing mutagenic effects of recombination through a sequence-level genetic map

Genetic diversity arises from recombination and de novo mutation (DNM). Using a combination of microarray genotype and whole-genome sequence data on parent-child pairs, we identified 4,531,535 crossover recombinations and 200,435 DNMs. The resulting genetic map has a resolution of 682 base pairs. Crossovers exhibit a mutagenic effect, with overrepresentation of DNMs within 1 kilobase of crossovers in males and females. In females, a higher mutation rate is observed up to 40 kilobases from crossovers, particularly for complex crossovers, which increase with maternal age. We identified 35 loci associated with the recombination rate or the location of crossovers, demonstrating extensive genetic control of meiotic recombination, and our results highlight genes linked to the formation of the synaptonemal complex as determinants of crossovers.

Old Trade, New Tricks: Insights into the Spontaneous Mutation Process from the Partnering of Classical Mutation Accumulation Experiments with High-Throughput Genomic Approaches

Mutations spawn genetic variation which, in turn, fuels evolution. Hence, experimental investigations into the rate and fitness effects of spontaneous mutations are central to the study of evolution. Mutation accumulation (MA) experiments have served as a cornerstone for furthering our understanding of spontaneous mutations for four decades. In the pregenomic era, phenotypic measurements of fitness-related traits in MA lines were used to indirectly estimate key mutational parameters, such as the genomic mutation rate, new mutational variance per generation, and the average fitness effect of mutations. Rapidly emerging next-generating sequencing technology has supplanted this phenotype-dependent approach, enabling direct empirical estimates of the mutation rate and a more nuanced understanding of the relative contributions of different classes of mutations to the standing genetic variation. Whole-genome sequencing of MA lines bears immense potential to provide a unified account of the evolutionary process at multiple levels-the genetic basis of variation, and the evolutionary dynamics of mutations under the forces of selection and drift. In this review, we have attempted to synthesize key insights into the spontaneous mutation process that are rapidly emerging from the partnering of classical MA experiments with high-throughput sequencing, with particular emphasis on the spontaneous rates and molecular properties of different mutational classes in nuclear and mitochondrial genomes of diverse taxa, the contribution of mutations to the evolution of gene expression, and the rate and stability of transgenerational epigenetic modifications. Future advances in sequencing technologies will enable greater species representation to further refine our understanding of mutational parameters and their functional consequences.

Error-prone bypass of DNA lesions during lagging-strand replication is a common source of germline and cancer mutations

Studies in experimental systems have identified a multitude of mutational mechanisms including DNA replication infidelity and DNA damage followed by inefficient repair or replicative bypass. However, the relative contributions of these mechanisms to human germline mutation remain unknown. Here, we show that error-prone damage bypass on the lagging strand plays a major role in human mutagenesis. Transcription-coupled DNA repair removes lesions on the transcribed strand; lesions on the non-transcribed strand are preferentially converted into mutations. In human polymorphism we detect a striking similarity between mutation types predominant on the non-transcribed strand and on the strand lagging during replication. Moreover, damage-induced mutations in cancers accumulate asymmetrically with respect to the direction of replication, suggesting that DNA lesions are resolved asymmetrically. We experimentally demonstrate that replication delay greatly attenuates the mutagenic effect of ultraviolet irradiation, confirming that replication converts DNA damage into mutations. We estimate that at least 10% of human mutations arise due to DNA damage.

Role of gene body methylation in acclimatization and adaptation in a basal metazoan

Gene body methylation (GBM) has been hypothesized to modulate responses to environmental change, including transgenerational plasticity, but the evidence thus far has been lacking. Here we show that coral fragments reciprocally transplanted between two distant reefs respond predominantly by increase or decrease in genome-wide GBM disparity: The range of methylation levels between lowly and highly methylated genes becomes either wider or narrower. Remarkably, at a broad functional level this simple adjustment correlated very well with gene expression change, reflecting a shifting balance between expressions of environmentally responsive and housekeeping genes. In our experiment, corals in a lower-quality habitat up-regulated genes involved in environmental responses, while corals in a higher-quality habitat invested more in housekeeping genes. Transplanted fragments showing closer GBM match to local corals attained higher fitness characteristics, which supports GBM's role in acclimatization. Fixed differences in GBM between populations did not align with plastic GBM changes and were mostly observed in genes with elevated F _ST, which suggests that they arose predominantly through genetic divergence. However, we cannot completely rule out transgenerational inheritance of acquired GBM states.

Mendelian Randomization Studies Promise to Shorten the Journey to FDA Approval

There has been a dearth of new drugs approved for cardiovascular disorders. The cost is prohibitive, averaging to $2.5 billion, and requiring 12.5 years. This is in large part due to the high failure rate, with only 5% approval by the Food and Drug Administration. Despite preclinical studies showing potential safety and efficacy, most fail when they go to clinical trials phase I to III. One cause for failure is the drug target, often discovered to be a biomarker rather than causative for the disease. Mendelian randomization (MR) studies would determine whether the drug target is causative and could save millions of dollars and time, and prevent unnecessary exposure to adverse drug effects. This was demonstrated in 3 clinical trials that were negative with 2 drugs, veraspladib and darapladib. MR studies during the trials showed the targets of secretory and lipoprotein-associated phospholipids A2 are not causative for coronary artery disease and predicted negative results. The requirement for MR studies is a genetic risk variant with altered function, randomized at conception that remains fixed throughout one’s lifetime. It is not confounded by dietary, lifestyle, or socioeconomic factors. It is more sensitive than randomized controlled trials because exposure to the risk factor is fixed for a lifetime. MR studies showed plasma high-density lipoprotein cholesterol is not a causative target of coronary artery disease, and neither is uric acid, C-reactive protein, and others. MR studies are highly sensitive in determining whether drug targets are causative, and are relatively easy, inexpensive, and not time consuming. It is recommended that drug targets undergo MR studies before proceeding to randomized controlled trials.

Species Concepts: What about Ticks?

Since ancient times, philosophers and taxonomists have tried to classify forms of life. This is what taxonomy is about: the science of identifying, naming, classifying, and describing organisms. In this article I address the issue of the species concept in tick taxonomy. While the typological species concept is still the most widely used, the biological and phylogenetic species concepts are growing in popularity among tick taxonomists. The integrative approach is increasingly being used, but the question is how to define a tick species when using this approach, particularly if data are incongruent. The adoption of an integrative species concept is discussed, in light of recent advances in our understanding of the genetics, morphology, and biology of ticks.

Family-based germline sequencing in children with cancer

The discovery of cancer-predisposing syndromes (CPSs) using next-generation sequencing (NGS) technologies is of increasing importance in pediatric oncology with regard to diagnosis, treatment, surveillance, family counselling and research. Recent studies indicate that a considerable percentage of childhood cancers are associated with CPSs. However, the ratio of CPSs that are caused by inherited vs. de novo mutations (DNMs), the risk of recurrence, and even the total number of genes, which should be considered as a true cancer-predisposing gene, are still unknown. In contrast to sequencing only single index patients, family-based NGS of the germline is a very powerful tool for providing unique insights into inheritance patterns (e.g., DNMs, parental mosaicism) and types of aberrations (e.g., SNV, CNV, indels, SV). Furthermore, functional perturbations of key cancer pathways (e.g., TP53, FA/BRCA) by at least two co-inherited heterozygous digenic mutations from each parent and currently unrecognized rare variants and unmeasured genetic interactions between common and rare variants may be a widespread genetic phenomenon in the germline of affected children. Therefore, family-based trio sequencing has the potential to reveal a striking new landscape of inheritance in childhood cancer and to facilitate the integration and efforts of individualized treatment strategies, including personalized and preventive medicine and cancer surveillance programs. Consequently, cancer genetics is becoming an increasingly common approach in modern oncology, so trio-sequencing should also be routinely integrated into pediatric oncology.

Reproductive Longevity Predicts Mutation Rates in Primates

Mutation rates vary between species across several orders of magnitude, with larger organisms having the highest per-generation mutation rates. Hypotheses for this pattern typically invoke physiological or population-genetic constraints imposed on the molecular machinery preventing mutations [1]. However, continuing germline cell division in multicellular eukaryotes means that organisms with longer generation times and of larger size will leave more mutations to their offspring simply as a byproduct of their increased lifespan [2, 3]. Here, we deeply sequence the genomes of 30 owl monkeys (Aotus nancymaae) from six multi-generation pedigrees to demonstrate that paternal age is the major factor determining the number of de novo mutations in this species. We find that owl monkeys have an average mutation rate of 0.81 × 10^-8 per site per generation, roughly 32% lower than the estimate in humans. Based on a simple model of reproductive longevity that does not require any changes to the mutational machinery, we show that this is the expected mutation rate in owl monkeys. We further demonstrate that our model predicts species-specific mutation rates in other primates, including study-specific mutation rates in humans based on the average paternal age. Our results suggest that variation in life history traits alone can explain variation in the per-generation mutation rate among primates, and perhaps among a wide range of multicellular organisms.

Prioritized High-Confidence Risk Genes for Intellectual Disability Reveal Molecular Convergence During Brain Development

Dissecting the genetic susceptibility to intellectual disability (ID) based on de novo mutations (DNMs) will aid our understanding of the neurobiological and genetic basis of ID. In this study, we identify 63 high-confidence ID genes with q-values < 0.1 based on four background DNM rates and coding DNM data sets from multiple sequencing cohorts. Bioinformatic annotations revealed a higher burden of these 63 ID genes in FMRP targets and CHD8 targets, and these genes show evolutionary constraint against functional genetic variation. Moreover, these ID risk genes were preferentially expressed in the cortical regions from the early fetal to late mid-fetal stages. In particular, a genome-wide weighted co-expression network analysis suggested that ID genes tightly converge onto two biological modules (M1 and M2) during human brain development. Functional annotations showed specific enrichment of chromatin modification and transcriptional regulation for M1 and synaptic function for M2, implying the divergent etiology of the two modules. In addition, we curated 12 additional strong ID risk genes whose molecular interconnectivity with known ID genes (q-values < 0.3) was greater than random. These findings further highlight the biological convergence of ID risk genes and help improve our understanding of the genetic architecture of ID.

Extremely rare variants reveal patterns of germline mutation rate heterogeneity in humans

A detailed understanding of the genome-wide variability of single-nucleotide germline mutation rates is essential to studying human genome evolution. Here, we use ~36 million singleton variants from 3560 whole-genome sequences to infer fine-scale patterns of mutation rate heterogeneity. Mutability is jointly affected by adjacent nucleotide context and diverse genomic features of the surrounding region, including histone modifications, replication timing, and recombination rate, sometimes suggesting specific mutagenic mechanisms. Remarkably, GC content, DNase hypersensitivity, CpG islands, and H3K36 trimethylation are associated with both increased and decreased mutation rates depending on nucleotide context. We validate these estimated effects in an independent dataset of ~46,000 de novo mutations, and confirm our estimates are more accurate than previously published results based on ancestrally older variants without considering genomic features. Our results thus provide the most refined portrait to date of the factors contributing to genome-wide variability of the human germline mutation rate.

Background selection and biased gene conversion affect more than 95% of the human genome and bias demographic inferences

Disentangling the effect on genomic diversity of natural selection from that of demography is notoriously difficult, but necessary to properly reconstruct the history of species. Here, we use high-quality human genomic data to show that purifying selection at linked sites (i.e. background selection, BGS) and GC-biased gene conversion (gBGC) together affect as much as 95% of the variants of our genome. We find that the magnitude and relative importance of BGS and gBGC are largely determined by variation in recombination rate and base composition. Importantly, synonymous sites and non-transcribed regions are also affected, albeit to different degrees. Their use for demographic inference can lead to strong biases. However, by conditioning on genomic regions with recombination rates above 1.5 cM/Mb and mutation types (C↔G, A↔T), we identify a set of SNPs that is mostly unaffected by BGS or gBGC, and that avoids these biases in the reconstruction of human history.

Insights Into de novo Mutation Variation in Lithuanian Exome

In the last decade, one of the biggest challenges in genomics research has been to distinguish definitive pathogenic variants from all likely pathogenic variants identified by next-generation sequencing. This task is particularly complex because of our lack of knowledge regarding overall genome variation and pathogenicity of the variants. Therefore, obtaining sufficient information about genome variants in the general population is necessary as such data could be used for the interpretation of de novo mutations (DNMs) in the context of patient's phenotype in cases of sporadic genetic disease. In this study, data from whole-exome sequencing of the general population in Lithuania were directly examined. In total, 84 (VarScan) and 95 (VarSeqTM) DNMs were identified and validated using different algorithms. Thirty-nine of these mutations were considered likely to be pathogenic based on gene function, evolutionary conservation, and mutation impact. The mutation rate estimated per position pair per generation was 2.74 × 10-8 [95% CI: 2.24 × 10-8-3.35 × 10-8] (VarScan) and 2.4 × 10-8 [95% CI: 1.96 × 10-8-2.99 × 10-8] (VarSeqTM), with 1.77 × 10-8 [95% CI: 6.03 × 10-9-5.2 × 10-8] de novo indels per position per generation. The rate of germline DNMs in the Lithuanian population and the effects of the genomic and epigenetic context on DNM formation were calculated for the first time in this study, providing a basis for further analysis of DNMs in individuals with genetic diseases. Considering these findings, additional studies in patient groups with genetic diseases with unclear etiology may facilitate our ability to distinguish certain pathogenic or adaptive DNMs from tolerated background DNMs and to reliably identify disease-causing DNMs by their properties through direct observation.