Human Single Nucleotide Variants Research Articles

Affinity-optimizing enhancer variants disrupt development

Enhancers control the location and timing of gene expression and contain the majority of variants associated with disease1–3. The ZRS is arguably the most well-studied vertebrate enhancer and mediates the expression of Shh in the developing limb4. Thirty-one human single-nucleotide variants (SNVs) within the ZRS are associated with polydactyly4–6. However, how this enhancer encodes tissue-specific activity, and the mechanisms by which SNVs alter the number of digits, are poorly understood. Here we show that the ETS sites within the ZRS are low affinity, and identify a functional ETS site, ETS-A, with extremely low affinity. Two human SNVs and a synthetic variant optimize the binding affinity of ETS-A subtly from 15% to around 25% relative to the strongest ETS binding sequence, and cause polydactyly with the same penetrance and severity. A greater increase in affinity results in phenotypes that are more penetrant and more severe. Affinity-optimizing SNVs in other ETS sites in the ZRS, as well as in ETS, interferon regulatory factor (IRF), HOX and activator protein 1 (AP-1) sites within a wide variety of enhancers, cause gain-of-function gene expression. The prevalence of binding sites with suboptimal affinity in enhancers creates a vulnerability in genomes whereby SNVs that optimize affinity, even slightly, can be pathogenic. Searching for affinity-optimizing SNVs in genomes could provide a mechanistic approach to identify causal variants that underlie enhanceropathies.

Increased mutation and gene conversion within human segmental duplications

Single-nucleotide variants (SNVs) in segmental duplications (SDs) have not been systematically assessed because of the limitations of mapping short-read sequencing data1,2. Here we constructed 1:1 unambiguous alignments spanning high-identity SDs across 102 human haplotypes and compared the pattern of SNVs between unique and duplicated regions3,4. We find that human SNVs are elevated 60% in SDs compared to unique regions and estimate that at least 23% of this increase is due to interlocus gene conversion (IGC) with up to 4.3 megabase pairs of SD sequence converted on average per human haplotype. We develop a genome-wide map of IGC donors and acceptors, including 498 acceptor and 454 donor hotspots affecting the exons of about 800 protein-coding genes. These include 171 genes that have ‘relocated’ on average 1.61 megabase pairs in a subset of human haplotypes. Using a coalescent framework, we show that SD regions are slightly evolutionarily older when compared to unique sequences, probably owing to IGC. SNVs in SDs, however, show a distinct mutational spectrum: a 27.1% increase in transversions that convert cytosine to guanine or the reverse across all triplet contexts and a 7.6% reduction in the frequency of CpG-associated mutations when compared to unique DNA. We reason that these distinct mutational properties help to maintain an overall higher GC content of SD DNA compared to that of unique DNA, probably driven by GC-biased conversion between paralogous sequences5,6.

The PRALINE database: protein and Rna humAn singLe nucleotIde variaNts in condEnsates.

Biological condensates are membraneless organelles with different material properties. Proteins and RNAs are the main components, but most of their interactions are still unknown. Here, we introduce PRALINE, a database for the interrogation of proteins and RNAs contained in stress granules, processing bodies and other assemblies including droplets and amyloids. PRALINE provides information about the predicted and experimentally validated protein-protein, protein-RNA and RNA-RNA interactions. For proteins, it reports the liquid-liquid phase separation and liquid-solid phase separation propensities. For RNAs, it provides information on predicted secondary structure content. PRALINE shows detailed information on human single-nucleotide variants, their clinical significance and presence in protein and RNA binding sites, and how they can affect condensates' physical properties. PRALINE is freely accessible on the web at http://praline.tartaglialab.com.

Discovery of human ACE2 variants with altered recognition by the SARS-CoV-2 spike protein.

Understanding how human ACE2 genetic variants differ in their recognition by SARS-CoV-2 can facilitate the leveraging of ACE2 as an axis for treating and preventing COVID-19. In this work, we experimentally interrogate thousands of ACE2 mutants to identify over one hundred human single-nucleotide variants (SNVs) that are likely to have altered recognition by the virus, and make the complementary discovery that ACE2 residues distant from the spike interface influence the ACE2-spike interaction. These findings illuminate new links between ACE2 sequence and spike recognition, and could find substantial utility in further fundamental research that augments epidemiological analyses and clinical trial design in the contexts of both existing strains of SARS-CoV-2 and novel variants that may arise in the future.

Role of a genetic variation in the microRNA-4421 binding site of ERP29 regarding risk of oropharynx cancer and prognosis

We conducted a two-stage association study on patients with oropharynx (OP) squamous cell carcinoma (SCC) and healthy controls to identify single nucleotide variants (SNVs) located at the microRNA (miR)-binding sites of carcinogenesis genes associated with risk and prognosis of the disease. In stage 1, 49 patients and 49 controls were analyzed using Genome-Wide Human SNV Arrays to identify variants in the 3′-untranslated region (3′-UTR) of carcinogenesis-related genes, and one SNV was selected for data validation in stage 2 by TaqMan assays in 250 OPSCC patients and 250 controls. The ERP29 c.*293A > G (rs7114) SNV located at miR-4421 binding site was selected for data validation among 46 SNVs. The ERp29 and miR-4421 levels were evaluated by quantitative-PCR and Western blotting. Interaction between miR-4421 with 3′-UTR of ERP29 was evaluated by luciferase reporter assay. Event-free survival (EFS) was calculated by Kaplan–Meier and Cox methods. ERP29 GG variant genotype was more common in OPSCC patients than in controls (6.4% vs 3.6%, p = 0.02; odds ratio: 5.67; 95% confidence interval (CI) 1.27–25.26). Shorter EFS were seen in the base of tongue (BT) SCC patients with GG genotype (0.0% vs 36.2%, p = 0.01; hazard ratio: 2.31; 95% CI: 1.03–5.15). Individuals with ERP29 AG or GG genotypes featured lower levels of ERP29 mRNA (p = 0.005), ERp29 protein (p < 0.001) and higher levels of miR-4421 (p = 0.02). The miR-4421 showed more efficient binding with 3′-UTR of the variant G allele when compared with wild-type allele A (p = 0.001). Our data suggest that ERP29 rs7114 SNV may alter the risk and prognosis of OPSCC due to variation in the ERp29 production possibly modulated by miR-4421.

Nucleotide sequence conservation of novel and established cis-regulatory sites within the tyrosine hydroxylase gene promoter.

Tyrosine hydroxylase (TH) is the rate-limiting enzyme in catecholamine biosynthesis and its gene proximal promoter ( < 1 kb upstream from the transcription start site) is essential for regulating transcription in both the developing and adult nervous systems. Several putative regulatory elements within the TH proximal promoter have been reported, but evolutionary conservation of these elements has not been thoroughly investigated. Since many vertebrate species are used to model development, function and disorders of human catecholaminergic neurons, identifying evolutionarily conserved transcription regulatory mechanisms is a high priority. In this study, we align TH proximal promoter nucleotide sequences from several vertebrate species to identify evolutionarily conserved motifs. This analysis identified three elements (a TATA box, cyclic AMP response element (CRE) and a 5'-GGTGG-3' site) that constitute the core of an ancient vertebrate TH promoter. Focusing on only eutherian mammals, two regions of high conservation within the proximal promoter were identified: a ∼250 bp region adjacent to the transcription start site and a ∼85 bp region located approximately 350 bp further upstream. Within both regions, conservation of previously reported cis-regulatory motifs and human single nucleotide variants was evaluated. Transcription reporter assays in a TH -expressing cell line demonstrated the functionality of highly conserved motifs in the proximal promoter regions and electromobility shift assays showed that brain-region specific complexes assemble on these motifs. These studies also identified a non-canonical CRE binding (CREB) protein recognition element in the proximal promoter. Together, these studies provide a detailed analysis of evolutionary conservation within the TH promoter and identify potential cis-regulatory motifs that underlie a core set of regulatory mechanisms in mammals.

A general framework for estimating the relative pathogenicity of human genetic variants

Our capacity to sequence human genomes has exceeded our ability to interpret genetic variation. Current genomic annotations tend to exploit a single information type (e.g. conservation) and/or are restricted in scope (e.g. to missense changes). Here, we describe Combined Annotation Dependent Depletion (CADD), a framework that objectively integrates many diverse annotations into a single, quantitative score. We implement CADD as a support vector machine trained to differentiate 14.7 million high-frequency human derived alleles from 14.7 million simulated variants. We pre-compute “C-scores” for all 8.6 billion possible human single nucleotide variants and enable scoring of short insertions/deletions. C-scores correlate with allelic diversity, annotations of functionality, pathogenicity, disease severity, experimentally measured regulatory effects, and complex trait associations, and highly rank known pathogenic variants within individual genomes. The ability of CADD to prioritize functional, deleterious, and pathogenic variants across many functional categories, effect sizes and genetic architectures is unmatched by any current annotation.

Base-pair neutral homozygotes can be discriminated by calibrated high-resolution melting of small amplicons

Genotyping by high-resolution melting analysis of small amplicons is homogeneous and simple. However, this approach can be limited by physical and chemical components of the system that contribute to intersample melting variation. It is challenging for this method to distinguish homozygous G::C from C::G or A::T from T::A base-pair neutral variants, which comprise ∼16% of all human single nucleotide polymorphisms (SNPs). We used internal oligonucleotide calibrators and custom analysis software to improve small amplicon (42–86 bp) genotyping on the LightScanner®. Three G/C (PAH c.1155C>G, CHK2 c.1-3850G>C and candidate gene BX647987 c.261+22,290C>G) and three T/A (CPS1 c.3405-29A>T, OTC c.299-8T>A and MSH2 c.1511-9A>T) human single nucleotide variants were analyzed. Calibration improved homozygote genotyping accuracy from 91.7 to 99.7% across 1105 amplicons from 141 samples for five of the six targets. The average Tm standard deviations of these targets decreased from 0.067°C before calibration to 0.022°C after calibration. We were unable to generate a small amplicon that could discriminate the BX647987 c.261+22,290C>G (rs1869458) SNP, despite reducing standard deviations from 0.086°C to 0.032°C. Two of the sites contained symmetric nearest neighbors adjacent to the SNPs. Unexpectedly, we were able to distinguish these homozygotes by Tm even though current nearest neighbor models predict that the two homozygous alleles would be identical.