Variation in genome-wide mutation rates within and between human families

The Evolutionary Interplay of Somatic and Germline Mutation Rates

Novel sequencing technologies are making it increasingly possible to measure the mutation rates of somatic cell lineages. Accurate germline mutation rate measurement technologies have also been available for a decade, making it possible to assess how this fundamental evolutionary parameter varies across the tree of life. Here, we review some classical theories about germline and somatic mutation rate evolution that were formulated using principles of population genetics and the biology of aging and cancer. We find that somatic mutation rate measurements, while still limited in phylogenetic diversity, seem consistent with the theory that selection to preserve the soma is proportional to life span. However, germline and somatic theories make conflicting predictions regarding which species should have the most accurate DNA repair. Resolving this conflict will require carefully measuring how mutation rates scale with time and cell division and achieving a better understanding of mutation rate pleiotropy among cell types.

Severity of effect considerations regarding the use of mutation as a toxicological endpoint for risk assessment: A report from the 8th International Workshop on Genotoxicity Testing (IWGT)

Exposure levels without appreciable human health risk may be determined by dividing a point of departure on a dose-response curve (e.g., benchmark dose) by a composite adjustment factor (AF). An "effect severity" AF (ESAF) is employed in some regulatory contexts. An ESAF of 10 may be incorporated in the derivation of a health-based guidance value (HBGV) when a "severe" toxicological endpoint, such as teratogenicity, irreversible reproductive effects, neurotoxicity, or cancer was observed in the reference study. Although mutation data have been used historically for hazard identification, this endpoint is suitable for quantitative dose-response modeling and risk assessment. As part of the 8th International Workshops on Genotoxicity Testing, a sub-group of the Quantitative Analysis Work Group (WG) explored how the concept of effect severity could be applied to mutation. To approach this question, the WG reviewed the prevailing regulatory guidance on how an ESAF is incorporated into risk assessments, evaluated current knowledge of associations between germline or somatic mutation and severe disease risk, and mined available data on the fraction of human germline mutations expected to cause severe disease. Based on this review and given that mutations are irreversible and some cause severe human disease, in regulatory settings where an ESAF is used, a majority of the WG recommends applying an ESAF value between 2 and 10 when deriving a HBGV from mutation data. This recommendation may need to be revisited in the future if direct measurement of disease-causing mutations by error-corrected next generation sequencing clarifies selection of ESAF values.

Massively parallel jumping assay decodes Alu retrotransposition activity

The human genome contains millions of retrotransposons, several of which could become active due to somatic mutations having phenotypic consequences, including disease. However, it is not thoroughly understood how nucleotide changes in retrotransposons affect their jumping activity. Here, we developed a novel massively parallel jumping assay (MPJA) that can test the jumping potential of thousands of transposons en masse. We generated nucleotide variant library of selected four Alu retrotransposons containing 165,087 different haplotypes and tested them for their jumping ability using MPJA. We found 66,821 unique jumping haplotypes, allowing us to pinpoint domains and variants vital for transposition. Mapping these variants to the Alu-RNA secondary structure revealed stem-loop features that contribute to jumping potential. Combined, our work provides a novel high-throughput assay that assesses the ability of retrotransposons to jump and identifies nucleotide changes that have the potential to reactivate them in the human genome.

Statistical methods for assessing the effects of de novo variants on birth defects

With the development of next-generation sequencing technology, de novo variants (DNVs) with deleterious effects can be identified and investigated for their effects on birth defects such as congenital heart disease (CHD). However, statistical power is still limited for such studies because of the small sample size due to the high cost of recruiting and sequencing samples and the low occurrence of DNVs. DNV analysis is further complicated by genetic heterogeneity across diseased individuals. Therefore, it is critical to jointly analyze DNVs with other types of genomic/biological information to improve statistical power to identify genes associated with birth defects. In this review, we discuss the general workflow, recent developments in statistical methods, and future directions for DNV analysis.

Estimating evolutionary and demographic parameters via ARG-derived IBD

Inference of demographic and evolutionary parameters from a sample of genome sequences often proceeds by first inferring identical-by-descent (IBD) genome segments. By exploiting efficient data encoding based on the ancestral recombination graph (ARG), we obtain three major advantages over current approaches: (i) no need to impose a length threshold on IBD segments, (ii) IBD can be defined without the hard-to-verify requirement of no recombination, and (iii) computation time can be reduced with little loss of statistical efficiency using only the IBD segments from a set of sequence pairs that scales linearly with sample size. We first demonstrate powerful inferences when true IBD information is available from simulated data. For IBD inferred from real data, we propose an approximate Bayesian computation inference algorithm and use it to show that poorly-inferred short IBD segments can improve estimation precision. We show estimation precision similar to a previously-published estimator despite a 4 000-fold reduction in data used for inference. Computational cost limits model complexity in our approach, but we are able to incorporate unknown nuisance parameters and model misspecification, still finding improved parameter inference.

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.

PMSeeker: A Scheme Based on the Greedy Algorithm and the Exhaustive Algorithm to Screen Low-Redundancy Marker Sets for Large-Scale Parentage Assignment with Full Parental Genotyping

Parentage assignment is a genetic test that utilizes genetic characteristics, such as molecular markers, to identify the parental relationships within populations, which, in commercial fish farming, are almost always large and where full information on potential parents is known. To accurately find the true parents, the genotypes of all loci in the parentage marker set (PMS) are required for each individual being tested. With the same accuracy, a PMS containing a smaller number of markers will undoubtedly save experimental costs. Thus, this study established a scheme to screen low-redundancy PMSs using the exhaustive algorithm and greedy algorithm. When screening PMSs, the greedy algorithm selects markers based on the parental dispersity index (PDI), a uniquely defined metric that outperforms the probability of exclusion (PE). With the conjunctive use of the two algorithms, non-redundant PMSs were found for more than 99.7% of solvable cases in three groups of random sample experiments in this study. Then, a low-redundancy PMS can be composed using two or more of these non-redundant PMSs. This scheme effectively reduces the number of markers in PMSs, thus conserving human and experimental resources and laying the groundwork for the widespread implementation of parentage assignment technology in economic species breeding.

The unusual gene architecture of polyubiquitin is created by dual-specific splice sites

BACKGROUND

The removal of introns occurs through the splicing of a 5' splice site (5'ss) with a 3' splice site (3'ss). These two elements are recognized by distinct components of the spliceosome. However, introns in higher eukaryotes contain many matches to the 5' and 3' splice-site motifs that are presumed not to be used.

RESULTS

Here, we find that many of these sites can be used. We also find occurrences of the AGGT motif that can function as either a 5'ss or a 3'ss-previously referred to as dual-specific splice sites (DSSs)-within introns. Analysis of the Sequence Read Archive reveals a 3.1-fold enrichment of DSSs relative to expectation, implying synergy between the ability to function as a 5'ss and 3'ss. Despite this suggested mechanistic advantage, DSSs are 2.7- and 4.7-fold underrepresented in annotated 5' and 3' splice sites. A curious exception is the polyubiquitin gene UBC, which contains a tandem array of DSSs that precisely delimit the boundary of each ubiquitin monomer. The resulting isoforms splice stochastically to include a variable number of ubiquitin monomers. We found no evidence of tissue-specific or feedback regulation but note the 8.4-fold enrichment of DSS-spliced introns in tandem repeat genes suggests a driving role in the evolution of genes like UBC.

CONCLUSIONS

We find an excess of unannotated splice sites and the utilization of DSSs in tandem repeats supports the role of splicing in gene evolution. These findings enhance our understanding of the diverse and complex nature of the splicing process.

A naturally occurring variant of MBD4 causes maternal germline hypermutation in primates

As part of an ongoing genome sequencing project at the Oregon National Primate Research Center, we identified a rhesus macaque with a rare homozygous frameshift mutation in the gene methyl-CpG binding domain 4, DNA glycosylase (MBD4). MBD4 is responsible for the repair of C > T deamination mutations at CpG dinucleotides and has been linked to somatic hypermutation and cancer predisposition in humans. We show here that MBD4-associated hypermutation also affects the germline: The six offspring of the MBD4-null dam have a fourfold to sixfold increase in de novo mutation burden. This excess burden was predominantly C > T mutations at CpG dinucleotides consistent with MBD4 loss of function in the dam. There was also a significant excess of C > T at CpA sites, indicating an important, unappreciated role for MBD4 to repair deamination in CpA contexts. The MBD4-null dam developed sustained eosinophilia later in life, but we saw no other signs of neoplastic processes associated with MBD4 loss of function in humans nor any obvious disease in the hypermutated offspring. This work provides the first evidence for a genetic factor causing hypermutation in the maternal germline of a mammal and adds to the very small list of naturally occurring variants known to modulate germline mutation rates in mammals.

Epigenetic aging of mammalian gametes

The process of aging refers to physiological changes that occur to an organism as time progresses and involves changes to DNA, proteins, metabolism, cells, and organs. Like the rest of the cells in the body, gametes age, and it is well established that there is a decline in reproductive capabilities in females and males with aging. One of the major pathways known to be involved in aging is epigenetic changes. The epigenome is the multitude of chemical modifications performed on DNA and chromatin that affect the ability of chromatin to be transcribed. In this review, we explore the effects of aging on female and male gametes with a focus on the epigenetic changes that occur in gametes throughout aging. Quality decline in oocytes occurs at a relatively early age. Epigenetic changes constitute an important part of oocyte aging. DNA methylation is reduced with age, along with reduced expression of DNA methyltransferases (DNMTs). Histone deacetylases (HDAC) expression is also reduced, and a loss of heterochromatin marks occurs with age. As a consequence of heterochromatin loss, retrotransposon expression is elevated, and aged oocytes suffer from DNA damage. In sperm, aging affects sperm number, motility and fecundity, and epigenetic changes may constitute a part of this process. 5 methyl-cytosine (5mC) methylation is elevated in sperm from aged men, but methylation on Long interspersed nuclear elements (LINE) elements is reduced. Di and trimethylation of histone 3 lysine 9 (H3K9me2/3) is reduced in sperm from aged men and trimethylation of histone 3 lysine 27 (H3K27me3) is elevated. The protamine makeup of sperm from aged men is also changed, with reduced protamine expression and a misbalanced ratio between protamine proteins protamine P1 and protamine P2. The study of epigenetic reproductive aging is recently gaining interest. The current status of the field suggests that many aspects of gamete epigenetic aging are still open for investigation. The clinical applications of these investigations have far-reaching consequences for fertility and sociological human behavior.

Quartet DNA reference materials and datasets for comprehensively evaluating germline variant calling performance

BACKGROUND

Genomic DNA reference materials are widely recognized as essential for ensuring data quality in omics research. However, relying solely on reference datasets to evaluate the accuracy of variant calling results is incomplete, as they are limited to benchmark regions. Therefore, it is important to develop DNA reference materials that enable the assessment of variant detection performance across the entire genome.

RESULTS

We established a DNA reference material suite from four immortalized cell lines derived from a family of parents and monozygotic twins. Comprehensive reference datasets of 4.2 million small variants and 15,000 structural variants were integrated and certified for evaluating the reliability of germline variant calls inside the benchmark regions. Importantly, the genetic built-in-truth of the Quartet family design enables estimation of the precision of variant calls outside the benchmark regions. Using the Quartet reference materials along with study samples, batch effects are objectively monitored and alleviated by training a machine learning model with the Quartet reference datasets to remove potential artifact calls. Moreover, the matched RNA and protein reference materials and datasets from the Quartet project enables cross-omics validation of variant calls from multiomics data.

CONCLUSIONS

The Quartet DNA reference materials and reference datasets provide a unique resource for objectively assessing the quality of germline variant calls throughout the whole-genome regions and improving the reliability of large-scale genomic profiling.

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.

The Role of Genetics in the Management of Heart Failure Patients

Over the last decades, the relevance of genetics in cardiovascular diseases has expanded, especially in the context of cardiomyopathies. Its relevance extends to the management of patients diagnosed with heart failure (HF), given its capacity to provide invaluable insights into the etiology of cardiomyopathies and identify individuals at a heightened risk of poor outcomes. Notably, the identification of an etiological genetic variant necessitates a comprehensive evaluation of the family lineage of the affected patients. In the future, these genetic variants hold potential as therapeutic targets with the capability to modify gene expression. In this complex setting, collaboration among cardiologists, specifically those specializing in cardiomyopathies and HF, and geneticists becomes paramount to improving individual and family health outcomes, as well as therapeutic clinical results. This review is intended to offer geneticists and cardiologists an updated perspective on the value of genetic research in HF and its implications in clinical practice.

Rare variant modifier analysis identifies variants in SEC24D associated with orofacial cleft subtypes

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p = 6.86 [Formula: see text] 10-7). Of the 15 variants within SEC24D, 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p = 0.0009), and created binding sites for 23 TFs, including Pax6 (p = 6.12 [Formula: see text] 10-5) and Pax9 (p = 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study indicates that rare genetic variation may contribute to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC.

Genomic inference of a severe human bottleneck during the Early to Middle Pleistocene transition

Population size history is essential for studying human evolution. However, ancient population size history during the Pleistocene is notoriously difficult to unravel. In this study, we developed a fast infinitesimal time coalescent process (FitCoal) to circumvent this difficulty and calculated the composite likelihood for present-day human genomic sequences of 3154 individuals. Results showed that human ancestors went through a severe population bottleneck with about 1280 breeding individuals between around 930,000 and 813,000 years ago. The bottleneck lasted for about 117,000 years and brought human ancestors close to extinction. This bottleneck is congruent with a substantial chronological gap in the available African and Eurasian fossil record. Our results provide new insights into our ancestry and suggest a coincident speciation event.

Extensive variation in germline de novo mutations in Poecilia reticulata

The rate of germline mutation is fundamental to evolutionary processes, as it generates the variation upon which selection acts. The guppy, Poecilia reticulata, is a model of rapid adaptation, however the relative contribution of standing genetic variation versus de novo mutation (DNM) to evolution in this species remains unclear. Here, we use pedigree-based approaches to quantify and characterize germline DNMs in three large guppy families. Our results suggest germline mutation rate in the guppy varies substantially across individuals and families. Most DNMs are shared across multiple siblings, suggesting they arose during early embryonic development. DNMs are randomly distributed throughout the genome, and male-biased mutation rate is low, as would be expected from the short guppy generation time. Overall, our study shows remarkable variation in germline mutation rate and provides insights into rapid evolution of guppies.

The contributions of rare inherited and polygenic risk to ASD in multiplex families

Autism spectrum disorder (ASD) has a complex genetic architecture involving contributions from both de novo and inherited variation. Few studies have been designed to address the role of rare inherited variation or its interaction with common polygenic risk in ASD. Here, we performed whole-genome sequencing of the largest cohort of multiplex families to date, consisting of 4,551 individuals in 1,004 families having two or more autistic children. Using this study design, we identify seven previously unrecognized ASD risk genes supported by a majority of rare inherited variants, finding support for a total of 74 genes in our cohort and a total of 152 genes after combined analysis with other studies. Autistic children from multiplex families demonstrate an increased burden of rare inherited protein-truncating variants in known ASD risk genes. We also find that ASD polygenic score (PGS) is overtransmitted from nonautistic parents to autistic children who also harbor rare inherited variants, consistent with combinatorial effects in the offspring, which may explain the reduced penetrance of these rare variants in parents. We also observe that in addition to social dysfunction, language delay is associated with ASD PGS overtransmission. These results are consistent with an additive complex genetic risk architecture of ASD involving rare and common variation and further suggest that language delay is a core biological feature of ASD.

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.

The origins and functional effects of postzygotic mutations throughout the human life span

Postzygotic mutations (PZMs) begin to accrue in the human genome immediately after fertilization, but how and when PZMs affect development and lifetime health remain unclear. To study the origins and functional consequences of PZMs, we generated a multitissue atlas of PZMs spanning 54 tissue and cell types from 948 donors. Nearly half the variation in mutation burden among tissue samples can be explained by measured technical and biological effects, and 9% can be attributed to donor-specific effects. Through phylogenetic reconstruction of PZMs, we found that their type and predicted functional impact vary during prenatal development, across tissues, and through the germ cell life cycle. Thus, methods for interpreting effects across the body and the life span are needed to fully understand the consequences of genetic variants.

Increased male investment in sperm competition results in reduced maintenance of gametes

Male animals often show higher mutation rates than their female conspecifics. A hypothesis for this male bias is that competition over fertilization of female gametes leads to increased male investment into reproduction at the expense of maintenance and repair, resulting in a trade-off between male success in sperm competition and offspring quality. Here, we provide evidence for this hypothesis by harnessing the power of experimental evolution to study effects of sexual selection on the male germline in the seed beetle Callosobruchus maculatus. We first show that 50 generations of evolution under strong sexual selection, coupled with experimental removal of natural selection, resulted in males that are more successful in sperm competition. We then show that these males produce progeny of lower quality if engaging in sociosexual interactions prior to being challenged to surveil and repair experimentally induced damage in their germline and that the presence of male competitors alone can be enough to elicit this response. We identify 18 candidate genes that showed differential expression in response to the induced germline damage, with several of these previously implicated in processes associated with DNA repair and cellular maintenance. These genes also showed significant expression changes across sociosexual treatments of fathers and predicted the reduction in quality of their offspring, with expression of one gene also being strongly correlated to male sperm competition success. Sex differences in expression of the same 18 genes indicate a substantially higher female investment in germline maintenance. While more work is needed to detail the exact molecular underpinnings of our results, our findings provide rare experimental evidence for a trade-off between male success in sperm competition and germline maintenance. This suggests that sex differences in the relative strengths of sexual and natural selection are causally linked to male mutation bias. The tenet advocated here, that the allocation decisions of an individual can affect plasticity of its germline and the resulting genetic quality of subsequent generations, has several interesting implications for mate choice processes.

Rare genetic variants in SEC24D modify orofacial cleft phenotypes

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p=6.86×10-7). Of the 15 variants within SEC24D, 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p=0.0009), and created binding sites for 23 TFs, including Pax6 (p=6.12×10-5) and Pax9 (p= 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study demonstrates that rare genetic variation contributes to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC.

Mutation Rate and Spectrum of the Silkworm in Normal and Temperature Stress Conditions

Mutation rate is a crucial parameter in evolutionary genetics. However, the mutation rate of most species as well as the extent to which the environment can alter the genome of multicellular organisms remain poorly understood. Here, we used parents–progeny sequencing to investigate the mutation rate and spectrum of the domestic silkworm (Bombyx mori) among normal and two temperature stress conditions (32 °C and 0 °C). The rate of single-nucleotide mutations in the normal temperature rearing condition was 0.41 × 10−8 (95% confidence interval, 0.33 × 10−8–0.49 × 10−8) per site per generation, which was up to 1.5-fold higher than in four previously studied insects. Moreover, the mutation rates of the silkworm under the stresses are significantly higher than in normal conditions. Furthermore, the mutation rate varies less in gene regions under normal and temperature stresses. Together, these findings expand the known diversity of the mutation rate among eukaryotes but also have implications for evolutionary analysis that assumes a constant mutation rate among species and environments.

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.

Anatomical variations of the flexor carpi ulnaris in the fetal period

Introduction: The Flexor Carpi Ulnaris (FCU) is a part of the palmar the forearm muscle group and one of the most important muscles for upper limb functioning - is responsible for flexion and adduc­tion of the hand at the radio-carpal joint. There are clinically significant but rare anatomical variations of FCU. The variability of the FCU has not been described up to now, and no typology of the muscle based on its more variable terminal attachment has been created. Aim of the study: Determination of FCU muscle typology based on available fetal material. Material and methods: A total of 114 human fetuses (53 female, 61 male) between 117 and 197 days of fetal life were eligible for the study. Preparations were carried out using classical anatomical techniques based on a previously published procedure. Thanks to that significant anthropometric landmarks were vis­ible for the gathering of metric measurements. Metric measurements were taken and statistically analysed using R-Project software. Results: A new typology was created based on variable muscle insertions. Additionally, the presence of an atypically located, additional, separated muscle belly was described. A comparison of measurements of the left upper limb in relation to the right upper limb showed significant differences for forearm length to the anthropometric point of the stylion radiale, limb length, total FCU length and FCU length which means that the left limb is longer than the right limb. A comparison of FCU insertion types between left and right upper limb showed there’s no significant difference between counts of each type. Conclusion: The FCU is a muscle that is easy to palpate and may therefore act as a topographical marker for healthcare professionals. Knowledge of its variability is not only of theoretical importance but also has clinical sig­nificance. The current publication demonstrates presence of variability in FCU terminal attachment. Certainly, this topic requires further research and continued work on a detailed understanding of forearm anatomy in the fetal period.

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.

Recombination and mutation shape variations in the major histocompatibility complex

The major histocompatibility complex (MHC) is closely associated with numerous diseases, but its high degree of polymorphism complicates the discovery of disease-associated variants. In principle, recombination and de novo mutations are two critical factors responsible for MHC polymorphisms. However, direct evidence for this hypothesis is lacking. Here, we report the generation of fine-scale MHC recombination and de novo mutation maps of ∼5 Mb by deep sequencing (> 100×) of the MHC genome for 17 MHC recombination and 30 non-recombination Han Chinese families (a total of 190 individuals). Recombination hotspots and Han-specific breakpoints are located in close proximity at haplotype block boundaries. The average MHC de novo mutation rate is higher than the genome-wide de novo mutation rate, particularly in MHC recombinant individuals. Notably, mutation and recombination generated polymorphisms are located within and outside linkage disequilibrium regions of the MHC, respectively, and evolution of the MHC locus was mainly controlled by positive selection. These findings provide insights on the evolutionary causes of the MHC diversity and may facilitate the identification of disease-associated genetic variants.

Sequencing individual genomes with recurrent genomic disorder deletions: an approach to characterize genes for autosomal recessive rare disease traits

BACKGROUND

In medical genetics, discovery and characterization of disease trait contributory genes and alleles depends on genetic reasoning, study design, and patient ascertainment; we suggest a segmental haploid genetics approach to enhance gene discovery and molecular diagnostics.

METHODS

We constructed a genome-wide map for nonallelic homologous recombination (NAHR)-mediated recurrent genomic deletions and used this map to estimate population frequencies of NAHR deletions based on large-scale population cohorts and region-specific studies. We calculated recessive disease carrier burden using high-quality pathogenic or likely pathogenic variants from ClinVar and gnomAD. We developed a NIRD (NAHR deletion Impact to Recessive Disease) score for recessive disorders by quantifying the contribution of NAHR deletion to the overall allele load that enumerated all pairwise combinations of disease-causing alleles; we used a Punnett square approach based on an assumption of random mating. Literature mining was conducted to identify all reported patients with defects in a gene with a high NIRD score; meta-analysis was performed on these patients to estimate the representation of NAHR deletions in recessive traits from contemporary human genomics studies. Retrospective analyses of extant clinical exome sequencing (cES) were performed for novel rare recessive disease trait gene and allele discovery from individuals with NAHR deletions.

RESULTS

We present novel genomic insights regarding the genome-wide impact of NAHR recurrent segmental variants on recessive disease burden; we demonstrate the utility of NAHR recurrent deletions to enhance discovery in the challenging context of autosomal recessive (AR) traits and biallelic variation. Computational results demonstrate new mutations mediated by NAHR, involving recurrent deletions at 30 genomic regions, likely drive recessive disease burden for over 74% of loci within these segmental deletions or at least 2% of loci genome-wide. Meta-analyses on 170 literature-reported patients implicate that NAHR deletions are depleted from the ascertained pool of AR trait alleles. Exome reanalysis of personal genomes from subjects harboring recurrent deletions uncovered new disease-contributing variants in genes including COX10, ERCC6, PRRT2, and OTUD7A.

CONCLUSIONS

Our results demonstrate that genomic sequencing of personal genomes with NAHR deletions could dramatically improve allele and gene discovery and enhance clinical molecular diagnosis. Moreover, results suggest NAHR events could potentially enable human haploid genetic screens as an approach to experimental inquiry into disease biology.

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.

A genotype likelihood function for DNA mixtures

The recent advent of genetic genealogy has brought about a renewed interest in genome-scale forensic analyses, of which kinship estimation is a critical component. Most genomic kinship estimators consider SNPs (single nucleotide polymorphisms), often leveraging the co-inheritance of shared alleles to inform their analyses. While current estimators cannot directly evaluate mixed samples, there exist well-established SNP-based kinship estimators tailored to considering challenged samples, including low-pass whole genome sequencing. As an example, several studies have shown remarkable success in imputing genotype posterior probabilities in low template samples when linked sites are considered. Critical to these approaches is the ability to account for genotype uncertainty; the lack of an expression for a genotype likelihood in imbalanced mixtures has prevented direct application. This work develops such an expression. The formulation is fully compatible with genotype imputation software, suggesting a genomic pipeline that estimates genotype likelihoods, performs imputation, and then estimates kinship when the sample is a mixture. Further, when framed as an imbalanced mixture, the problem of mixture deconvolution is reducible to the problem of genotyping mixed samples. Herein, the ability to genotype two-person mixtures is assessed through example and in silico settings. While certain mixture scenarios and classes of sites are inherently inseparable, simulations of read depths between 60 and 190 appear to produce likelihoods of sufficient magnitude to deconvolve two-person mixtures whenever the mixture fraction is moderately imbalanced. The described approach and results suggest a path forward for estimating the kinship coefficient (and similar inferences on relatedness) when the sample is a mixture.

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.

No evidence of increased mutations in the germline of a group of British nuclear test veterans

The potential germline effects of radiation exposure to military veterans present at British nuclear tests in Australia and the South Pacific is of considerable interest. We analyzed germline mutations in 60 families of UK military personnel comprising 30 control and 30 nuclear test veterans (NTV). Using whole-genome sequencing we studied the frequency and spectra of de novo mutations to investigate the transgenerational effect of veterans' (potential) exposure to radiation at nuclear bomb test sites. We find no elevation in total de novo single nucleotide variants, small insertion-deletions, structural variants or clustered mutations among the offspring of nuclear test veterans compared to those of control personnel. We did observe an elevated occurrence of single base substitution mutations within mutation signature SBS16, due to a subset of NTV offspring. The relevance of this elevation to potential exposure of veteran fathers and, future health risks, require further investigation. Overall, we find no evidence of increased mutations in the germline of a group of British nuclear test veterans. ISRCTN Registry 17461668.

Germline De Novo Mutations as a Cause of Childhood Cancer

Germline de novo mutations (DNMs) represent one of the important topics that need extensive attention from epidemiologists, geneticists, and other relevant stakeholders. Advances in next-generation sequencing technologies allowed examination of parent-offspring trios to ascertain the frequency of germline DNMs. Many epidemiological risk factors for childhood cancer are indicative of DNMs as a mechanism. The aim of this review was to give an overview of germline DNMs, their causes in general, and to discuss their relation to childhood cancer risk. In addition, we highlighted existing gaps in knowledge in many topics of germline DNMs in childhood cancer that need exploration and collaborative efforts.

Cancer Risk and Mutational Patterns Following Organ Transplantation

The rapid development of medical technology and widespread application of immunosuppressive drugs have improved the success rate of organ transplantation significantly. However, the use of immunosuppressive agents increases the frequency of malignancy greatly. With the prospect of “precision medicine” for tumors and development of next-generation sequencing technology, more attention has been paid to the application of high-throughput sequencing technology in clinical oncology research, which is mainly applied to the early diagnosis of tumors and analysis of tumor-related genes. All generations of cancers carry somatic mutations, meanwhile, significant differences were observed in mutational signatures across tumors. Systematic sequencing of cancer genomes from patients after organ transplantation can reveal DNA damage and repair processes in exposed cancer cells and their precursors. In this review, we summarize the application of high-throughput sequencing and organoids in the field of organ transplantation, the mutational patterns of cancer genomes, and propose a new research strategy for understanding the mechanism of cancer following organ transplantation.

[Exploring the association between de novo mutations and non-syndromic cleft lip with or without palate based on whole exome sequencing of case-parent trios]

OBJECTIVE

To explore the association between de novo mutations (DNM) and non-syndromic cleft lip with or without palate (NSCL/P) using case-parent trio design.

METHODS

Whole-exome sequencing was conducted for twenty-two NSCL/P trios and Genome Analysis ToolKit (GATK) was used to identify DNM by comparing the alleles of the cases and their parents. Information of predictable functions was annotated to the locus with SnpEff. Enrichment analysis for DNM was conducted to test the difference between the actual number and the expected number of DNM, and to explore whether there were genes with more DNM than expected. NSCL/P-related genes indicated by previous studies with solid evidence were selected by literature reviewing. Protein-protein interactions analysis was conducted among the genes with protein-altering DNM and NSCL/P-related genes. R package "denovolyzeR" was used for the enrichment analysis (Bonferroni correction: P=0.05/n, n is the number of genes in the whole genome range). Protein-protein interactions among genes with DNM and genes with solid evidence on the risk factors of NSCL/P were predicted depending on the information provided by STRING database.

RESULTS

A total of 339 908 SNPs were qualified for the subsequent analysis after quality control. The number of high confident DNM identified by GATK was 345. Among those DNM, forty-four DNM were missense mutations, one DNM was nonsense mutation, two DNM were splicing site mutations, twenty DNM were synonymous mutations and others were located in intron or intergenic regions. The results of enrichment analysis showed that the number of protein-altering DNM on the exome regions was larger than expected (P < 0.05), and five genes (KRTCAP2, HMCN2, ANKRD36C, ADGRL2 and DIPK2A) had more DNM than expected (P < 0.05/(2×19 618)). Protein-protein interaction analysis was conducted among forty-six genes with protein-altering DNM and thirteen genes associated with NSCL/P selected by literature reviewing. Six pairs of interactions occurred between the genes with DNM and known NSCL/P-related genes. The score measuring the confidence level of the predicted interaction between RGPD4 and SUMO1 was 0.868, which was higher than the scores for other pairs of genes.

CONCLUSION

Our study provided novel insights into the development of NSCL/P and demonstrated that functional analyses of genes carrying DNM were warranted to understand the genetic architecture of complex diseases.

Decoding signal transducer and activator of transcription 1 across various cancers through data mining and integrative analysis

OBJECTIVES

Using different online available databases and Bioinformatics tools, we extensively studied the role STAT1 across different cancers.

METHODS

STAT1 mRNA, protein expression, and promoter methylation were analyzed and validated using UALCAN, GENT2, Human Protein Atlas (HPA), and MEXPRESS. Furthermore, the potential prognostic values were evaluated through KM plotter. Then, cBioPortal was utilized to examine the STAT1-related genetic mutations, while pathway enrichment analysis was performed using DAVID. To identify STAT1 targeted microRNAs (miRNAs) and transcription factors (TFs) we used Enricher. Moreover, a correlational analysis between STAT1 expression tumor purity and CD8+ T immune cells and a gene-drug interaction network analysis was performed using TIMER, CTD, and Cytoscape.

RESULTS

In 23 major human cancers, STAT1 expression was notably up-regulated relative to corresponding controls. As well, the elevated expression of STAT1 was exclusively found to be associated with the reduced overall survival (OS) of Esophageal Carcinoma (ESCA), Kidney Renal Clear Cell Carcinoma (KIRC), and Lung adenocarcinoma (LUAD) patients. This implies that STAT1 plays a significant role in the development and progression of these three cancers. Further pathway analysis indicated that STAT1 enriched genes were involved in six critical pathways, while a few interesting correlations were also documented between STAT1 expression and promoter methylation level, tumor purity, CD8+ T immune cells infiltration, and genetic alteration. In addition, we have also predicted a few miRNAs, TFs, and chemotherapeutic drugs that could regulate the STAT1 expression.

CONCLUSION

The current study revealed the shared oncogenic, diagnostic, and prognostic role of STAT1 in ESCA, KIRC, and LUAD.

Genetic and chemotherapeutic influences on germline hypermutation

Mutations in the germline generates all evolutionary genetic variation and is a cause of genetic disease. Parental age is the primary determinant of the number of new germline mutations in an individual's genome1,2. Here we analysed the genome-wide sequences of 21,879 families with rare genetic diseases and identified 12 individuals with a hypermutated genome with between two and seven times more de novo single-nucleotide variants than expected. In most families (9 out of 12), the excess mutations came from the father. Two families had genetic drivers of germline hypermutation, with fathers carrying damaging genetic variation in DNA-repair genes. For five of the families, paternal exposure to chemotherapeutic agents before conception was probably a key driver of hypermutation. Our results suggest that the germline is well protected from mutagenic effects, hypermutation is rare, the number of excess mutations is relatively modest and most individuals with a hypermutated genome will not have a genetic disease.

Familial long-read sequencing increases yield of de novo mutations

Studies of de novo mutation (DNM) have typically excluded some of the most repetitive and complex regions of the genome because these regions cannot be unambiguously mapped with short-read sequencing data. To better understand the genome-wide pattern of DNM, we generated long-read sequence data from an autism parent-child quad with an affected female where no pathogenic variant had been discovered in short-read Illumina sequence data. We deeply sequenced all four individuals by using three sequencing platforms (Illumina, Oxford Nanopore, and Pacific Biosciences) and three complementary technologies (Strand-seq, optical mapping, and 10X Genomics). Using long-read sequencing, we initially discovered and validated 171 DNMs across two children-a 20% increase in the number of de novo single-nucleotide variants (SNVs) and indels when compared to short-read callsets. The number of DNMs further increased by 5% when considering a more complete human reference (T2T-CHM13) because of the recovery of events in regions absent from GRCh38 (e.g., three DNMs in heterochromatic satellites). In total, we validated 195 de novo germline mutations and 23 potential post-zygotic mosaic mutations across both children; the overall true substitution rate based on this integrated callset is at least 1.41 × 10-8 substitutions per nucleotide per generation. We also identified six de novo insertions and deletions in tandem repeats, two of which represent structural variants. We demonstrate that long-read sequencing and assembly, especially when combined with a more complete reference genome, increases the number of DNMs by >25% compared to previous studies, providing a more complete catalog of DNM compared to short-read data alone.

Neurofibromatosis type 1 families with first-degree relatives harbouring distinct NF1 pathogenic variants. Genetic counselling and familial diagnosis: what should be offered?

Neurofibromatosis type 1 (NF1) is an autosomal dominant disorder caused by pathogenic variants in NF1 Recently, NF1 testing has been included as a clinical criterion for NF1 diagnosis. Additionally, preconception genetic counselling in patients with NF1 focuses on a 50% risk of transmitting the familial variant as the risk of having a sporadic NF1 is considered the same as the general population.

METHODS

829 individuals, 583 NF1 sporadic cases and 246 patients with NF1 with documented family history, underwent genetic testing for NF1. Genotyping and segregation analysis of NF1 familial variants was determined by microsatellite analysis and NF1 sequencing.

RESULTS

The mutational analysis of NF1 in 154 families with two or more affected cases studied showed the co-occurrence of two different NF1 germline pathogenic variants in four families. The estimated mutation rate in those families was 3.89×10^-3, 20 times higher than the NF1 mutation rate (~2×10^-4) (p=0.0008). Furthermore, the co-occurrence of two different NF1 germline pathogenic variants in these families was 1:39, 60 times the frequency of sporadic NF1 (1:2500) (p=0.003). In all cases, the de novo NF1 pathogenic variant was present in a descendant of an affected male. In two cases, variants were detected in the inherited paternal wild-type allele.

CONCLUSIONS

Our results, together with previous cases reported, suggest that the offspring of male patients with NF1 could have an increased risk of experiencing de novo NF1 pathogenic variants. This observation, if confirmed in additional cohorts, could have relevant implications for NF1 genetic counselling, family planning and NF1 genetic testing.

Computational Genomics in the Era of Precision Medicine: Applications to Variant Analysis and Gene Therapy

Rapid methodological advances in statistical and computational genomics have enabled researchers to better identify and interpret both rare and common variants responsible for complex human diseases. As we continue to see an expansion of these advances in the field, it is now imperative for researchers to understand the resources and methodologies available for various data types and study designs. In this review, we provide an overview of recent methods for identifying rare and common variants and understanding their roles in disease etiology. Additionally, we discuss the strategy, challenge, and promise of gene therapy. As computational and statistical approaches continue to improve, we will have an opportunity to translate human genetic findings into personalized health care.

De novo mutation rates at the single-mutation resolution in a human HBB gene-region associated with adaptation and genetic disease

Although it is known that the mutation rate varies across the genome, previous estimates were based on averaging across various numbers of positions. Here, we describe a method to measure the origination rates of target mutations at target base positions and apply it to a 6-bp region in the human hemoglobin subunit beta (HBB) gene and to the identical, paralogous hemoglobin subunit delta (HBD) region in sperm cells from both African and European donors. The HBB region of interest (ROI) includes the site of the hemoglobin S (HbS) mutation, which protects against malaria, is common in Africa, and has served as a classic example of adaptation by random mutation and natural selection. We found a significant correspondence between de novo mutation rates and past observations of alleles in carriers, showing that mutation rates vary substantially in a mutation-specific manner that contributes to the site frequency spectrum. We also found that the overall point mutation rate is significantly higher in Africans than in Europeans in the HBB region studied. Finally, the rate of the 20A→T mutation, called the “HbS mutation” when it appears in HBB, is significantly higher than expected from the genome-wide average for this mutation type. Nine instances were observed in the African HBB ROI, where it is of adaptive significance, representing at least three independent originations; no instances were observed elsewhere. Further studies will be needed to examine mutation rates at the single-mutation resolution across these and other loci and organisms and to uncover the molecular mechanisms responsible.

The mutationathon highlights the importance of reaching standardization in estimates of pedigree-based germline mutation rates

In the past decade, several studies have estimated the human per-generation germline mutation rate using large pedigrees. More recently, estimates for various nonhuman species have been published. However, methodological differences among studies in detecting germline mutations and estimating mutation rates make direct comparisons difficult. Here, we describe the many different steps involved in estimating pedigree-based mutation rates, including sampling, sequencing, mapping, variant calling, filtering, and appropriately accounting for false-positive and false-negative rates. For each step, we review the different methods and parameter choices that have been used in the recent literature. Additionally, we present the results from a ‘Mutationathon,’ a competition organized among five research labs to compare germline mutation rate estimates for a single pedigree of rhesus macaques. We report almost a twofold variation in the final estimated rate among groups using different post-alignment processing, calling, and filtering criteria, and provide details into the sources of variation across studies. Though the difference among estimates is not statistically significant, this discrepancy emphasizes the need for standardized methods in mutation rate estimations and the difficulty in comparing rates from different studies. Finally, this work aims to provide guidelines for computational and statistical benchmarks for future studies interested in identifying germline mutations from pedigrees.

Spontaneous mutation rate estimates for the principal malaria vectors Anopheles coluzzii and Anopheles stephensi

Using high-depth whole genome sequencing of F0 mating pairs and multiple individual F1 offspring, we estimated the nuclear mutation rate per generation in the malaria vectors Anopheles coluzzii and Anopheles stephensi by detecting de novo genetic mutations. A purpose-built computer program was employed to filter actual mutations from a deep background of superficially similar artifacts resulting from read misalignment. Performance of filtering parameters was determined using software-simulated mutations, and the resulting estimate of false negative rate was used to correct final mutation rate estimates. Spontaneous mutation rates by base substitution were estimated at 1.00 × 10⁻⁹ (95% confidence interval, 2.06 × 10⁻¹⁰—2.91 × 10⁻⁹) and 1.36 × 10⁻⁹ (95% confidence interval, 4.42 × 10⁻¹⁰—3.18 × 10⁻⁹) per site per generation in A. coluzzii and A. stephensi respectively. Although similar studies have been performed on other insect species including dipterans, this is the first study to empirically measure mutation rates in the important genus Anopheles, and thus provides an estimate of µ that will be of utility for comparative evolutionary genomics, as well as for population genetic analysis of malaria vector mosquito species.

Resolving inherited and de novo germline predisposing sequence variants by means of whole exome trio analyses in childhood hematological malignancies

Molecular screening tools have significantly eased the assessment of potential germline susceptibility factors that may underlie the development of pediatric malignancies. Most of the hitherto published studies utilize the comparative analyses of the respective patients' germline and tumor tissues for this purpose. Since this approach is not able to discriminate between de novo and inherited sequence variants, we performed whole exome trio analyses in a consecutive series of 131 children with various forms of hematologic malignancies and their parents. In total, we identified 458 de novo variants with a range from zero to 28 (median value = 3) per patient, although most of them (58%) had only up to three per exome. Overall, we identified bona fide cancer predisposing alterations in five of the investigated 131 (3.8%) patients. Three of them had de novo pathogenic lesions in the SOS1, PTPN11 and TP53 genes and two of them parentally inherited ones in the STK11 and PMS2 genes that are specific for a Peutz-Jeghers and a constitutional mismatch repair deficiency (CMMRD) syndrome, respectively. Notwithstanding that we did not identify a disease-specific alteration in the two cases with the highest number of de novo variants, one of them developed two almost synchronous malignancies: a myelodysplastic syndrome and successively within two months a cerebral astrocytoma. Moreover, we also found that the rate of de novo sequence variants in the offspring increased especially with the age of the father, but less so with that of the mother. We therefore conclude that trio analyses deliver an immediate overview about the inheritance pattern of the entire spectrum of sequence variants, which not only helps to securely identify the de novo or inherited nature of genuinely disease-related lesions, but also of all other less obvious variants that in one or the other way may eventually advance our understanding of the disease process.