Genome-wide identification of dominant polyadenylation hexamers for use in variant classification.

Polyadenylation is an essential process for the stabilization and export of mRNAs to the cytoplasm and the polyadenylation signal hexamer (herein referred to as hexamer) plays a key role in this process. Yet, only 14 Mendelian disorders have been associated with hexamer variants. This is likely an under-ascertainment as hexamers are not well defined and not routinely examined in molecular analysis. To facilitate the interrogation of putatively pathogenic hexamer variants, we set out to define functionally important hexamers genome-wide as a resource for research and clinical testing interrogation. We identified predominant polyA sites (herein referred to as pPAS) and putative predominant hexamers across protein coding genes (PAS usage >50% per gene). As a measure of the validity of these sites, the population constraint of 4532 predominant hexamers were measured. The predominant hexamers had fewer observed variants compared to non-predominant hexamers and trimer controls, and CADD scores for variants in these hexamers were significantly higher than controls. Exome data for 1477 individuals were interrogated for hexamer variants and transcriptome data were generated for 76 individuals with 65 variants in predominant hexamers. 3' RNA-seq data showed these variants resulted in alternate polyadenylation events (38%) and in elongated mRNA transcripts (12%). Our list of pPAS and predominant hexamers are available in the UCSC genome browser and on GitHub. We suggest this list of predominant hexamers can be used to interrogate exome and genome data. Variants in these predominant hexamers should be considered candidates for pathogenic variation in human disease, and to that end we suggest pathogenicity criteria for classifying hexamer variants.

Low-level variant calling for non-matched samples using a position-based and nucleotide-specific approach.

BACKGROUND: The widespread use of next-generation sequencing has identified an important role for somatic mosaicism in many diseases. However, detecting low-level mosaic variants from next-generation sequencing data remains challenging. RESULTS: Here, we present a method for Position-Based Variant Identification (PBVI) that uses empirically-derived distributions of alternate nucleotides from a control dataset. We modeled this approach on 11 segmental overgrowth genes. We show that this method improves detection of single nucleotide mosaic variants of 0.01-0.05 variant allele fraction compared to other low-level variant callers. At depths of 600 × and 1200 ×, we observed > 85% and > 95% sensitivity, respectively. In a cohort of 26 individuals with somatic overgrowth disorders PBVI showed improved signal to noise, identifying pathogenic variants in 17 individuals. CONCLUSION: PBVI can facilitate identification of low-level mosaic variants thus increasing the utility of next-generation sequencing data for research and diagnostic purposes.

A mouse model of Proteus syndrome.

Proteus syndrome is a mosaic, progressive overgrowth disorder caused by a somatic activating variant c.49G > A p.(E17K) in AKT1. The presentation in affected individuals is variable, with a diversity of tissues demonstrating abnormalities. Common manifestations include skin and bony overgrowth, vascular malformations (VMs), cysts and benign tumors. We used two methods to create mouse models that had endogenously-regulated mosaic expression of the Proteus syndrome variant. Variant allele fractions (VAFs) ranged from 0% to 50% across numerous tissues in 44 Proteus syndrome mice. Mice were phenotypically heterogeneous with lesions rarely observed before 12 months of age. VMs were the most frequent finding with a total of 69 found in 29 of 44 Proteus syndrome mice. Twenty-eight cysts and ectasia, frequently biliary, were seen in 22 of 44 Proteus syndrome mice. Varying levels of mammary hyperplasia were seen in 10 of 16 female Proteus syndrome mice with other localized regions of hyperplasia and stromal expansion noted in several additional animals. Interestingly, 27 of 31 Proteus syndrome animals had non-zero blood VAF that is in contrast to the human disorder where it is rarely seen in peripheral blood. Identification of variant-positive cells by green fluorescent protein (GFP) staining in chimeric Proteus syndrome mice showed that in some lesions, hyperplastic cells were predominantly GFP/Akt1E17K-positive and showed increased pAKT signal compared to GFP-negative cells. However, hyperplastic mammary epithelium was a mixture of GFP/Akt1E17K-positive and negative cells with some GFP/Akt1E17K-negative cells also having increased pAKT signal suggesting that the variant-positive cells can induce lesion formation in a non-cell autonomous manner.