DNA binding analysis of rare variants in homeodomains reveals homeodomain specificity-determining residues.

Homeodomains (HDs) are the second largest class of DNA binding domains (DBDs) among eukaryotic sequence-specific transcription factors (TFs) and are the TF structural class with the largest number of disease-associated mutations in the Human Gene Mutation Database (HGMD). Despite numerous structural studies and large-scale analyses of HD DNA binding specificity, HD-DNA recognition is still not fully understood. Here, we analyze 92 human HD mutants, including disease-associated variants and variants of uncertain significance (VUS), for their effects on DNA binding activity. Many of the variants alter DNA binding affinity and/or specificity. Detailed biochemical analysis and structural modeling identifies 14 previously unknown specificity-determining positions, 5 of which do not contact DNA. The same missense substitution at analogous positions within different HDs often exhibits different effects on DNA binding activity. Variant effect prediction tools perform moderately well in distinguishing variants with altered DNA binding affinity, but poorly in identifying those with altered binding specificity. Our results highlight the need for biochemical assays of TF coding variants and prioritize dozens of variants for further investigations into their pathogenicity and the development of clinical diagnostics and precision therapies.

Single-nucleotide variant calling in single-cell sequencing data with Monopogen.

Single-cell omics technologies enable molecular characterization of diverse cell types and states, but how the resulting transcriptional and epigenetic profiles depend on the cell's genetic background remains understudied. We describe Monopogen, a computational tool to detect single-nucleotide variants (SNVs) from single-cell sequencing data. Monopogen leverages linkage disequilibrium from external reference panels to identify germline SNVs and detects putative somatic SNVs using allele cosegregating patterns at the cell population level. It can identify 100 K to 3 M germline SNVs achieving a genotyping accuracy of 95%, together with hundreds of putative somatic SNVs. Monopogen-derived genotypes enable global and local ancestry inference and identification of admixed samples. It identifies variants associated with cardiomyocyte metabolic levels and epigenomic programs. It also improves putative somatic SNV detection that enables clonal lineage tracing in primary human clonal hematopoiesis. Monopogen brings together population genetics, cell lineage tracing and single-cell omics to uncover genetic determinants of cellular processes.

Common variants in signaling transcription-factor-binding sites drive phenotypic variability in red blood cell traits.

Genome-wide association studies identify genomic variants associated with human traits and diseases. Most trait-associated variants are located within cell-type-specific enhancers, but the molecular mechanisms governing phenotypic variation are less well understood. Here, we show that many enhancer variants associated with red blood cell (RBC) traits map to enhancers that are co-bound by lineage-specific master transcription factors (MTFs) and signaling transcription factors (STFs) responsive to extracellular signals. The majority of enhancer variants reside on STF and not MTF motifs, perturbing DNA binding by various STFs (BMP/TGF-ß-directed SMADs or WNT-induced TCFs) and affecting target gene expression. Analyses of engineered human blood cells and expression quantitative trait loci verify that disrupted STF binding leads to altered gene expression. Our results propose that the majority of the RBC-trait-associated variants that reside on transcription-factor-binding sequences fall in STF target sequences, suggesting that the phenotypic variation of RBC traits could stem from altered responsiveness to extracellular stimuli.

Transcription factor-DNA binding: beyond binding site motifs.

Sequence-specific transcription factors (TFs) regulate gene expression by binding to cis-regulatory elements in promoter and enhancer DNA. While studies of TF-DNA binding have focused on TFs' intrinsic preferences for primary nucleotide sequence motifs, recent studies have elucidated additional layers of complexity that modulate TF-DNA binding. In this review, we discuss technological developments for identifying TF binding preferences and highlight recent discoveries that elaborate how TF interactions, local DNA structure, and genomic features influence TF-DNA binding. We highlight novel approaches for characterizing functional binding site motifs that promise to inform our understanding of how TF binding controls gene expression and ultimately contributes to phenotype.

Survey of variation in human transcription factors reveals prevalent DNA binding changes.

Sequencing of exomes and genomes has revealed abundant genetic variation affecting the coding sequences of human transcription factors (TFs), but the consequences of such variation remain largely unexplored. We developed a computational, structure-based approach to evaluate TF variants for their impact on DNA binding activity and used universal protein-binding microarrays to assay sequence-specific DNA binding activity across 41 reference and 117 variant alleles found in individuals of diverse ancestries and families with Mendelian diseases. We found 77 variants in 28 genes that affect DNA binding affinity or specificity and identified thousands of rare alleles likely to alter the DNA binding activity of human sequence-specific TFs. Our results suggest that most individuals have unique repertoires of TF DNA binding activities, which may contribute to phenotypic variation.

HMGB1 binds to the rs7903146 locus in TCF7L2 in human pancreatic islets.

The intronic SNP rs7903146 in the T-cell factor 7-like 2 gene (TCF7L2) is the common genetic variant most highly associated with Type 2 diabetes known to date. The risk T-allele is located in an open chromatin region specific to human pancreatic islets of Langerhans, thereby accessible for binding of regulatory proteins. The risk T-allele locus exhibits stronger enhancer activity compared to the non-risk C-allele. The aim of this study was to identify transcriptional regulators that bind the open chromatin region in the rs7903146 locus and thereby potentially regulate TCF7L2 expression and activity. Using affinity chromatography followed by Edman sequencing, we identified one candidate regulatory protein, i.e. high-mobility group protein B1 (HMGB1). The binding of HMGB1 to the rs7903146 locus was confirmed in pancreatic islets from human deceased donors, in HCT116 and in HEK293 cell lines using: (i) protein purification on affinity columns followed by Western blot, (ii) chromatin immunoprecipitation followed by qPCR and (iii) electrophoretic mobility shift assay. The results also suggested that HMGB1 might have higher binding affinity to the C-allele of rs7903146 compared to the T-allele, which was supported in vitro using Dynamic Light Scattering, possibly in a tissue-specific manner. The functional consequence of HMGB1 depletion in HCT116 and INS1 cells was reduced insulin and TCF7L2 mRNA expression, TCF7L2 transcriptional activity and glucose stimulated insulin secretion. These findings suggest that the rs7903146 locus might exert its enhancer function by interacting with HMGB1 in an allele dependent manner.

Stable Formation of Gold Nanoparticles onto Redox-Active Solid Biosubstrates Made of Squid Suckerin Proteins.

The use of biomolecules to synthesize inorganic nanomaterials, including metallic nanoparticles, offers the ability to induce controlled growth under mild environmental conditions. Here, recently discovered silk-like "suckerin" proteins are used to induce the formation of gold nanoparticles (AuNPs). Advantage is taken of the distinctive biological and physico-chemical characteristics of suckerins, namely their facile recombinant expression, their solubility in aqueous solutions, and their modular primary structure with high molar content of redox-active tyrosine (Tyr) residues to induce the formation of AuNPs not only in solution, but also from nanostructured solid substrates fabricated from suckerins.

Functional VEGFA knockdown with artificial 3'-tailed mirtrons defined by 5' splice site and branch point.

Mirtrons are introns that form pre-miRNA hairpins after splicing to produce RNA interference (RNAi) effectors distinct from Drosha-dependent intronic miRNAs, and will be especially useful for co-delivery of coding genes and RNAi. A specific family of mirtrons - 3'-tailed mirtrons - has hairpins precisely defined on the 5' end by the 5' splice site and 3' end by the branch point. Here, we present design principles for artificial 3'-tailed mirtrons and demonstrate, for the first time, efficient gene knockdown with tailed mirtrons within eGFP coding region. These artificial tailed mirtrons, unlike canonical mirtrons, have very few sequence design restrictions. Tailed mirtrons targeted against VEGFA mRNA, the misregulation of which is causative of several disorders including cancer, achieved significant levels of gene knockdown. Tailed mirtron-mediated knockdown was further shown to be splicing-dependent, and at least as effective as equivalent artificial intronic miRNAs, with the added advantage of very defined cleavage sites for generation of mature miRNA guide strands. Further development and exploitation of this unique mirtron biogenesis pathway for therapeutic RNAi coupled into protein-expressing genes can potentially enable the development of precisely controlled combinatorial gene therapy.