Missense Variants Reveal Functional Insights Into the Human ARID Family of Gene Regulators.

Missense variants are alterations to protein coding sequences that result in amino acid substitutions. They can be deleterious if the amino acid is required for maintaining structure or/and function, but are likely to be tolerated at other sites. Consequently, missense variation within a healthy population can mirror the effects of negative selection on protein structure and function, such that functional sites on proteins are often depleted of missense variants. Advances in high-throughput sequencing have dramatically increased the sample size of available human variation data, allowing for population-wide analysis of selective pressures. In this study, we developed a convenient set of tools, called 1D-to-3D, for visualizing the positions of missense variants on protein sequences and structures. We used these tools to characterize human homologues of the ARID family of gene regulators. ARID family members are implicated in multiple cancer types, developmental disorders, and immunological diseases but current understanding of their mechanistic roles is incomplete. Combined with phylogenetic and structural analyses, our approach allowed us to characterise sites important for protein-protein interactions, histone modification recognition, and DNA binding by the ARID proteins. We find that comparing missense depletion patterns among paralogs can reveal sub-functionalization at the level of domains. We propose that visualizing missense variants and their depletion on structures can serve as a valuable tool for complementing evolutionary and experimental findings.

Classification of MSH6 Variants of Uncertain Significance Using Functional Assays.

Lynch syndrome (LS) is one of the most common hereditary cancer predisposition syndromes worldwide. Individuals with LS have a high risk of developing colorectal or endometrial cancer, as well as several other cancers. LS is caused by autosomal dominant pathogenic variants in one of the DNA mismatch repair (MMR) genes MLH1, MSH2, PMS2 or MSH6, and typically include truncating variants, such as frameshift, nonsense or splicing variants. However, a significant number of missense, intronic, or silent variants, or small in-frame insertions/deletions, are detected during genetic screening of the MMR genes. The clinical effects of these variants are often more difficult to predict, and a large fraction of these variants are classified as variants of uncertain significance (VUS). It is pivotal for the clinical management of LS patients to have a clear genetic diagnosis, since patients benefit widely from screening, preventive and personal therapeutic measures. Moreover, in families where a pathogenic variant is identified, testing can be offered to family members, where non-carriers can be spared frequent surveillance, while carriers can be included in cancer surveillance programs. It is therefore important to reclassify VUSs, and, in this regard, functional assays can provide insight into the effect of a variant on the protein or mRNA level. Here, we briefly describe the disorders that are related to MMR deficiency, as well as the structure and function of MSH6. Moreover, we review the functional assays that are used to examine VUS identified in MSH6 and discuss the results obtained in relation to the ACMG/AMP PS3/BS3 criterion. We also provide a compiled list of the MSH6 variants examined by these assays. Finally, we provide a future perspective on high-throughput functional analyses with specific emphasis on the MMR genes.

A phylogenetically conserved APETALA2/ETHYLENE RESPONSE FACTOR, ERF12, regulates Arabidopsis floral development.

KEY MESSAGE: Arabidopsis ETHYLENE RESPONSE FACTOR12 (ERF12), the rice MULTIFLORET SPIKELET1 orthologue pleiotropically affects meristem identity, floral phyllotaxy and organ initiation and is conserved among angiosperms. Reproductive development necessitates the coordinated regulation of meristem identity and maturation and lateral organ initiation via positive and negative regulators and network integrators. We have identified ETHYLENE RESPONSE FACTOR12 (ERF12) as the Arabidopsis orthologue of MULTIFLORET SPIKELET1 (MFS1) in rice. Loss of ERF12 function pleiotropically affects reproductive development, including defective floral phyllotaxy and increased floral organ merosity, especially supernumerary sepals, at incomplete penetrance in the first-formed flowers. Wildtype floral organ number in early formed flowers is labile, demonstrating that floral meristem maturation involves the stabilisation of positional information for organogenesis, as well as appropriate identity. A subset of erf12 phenotypes partly defines a narrow developmental time window, suggesting that ERF12 functions heterochronically to fine-tune stochastic variation in wild type floral number and similar to MFS1, promotes meristem identity. ERF12 expression encircles incipient floral primordia in the inflorescence meristem periphery and is strong throughout the floral meristem and intersepal regions. ERF12 is a putative transcriptional repressor and genetically opposes the function of its relatives DORNRÖSCHEN, DORNRÖSCHEN-LIKE and PUCHI and converges with the APETALA2 pathway. Phylogenetic analysis suggests that ERF12 is conserved among all eudicots and appeared in angiosperm evolution concomitant with the generation of floral diversity.

Mechanistic insights into the evolution of DUF26-containing proteins in land plants.

Large protein families are a prominent feature of plant genomes and their size variation is a key element for adaptation. However, gene and genome duplications pose difficulties for functional characterization and translational research. Here we infer the evolutionary history of the DOMAIN OF UNKNOWN FUNCTION (DUF) 26-containing proteins. The DUF26 emerged in secreted proteins. Domain duplications and rearrangements led to the appearance of CYSTEINE-RICH RECEPTOR-LIKE PROTEIN KINASES (CRKs) and PLASMODESMATA-LOCALIZED PROTEINS (PDLPs). The DUF26 is land plant-specific but structural analyses of PDLP ectodomains revealed strong similarity to fungal lectins and thus may constitute a group of plant carbohydrate-binding proteins. CRKs expanded through tandem duplications and preferential retention of duplicates following whole genome duplications, whereas PDLPs evolved according to the dosage balance hypothesis. We propose that new gene families mainly expand through small-scale duplications, while fractionation and genetic drift after whole genome multiplications drive families towards dosage balance.

Automated in situ chromatin profiling efficiently resolves cell types and gene regulatory programs.

BACKGROUND: Our understanding of eukaryotic gene regulation is limited by the complexity of protein-DNA interactions that comprise the chromatin landscape and by inefficient methods for characterizing these interactions. We recently introduced CUT&RUN, an antibody-targeted nuclease cleavage method that profiles DNA-binding proteins, histones and chromatin-modifying proteins in situ with exceptional sensitivity and resolution. RESULTS: Here, we describe an automated CUT&RUN platform and apply it to characterize the chromatin landscapes of human cells. We find that automated CUT&RUN profiles of histone modifications crisply demarcate active and repressed chromatin regions, and we develop a continuous metric to identify cell-type-specific promoter and enhancer activities. We test the ability of automated CUT&RUN to profile frozen tumor samples and find that our method readily distinguishes two pediatric glioma xenografts by their subtype-specific gene expression programs. CONCLUSIONS: The easy, cost-effective workflow makes automated CUT&RUN an attractive tool for high-throughput characterization of cell types and patient samples.

Caught with One's Zinc Fingers in the Genome Integrity Cookie Jar.

Zinc finger (ZnF) domains are present in at least 5% of human proteins. First characterized as binding to DNA, ZnFs display extraordinary binding plasticity and can bind to RNA, lipids, proteins, and protein post-translational modifications (PTMs). The diverse binding properties of ZnFs have made their functional characterization challenging. While once confined to large and poorly characterized protein families, proteomic, cellular, and molecular studies have begun to shed light on their involvement as protectors of the genome. We focus here on the emergent roles of ZnF domain-containing proteins in promoting genome integrity, including their involvement in telomere maintenance and DNA repair. These findings have highlighted the need for further characterization of ZnF proteins, which can reveal the functions of this large gene class in normal cell function and human diseases, including those involving genome instability such as aging and cancer.

Cloning, characterization and subcellular localization of Nuclear LIM interactor interacting factor gene from Leishmania donovani.

LIM domains are zinc-binding motifs that mediate protein-protein interactions and are found in a wide variety of cytoplasmic and nuclear proteins. The nuclear LIM domain family members have a number of different functions including transcription factors, gene regulation, cell fate determination, organization of the cytoskeleton and tumour formation exerting their function through various LIM domain interacting protein partners/cofactors. Nuclear LIM domain interacting proteins/factors have not been reported in any protozoan parasites including Leishmania. Here, we report for the first time cloning, characterization and subcellular localization of nuclear LIM interactor-interacting factor (NLI) like protein from Leishmania donovani, the causative agent of Indian Kala-azar. Primary sequence analysis of LdNLI revealed presence of characteristic features of nuclear LIM interactor-interacting factor. However, leishmanial NLI represents a distinct kinetoplastid group, clustered in a separate branch of the phylogenic tree. The sub-cellular distribution of LdNLI revealed the discreet localization in nucleus and kinetoplast only, suggesting that the gene may have a role in parasite gene expression.

Embryo development. A cysteine-clamp gene drives embryo polarity in the midge Chironomus.

In the fruit fly Drosophila, head formation is driven by a single gene, bicoid, which generates head-to-tail polarity of the main embryonic axis. Bicoid deficiency results in embryos with tail-to-tail polarity and no head. However, most insects lack bicoid, and the molecular mechanism for establishing head-to-tail polarity is poorly understood. We have identified a gene that establishes head-to-tail polarity of the mosquito-like midge, Chironomus riparius. This gene, named panish, encodes a cysteine-clamp DNA binding domain and operates through a different mechanism than bicoid. This finding, combined with the observation that the phylogenetic distributions of panish and bicoid are limited to specific families of flies, reveals frequent evolutionary changes of body axis determinants and a remarkable opportunity to study gene regulatory network evolution.

CpG island-mediated global gene regulatory modes in mouse embryonic stem cells.

Both transcriptional and epigenetic regulations are fundamental for the control of eukaryotic gene expression. Here we perform a compendium analysis of >200 large sequencing data sets to elucidate the regulatory logic of global gene expression programs in mouse embryonic stem (ES) cells. We define four major classes of DNA-binding proteins (Core, PRC, MYC and CTCF) based on their target co-occupancy, and discover reciprocal regulation between the MYC and PRC classes for the activity of nearly all genes under the control of the CpG island (CGI)-containing promoters. This CGI-dependent regulatory mode explains the functional segregation between CGI-containing and CGI-less genes during early development. By defining active enhancers based on the co-occupancy of the Core class, we further demonstrate their additive roles in CGI-containing gene expression and cell type-specific roles in CGI-less gene expression. Altogether, our analyses provide novel insights into previously unknown CGI-dependent global gene regulatory modes.

Genome-wide screens for sensitivity to ionizing radiation identify the fission yeast nonhomologous end joining factor Xrc4.

Nonhomologous end joining (NHEJ) is the main means for repairing DNA double-strand breaks (DSBs) in human cells. Molecular understanding of NHEJ has benefited from analyses in the budding yeast Saccharomyces cerevisiae and the fission yeast Schizosaccharomyces pombe. In human cells, the DNA ligation reaction of the classical NHEJ pathway is carried out by a protein complex composed of DNA ligase IV (LigIV) and XRCC4. In S. cerevisiae, this reaction is catalyzed by a homologous complex composed of Dnl4 and Lif1. Intriguingly, no homolog of XRCC4 has been found in S. pombe, raising the possibility that such a factor may not always be required for classical NHEJ. Here, through screening the ionizing radiation (IR) sensitivity phenotype of a genome-wide fission yeast deletion collection in both the vegetative growth state and the spore state, we identify Xrc4, a highly divergent homolog of human XRCC4. Like other fission yeast NHEJ factors, Xrc4 is critically important for IR resistance of spores, in which no homologous recombination templates are available. Using both extrachromosomal and chromosomal DSB repair assays, we show that Xrc4 is essential for classical NHEJ. Exogenously expressed Xrc4 colocalizes with the LigIV homolog Lig4 at the chromatin region of the nucleus in a mutually dependent manner. Furthermore, like their human counterparts, Xrc4 and Lig4 interact with each other and this interaction requires the inter-BRCT linker and the second BRCT domain of Lig4. Our discovery of Xrc4 suggests that an XRCC4 family protein is universally required for classical NHEJ in eukaryotes.

Genome-wide identification, classification and expression analysis of the heat shock transcription factor family in Chinese cabbage.

The Hsf gene family, one of the most important transcription factor families, plays crucial roles in regulating heat resistance. However, a systematic and comprehensive analysis of this gene family has not been reported in Chinese cabbage. Therefore, systematic analysis of the Hsf gene family in Chinese cabbage has profound significance. In this study, 35 BrHsf genes were identified from Chinese cabbage, which could be classified into three groups according to their structural characteristics and phylogenetic comparisons with Arabidopsis and rice. Thirty-three BrHsf genes mapped on chromosomes were further assigned to three subgenomes and eight ancestral karyotypes. Distribution mapping showed that BrHsf genes were non-randomly localized on chromosomes. Chinese cabbage and Arabidopsis shared 22 orthologous gene pairs. The expansion of BrHsf genes mainly resulted from genome triplication. Comparative analysis showed that the most Hsf genes were in Chinese cabbage among the five species analyzed. Interestingly, the number of Hsf genes of heat-resistant plants (Theobroma cacao and Musa acuminata) was fewer than that in Chinese cabbage. The expression patterns of BrHsf genes were different in six tissues, based on RNA-seq. Quantitative real-time-PCR analysis showed that the expression level of BrHsf genes varied under various abiotic stresses. In conclusion, this comprehensive analysis of BrHsf genes will provide rich resources, aiding the determination of Hsfs functions in plant heat resistance. Furthermore, the comparative genomics analysis deepened our understanding of Hsf genes' evolution accompanied by the polyploidy event of Chinese cabbage.

Composition and interrelationships of a large Neotropical freshwater fish group, the subfamily Cheirodontinae (Characiformes: Characidae): a case study based on mitochondrial and nuclear DNA sequences.

Characidae is the most species-rich family of freshwater fishes in the order Characiformes, with more than 1000 valid species that correspond to approximately 55% of the order. Few hypotheses about the composition and internal relationships within this family are available and most fail to reach an agreement. Among Characidae, Cheirodontinae is an emblematic group that includes 18 genera (1 fossil) and approximately 60 described species distributed throughout the Neotropical region. The taxonomic and systematic history of Cheirodontinae is complex, and only two hypotheses about the internal relationships in this subfamily have been reported to date. In the present study, we test the composition and relationships of fishes assigned to Cheirodontinae based on a broad taxonomic sample that also includes some characid incertae sedis taxa that were previously considered to be part of Cheirodontinae. We present phylogenetic analyses of a large molecular dataset of mitochondrial and nuclear DNA sequences. Our results reject the monophyly of Cheirodontinae as previously conceived, as well as the tribes Cheirodontini and Compsurini, and the genera Cheirodon, Compsura, Leptagoniates, Macropsobrycon, Odontostilbe, and Serrapinnus. On the basis of these results we propose: (1) the exclusion of Amazonspinther and Spintherobolus from the subfamily Cheirodontinae since they are the sister-group of all remaining Characidae; (2) the removal of Macropsobrycon xinguensis of the genus Macropsobrycon; (3) the removal of Leptagoniates pi of the genus Leptagoniates; (4) the inclusion of Leptagoniates pi in the subfamily Cheirodontinae; (5) the removal of Cheirodon stenodon of the genus Cheirodon and its inclusion in the subfamily Cheirodontinae under a new genus name; (6) the need to revise the polyphyletic genera Compsura, Odontostilbe, and Serrapinnus; and (7) the division of Cheirodontinae in three newly defined monophyletic tribes: Cheirodontini, Compsurini, and Pseudocheirodontini. Our results suggest that our knowledge about the largest Neotropical fish family, Characidae, still is incipient.

[Classification of amino acids based on a comparative analysis of contacts in DNA-protein complexes and specific DNA-protein interactions].

The classification of amino acid residues based on the events of contact formation between distinct amino acid and selected nucleotides was constructed. Thus, the most integral properties, that characterize interactions in organization of DNA-protein complexes, were used. We applied the Voronoi-Delaunay tessellation to draw statistics of contacts and area of contacts for the set included 1937 DNA-protein complexes. Similarities of amino acid residues have been searched for based on the comparison of corresponded rows and matrixes of contacts and areas of contacts. Nine measures of distance were used for estimation of rows similarity degree. The procedure of clustering amino acids in groups included three hierarchical and two nonhierarchical methods. A total tree was built using nine techniques of estimating distance with three hierarchical clustering methods. It was shown that clustering centers in the main groups are always constant while other relationships between objects vary. Clustering of binary associations was found for the most amino acids. Major classes of up to six amino acids correspond to the certain local structures of the polypeptide chain in the context of amino acid composition. These data should be taken into account when designing DNA-protein ligands.

c-Ski in health and disease.

c-Ski is an evolutionary conserved protein that is involved in diverse cellular processes such as proliferation, differentiation, transformation, and tumor progression. A large range of cellular partners of c-Ski, including transcription factors, chromatin-remodeling molecules, tumor suppressors, and nuclear hormone receptors, has been identified. Moreover, numerous mechanisms have been described by which c-Ski regulates essential signaling pathways, e.g., the TGFß pathway. In this review, we summarize the diverse roles attributed to c-Ski during normal development and in cancer progression and discuss future strategies to unravel further the complex nature of c-Ski actions in a context-dependent manner.

The rhomboid protease family: a decade of progress on function and mechanism.

Rhomboid proteases are the largest family of enzymes that hydrolyze peptide bonds within the cell membrane. Although discovered to be serine proteases only a decade ago, rhomboid proteases are already considered to be the best understood intramembrane proteases. The presence of rhomboid proteins in all domains of life emphasizes their importance but makes their evolutionary history difficult to chart with confidence. Phylogenetics nevertheless offers three guiding principles for interpreting rhomboid function. The near ubiquity of rhomboid proteases across evolution suggests broad, organizational roles that are not directly essential for cell survival. Functions have been deciphered in only about a dozen organisms and fall into four general categories: initiating cell signaling in animals, facilitating bacterial quorum sensing, regulating mitochondrial homeostasis, and dismantling adhesion complexes of parasitic protozoa. Although in no organism has the full complement of rhomboid function yet been elucidated, links to devastating human disease are emerging rapidly, including to Parkinson's disease, type II diabetes, cancer, and bacterial and malaria infection. Rhomboid proteases are unlike most proteolytic enzymes, because they are membrane-immersed; understanding how the membrane immersion affects their function remains a key challenge.

Characterisation of the heat shock factor of the human thermodimorphic pathogen Paracoccidioides lutzii.

Thermodimorphic fungi include most causative agents of systemic mycoses, but the molecular mechanisms that underlie their defining trait, i.e. the ability to shift between mould and yeast on temperature change alone, remain poorly understood. We hypothesised that the heat shock factor (Hsf), a protein that evolved to sense thermal stimuli quickly, might play a role in this process in addition to the known regulator Drk1 and the Ryp proteins. To test this hypothesis, we characterised the Hsf from the thermodimorph Paracoccidioides lutzii (formerly Paracoccidioides brasiliensis isolate 01). We show in the present work that PlHsf possesses regulatory domains that are exclusive of the Eurotiomycetidae family, suggesting evolutionary specialisation; that it can successfully rescue the otherwise lethal loss of the native protein of Saccharomyces cerevisiae; and that its DNA-binding domain is able to recognise regulatory elements from the promoters of both Drk1 and Ryp1. An in silico screening of all 1 kb sequences upstream of P. lutzii ORFs revealed that 7% of them possess a heat shock element. This is the first description of a heat shock factor in a thermodimorphic fungus.

ATM activation in the presence of oxidative stress.

The Ataxia-Telangiectasia mutated (ATM) kinase is regarded as the major regulator of the cellular response to DNA double strand breaks (DSBs). In response to DSBs, ATM dimers dissociate into active monomers in a process promoted by the Mre11-Rad50-Nbs1 (MRN) complex. ATM can also be activated by oxidative stress directly in the form of exposure to H2O2. The active ATM in this case is a disulfide-crosslinked dimer containing 2 or more disulfide bonds. Mutation of a critical cysteine residue in the FATC domain involved in disulfide bond formation specifically blocks ATM activation by oxidative stress. Here we show that ATM activation by DSBs is inhibited in the presence of H2O2 because oxidation blocks the ability of MRN to bind to DNA. However, ATM activation via direct oxidation by H2O2 complements the loss of MRN/DSB-dependent activation and contributes significantly to the overall level of ATM activity in the presence of both DSBs and oxidative stress.

Structural basis for the specialization of Nur, a nickel-specific Fur homolog, in metal sensing and DNA recognition.

Nur, a member of the Fur family, is a nickel-responsive transcription factor that controls nickel homeostasis and anti-oxidative response in Streptomyces coelicolor. Here we report the 2.4-A resolution crystal structure of Nur. It contains a unique nickel-specific metal site in addition to a nonspecific common metal site. The identification of the 6-5-6 motif of the Nur recognition box and a Nur/DNA complex model reveals that Nur mainly interacts with terminal bases of the palindrome on complex formation. This contrasts with more distributed contacts between Fur and the n-1-n type of the Fur-binding motif. The disparity between Nur and Fur in the conformation of the S1-S2 sheet in the DNA-binding domain can explain their different DNA-recognition patterns. Furthermore, the fact that the specificity of Nur in metal sensing and DNA recognition is conferred by the specific metal site suggests that its introduction drives the evolution of Nur orthologs in the Fur family.