Analysis of genetic inheritance in a family quartet by whole-genome sequencing

Muons, mutations, and planetary shielding

Life on earth is protected from astrophysical cosmic rays by the heliospheric magnetic and slowly varying geomagnetic fields, and by collisions with oxygen and nitrogen molecules in the atmosphere. The collisions generate showers of particles of lesser energy; only muons, a charged particle with a mass between that of an electron and a proton, can reach earth's surface in substantial quantities. Muons are easily detected, used to image interior spaces of pyramids, and known to limit the stability of qubits in quantum computing; yet, despite their charge, average energy of 4 GeV and ionizing properties, muons are not considered to affect chemical reactions or biology. In this Perspective the potential damaging effects of muons on DNA, and hence the repercussions for evolution and disease, are examined. It is argued here that the effect of muons on life through DNA mutations should be considered when investigating the protection provided by the magnetic environment and atmosphere from cosmic rays on earth and exoplanets.

Duet: SNP-assisted structural variant calling and phasing using Oxford nanopore sequencing

BACKGROUND

Whole genome sequencing using the long-read Oxford Nanopore Technologies (ONT) MinION sequencer provides a cost-effective option for structural variant (SV) detection in clinical applications. Despite the advantage of using long reads, however, accurate SV calling and phasing are still challenging.

RESULTS

We introduce Duet, an SV detection tool optimized for SV calling and phasing using ONT data. The tool uses novel features integrated from both SV signatures and single-nucleotide polymorphism signatures, which can accurately distinguish SV haplotype from a false signal. Duet was benchmarked against state-of-the-art tools on multiple ONT sequencing datasets of sequencing coverage ranging from 8× to 40×. At low sequencing coverage of 8×, Duet performs better than all other tools in SV calling, SV genotyping and SV phasing. When the sequencing coverage is higher (20× to 40×), the F1-score for SV phasing is further improved in comparison to the performance of other tools, while its performance of SV genotyping and SV calling remains higher than other tools.

CONCLUSION

Duet can perform accurate SV calling, SV genotyping and SV phasing using low-coverage ONT data, making it very useful for low-coverage genomes. It has great performance when scaled to high-coverage genomes, which is adaptable to various clinical applications. Duet is open source and is available at https://github.com/yekaizhou/duet .

Genetic drift drives rapid speciation of an Arctic insular endemic shrew (Sorex pribilofensis)

Episodes of Quaternary environmental change shaped the genomes of extant species, influencing their response to contemporary environments, which are changing rapidly. Island endemics are among the most vulnerable to such change, accounting for a disproportionate number of recent extinctions. To prevent extinctions and conserve island biodiversity it is vital to combine knowledge of species' ecologies with their complex evolutionary histories. The Bering Sea has a history of cyclical island isolation and reconnection, coupled with modern rates of climate change that exceed global averages. The endangered Pribilof Island shrew (Sorex pribilofensis) is endemic to St. Paul Island, Alaska, which was isolated from mainland Beringia ~14,000 years ago by rising sea levels. Using ~11,000 single nucleotide polymorphisms, 17 microsatellites and mitochondrial sequence data, we test predictions about the evolutionary processes driving shrew speciation across Beringia. Our data show considerable differentiation of S. pribilofensis from mainland sibling species, relative to levels of divergence between mainland shrews. We also find a genome-wide loss of diversity and extremely low N_e for S. pribilofensis. We then show that intraspecific genetic diversity is significantly related to interspecific divergence, and that differentiation between S. pribilofensis and other Beringian shrews is highest across loci that are fixed in S. pribilofensis, indicating that strong drift has driven differentiation of this island species. Our findings show that drift as a consequence of Arctic climate cycling can rapidly reshape insular biodiversity. Arctic island species that lack genomic diversity and have evolved in response to past climate may have limited ability to respond to modern environmental changes.

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.

Rare variants implicate NMDA receptor signaling and cerebellar gene networks in risk for bipolar disorder

Bipolar disorder is an often-severe mental health condition characterized by alternation between extreme mood states of mania and depression. Despite strong heritability and the recent identification of 64 common variant risk loci of small effect, pathophysiological mechanisms remain unknown. Here, we analyzed genome sequences from 41 multiply-affected pedigrees and identified variants in 741 genes with nominally significant linkage or association with bipolar disorder. These 741 genes overlapped known risk genes for neurodevelopmental disorders and clustered within gene networks enriched for synaptic and nuclear functions. The top variant in this analysis - prioritized by statistical association, predicted deleteriousness, and network centrality - was a missense variant in the gene encoding D-amino acid oxidase (DAOG131V). Heterologous expression of DAOG131V in human cells resulted in decreased DAO protein abundance and enzymatic activity. In a knock-in mouse model of DAOG131, DaoG130V/+, we similarly found decreased DAO protein abundance in hindbrain regions, as well as enhanced stress susceptibility and blunted behavioral responses to pharmacological inhibition of N-methyl-D-aspartate receptors (NMDARs). RNA sequencing of cerebellar tissue revealed that DaoG130V resulted in decreased expression of two gene networks that are enriched for synaptic functions and for genes expressed, respectively, in Purkinje neurons or granule neurons. These gene networks were also down-regulated in the cerebellum of patients with bipolar disorder compared to healthy controls and were enriched for additional rare variants associated with bipolar disorder risk. These findings implicate dysregulation of NMDAR signaling and of gene expression in cerebellar neurons in bipolar disorder pathophysiology and provide insight into its genetic architecture.

The Coaching for Cognition in Alzheimer's (COCOA) trial: Study design

Comprehensive treatment of Alzheimer's disease (AD) requires not only pharmacologic treatment but also management of existing medical conditions and lifestyle modifications including diet, cognitive training, and exercise. We present the design and methodology for the Coaching for Cognition in Alzheimer's (COCOA) trial. AD and other dementias result from the interplay of multiple interacting dysfunctional biological systems. Monotherapies have had limited success. More interventional studies are needed to test the effectiveness of multimodal multi-domain therapies for dementia prevention and treatment. Multimodal therapies use multiple interventions to address multiple systemic causes and potentiators of cognitive decline and functional loss; they can be personalized, as different sets of etiologies and systems responsive to therapy may be present in different individuals. COCOA is designed to test the hypothesis that coached multimodal interventions beneficially alter the trajectory of cognitive decline for individuals on the spectrum of AD and related dementias (ADRD). COCOA is a two-arm prospective randomized controlled trial (RCT). COCOA collects psychometric, clinical, lifestyle, genomic, proteomic, metabolomic, and microbiome data at multiple timepoints across 2 years for each participant. These data enable systems biology analyses. One arm receives standard of care and generic healthy aging recommendations. The other arm receives standard of care and personalized data-driven remote coaching. The primary outcome measure is the Memory Performance Index (MPI), a measure of cognition. The MPI is a summary statistic of the MCI Screen (MCIS). Secondary outcome measures include the Functional Assessment Staging Test (FAST), a measure of function. COCOA began enrollment in January 2018. We hypothesize that multimodal interventions will ameliorate cognitive decline and that data-driven health coaching will increase compliance, assist in personalizing multimodal interventions, and improve outcomes for patients, particularly for those in the early stages of the AD spectrum.

Highlights

The Coaching for Cognition in Alzheimer's (COCOA) trial tests personalized multimodal lifestyle interventions for Alzheimer's disease and related dementias.Dense longitudinal molecular data will be useful for future studies.Increased use of Hill's criteria in analyses may advance knowledge generation.Remote coaching may be an effective intervention.Because lifestyle interventions are inexpensive, they may be particularly valuable in reducing global socioeconomic disparities in dementia care.

Information enhanced model selection for Gaussian graphical model with application to metabolomic data

In light of the low signal-to-noise nature of many large biological data sets, we propose a novel method to learn the structure of association networks using Gaussian graphical models combined with prior knowledge. Our strategy includes two parts. In the first part, we propose a model selection criterion called structural Bayesian information criterion, in which the prior structure is modeled and incorporated into Bayesian information criterion. It is shown that the popular extended Bayesian information criterion is a special case of structural Bayesian information criterion. In the second part, we propose a two-step algorithm to construct the candidate model pool. The algorithm is data-driven and the prior structure is embedded into the candidate model automatically. Theoretical investigation shows that under some mild conditions structural Bayesian information criterion is a consistent model selection criterion for high-dimensional Gaussian graphical model. Simulation studies validate the superiority of the proposed algorithm over the existing ones and show the robustness to the model misspecification. Application to relative concentration data from infant feces collected from subjects enrolled in a large molecular epidemiological cohort study validates that metabolic pathway involvement is a statistically significant factor for the conditional dependence between metabolites. Furthermore, new relationships among metabolites are discovered which can not be identified by the conventional methods of pathway analysis. Some of them have been widely recognized in biological literature.

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.

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.

Parental mosaicism for apparent de novo genetic variants: Scope, detection, and counseling challenges

The disease burden of de novo mutations (DNMs) has been evidenced only recently when the common application of next-generation sequencing technologies enabled their reliable and affordable detection through family-based clinical exome or genome sequencing. Implementation of exome sequencing into prenatal diagnostics revealed that up to 63% of pathogenic or likely pathogenic variants associated with fetal structural anomalies are apparently de novo, primarily for autosomal dominant disorders. Apparent DNMs have been considered to primarily occur as germline or zygotic events, with consequently negligible recurrence risks. However, there is now evidence that a considerable proportion of them are in fact inherited from a parent mosaic for the variant. Here, we review the burden of DNMs in prenatal diagnostics and the influence of parental mosaicism on the interpretation of apparent DNMs and discuss the challenges with detecting and quantifying parental mosaicism and its effect on recurrence risk. We also describe new bioinformatic and technological tools developed to assess mosaicism and discuss how they improve the accuracy of reproductive risk counseling when parental mosaicism is detected.

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.

A complete pedigree-based graph workflow for rare candidate variant analysis

Methods that use a linear genome reference for genome sequencing data analysis are reference-biased. In the field of clinical genetics for rare diseases, a resulting reduction in genotyping accuracy in some regions has likely prevented the resolution of some cases. Pangenome graphs embed population variation into a reference structure. Although pangenome graphs have helped to reduce reference mapping bias, further performance improvements are possible. We introduce VG-Pedigree, a pedigree-aware workflow based on the pangenome-mapping tool of Giraffe and the variant calling tool DeepTrio using a specially trained model for Giraffe-based alignments. We demonstrate mapping and variant calling improvements in both single-nucleotide variants (SNVs) and insertion and deletion (indel) variants over those produced by alignments created using BWA-MEM to a linear-reference and Giraffe mapping to a pangenome graph containing data from the 1000 Genomes Project. We have also adapted and upgraded deleterious-variant (DV) detecting methods and programs into a streamlined workflow. We used these workflows in combination to detect small lists of candidate DVs among 15 family quartets and quintets of the Undiagnosed Diseases Program (UDP). All candidate DVs that were previously diagnosed using the Mendelian models covered by the previously published methods were recapitulated by these workflows. The results of these experiments indicate that a slightly greater absolute count of DVs are detected in the proband population than in their matched unaffected siblings.

Rocks and clocks revised: New promises and challenges in dating the primate tree of life

In recent years, multiple technological and methodological advances have increased our ability to estimate phylogenies, leading to more accurate dating of the primate tree of life. Here we provide an overview of the limitations and potentials of some of these advancements and discuss how dated phylogenies provide the crucial temporal scale required to understand primate evolution. First, we review new methods, such as the total-evidence dating approach, that promise a better integration between the fossil record and molecular data. We then explore how the ever-increasing availability of genomic-level data for more primate species can impact our ability to accurately estimate timetrees. Finally, we discuss more recent applications of mutation rates to date divergence times. We highlight example studies that have applied these approaches to estimate divergence dates within primates. Our goal is to provide a critical overview of these new developments and explore the promises and challenges of their application in evolutionary anthropology.

Evolutionary Landscapes of Host-Virus Arms Races

Vertebrate immune systems suppress viral infection using both innate restriction factors and adaptive immunity. Viruses mutate to escape these defenses, driving hosts to counterevolve to regain fitness. This cycle recurs repeatedly, resulting in an evolutionary arms race whose outcome depends on the pace and likelihood of adaptation by host and viral genes. Although viruses evolve faster than their vertebrate hosts, their proteins are subject to numerous functional constraints that impact the probability of adaptation. These constraints are globally defined by evolutionary landscapes, which describe the fitness and adaptive potential of all possible mutations. We review deep mutational scanning experiments mapping the evolutionary landscapes of both host and viral proteins engaged in arms races. For restriction factors and some broadly neutralizing antibodies, landscapes favor the host, which may help to level the evolutionary playing field against rapidly evolving viruses. We discuss the biophysical underpinnings of these landscapes and their therapeutic implications.

Molecular Evolutionary Rate Predicts Intraspecific Genetic Polymorphism and Species-Specific Selection

It is unknown what determines genetic diversity and how genetic diversity is associated with various biological traits. In this work, we provide insight into these issues. By comparing genetic variation of 14,671 mammalian gene trees with thousands of individual human, chimpanzee, gorilla, mouse, and dog/wolf genomes, we found that intraspecific genetic diversity can be predicted by long-term molecular evolutionary rates rather than de novo mutation rates. This relationship was established during the early stage of mammalian evolution. Moreover, we developed a method to detect fluctuations of species-specific selection on genes based on the deviations of intraspecific genetic diversity predicted from long-term rates. We showed that the evolution of epithelial cells, rather than connective tissue, mainly contributed to morphological evolution of different species. For humans, evolution of the immune system and selective sweeps caused by infectious diseases are the most representative examples of adaptive evolution.

Detection of germline mosaicism in fathers of children with intellectual disability syndromes caused by de novo variants

Background

De novo variants are a common cause to rare intellectual disability syndromes, associated with low recurrence risk. However, when such variants occur pre‐zygotically in parental germ cells, the recurrence risk might be higher. Still, the recurrence risk estimates are mainly based on empirical data and the prevalence of germline mosaicism is often unknown.

Methods

To establish the prevalence of mosaicism in parents of children with intellectual disability syndromes caused by de novo variants, we performed droplet digital PCR on DNA extracted from blood (43 trios), and sperm (31 fathers).

Results

We detected low‐level mosaicism in sperm‐derived DNA but not in blood in the father of a child with Kleefstra syndrome caused by an EHMT1 variant. Additionally, we found a higher level of paternal mosaicism in sperm compared to blood in the father of a child with Gillespie syndrome caused by an ITPR1 variant.

Conclusion

By employing droplet digital PCR, we detected paternal germline mosaicism in two intellectual disability syndromes. In both cases, the mosaicism level was higher in sperm than blood, indicating that analysis of blood alone may underestimate germline mosaicism. Therefore, sperm analysis can be clinically useful to establish the recurrence risk for parents and improve genetic counselling.

Modeling dopamine dysfunction in autism spectrum disorder: From invertebrates to vertebrates

Autism Spectrum Disorder (ASD) is a highly heterogeneous neurodevelopmental disorder characterized by deficits in social communication and by patterns of restricted interests and/or repetitive behaviors. The Simons Foundation Autism Research Initiative's Human Gene and CNV Modules now list over 1000 genes implicated in ASD and over 2000 copy number variant loci reported in individuals with ASD. Given this ever-growing list of genetic changes associated with ASD, it has become evident that there is likely not a single genetic cause of this disorder nor a single neurobiological basis of this disorder. Instead, it is likely that many different neurobiological perturbations (which may represent subtypes of ASD) can result in the set of behavioral symptoms that we called ASD. One such of possible subtype of ASD may be associated with dopamine dysfunction. Precise regulation of synaptic dopamine (DA) is required for reward processing and behavioral learning, behaviors which are disrupted in ASD. Here we review evidence for DA dysfunction in ASD and in animal models of ASD. Further, we propose that these studies provide a scaffold for scientists and clinicians to consider subcategorizing the ASD diagnosis based on the genetic changes, neurobiological difference, and behavioral features identified in individuals with ASD.

Long-read genome sequencing provides molecular insights into scavenging and societal complexity in spotted hyena Crocuta crocuta

The spotted hyena (Crocuta crocuta) is a large and unique terrestrial carnivore. It is a particularly fascinating species due to its distinct phenotypic traits, especially its complex social structure and scavenging lifestyle, with associated high dietary exposure to microbial pathogens. However, the underlying molecular mechanisms related to these phenotypes remain elusive. Here, we sequenced and assembled a high-quality long-read genome of the spotted hyena, with a contig N50 length of ∼13.75 Mb. Based on comparative genomics, immunoglobulin family members (e.g., IGKV4-1) showed significant adaptive duplications in the spotted hyena and striped hyena. Furthermore, immune-related genes (e.g., CD8A, LAG3, and TLR3) experienced species-specific positive selection in the spotted hyena lineage. These results suggest that immune tolerance between the spotted hyena and closely related striped hyena has undergone adaptive divergence to cope with prolonged dietary exposure to microbial pathogens from scavenging. Furthermore, we provided the potential genetic insights underlying social complexity, hinting at social behavior and cognition. Specifically, the RECNE-associated genes (e.g., UGP2 and ACTR2) in the spotted hyena genome are involved in regulation of social communication. Taken together, our genomic analyses provide molecular insights into the scavenging lifestyle and societal complexity of spotted hyenas.

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.

Comparing Clinical and Genetic Characteristics of De Novo and Inherited COL1A1/COL1A2 Variants in a Large Chinese Cohort of Osteogenesis Imperfecta

PURPOSE

Nearly 85%-90% of osteogenesis imperfecta (OI) cases are caused by autosome dominant mutations of COL1A1 and COL1A2 genes, of which de novo mutations cover a large proportion, whereas their characteristics remain to be elucidated. This study aims to compare the differences in clinical and genetic characteristics of de novo and inherited COL1A1/COL1A2 mutations of OI, assess the average paternal and maternal age at conception in de novo mutations, and research the rate of nonpenetrance in inherited mutations.

MATERIALS AND METHODS

A retrospective comparison between de novo and inherited mutations was performed among 135 OI probands with COL1A1/COL1A2 mutations. Mutational analyses of all probands and their family members were completed by Sanger sequencing. A new clinical scoring system was developed to assess the clinical severity of OI quantitatively.

RESULTS

A total of 51 probands (37.78%) with de novo mutations and 84 probands (62.22%) with inherited mutations were grouped by the results of the parental gene verification. The proportion of clinical type III (P<0.001) and clinical scores (P<0.001) were significantly higher in de novo mutations. Missense mutations covered a slightly higher proportion of de novo COL1A1 mutations (46.34%) compared with inherited COL1A1 mutations (33.33%), however, lacking a significant difference (P=0.1923). The mean BMD Z/T-score at the lumbar spine in de novo mutations was -2.3 ± 1.5, lower than inherited mutations (-1.7 ± 1.8), but lacking statistical significance (P=0.0742). There was no significant difference between the two groups in OI-related phenotypes (like fracture frequency, blue sclera, and hearing loss) and biochemical indexes. In de novo mutations, the average paternal and maternal age at conception was 29.2 (P<0.05) and 26.8 (P<0.0001), respectively, which were significantly younger than the average gestational age of the population. Additionally, 98.04% of pedigrees (50/51) with de novo mutations were spontaneous conception. The rate of nonpenetrance of parents with pathogenic variants in the inherited mutation group was 25.64% (20/78).

CONCLUSIONS

Our data revealed that the proportion of clinical type III and clinical scores were significantly higher in de novo mutations than in inherited mutations, demonstrating that de novo mutations are more damaging because they have not undergone purifying selection.

Occam's razor dulled: the occurrence of multiple genetic diagnoses

PURPOSE OF REVIEW

A single genetic diagnosis, especially from the analysis of a limited number of genes, may not signal the end of a diagnostic odyssey. When a patient with a genetic syndrome presents with symptoms that are not usually associated with their disease phenotype, additional genetic testing is warranted.

RECENT FINDINGS

Although multiple co-existing genetic diagnoses may sound unlikely, many recent studies and case reports have demonstrated that this scenario is more common than expected. Studies involving whole exome and genome sequencing have identified a frequency of multiple genetic diagnoses and have identified clinical findings that make a second diagnosis more likely, which we have seen reflected in recent cases from our own clinic and consult service. These include multisystem disease, consanguinity, well described aneuploidies with rare or new symptoms, and complex structural chromosomal anomalies which may include multiple chromosomes and breakpoints that disrupt gene function.

SUMMARY

Identifying a second diagnosis can have vast implications for patient management and counseling. Patients can be followed with appropriate medical screening and early interventions to support optimal child development. Furthermore, the patient's family can be impacted by ending the diagnostic odyssey, providing testing for other at-risk family members, and offering prenatal options.

HaplotypeTools: a toolkit for accurately identifying recombination and recombinant genotypes

Background

Identifying haplotypes is central to sequence analysis in diploid or polyploid genomes. Despite this, there remains a lack of research and tools designed for physical phasing and its downstream analysis.

Results

HaplotypeTools is a new toolset to phase variant sites using VCF and BAM files and to analyse phased VCFs. Phasing is achieved via the identification of reads overlapping ≥ 2 heterozygous positions and then extended by additional reads, a process that can be parallelized across a computer cluster. HaplotypeTools includes various utility scripts for downstream analysis including crossover detection and phylogenetic placement of haplotypes to other lineages or species. HaplotypeTools was assessed for accuracy against WhatsHap using simulated short and long reads, demonstrating higher accuracy, albeit with reduced haplotype length. HaplotypeTools was also tested on real Illumina data to determine the ancestry of hybrid fungal isolate Batrachochytrium dendrobatidis (Bd) SA-EC3, finding 80% of haplotypes across the genome phylogenetically cluster with parental lineages BdGPL (39%) and BdCAPE (41%), indicating those are the parental lineages. Finally, ~ 99% of phasing was conserved between overlapping phase groups between SA-EC3 and either parental lineage, indicating mitotic gene conversion/parasexuality as the mechanism of recombination for this hybrid isolate. HaplotypeTools is open source and freely available from https://github.com/rhysf/HaplotypeTools under the MIT License.

Conclusions

HaplotypeTools is a powerful resource for analyzing hybrid or recombinant diploid or polyploid genomes and identifying parental ancestry for sub-genomic regions.

GC-biased gene conversion in X-chromosome palindromes conserved in human, chimpanzee, and rhesus macaque

Gene conversion is GC-biased across a wide range of taxa. Large palindromes on mammalian sex chromosomes undergo frequent gene conversion that maintains arm-to-arm sequence identity greater than 99%, which may increase their susceptibility to the effects of GC-biased gene conversion. Here, we demonstrate a striking history of GC-biased gene conversion in 12 palindromes conserved on the X chromosomes of human, chimpanzee, and rhesus macaque. Primate X-chromosome palindrome arms have significantly higher GC content than flanking single-copy sequences. Nucleotide replacements that occurred in human and chimpanzee palindrome arms over the past 7 million years are one-and-a-half times as GC-rich as the ancestral bases they replaced. Using simulations, we show that our observed pattern of nucleotide replacements is consistent with GC-biased gene conversion with a magnitude of 70%, similar to previously reported values based on analyses of human meioses. However, GC-biased gene conversion since the divergence of human and rhesus macaque explains only a fraction of the observed difference in GC content between palindrome arms and flanking sequence, suggesting that palindromes are older than 29 million years and/or had elevated GC content at the time of their formation. This work supports a greater than 2:1 preference for GC bases over AT bases during gene conversion and demonstrates that the evolution and composition of mammalian sex chromosome palindromes is strongly influenced by GC-biased gene conversion.

De Novo Germline and Somatic Variants Convergently Promote Endothelial-to-Mesenchymal Transition in Simplex Brain Arteriovenous Malformation

[Figure: see text].

Toward a methodology for evaluating DNA variants in nuclear families

The genetic underpinnings of most pediatric-cancer cases are unknown. Population-based studies use large sample sizes but have accounted for only a small proportion of the estimated heritability of pediatric cancers. Pedigree-based studies are infeasible for most human populations. One alternative is to collect genetic data from a single nuclear family and use inheritance patterns within the family to filter candidate variants. This approach can be applied to common and rare variants, including those that are private to a given family or to an affected individual. We evaluated this approach using genetic data from three nuclear families with 5, 4, and 7 children, respectively. Only one child in each nuclear family had been diagnosed with cancer, and neither parent had been affected. Diagnoses for the affected children were benign low-grade astrocytoma, Wilms tumor (stage 2), and Burkitt's lymphoma, respectively. We used whole-genome sequencing to profile normal cells from each family member and a linked-read technology for genomic phasing. For initial variant filtering, we used global minor allele frequencies, deleteriousness scores, and functional-impact annotations. Next, we used genetic variation in the unaffected siblings as a guide to filter the remaining variants. As a way to evaluate our ability to detect variant(s) that may be relevant to disease status, the corresponding author blinded the primary author to affected status; the primary author then assigned a risk score to each child. Based on this evidence, the primary author predicted which child had been affected in each family. The primary author's prediction was correct for the child who had been diagnosed with a Wilms tumor; the child with Burkitt's lymphoma had the second-highest risk score among the seven children in that family. This study demonstrates a methodology for filtering and evaluating candidate genomic variants and genes within nuclear families that may merit further exploration.

Genome-wide mutational signatures revealed distinct developmental paths for human B cell lymphomas

Both somatic hypermutation (SHM) and class switch recombination (CSR) are initiated by activation-induced cytidine deaminase (AID). Dysregulation of these processes has been linked to B cell lymphomagenesis. Here we performed an in-depth analysis of diffuse large B cell lymphoma (DLBCL) and follicular lymphoma (FL) genomes. We characterized seven genomic mutational signatures, including two B cell tumor-specific signatures, one of which is novel and associated with aberrant SHM. We further identified two major mutational signatures (K1 and K2) of clustered mutations (kataegis) resulting from the activities of AID or error-prone DNA polymerase η, respectively. K1 was associated with the immunoglobulin (Ig) switch region mutations/translocations and the ABC subtype of DLBCL, whereas K2 was related to the Ig variable region mutations and the GCB subtype of DLBCL and FL. Similar patterns were also observed in chronic lymphocytic leukemia subtypes. Thus, alterations associated with aberrant CSR and SHM activities can be linked to distinct developmental paths for different subtypes of B cell lymphomas.

Packpred: Predicting the Functional Effect of Missense Mutations

Predicting the functional consequences of single point mutations has relevance to protein function annotation and to clinical analysis/diagnosis. We developed and tested Packpred that makes use of a multi-body clique statistical potential in combination with a depth-dependent amino acid substitution matrix (FADHM) and positional Shannon entropy to predict the functional consequences of point mutations in proteins. Parameters were trained over a saturation mutagenesis data set of T4-lysozyme (1,966 mutations). The method was tested over another saturation mutagenesis data set (CcdB; 1,534 mutations) and the Missense3D data set (4,099 mutations). The performance of Packpred was compared against those of six other contemporary methods. With MCC values of 0.42, 0.47, and 0.36 on the training and testing data sets, respectively, Packpred outperforms all methods in all data sets, with the exception of marginally underperforming in comparison to FADHM in the CcdB data set. A meta server analysis was performed that chose best performing methods of wild-type amino acids and for wild-type mutant amino acid pairs. This led to an increase in the MCC value of 0.40 and 0.51 for the two meta predictors, respectively, on the Missense3D data set. We conjecture that it is possible to improve accuracy with better meta predictors as among the seven methods compared, at least one method or another is able to correctly predict ∼99% of the data.

Evidence of local adaptation despite strong drift in a Neotropical patchily distributed bromeliad

Both genetic drift and divergent selection are predicted to be drivers of population differentiation across patchy habitats, but the extent to which these forces act on natural populations to shape traits is strongly affected by species' ecological features. In this study, we infer the genomic structure of Pitcairnia lanuginosa, a widespread herbaceous perennial plant with a patchy distribution. We sampled populations in the Brazilian Cerrado and the Central Andean Yungas and discovered and genotyped SNP markers using double-digest restriction-site associated DNA sequencing. In addition, we analyzed ecophysiological traits obtained from a common garden experiment and compared patterns of phenotypic and genetic divergence (P_ST-F_ST comparisons) in a subset of populations from the Cerrado. Our results from molecular analyses pointed to extremely low genetic diversity and a remarkable population differentiation, supporting a major role of genetic drift. Approximately 0.3% of genotyped SNPs were flagged as differentiation outliers by at least two distinct methods, and Bayesian generalized linear mixed models revealed a signature of isolation by environment in addition to isolation by distance for high-differentiation outlier SNPs among the Cerrado populations. P_ST-F_ST comparisons suggested divergent selection on two ecophysiological traits linked to drought tolerance. We showed that these traits vary among populations, although without any particular macro-spatial pattern, suggesting local adaptation to differences in micro-habitats. Our study shows that selection might be a relevant force, particularly for traits involved in drought stress, even for populations experiencing strong drift, which improves our knowledge on eco-evolutionary processes acting on non-continuously distributed species.

Computational and Mass Spectrometry-Based Approach Identify Deleterious Non-Synonymous Single Nucleotide Polymorphisms (nsSNPs) in JMJD6

The jumonji domain-containing protein 6 (JMJD6) gene catalyzes the arginine demethylation and lysine hydroxylation of histone and a growing list of its known substrate molecules, including p53 and U2AF65, suggesting a possible role in mRNA splicing and transcription in cancer progression. Mass spectrometry-based technology offers the opportunity to detect SNP variants accurately and effectively. In our study, we conducted a combined computational and filtration workflow to predict the nonsynonymous single nucleotide polymorphisms (nsSNPs) present in JMJD6, followed by a liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis and validation. The computational approaches SIFT, PolyPhen-2, SNAP, I-Mutant 2.0, PhD-SNP, PANTHER, and SNPS&GO were integrated to screen out the predicted damaging/deleterious nsSNPs. Through the three-dimensional structure of JMJD6, H187R (rs1159480887) was selected as a candidate for validation. The validation experiments showed that the mutation of this nsSNP in JMJD6 obviously affected mRNA splicing or the transcription of downstream genes through the reduced lysyl-hydroxylase activity of its substrates, U2AF65 and p53, further indicating the accuracy of this prediction method. This research provides an effective computational workflow for researchers with an opportunity to select prominent deleterious nsSNPs and, thus, remains promising for examining the dysfunction of proteins.

Discovery of genomic variation across a generation

Over the past 30 years (the timespan of a generation), advances in genomics technologies have revealed tremendous and unexpected variation in the human genome and have provided increasingly accurate answers to long-standing questions of how much genetic variation exists in human populations and to what degree the DNA complement changes between parents and offspring. Tracking the characteristics of these inherited and spontaneous (or de novo) variations has been the basis of the study of human genetic disease. From genome-wide microarray and next-generation sequencing scans, we now know that each human genome contains over 3 million single nucleotide variants when compared with the ~ 3 billion base pairs in the human reference genome, along with roughly an order of magnitude more DNA—approximately 30 megabase pairs (Mb)—being ‘structurally variable’, mostly in the form of indels and copy number changes. Additional large-scale variations include balanced inversions (average of 18 Mb) and complex, difficult-to-resolve alterations. Collectively, ~1% of an individual’s genome will differ from the human reference sequence. When comparing across a generation, fewer than 100 new genetic variants are typically detected in the euchromatic portion of a child’s genome. Driven by increasingly higher-resolution and higher-throughput sequencing technologies, newer and more accurate databases of genetic variation (for instance, more comprehensive structural variation data and phasing of combinations of variants along chromosomes) of worldwide populations will emerge to underpin the next era of discovery in human molecular genetics.

Contrasting Gene Decay in Subterranean Vertebrates: Insights from Cavefishes and Fossorial Mammals

Evolution sometimes proceeds by loss, especially when structures and genes become dispensable after an environmental shift relaxes functional constraints. Subterranean vertebrates are outstanding models to analyze this process, and gene decay can serve as a readout. We sought to understand some general principles on the extent and tempo of the decay of genes involved in vision, circadian clock, and pigmentation in cavefishes. The analysis of the genomes of two Cuban species belonging to the genus Lucifuga provided evidence for the largest loss of eye-specific genes and nonvisual opsin genes reported so far in cavefishes. Comparisons with a recently evolved cave population of Astyanax mexicanus and three species belonging to the Chinese tetraploid genus Sinocyclocheilus revealed the combined effects of the level of eye regression, time, and genome ploidy on eye-specific gene pseudogenization. The limited extent of gene decay in all these cavefishes and the very small number of loss-of-function mutations per pseudogene suggest that their eye degeneration may not be very ancient, ranging from early to late Pleistocene. This is in sharp contrast with the identification of several vision genes carrying many loss-of-function mutations in ancient fossorial mammals, further suggesting that blind fishes cannot thrive more than a few million years in cave ecosystems.

The challenge and promise of estimating the de novo mutation rate from whole-genome comparisons among closely related individuals

Germline mutations are the raw material for natural selection, driving species evolution and the generation of earth's biodiversity. Without this driver of genetic diversity, life on earth would stagnate. Yet, it is a double-edged sword. An excess of mutations can have devastating effects on fitness and population viability. It is therefore one of the great challenges of molecular ecology to determine the rate and mechanisms by which these mutations accrue across the tree of life. Advances in high-throughput sequencing technologies are providing new opportunities for characterizing the rates and mutational spectra within species and populations thus informing essential evolutionary parameters such as the timing of speciation events, the intricacies of historical demography, and the degree to which lineages are subject to the burdens of mutational load. Here, we will focus on both the challenge and promise of whole-genome comparisons among parents and their offspring from known pedigrees for the detection of germline mutations as they arise in a single generation. The potential of these studies is high, but the field is still in its infancy and much uncertainty remains. Namely, the technical challenges are daunting given that pedigree-based genome comparisons are essentially searching for needles in a haystack given the very low signal to noise ratio. Despite the challenges, we predict that rapidly developing methods for whole-genome comparisons hold great promise for integrating empirically derived estimates of de novo mutation rates and mutation spectra across many molecular ecological applications.

Cataloguing experimentally confirmed 80.7 kb-long ACKR1 haplotypes from the 1000 Genomes Project database

Background

Clinically effective and safe genotyping relies on correct reference sequences, often represented by haplotypes. The 1000 Genomes Project recorded individual genotypes across 26 different populations and, using computerized genotype phasing, reported haplotype data. In contrast, we identified long reference sequences by analyzing the homozygous genomic regions in this online database, a concept that has rarely been reported since next generation sequencing data became available.

Study design and methods

Phased genotype data for a 80.6 kb region of chromosome 1 was downloaded for all 2,504 unrelated individuals of the 1000 Genome Project Phase 3 cohort. The data was centered on the ACKR1 gene and bordered by the CADM3 and FCER1A genes. Individuals with heterozygosity at a single site or with complete homozygosity allowed unambiguous assignment of an ACKR1 haplotype. A computer algorithm was developed for extracting these haplotypes from the 1000 Genome Project in an automated fashion. A manual analysis validated the data extracted by the algorithm.

Results

We confirmed 902 ACKR1 haplotypes of varying lengths, the longest at 80,584 nucleotides and shortest at 1,901 nucleotides. The combined length of haplotype sequences comprised 19,895,388 nucleotides with a median of 16,014 nucleotides. Based on our approach, all haplotypes can be considered experimentally confirmed and not affected by the known errors of computerized genotype phasing.

Conclusions

Tracts of homozygosity can provide definitive reference sequences for any gene. They are particularly useful when observed in unrelated individuals of large scale sequence databases. As a proof of principle, we explored the 1000 Genomes Project database for ACKR1 gene data and mined long haplotypes. These haplotypes are useful for high throughput analysis with next generation sequencing. Our approach is scalable, using automated bioinformatics tools, and can be applied to any gene.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04169-6.

Stability across the Whole Nuclear Genome in the Presence and Absence of DNA Mismatch Repair

We describe the contribution of DNA mismatch repair (MMR) to the stability of the eukaryotic nuclear genome as determined by whole-genome sequencing. To date, wild-type nuclear genome mutation rates are known for over 40 eukaryotic species, while measurements in mismatch repair-defective organisms are fewer in number and are concentrated on Saccharomyces cerevisiae and human tumors. Well-studied organisms include Drosophila melanogaster and Mus musculus, while less genetically tractable species include great apes and long-lived trees. A variety of techniques have been developed to gather mutation rates, either per generation or per cell division. Generational rates are described through whole-organism mutation accumulation experiments and through offspring–parent sequencing, or they have been identified by descent. Rates per somatic cell division have been estimated from cell line mutation accumulation experiments, from systemic variant allele frequencies, and from widely spaced samples with known cell divisions per unit of tissue growth. The latter methods are also used to estimate generational mutation rates for large organisms that lack dedicated germlines, such as trees and hyphal fungi. Mechanistic studies involving genetic manipulation of MMR genes prior to mutation rate determination are thus far confined to yeast, Arabidopsis thaliana, Caenorhabditis elegans, and one chicken cell line. A great deal of work in wild-type organisms has begun to establish a sound baseline, but far more work is needed to uncover the variety of MMR across eukaryotes. Nonetheless, the few MMR studies reported to date indicate that MMR contributes 100-fold or more to genome stability, and they have uncovered insights that would have been impossible to obtain using reporter gene assays.

Identification of two pathogenic mutations in SORL1 in early-onset Alzheimer's disease

The sortilin-related receptor 1 (SORL1) gene has been the subject of many studies focusing on frequent polymorphisms, which is associated with increased risk for Alzheimer's Disease (AD). By whole-exome sequencing (WES), we identified two pathogenic missense mutations c.579C > G (p.F193L) and c.1397A > G (p.N466S) in SORL1. The two mutations were located in the same protein domain, and the two unrelated probands both had an onset of memory problems at less than 65 years of age, but their clinical manifestations and cranial imaging are different. The protein structure and function affected by these mutations were predicted using bioinformatics analysis, which suggested they were pathogenic. 3D protein structural analysis revealed that these amino acid substitutions might result in instability of protein structure and adverse intramolecular interactions. These findings suggest that both F193L and N466S should be thought as potential causative mutations in early-onset Alzheimer's disease (EOAD) patients. Further functional studies are warranted to evaluate their roles in the pathogenesis of AD.

Factors Associated with Mutations: Their Matching Rates to Cardiovascular and Neurological Diseases

Monogenic hypertension is rare and caused by genetic mutations, but whether factors associated with mutations are disease-specific remains uncertain. Given two factors associated with high mutation rates, we tested how many previously known genes match with (i) proximity to telomeres or (ii) high adenine and thymine content in cardiovascular diseases (CVDs) related to vascular stiffening. We extracted genomic information using a genome data viewer. In human chromosomes, 64 of 79 genetic loci involving >25 rare mutations and single nucleotide polymorphisms satisfied (i) or (ii), resulting in an 81% matching rate. However, this high matching rate was no longer observed as we checked the two factors in genes associated with essential hypertension (EH), thoracic aortic aneurysm (TAA), and congenital heart disease (CHD), resulting in matching rates of 53%, 70%, and 75%, respectively. A matching of telomere proximity or high adenine and thymine content projects the list of loci involving rare mutations of monogenic hypertension better than those of other CVDs, likely due to adoption of rigorous criteria for true-positive signals. Our data suggest that the factor–disease matching rate is an accurate tool that can explain deleterious mutations of monogenic hypertension at a >80% match—unlike the relatively lower matching rates found in human genes of EH, TAA, CHD, and familial Parkinson’s disease.

The germline mutational process in rhesus macaque and its implications for phylogenetic dating

BACKGROUND

Understanding the rate and pattern of germline mutations is of fundamental importance for understanding evolutionary processes.

RESULTS

Here we analyzed 19 parent-offspring trios of rhesus macaques (Macaca mulatta) at high sequencing coverage of ∼76× per individual and estimated a mean rate of 0.77 × 10-8de novo mutations per site per generation (95% CI: 0.69 × 10-8 to 0.85 × 10-8). By phasing 50% of the mutations to parental origins, we found that the mutation rate is positively correlated with the paternal age. The paternal lineage contributed a mean of 81% of the de novo mutations, with a trend of an increasing male contribution for older fathers. Approximately 3.5% of de novo mutations were shared between siblings, with no parental bias, suggesting that they arose from early development (postzygotic) stages. Finally, the divergence times between closely related primates calculated on the basis of the yearly mutation rate of rhesus macaque generally reconcile with divergence estimated with molecular clock methods, except for the Cercopithecoidea/Hominoidea molecular divergence dated at 58 Mya using our new estimate of the yearly mutation rate.

CONCLUSIONS

When compared to the traditional molecular clock methods, new estimated rates from pedigree samples can provide insights into the evolution of well-studied groups such as primates.

Genetic variation in the Mauritian cynomolgus macaque population reflects variation in the human population

The cynomolgus macaque is an important species for preclinical research, however the extent of genetic variation in this population and its similarity to the human population is not well understood. Exome sequencing was conducted for 101 cynomolgus macaques to characterize genetic variation. The variant distribution frequency was 7.81 variants per kilobase across the sequenced regions, with a total of 2,770,009 single nucleotide variants identified from 2,996,041 loci. A large portion (85.6%) had minor allele frequencies greater than 5%. Enriched pathways for genes with high genetic diversity (≥10 variants per kilobase) were those involving signaling peptides and immune response. Compared to human, the variant distribution frequency and nucleotide diversity in the macaque exome was approximately 4 times greater; however the ratio of non-synonymous to synonymous variants was similar (0.735 and 0.831, respectively). Understanding genetic variability in cynomolgus macaques will enable better interpretation and human translation of phenotypic variability in this species.

CRISPR/Cas9-mediated correction of MITF homozygous point mutation in a Waardenburg syndrome 2A pig model

Gene therapy for curing congenital human diseases is promising, but the feasibility and safety need to be further evaluated. In this study, based on a pig model that carries the c.740T>C (L247S) mutation in MITF with an inheritance pattern and clinical pathology that mimics Waardenburg syndrome 2A (WS2A), we corrected the point mutation by the CRISPR-Cas9 system in the mutant fibroblast cells using single-stranded oligodeoxynucleotide (ssODN) and long donor plasmid DNA as the repair template. By using long donor DNA, precise correction of this point mutation was achieved. The corrected cells were then used as the donor cell for somatic cell nuclear transfer (SCNT) to produce piglets, which exhibited a successfully rescued phenotype of WS2A, including anophthalmia and hearing loss. Furthermore, engineered base editors (BEs) were exploited to make the correction in mutant porcine fibroblast cells and early embryos. The correction efficiency was greatly improved, whereas substantial off-targeting mutations were detected, raising a safety concern for their potential applications in gene therapy. Thus, we explored the possibility of precise correction of WS2A-causing gene mutation by the CRISPR-Cas9 system in a large-animal model, suggesting great prospects for its future applications in treating human genetic diseases.

De novo structural mutation rates and gamete-of-origin biases revealed through genome sequencing of 2,396 families

Each human genome includes de novo mutations that arose during gametogenesis. While these germline mutations represent a fundamental source of new genetic diversity, they can also create deleterious alleles that impact fitness. Whereas the rate and patterns of point mutations in the human germline are now well understood, far less is known about the frequency and features that impact de novo structural variants (dnSVs). We report a family-based study of germline mutations among 9,599 human genomes from 33 multigenerational CEPH-Utah families and 2,384 families from the Simons Foundation Autism Research Initiative. We find that de novo structural mutations detected by alignment-based, short-read WGS occur at an overall rate of at least 0.160 events per genome in unaffected individuals, and we observe a significantly higher rate (0.206 per genome) in ASD-affected individuals. In both probands and unaffected samples, nearly 73% of de novo structural mutations arose in paternal gametes, and we predict most de novo structural mutations to be caused by mutational mechanisms that do not require sequence homology. After multiple testing correction, we did not observe a statistically significant correlation between parental age and the rate of de novo structural variation in offspring. These results highlight that a spectrum of mutational mechanisms contribute to germline structural mutations and that these mechanisms most likely have markedly different rates and selective pressures than those leading to point mutations.

Whole Genome Sequencing in the Evaluation of Fetal Structural Anomalies: A Parallel Test with Chromosomal Microarray Plus Whole Exome Sequencing

Whole genome sequencing (WGS) is a powerful tool for postnatal genetic diagnosis, but relevant clinical studies in the field of prenatal diagnosis are limited. The present study aimed to prospectively evaluate the utility of WGS compared with chromosomal microarray (CMA) and whole exome sequencing (WES) in the prenatal diagnosis of fetal structural anomalies. We performed trio WGS (≈40-fold) in parallel with CMA in 111 fetuses with structural or growth anomalies, and sequentially performed WES when CMA was negative (CMA plus WES). In comparison, WGS not only detected all pathogenic genetic variants in 22 diagnosed cases identified by CMA plus WES, yielding a diagnostic rate of 19.8% (22/110), but also provided additional and clinically significant information, including a case of balanced translocations and a case of intrauterine infection, which might not be detectable by CMA or WES. WGS also required less DNA (100 ng) as input and could provide a rapid turnaround time (TAT, 18 ± 6 days) compared with that (31 ± 8 days) of the CMA plus WES. Our results showed that WGS provided more comprehensive and precise genetic information with a rapid TAT and less DNA required than CMA plus WES, which enables it as an alternative prenatal diagnosis test for fetal structural anomalies.

Human molecular evolutionary rate, time dependency and transient polymorphism effects viewed through ancient and modern mitochondrial DNA genomes

Human evolutionary genetics gives a chronological framework to interpret the human history. It is based on the molecular clock hypothesis that suppose a straightforward relationship between the mutation rate and the substitution rate with independence of other factors as demography dynamics. Analyzing ancient and modern human complete mitochondrial genomes we show here that, along the time, the substitution rate can be significantly slower or faster than the average germline mutation rate confirming a time dependence effect mainly attributable to changes in the effective population size of the human populations, with an exponential growth in recent times. We also detect that transient polymorphisms play a slowdown role in the evolutionary rate deduced from haplogroup intraspecific trees. Finally, we propose the use of the most divergent lineages within haplogroups as a practical approach to correct these molecular clock mismatches.

SequencErr: measuring and suppressing sequencer errors in next-generation sequencing data

Background

There is currently no method to precisely measure the errors that occur in the sequencing instrument/sequencer, which is critical for next-generation sequencing applications aimed at discovering the genetic makeup of heterogeneous cellular populations.

Results

We propose a novel computational method, SequencErr, to address this challenge by measuring the base correspondence between overlapping regions in forward and reverse reads. An analysis of 3777 public datasets from 75 research institutions in 18 countries revealed the sequencer error rate to be ~ 10 per million (pm) and 1.4% of sequencers and 2.7% of flow cells have error rates > 100 pm. At the flow cell level, error rates are elevated in the bottom surfaces and > 90% of HiSeq and NovaSeq flow cells have at least one outlier error-prone tile. By sequencing a common DNA library on different sequencers, we demonstrate that sequencers with high error rates have reduced overall sequencing accuracy, and removal of outlier error-prone tiles improves sequencing accuracy. We demonstrate that SequencErr can reveal novel insights relative to the popular quality control method FastQC and achieve a 10-fold lower error rate than popular error correction methods including Lighter and Musket.

Conclusions

Our study reveals novel insights into the nature of DNA sequencing errors incurred on DNA sequencers. Our method can be used to assess, calibrate, and monitor sequencer accuracy, and to computationally suppress sequencer errors in existing datasets.

Family-Based Whole Genome Sequencing Identified Novel Variants in ABCA5 Gene in a Patient with Idiopathic Ventricular Tachycardia

Idiopathic ventricular tachycardia (IVT) is the major cause of sudden cardiac death. Patients with IVT were usually manifested without structural heart disease. In this present study, we performed family-based whole genome sequencing (WGS) and Sanger sequencing for a 5-year-old Chinese boy with IVT and all the unaffected family members in order to identify the candidate gene and disease-causing mutation underlying the disease phenotype. Results showed that a novel heterozygous single-nucleotide duplication (c.128dup) and a novel heterozygous missense (c.3328A > G) variant in ABCA5 gene were identified in the proband. The single-nucleotide duplication (c.128dupT), inherited from his father and patrilineal grandfather, leads to a frameshift which results into the formation of a truncated ABCA5 protein of 50 (p.Leu43Phefs*8) amino acids. Hence, it is a loss-of-function mutation. The missense (c.3328A > G) variant, inherited from his mother, leads to the replacement of isoleucine by valine at the position of 1110 (p.Ile1110Val) of the ABCA5 protein. Multiple sequence alignment showed that p.Ile1110 is evolutionarily conserved among several species indicating both the structural and functional significance of the p.Ile1110 residue in the wild-type ABCA5 protein. Quantitative RT-PCR showed that the ABCA5 mRNA expression levels were decreased in the proband. These two novel variants of ABCA5 gene were co-segregated well among all the members of this family. Our present study also strongly supports the importance of using family-based whole genome sequencing for identifying novel candidate genes associated with IVT.

Accuracy and efficiency of germline variant calling pipelines for human genome data

Advances in next-generation sequencing technology have enabled whole genome sequencing (WGS) to be widely used for identification of causal variants in a spectrum of genetic-related disorders, and provided new insight into how genetic polymorphisms affect disease phenotypes. The development of different bioinformatics pipelines has continuously improved the variant analysis of WGS data. However, there is a necessity for a systematic performance comparison of these pipelines to provide guidance on the application of WGS-based scientific and clinical genomics. In this study, we evaluated the performance of three variant calling pipelines (GATK, DRAGEN and DeepVariant) using the Genome in a Bottle Consortium, "synthetic-diploid" and simulated WGS datasets. DRAGEN and DeepVariant show better accuracy in SNP and indel calling, with no significant differences in their F1-score. DRAGEN platform offers accuracy, flexibility and a highly-efficient execution speed, and therefore superior performance in the analysis of WGS data on a large scale. The combination of DRAGEN and DeepVariant also suggests a good balance of accuracy and efficiency as an alternative solution for germline variant detection in further applications. Our results facilitate the standardization of benchmarking analysis of bioinformatics pipelines for reliable variant detection, which is critical in genetics-based medical research and clinical applications.

Accuracy and efficiency of germline variant calling pipelines for human genome data

Advances in next-generation sequencing technology have enabled whole genome sequencing (WGS) to be widely used for identification of causal variants in a spectrum of genetic-related disorders, and provided new insight into how genetic polymorphisms affect disease phenotypes. The development of different bioinformatics pipelines has continuously improved the variant analysis of WGS data. However, there is a necessity for a systematic performance comparison of these pipelines to provide guidance on the application of WGS-based scientific and clinical genomics. In this study, we evaluated the performance of three variant calling pipelines (GATK, DRAGEN and DeepVariant) using the Genome in a Bottle Consortium, "synthetic-diploid" and simulated WGS datasets. DRAGEN and DeepVariant show better accuracy in SNP and indel calling, with no significant differences in their F1-score. DRAGEN platform offers accuracy, flexibility and a highly-efficient execution speed, and therefore superior performance in the analysis of WGS data on a large scale. The combination of DRAGEN and DeepVariant also suggests a good balance of accuracy and efficiency as an alternative solution for germline variant detection in further applications. Our results facilitate the standardization of benchmarking analysis of bioinformatics pipelines for reliable variant detection, which is critical in genetics-based medical research and clinical applications.

Next Generation Sequencing and Bioinformatics Analysis of Family Genetic Inheritance

Mendelian and complex genetic trait diseases continue to burden and affect society both socially and economically. The lack of effective tests has hampered diagnosis thus, the affected lack proper prognosis. Mendelian diseases are caused by genetic mutations in a singular gene while complex trait diseases are caused by the accumulation of mutations in either linked or unlinked genomic regions. Significant advances have been made in identifying novel diseases associated mutations especially with the introduction of next generation and third generation sequencing. Regardless, some diseases are still without diagnosis as most tests rely on SNP genotyping panels developed from population based genetic analyses. Analysis of family genetic inheritance using whole genomes, whole exomes or a panel of genes has been shown to be effective in identifying disease-causing mutations. In this review, we discuss next generation and third generation sequencing platforms, bioinformatic tools and genetic resources commonly used to analyze family based genomic data with a focus on identifying inherited or novel disease-causing mutations. Additionally, we also highlight the analytical, ethical and regulatory challenges associated with analyzing personal genomes which constitute the data used for family genetic inheritance.

Elevated DHODH expression promotes cell proliferation via stabilizing β-catenin in esophageal squamous cell carcinoma

As a key enzyme in de novo pyrimidine biosynthesis, the expression level of dihydroorotate dehydrogenase (DHODH) has been reported to be elevated in various types of malignant tumors and its tumor-promoting effect was considered to relate to its pyrimidine synthesis function. Here, we revealed one intriguing potential mechanism that DHODH modulated β-catenin signaling in esophageal squamous cell carcinoma (ESCC). We demonstrated that DHODH directly bound to the NH2 terminal of β-catenin, thereby, interrupting the interaction of GSK3β with β-catenin and leading to the abrogation of β-catenin degradation and accumulation of β-catenin in the nucleus, which in turn, resulted in the activation of β-catenin downstream genes, including CCND1, E2F3, Nanog, and OCT4. We further demonstrated that the regulation of β-catenin by DHODH was independent of DHODH catalyzing activity. Univariate and multivariate analyses suggested that DHODH expression might be an independent prognostic factor for ESCC patients. Collectively, our study highlights the pivotal role of DHODH mediated β-catenin signaling and indicates that DHODH may act as a multi-functional switcher from catalyzing pyrimidine metabolism to regulating tumor-related signaling pathways in ESCC.