Variable Glutamine-Rich Repeats Modulate Transcription Factor Activity

Polyglutamine expansion induced dynamic misfolding of androgen receptor

Spinal bulbar muscular atrophy (SBMA) is caused by a polyglutamine expansion (pQe) in the N-terminal transactivation domain of the human androgen receptor (AR-NTD), resulting in a combination of toxic gain- and loss-of-function mechanisms. The structural basis of these processes has not been resolved due to the disordered nature of the NTD, which hinders experimental analyses of its detailed conformations. Here, using extensive computational modeling, we show that AR-NTD forms dynamic compact regions, which upon pQe re-organize dynamically, mediated partly by direct pQ interaction with the Androgen N-Terminal Signature (ANTS) motif. The altered dynamics of the NTD result in a perturbation of interdomain interactions, with potential implications for the binding of the receptor protein to its response element. Oligomeric aggregation of the dynamic misfolded NTD exposes pQe, but blocks tau-5 and the FQNLF motif, which could lead to aberrant receptor transcriptional activity. These observations suggest a structural mechanism for AR dysfunction in SBMA.

DNA binding and mitotic phosphorylation protect polyglutamine proteins from assembly formation

Polyglutamine (polyQ) expansion is associated with pathogenic protein aggregation in neurodegenerative disorders. However, long polyQ tracts are also found in many transcription factors (TFs), such as FOXP2, a TF implicated in human speech. Here, we explore how FOXP2 and other glutamine-rich TFs avoid unscheduled assembly. Throughout interphase, DNA binding, irrespective of sequence specificity, has a solubilizing effect. During mitosis, multiple phosphorylation events promote FOXP2's eviction from chromatin and supplant the solubilizing function of DNA. Further, human-specific amino acid substitutions linked to the evolution of speech map to a mitotic phospho-patch, the "EVO patch," and reduce the propensity of the human FOXP2 to assemble. Fusing the pathogenic form of Huntingtin to either a DNA-binding domain, a phosphomimetic variant of this EVO patch, or a negatively charged peptide is sufficient to diminish assembly formation, suggesting that hijacking mechanisms governing solubility of glutamine-rich TFs may offer new strategies for treatment of polyQ expansion diseases.

PolyQ Length of the Clock Gene Is Correlated With Pelagic Larval Duration in the Damselfishes (Pomacentridae), but Within a Species Habitat Availability Counts

Immediately after hatching, larvae of coral reef fish leave their natal reef environment and begin their pelagic larval phase probably to avoid high predation on the reef. The time they spend in the open ocean (pelagic larval duration, PLD), before settlement varies from species to species and depends partly on developmental processes that eventually require re-settlement to a reef. The polyglutamine region (PolyQ) as part of the clock gene has been suggested as a possible candidate that could control developmental processes and potentially the time until settlement, which can be determined by counting the rings of the otoliths. We studied the potential relationship between the number of glutamine repeats in 20 species of pomacentrids and their PLDs. Most pomacentrids came from similar locations, so we avoided the impact of latitudinal clines on PLD. Within the clock gene, we found two main distinguishable, variable glutamine-rich regions (PolyQ and Qrich). Considering phylogenetic relationships, PolyQ/Qrich repeat length and pelagic larval duration were significantly positively correlated. However, when analyzing this relationship in a single species, the neon damselfish (Pomacentrus coelestis), we did not find a significant correlation between PolyQ length variation and PLD. Instead, we found a significant reduction of PLD in years with increased habitat availability. Our results show that glutamine-rich regions can influence the timing of settlement on a broader scale, but that ecological factors-such as habitat availability-can also have a significant impact.

Genome-wide profiling of polymorphic short tandem repeats and their influence on gene expression and trait variation in diverse rice populations

Short tandem repeats (STRs) modulate gene expression and contribute to trait variation. However, a systematic evaluation of the genomic characteristics of STRs has not been conducted, and their influence on gene expression in rice remains unclear. Here, we construct a map of 137,629 polymorphic STRs in the rice (Oryza sativa L.) genome using a population-scale resequencing dataset. A genome-wide survey encompassing 4726 accessions shows that the occurrence frequency, mutational patterns, chromosomal distribution, and functional properties of STRs are correlated with the sequences and lengths of repeat motifs. Leveraging a transcriptome dataset from 127 rice accessions, we identify 44,672 expression STRs (eSTRs) by modeling gene expression in response to the length variation of STRs. These eSTRs are notably enriched in the regulatory regions of genes with active transcriptional signatures. Population analysis identifies numerous STRs that have undergone genetic divergence among different rice groups and 1726 tagged STRs that may be associated with agronomic traits. By editing the (ACT)7 STR in OsFD1 promoter, we further experimentally validate its role in regulating gene expression and phenotype. Our study highlights the contribution of STRs to transcriptional regulation in plants and establishes the foundation for their potential use as alternative targets for genetic improvement.

The Role of Med15 Sequence Features in Transcription Factor Interactions

Med15 is a general transcriptional regulator and tail module subunit within the RNA Pol II mediator complex. The Saccharomyces cerevisiae Med15 protein has a well-structured N-terminal KIX domain, three activator binding domains (ABDs) and several naturally variable polyglutamine (poly-Q) tracts (Q1, Q2, Q3) embedded in an intrinsically disordered central region, and a C-terminal mediator association domain (MAD). We investigated how the presence of ABDs and changes in length and composition of poly-Q tracts influences Med15 activity using phenotypic, gene expression, transcription factor interaction and phase separation assays of truncation, deletion, and synthetic alleles. We found that individual Med15 activities were influenced by the number of activator binding domains (ABDs) and adjacent polyglutamine tract composition. Robust Med15 activity required at least the Q1 tract and the length of that tract modulated activity in a context-dependent manner. Reduced Msn2-dependent transcriptional activation due to Med15 Q1 tract variation correlated with reduced Msn2:Med15 interaction strength, but interaction strength did not always mirror phase separation propensity. We also observed that distant glutamine tracts and Med15 phosphorylation affected the activities of the KIX domain, and interaction studies revealed that intramolecular interactions may affect some Med15-transcription factor interactions.

Intrinsically Disordered Compositional Bias in Proteins: Sequence Traits, Region Clustering, and Generation of Hypothetical Functional Associations

Compositionally biased regions (CBRs), ie, tracts that are dominated by a subset of residue types, are common features of eukaryotic proteins. These are often found bounded within or almost coterminous with intrinsically disordered or 'natively unfolded' parts. Here, it is investigated how the function of such intrinsically disordered compositionally biased regions (ID-CBRs) is directly linked to their compositional traits, focusing on the well-characterized yeast (Saccharomyces cerevisiae) proteome as a test case. The ID-CBRs that are clustered together using compositional distance are discovered to have clear functional linkages at various levels of diversity. The specific case of the Sup35p and Rnq1p proteins that underlie causally linked prion phenomena ([PSI+] and [RNQ+]) is highlighted. Their prion-forming ID-CBRs are typically clustered very close together indicating some compositional engendering for [RNQ+] seeding of [PSI+] prions. Delving further, ID-CBRs with distinct types of residue patterning such as 'blocking' or relative segregation of residues into homopeptides are found to have significant functional trends. Specific examples of such ID-CBR functional linkages that are discussed are: Q/N-rich ID-CBRs linked to transcriptional coactivation, S-rich to transcription-factor binding, R-rich to DNA-binding, S/E-rich to protein localization, and D-rich linked to chromatin remodelling. These data may be useful in informing experimental hypotheses for proteins containing such regions.

The molecular basis for cellular function of intrinsically disordered protein regions

Intrinsically disordered protein regions exist in a collection of dynamic interconverting conformations that lack a stable 3D structure. These regions are structurally heterogeneous, ubiquitous and found across all kingdoms of life. Despite the absence of a defined 3D structure, disordered regions are essential for cellular processes ranging from transcriptional control and cell signalling to subcellular organization. Through their conformational malleability and adaptability, disordered regions extend the repertoire of macromolecular interactions and are readily tunable by their structural and chemical context, making them ideal responders to regulatory cues. Recent work has led to major advances in understanding the link between protein sequence and conformational behaviour in disordered regions, yet the link between sequence and molecular function is less well defined. Here we consider the biochemical and biophysical foundations that underlie how and why disordered regions can engage in productive cellular functions, provide examples of emerging concepts and discuss how protein disorder contributes to intracellular information processing and regulation of cellular function.

The Potential Role of Genic-SSRs in Driving Ecological Adaptation Diversity in Caragana Plants

Caragana, a xerophytic shrub genus widely distributed in northern China, exhibits distinctive geographical substitution patterns and ecological adaptation diversity. This study employed transcriptome sequencing technology to investigate 12 Caragana species, aiming to explore genic-SSR variations in the Caragana transcriptome and identify their role as a driving force for environmental adaptation within the genus. A total of 3666 polymorphic genic-SSRs were identified across different species. The impact of these variations on the expression of related genes was analyzed, revealing a significant linear correlation (p < 0.05) between the length variation of 264 polymorphic genic-SSRs and the expression of associated genes. Additionally, 2424 polymorphic genic-SSRs were located in differentially expressed genes among Caragana species. Through weighted gene co-expression network analysis, the expressions of these genes were correlated with 19 climatic factors and 16 plant functional traits in various habitats. This approach facilitated the identification of biological processes associated with habitat adaptations in the studied Caragana species. Fifty-five core genes related to functional traits and climatic factors were identified, including various transcription factors such as MYB, TCP, ARF, and structural proteins like HSP90, elongation factor TS, and HECT. The roles of these genes in the ecological adaptation diversity of Caragana were discussed. Our study identified specific genomic components and genes in Caragana plants responsive to heterogeneous habitats. The results contribute to advancements in the molecular understanding of their ecological adaptation, lay a foundation for the conservation and development of Caragana germplasm resources, and provide a scientific basis for plant adaptation to global climate change.

Structure-function relationships in protein homorepeats

Homorepeats (or polyX), protein segments containing repetitions of the same amino acid, are abundant in proteomes from all kingdoms of life and are involved in crucial biological functions as well as several neurodegenerative and developmental diseases. Mainly inserted in disordered segments of proteins, the structure/function relationships of homorepeats remain largely unexplored. In this review, we summarize present knowledge for the most abundant homorepeats, highlighting the role of the inherent structure and the conformational influence exerted by their flanking regions. Recent experimental and computational methods enable residue-specific investigations of these regions and promise novel structural and dynamic information for this elusive group of proteins. This information should increase our knowledge about the structural bases of phenomena such as liquid-liquid phase separation and trinucleotide repeat disorders.

PolyQ length-based molecular encoding of vocalization frequency in FOXP2

The transcription factor FOXP2, a regulator of vocalization- and speech/language-related phenotypes, contains two long polyQ repeats (Q1 and Q2) displaying marked, still enigmatic length variation across mammals. We found that the Q1/Q2 length ratio quantitatively encodes vocalization frequency ranges, from the infrasonic to the ultrasonic, displaying striking convergent evolution patterns. Thus, species emitting ultrasonic vocalizations converge with bats in having a low ratio, whereas species vocalizing in the low-frequency/infrasonic range converge with elephants and whales, which have higher ratios. Similar, taxon-specific patterns were observed for the FOXP2-related protein FOXP1. At the molecular level, we observed that the FOXP2 polyQ tracts form coiled coils, assembling into condensates and fibrils, and drive liquid-liquid phase separation (LLPS). By integrating evolutionary and molecular analyses, we found that polyQ length variation related to vocalization frequency impacts FOXP2 structure, LLPS, and transcriptional activity, thus defining a novel form of polyQ length-based molecular encoding of vocalization frequency.

Teleost genomic repeat landscapes in light of diversification rates and ecology

Repetitive DNA make up a considerable fraction of most eukaryotic genomes. In fish, transposable element (TE) activity has coincided with rapid species diversification. Here, we annotated the repetitive content in 100 genome assemblies, covering the major branches of the diverse lineage of teleost fish. We investigated if TE content correlates with family level net diversification rates and found support for a weak negative correlation. Further, we demonstrated that TE proportion correlates with genome size, but not to the proportion of short tandem repeats (STRs), which implies independent evolutionary paths. Marine and freshwater fish had large differences in STR content, with the most extreme propagation detected in the genomes of codfish species and Atlantic herring. Such a high density of STRs is likely to increase the mutational load, which we propose could be counterbalanced by high fecundity as seen in codfishes and herring.

Proteins with amino acid repeats constitute a rapidly evolvable and human-specific essentialome

Protein products of essential genes, indispensable for organismal survival, are highly conserved and bring about fundamental functions. Interestingly, proteins that contain amino acid homorepeats that tend to evolve rapidly are enriched in eukaryotic essentialomes. Why are proteins with hypermutable homorepeats enriched in conserved and functionally vital essential proteins? We solve this functional versus evolutionary paradox by demonstrating that human essential proteins with homorepeats bring about crosstalk across biological processes through high interactability and have distinct regulatory functions affecting expansive global regulation. Importantly, essential proteins with homorepeats rapidly diverge with the amino acid substitutions frequently affecting functional sites, likely facilitating rapid adaptability. Strikingly, essential proteins with homorepeats influence human-specific embryonic and brain development, implying that the presence of homorepeats could contribute to the emergence of human-specific processes. Thus, we propose that homorepeat-containing essential proteins affecting species-specific traits can be potential intervention targets across pathologies, including cancers and neurological disorders.

Functional coordination between transcription factor clustering and gene activity

The prevailing view of metazoan gene regulation is that transcription is facilitated through the formation of static activator complexes at distal regulatory regions. Here, we employed quantitative single-cell live-imaging and computational analysis to provide evidence that the dynamic assembly and disassembly process of transcription factor (TF) clusters at enhancers is a major source of transcriptional bursting in developing Drosophila embryos. We further show that the regulatory connectivity between TF clustering and burst induction is highly regulated through intrinsically disordered regions (IDRs). Addition of a poly-glutamine tract to the maternal morphogen Bicoid demonstrated that extended IDR length leads to ectopic TF clustering and burst induction from its endogenous target genes, resulting in defects in body segmentation during embryogenesis. Moreover, we successfully visualized the presence of "shared" TF clusters during the co-activation of two distant genes, which provides a concrete molecular explanation for the newly proposed "topological operon" hypothesis in metazoan gene regulation.

Adaptive protein evolution through length variation of short tandem repeats in Arabidopsis

Intrinsically disordered protein regions are of high importance for biotic and abiotic stress responses in plants. Tracts of identical amino acids accumulate in these regions and can vary in length over generations because of expansions and retractions of short tandem repeats at the genomic level. However, little attention has been paid to what extent length variation is shaped by natural selection. By environmental association analysis on 2514 length variable tracts in 770 whole-genome sequenced Arabidopsis thaliana, we show that length variation in glutamine and asparagine amino acid homopolymers, as well as in interaction hotspots, correlate with local bioclimatic habitat. We determined experimentally that the promoter activity of a light-stress gene depended on polyglutamine length variants in a disordered transcription factor. Our results show that length variations affect protein function and are likely adaptive. Length variants modulating protein function at a global genomic scale has implications for understanding protein evolution and eco-evolutionary biology.

Genome-wide survey of D/E repeats in human proteins uncovers their instability and aids in identifying their role in the chromatin regulator ATAD2

D/E repeats are stretches of aspartic and/or glutamic acid residues found in over 150 human proteins. We examined genomic stability of D/E repeats and functional characteristics of D/E repeat-containing proteins vis-à-vis the proteins with poly-Q or poly-A repeats, which are known to undergo pathologic expansions. Mining of tumor sequencing data revealed that D/E repeat-coding regions are similar to those coding poly-Qs and poly-As in increased incidence of trinucleotide insertions/deletions but differ in types and incidence of substitutions. D/E repeat-containing proteins preferentially function in chromatin metabolism and are the more likely to be nuclear and interact with core histones, the longer their repeats are. One of the longest D/E repeats of unknown function is in ATAD2, a bromodomain family ATPase frequently overexpressed in tumors. We demonstrate that D/E repeat deletion in ATAD2 suppresses its binding to nascent and mature chromatin and to the constitutive pericentromeric heterochromatin, where ATAD2 represses satellite transcription.

•

Many human proteins contain runs of aspartic/glutamic acid residues (D/E repeats)

•

D/E repeats show increased incidence of in-frame insertions/deletions in tumors

•

Nuclear and histone-interacting proteins often have long D/E repeats

•

D/E repeat of the oncogene ATAD2 controls its binding to pericentric chromatin

The evolution and polymorphism of mono-amino acid repeats in androgen receptor and their regulatory role in health and disease

Androgen receptor (AR) is a key member of nuclear hormone receptors with the longest intrinsically disordered N-terminal domain (NTD) in its protein family. There are four mono-amino acid repeats (polyQ1, polyQ2, polyG, and polyP) located within its NTD, of which two are polymorphic (polyQ1 and polyG). The length of both polymorphic repeats shows clinically important correlations with disease, especially with cancer and neurodegenerative diseases, as shorter and longer alleles exhibit significant differences in expression, activity and solubility. Importantly, AR has also been shown to undergo condensation in the nucleus by liquid-liquid phase separation, a process highly sensitive to protein solubility and concentration. Nonetheless, in prostate cancer cells, AR variants also partition into transcriptional condensates, which have been shown to alter the expression of target gene products. In this review, we summarize current knowledge on the link between AR repeat polymorphisms and cancer types, including mechanistic explanations and models comprising the relationship between condensate formation, polyQ1 length and transcriptional activity. Moreover, we outline the evolutionary paths of these recently evolved amino acid repeats across mammalian species, and discuss new research directions with potential breakthroughs and controversies in the literature.

Regions with two amino acids in protein sequences: a step forward from homorepeats into the low complexity landscape

Low complexity regions (LCRs) differ in amino acid composition from the background provided by the corresponding proteomes. The simplest LCRs are homorepeats (or polyX), regions composed of mostly-one amino acid type. Extensive research has been done to characterize homorepeats, and their taxonomic, functional and structural features depend on the amino acid type and sequence context. From them, the next step towards the study of LCRs are the regions composed of two types of amino acids, which we call polyXY. We classify polyXY in three categories based on the arrangement of the two amino acid types ‘X’ and ‘Y’: direpeats (e.g. ‘XYXYXY’), joined (e.g. ‘XXXYYY’) and shuffled (e.g. ‘XYYXXY’). We developed a script to search for polyXY, and located them in a comprehensive set of 20,340 reference proteomes. These results are available in a dedicated web server called XYs, in which the user can also submit their own protein datasets to detect polyXY. We studied the distribution of polyXY types by amino acid pair XY and category, and show that polyXY in Eukaryota are mainly located within intrinsically disordered regions. Our study provides a first step towards the characterization of polyXY as protein motifs.

SIN-3 functions through multi-protein interaction to regulate apoptosis, autophagy, and longevity in Caenorhabditis elegans

SIN3/HDAC is a multi-protein complex that acts as a regulatory unit and functions as a co-repressor/co-activator and a general transcription factor. SIN3 acts as a scaffold in the complex, binding directly to HDAC1/2 and other proteins and plays crucial roles in regulating apoptosis, differentiation, cell proliferation, development, and cell cycle. However, its exact mechanism of action remains elusive. Using the Caenorhabditis elegans (C. elegans) model, we can surpass the challenges posed by the functional redundancy of SIN3 isoforms. In this regard, we have previously demonstrated the role of SIN-3 in uncoupling autophagy and longevity in C. elegans. In order to understand the mechanism of action of SIN3 in these processes, we carried out a comparative analysis of the SIN3 protein interactome from model organisms of different phyla. We identified conserved, expanded, and contracted gene classes. The C. elegans SIN-3 interactome -revealed the presence of  well-known proteins, such as DAF-16, SIR-2.1, SGK-1, and AKT-1/2, involved in autophagy, apoptosis, and longevity. Overall, our analyses propose  potential mechanisms by which SIN3 participates in multiple biological processes and their conservation across species and identifies candidate genes for further experimental analysis.

A novel transcriptional activation mechanism of inulinase gene in Kluyveromyces marxianus involving a glycolysis regulator KmGcr1p with unique and functional Q-rich repeats

Kluyveromyces marxianus is the most suitable fungus for inulinase industrial production. However, the underlying transcriptional activation mechanism of the inulinase gene (INU1) is hitherto unclear. Here, we undertook genetic and biochemical analyses to elucidate that a glycolysis regulator KmGcr1p with unique Q-rich repeats is the key transcriptional activator of INU1. We determined that INU1 and glycolytic genes share similar transcriptional activation patterns, and that inulinase activity is induced by fermentable carbon sources including the hydrolysis products of inulin (fructose and glucose), which suggests a novel model of product feedback activation. Furthermore, all four CT-boxes in the INU1 promoter are important for KmGcr1p DNA binding in vitro, but the most downstream CT-box 1 primarily confers upstream activating sequence activity in vivo. More intriguingly, the use of artificial and natural GCR1 mutants suggests that the Q-rich repeats act as a functional module to maintain KmGcr1p transcriptional activity by contributing to its solubility and DNA binding affinity. Altogether, this study uncovers a novel transcriptional activation mechanism for the inulinase gene that is different from the previous understanding for filamentous fungi, but might have universal significance among inulinase-producing yeasts, thereby leading to a better understanding of the regulation mechanism of yeast inulinase genes.

Decoding the role of coiled-coil motifs in human prion-like proteins

Prions are self-propagating proteins that cause fatal neurodegenerative diseases in humans. However, increasing evidence suggests that eukaryotic cells exploit prion conformational conversion for functional purposes. A recent study delineated a group of twenty prion-like proteins in humans, characterized by the presence of low-complexity glutamine-rich sequences with overlapping coiled-coil (CCs) motifs. This is the case of Mediator complex subunit 15 (MED15), which is overexpressed in a wide range of human cancers. Biophysical studies demonstrated that the prion-like domain (PrLD) of MED15 forms homodimers in solution, sustained by CCs interactions. Furthermore, the same coiled-coil (CC) region plays a crucial role in the PrLD structural transition to a transmissible β-sheet amyloid state. In this review, we discuss the role of CCs motifs and their contribution to amyloid transitions in human prion-like domains (PrLDs), while providing a comprehensive overview of six predicted human prion-like proteins involved in transcription, gene expression, or DNA damage response and associated with human disease, whose PrLDs contain or overlap with CCs sequences. Finally, we try to rationalize how these molecular signatures might relate to both their function and involvement in disease.

Possible Role for Allelic Variation in Yeast MED15 in Ecological Adaptation

The propensity for Saccharomyces cerevisiae yeast to ferment sugars into ethanol and CO₂ has long been useful in the production of a wide range of food and drink. In the production of alcoholic beverages, the yeast strain selected for fermentation is crucial because not all strains are equally proficient in tolerating fermentation stresses. One potential mechanism by which domesticated yeast may have adapted to fermentation stresses is through changes in the expression of stress response genes. MED15 is a general transcriptional regulator and RNA Pol II Mediator complex subunit which modulates the expression of many metabolic and stress response genes. In this study, we explore the role of MED15 in alcoholic fermentation. In addition, we ask whether MED15 alleles from wine, sake or palm wine yeast improve fermentation activity and grape juice fermentation stress responses. And last, we investigate to what extent any differences in activity are due to allelic differences in the lengths of three polyglutamine tracts in MED15. We find that strains lacking MED15 are deficient in fermentation and fermentation stress responses and that MED15 alleles from alcoholic beverage yeast strains can improve both the fermentation capacity and the response to ethanol stresses when transplanted into a standard laboratory strain. Finally, we find that polyglutamine tract length in the Med15 protein is one determinant in the efficiency of the alcoholic fermentation process. These data lead to a working model in which polyglutamine tract length and other types of variability within transcriptional hubs like the Mediator subunit, Med15, may contribute to a reservoir of transcriptional profiles that may provide a fitness benefit in the face of environmental fluctuations.

Between Interactions and Aggregates: The PolyQ Balance

Polyglutamine (polyQ) regions are highly abundant consecutive runs of glutamine residues. They have been generally studied in relation to the so-called polyQ-associated diseases, characterized by protein aggregation caused by the expansion of the polyQ tract via a CAG-slippage mechanism. However, more than 4,800 human proteins contain a polyQ, and only nine of these regions are known to be associated with disease. Computational sequence studies and experimental structure determinations are completing a more interesting picture in which polyQ emerge as a motif for modulation of protein-protein interactions. But long polyQ regions may lead to an excess of interactions, and produce aggregates. Within this mechanistic perspective of polyQ function and malfunction, we discuss polyQ definition and properties such as variable codon usage, sequence and context structure imposition, functional relevance, evolutionary patterns in species-centered analyses, and open resources.

Stop codon readthrough alters the activity of a POU/Oct transcription factor during Drosophila development

Background

A number of cellular processes have evolved in metazoans that increase the proteome repertoire in relation to the genome, such as alternative splicing and translation recoding. Another such process, translational stop codon readthrough (SCR), generates C-terminally extended protein isoforms in many eukaryotes, including yeast, plants, insects, and humans. While comparative genome analyses have predicted the existence of programmed SCR in many species including humans, experimental proof of its functional consequences are scarce.

Results

We show that SCR of the Drosophila POU/Oct transcription factor Ventral veins lacking/Drifter (Vvl/Dfr) mRNA is prevalent in certain tissues in vivo, reaching a rate of 50% in the larval prothoracic gland. Phylogenetically, the C-terminal extension is conserved and harbors intrinsically disordered regions and amino acid stretches implied in transcriptional activation. Elimination of Vvl/Dfr translational readthrough by CRISPR/Cas9 mutagenesis changed the expression of a large number of downstream genes involved in processes such as chromatin regulation, neurogenesis, development, and immune response. As a proof-of-principle, we demonstrate that the C-terminal extension of Vvl/Dfr is necessary for correct timing of pupariation, by increasing the capacity to regulate its target genes. The extended Vvl/Dfr isoform acts in synergy with the transcription factor Molting defective (Mld) to increase the expression and biosynthesis of the steroid hormone ecdysone, thereby advancing pupariation. Consequently, late-stage larval development was prolonged and metamorphosis delayed in vvl/dfr readthrough mutants.

Conclusions

We demonstrate that translational recoding of a POU/Oct transcription factor takes place in a highly tissue-specific and temporally controlled manner. This dynamic and regulated recoding is necessary for normal expression of a large number of genes involved in many cellular and developmental processes. Loss of Vvl/Dfr translational readthrough negatively affects steroid hormone biosynthesis and delays larval development and progression into metamorphosis. Thus, this study demonstrates how SCR of a transcription factor can act as a developmental switch in a spatiotemporal manner, feeding into the timing of developmental transitions between different life-cycle stages.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01106-0.

PolyQ length co-evolution in neural proteins

Intermolecular co-evolution optimizes physiological performance in functionally related proteins, ultimately increasing molecular co-adaptation and evolutionary fitness. Polyglutamine (polyQ) repeats, which are over-represented in nervous system-related proteins, are increasingly recognized as length-dependent regulators of protein function and interactions, and their length variation contributes to intraspecific phenotypic variability and interspecific divergence. However, it is unclear whether polyQ repeat lengths evolve independently in each protein or rather co-evolve across functionally related protein pairs and networks, as in an integrated regulatory system. To address this issue, we investigated here the length evolution and co-evolution of polyQ repeats in clusters of functionally related and physically interacting neural proteins in Primates. We observed function-/disease-related polyQ repeat enrichment and evolutionary hypervariability in specific neural protein clusters, particularly in the neurocognitive and neuropsychiatric domains. Notably, these analyses detected extensive patterns of intermolecular polyQ length co-evolution in pairs and clusters of functionally related, physically interacting proteins. Moreover, they revealed both direct and inverse polyQ length co-variation in protein pairs, together with complex patterns of coordinated repeat variation in entire polyQ protein sets. These findings uncover a whole system of co-evolving polyQ repeats in neural proteins with direct implications for understanding polyQ-dependent phenotypic variability, neurocognitive evolution and neuropsychiatric disease pathogenesis.

LCD-Composer: an intuitive, composition-centric method enabling the identification and detailed functional mapping of low-complexity domains

Low complexity domains (LCDs) in proteins are regions predominantly composed of a small subset of the possible amino acids. LCDs are involved in a variety of normal and pathological processes across all domains of life. Existing methods define LCDs using information-theoretical complexity thresholds, sequence alignment with repetitive regions, or statistical overrepresentation of amino acids relative to whole-proteome frequencies. While these methods have proven valuable, they are all indirectly quantifying amino acid composition, which is the fundamental and biologically-relevant feature related to protein sequence complexity. Here, we present a new computational tool, LCD-Composer, that directly identifies LCDs based on amino acid composition and linear amino acid dispersion. Using LCD-Composer's default parameters, we identified simple LCDs across all organisms available through UniProt and provide the resulting data in an accessible form as a resource. Furthermore, we describe large-scale differences between organisms from different domains of life and explore organisms with extreme LCD content for different LCD classes. Finally, we illustrate the versatility and specificity achievable with LCD-Composer by identifying diverse classes of LCDs using both simple and multifaceted composition criteria. We demonstrate that the ability to dissect LCDs based on these multifaceted criteria enhances the functional mapping and classification of LCDs.

Defining the role of the polyasparagine repeat domain of the S. cerevisiae transcription factor Azf1p

Across eukaryotes, homopolymeric repeats of amino acids are enriched in regulatory proteins such as transcription factors and chromatin remodelers. These domains play important roles in signaling, binding, prion formation, and functional phase separation. Azf1p is a prion-forming yeast transcription factor that contains two homorepeat domains, a polyglutamine and a polyasparagine domain. In this work, we report a new phenotype for Azf1p and identify a large set of genes that are regulated by Azf1p during growth in glucose. We show that the polyasparagine (polyN) domain plays a subtle role in transcription but is dispensable for Azf1p localization and prion formation. Genes upregulated upon deletion of the polyN domain are enriched in functions related to carbon metabolism and storage. This domain may therefore be a useful target for engineering yeast strains for fermentation applications and small molecule production. We also report that both the polyasparagine and polyglutamine domains vary in length across strains of S. cerevisiae and propose a model for how this variation may impact protein function.

Regulation of mating and mating-type-specific genes in Zygosaccharomyces sp. yeast

Zygosaccharomyces sp. is an industrially important yeast for the production traditional fermented foods in Japan. At present, however, there is no easy method for mating Zygosaccharomyces sp. strains in the laboratory; furthermore, little is known about the expression of mating-type-specific genes in this yeast. Here, mating was observed when Zygosaccharomyces sp. was subjected to nitrogen-starvation conditions. The expression of mating-type-specific genes, Zygo STE6 and Zygo MFα1, was induced under nitrogen-starvation conditions, as confirmed by lacZ reporter assay. This expression was mating-type-specific: Zygo STE6 was expressed specifically for mating-type a, whereas and Zygo MFα1 was expressed specifically for mating-type α. Yeast strains Zygosaccharomyces rouxii DL2 and DA2, derived from type strain Z. rouxii CBS732, did not show mating even under nitrogen-starvation conditions. Gene sequencing revealed that the Zygo STE12 in Z. rouxii CBS732 has a frameshift mutation. Under nitrogen starvation, mating was observed in both DL2 and DA2 transformed with the wild-type Zygo STE12. The expression of Zygo STE6 in Z. rouxii DL2 transformed with wild-type Zygo STE12 under nitrogen-starvation conditions was confirmed by lacZ reporter assay. Collectively, these results revealed that, under nitrogen-starvation conditions, Zygosaccharomyces sp. can mate and mating-type-specific genes are expressed. Furthermore, Zygo Ste12 is essential for both mating and the expression of mating-type-specific genes in Zygosaccharomyces sp.

MED15 prion-like domain forms a coiled-coil responsible for its amyloid conversion and propagation

A disordered to β-sheet transition was thought to drive the functional switch of Q/N-rich prions, similar to pathogenic amyloids. However, recent evidence indicates a critical role for coiled-coil (CC) regions within yeast prion domains in amyloid formation. We show that many human prion-like domains (PrLDs) contain CC regions that overlap with polyQ tracts. Most of the proteins bearing these domains are transcriptional coactivators, including the Mediator complex subunit 15 (MED15) involved in bridging enhancers and promoters. We demonstrate that the human MED15-PrLD forms homodimers in solution sustained by CC interactions and that it is this CC fold that mediates the transition towards a β-sheet amyloid state, its chemical or genetic disruption abolishing aggregation. As in functional yeast prions, a GFP globular domain adjacent to MED15-PrLD retains its structural integrity in the amyloid state. Expression of MED15-PrLD in human cells promotes the formation of cytoplasmic and perinuclear inclusions, kidnapping endogenous full-length MED15 to these aggregates in a prion-like manner. The prion-like properties of MED15 are conserved, suggesting novel mechanisms for the function and malfunction of this transcription coactivator.

Coxsackievirus B3 targets TFEB to disrupt lysosomal function

Coxsackievirus B3 (CVB3) is a prevalent etiological agent for viral myocarditis and neurological disorders, particularly in infants and young children. Virus-encoded proteinases have emerged as cytopathic factors that contribute to disease pathogenesis in part through targeting the cellular recycling machinery of autophagy. Although it is appreciated that CVB3 can usurp cellular macroautophagy/autophagy for pro-viral functions, the precise mechanisms by which viral proteinases disrupt autophagy remain incompletely understood. Here we identified TFEB (transcription factor EB), a master regulator of autophagy and lysosome biogenesis, as a novel target of CVB3 proteinase 3 C. Time-course infections uncovered a significant loss of full-length TFEB and the emergence of a lower-molecular mass (~63 kDa) fragment. Cellular and in vitro cleavage assays revealed the involvement of viral proteinase 3 C in the proteolytic processing of TFEB, while site-directed mutagenesis identified the site of cleavage after glutamine 60. Assessment of TFEB transcriptional activity using a reporter construct discovered a loss of function of the cleavage fragment despite nuclear localization and retaining of the ability of DNA and protein binding. Furthermore, we showed that CVB3 infection was also able to trigger cleavage-independent nuclear translocation of TFEB that relied on the serine-threonine phosphatase PPP3/calcineurin. Finally, we demonstrated that both TFEB and TFEB [Δ60] serve roles in viral egress albeit through differing mechanisms. Collectively, this study reveals that CVB3 targets TFEB for proteolytic processing to disrupt host lysosomal function and enhance viral infection.

Abbreviations:ACTB: actin beta; CLEAR: coordinated lysosomal enhancement and regulation; CVB3: coxsackievirus B3; DAPI: 4′,6-diamidino-2-phenylindole; GFP: green fluorescent protein; LAMP1: lysosomal associated membrane protein 1; LTR: LysoTracker Red; PPP3/calcineurin: protein phosphatase 3; PPP3CA: protein phosphatase 3 catalytic subunit A; p-TFEB: phospho-Ser211 TFEB; si-CON: scramble control siRNA; TFEB: transcription factor EB; TFEB [Δ60]: TFEB cleavage fragment that lacks the first 60 amino acids; VP1: viral capsid protein 1

Identifying molecular features that are associated with biological function of intrinsically disordered protein regions

In previous work, we showed that intrinsically disordered regions (IDRs) of proteins contain sequence-distributed molecular features that are conserved over evolution, despite little sequence similarity that can be detected in alignments (Zarin et al., 2019). Here, we aim to use these molecular features to predict specific biological functions for individual IDRs and identify the molecular features within them that are associated with these functions. We find that the predictable functions are diverse. Examining the associated molecular features, we note some that are consistent with previous reports and identify others that were previously unknown. We experimentally confirm that elevated isoelectric point and hydrophobicity, features that are positively associated with mitochondrial localization, are necessary for mitochondrial targeting function. Remarkably, increasing isoelectric point in a synthetic IDR restores weak mitochondrial targeting. We believe feature analysis represents a new systematic approach to understand how biological functions of IDRs are specified by their protein sequences.

Moving beyond disease to function: Physiological roles for polyglutamine-rich sequences in cell decisions

Glutamine-rich tracts, also known as polyQ domains, have received a great deal of attention for their role in multiple neurodegenerative diseases, including Huntington's disease (HD), spinocerebellar ataxia (SCA), and others [22], [27]. Expansions in the normal polyQ tracts are thus commonly linked to disease, but polyQ domains themselves play multiple important functional roles in cells that are being increasingly appreciated. The biochemical nature of these domains allows them to adopt a number of different structures and form large assemblies that enable environmental responsiveness, localized signaling, and cellular memory. In many cases, these involve the formation of condensates that have varied material states. In this review, we highlight known and emerging functional roles for polyQ tracts in normal cell physiology.

Insights into gene regulation: From regulatory genomic elements to DNA-protein and protein-protein interactions

Transcription is orchestrated by non-coding regulatory elements embedded in chromatin, which exist within the larger context of chromosome topology. Here, we review recent insights into the functions of non-coding regulatory elements and their protein interactors during transcription control. A picture emerges in which the topological environment constraints enhancer-promoter interactions and specific enhancer-bound proteins with distinct promoter-compatibilities refine target promoter choice. Such compatibilities are encoded within the sequences of enhancers and promoters and realized by diverse transcription factors and cofactors with distinct biochemical activities. An emerging property of transcription factors and cofactors is the formation of nuclear microenvironments or membraneless compartments that can have properties of phase-separated liquids. These environments are able to selectively enrich certain proteins and small molecules over others. Further investigation into the interaction of transcriptional regulators with themselves and regulatory DNA elements will help reveal the complexities of gene regulation within the context of the nucleus.

β-catenin regulates FOXP2 transcriptional activity via multiple binding sites

The transcription factor forkhead box protein P2 (FOXP2) is a highly conserved key regulator of embryonal development. The molecular mechanisms of how FOXP2 regulates embryonal development, however, remain elusive. Using RNA sequencing, we identified the Wnt signaling pathway as key target of FOXP2‐dependent transcriptional regulation. Using cell‐based assays, we show that FOXP2 transcriptional activity is regulated by the Wnt coregulator β‐catenin and that β‐catenin contacts multiple regions within FOXP2. Using nuclear magnetic resonance spectroscopy, we uncovered the molecular details of these interactions. β‐catenin contacts a disordered FOXP2 region with α‐helical propensity via its folded armadillo domain, whereas the intrinsically disordered β‐catenin N terminus and C terminus bind to the conserved FOXP2 DNA‐binding domain. Using RNA sequencing, we confirmed that β‐catenin indeed regulates transcriptional activity of FOXP2 and that the FOXP2 α‐helical motif acts as a key regulatory element of FOXP2 transcriptional activity. Taken together, our findings provide first insight into novel regulatory interactions and help to understand the intricate mechanisms of FOXP2 function and (mis)‐regulation in embryonal development and human diseases.

Database

Expression data are available in the GEO database under the accession number GSE138938.

Evolution and expansion of the RUNX2 QA repeat corresponds with the emergence of vertebrate complexity

Runt-related transcription factor 2 (RUNX2) is critical for the development of the vertebrate bony skeleton. Unlike other RUNX family members, RUNX2 possesses a variable poly-glutamine, poly-alanine (QA) repeat domain. Natural variation within this repeat is able to alter the transactivation potential of RUNX2, acting as an evolutionary 'tuning knob' suggested to influence mammalian skull shape. However, the broader role of the RUNX2 QA repeat throughout vertebrate evolution is unknown. In this perspective, we examine the role of the RUNX2 QA repeat during skeletal development and discuss how its emergence and expansion may have facilitated the evolution of morphological novelty in vertebrates.

Crosstalk of promoter and terminator during RNA polymerase II transcription cycle

The transcription cycle of RNAPII is comprised of three consecutive steps; initiation, elongation and termination. It has been assumed that the initiation and termination steps occur in spatial isolation, essentially as independent events. A growing body of evidence, however, has challenged this dogma. First, factors involved in initiation and termination exhibit both a genetic and a physical interaction during transcription. Second, the initiation and termination factors have been found to occupy both ends of a transcribing gene. Third, physical interaction of initiation and termination factors occupying distal ends of a gene sometime results in the entire terminator region of a genes looping back and contact its cognate promoter, thereby forming a looped gene architecture during transcription. A logical interpretation of these findings is that the initiation and termination steps of transcription do not occur in isolation. There is extensive communication of factors occupying promoter and terminator ends of a gene during transcription cycle. This review entails a discussion of the promoter-terminator crosstalk and its implication in the context of transcription.

Amino acid homorepeats in proteins

Amino acid homorepeats, or homorepeats, are polypeptide segments found in proteins that contain stretches of identical amino acid residues. Although abnormal homorepeat expansions are linked to pathologies such as neurodegenerative diseases, homorepeats are prevalent in eukaryotic proteomes, suggesting that they are important for normal physiology. In this Review, we discuss recent advances in our understanding of the biological functions of homorepeats, which range from facilitating subcellular protein localization to mediating interactions between proteins across diverse cellular pathways. We explore how the functional diversity of homorepeat-containing proteins could be linked to the ability of homorepeats to adopt different structural conformations, an ability influenced by repeat composition, repeat length and the nature of flanking sequences. We conclude by highlighting how an understanding of homorepeats will help us better characterize and develop therapeutics against the human diseases to which they contribute.

Contractions of the C-Terminal Domain of Saccharomyces cerevisiae Rpb1p Are Mediated by Rad5p

The C-terminal domain (CTD) is an essential domain of the largest subunit of RNA polymerase II, Rpb1p, and is composed of 26 tandem repeats of a seven-amino acid sequence, YSPTSPS. Despite being an essential domain within an essential gene, we have previously demonstrated that the CTD coding region is genetically unstable. Furthermore, yeast with a truncated or mutated CTD sequence are capable of promoting spontaneous genetic expansion or contraction of this coding region to improve fitness. We investigated the mechanism by which the CTD contracts using a tet-off reporter system for RPB1 to monitor genetic instability within the CTD coding region. We report that contractions require the post-replication repair factor Rad5p but, unlike expansions, not the homologous recombination factors Rad51p and Rad52p. Sequence analysis of contraction events reveals that deleted regions are flanked by microhomologies. We also find that G-quadruplex forming sequences predicted by the QGRS Mapper are enriched on the noncoding strand of the CTD compared to the body of RPB1. Formation of G-quadruplexes in the CTD coding region could block the replication fork, necessitating post-replication repair. We propose that contractions of the CTD result when microhomologies misalign during Rad5p-dependent template switching via fork reversal.

MLLT10 in benign and malignant hematopoiesis

Non-random chromosomal translocations involving the putative transcription factor Mixed Lineage Leukemia Translocated to 10 (MLLT10, also known as AF10) are commonly observed in both acute myeloid and lymphoid leukemias and are indicative of a poor prognosis. Despite the well-described actions of oncogenic MLLT10 fusion proteins, the role of wild-type MLLT10 in hematopoiesis is not well characterized. The protein structure and several interacting partners have been described and provide indications as to the potential functions of MLLT10. This review examines these aspects of MLLT10, contextualizing its function in benign and malignant hematopoiesis.

Tuning up Transcription Factors for Therapy

The recent developments in the delivery and design of transcription factors put their therapeutic applications within reach, exemplified by cell replacement, cancer differentiation and T-cell based cancer therapies. The success of such applications depends on the efficacy and precision in the action of transcription factors. The biophysical and genetic characterization of the paradigmatic prokaryotic repressors, LacI and TetR and the designer transcription factors, transcription activator-like effector (TALE) and CRISPR-dCas9 revealed common principles behind their efficacy, which can aid the optimization of transcriptional activators and repressors. Further studies will be required to analyze the linkage between dissociation constants and enzymatic activity, the role of phase separation and squelching in activation and repression and the long-range interaction of transcription factors with epigenetic regulators in the context of the chromosomes. Understanding these mechanisms will help to tailor natural and synthetic transcription factors to the needs of specific applications.

Transcribed microsatellite allele lengths are often correlated with gene expression in natural sunflower populations

Microsatellites are common in genomes of most eukaryotic species. Due to their high mutability, an adaptive role for microsatellites has been considered. However, little is known concerning the contribution of microsatellites towards phenotypic variation. We used populations of the common sunflower (Helianthus annuus) at two latitudes to quantify the effect of microsatellite allele length on phenotype at the level of gene expression. We conducted a common garden experiment with seed collected from sunflower populations in Kansas and Oklahoma followed by an RNA-Seq experiment on 95 individuals. The effect of microsatellite allele length on gene expression was assessed across 3,325 microsatellites that could be consistently scored. Our study revealed 479 microsatellites at which allele length significantly correlates with gene expression (eSTRs). When irregular allele sizes not conforming to the motif length were removed, the number of eSTRs rose to 2,379. The percentage of variation in gene expression explained by eSTRs ranged from 1%-86% when controlling for population and allele-by-population interaction effects at the 479 eSTRs. Of these eSTRs, 70.4% are in untranslated regions (UTRs). A gene ontology (GO) analysis revealed that eSTRs are significantly enriched for GO terms associated with cis- and trans-regulatory processes. Our findings suggest that a substantial number of transcribed microsatellites can influence gene expression.

Compound Dynamics and Combinatorial Patterns of Amino Acid Repeats Encode a System of Evolutionary and Developmental Markers

Homopolymeric amino acid repeats (AARs) like polyalanine (polyA) and polyglutamine (polyQ) in some developmental proteins (DPs) regulate certain aspects of organismal morphology and behavior, suggesting an evolutionary role for AARs as developmental "tuning knobs." It is still unclear, however, whether these are occasional protein-specific phenomena or hints at the existence of a whole AAR-based regulatory system in DPs. Using novel approaches to trace their functional and evolutionary history, we find quantitative evidence supporting a generalized, combinatorial role of AARs in developmental processes with evolutionary implications. We observe nonrandom AAR distributions and combinations in HOX and other DPs, as well as in their interactomes, defining elements of a proteome-wide combinatorial functional code whereby different AARs and their combinations appear preferentially in proteins involved in the development of specific organs/systems. Such functional associations can be either static or display detectable evolutionary dynamics. These findings suggest that progressive changes in AAR occurrence/combination, by altering embryonic development, may have contributed to taxonomic divergence, leaving detectable traces in the evolutionary history of proteomes. Consistent with this hypothesis, we find that the evolutionary trajectories of the 20 AARs in eukaryotic proteomes are highly interrelated and their individual or compound dynamics can sharply mark taxonomic boundaries, or display clock-like trends, carrying overall a strong phylogenetic signal. These findings provide quantitative evidence and an interpretive framework outlining a combinatorial system of AARs whose compound dynamics mark at the same time DP functions and evolutionary transitions.

Unique near-complete deletion of GLI2 in a patient with combined pituitary hormone deficiency and post-axial polydactyly

INTRODUCTION

Combined pituitary hormone deficiency (CPHD) can cause a broad spectrum of health problems, ranging from short stature only, to convulsions or even death. In the majority of patients, the cause is unknown.

METHODS

The idex case had unexplained CPHD, pituitary anomalies on MRI and polydactyly. In the patients and her unaffected parents, we performed SNP array analysis and Whole Exome Sequencing, after candidate gene analysis turned out negative.

RESULTS

We found a unique de novo heterozygous 229.9 kb deletion in the index case on chr. 2q14.2. This deletion covered 12 out of the 13 coding exons of the GLI2 gene, a transcription factor involved in midline formation and previously associated with CPHD. As reported GLI2 deletions and mutations show a large phenotypic variability, we performed a genotype-phenotype analysis. This revealed that GLI2 missense mutations usually present with a 'ppp-only' phenotype (pituitary anomalies ± postaxial polydactyly without brain phenotype), whereas the 'ppp-plus' phenotype (with major brain malformations and/or intellectual disabilities) is more frequent in patients with larger deletions, and those with frameshift mutations/point mutations or splice variants resulting in a stop codon (p < .001).

CONCLUSION

The present case shows that a deletion of the GLI2 gene only (not affecting any of the adjacent genes) causes pituitary anomalies without brain phenotype. This suggests that brain phenotype only occurs when additional genes adjacent to GLI2 are deleted, or when mutations result in truncated GLI2 mRNA/protein. However, due to the lack of functional data for many GLI2 mutations and based on the available information regarding variable penetrance, phenotype-genotype correlations need to be made with caution.

Protein Abundance Biases the Amino Acid Composition of Disordered Regions to Minimize Non-functional Interactions

In eukaryotes, disordered regions cover up to 50% of proteomes and mediate fundamental cellular processes. In contrast to globular domains, where about half of the amino acids are buried in the protein interior, disordered regions show higher solvent accessibility, which makes them prone to engage in non-functional interactions. Such interactions are exacerbated by the law of mass action, prompting the question of how they are minimized in abundant proteins. We find that interaction propensity or "stickiness" of disordered regions negatively correlates with their cellular abundance, both in yeast and human. Strikingly, considering yeast proteins where a large fraction of the sequence is disordered, the correlation between stickiness and abundance reaches R=-0.55. Beyond this global amino-acid composition bias, we identify three rules by which amino-acid composition of disordered regions adjusts with high abundance. First, lysines are preferred over arginines, consistent with the latter amino acid being stickier than the former. Second, compensatory effects exist, whereby a sticky region can be tolerated if it is compensated by a distal non-sticky region. Third, such compensation requires a lower average stickiness at the same abundance when compared to a scenario where stickiness is homogeneous throughout the sequence. We validate these rules experimentally, employing them as different strategies to rescue an otherwise sticky protein fragment from aggregation. Our results highlight that non-functional interactions represent a significant constraint in cellular systems and reveal simple rules by which protein sequences adapt to that constraint. Data from this work are deposited in Figshare, at https://doi.org/10.6084/m9.figshare.8068937.v3.

No evidence for an association between Clock gene allelic variation and migration timing in a long-distance migratory shorebird (Limosa lapponica baueri)

The gene Clock is a key part of the Core Circadian Oscillator, and the length of the polyglutamine (poly-Q) repeat sequence in Clock (ClkpolyQcds) has been proposed to be associated with the timing of annual cycle events in birds. We tested whether variation in ClkpolyQcds corresponds to variation in migration timing in the bar-tailed godwit (Limosa lapponica baueri), a species in which individuals show strong annual consistency in their migration timing despite the New Zealand population migrating across a 5-week period. We describe allelic variation of the ClkpolyQcds in 135 godwits over-wintering in New Zealand (N.Z.) and investigate whether polymorphism in this region is associated with northward migration timing (chronophenotype) from N.Z. or (for 32 birds tracked by geolocator) after the primary stopover in Asia. Six Clock alleles were detected (Q₇‒Q₁₂) and there was substantial variation between individuals (heterozygosity of 0.79). There was no association between ClkpolyQcds polymorphism and migration timing from N.Z. The length of the shorter Clock allele was related to migration timing from Asia, though this relationship arose largely from just a few northern-breeding birds with longer alleles. Other studies show no consistent associations between ClkpolyQcds and migration timing in birds, although Clock may be associated with breeding latitude in some species (as an adaptation to photoperiodic regime). Apparent relationships with migration timing could reflect latitude-related variation in migration timing, rather than Clock directly affecting migration timing. On current evidence, ClkpolyQcds is not a strong candidate for driving migration timing in migratory birds generally.

CAGI SickKids challenges: Assessment of phenotype and variant predictions derived from clinical and genomic data of children with undiagnosed diseases

Whole-genome sequencing (WGS) holds great potential as a diagnostic test. However, the majority of patients currently undergoing WGS lack a molecular diagnosis, largely due to the vast number of undiscovered disease genes and our inability to assess the pathogenicity of most genomic variants. The CAGI SickKids challenges attempted to address this knowledge gap by assessing state-of-the-art methods for clinical phenotype prediction from genomes. CAGI4 and CAGI5 participants were provided with WGS data and clinical descriptions of 25 and 24 undiagnosed patients from the SickKids Genome Clinic Project, respectively. Predictors were asked to identify primary and secondary causal variants. In addition, for CAGI5, groups had to match each genome to one of three disorder categories (neurologic, ophthalmologic, and connective), and separately to each patient. The performance of matching genomes to categories was no better than random but two groups performed significantly better than chance in matching genomes to patients. Two of the ten variants proposed by two groups in CAGI4 were deemed to be diagnostic, and several proposed pathogenic variants in CAGI5 are good candidates for phenotype expansion. We discuss implications for improving in silico assessment of genomic variants and identifying new disease genes.

Copper(II)-Binding Induces a Unique Polyproline Type II Helical Structure within the Ion-Binding Segment in the Intrinsically Disordered F-Domain of Ecdysteroid Receptor from Aedes aegypti

Reproduction of the dominant vector of Zika and dengue diseases, Aedes aegypti mosquito, is controlled by an active heterodimer complex composed of the 20-hydroxyecdysone receptor (EcR) and ultraspiracle protein. Although A. aegypti EcR shares the structural and functional organization with other nuclear receptors, its C-terminus has an additional long F domain (AaFEcR). Recently, we showed that the full length AaFEcR is intrinsically disordered with the ability to specifically bind divalent metal ions. Here, we describe the details of the exhaustive structural and thermodynamic properties of Zn²⁺- and Cu²⁺-complexes with the AaFEcR domain, based on peptide models of its two putative metal binding sites (Ac-HGPHPHPHG-NH₂ and Ac-QQLTPNQQQHQQQHSQLQQVHANGS-NH₂). Unexpectedly, only in the presence of increasing concentrations of Cu²⁺ ions, the Ac-HGPHPHPHG-NH₂ peptide gained a metal ion-induced poly-l-proline type II helical structure, which is unique for members of the family of nuclear receptors.

Mutations in Cytosine-5 tRNA Methyltransferases Impact Mobile Element Expression and Genome Stability at Specific DNA Repeats

The maintenance of eukaryotic genome stability is ensured by the interplay of transcriptional as well as post-transcriptional mechanisms that control recombination of repeat regions and the expression and mobility of transposable elements. We report here that mutations in two (cytosine-5) RNA methyltransferases, Dnmt2 and NSun2, impact the accumulation of mobile element-derived sequences and DNA repeat integrity in Drosophila. Loss of Dnmt2 function caused moderate effects under standard conditions, while heat shock exacerbated these effects. In contrast, NSun2 function affected mobile element expression and genome integrity in a heat shock-independent fashion. Reduced tRNA stability in both RCMT mutants indicated that tRNA-dependent processes affected mobile element expression and DNA repeat stability. Importantly, further experiments indicated that complex formation with RNA could also contribute to the impact of RCMT function on gene expression control. These results thus uncover a link between tRNA modification enzymes, the expression of repeat DNA, and genomic integrity.

TcaB, an insecticidal protein from Photorhabdus akhurstii causes cytotoxicity in the greater wax moth, Galleria mellonella

Photorhabdus akhurstii can produce a variety of proteins that aid this bacterium and its mutualistic nematode vector, Heterorhabditis indica to kill the insect host. Herein, we characterized (by heterologously expressing in E. coli) an open reading frame (1713 bp) of the toxin complex protein, TcaB from P. akhurstii strains IARI-SGHR2 and IARI-SGMS1 and assessed its toxic effect on G. mellonella larvae. The intra-hemocoel injection of purified TcaB (molecular weight-63 kDa) caused fourth instar larval bodies to blacken and die with LD₅₀ values of 67.25 (IARI-SGHR2) and 52.08 (IARI-SGMS1) ng per larva at 12 h. Additionally, oral administration of the toxin caused larval mortality with LD₅₀ values of 709.55 (IARI-SGHR2) and 598.44 (IARI-SGMS1) ng per g diet per larva at 7 days post feeding. Injection of purified TcaB caused loss of viability of fourth instar G. mellonella hemocytes at 6 h post incubation; cells displayed morphological changes typical of apoptosis, including cell shrinkage, membrane blebbing, nuclear condensation and disintegration. Injection of TcaB also elevated the phenoloxidase activity in insect hemolymph which triggers an extensive immune response that potentially leads to larval death. Similar to other bacterial toxins TcaB possesses potent biological activity which may enable it to be used as an efficient agent for pest management.

Function determinants of TET proteins: the arrangements of sequence motifs with specific codes

The ten-eleven translocation (TET) proteins play a crucial role in promoting locus-specific reversal of DNA methylation, a type of chromatin modification. Considerable evidences have demonstrated that the sequence motifs with specific codes are important to determine the functions of domains and active sites. Here, we surveyed major studies and reviews regarding the multiple functions of the TET proteins and established the patterns of the motif arrangements that determine the functions of TET proteins. First, we summarized the functional sequence basis of TET proteins and identified the new functional motifs based on the phylogenetic relationship. Next, we described the sequence characteristics of the functional motifs in detail and provided an overview of the relationship between the sequence motifs and the functions of TET proteins, including known functions and potential functions. Finally, we highlighted that sequence motifs with diverse post-translational modifications perform unique functions, and different selection pressures lead to different arrangements of sequence motifs, resulting in different paralogs and isoforms.

Side chain to main chain hydrogen bonds stabilize a polyglutamine helix in a transcription factor

Polyglutamine (polyQ) tracts are regions of low sequence complexity frequently found in transcription factors. Tract length often correlates with transcriptional activity and expansion beyond specific thresholds in certain human proteins is the cause of polyQ disorders. To study the structural basis of the association between tract length, transcriptional activity and disease, we addressed how the conformation of the polyQ tract of the androgen receptor, associated with spinobulbar muscular atrophy (SBMA), depends on its length. Here we report that this sequence folds into a helical structure stabilized by unconventional hydrogen bonds between glutamine side chains and main chain carbonyl groups, and that its helicity directly correlates with tract length. These unusual hydrogen bonds are bifurcate with the conventional hydrogen bonds stabilizing α-helices. Our findings suggest a plausible rationale for the association between polyQ tract length and androgen receptor transcriptional activity and have implications for establishing the mechanistic basis of SBMA.

Polyglutamine (polyQ) tracts are low-complexity regions and their expansion is linked to certain neurodegenerative diseases. Here the authors combine experimental and computational approaches to find that the length of the androgen receptor polyQ tract correlates with its helicity and show that the polyQ helical structure is stabilized by hydrogen bonds between the Gln side chains and main chain carbonyl groups.