Support for down-tuning of the calreticulin gene in the process of human evolution

Distinct structural bases for sequence-specific DNA binding by mammalian BEN domain proteins

The BEN domain is a recently recognized DNA binding module that is present in diverse metazoans and certain viruses. Several BEN domain factors are known as transcriptional repressors, but, overall, relatively little is known of how BEN factors identify their targets in humans. In particular, X-ray structures of BEN domain:DNA complexes are only known for Drosophila factors bearing a single BEN domain, which lack direct vertebrate orthologs. Here, we characterize several mammalian BEN domain (BD) factors, including from two NACC family BTB-BEN proteins and from BEND3, which has four BDs. In vitro selection data revealed sequence-specific binding activities of isolated BEN domains from all of these factors. We conducted detailed functional, genomic, and structural studies of BEND3. We show that BD4 is a major determinant for in vivo association and repression of endogenous BEND3 targets. We obtained a high-resolution structure of BEND3-BD4 bound to its preferred binding site, which reveals how BEND3 identifies cognate DNA targets and shows differences with one of its non-DNA-binding BEN domains (BD1). Finally, comparison with our previous invertebrate BEN structures, along with additional structural predictions using AlphaFold2 and RoseTTAFold, reveal distinct strategies for target DNA recognition by different types of BEN domain proteins. Together, these studies expand the DNA recognition activities of BEN factors and provide structural insights into sequence-specific DNA binding by mammalian BEN proteins.

Evolving Evidence of Calreticulin as a Pharmacological Target in Neurological Disorders

Calreticulin (CALR), a lectin-like ER chaperone, was initially known only for its housekeeping function, but today it is recognized for many versatile roles in different compartments of a cell. Apart from canonical roles in protein folding and calcium homeostasis, it performs a variety of noncanonical roles, mostly in CNS development. In the past, studies have linked Calreticulin with various other biological components which are detrimental in deciding the fate of neurons. Many neurological disorders that differ in their etiology are commonly associated with aberrant levels of Calreticulin, that lead to modulation of apoptosis and phagocytosis, and impact on transcriptional pathways, impairment in proteostatis, and calcium imbalances. Such multifaceted properties of Calreticulin are the reason why it has been implicated in vital roles of the nervous system in recent years. Hence, understanding its role in the physiology of neurons would help to unearth its involvement in the spectrum of neurological disorders. This Review aims toward exploring the interplay of Calreticulin in neurological disorders which would aid in targeting Calreticulin for developing novel neurotherapeutics.

Decreased gene expression activity as a result of a mutation in the calreticulin gene promoter in a family case of schizoaffective disorder

Accumulating evidence of population association studies support the hypothesis that the high heritability of major psychiatric disorders is a combination of relatively common alleles of modest effect, and rare alleles some with relatively larger effects. We have previously reported low frequency mutations in the proximal promoter of the human calreticulin (CALR) gene that co-occur with the spectrum of major psychiatric disorders. One of those mutations at -205C>T (rs556992558) was detected in an isolate case of schizoaffective disorder. In the current study, the functional implication of mutation -205T is studied in the human neuronal cell lines LAN-5, BE(2)-C and HEK-293. In contrast with other mutations in the promoter region which increase gene expression activity, the -205T mutation significantly decreased gene expression in those cell lines in comparison with the wild-type -205C nucleotide (p < 0.000001, p < 0.0005, and p < 0.017, respectively). Treatment of the cell lines with the mood-stabilizing drug, valproic acid (VPA) resulted in differential gene expression activity in the mutant -205T versus the wild-type -205C construct. VPA increased gene expression activity in both constructs, while a significantly higher expression activity was observed in the mutant construct (p < 0.01), indicative of the creation of a positive effector binding site for VPA as a result of the -205T mutation. We conclude that deviation from normalcy in the level of CALR in either direction is associated with major psychiatric disorders.

Genome-wide identification of human- and primate-specific core promoter short tandem repeats

Recent reports of a link between human- and primate-specific genetic factors and human/primate-specific characteristics and diseases necessitate genome-wide identification of those factors. We have previously reported core promoter short tandem repeats (STRs) of extreme length (≥6-repeats) that have expanded exceptionally in primates vs. non-primates, and may have a function in adaptive evolution. In the study reported here, we extended our study to the human STRs of ≥3-repeats in the category of penta and hexaucleotide STRs, across the entire human protein coding gene core promoters, and analyzed their status in several superorders and orders of vertebrates, using the Ensembl database. The ConSite software was used to identify the transcription factor (TF) sets binding to those STRs. STR specificity was observed at different levels of human and non-human primate (NHP) evolution. 73% of the pentanucleotide STRs and 68% of the hexanucleotide STRs were found to be specific to human and NHPs. AP-2alpha, Sp1, and MZF were the predominantly selected TFs (90%) binding to the human-specific STRs. Furthermore, the number of TF sets binding to a given STR was found to be a selection factor for that STR. Our findings indicate that selected STRs, the cognate binding TFs, and the number of TF set binding to those STRs function as switch codes at different levels of human and NHP evolution and speciation.

BEND3 is involved in the human-specific repression of calreticulin: Implication for the evolution of higher brain functions in human

Recent emerging evidence indicates that changes in gene expression levels are linked to human evolution. We have previously reported a human-specific nucleotide in the promoter sequence of the calreticulin (CALR) gene at position -220C, which is the site of action of valproic acid. Reversion of this nucleotide to the ancestral A-allele has been detected in patients with degrees of deficit in higher brain cognitive functions. This mutation has since been reported in the 1000 genomes database at an approximate frequency of <0.0004 in humans (rs138452745). In the study reported here, we present update on the status of rs138452745 across evolution, based on the Ensembl and NCBI databases. The DNA pulldown assay was also used to identify the proteins binding to the C- and A-alleles, using two cell lines, SK-N-BE and HeLa. Consistent with our previous findings, the C-allele is human-specific, and the A-allele is the rule across all other species (N=38). This nucleotide resides in a block of 12-nucleotides that is strictly conserved across evolution. The DNA pulldown experiments revealed that in both SK-N-BE and HeLa cells, the transcription repressor BEN domain containing 3 (BEND3) binds to the human-specific C-allele, whereas the nuclear factor I (NFI) family members, NF1A, B, C, and X, specifically bind to the ancestral A-allele. This binding pattern is consistent with a previously reported decreased promoter activity of the C-allele vs. the A-allele. We propose that there is a link between binding of BEND3 to the CALR rs138452745 C-allele and removal of NFI binding site from this nucleotide, and the evolution of human-specific higher brain functions. To our knowledge, CALR rs138452745 is the first instance of enormous nucleotide conservation across evolution, except in the human species.

An exceptionally long CA-repeat in the core promoter of SCGB2B2 links with the evolution of apes and Old World monkeys

We have recently reported a genome-scale catalog of human protein-coding genes that contain "exceptionally long" STRs (≥6-repeats) in their core promoter, which may be of selective advantage in this species. At the top of that list, SCGB2B2 (also known as SCGBL), contains one of the longest CA-repeat STRs identified in a human gene core promoter, at 25-repeats. In the study reported here, we analyzed the conservation status of this CA-STR across evolution. The functional implication of this STR to alter gene expression activity was also analyzed in the HEK-293 cell line. We report that the SCGB2B2 core promoter CA-repeat reaches exceptional lengths, ranging from 9- to 25-repeats, across Apes (Hominoids) and the Old World monkeys (CA>2-repeats were not detected in any other species). The longest CA-repeats and highest identity in the SCGB2B2 protein sequence were observed between human and bonobo. A trend for increased gene expression activity was observed from the shorter to the longer CA-repeats (p<0.009), and the CA-repeat increased gene expression activity, per se (p<0.02). We propose that the SCGB2B2 gene core promoter CA-repeat functions as an expression code for the evolution of Apes and the Old World monkeys.

Exceptionally long 5' UTR short tandem repeats specifically linked to primates

We have previously reported genome-scale short tandem repeats (STRs) in the core promoter interval (i.e. -120 to +1 to the transcription start site) of protein-coding genes that have evolved identically in primates vs. non-primates. Those STRs may function as evolutionary switch codes for primate speciation. In the current study, we used the Ensembl database to analyze the 5' untranslated region (5' UTR) between +1 and +60 of the transcription start site of the entire human protein-coding genes annotated in the GeneCards database, in order to identify "exceptionally long" STRs (≥5-repeats), which may be of selective/adaptive advantage. The importance of this critical interval is its function as core promoter, and its effect on transcription and translation. In order to minimize ascertainment bias, we analyzed the evolutionary status of the human 5' UTR STRs of ≥5-repeats in several species encompassing six major orders and superorders across mammals, including primates, rodents, Scandentia, Laurasiatheria, Afrotheria, and Xenarthra. We introduce primate-specific STRs, and STRs which have expanded from mouse to primates. Identical co-occurrence of the identified STRs of rare average frequency between 0.006 and 0.0001 in primates supports a role for those motifs in processes that diverged primates from other mammals, such as neuronal differentiation (e.g. APOD and FGF4), and craniofacial development (e.g. FILIP1L). A number of the identified STRs of ≥5-repeats may be human-specific (e.g. ZMYM3 and DAZAP1). Future work is warranted to examine the importance of the listed genes in primate/human evolution, development, and disease.

Dominant and Protective Role of the CYTH4 Primate-Specific GTTT-Repeat Longer Alleles Against Neurodegeneration

Primate-specific genes and regulatory mechanisms could provide insight into human brain functioning and disease. In a genome-scale analysis of the entire protein-coding genes listed in the GeneCards database, we have recently reported human genes that contain "exceptionally long" short tandem repeats (STRs) in their core promoter, which may be of adaptive/selective evolutionary advantage in this species. The longest tetra-nucleotide repeat identified in a human gene core promoter belongs to the CYTH4 gene. This GTTT-repeat is specific to Hominidae and Old World monkeys, and the shortest allele of this repeat, (GTTT)6, is linked with neural dysfunction and type I bipolar disorder in human. In the present study, we sought a possibly broader role for the CYTH4 gene core promoter GTTT-repeat in neural functioning and investigated its allelic distribution in a total of 949 human subjects, consisting of two neurodegenerative disorders, multiple sclerosis (MS) (n = 272) and Alzheimer's disease (AD) (n = 257), and controls (n = 420). The range of the alleles of this GTTT-repeat in the human sample studied was between 6- and 9-repeats. The shortest allele, (GTTT)6, was significantly in excess in the MS and AD patients in comparison with the controls (p < 0.004). The 6/6, 6/7, and 7/7 genotypes were in excess in the MS and AD patients, whereas the overall frequency of all other genotypes (consisting of at least one longer allele, i.e., 8- or 9-repeat) was higher in the controls (p < 0.005), indicating a dominant and protective effect for the longer alleles against neurodegeneration. This is the first indication of the involvement of a primate-specific STR in neurodegeneration in humans. We propose an adaptive evolutionary role for the expansion of the CYTH4 gene core promoter GTTT-repeat in the human brain, which is supported by a link between the shortest allele of this repeat with neuropsychiatric disorders.

A primate-specific functional GTTT-repeat in the core promoter of CYTH4 is linked to bipolar disorder in human

Evidence of primate-specific genes and gene regulatory mechanisms linked to bipolar disorder (BD) lend support to evolutionary/adaptive processes in the pathogenesis of this disorder. Following a genome-scale analysis of the entire protein coding genes annotated in the GeneCards database, we have recently reported that cytohesin-4 (CYTH4) contains the longest tetra-nucleotide short tandem repeat (STR) identified in a human protein-coding gene core promoter, which may be of adaptive advantage to this species. In the current study, we analyzed the evolutionary trend of this STR across evolution. We also analyzed the functional implication and distribution of this STR in a group of patients with type 1 BD (n=233) and controls (n=262). We found that this STR is exceptionally expanded in primates (Fisher exact p<0.00003). Association was observed between type I BD and the 6-repeat allele of this STR, (GTTT)₆ (Yates corrected Χ(2)=12.68, p<0.0001, OR: 1.68). This allele is the shortest length of the GTTT-repeat identified in the human subjects studied. Consistent with that finding, excess homozygosity was observed for the shorter alleles, (GTTT)₆ and (GTTT)₇, vs. the longer alleles, (GTTT)₈ and (GTTT)₉ in the BD group (Yates corrected Χ(2)=5.18, p<0.01, 1 df, OR: 1.96). Using Dual Glo luciferase system in HEK-293 cells, a trend for gene expression repression was observed from the 6- to the 9-repeat allele (p<0.003), and the GTTT-repeat significantly down-regulated gene expression, per se (p<0.0006). This is the first evidence of a link between a primate-specific STR and a major psychiatric disorder in human. It may be speculated that the CYTH4 GTTT-repeat in primates may have conferred selective advantage to this order, reflected in neural function and neurophenotypes. The role of the CYTH4 gene in the pathogenesis of type I BD remains to be clarified in the future studies.

Core promoter short tandem repeats as evolutionary switch codes for primate speciation

Alteration in gene expression levels underlies many of the phenotypic differences across species. Because of their highly mutable nature, proximity to the +1 transcription start site (TSS), and the emerging evidence of functional impact on gene expression, core promoter short tandem repeats (STRs) may be considered an ideal source of variation across species. In a genome-scale analysis of the entire Homo sapiens protein-coding genes, we have previously identified core promoters with at least one STR of ≥ 6-repeats, with possible selective advantage in this species. In the current study, we performed reverse analysis of the entire Homo sapiens orthologous genes in mouse in the Ensembl database, in order to identify conserved STRs that have shrunk as an evolutionary advantage to humans. Two protocols were used to minimize ascertainment bias. Firstly, two species sharing a more recent ancestor with Homo sapiens (i.e. Pan troglodytes and Gorilla gorilla gorilla) were also included in the study. Secondly, four non-primate species encompassing the major orders across Mammals, including Scandentia, Laurasiatheria, Afrotheria, and Xenarthra were analyzed as out-groups. We introduce STR evolutionary events specifically identical in primates (i.e. Homo sapiens, Pan troglodytes, and Gorilla gorilla gorilla) vs. non-primate out-groups. The average frequency of the identically shared STR motifs across those primates ranged between 0.00005 and 0.06. The identified genes are involved in important evolutionary and developmental processes, such as normal craniofacial development (TFAP2B), regulation of cell shape (PALMD), learning and long-term memory (RGS14), nervous system development (GFRA2), embryonic limb morphogenesis (PBX2), and forebrain development (APAF1). We provide evidence of core promoter STRs as evolutionary switch codes for primate speciation, and the first instance of identity-by-descent for those motifs at the interspecies level.

Exceptional expansion and conservation of a CT-repeat complex in the core promoter of PAXBP1 in primates

Adaptive evolution may be linked with the genomic distribution and function of short tandem repeats (STRs). Proximity of the core promoter STRs to the +1 transcription start site (TSS), and their mutable nature are characteristics that highlight those STRs as a novel source of interspecies variation. The PAXBP1 gene (alternatively known as GCFC1) core promoter contains the longest STR identified in a Homo sapiens gene core promoter. Indeed, this core promoter is a stretch of four consecutive CT-STRs. In the current study, we used the Ensembl, NCBI, and UCSC databases to analyze the evolutionary trend and functional implication of this CT-STR complex in six major lineages across vertebrates, including primates, non-primate mammals, birds, reptiles, amphibians, and fish. We observed exceptional expansion (≥4-repeats) and conservation of this CT-STR complex across primates, except prosimians, Microcebus murinus and Otolemur garnettii (Fisher exact P<4.1×10(-7)). H. sapiens has the most complex STR formula, and longest repeats. Macaca mulatta and Callithrix jacchus monkeys have the simplest STR formulas, and shortest repeat numbers. CT≥4-repeats were not detected in non-primate lineages. Different length alleles across the PAXBP1 core promoter CT-STRs significantly altered gene expression in vitro (P<0.001, t-test). PAXBP1 has a crucial role in craniofacial development, myogenesis, and spine morphogenesis, properties that have been diverged between primates and non-primates. To our knowledge, this is the first instance of expansion and conservation of a STR complex co-occurring specifically with the primate lineage.

Novel evidence of the involvement of calreticulin in major psychiatric disorders

Calreticulin (CALR) is a multi-functional protein that is strictly conserved across species. Two mRNA transcripts have been recognized for the CALR gene in humans, which use a common promoter sequence. We have recently reported mutations in the CALR promoter that co-occur with psychosis. One of those mutations at -220A increases gene expression in human BE(2)-C and HEK-293 cell lines. This mutation is the first instance of a functional cognition-deficit mutation reversing a human gene promoter to the primitive type. In the current study, we analyzed the effect of the most widely-used mood-stabilizing drug, valproic acid (VPA), on nucleotide -220 in two neuronal cell lines, LAN-5 and N2A. Remarkably, VPA increased gene expression in the cells with the wild-type -220C construct, whereas a dramatic decrease in gene expression was observed in the cell lines with the mutant construct (p<0.000004 and p<0.016, respectively). We also sequenced the 600-bp CALR promoter, and the highly conserved intron 1 sequence in an independent sample of patients afflicted with major psychiatric disorders and controls. A new case of major depressive disorder with psychotic features with the -220A mutation was identified. A novel 1-bp insertion was also detected in intron 1 at IVSI-310, in a case of amphetamine-induced psychosis. As for the psychosis-linked CALR promoter mutations identified to-date, the IVSI mutation was not detected in the control pool. This mutation creates a RREB-1 transcription factor binding site within the first intron. Our present findings identify the site of action of VPA in the CALR promoter, and introduce a novel mutation in a case of substance-induced psychosis in the first intron of CALR.