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

A familial, telomere-to-telomere reference for human de novo mutation and recombination from a four-generation pedigree

Using five complementary short- and long-read sequencing technologies, we phased and assembled >95% of each diploid human genome in a four-generation, 28-member family (CEPH 1463) allowing us to systematically assess de novo mutations (DNMs) and recombination. From this family, we estimate an average of 192 DNMs per generation, including 75.5 de novo single-nucleotide variants (SNVs), 7.4 non-tandem repeat indels, 79.6 de novo indels or structural variants (SVs) originating from tandem repeats, 7.7 centromeric de novo SVs and SNVs, and 12.4 de novo Y chromosome events per generation. STRs and VNTRs are the most mutable with 32 loci exhibiting recurrent mutation through the generations. We accurately assemble 288 centromeres and six Y chromosomes across the generations, documenting de novo SVs, and demonstrate that the DNM rate varies by an order of magnitude depending on repeat content, length, and sequence identity. We show a strong paternal bias (75-81%) for all forms of germline DNM, yet we estimate that 17% of de novo SNVs are postzygotic in origin with no paternal bias. We place all this variation in the context of a high-resolution recombination map (∼3.5 kbp breakpoint resolution). We observe a strong maternal recombination bias (1.36 maternal:paternal ratio) with a consistent reduction in the number of crossovers with increasing paternal (r=0.85) and maternal (r=0.65) age. However, we observe no correlation between meiotic crossover locations and de novo SVs, arguing against non-allelic homologous recombination as a predominant mechanism. The use of multiple orthogonal technologies, near-telomere-to-telomere phased genome assemblies, and a multi-generation family to assess transmission has created the most comprehensive, publicly available "truth set" of all classes of genomic variants. The resource can be used to test and benchmark new algorithms and technologies to understand the most fundamental processes underlying human genetic variation.

Whole genome sequencing study of identical twins discordant for psychosis

Monozygotic (MZ) twins are often thought to have identical genomes, but recent work has shown that early post-zygotic events can result in a spectrum of DNA variants that are different between MZ twins. Such variants may explain phenotypic discordance and contribute to disease etiology. Here we performed whole genome sequencing in 17 pairs of MZ twins discordant for a psychotic disorder (schizophrenia, schizoaffective disorder or bipolar disorder). We examined various classes of rare variants that are discordant within a twin pair. We identified four genes harboring rare, predicted deleterious missense variants that were private to an affected individual in the cohort. Variants in FOXN1 and FLOT2 would have been categorized as damaging from recent schizophrenia and bipolar exome sequencing studies. Additionally, we identified four rare genic copy number variants (CNVs) private to an affected sample, two of which overlapped genes that have shown evidence for association with schizophrenia or bipolar disorder. One such CNV was a 3q29 duplication previously implicated in autism and developmental delay. We have performed the largest MZ twin study for discordant psychotic phenotypes to date. These findings warrant further investigation using other analytical approaches.

SMAD4 mutations causing Myhre syndrome are under positive selection in the male germline

While it is widely thought that de novo mutations (DNMs) occur randomly, we previously showed that some DNMs are enriched because they are positively selected in the testes of aging men. These "selfish" mutations cause disorders with a shared presentation of features, including exclusive paternal origin, significant increase of the father's age, and high apparent germline mutation rate. To date, all known selfish mutations cluster within the components of the RTK-RAS-MAPK signaling pathway, a critical modulator of testicular homeostasis. Here, we demonstrate the selfish nature of the SMAD4 DNMs causing Myhre syndrome (MYHRS). By analyzing 16 informative trios, we show that MYHRS-causing DNMs originated on the paternally derived allele in all cases. We document a statistically significant epidemiological paternal age effect of 6.3 years excess for fathers of MYHRS probands. We developed an ultra-sensitive assay to quantify spontaneous MYHRS-causing SMAD4 variants in sperm and show that pathogenic variants at codon 500 are found at elevated level in sperm of most men and exhibit a strong positive correlation with donor's age, indicative of a high apparent germline mutation rate. Finally, we performed in vitro assays to validate the peculiar functional behavior of the clonally selected DNMs and explored the basis of the pathophysiology of the different SMAD4 sperm-enriched variants. Taken together, these data provide compelling evidence that SMAD4, a gene operating outside the canonical RAS-MAPK signaling pathway, is associated with selfish spermatogonial selection and raises the possibility that other genes/pathways are under positive selection in the aging human testis.

Elucidating the Transcriptional States of Spermatogenesis-Joint Analysis of Germline and Supporting Cell, Mice and Human, Normal and Perturbed, Bulk and Single-Cell RNA-Seq

In studying the molecular underpinning of spermatogenesis, we expect to understand the fundamental biological processes better and potentially identify genes that may lead to novel diagnostic and therapeutic strategies toward precision medicine in male infertility. In this review, we emphasized our perspective that the path forward necessitates integrative studies that rely on complementary approaches and types of data. To comprehensively analyze spermatogenesis, this review proposes four axes of integration. First, spanning the analysis of spermatogenesis in the healthy state alongside pathologies. Second, the experimental analysis of model systems (in which we can deploy treatments and perturbations) alongside human data. Third, the phenotype is measured alongside its underlying molecular profiles using known markers augmented with unbiased profiles. Finally, the testicular cells are studied as ecosystems, analyzing the germ cells alongside the states observed in the supporting somatic cells. Recently, the study of spermatogenesis has been advancing using single-cell RNA sequencing, where scientists have uncovered the unique stages of germ cell development in mice, revealing new regulators of spermatogenesis and previously unknown cell subtypes in the testis. An in-depth analysis of meiotic and postmeiotic stages led to the discovery of marker genes for spermatogonia, Sertoli and Leydig cells and further elucidated all the other germline and somatic cells in the testis microenvironment in normal and pathogenic conditions. The outcome of an integrative analysis of spermatogenesis using advanced molecular profiling technologies such as scRNA-seq has already propelled our biological understanding, with additional studies expected to have clinical implications for the study of male fertility. By uncovering new genes and pathways involved in abnormal spermatogenesis, we may gain insights into subfertility or sterility.

Transgenerational transmission of post-zygotic mutations suggests symmetric contribution of first two blastomeres to human germline

Little is known about the origin of germ cells in humans. We previously leveraged post-zygotic mutations to reconstruct zygote-rooted cell lineage ancestry trees in a phenotypically normal woman, termed NC0. Here, by sequencing the genome of her children and their father, we analyzed the transmission of early pregastrulation lineages and corresponding mutations across human generations. We found that the germline in NC0 is polyclonal and is founded by at least two cells likely descending from the two blastomeres arising from the first zygotic cleavage. Analyses of public data from several multi-children families and from 1,934 familial quads confirmed this finding in larger cohorts, revealing that known imbalances of up to 90:10 in early lineages allocation in somatic tissues are not reflected in transmission to offspring, establishing a fundamental difference in lineage allocation between the soma and the germline. Analyses of all the data consistently suggest that germline has a balanced 50:50 lineage allocation from the first two blastomeres.

The clock-like accumulation of germline and somatic mutations can arise from the interplay of DNA damage and repair

The rates at which mutations accumulate across human cell types vary. To identify causes of this variation, mutations are often decomposed into a combination of the single-base substitution (SBS) "signatures" observed in germline, soma, and tumors, with the idea that each signature corresponds to one or a small number of underlying mutagenic processes. Two such signatures turn out to be ubiquitous across cell types: SBS signature 1, which consists primarily of transitions at methylated CpG sites thought to be caused by spontaneous deamination, and the more diffuse SBS signature 5, which is of unknown etiology. In cancers, the number of mutations attributed to these 2 signatures accumulates linearly with age of diagnosis, and thus the signatures have been termed "clock-like." To better understand this clock-like behavior, we develop a mathematical model that includes DNA replication errors, unrepaired damage, and damage repaired incorrectly. We show that mutational signatures can exhibit clock-like behavior because cell divisions occur at a constant rate and/or because damage rates remain constant over time, and that these distinct sources can be teased apart by comparing cell lineages that divide at different rates. With this goal in mind, we analyze the rate of accumulation of mutations in multiple cell types, including soma as well as male and female germline. We find no detectable increase in SBS signature 1 mutations in neurons and only a very weak increase in mutations assigned to the female germline, but a significant increase with time in rapidly dividing cells, suggesting that SBS signature 1 is driven by rounds of DNA replication occurring at a relatively fixed rate. In contrast, SBS signature 5 increases with time in all cell types, including postmitotic ones, indicating that it accumulates independently of cell divisions; this observation points to errors in DNA repair as the key underlying mechanism. Thus, the two "clock-like" signatures observed across cell types likely have distinct origins, one set by rates of cell division, the other by damage rates.

[Sperm Mosaic Variants and Their Influence on the Offspring]

Genomic mosaicism arising from mosaic variants is a phenomenon that describes the presence of a cell or cell populations with different genome compositions from the germline cells of an individual. It comprises all types of genetic variants. A large proportion of childhood genetic disorders are defined as being de novo, meaning that the disease-causing mutations are only detected in the proband, not in any of the parents. Population studies show that 80% of the de novo mutations arise from the paternal haplotype, that is, from paternal sperm mosaicism. This review provides a summary of the types and detection strategies of sperm mosaicism. In addition, it provides discussions on how recent studies demonstrated that genomic mosaic mutations in parents, especially those in the paternal sperms, could be inherited by the offspring and cause childhood disorders. According to the previous findings of the author's research team, sperm mosaicism derived from early embryogenesis and primordial germ cell stages can explain 5% to 20% of the de novo mutations related to clinical phenotypes and can serve as an important predictor of both rare and complex disorders. Sperm mosaicism shows great potential for clinical genetic diagnosis and consultations. Based on the published literature, the author suggests that, large-scale screening for de novo sperm mosaic mutations and population-based genetic screening should be conducted in future studies, which will greatly enhance the risk assessment in the offspring and effectively improve the genetic health at the population level. Implementation of direct sperm detection for de novo mutations will significantly increase the efficiency of the stratification of patient cohorts and improve recurrence risk assessment for future births. Future research in the field should be focused on the impact of environmental and lifestyle factors on the health of the offspring through sperms and their modeling of mutation signatures. In addition, targeted in vitro modeling of sperm mutations will also be a promising direction.

Human embryonic genetic mosaicism and its effects on development and disease

Nearly every mammalian cell division is accompanied by a mutational event that becomes fixed in a daughter cell. When carried forward to additional cell progeny, a clone of variant cells can emerge. As a result, mammals are complex mosaics of clones that are genetically distinct from one another. Recent high-throughput sequencing studies have revealed that mosaicism is common, clone sizes often increase with age and specific variants can affect tissue function and disease development. Variants that are acquired during early embryogenesis are shared by multiple cell types and can affect numerous tissues. Within tissues, variant clones compete, which can result in their expansion or elimination. Embryonic mosaicism has clinical implications for genetic disease severity and transmission but is likely an under-recognized phenomenon. To better understand its implications for mosaic individuals, it is essential to leverage research tools that can elucidate the mechanisms by which expanded embryonic variants influence development and disease.

Effects of parental age and polymer composition on short tandem repeat de novo mutation rates

Short tandem repeats (STRs) are hotspots of genomic variability in the human germline because of their high mutation rates, which have long been attributed largely to polymerase slippage during DNA replication. This model suggests that STR mutation rates should scale linearly with a father's age, as progenitor cells continually divide after puberty. In contrast, it suggests that STR mutation rates should not scale with a mother's age at her child's conception, since oocytes spend a mother's reproductive years arrested in meiosis II and undergo a fixed number of cell divisions that are independent of the age at ovulation. Yet, mirroring recent findings, we find that STR mutation rates covary with paternal and maternal age, implying that some STR mutations are caused by DNA damage in quiescent cells rather than polymerase slippage in replicating progenitor cells. These results echo the recent finding that DNA damage in oocytes is a significant source of de novo single nucleotide variants and corroborate evidence of STR expansion in postmitotic cells. However, we find that the maternal age effect is not confined to known hotspots of oocyte mutagenesis, nor are postzygotic mutations likely to contribute significantly. STR nucleotide composition demonstrates divergent effects on de novo mutation (DNM) rates between sexes. Unlike the paternal lineage, maternally derived DNMs at A/T STRs display a significantly greater association with maternal age than DNMs at G/C-containing STRs. These observations may suggest the mechanism and developmental timing of certain STR mutations and contradict prior attribution of replication slippage as the primary mechanism of STR mutagenesis.

Automated Identification of Germline de novo Mutations in Family Trios: A Consensus-Based Informatic Approach

Accurate identification of germline de novo variants (DNVs) remains a challenging problem despite rapid advances in sequencing technologies as well as methods for the analysis of the data they generate, with putative solutions often involving ad hoc filters and visual inspection of identified variants. Here, we present a purely informatic method for the identification of DNVs by analyzing short-read genome sequencing data from proband-parent trios. Our method evaluates variant calls generated by three genome sequence analysis pipelines utilizing different algorithms-GATK HaplotypeCaller, DeepTrio and Velsera GRAF-exploring the assumption that a requirement of consensus can serve as an effective filter for high-quality DNVs. We assessed the efficacy of our method by testing DNVs identified using a previously established, highly accurate classification procedure that partially relied on manual inspection and used Sanger sequencing to validate a DNV subset comprising less confident calls. The results show that our method is highly precise and that applying a force-calling procedure to putative variants further removes false-positive calls, increasing precision of the workflow to 99.6%. Our method also identified novel DNVs, 87% of which were validated, indicating it offers a higher recall rate without compromising accuracy. We have implemented this method as an automated bioinformatics workflow suitable for large-scale analyses without need for manual intervention.

Meta-analysis of 46,000 germline de novo mutations linked to human inherited disease

BACKGROUND

De novo mutations (DNMs) are variants that occur anew in the offspring of noncarrier parents. They are not inherited from either parent but rather result from endogenous mutational processes involving errors of DNA repair/replication. These spontaneous errors play a significant role in the causation of genetic disorders, and their importance in the context of molecular diagnostic medicine has become steadily more apparent as more DNMs have been reported in the literature. In this study, we examined 46,489 disease-associated DNMs annotated by the Human Gene Mutation Database (HGMD) to ascertain their distribution across gene and disease categories.

RESULTS

Most disease-associated DNMs reported to date are found to be associated with developmental and psychiatric disorders, a reflection of the focus of sequencing efforts over the last decade. Of the 13,277 human genes in which DNMs have so far been found, the top-10 genes with the highest proportions of DNM relative to gene size were H3-3 A, DDX3X, CSNK2B, PURA, ZC4H2, STXBP1, SCN1A, SATB2, H3-3B and TUBA1A. The distribution of CADD and REVEL scores for both disease-associated DNMs and those mutations not reported to be de novo revealed a trend towards higher deleteriousness for DNMs, consistent with the likely lower selection pressure impacting them. This contrasts with the non-DNMs, which are presumed to have been subject to continuous negative selection over multiple generations.

CONCLUSION

This meta-analysis provides important information on the occurrence and distribution of disease-associated DNMs in association with heritable disease and should make a significant contribution to our understanding of this major type of mutation.

Epistasis between mutator alleles contributes to germline mutation spectrum variability in laboratory mice

Maintaining germline genome integrity is essential and enormously complex. Although many proteins are involved in DNA replication, proofreading, and repair, mutator alleles have largely eluded detection in mammals. DNA replication and repair proteins often recognize sequence motifs or excise lesions at specific nucleotides. Thus, we might expect that the spectrum of de novo mutations - the frequencies of C>T, A>G, etc. - will differ between genomes that harbor either a mutator or wild-type allele. Previously, we used quantitative trait locus mapping to discover candidate mutator alleles in the DNA repair gene Mutyh that increased the C>A germline mutation rate in a family of inbred mice known as the BXDs (Sasani et al., 2022, Ashbrook et al., 2021). In this study we developed a new method to detect alleles associated with mutation spectrum variation and applied it to mutation data from the BXDs. We discovered an additional C>A mutator locus on chromosome 6 that overlaps Ogg1, a DNA glycosylase involved in the same base-excision repair network as Mutyh (David et al., 2007). Its effect depends on the presence of a mutator allele near Mutyh, and BXDs with mutator alleles at both loci have greater numbers of C>A mutations than those with mutator alleles at either locus alone. Our new methods for analyzing mutation spectra reveal evidence of epistasis between germline mutator alleles and may be applicable to mutation data from humans and other model organisms.

Familial factors rather than paternal age contribute to the aetiology of epilepsy

BACKGROUND

Whether paternal age associated with offspring's epilepsy risk is a cause of de novo mutation as men age, or just an association due to confounding factors, is still unclear.

METHODS

We performed a population-based, multi-generation and sibling comparison study in Taiwan, which included 2 751 232 singletons born in 2001-17 who were followed until 2020. Of these, 819 371/826 087 with information on paternal/maternal grandparents were selected for multi-generation analyses and 1 748 382 with sibling(s) were selected for sibling comparison. Cox proportional hazard regression was used to estimate the hazard ratio (HR) and 95% confidence interval (CI).

RESULTS

In the total cohort, there was an increased risk of epilepsy in individuals with advanced paternal age, e.g. the HR for paternal age ≥50 was1.36 (95% CI: 1.15-1.61) compared with paternal age 25-29, and fathers older than mothers, e.g. the HR for parental age difference ≥15 years was 1.29 (95% CI: 1.16-1.43). When accounting for parental age difference, the association between paternal age and epilepsy in offspring was attenuated (HR for paternal age ≥50 was 1.11, 95% CI: 0.93-1.34). Multi-generation analyses did not support the association of advanced grand-paternal age at childbirth of the parent with offspring's risk of epilepsy. Sibling comparison analyses did not support the association of older paternal age with increased risk of epilepsy (HR was 0.96 for per year increase in paternal age, 95% CI: 0.96-0.97).

CONCLUSIONS

These results do not support the hypothesis that advanced paternal age is associated with epilepsy in offspring. Instead, familial factors may explain the observed paternal age association with the offspring's risk of epilepsy.

Exploring FGFR3 Mutations in the Male Germline: Implications for Clonal Germline Expansions and Paternal Age-Related Dysplasias

Delayed fatherhood results in a higher risk of inheriting a new germline mutation that might result in a congenital disorder in the offspring. In particular, some FGFR3 mutations increase in frequency with age, but there are still a large number of uncharacterized FGFR3 mutations that could be expanding in the male germline with potentially early- or late-onset effects in the offspring. Here, we used digital polymerase chain reaction to assess the frequency and spatial distribution of 10 different FGFR3 missense substitutions in the sexually mature male germline. Our functional assessment of the receptor signaling of the variants with biophysical methods showed that 9 of these variants resulted in a higher activation of the receptor´s downstream signaling, resulting in 2 different expansion behaviors. Variants that form larger subclonal expansions in a dissected postmortem testis also showed a positive correlation of the substitution frequency with the sperm donor's age, and a high and ligand-independent FGFR3 activation. In contrast, variants that measured high FGFR3 signaling and elevated substitution frequencies independent of the donor's age did not result in measurable subclonal expansions in the testis. This suggests that promiscuous signal activation might also result in an accumulation of mutations before the sexual maturation of the male gonad with clones staying relatively constant in size throughout time. Collectively, these results provide novel insights into our understanding of the mutagenesis of driver mutations and their resulting mosaicism in the male germline with important consequences for the transmission and recurrence of associated disorders.

Hotspot propensity across mutational processes

The sparsity of mutations observed across tumours hinders our ability to study mutation rate variability at nucleotide resolution. To circumvent this, here we investigated the propensity of mutational processes to form mutational hotspots as a readout of their mutation rate variability at single base resolution. Mutational signatures 1 and 17 have the highest hotspot propensity (5-78 times higher than other processes). After accounting for trinucleotide mutational probabilities, sequence composition and mutational heterogeneity at 10 Kbp, most (94-95%) signature 17 hotspots remain unexplained, suggesting a significant role of local genomic features. For signature 1, the inclusion of genome-wide distribution of methylated CpG sites into models can explain most (80-100%) of the hotspot propensity. There is an increased hotspot propensity of signature 1 in normal tissues and de novo germline mutations. We demonstrate that hotspot propensity is a useful readout to assess the accuracy of mutation rate models at nucleotide resolution. This new approach and the findings derived from it open up new avenues for a range of somatic and germline studies investigating and modelling mutagenesis.

Effects of parental age and polymer composition on short tandem repeat de novo mutation rates

Short tandem repeats (STRs) are hotspots of genomic variability in the human germline because of their high mutation rates, which have long been attributed largely to polymerase slippage during DNA replication. This model suggests that STR mutation rates should scale linearly with a father's age, as progenitor cells continually divide after puberty. In contrast, it suggests that STR mutation rates should not scale with a mother's age at her child's conception, since oocytes spend a mother's reproductive years arrested in meiosis II and undergo a fixed number of cell divisions that are independent of the age at ovulation. Yet, mirroring recent findings, we find that STR mutation rates covary with paternal and maternal age, implying that some STR mutations are caused by DNA damage in quiescent cells rather than the classical mechanism of polymerase slippage in replicating progenitor cells. These results also echo the recent finding that DNA damage in quiescent oocytes is a significant source of de novo SNVs and corroborate evidence of STR expansion in postmitotic cells. However, we find that the maternal age effect is not confined to previously discovered hotspots of oocyte mutagenesis, nor are post-zygotic mutations likely to contribute significantly. STR nucleotide composition demonstrates divergent effects on DNM rates between sexes. Unlike the paternal lineage, maternally derived DNMs at A/T STRs display a significantly greater association with maternal age than DNMs at GC-containing STRs. These observations may suggest the mechanism and developmental timing of certain STR mutations and are especially surprising considering the prior belief in replication slippage as the dominant mechanism of STR mutagenesis.

Characterizing genetic variation in the regulation of the ER stress response through computational and cis-eQTL analyses

Misfolded proteins in the endoplasmic reticulum (ER) elicit the ER stress response, a large transcriptional response driven by 3 well-characterized transcription factors (TFs). This transcriptional response is variable across different genetic backgrounds. One mechanism in which genetic variation can lead to transcriptional variability in the ER stress response is through altered binding and activity of the 3 main TFs: XBP1, ATF6, and ATF4. This work attempts to better understand this mechanism by first creating a computational pipeline to identify potential binding sites throughout the human genome. We utilized GTEx data sets to identify cis-eQTLs that fall within predicted TF binding sites (TFBSs). We also utilized the ClinVar database to compare the number of pathogenic vs benign variants at different positions of the binding motifs. Finally, we performed a cis-eQTL analysis on human cell lines experiencing ER stress to identify cis-eQTLs that regulate the variable ER stress response. The majority of these cis-eQTLs are unique to a given condition: control or ER stress. Some of these stress-specific cis-eQTLs fall within putative binding sites of the 3 main ER stress response TFs, providing a potential mechanism by which these cis-eQTLs might be impacting gene expression under ER stress conditions through altered TF binding. This study represents the first cis-eQTL analysis on human samples experiencing ER stress and is a vital step toward identifying the genetic components responsible for the variable ER stress response.

Nutrigenomics and redox regulation: Concepts relating to the Special Issue on nutrigenomics

During our whole lifespan, from conception to death, the epigenomes of all tissues and cell types of our body integrate signals from the environment. This includes signals derived from our diet and the uptake of macro- and micronutrients. In most cases, this leads only to transient changes, but some effects of this epigenome programming process are persistent and can even be transferred to the next generation. Both epigenetic programming and redox processes are affected by the individual choice of diet and other lifestyle decisions like physical activity. The nutrient-gene communication pathways have adapted during human evolution and are essential for maintaining health. However, when they are maladaptive, such as in long-term obesity, they significantly contribute to diseases like type 2 diabetes and cancer. The field of nutrigenomics investigates nutrition-related signal transduction pathways and their effect on gene expression involving interactions both with the genome and the epigenomes. Several of these diet-(epi)genome interactions and the involved signal transduction cascades are redox-regulated. Examples include the effects of the NAD+/NADH ratio, vitamin C levels and secondary metabolites of dietary molecules from plants on the acetylation and methylation state of the epigenome as well as on gene expression through redox-sensitive pathways via the transcription factors NFE2L2 and FOXO. In this review, we summarize and extend on these topics as well as those discussed in the articles of this Special Issue and take them into the context of redox biology.

Germline mutation rate is elevated in young and old parents in Caenorhabditis remanei

The effect of parental age on germline mutation rate across generations is not fully understood. While some studies report a positive linear relationship of mutation rate with increasing age, others suggest that mutation rate varies with age but not in a linear fashion. We investigated the effect of parental age on germline mutations by generating replicated mutation accumulation lines in Caenorhabditis remanei at three parental ages ("Young T1" [Day 1], "Peak T2" [Day 2], and "Old T5" [Day 5] parents). We conducted whole-genome resequencing and variant calling to compare differences in mutation rates after three generations of mutation accumulation. We found that Peak T2 lines had an overall reduced mutation rate compared to Young T1 and Old T5 lines, but this pattern of the effect varied depending on the variant impact. Specifically, we found no high-impact variants in Peak T2 lines, and modifiers and up- and downstream gene variants were less frequent in these lines. These results suggest that animals at the peak of reproduction have better DNA maintenance and repair compared to young and old animals. We propose that C. remanei start to reproduce before they optimize their DNA maintenance and repair, trading the benefits of earlier onset of reproduction against offspring mutation load. The increase in offspring mutation load with age likely represents germline senescence.

The 'faulty male' hypothesis for sex-biased mutation and disease

Biological differences between males and females lead to many differences in physiology, disease, and overall health. One of the most prominent disparities is in the number of germline mutations passed to offspring: human males transmit three times as many mutations as do females. While the classic explanation for this pattern invokes differences in post-puberty germline replication between the sexes, recent whole-genome evidence in humans and other mammals has cast doubt on this mechanism. Here, we review recent work that is inconsistent with a replication-driven model of male-biased mutation, and propose an alternative, 'faulty male' hypothesis. This model proposes that males are less able to repair and/or protect DNA from damage compared to females. Importantly, we suggest that this new model for male-biased mutation may also help to explain several pronounced differences between the sexes in cancer, aging, and DNA repair. Although the detailed contributions of genetic, epigenetic, and hormonal influences of biological sex on mutation remain to be fully understood, a reconsideration of the mechanisms underlying these differences will lead to a deeper understanding of evolution and disease.

Epistasis between mutator alleles contributes to germline mutation spectra variability in laboratory mice

Maintaining germline genome integrity is essential and enormously complex. Although many proteins are involved in DNA replication, proofreading, and repair [1], mutator alleles have largely eluded detection in mammals. DNA replication and repair proteins often recognize sequence motifs or excise lesions at specific nucleotides. Thus, we might expect that the spectrum of de novo mutations - the frequencies of C>T, A>G, etc. - will differ between genomes that harbor either a mutator or wild-type allele. Previously, we used quantitative trait locus mapping to discover candidate mutator alleles in the DNA repair gene Mutyh that increased the C>A germline mutation rate in a family of inbred mice known as the BXDs [2,3]. In this study we developed a new method to detect alleles associated with mutation spectrum variation and applied it to mutation data from the BXDs. We discovered an additional C>A mutator locus on chromosome 6 that overlaps Ogg1, a DNA glycosylase involved in the same base-excision repair network as Mutyh [4]. Its effect depended on the presence of a mutator allele near Mutyh, and BXDs with mutator alleles at both loci had greater numbers of C>A mutations than those with mutator alleles at either locus alone. Our new methods for analyzing mutation spectra reveal evidence of epistasis between germline mutator alleles and may be applicable to mutation data from humans and other model organisms.

Validation of machine learning approach for direct mutation rate estimation

Mutations are the primary source of all genetic variation. Knowledge about their rates is critical for any evolutionary genetic analyses, but for a long time, that knowledge has remained elusive and indirectly inferred. In recent years, parent-offspring comparisons have yielded the first direct mutation rate estimates. The analyses are, however, challenging due to high rate of false positives and no consensus regarding standardized filtering of candidate de novo mutations. Here, we validate the application of a machine learning approach for such a task and estimate the mutation rate for the guppy (Poecilia reticulata), a model species in eco-evolutionary studies. We sequenced 4 parents and 20 offspring, followed by screening their genomes for de novo mutations. The initial large number of candidate de novo mutations was hard-filtered to remove false-positive results. These results were compared with mutation rate estimated with a supervised machine learning approach. Both approaches were followed by molecular validation of all candidate de novo mutations and yielded similar results. The ML method uniquely identified three mutations, but overall required more hands-on curation and had higher rates of false positives and false negatives. Both methods concordantly showed no difference in mutation rates between families. Estimated here the guppy mutation rate is among the lowest directly estimated mutation rates in vertebrates; however, previous research has also found low estimated rates in other teleost fishes. We discuss potential explanations for such a pattern, as well as future utility and limitations of machine learning approaches.

Scaling the discrete-time Wright-Fisher model to biobank-scale datasets

The discrete-time Wright-Fisher (DTWF) model and its diffusion limit are central to population genetics. These models can describe the forward-in-time evolution of allele frequencies in a population resulting from genetic drift, mutation, and selection. Computing likelihoods under the diffusion process is feasible, but the diffusion approximation breaks down for large samples or in the presence of strong selection. Existing methods for computing likelihoods under the DTWF model do not scale to current exome sequencing sample sizes in the hundreds of thousands. Here, we present a scalable algorithm that approximates the DTWF model with provably bounded error. Our approach relies on two key observations about the DTWF model. The first is that transition probabilities under the model are approximately sparse. The second is that transition distributions for similar starting allele frequencies are extremely close as distributions. Together, these observations enable approximate matrix-vector multiplication in linear (as opposed to the usual quadratic) time. We prove similar properties for Hypergeometric distributions, enabling fast computation of likelihoods for subsamples of the population. We show theoretically and in practice that this approximation is highly accurate and can scale to population sizes in the tens of millions, paving the way for rigorous biobank-scale inference. Finally, we use our results to estimate the impact of larger samples on estimating selection coefficients for loss-of-function variants. We find that increasing sample sizes beyond existing large exome sequencing cohorts will provide essentially no additional information except for genes with the most extreme fitness effects.

Disentangling sources of clock-like mutations in germline and soma

The rates of mutations vary across cell types. To identify causes of this variation, mutations are often decomposed into a combination of the single base substitution (SBS) "signatures" observed in germline, soma and tumors, with the idea that each signature corresponds to one or a small number of underlying mutagenic processes. Two such signatures turn out to be ubiquitous across cell types: SBS signature 1, which consists primarily of transitions at methylated CpG sites caused by spontaneous deamination, and the more diffuse SBS signature 5, which is of unknown etiology. In cancers, the number of mutations attributed to these two signatures accumulates linearly with age of diagnosis, and thus the signatures have been termed "clock-like." To better understand this clock-like behavior, we develop a mathematical model that includes DNA replication errors, unrepaired damage, and damage repaired incorrectly. We show that mutational signatures can exhibit clock-like behavior because cell divisions occur at a constant rate and/or because damage rates remain constant over time, and that these distinct sources can be teased apart by comparing cell lineages that divide at different rates. With this goal in mind, we analyze the rate of accumulation of mutations in multiple cell types, including soma as well as male and female germline. We find no detectable increase in SBS signature 1 mutations in neurons and only a very weak increase in mutations assigned to the female germline, but a significant increase with time in rapidly-dividing cells, suggesting that SBS signature 1 is driven by rounds of DNA replication occurring at a relatively fixed rate. In contrast, SBS signature 5 increases with time in all cell types, including post-mitotic ones, indicating that it accumulates independently of cell divisions; this observation points to errors in DNA repair as the key underlying mechanism. Thus, the two "clock-like" signatures observed across cell types likely have distinct origins, one set by rates of cell division, the other by damage rates.

The rate of de novo structural variation is increased in in vitro-produced offspring and preferentially affects the paternal genome

Assisted reproductive technologies (ARTs), including in vitro maturation and fertilization (IVF), are increasingly used in human and animal reproduction. Whether these technologies directly affect the rate of de novo mutation (DNM), and to what extent, has been a matter of debate. Here we take advantage of domestic cattle, characterized by complex pedigrees that are ideally suited to detect DNMs and by the systematic use of ART, to study the rate of de novo structural variation (dnSV) in this species and how it is impacted by IVF. By exploiting features of associated de novo point mutations (dnPMs) and dnSVs in clustered DNMs, we provide strong evidence that (1) IVF increases the rate of dnSV approximately fivefold, and (2) the corresponding mutations occur during the very early stages of embryonic development (one- and two-cell stage), yet primarily affect the paternal genome.

De Novo Mutation Rates in Sticklebacks

Mutation rate is a fundamental parameter in population genetics. Apart from being an important scaling parameter for demographic and phylogenetic inference, it allows one to understand at what rate new genetic diversity is generated and what the expected level of genetic diversity is in a population at equilibrium. However, except for well-established model organisms, accurate estimates of de novo mutation rates are available for a very limited number of organisms from the wild. We estimated mutation rates (µ) in two marine populations of the nine-spined stickleback (Pungitius pungitius) with the aid of several 2- and 3-generational family pedigrees, deep (>50×) whole-genome resequences and a high-quality reference genome. After stringent filtering, we discovered 308 germline mutations in 106 offspring translating to µ = 4.83 × 10-9 and µ = 4.29 × 10-9 per base per generation in the two populations, respectively. Up to 20% of the mutations were shared by full-sibs showing that the level of parental mosaicism was relatively high. Since the estimated µ was 3.1 times smaller than the commonly used substitution rate, recalibration with µ led to substantial increase in estimated divergence times between different stickleback species. Our estimates of the de novo mutation rate should provide a useful resource for research focused on fish population genetics and that of sticklebacks in particular.

Emergence and influence of sequence bias in evolutionarily malleable, mammalian tandem arrays

BACKGROUND

The radiation of mammals at the extinction of the dinosaurs produced a plethora of new forms-as diverse as bats, dolphins, and elephants-in only 10-20 million years. Behind the scenes, adaptation to new niches is accompanied by extensive innovation in large families of genes that allow animals to contact the environment, including chemosensors, xenobiotic enzymes, and immune and barrier proteins. Genes in these "outward-looking" families are allelically diverse among humans and exhibit tissue-specific and sometimes stochastic expression.

RESULTS

Here, we show that these tandem arrays of outward-looking genes occupy AT-biased isochores and comprise the "tissue-specific" gene class that lack CpG islands in their promoters. Models of mammalian genome evolution have not incorporated the sharply different functions and transcriptional patterns of genes in AT- versus GC-biased regions. To examine the relationship between gene family expansion, sequence content, and allelic diversity, we use population genetic data and comparative analysis. First, we find that AT bias can emerge during evolutionary expansion of gene families in cis. Second, human genes in AT-biased isochores or with GC-poor promoters experience relatively low rates of de novo point mutation today but are enriched for non-synonymous variants. Finally, we find that isochores containing gene clusters exhibit low rates of recombination.

CONCLUSIONS

Our analyses suggest that tolerance of non-synonymous variation and low recombination are two forces that have produced the depletion of GC bases in outward-facing gene arrays. In turn, high AT content exerts a profound effect on their chromatin organization and transcriptional regulation.

Direct additive genetics and maternal effect contribute to the risk of Tourette disorder

BACKGROUND

Risk for Tourette disorder, and chronic motor or vocal tic disorders (referenced here inclusively as CTD), arise from a combination of genetic and environmental factors. While multiple studies have demonstrated the importance of direct additive genetic variation for CTD risk, little is known about the role of cross-generational transmission of genetic risk, such as maternal effect, which is not transmitted via the inherited parental genomes. Here, we partition sources of variation on CTD risk into direct additive genetic effect (narrow-sense heritability) and maternal effect.

METHODS

The study population consists of 2 522 677 individuals from the Swedish Medical Birth Register, who were born in Sweden between 1 January 1973 and 31 December 2000, and followed for a diagnosis of CTD through 31 December, 2013. We used generalised linear mixed models to partition the liability of CTD into: direct additive genetic effect, genetic maternal effect and environmental maternal effect.

RESULTS

We identified 6227 (0.2%) individuals in the birth cohort with a CTD diagnosis. A study of half-siblings showed that maternal half-siblings had twice higher risk of developing a CTD compared with paternal ones. We estimated 60.7% direct additive genetic effect (95% credible interval, 58.5% to 62.4%), 4.8% genetic maternal effect (95% credible interval, 4.4% to 5.1%) and 0.5% environmental maternal effect (95% credible interval, 0.2% to 7%).

CONCLUSIONS

Our results demonstrate genetic maternal effect contributes to the risk of CTD. Failure to account for maternal effect results in an incomplete understanding of the genetic risk architecture of CTD, as the risk for CTD is impacted by maternal effect which is above and beyond the risk from transmitted genetic effect.

Germline de novo mutations in families with Mendelian cancer syndromes caused by defects in DNA repair

DNA repair defects underlie many cancer syndromes. We tested whether de novo germline mutations (DNMs) are increased in families with germline defects in polymerase proofreading or base excision repair. A parent with a single germline POLE or POLD1 mutation, or biallelic MUTYH mutations, had 3-4 fold increased DNMs over sex-matched controls. POLE had the largest effect. The DNMs carried mutational signatures of the appropriate DNA repair deficiency. No DNM increase occurred in offspring of MUTYH heterozygous parents. Parental DNA repair defects caused about 20-150 DNMs per child, additional to the ~60 found in controls, but almost all extra DNMs occurred in non-coding regions. No increase in post-zygotic mutations was detected, excepting a child with bi-allelic MUTYH mutations who was excluded from the main analysis; she had received chemotherapy and may have undergone oligoclonal haematopoiesis. Inherited DNA repair defects associated with base pair-level mutations increase DNMs, but phenotypic consequences appear unlikely.

The landscape of somatic mutations in lymphoblastoid cell lines

Somatic mutations have important biological ramifications while exerting substantial rate, type, and genomic location heterogeneity. Yet, their sporadic occurrence makes them difficult to study at scale and across individuals. Lymphoblastoid cell lines (LCLs), a model system for human population and functional genomics, harbor large numbers of somatic mutations and have been extensively genotyped. By comparing 1,662 LCLs, we report that the mutational landscape of the genome varies across individuals in terms of the number of mutations, their genomic locations, and their spectra; this variation may itself be modulated by somatic trans-acting mutations. Mutations attributed to the translesion DNA polymerase η follow two different modes of formation, with one mode accounting for the hypermutability of the inactive X chromosome. Nonetheless, the distribution of mutations along the inactive X chromosome appears to follow an epigenetic memory of the active form.

MECP2 germline mosaicism plays an important part in the inheritance of Rett syndrome: a study of MECP2 germline mosaicism in males

BACKGROUND

Germline mosaicisms could be inherited to offspring, which considered as "de novo" in most cases. Paternal germline MECP2 mosaicism has been reported in fathers of girls with Rett syndrome (RTT) previously. For further study, we focused on MECP2 germline mosaicism in males, not only RTT fathers.

METHODS

Thirty-two fathers of RTT girls with MECP2 pathogenic mutations and twenty-five healthy adult males without history and family history of RTT or other genetic disorders were recruited. Sperm samples were collected and ten MECP2 hotspot mutations were detected by micro-droplet digital PCR (mDDPCR). And routine semen test was performed at the same time if the sample was sufficient. Additionally, blood samples were also detected for those with sperm MECP2 mosaicisms.

RESULTS

Nine fathers with RTT daughters (28.1%, 9/32) were found to have MECP2 mosaicism in their sperm samples, with the mutant allele fractions (MAFs) ranging from 0.05% to 7.55%. Only one father with MECP2 c.806delG germline mosaicism (MAF 7.55%) was found to have mosaicism in the blood sample, with the MAF was 0.28%. In the group of healthy adult males, MECP2 mosaicism was found in 7 sperm samples (28.0%, 7/25), with the MAFs ranging from 0.05% to 0.18%. None of the healthy adult males with MECP2 germline mosaicisms were found with MECP2 mosaicism in blood samples. There were no statistical differences in age, or the incidence of asthenospermia between fathers with RTT daughters and healthy adult males with MECP2 germline mosaicisms. Additionally, there was no linear correlation between MAFs of MECP2 mosaicisms and the age of males with germline MECP2 mosaicisms.

CONCLUSIONS

Germline MECP2 mosaicism could be found not only in fathers with RTT daughters but also in healthy adult males without family history of RTT. As germline mosaic mutations may be passed on to offspring which commonly known as "de novo", more attention should be paid to germline mosaicism, especially in families with a proband diagnosed with genetic disorders.

Identification and characterisation of de novo germline structural variants in two commercial pig lines using trio-based whole genome sequencing

BACKGROUND

De novo mutations arising in the germline are a source of genetic variation and their discovery broadens our understanding of genetic disorders and evolutionary patterns. Although the number of de novo single nucleotide variants (dnSNVs) has been studied in a number of species, relatively little is known about the occurrence of de novo structural variants (dnSVs). In this study, we investigated 37 deeply sequenced pig trios from two commercial lines to identify dnSVs present in the offspring. The identified dnSVs were characterised by identifying their parent of origin, their functional annotations and characterizing sequence homology at the breakpoints.

RESULTS

We identified four swine germline dnSVs, all located in intronic regions of protein-coding genes. Our conservative, first estimate of the swine germline dnSV rate is 0.108 (95% CI 0.038-0.255) per generation (one dnSV per nine offspring), detected using short-read sequencing. Two detected dnSVs are clusters of mutations. Mutation cluster 1 contains a de novo duplication, a dnSNV and a de novo deletion. Mutation cluster 2 contains a de novo deletion and three de novo duplications, of which one is inverted. Mutation cluster 2 is 25 kb in size, whereas mutation cluster 1 (197 bp) and the other two individual dnSVs (64 and 573 bp) are smaller. Only mutation cluster 2 could be phased and is located on the paternal haplotype. Mutation cluster 2 originates from both micro-homology as well as non-homology mutation mechanisms, where mutation cluster 1 and the other two dnSVs are caused by mutation mechanisms lacking sequence homology. The 64 bp deletion and mutation cluster 1 were validated through PCR. Lastly, the 64 bp deletion and the 573 bp duplication were validated in sequenced offspring of probands with three generations of sequence data.

CONCLUSIONS

Our estimate of 0.108 dnSVs per generation in the swine germline is conservative, due to our small sample size and restricted possibilities of dnSV detection from short-read sequencing. The current study highlights the complexity of dnSVs and shows the potential of breeding programs for pigs and livestock species in general, to provide a suitable population structure for identification and characterisation of dnSVs.

Genome-wide identification of copy neutral loss of heterozygosity reveals its possible association with spatial positioning of chromosomes

Loss of heterozygosity (LOH) is a genetic alteration that results from the loss of one allele at a heterozygous locus. In particular, copy neutral LOH (CN-LOH) events are generated, for example, by mitotic homologous recombination after monoallelic defection or gene conversion, resulting in novel homozygous locus having two copies of the normal counterpart allele. This phenomenon can serve as a source of genome diversity and is associated with various diseases. To clarify the nature of the CN-LOH such as the frequency, genomic distribution and inheritance pattern, we made use of whole-genome sequencing data of the three-generation CEPH/Utah family cohort, with the pedigree consisting of grandparents, parents and offspring. We identified an average of 40.7 CN-LOH events per individual taking advantage of 285 healthy individuals from 33 families in the cohort. On average 65% of them were classified as gonosomal-mosaicism-associated CN-LOH, which exists in both germline and somatic cells. We also confirmed that the incidence of the CN-LOH has little to do with the parents' age and sex. Furthermore, through the analysis of the genomic region including the CN-LOH, we found that the chance of the occurrence of the CN-LOH tends to increase at the GC-rich locus and/or on the chromosome having a relatively close inter-homolog distance. We expect that these results provide significant insights into the association between genetic alteration and spatial position of chromosomes as well as the intrinsic genetic property of the CN-LOH.

Revisiting mutagenesis at non-B DNA motifs in the human genome

Non-B DNA structures formed by repetitive sequence motifs are known instigators of mutagenesis in experimental systems. Analyzing this phenomenon computationally in the human genome requires careful disentangling of intrinsic confounding factors, including overlapping and interrupted motifs and recurrent sequencing errors. Here, we show that accounting for these factors eliminates all signals of repeat-induced mutagenesis that extend beyond the motif boundary, and eliminates or dramatically shrinks the magnitude of mutagenesis within some motifs, contradicting previous reports. Mutagenesis not attributable to artifacts revealed several biological mechanisms. Polymerase slippage generates frequent indels within every variety of short tandem repeat motif, implicating slipped-strand structures. Interruption-correcting single nucleotide variants within short tandem repeats may originate from error-prone polymerases. Secondary-structure formation promotes single nucleotide variants within palindromic repeats and duplications within direct repeats. G-quadruplex motifs cause recurrent sequencing errors, whereas mutagenesis at Z-DNAs is conspicuously absent.

Evolution of the germline mutation rate across vertebrates

The germline mutation rate determines the pace of genome evolution and is an evolving parameter itself¹. However, little is known about what determines its evolution, as most studies of mutation rates have focused on single species with different methodologies². Here we quantify germline mutation rates across vertebrates by sequencing and comparing the high-coverage genomes of 151 parent-offspring trios from 68 species of mammals, fishes, birds and reptiles. We show that the per-generation mutation rate varies among species by a factor of 40, with mutation rates being higher for males than for females in mammals and birds, but not in reptiles and fishes. The generation time, age at maturity and species-level fecundity are the key life-history traits affecting this variation among species. Furthermore, species with higher long-term effective population sizes tend to have lower mutation rates per generation, providing support for the drift barrier hypothesis³. The exceptionally high yearly mutation rates of domesticated animals, which have been continually selected on fecundity traits including shorter generation times, further support the importance of generation time in the evolution of mutation rates. Overall, our comparative analysis of pedigree-based mutation rates provides ecological insights on the mutation rate evolution in vertebrates.

Personalized recurrence risk assessment following the birth of a child with a pathogenic de novo mutation

Following the diagnosis of a paediatric disorder caused by an apparently de novo mutation, a recurrence risk of 1-2% is frequently quoted due to the possibility of parental germline mosaicism; but for any specific couple, this figure is usually incorrect. We present a systematic approach to providing individualized recurrence risk. By combining locus-specific sequencing of multiple tissues to detect occult mosaicism with long-read sequencing to determine the parent-of-origin of the mutation, we show that we can stratify the majority of couples into one of seven discrete categories associated with substantially different risks to future offspring. Among 58 families with a single affected offspring (representing 59 de novo mutations in 49 genes), the recurrence risk for 35 (59%) was decreased below 0.1%, but increased owing to parental mixed mosaicism for 5 (9%)-that could be quantified in semen for paternal cases (recurrence risks of 5.6-12.1%). Implementation of this strategy offers the prospect of driving a major transformation in the practice of genetic counselling.

The evolution of personality disorders: A review of proposals

Personality disorders (PDs) are currently considered dysfunctions. However, personality differences are older than humanity and are ubiquitous in nature, from insects to higher primates. This suggests that a number of evolutionary mechanisms-other than dysfunctions-may be able to maintain stable behavioral variation in the gene pool. First of all, apparently maladaptive traits may actually improve fitness by enabling better survival or successful mating or reproduction, as exemplified by neuroticism, psychopathy, and narcissism. Furthermore, some PDs may harm important biological goals while facilitating others, or may be globally beneficial or detrimental depending on environmental circumstances or body condition. Alternatively, certain traits may form part of life history strategies: Coordinated suites of morphological, physiological and behavioral characters that optimize fitness through alternative routes and respond to selection as a whole. Still others may be vestigial adaptations that are no longer beneficial in present times. Finally, variation may be adaptative in and by itself, as it reduces competition for finite resources. These and other evolutionary mechanisms are reviewed and illustrated through human and non-human examples. Evolutionary theory is the best-substantiated explanatory framework across the life sciences, and may shed light on the question of why harmful personalities exist at all.

A method to build extended sequence context models of point mutations and indels

The mutation rate of a specific position in the human genome depends on the sequence context surrounding it. Modeling the mutation rate by estimating a rate for each possible k-mer, however, only works for small values of k since the data becomes too sparse for larger values of k. Here we propose a new method that solves this problem by grouping similar k-mers. We refer to the method as k-mer pattern partition and have implemented it in a software package called kmerPaPa. We use a large set of human de novo mutations to show that this new method leads to improved prediction of mutation rates and makes it possible to create models using wider sequence contexts than previous studies. As the first method of its kind, it does not only predict rates for point mutations but also insertions and deletions. We have additionally created a software package called Genovo that, given a k-mer pattern partition model, predicts the expected number of synonymous, missense, and other functional mutation types for each gene. Using this software, we show that the created mutation rate models increase the statistical power to detect genes containing disease-causing variants and to identify genes under strong selective constraint.

The mutational dynamics of short tandem repeats in large, multigenerational families

BACKGROUND

Short tandem repeats (STRs) compose approximately 3% of the genome, and mutations at STR loci have been linked to dozens of human diseases including amyotrophic lateral sclerosis, Friedreich ataxia, Huntington disease, and fragile X syndrome. Improving our understanding of these mutations would increase our knowledge of the mutational dynamics of the genome and may uncover additional loci that contribute to disease. To estimate the genome-wide pattern of mutations at STR loci, we analyze blood-derived whole-genome sequencing data for 544 individuals from 29 three-generation CEPH pedigrees. These pedigrees contain both sets of grandparents, the parents, and an average of 9 grandchildren per family.

RESULTS

We use HipSTR to identify de novo STR mutations in the 2nd generation of these pedigrees and require transmission to the third generation for validation. Analyzing approximately 1.6 million STR loci, we estimate the empirical de novo STR mutation rate to be 5.24 × 10^-5 mutations per locus per generation. Perfect repeats mutate about 2 × more often than imperfect repeats. De novo STRs are significantly enriched in Alu elements.

CONCLUSIONS

Approximately 30% of new STR mutations occur within Alu elements, which compose only 11% of the genome, but only 10% are found in LINE-1 insertions, which compose 17% of the genome. Phasing these mutations to the parent of origin shows that parental transmission biases vary among families. We estimate the average number of de novo genome-wide STR mutations per individual to be approximately 85, which is similar to the average number of observed de novo single nucleotide variants.

de novo variant calling identifies cancer mutation signatures in the 1000 Genomes Project

Detection of de novo variants (DNVs) is critical for studies of disease-related variation and mutation rates. To accelerate DNV calling, we developed a graphics processing units-based workflow. We applied our workflow to whole-genome sequencing data from three parent-child sequenced cohorts including the Simons Simplex Collection (SSC), Simons Foundation Powering Autism Research (SPARK), and the 1000 Genomes Project (1000G) that were sequenced using DNA from blood, saliva, and lymphoblastoid cell lines (LCLs), respectively. The SSC and SPARK DNV callsets were within expectations for number of DNVs, percent at CpG sites, phasing to the paternal chromosome of origin, and average allele balance. However, the 1000G DNV callset was not within expectations and contained excessive DNVs that are likely cell line artifacts. Mutation signature analysis revealed 30% of 1000G DNV signatures matched B-cell lymphoma. Furthermore, we found variants in DNA repair genes and at Clinvar pathogenic or likely-pathogenic sites and significant excess of protein-coding DNVs in IGLL5; a gene known to be involved in B-cell lymphomas. Our study provides a new rapid DNV caller for the field and elucidates important implications of using sequencing data from LCLs for reference building and disease-related projects.

Fertility in the aging male: a systematic review

The negative effect of advanced female age on fertility and offspring health is well understood. In comparison, much less is known about the implications of male age on fertility, with many studies showing conflicting results. Nevertheless, increasing evidence suggests that advanced paternal age has negative effects on sperm parameters, reproductive success, and offspring health. Herein, we summarize the current body of knowledge on this controversial topic, with the belief that this review will serve as a resource for the clinicians providing fertility counseling to couples with older male partners.

The impact of paternal age on new mutations and disease in the next generation

Advanced paternal age is associated with an increased risk of fathering children with genetic disorders and other adverse reproductive consequences. However, the mechanisms underlying this phenomenon remain largely unexplored. In this review, we focus on the impact of paternal age on de novo mutations that are an important contributor to genetic disease and can be studied both indirectly through large-scale sequencing studies and directly in the tissue in which they predominantly arise-the aging testis. We discuss the recent data that have helped establish the origins and frequency of de novo mutations, and highlight experimental evidence about the close link between new mutations, parental age, and genetic disease. We then focus on a small group of rare genetic conditions, the so-called "paternal age effect" disorders that show a strong association between paternal age and disease prevalence, and discuss the underlying mechanism ("selfish selection") and implications of this process in more detail. More broadly, understanding the causes and consequences of paternal age on genetic risk has important implications both for individual couples and for public health advice given that the average age of fatherhood is steadily increasing in many developed nations.

Paternal age and 13 psychiatric disorders in the offspring: a population-based cohort study of 7 million children in Taiwan

Although paternal age has been linked to certain psychiatric disorders, the nature of any causal relationship remains elusive. Here, we aimed to comprehensively assess the magnitude of a wide range of offspring's psychiatric risk conferred by paternal age, leveraging a pedigree inferred from covered-insurance relationship (accuracy >98%) in Taiwan's single-payer compulsory insurance program. We also examined whether there is an independent role of paternal age and explored the potential effect of parental age difference. A total cohort of 7,264,788 individuals born between 1980 and 2018 were included; 5,572,232 with sibling(s) were selected for sibling-comparison analyses and 1,368,942 and 1,044,420 children with information of paternal-grandparents and maternal-grandparents, respectively, were selected for multi-generation analyses. Using inpatient/outpatient claims data (1997-2018), we identified schizophrenia, autism, bipolar disorder (BPD), attention deficit-hyperactivity disorder (ADHD), major depressive disorder (MDD), eating disorder (ED), substance use disorder (SUD), mental retardation (MR), tic disorder, obsessive-compulsive disorder (OCD), anxiety, and somatoform disorder. We identified suicides using death certificates. Logistic regression analysis was used to estimate the paternal/maternal/grand-paternal age association with psychiatric risk in the offspring. The total cohort and sibling-comparison cohort resulted in similar estimates. Paternal age had a U-shaped relationship with offspring's MDD, ED, SUD, and anxiety. A very young maternal age (<20 years) was associated with markedly higher risk in offspring's SUD, MR, and suicide. Older paternal age (>25 years) was linearly associated with offspring's schizophrenia, autism, BPD, ADHD, MDD, ED, SUD, MR, OCD, anxiety, and suicide. Older grand-paternal age was linearly associated with offspring's schizophrenia, autism, ADHD, and MR. Dissimilar parental age was positively associated with offspring's ADHD, MDD, SUD, MR, anxiety, and suicide, and negatively associated with offspring's OCD. This comprehensive assessment provides solid evidence for the independent role of paternal age in psychiatric risk in the offspring and clarifies the significance of both early parenthood and delayed paternity.

De novo variants are a common cause of genetic hearing loss

PURPOSE

De novo variants (DNVs) are a well-recognized cause of genetic disorders. The contribution of DNVs to hearing loss (HL) is poorly characterized. We aimed to evaluate the rate of DNVs in HL-associated genes and assess their contribution to HL.

METHODS

Targeted genomic enrichment and massively parallel sequencing were used for molecular testing of all exons and flanking intronic sequences of known HL-associated genes, with no exclusions on the basis of type of HL or clinical features. Segregation analysis was performed, and previous reports of DNVs in PubMed and ClinVar were reviewed to characterize the rate, distribution, and spectrum of DNVs in HL.

RESULTS

DNVs were detected in 10% (24/238) of trios for whom segregation analysis was performed. Overall, DNVs were causative in at least ∼1% of probands for whom a genetic diagnosis was resolved, with marked variability based on inheritance mode and phenotype. DNVs of MITF were most common (21% of DNVs), followed by GATA3 (13%), STRC (13%), and ACTG1 (8%). Review of reported DNVs revealed gene-specific variability in contribution of DNV to the mutational spectrum of HL-associated genes.

CONCLUSION

DNVs are a relatively common cause of genetic HL and must be considered in all cases of sporadic HL.

Combining genetic constraint with predictions of alternative splicing to prioritize deleterious splicing in rare disease studies

Background

Despite numerous molecular and computational advances, roughly half of patients with a rare disease remain undiagnosed after exome or genome sequencing. A particularly challenging barrier to diagnosis is identifying variants that cause deleterious alternative splicing at intronic or exonic loci outside of canonical donor or acceptor splice sites.

Results

Several existing tools predict the likelihood that a genetic variant causes alternative splicing. We sought to extend such methods by developing a new metric that aids in discerning whether a genetic variant leads to deleterious alternative splicing. Our metric combines genetic variation in the Genome Aggregate Database with alternative splicing predictions from SpliceAI to compare observed and expected levels of splice-altering genetic variation. We infer genic regions with significantly less splice-altering variation than expected to be constrained. The resulting model of regional splicing constraint captures differential splicing constraint across gene and exon categories, and the most constrained genic regions are enriched for pathogenic splice-altering variants. Building from this model, we developed ConSpliceML. This ensemble machine learning approach combines regional splicing constraint with multiple per-nucleotide alternative splicing scores to guide the prediction of deleterious splicing variants in protein-coding genes. ConSpliceML more accurately distinguishes deleterious and benign splicing variants than state-of-the-art splicing prediction methods, especially in “cryptic” splicing regions beyond canonical donor or acceptor splice sites.

Conclusion

Integrating a model of genetic constraint with annotations from existing alternative splicing tools allows ConSpliceML to prioritize potentially deleterious splice-altering variants in studies of rare human diseases.

De novo mutations, genetic mosaicism and human disease

Mosaicism—the existence of genetically distinct populations of cells in a particular organism—is an important cause of genetic disease. Mosaicism can appear as de novo DNA mutations, epigenetic alterations of DNA, and chromosomal abnormalities. Neurodevelopmental or neuropsychiatric diseases, including autism—often arise by de novo mutations that usually not present in either of the parents. De novo mutations might occur as early as in the parental germline, during embryonic, fetal development, and/or post-natally, through ageing and life. Mutation timing could lead to mutation burden of less than heterozygosity to approaching homozygosity. Developmental timing of somatic mutation attainment will affect the mutation load and distribution throughout the body. In this review, we discuss the timing of de novo mutations, spanning from mutations in the germ lineage (all ages), to post-zygotic, embryonic, fetal, and post-natal events, through aging to death. These factors can determine the tissue specific distribution and load of de novo mutations, which can affect disease. The disease threshold burden of somatic de novo mutations of a particular gene in any tissue will be important to define.

Thousands of human mutation clusters are explained by short-range template switching

Variation within human genomes is unevenly distributed, and variants show spatial clustering. DNA replication-related template switching is a poorly known mutational mechanism capable of causing major chromosomal rearrangements as well as creating short inverted sequence copies that appear as local mutation clusters in sequence comparisons. In this study, haplotype-resolved genome assemblies representing 25 human populations and multinucleotide variants aggregated from 140,000 human sequencing experiments were reanalyzed. Local template switching could explain thousands of complex mutation clusters across the human genome, the loci segregating within and between populations. During the study, computational tools were developed for identification of template switch events using both short-read sequencing data and genotype data, and for genotyping candidate loci using short-read data. The characteristics of template-switch mutations complicate their detection, and widely used analysis pipelines for short-read sequencing data, normally capable of identifying single nucleotide changes, were found to miss template-switch mutations of tens of base pairs, potentially invalidating medical genetic studies searching for a causative allele behind genetic diseases. Combined with the massive sequencing data now available for humans, the novel tools described here enable building catalogs of affected loci and studying the cellular mechanisms behind template switching in both healthy organisms and disease.

A paternal bias in germline mutation is widespread in amniotes and can arise independently of cell divisions

In humans and other mammals, germline mutations are more likely to arise in fathers than in mothers. Although this sex bias has long been attributed to DNA replication errors in spermatogenesis, recent evidence from humans points to the importance of mutagenic processes that do not depend on cell division, calling into question our understanding of this basic phenomenon. Here, we infer the ratio of paternal-to-maternal mutations, α, in 42 species of amniotes, from putatively neutral substitution rates of sex chromosomes and autosomes. Despite marked differences in gametogenesis, physiologies and environments across species, fathers consistently contribute more mutations than mothers in all the species examined, including mammals, birds, and reptiles. In mammals, α is as high as 4 and correlates with generation times; in birds and snakes, α appears more stable around 2. These observations are consistent with a simple model, in which mutations accrue at equal rates in both sexes during early development and at a higher rate in the male germline after sexual differentiation, with a conserved paternal-to-maternal ratio across species. Thus, α may reflect the relative contributions of two or more developmental phases to total germline mutations, and is expected to depend on generation time even if mutations do not track cell divisions.

De novo Mutations in Domestic Cat are Consistent with an Effect of Reproductive Longevity on Both the Rate and Spectrum of Mutations

The mutation rate is a fundamental evolutionary parameter with direct and appreciable effects on the health and function of individuals. Here, we examine this important parameter in the domestic cat, a beloved companion animal as well as a valuable biomedical model. We estimate a mutation rate of 0.86 × 10-8 per bp per generation for the domestic cat (at an average parental age of 3.8 years). We find evidence for a significant paternal age effect, with more mutations transmitted by older sires. Our analyses suggest that the cat and the human have accrued similar numbers of mutations in the germline before reaching sexual maturity. The per-generation mutation rate in the cat is 28% lower than what has been observed in humans, but is consistent with the shorter generation time in the cat. Using a model of reproductive longevity, which takes into account differences in the reproductive age and time to sexual maturity, we are able to explain much of the difference in per-generation rates between species. We further apply our reproductive longevity model in a novel analysis of mutation spectra and find that the spectrum for the cat resembles the human mutation spectrum at a younger age of reproduction. Together, these results implicate changes in life-history as a driver of mutation rate evolution between species. As the first direct observation of the paternal age effect outside of rodents and primates, our results also suggest a phenomenon that may be universal among mammals.

Patterns and distribution of de novo mutations in multiplex Middle Eastern families

While de novo mutations (DNMs) are key to genetic diversity, they are also responsible for a high number of rare disorders. To date, no study has systematically examined the rate and distribution of DNMs in multiplex families in highly consanguineous populations. Leveraging WGS profiles of 645 individuals in 146 families, we implemented a combinatorial approach using 3 complementary tools for DNM discovery in 353 unique trio combinations. We found a total of 27,168 DNMs (median: 70 single-nucleotide and 6 insertion-deletions per individual). Phasing revealed around 80% of DNMs were paternal in origin. Notably, using whole-genome methylation data of spermatogonial stem cells, these DNMs were significantly more likely to occur at highly methylated CpGs (OR: 2.03; p value = 6.62 × 10⁻¹¹). We then examined the effects of consanguinity and ethnicity on DNMs, and found that consanguinity does not seem to correlate with DNM rate, and special attention has to be considered while measuring such a correlation. Additionally, we found that Middle-Eastern families with Arab ancestry had fewer DNMs than African families, although not significant (p value = 0.16). Finally, for families with diseased probands, we examined the difference in DNM counts and putative impact across affected and unaffected siblings, but did not find significant differences between disease groups, likely owing to the enrichment for recessive disorders in this part of the world, or the small sample size per clinical condition. This study serves as a reference for DNM discovery in multiplex families from the globally under-represented populations of the Middle-East.