YBX1-interacting small RNAs and RUNX2 can be blocked in primary bone cancer using CADD522.

Primary bone cancer (PBC) comprises several subtypes each underpinned by distinctive genetic drivers. This driver diversity produces novel morphological features and clinical behaviour that serendipitously makes PBC an excellent metastasis model. Here, we report that some transfer RNA-derived small RNAs termed tRNA fragments (tRFs) perform as a constitutive tumour suppressor mechanism by blunting a potential pro-metastatic protein-RNA interaction. This mechanism is reduced in PBC progression with a gradual loss of tRNAGlyTCC cleavage into 5' end tRF-GlyTCC when comparing low-grade, intermediate-grade and high-grade patient tumours. We detected recurrent activation of miR-140 leading to upregulated RUNX2 expression in high-grade patient tumours. Both tRF-GlyTCC and RUNX2 share a sequence motif in their 3' ends that matches the YBX1 recognition site known to stabilise pro-metastatic mRNAs. Investigating some aspects of this interaction network, gain- and loss-of-function experiments using small RNA mimics and antisense LNAs, respectively, showed that ectopic tRF-GlyTCC reduced RUNX2 expression and dispersed 3D micromass architecture in vitro. iCLIP sequencing revealed YBX1 physical binding to the 3' UTR of RUNX2. The interaction between YBX1, tRF-GlyTCC and RUNX2 led to the development of the RUNX2 inhibitor CADD522 as a PBC treatment. CADD522 assessment in vitro revealed significant effects on PBC cell behaviour. In xenograft mouse models, CADD522 as a single agent without surgery significantly reduced tumour volume, increased overall and metastasis-free survival and reduced cancer-induced bone disease. Our results provide insight into PBC molecular abnormalities that have led to the identification of new targets and a new therapeutic.

3' Untranslated Region Structural Elements in CYP24A1 Are Associated With Infantile Hypercalcemia Type 1.

Loss-of-function mutations in the CYP24A1 protein-coding region causing reduced 25 hydroxyvitamin D (25OHD) and 1,25 dihydroxyvitamin D (1,25(OH)2 D) catabolism have been observed in some cases of infantile hypercalcemia type 1 (HCINF1), which can manifest as nephrocalcinosis, hypercalcemia and adult-onset hypercalciuria, and renal stone formation. Some cases present with apparent CYP24A1 phenotypes but do not exhibit pathogenic mutations. Here, we assessed the molecular mechanisms driving apparent HCINF1 where there was a lack of CYP24A1 mutation. We obtained blood samples from 47 patients with either a single abnormality of no obvious cause or a combination of hypercalcemia, hypercalciuria, and nephrolithiasis as part of our metabolic and stone clinics. We used liquid chromatography tandem mass spectrometry (LC-MS/MS) to determine serum vitamin D metabolites and direct sequencing to confirm CYP24A1 genotype. Six patients presented with profiles characteristic of altered CYP24A1 function but lacked protein-coding mutations in CYP24A1. Analysis upstream and downstream of the coding sequence showed single nucleotide variants (SNVs) in the CYP24A1 3' untranslated region (UTR). Bioinformatics approaches revealed that these 3' UTR abnormalities did not result in microRNA silencing but altered the CYP24A1 messenger RNA (mRNA) secondary structure, which negatively impacted translation. Our experiments showed that mRNA misfolding driven by these 3' UTR sequence-dependent structural elements was associated with normal 25OHD but abnormal 1,25(OH)2 D catabolism. Using CRISPR-Cas9 gene editing, we developed an in vitro mutant model for future CYP24A1 studies. Our results form a basis for future studies investigating structure-function relationships and novel CYP24A1 mutations producing a semifunctional protein. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Mechanistic insights into non-coding Y RNA processing.

Y RNAs (84-112 nt) are non-coding RNAs transcribed by RNA polymerase III and are characterized by a distinctive secondary structure. Human Y RNAs interact with the autoimmune proteins SSB and RO60 that together form a ribonucleoprotein (RNP) complex termed RoRNP and Y RNAs also perform regulatory roles in DNA and RNA replication and stability, which has major implications for diseases including cancer. During cellular stress and apoptosis, Y RNAs are cleaved into 3' and 5' end fragments termed Y RNA-derived small RNAs (ysRNAs). Although some ysRNA functions in stress, apoptosis and cancer have been reported, their fundamental biogenesis has not been described. Here we report that 3' end RNY5 cleavage is structure dependent. In high throughput mutagenesis experiments, cleavage occurred between the 2nd and 3rd nt above a double stranded stem comprising high GC content. We demonstrate that an internal loop above stem S3 is critical for producing 3' end ysRNAs (31 nt) with mutants resulting in longer or no ysRNAs. We show a UGGGU sequence motif at position 22 of RNY5 is critical for producing 5' end ysRNAs (22-25 nt). We show that intact RO60 is critical for ysRNA biogenesis. We conclude that ribonuclease L (RNASEL) contributes to Y RNA cleavage in mouse embryonic fibroblasts but is not the only endoribonuclease important in human cells.

miR-7b-3p Exerts a Dual Role After Spinal Cord Injury, by Supporting Plasticity and Neuroprotection at Cortical Level.

Spinal cord injury (SCI) affects 6 million people worldwide with no available treatment. Despite research advances, the inherent poor regeneration potential of the central nervous system remains a major hurdle. Small RNAs (sRNAs) 19-33 nucleotides in length are a set of non-coding RNA molecules that regulate gene expression and have emerged as key players in regulating cellular events occurring after SCI. Here we profiled a class of sRNA known as microRNAs (miRNAs) following SCI in the cortex where the cell bodies of corticospinal motor neurons are located. We identified miR-7b-3p as a candidate target given its significant upregulation after SCI in vivo and we screened by miRWalk PTM the genes predicted to be targets of miR-7b-3p (among which we identified Wipf2, a gene regulating neurite extension). Moreover, 16 genes, involved in neural regeneration and potential miR-7b-3p targets, were found to be downregulated in the cortex following SCI. We also analysed miR-7b-3p function during cortical neuron development in vitro: we observed that the overexpression of miR-7b-3p was important (1) to maintain neurons in a more immature and, likely, plastic neuronal developmental phase and (2) to contrast the apoptotic pathway; however, in normal conditions it did not affect the Wipf2 expression. On the contrary, the overexpression of miR-7b-3p upon in vitro oxidative stress condition (mimicking the SCI environment) significantly reduced the expression level of Wipf2, as observed in vivo, confirming it as a direct miR-7b-3p target. Overall, these data suggest a dual role of miR-7b-3p: (i) the induction of a more plastic neuronal condition/phase, possibly at the expense of the axon growth, (ii) the neuroprotective role exerted through the inhibition of the apoptotic cascade. Increasing the miR-7b-3p levels in case of SCI could reactivate in adult neurons silenced developmental programmes, supporting at the same time the survival of the axotomised neurons.

The P2RX7B splice variant modulates osteosarcoma cell behaviour and metastatic properties.

Background: Osteosarcoma (OS) is the most common type of primary bone cancer affecting children and adolescents. OS has a high propensity to spread meaning the disease is often incurable and fatal. There have been no improvements in survival rates for decades. This highlights an urgent need for the development of novel therapeutic strategies. Here, we report in vitro and in vivo data that demonstrates the role of purinergic signalling, specifically, the B isoform of the purinergic receptor P2RX7 (P2RX7B), in OS progression and metastasis. Methods: TE85 and MNNG-HOS OS cells were transfected with P2RX7B. These cell lines were then characterised and assessed for proliferation, cell adhesion, migration and invasion in vitro. We used these cells to perform both paratibial and tail vein injected mouse studies where the primary tumour, bone and lungs were analysed. We used RNA-seq to identify responsive pathways relating to P2RX7B. Results: Our data shows that P2RX7B expression confers a survival advantage in TE85 + P2RX7B and MNNG-HOS + P2RX7B human OS cell lines in vitro that is minimised following treatment with A740003, a specific P2RX7 antagonist. P2RX7B expression reduced cell adhesion and P2RX7B activation promoted invasion and migration in vitro, demonstrating a metastatic phenotype. Using an in vivo OS xenograft model, MNNG-HOS + P2RX7B tumours exhibited cancer-associated ectopic bone formation that was abrogated with A740003 treatment. A pro-metastatic phenotype was further demonstrated in vivo as expression of P2RX7B in primary tumour cells increased the propensity of tumour cells to metastasise to the lungs. RNA-seq identified a novel gene axis, FN1/LOX/PDGFB/IGFBP3/BMP4, downregulated in response to A740003 treatment. Conclusion: Our data illustrates a role for P2RX7B in OS tumour growth, progression and metastasis. We show that P2RX7B is a future therapeutic target in human OS.

Targeting the MAPK7/MMP9 axis for metastasis in primary bone cancer.

Metastasis is the leading cause of cancer-related death. This multistage process involves contribution from both tumour cells and the tumour stroma to release metastatic cells into the circulation. Circulating tumour cells (CTCs) survive circulatory cytotoxicity, extravasate and colonise secondary sites effecting metastatic outcome. Reprogramming the transcriptomic landscape is a metastatic hallmark, but detecting underlying master regulators that drive pathological gene expression is a key challenge, especially in childhood cancer. Here we used whole tumour plus single-cell RNA-sequencing in primary bone cancer and CTCs to perform weighted gene co-expression network analysis to systematically detect coordinated changes in metastatic transcript expression. This approach with comparisons applied to data collected from cell line models, clinical samples and xenograft mouse models revealed mitogen-activated protein kinase 7/matrix metallopeptidase 9 (MAPK7/MMP9) signalling as a driver for primary bone cancer metastasis. RNA interference knockdown of MAPK7 reduces proliferation, colony formation, migration, tumour growth, macrophage residency/polarisation and lung metastasis. Parallel to these observations were reduction of activated interleukins IL1B, IL6, IL8 plus mesenchymal markers VIM and VEGF in response to MAPK7 loss. Our results implicate a newly discovered, multidimensional MAPK7/MMP9 signalling hub in primary bone cancer metastasis that is clinically actionable.

Artificially induced phased siRNAs promote virus resistance in transgenic plants.

We previously developed transgenic tobacco plants that were resistant to two geminiviruses. We generated resistance using RNAi constructs that produced trans-acting siRNA (tasiRNA) like secondary siRNAs known as phased siRNA (phasiRNA) that targeted several regions of Tomato Leaf Curl New Delhi Virus (ToLCNDV) and Tomato Leaf Curl Gujarat Virus (ToLCGV) transcripts encoding the RNA silencing suppressor proteins AC2 and AC4. Here, we performed degradome analysis to determine the precise cleavage sites of RNA-RNA interaction between phasiRNA and viral transcripts. We then applied our RNAi technology in tomato, which is the natural host for ToLCNDV and ToLCGV. The relative ease of developing and using phasiRNA constructs represents a significant technical advance in imparting virus resistance in crops and/or important model systems.

Molecular insights into an ancient form of Paget's disease of bone.

Paget's disease of bone (PDB) is a chronic skeletal disorder that can affect one or several bones in individuals older than 55 y of age. PDB-like changes have been reported in archaeological remains as old as Roman, although accurate diagnosis and natural history of the disease is lacking. Six skeletons from a collection of 130 excavated at Norton Priory in the North West of England, which dates to medieval times, show atypical and extensive pathological changes resembling contemporary PDB affecting as many as 75% of individual skeletons. Disease prevalence in the remaining collection is high, at least 16% of adults, with age at death estimations as low as 35 y. Despite these atypical features, paleoproteomic analysis identified sequestosome 1 (SQSTM1) or p62, a protein central to the pathological milieu of PDB, as one of the few noncollagenous human sequences preserved in skeletal samples. Targeted proteomic analysis detected >60% of the ancient p62 primary sequence, with Western blotting indicating p62 abnormalities, including in dentition. Direct sequencing of ancient DNA excluded contemporary PDB-associated SQSTM1 mutations. Our observations indicate that the ancient p62 protein is likely modified within its C-terminal ubiquitin-associated domain. Ancient miRNAs were remarkably preserved in an osteosarcoma from a skeleton with extensive disease, with miR-16 expression consistent with that reported in contemporary PDB-associated bone tumors. Our work displays the use of proteomics to inform diagnosis of ancient diseases such as atypical PDB, which has unusual features presumably potentiated by yet-unidentified environmental or genetic factors.

microRNAs associated with early neural crest development in Xenopus laevis.

BACKGROUND: The neural crest (NC) is a class of transitory stem cell-like cells unique to vertebrate embryos. NC cells arise within the dorsal neural tube where they undergo an epithelial to mesenchymal transition in order to migrate and differentiate throughout the developing embryo. The derivative cell types give rise to multiple tissues, including the craniofacial skeleton, peripheral nervous system and skin pigment cells. Several well-studied gene regulatory networks underpin NC development, which when disrupted can lead to various neurocristopathies such as craniofrontonasal dysplasia, DiGeorge syndrome and some forms of cancer. Small RNAs, such as microRNAs (miRNAs) are non-coding RNA molecules important in post-transcriptional gene silencing and critical for cellular regulation of gene expression. RESULTS: To uncover novel small RNAs in NC development we used high definition adapters and next generation sequencing of libraries derived from ectodermal explants of Xenopus laevis embryos induced to form neural and NC tissue. Ectodermal and blastula animal pole (blastula) stage tissues were also sequenced. We show that miR-427 is highly abundant in all four tissue types though in an isoform specific manner and we define a set of 11 miRNAs that are enriched in the NC. In addition, we show miR-301a and miR-338 are highly expressed in both the NC and blastula suggesting a role for these miRNAs in maintaining the stem cell-like phenotype of NC cells. CONCLUSION: We have characterised the miRNAs expressed in Xenopus embryonic explants treated to form ectoderm, neural or NC tissue. This has identified novel tissue specific miRNAs and highlighted differential expression of miR-427 isoforms.

MicroRNA expression in a phosphaturic mesenchymal tumour.

Phosphaturic mesenchymal tumours are a heterogeneous set of bone and soft tissue neoplasms that can cause a number of paraneoplastic syndromes such as tumour induced osteomalacia. The term phosphaturic comes from the common finding that these tumours secrete high levels of fibroblast growth factor 23 which causes renal phosphate wasting leading to hypophosphatemia. Phosphaturic mesenchymal tumours are rare and diagnosis is difficult. A very active 68 year old male presented with bone pain and muscle weakness. He was hypophosphataemic and total alkaline phosphatase was markedly elevated. The patient was placed on vitamin D supplementation but his condition progressed. In the fifth year of presentation the patient required the use of a wheelchair and described "explosive" bone pain on physical contact. Serum 1,25 dihydroxyvitamin D was low and serum fibroblast growth factor 23 was significantly elevated, raising suspicion of a phosphaturic mesenchymal tumour. A lesion was detected in his left femoral head and the patient underwent a total hip replacement. The patient displayed a rapid improvement to his condition and during a three year follow up period he returned to an active lifestyle. As molecular testing may help provide a robust diagnosis and is particularly useful in rare diseases we took a next generation sequencing approach to identify a differential expression of small RNAs in the resected tumour. Small RNAs are non-coding RNA molecules that play a key role in regulation of gene expression and can be used as specific biomarkers. We found an upregulation of miR-197. We also found a downregulation of miR-20b, miR-144 and miR-335 which is a small RNA profile typical of osteosarcoma. MiR-21, the most frequently upregulated microRNA in cancer, was downregulated. We conclude that the specific small RNA profile is typical of osteosarcoma except for the downregulation of oncogenic miR-21. Transcriptional plasticity of miR-197, which is computationally predicted to target fibroblast growth factor 23 messenger RNA, may be upregulated in a cellular effort to correct the ectopic expression of the protein.

Correction: Chromosomal-Level Assembly of the Asian Seabass Genome Using Long Sequence Reads and Multi-layered Scaffolding.

[This corrects the article DOI: 10.1371/journal.pgen.1005954.].

Transfer RNA-derived small RNAs in the cancer transcriptome.

The cellular lifetime includes stages such as differentiation, proliferation, division, senescence and apoptosis. These stages are driven by a strictly ordered process of transcription dynamics. Molecular disruption to RNA polymerase assembly, chromatin remodelling and transcription factor binding through to RNA editing, splicing, post-transcriptional regulation and ribosome scanning can result in significant costs arising from genome instability. Cancer development is one example of when such disruption takes place. RNA silencing is a term used to describe the effects of post-transcriptional gene silencing mediated by a diverse set of small RNA molecules. Small RNAs are crucial for regulating gene expression and microguarding genome integrity. RNA silencing studies predominantly focus on small RNAs such as microRNAs, short-interfering RNAs and piwi-interacting RNAs. We describe an emerging renewal of interest in a 'larger' small RNA, the transfer RNA (tRNA). Precisely generated tRNA-derived small RNAs, named tRNA halves (tiRNAs) and tRNA fragments (tRFs), have been reported to be abundant with dysregulation associated with cancer. Transfection of tiRNAs inhibits protein translation by displacing eukaryotic initiation factors from messenger RNA (mRNA) and inaugurating stress granule formation. Knockdown of an overexpressed tRF inhibits cancer cell proliferation. Recovery of lacking tRFs prevents cancer metastasis. The dual oncogenic and tumour-suppressive role is typical of functional small RNAs. We review recent reports on tiRNA and tRF discovery and biogenesis, identification and analysis from next-generation sequencing data and a mechanistic animal study to demonstrate their physiological role in cancer biology. We propose tRNA-derived small RNA-mediated RNA silencing is an innate defence mechanism to prevent oncogenic translation. We expect that cancer cells are percipient to their ablated control of transcription and attempt to prevent loss of genome control through RNA silencing.

Chromosomal-Level Assembly of the Asian Seabass Genome Using Long Sequence Reads and Multi-layered Scaffolding.

We report here the ~670 Mb genome assembly of the Asian seabass (Lates calcarifer), a tropical marine teleost. We used long-read sequencing augmented by transcriptomics, optical and genetic mapping along with shared synteny from closely related fish species to derive a chromosome-level assembly with a contig N50 size over 1 Mb and scaffold N50 size over 25 Mb that span ~90% of the genome. The population structure of L. calcarifer species complex was analyzed by re-sequencing 61 individuals representing various regions across the species' native range. SNP analyses identified high levels of genetic diversity and confirmed earlier indications of a population stratification comprising three clades with signs of admixture apparent in the South-East Asian population. The quality of the Asian seabass genome assembly far exceeds that of any other fish species, and will serve as a new standard for fish genomics.

Role of miR-140 in embryonic bone development and cancer.

Bone is increasingly viewed as an endocrine organ with key biological functions. The skeleton produces hormones and cytokines, such as FGF23 and osteocalcin, which regulate an extensive list of homoeostatic functions. Some of these functions include glucose metabolism, male fertility, blood cell production and calcium/phosphate metabolism. Many of the genes regulating these functions are specific to bone cells. Some of these genes can be wrongly expressed by other malfunctioning cells, driving the generation of disease. The miRNAs are a class of non-coding RNA molecules that are powerful regulators of gene expression by suppressing and fine-tuning target mRNAs. Expression of one such miRNA, miR-140, is ubiquitous in chondrocyte cells during embryonic bone development. Activity in cells found in the adult breast, colon and lung tissue can silence genes required for tumour suppression. The realization that the same miRNA can be both normal and detrimental, depending on the cell, tissue and time point, provides a captivating twist to the study of whole-organism functional genomics. With the recent interest in miRNAs in bone biology and RNA-based therapeutics on the horizon, we present a review on the role of miR-140 in the molecular events that govern bone formation in the embryo. Cellular pathways involving miR-140 may be reactivated or inhibited when treating skeletal injury or disorder in adulthood. These pathways may also provide a novel model system when studying cancer biology of other cells and tissues.

MicroRNAs influence reproductive responses by females to male sex peptide in Drosophila melanogaster.

Across taxa, female behavior and physiology change significantly following the receipt of ejaculate molecules during mating. For example, receipt of sex peptide (SP) in female Drosophila melanogaster significantly alters female receptivity, egg production, lifespan, hormone levels, immunity, sleep, and feeding patterns. These changes are underpinned by distinct tissue- and time-specific changes in diverse sets of mRNAs. However, little is yet known about the regulation of these gene expression changes, and hence the potential role of microRNAs (miRNAs), in female postmating responses. A preliminary screen of genomic responses in females to receipt of SP suggested that there were changes in the expression of several miRNAs. Here we tested directly whether females lacking four of the candidate miRNAs highlighted (miR-279, miR-317, miR-278, and miR-184) showed altered fecundity, receptivity, and lifespan responses to receipt of SP, when mated once or continually to SP null or control males. The results showed that miRNA-lacking females mated to SP null males exhibited altered receptivity, but not reproductive output, in comparison to controls. However, these effects interacted significantly with the genetic background of the miRNA-lacking females. No significant survival effects were observed in miRNA-lacking females housed continually with SP null or control males. However, continual exposure to control males that transferred SP resulted in significantly higher variation in miRNA-lacking female lifespan than did continual exposure to SP null males. The results provide the first insight into the effects and importance of miRNAs in regulating postmating responses in females.

Age-dependent female responses to a male ejaculate signal alter demographic opportunities for selection.

A central tenet of evolutionary explanations for ageing is that the strength of selection wanes with age. However, data on age-specific expression and benefits of sexually selected traits are lacking-particularly for traits subject to sexual conflict. We addressed this by using as a model the responses of Drosophila melanogaster females of different ages to receipt of sex peptide (SP), a seminal fluid protein transferred with sperm during mating. SP can mediate sexual conflict, benefitting males while causing fitness costs in females. Virgin and mated females of all ages showed significantly reduced receptivity in response to SP. However, only young virgin females also showed increased egg laying; hence, there was a narrow demographic window of maximal responses to SP. Males gained significant 'per mating' fitness benefits only when mating with young females. The pattern completely reversed in matings with older females, where SP transfer was costly. The overall benefits of SP transfer (hence opportunity for selection) therefore reversed with female age. The data reveal a new example of demographic variation in the strength of selection, with convergence and conflicts of interest between males and ageing females occurring over different facets of responses to a sexually antagonistic trait.