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

Reassessing the association: Evaluation of a polyalanine deletion variant of RUNX2 in non-syndromic sagittal and metopic craniosynostosis

The RUNT-related transcription factor RUNX2 plays a critical role in osteoblast differentiation, and alterations to gene dosage cause distinct craniofacial anomalies. Uniquely amongst the RUNT-related family, vertebrate RUNX2 encodes a polyglutamine/polyalanine repeat (Gln23-Glu-Ala17 in humans), with the length of the polyalanine component completely conserved in great apes. Surprisingly, a frequent 6-amino acid deletion polymorphism, p.(Ala84_Ala89)del, occurs in humans (termed 11A allele), and a previous association study (Cuellar et al. Bone 137:115395;2020) reported that the 11A variant was significantly more frequent in non-syndromic sagittal craniosynostosis (nsSag; allele frequency [AF] = 0.156; 95% confidence interval [CI] 0.126-0.189) compared to non-syndromic metopic craniosynostosis (nsMet; AF = 0.068; 95% CI 0.045-0.098). However, the gnomAD v.2.1.1 control population used by Cuellar et al. did not display Hardy-Weinberg equilibrium, hampering interpretation. To re-examine this association, we genotyped the RUNX2 11A polymorphism in 225 individuals with sporadic nsSag as parent-child trios and 164 singletons with sporadic nsMet, restricting our analysis to individuals of European ancestry. We compared observed allele frequencies to the non-transmitted alleles in the parent-child trios, and to the genome sequencing data from gnomAD v.4, which display Hardy-Weinberg equilibrium. Observed AFs (and 95% CI) were 0.076 (0.053-0.104) in nsSag and 0.082 (0.055-0.118) in nsMet, compared with 0.062 (0.042-0.089) in non-transmitted parental alleles and 0.065 (0.063-0.067) in gnomAD v.4.0.0 non-Finnish European control genomes. In summary, we observed a non-significant excess, compared to gnomAD data, of 11A alleles in both nsSag (relative risk 1.18, 95% CI 0.83-1.67) and nsMet (relative risk 1.29, 95% CI 0.87-1.92), but we did not replicate the much higher excess of RUNX2 11A alleles in nsSag previously reported (p = 0.0001).

Evolutionary mechanisms modulating the mammalian skull development

Mammals possess impressive craniofacial variation that mirrors their adaptation to diverse ecological niches, feeding behaviour, physiology and overall lifestyle. The spectrum of craniofacial geometries is established mainly during embryonic development. The formation of the head represents a sequence of events regulated on genomic, molecular, cellular and tissue level, with each step taking place under tight spatio-temporal control. Even minor variations in timing, position or concentration of the molecular drivers and the resulting events can affect the final shape, size and position of the skeletal elements and the geometry of the head. Our knowledge of craniofacial development increased substantially in the last decades, mainly due to research using conventional vertebrate model organisms. However, how developmental differences in head formation arise specifically within mammals remains largely unexplored. This review highlights three evolutionary mechanisms acknowledged to modify ontogenesis: heterochrony, heterotopy and heterometry. We present recent research that links changes in developmental timing, spatial organization or gene expression levels to the acquisition of species-specific skull morphologies. We highlight how these evolutionary modifications occur on the level of the genes, molecules and cellular processes, and alter conserved developmental programmes to generate a broad spectrum of skull shapes characteristic of the class Mammalia. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

Accumulation of NCOA1 dependent on HERC3 deficiency transactivates matrix metallopeptidases and promotes extracellular matrix degradation in intervertebral disc degeneration

BACKGROUND

Matrix metallopeptidases (MMPs) are critical matrix-degrading molecules and they are frequently overexpressed in degenerative discs. This study aimed to investigate the mechanism for MMP upregulation.

METHODS

Immunoblot and RT-qPCR were used for detecting protein and gene expression levels. 4-month-old and 24-month-old C57BL/6 mice were used for evaluating intervertebral disc degeneration (IDD). An ubiquitination assay was used to determine protein modification. Immunoprecipitation and mass spectrometry were used for identifying protein complex members.

RESULTS

We identified the elevation of 14 MMPs among 23 members in aged mice with IDD. Eleven of these 14 MMP gene promoters contained a Runx2 (runt-related transcription factor 2) binding site. Biochemical analyses revealed that Runx2 recruited a histone acetyltransferase p300 and a coactivator NCOA1 (nuclear receptor coactivator 1) to assemble a complex, transactivating MMP expression. The deficiency of an E3 ligase called HERC3 (HECT and RLD domain containing E3 ubiquitin-protein ligase 3) resulted in the accumulation of NCOA1 in the inflammatory microenvironment. High throughput screening of small molecules that specifically target the NCOA1-p300 interaction identified a compound SMTNP-191, which showed an inhibitory effect on suppressing MMP expression and attenuating the IDD process in aged mice.

CONCLUSION

Our data support a model in which deficiency of HERC3 fails to ubiquitinate NCOA1, leading to the assembly of NCOA1-p300-Runx2 and causing the transactivation of MMPs. These findings offer new insight into inflammation-mediated MMP accumulation and also provide a new therapeutic strategy to retard the IDD process.

Mutation and selection processes regulating short tandem repeats give rise to genetic and phenotypic diversity across species

Short tandem repeats (STRs) are units of 1-6 bp that repeat in a tandem fashion in DNA. Along with single nucleotide polymorphisms and large structural variations, they are among the major genomic variants underlying genetic, and likely phenotypic, divergence. STRs experience mutation rates that are orders of magnitude higher than other well-studied genotypic variants. Frequent copy number changes result in a wide range of alleles, and provide unique opportunities for modulating complex phenotypes through variation in repeat length. While classical studies have identified key roles of individual STR loci, the advent of improved sequencing technology, high-quality genome assemblies for diverse species, and bioinformatics methods for genome-wide STR analysis now enable more systematic study of STR variation across wide evolutionary ranges. In this review, we explore mutation and selection processes that affect STR copy number evolution, and how these processes give rise to varying STR patterns both within and across species. Finally, we review recent examples of functional and adaptive changes linked to STRs.

Lineage-specific protein repeat expansions and contractions reveal malleable regions of immune genes

Functional diversification, a higher evolutionary rate, and intense positive selection help a limited number of immune genes interact with many pathogens. Repeats in protein-coding regions are a well-known source of functional diversification, adaptive variation, and evolutionary novelty in a short time. Repeats play a crucial role in biochemical functions like functional diversification of transcription regulation, protein kinases, cell adhesion, signaling pathways, morphogenesis, DNA repair, recombination, and RNA processing. Repeat length variation can change the associated protein's interaction, efficacy, and overall protein network. Repeats have an intrinsic unstable nature and can potentially evolve rapidly and expedite the acquisition of complex phenotypic traits and functions. Because of their ability to generate rapid, adaptive variations over short evolutionary distances, repeats are considered "tuning knobs." Repeat length variation in specific genes, like RUNX2 and ALX4, is associated with morphological and physiological changes across vertebrates. Here we study repeat length variation as a potent source of species-specific immune diversification across several clades of tetrapods. Moreover, we provide a clade-wise comprehensive list of immune genes with repeat types for future studies of morphological/evolutionary changes within species groups. We observe significant repeat length variation of FASLG and C1QC in Rodentia and Primates' contrasting species groups, respectively.

Marsupials and Multi-Omics: Establishing New Comparative Models of Neural Crest Patterning and Craniofacial Development

Studies across vertebrates have revealed significant insights into the processes that drive craniofacial morphogenesis, yet we still know little about how distinct facial morphologies are patterned during development. Studies largely point to evolution in GRNs of cranial progenitor cell types such as neural crest cells, as the major driver underlying adaptive cranial shapes. However, this hypothesis requires further validation, particularly within suitable models amenable to manipulation. By utilizing comparative models between related species, we can begin to disentangle complex developmental systems and identify the origin of species-specific patterning. Mammals present excellent evolutionary examples to scrutinize how these differences arise, as sister clades of eutherians and marsupials possess suitable divergence times, conserved cranial anatomies, modular evolutionary patterns, and distinct developmental heterochrony in their NCC behaviours and craniofacial patterning. In this review, I lend perspectives into the current state of mammalian craniofacial biology and discuss the importance of establishing a new marsupial model, the fat-tailed dunnart, for comparative research. Through detailed comparisons with the mouse, we can begin to decipher mammalian conserved, and species-specific processes and their contribution to craniofacial patterning and shape disparity. Recent advances in single-cell multi-omics allow high-resolution investigations into the cellular and molecular basis of key developmental processes. As such, I discuss how comparative evolutionary application of these tools can provide detailed insights into complex cellular behaviours and expression dynamics underlying adaptive craniofacial evolution. Though in its infancy, the field of "comparative evo-devo-omics" presents unparalleled opportunities to precisely uncover how phenotypic differences arise during development.

Shaping modern human skull through epigenetic, transcriptional and post-transcriptional regulation of the RUNX2 master bone gene

RUNX2 encodes the master bone transcription factor driving skeletal development in vertebrates, and playing a specific role in craniofacial and skull morphogenesis. The anatomically modern human (AMH) features sequence changes in the RUNX2 locus compared with archaic hominins' species. We aimed to understand how these changes may have contributed to human skull globularization occurred in recent evolution. We compared in silico AMH and archaic hominins' genomes, and used mesenchymal stromal cells isolated from skull sutures of craniosynostosis patients for in vitro functional assays. We detected 459 and 470 nucleotide changes in noncoding regions of the AMH RUNX2 locus, compared with the Neandertal and Denisovan genomes, respectively. Three nucleotide changes in the proximal promoter were predicted to alter the binding of the zinc finger protein Znf263 and long-distance interactions with other cis-regulatory regions. By surface plasmon resonance, we selected nucleotide substitutions in the 3'UTRs able to affect miRNA binding affinity. Specifically, miR-3150a-3p and miR-6785-5p expression inversely correlated with RUNX2 expression during in vitro osteogenic differentiation. The expression of two long non-coding RNAs, AL096865.1 and RUNX2-AS1, within the same locus, was modulated during in vitro osteogenic differentiation and correlated with the expression of specific RUNX2 isoforms. Our data suggest that RUNX2 may have undergone adaptive phenotypic evolution caused by epigenetic and post-transcriptional regulatory mechanisms, which may explain the delayed suture fusion leading to the present-day globular skull shape.