Dynamics of microRNA expression during mouse prenatal development

Transcriptomic profiling of gastrointestinal tracts in dairy cattle during lactation reveals molecular adaptations for milk synthesis

During lactation, dairy cattle's digestive tract requires significant adaptations to meet the increased nutrient demands for milk production. As we attempt to improve milk-related traits through selective pressure, it is crucial to understand the biological functions of the epithelia of the rumen, small intestine, and colonic tissues in response to changes in physiological state driven by changes in nutrient demands for milk synthesis. In this study, we obtained a total of 108 transcriptome profiles from three tissues (epithelia of the colon, duodenum, and rumen) of five Holstein cows, spanning eight time points from the early, mid, late lactation periods to the dry period. On average 97.06% of reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We analyzed 27,607 gene expression patterns at multiple periods, enabling direct comparisons within and among tissues during different lactation stages, including early and peak lactation. We identified 1645, 813, and 2187 stage-specific genes in the colon, duodenum, and rumen, respectively, which were enriched for common or specific biological functions among different tissues. Time series analysis categorized the expressed genes within each tissue into four clusters. Furthermore, when the three tissues were analyzed collectively, 36 clusters of similarly expressed genes were identified. By integrating other comprehensive approaches such as gene co-expression analyses, functional enrichment, and cell type deconvolution, we gained profound insights into cattle lactation, revealing tissue-specific characteristics of the gastrointestinal tract and shedding light on the intricate molecular adaptations involved in nutrient absorption, immune regulation, and cellular processes for milk synthesis during lactation.

Temporal Expression Analysis to Unravel Gene Regulatory Dynamics by microRNAs

MicroRNAs (miRNAs) are a class of small non-coding RNAs (sncRNAs) of length 21-25 nucleotides. These sncRNAs hybridize to repress their target genes and inhibit protein translation, thereby controlling regulatory functions in the cell. Integration of time-series matched small and RNA-seq data enables investigation of dynamic gene regulation through miRNAs during development or in response to a stimulus, such as stress. Here we summarize analysis strategies, such as probabilistic and regression-based models, that take advantage of the temporal dimension to investigate the complexity of miRNA regulation.

Identification and developmental profiling of microRNAs in the acoel worm Hofstenia miamia

The acoel worm Hofstenia miamia (H. miamia) has recently emerged as a model organism for studying whole-body regeneration and embryonic development. Previous studies suggest that post-transcriptional mechanisms likely play important roles in whole-body regeneration. Here, we establish a resource for studying H. miamia microRNA-mediated gene regulation, a major aspect of post-transcriptional control in animals. Using small RNA-sequencing samples spanning key developmental stages, we annotated H. miamia microRNAs. Our analysis uncovered a total of 1,050 microRNA loci, including 479 high-confidence loci based on structural and read abundance criteria. Comparison of microRNA seed sequences with those in other bilaterian species revealed that H. miamia encodes the majority of known conserved bilaterian microRNA families and that several microRNA families previously reported only in protostomes or deuterostomes likely have ancient bilaterian origins. We profiled the expression dynamics of the H. miamia miRNAs across embryonic and post-embryonic development. We observed that the let-7 and mir-125 microRNAs are unconventionally enriched at early embryonic stages. To generate hypotheses for miRNA function, we annotated the 3' UTRs of H. miamia protein-coding genes and performed miRNA target site predictions. Focusing on genes that are known to function in the wound response, posterior patterning, and neural differentiation in H. miamia , we found that these processes may be under substantial miRNA regulation. Notably, we found that miRNAs in MIR-7 and MIR-9 families which have target sites in the posterior genes fz-1 , wnt-3 , and sp5 are indeed expressed in the anterior of the animal, consistent with a repressive effect on their corresponding target genes. Our annotation offers candidate miRNAs for further functional investigation, providing a resource for future studies of post-transcriptional control during development and regeneration.

A miR-383-5p Signaling Hub Coordinates the Axon Regeneration Response to Inflammation

Neuroinflammation can positively influence axon regeneration following injury in the central nervous system. Inflammation promotes the release of neurotrophic molecules and stimulates intrinsic proregenerative molecular machinery in neurons, but the detailed mechanisms driving this effect are not fully understood. We evaluated how microRNAs are regulated in retinal neurons in response to intraocular inflammation to identify their potential role in axon regeneration. We found that miR-383-5p is downregulated in retinal ganglion cells in response to zymosan-induced intraocular inflammation. MiR-383-5p downregulation in neurons is sufficient to promote axon growth in vitro, and the intravitreal injection of a miR-383-5p inhibitor into the eye promotes axon regeneration following optic nerve crush. MiR-383-5p directly targets ciliary neurotrophic factor (CNTF) receptor components, and miR-383-5p inhibition sensitizes adult retinal neurons to the outgrowth-promoting effects of CNTF. Interestingly, we also demonstrate that CNTF treatment is sufficient to reduce miR-383-5p levels in neurons, constituting a positive-feedback module, whereby initial CNTF treatment reduces miR-383-5p levels, which then disinhibits CNTF receptor components to sensitize neurons to the ligand. Additionally, miR-383-5p inhibition derepresses the mitochondrial antioxidant protein peroxiredoxin-3 (PRDX3) which was required for the proregenerative effects associated with miR-383-5p loss-of-function in vitro. We have thus identified a positive-feedback mechanism that facilitates neuronal CNTF sensitivity in neurons and a new molecular signaling module that promotes inflammation-induced axon regeneration.

A gene desert required for regulatory control of pleiotropic Shox2 expression and embryonic survival

Approximately a quarter of the human genome consists of gene deserts, large regions devoid of genes often located adjacent to developmental genes and thought to contribute to their regulation. However, defining the regulatory functions embedded within these deserts is challenging due to their large size. Here, we explore the cis-regulatory architecture of a gene desert flanking the Shox2 gene, which encodes a transcription factor indispensable for proximal limb, craniofacial, and cardiac pacemaker development. We identify the gene desert as a regulatory hub containing more than 15 distinct enhancers recapitulating anatomical subdomains of Shox2 expression. Ablation of the gene desert leads to embryonic lethality due to Shox2 depletion in the cardiac sinus venosus, caused in part by the loss of a specific distal enhancer. The gene desert is also required for stylopod morphogenesis, mediated via distributed proximal limb enhancers. In summary, our study establishes a multi-layered role of the Shox2 gene desert in orchestrating pleiotropic developmental expression through modular arrangement and coordinated dynamics of tissue-specific enhancers.

[Micro-RNAs as biomarkers of radiation-induced injuries]

In the event of a radiological or nuclear emergency following an accident or malicious act, potentially involving many victims, medical care requires the identification and diagnosis of individuals exposed to high doses of ionizing radiation as quickly as possible. While an initial screening can be carried out directly in the field, additional biological in-lab analyses are required to refine the diagnosis and optimize the therapeutic management of victims. The fast and simultaneous management of many patients is limited by currently established techniques. To overcome these constraints, the use of new biomarkers to predict the risk and severity of radiation-induced injuries is under investigation. This synthesis summarizes the latest scientific advances demonstrating the potential of microRNAs as biomarkers of radiationinduced injuries, highlighting their relevance for human health care and radioprotection.

The ENCODE Uniform Analysis Pipelines

The Encyclopedia of DNA elements (ENCODE) project is a collaborative effort to create a comprehensive catalog of functional elements in the human genome. The current database comprises more than 19000 functional genomics experiments across more than 1000 cell lines and tissues using a wide array of experimental techniques to study the chromatin structure, regulatory and transcriptional landscape of the Homo sapiens and Mus musculus genomes. All experimental data, metadata, and associated computational analyses created by the ENCODE consortium are submitted to the Data Coordination Center (DCC) for validation, tracking, storage, and distribution to community resources and the scientific community. The ENCODE project has engineered and distributed uniform processing pipelines in order to promote data provenance and reproducibility as well as allow interoperability between genomic resources and other consortia. All data files, reference genome versions, software versions, and parameters used by the pipelines are captured and available via the ENCODE Portal. The pipeline code, developed using Docker and Workflow Description Language (WDL; https://openwdl.org/) is publicly available in GitHub, with images available on Dockerhub (https://hub.docker.com), enabling access to a diverse range of biomedical researchers. ENCODE pipelines maintained and used by the DCC can be installed to run on personal computers, local HPC clusters, or in cloud computing environments via Cromwell. Access to the pipelines and data via the cloud allows small labs the ability to use the data or software without access to institutional compute clusters. Standardization of the computational methodologies for analysis and quality control leads to comparable results from different ENCODE collections - a prerequisite for successful integrative analyses.

The ENCODE Uniform Analysis Pipelines

The Encyclopedia of DNA elements (ENCODE) project is a collaborative effort to create a comprehensive catalog of functional elements in the human genome. The current database comprises more than 19000 functional genomics experiments across more than 1000 cell lines and tissues using a wide array of experimental techniques to study the chromatin structure, regulatory and transcriptional landscape of the Homo sapiens and Mus musculus genomes. All experimental data, metadata, and associated computational analyses created by the ENCODE consortium are submitted to the Data Coordination Center (DCC) for validation, tracking, storage, and distribution to community resources and the scientific community. The ENCODE project has engineered and distributed uniform processing pipelines in order to promote data provenance and reproducibility as well as allow interoperability between genomic resources and other consortia. All data files, reference genome versions, software versions, and parameters used by the pipelines are captured and available via the ENCODE Portal. The pipeline code, developed using Docker and Workflow Description Language (WDL; https://openwdl.org/) is publicly available in GitHub, with images available on Dockerhub (https://hub.docker.com), enabling access to a diverse range of biomedical researchers. ENCODE pipelines maintained and used by the DCC can be installed to run on personal computers, local HPC clusters, or in cloud computing environments via Cromwell. Access to the pipelines and data via the cloud allows small labs the ability to use the data or software without access to institutional compute clusters. Standardization of the computational methodologies for analysis and quality control leads to comparable results from different ENCODE collections - a prerequisite for successful integrative analyses.

Promiscuous splicing-derived hairpins are dominant substrates of tailing-mediated defense of miRNA biogenesis in mammals

Canonical microRNA (miRNA) hairpins are processed by the RNase III enzymes Drosha and Dicer into ∼22 nt RNAs loaded into an Argonaute (Ago) effector. In addition, splicing generates numerous intronic hairpins that bypass Drosha (mirtrons) to yield mature miRNAs. Here, we identify hundreds of previously unannotated, splicing-derived hairpins in intermediate-length (∼50-100 nt) but not small (20-30 nt) RNA data. Since we originally defined mirtrons from small RNA duplexes, we term this larger set as structured splicing-derived RNAs (ssdRNAs). These associate with Dicer and/or Ago complexes, but generally accumulate modestly and are poorly conserved. We propose they contaminate the canonical miRNA pathway, which consequently requires defense against the siege of splicing-derived substrates. Accordingly, ssdRNAs/mirtrons comprise dominant hairpin substrates for 3' tailing by multiple terminal nucleotidyltransferases, notably TUT4/7 and TENT2. Overall, the rampant proliferation of young mammalian mirtrons/ssdRNAs, coupled with an inhibitory molecular defense, comprises a Red Queen's race of intragenomic conflict.

Circulating miRNAs in hand osteoarthritis

OBJECTIVE

Hand osteoarthritis (OA) is a frequent musculoskeletal disorder with an increasing prevalence during ageing. This study aimed to evaluate circulating microRNAs (miRNAs) in the plasma of patients with hand OA compared with age- and sex-matched healthy control subjects.

METHODS

In total, 238 participants (96 with erosive and 73 with non-erosive hand OA patients and 69 healthy control subjects) were included in this study. All patients underwent clinical examinations, including self-reported measures (AUSCAN and Algofunctional index). Radiographs of both hands were scored with the Kallman scale. The profile of miRNAs in plasma was screened using TaqMan™ Low-Density Array, and candidate miRNAs were validated on two quantitative real-time PCR (qRT-PCR) systems (QuantStudio and SmartChip).

RESULTS

Of all the 754 miRNAs, 40 miRNAs were different between hand OA patients and healthy control subjects in the screening cohort. Following the two-phase validation process, three miRNAs (miR-23a-3p, miR-146a-5p, and miR-652-3p) were increased in patients with hand OA compared with healthy control subjects and were associated with the AUSCAN sum score and AUSCAN pain. Furthermore, an inverse correlation of miR-222-3p with the Kallman radiographic score was found. The expression of miRNAs did not differ between erosive and non-erosive hand OA.

CONCLUSION

The profile of circulating miRNAs could unveil candidate biomarkers associated with hand OA symptoms. Longitudinal studies are required to determine the role of miRNAs in hand OA.

Canalization of Phenotypes-When the Transcriptome is Constantly but Weakly Perturbed

Recent studies have increasingly pointed to microRNAs (miRNAs) as the agent of gene regulatory network (GRN) stabilization as well as developmental canalization against constant but small environmental perturbations. To analyze mild perturbations, we construct a Dicer-1 knockdown line (dcr-1 KD) in Drosophila that modestly reduces all miRNAs by, on average, ∼20%. The defining characteristic of stabilizers is that, when their capacity is compromised, GRNs do not change their short-term behaviors. Indeed, even with such broad reductions across all miRNAs, the changes in the transcriptome are very modest during development in stable environment. By comparison, broad knockdowns of other regulatory genes (esp. transcription factors) by the same method should lead to drastic changes in the GRNs. The consequence of destabilization may thus be in long-term development as postulated by the theory of canalization. Flies with modest miRNA reductions may gradually deviate from the developmental norm, resulting in late-stage failures such as shortened longevity. In the optimal culture condition, the survival to adulthood is indeed normal in the dcr-1 KD line but, importantly, adult longevity is reduced by ∼90%. When flies are stressed by high temperature, dcr-1 KD induces lethality earlier in late pupation and, as the perturbations are shifted earlier, the affected stages are shifted correspondingly. Hence, in late stages of development with deviations piling up, GRN would be increasingly in need of stabilization. In conclusion, miRNAs appear to be a solution to weak but constant environmental perturbations.

The embryonic transcriptome of Parhyale hawaiensis reveals different dynamics of microRNAs and mRNAs during the maternal-zygotic transition

Parhyale hawaiensis has emerged as the crustacean model of choice due to its tractability, ease of imaging, sequenced genome, and development of CRISPR/Cas9 genome editing tools. However, transcriptomic datasets spanning embryonic development are lacking, and there is almost no annotation of non-protein-coding RNAs, including microRNAs. We have sequenced microRNAs, together with mRNAs and long non-coding RNAs, in Parhyale using paired size-selected RNA-seq libraries at seven time-points covering important transitions in embryonic development. Focussing on microRNAs, we annotate 175 loci in Parhyale, 88 of which have no known homologs. We use these data to annotate the microRNAome of 37 crustacean genomes, and suggest a core crustacean microRNA set of around 61 sequence families. We examine the dynamic expression of microRNAs and mRNAs during the maternal-zygotic transition. Our data suggest that zygotic genome activation occurs in two waves in Parhyale with microRNAs transcribed almost exclusively in the second wave. Contrary to findings in other arthropods, we do not predict a general role for microRNAs in clearing maternal transcripts. These data significantly expand the available transcriptomics resources for Parhyale, and facilitate its use as a model organism for the study of small RNAs in processes ranging from embryonic development to regeneration.

BET Proteins Regulate Expression of Osr1 in Early Kidney Development

In utero renal development is subject to maternal metabolic and environmental influences affecting long-term renal function and the risk of developing chronic kidney failure and cardiovascular disease. Epigenetic processes have been implicated in the orchestration of renal development and prenatal programming of nephron number. However, the role of many epigenetic modifiers for kidney development is still unclear. Bromodomain and extra-terminal domain (BET) proteins act as histone acetylation reader molecules and promote gene transcription. BET family members Brd2, Brd3 and Brd4 are expressed in the nephrogenic zone during kidney development. Here, the effect of the BET inhibitor JQ1 on renal development is evaluated. Inhibition of BET proteins via JQ1 leads to reduced growth of metanephric kidney cultures, loss of the nephron progenitor cell population, and premature and disturbed nephron differentiation. Gene expression of key nephron progenitor transcription factor Osr1 is downregulated after 24 h BET inhibition, while Lhx1 and Pax8 expression is increased. Mining of BRD4 ChIP-seq and gene expression data identify Osr1 as a key factor regulated by BRD4-controlled gene activation. Inhibition of BRD4 by BET inhibitor JQ1 leads to downregulation of Osr1, thereby causing a disturbance in the balance of nephron progenitor cell self-renewal and premature differentiation of the nephron, which ultimately leads to kidney hypoplasia and disturbed nephron development. This raises questions about the potential teratogenic effects of BET inhibitors for embryonic development. In summary, our work highlights the role of BET proteins for prenatal programming of nephrogenesis and identifies Osr1 as a potential target of BET proteins.

Signaling Modulation by miRNA-221-3p During Tooth Morphogenesis in Mice

miRNAs are conserved short non-coding RNAs that play a role in the modulation of various biological pathways during tissue and organ morphogenesis. In this study, the function of miRNA-221-3p in tooth development, through its loss or gain in function was evaluated. A variety of techniques were utilized to evaluate detailed functional roles of miRNA-221-3p during odontogenesis, including in vitro tooth cultivation, renal capsule transplantation, in situ hybridization, real-time PCR, and immunohistochemistry. Two-day in vitro tooth cultivation at E13 identified altered cellular events, including cellular proliferation, apoptosis, adhesion, and cytoskeletal arrangement, with the loss and gain of miRNA-221-3p. qPCR analysis revealed alterations in gene expression of tooth-related signaling molecules, including β-catenin, Bmp2, Bmp4, Fgf4, Ptch1, and Shh, when inhibited with miRNA-221-3p and mimic. Also, the inhibition of miRNA-221-3p demonstrated increased mesenchymal localizations of pSMAD1/5/8, alongside decreased expression patterns of Shh and Fgf4 within inner enamel epithelium (IEE) in E13 + 2 days in vitro cultivated teeth. Moreover, 1-week renal transplantation of in vitro cultivated teeth had smaller tooth size with reduced enamel and dentin matrices, along with increased cellular proliferation and Shh expression along the Hertwig epithelial root sheath (HERS), within the inhibitor group. Similarly, in 3-week renal calcified teeth, the overexpression of miRNA-221-3p did not affect tooth phenotype, while the loss of function resulted in long and slender teeth with short mesiodistal length. This study provides evidence that a suitable level of miRNA-221-3p is required for the modulation of major signaling pathways, including Wnt, Bmp, and Shh, during tooth morphogenesis.

The lncRNA XIST/miR-150-5p/c-Fos axis regulates sepsis-induced myocardial injury via TXNIP-modulated pyroptosis

Myocardial injury is a severe complication of sepsis and contributes substantially to the death of critically ill patients. Long non-coding RNAs (lncRNAs) participate in the pathogenesis of sepsis-induced myocardial injury. In this study, we investigated the role of lncRNA X-inactive specific transcript (XIST) in septic myocardial injury and explored its mechanism. Lipopolysaccharide (LPS)-stimulated H9C2 cells and rats subjected to cecal ligation and puncture (CLP) were used as the in vitro and in vivo models. After exposure to LPS, XIST and c-Fos levels were upregulated, but miR-150-5p was downregulated in H9C2 cardiomyocytes and myocardial tissues. XIST affected viability, apoptosis, and pyroptosis in LPS-challenged H9C2 cells. Moreover, XIST knockdown attenuated LPS-induced injury in H9C2 cells by targeting the miR-150-5p/c-Fos axis. c-Fos could bound to the promoter of the TXNIP/XIST gene and enhanced TXNIP/XIST expression. Silencing of XIST improved cardiac function and survival rate and reduced apoptosis and pyroptosis by regulating the miR-150-5p/c-Fos axis in septic rats in vivo. Taken together, our data show that XIST/miR-150-5p/c-Fos axis affected septic myocardial injury, which may indicate a novel therapeutic strategy for sepsis-induced myocardial injury.

Distinct and Coordinated Regulation of Small Non-coding RNAs by E2f1 and p53 During Drosophila Development and in Response to DNA Damage

Small non-coding RNAs (ncRNAs), including microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs), play a pivotal role in biological processes. A comprehensive quantitative reference of small ncRNAs expression during development and in DNA damage response (DDR) would significantly advance our understanding of their roles. In this study, we systemically analyzed the expression profile of miRNAs and piRNAs in wild-type flies, e2f1 mutant, p53 mutant and e2f1 p53 double mutant during development and after X-ray irradiation. By using small RNA sequencing and bioinformatic analysis, we found that both miRNAs and piRNAs were expressed in a dynamic mode and formed 4 distinct clusters during development. Notably, the expression pattern of miRNAs and piRNAs was changed in e2f1 mutant at multiple developmental stages, while retained in p53 mutant, indicating a critical role of E2f1 played in mediating small ncRNAs expression. Moreover, we identified differentially expressed (DE) small ncRNAs in e2f1 mutant and p53 mutant after X-ray irradiation. Furthermore, we mapped the binding motif of E2f1 and p53 around the small ncRNAs. Our data suggested that E2f1 and p53 work differently yet coordinately to regulate small ncRNAs expression, and E2f1 may play a major role to regulate miRNAs during development and after X-ray irradiation. Collectively, our results provide comprehensive characterization of small ncRNAs, as well as the regulatory roles of E2f1 and p53 in small ncRNAs expression, during development and in DNA damage response, which reveal new insights into the small ncRNAs biology.

Novel insights into the interplay between m6A modification and noncoding RNAs in cancer

N6-methyladenosine (m⁶A) is one of the most common RNA modifications in eukaryotes, mainly in messenger RNA (mRNA). Increasing evidence shows that m⁶A methylation modification acts an essential role in various physiological and pathological bioprocesses. Noncoding RNAs (ncRNAs), including miRNAs, lncRNAs and circRNAs, are known to participate in regulating cell differentiation, angiogenesis, immune response, inflammatory response and carcinogenesis. m⁶A regulators, such as METTL3, ALKBH5 and IGF2BP1 have been reported to execute a m⁶A-dependent modification of ncRNAs involved in carcinogenesis. Meanwhile, ncRNAs can target or modulate m⁶A regulators to influence cancer development. In this review, we provide an insight into the interplay between m⁶A modification and ncRNAs in cancer.

Expanded encyclopaedias of DNA elements in the human and mouse genomes

The human and mouse genomes contain instructions that specify RNAs and proteins and govern the timing, magnitude, and cellular context of their production. To better delineate these elements, phase III of the Encyclopedia of DNA Elements (ENCODE) Project has expanded analysis of the cell and tissue repertoires of RNA transcription, chromatin structure and modification, DNA methylation, chromatin looping, and occupancy by transcription factors and RNA-binding proteins. Here we summarize these efforts, which have produced 5,992 new experimental datasets, including systematic determinations across mouse fetal development. All data are available through the ENCODE data portal (https://www.encodeproject.org), including phase II ENCODE1 and Roadmap Epigenomics2 data. We have developed a registry of 926,535 human and 339,815 mouse candidate cis-regulatory elements, covering 7.9 and 3.4% of their respective genomes, by integrating selected datatypes associated with gene regulation, and constructed a web-based server (SCREEN; http://screen.encodeproject.org) to provide flexible, user-defined access to this resource. Collectively, the ENCODE data and registry provide an expansive resource for the scientific community to build a better understanding of the organization and function of the human and mouse genomes.

The changing mouse embryo transcriptome at whole tissue and single-cell resolution

During mammalian embryogenesis, differential gene expression gradually builds the identity and complexity of each tissue and organ system1. Here we systematically quantified mouse polyA-RNA from day 10.5 of embryonic development to birth, sampling 17 tissues and organs. The resulting developmental transcriptome is globally structured by dynamic cytodifferentiation, body-axis and cell-proliferation gene sets that were further characterized by the transcription factor motif codes of their promoters. We decomposed the tissue-level transcriptome using single-cell RNA-seq (sequencing of RNA reverse transcribed into cDNA) and found that neurogenesis and haematopoiesis dominate at both the gene and cellular levels, jointly accounting for one-third of differential gene expression and more than 40% of identified cell types. By integrating promoter sequence motifs with companion ENCODE epigenomic profiles, we identified a prominent promoter de-repression mechanism in neuronal expression clusters that was attributable to known and novel repressors. Focusing on the developing limb, single-cell RNA data identified 25 candidate cell types that included progenitor and differentiating states with computationally inferred lineage relationships. We extracted cell-type transcription factor networks and complementary sets of candidate enhancer elements by using single-cell RNA-seq to decompose integrative cis-element (IDEAS) models that were derived from whole-tissue epigenome chromatin data. These ENCODE reference data, computed network components and IDEAS chromatin segmentations are companion resources to the matching epigenomic developmental matrix, and are available for researchers to further mine and integrate.