Quantifying negative selection on synonymous variants.

Widespread adoption of DNA sequencing has resulted in large numbers of genetic variants, whose contribution to disease is not easily determined. Although many types of variation are known to disrupt cellular processes in predictable ways, for some categories of variants, the effects may not be directly detectable. A particular example is synonymous variants, that is, those single-nucleotide variants that create a codon substitution, such that the produced amino acid sequence is unaffected. Contrary to the original theory suggesting that synonymous variants are benign, there is a growing volume of research showing that, despite their "silent" mechanism of action, some synonymous variation may be deleterious. Here, we studied the extent of the negative selective pressure acting on different classes of synonymous variants by analyzing the relative enrichment of synonymous singleton variants in the human exomes provided by gnomAD. Using a modification of the mutability-adjusted proportion of singletons (MAPS) metric as a measure of purifying selection, we found that some classes of synonymous variants are subject to stronger negative selection than others. For instance, variants that reduce codon optimality undergo stronger selection than optimality-increasing variants. Besides, selection affects synonymous variants implicated in splice-site-loss or splice-site-gain events. To understand what drives this negative selection, we tested a number of predictors in the aim to explain the variability in the selection scores. Our findings provide insights into the effects of synonymous variants at the population level, highlighting the specifics of the role that these variants play in health and disease.

ConanVarvar: a versatile tool for the detection of large syndromic copy number variation from whole-genome sequencing data.

BACKGROUND: A wide range of tools are available for the detection of copy number variants (CNVs) from whole-genome sequencing (WGS) data. However, none of them focus on clinically-relevant CNVs, such as those that are associated with known genetic syndromes. Such variants are often large in size, typically 1-5 Mb, but currently available CNV callers have been developed and benchmarked for the discovery of smaller variants. Thus, the ability of these programs to detect tens of real syndromic CNVs remains largely unknown. RESULTS: Here we present ConanVarvar, a tool which implements a complete workflow for the targeted analysis of large germline CNVs from WGS data. ConanVarvar comes with an intuitive R Shiny graphical user interface and annotates identified variants with information about 56 associated syndromic conditions. We benchmarked ConanVarvar and four other programs on a dataset containing real and simulated syndromic CNVs larger than 1 Mb. In comparison to other tools, ConanVarvar reports 10-30 times less false-positive variants without compromising sensitivity and is quicker to run, especially on large batches of samples. CONCLUSIONS: ConanVarvar is a useful instrument for primary analysis in disease sequencing studies, where large CNVs could be the cause of disease.

Whole genome sequencing in transposition of the great arteries and associations with clinically relevant heart, brain and laterality genes.

BACKGROUND: The most common cyanotic congenital heart disease (CHD) requiring management as a neonate is transposition of great arteries (TGA). Clinically, up to 50% of TGA patients develop some form of neurodevelopmental disability (NDD), thought to have a significant genetic component. A "ciliopathy" and links with laterality disorders have been proposed. This first report of whole genome sequencing in TGA, sought to identify clinically relevant variants contributing to heart, brain and laterality defects. METHODS: Initial whole genome sequencing analyses on 100 TGA patients focussed on established disease genes related to CHD (n = 107), NDD (n = 659) and heterotaxy (n = 74). Single variant as well as copy number variant analyses were conducted. Variant pathogenicity was assessed using the American College of Medical Genetics and Genomics-Association for Molecular Pathology guidelines. RESULTS: Fifty-five putatively damaging variants were identified in established disease genes associated with CHD, NDD and heterotaxy; however, no clinically relevant variants could be attributed to disease. Notably, case-control analyses identified significantly more predicted-damaging, silent and total variants in TGA cases than healthy controls in established CHD genes (P < .001), NDD genes (P < .001) as well as across the three gene panels (P < .001). CONCLUSION: We present compelling evidence that the majority of TGA is not caused by monogenic rare variants and is most likely oligogenic and/or polygenic in nature, highlighting the complex genetic architecture and multifactorial influences on this CHD sub-type and its long-term sequelae. Assessment of variant burden in key heart, brain and/or laterality genes may be required to unravel the genetic contributions to TGA and related disabilities.