Dissecting gene expression networks in the developing hippocampus through the lens of NEIL3 depletion.

Gene regulation in the hippocampus is fundamental for its development, synaptic plasticity, memory formation, and adaptability. Comparisons of gene expression among different developmental stages, distinct cell types, and specific experimental conditions have identified differentially expressed genes contributing to the organization and functionality of hippocampal circuits. The NEIL3 DNA glycosylase, one of the DNA repair enzymes, plays an important role in hippocampal maturation and neuron functionality by shaping transcription. While differential gene expression (DGE) analysis has identified key genes involved, broader gene expression patterns crucial for high-order hippocampal functions remain uncharted. By utilizing the weighted gene co-expression network analysis (WGCNA), we mapped gene expression networks in immature (p8-neonatal) and mature (3 m-adult) hippocampal circuits in wild-type and NEIL3-deficient mice. Our study unveiled intricate gene network structures underlying hippocampal maturation, delineated modules of co-expressed genes, and pinpointed highly interconnected hub genes specific to the maturity of hippocampal subregions. We investigated variations within distinct gene network modules following NEIL3 depletion, uncovering NEIL3-targeted hub genes that influence module connectivity and specificity. By integrating WGCNA with DGE, we delve deeper into the NEIL3-dependent molecular intricacies of hippocampal maturation. This study provides a comprehensive systems-level analysis for assessing the potential correlation between gene connectivity and functional connectivity within the hippocampal network, thus shaping hippocampal function throughout development.

Enhanced glutathione levels confer resistance to apoptotic and ferroptotic programmed cell death in NEIL DNA glycosylase deficient HAP1 cells.

The NTHL1 and NEIL1-3 DNA glycosylases are major enzymes in the removal of oxidative DNA base lesions, via the base excision repair (BER) pathway. It is expected that lack of these DNA glycosylases activities would render cells vulnerable to oxidative stress, promoting cell death. Intriguingly, we found that single, double, triple, and quadruple DNA glycosylase knockout HAP1 cells are, however, more resistant to oxidative stress caused by genotoxic agents than wild type cells. Furthermore, glutathione depletion in NEIL deficient cells further enhances resistance to cell death induced via apoptosis and ferroptosis. Finally, we observed higher basal level of glutathione and differential expression of NRF2-regulated genes associated with glutathione homeostasis in the NEIL triple KO cells. We propose that lack of NEIL DNA glycosylases causes aberrant transcription and subsequent errors in protein synthesis. This leads to increased endoplasmic reticulum stress and proteotoxic stress. To counteract the elevated intracellular stress, an adaptive response mediated by increased glutathione basal levels, rises in these cells. This study reveals an unforeseen role of NEIL glycosylases in regulation of resistance to oxidative stress, suggesting that modulation of NEIL glycosylase activities is a potential approach to improve the efficacy of e.g. anti-inflammatory therapies.

Loss of alkyladenine DNA glycosylase alters gene expression in the developing mouse brain and leads to reduced anxiety and improved memory.

Neurodevelopment is a tightly coordinated process, during which the genome is exposed to spectra of endogenous agents at different stages of differentiation. Emerging evidence indicates that DNA damage is an important feature of developing brain, tightly linked to gene expression and neuronal activity. Some of the most frequent DNA damage includes changes to DNA bases, which are recognized by DNA glycosylases and repaired through base excision repair (BER) pathway. The only mammalian DNA glycosylase able to remove frequent alkylated DNA based is alkyladenine DNA glycosylase (Aag, aka Mpg). We recently demonstrated that, besides its role in DNA repair, AAG affects expression of neurodevelopmental genes in human cells. Aag was further proposed to act as reader of epigenetic marks, including 5-hydroxymethylcytosine (5hmC), in the mouse brain. Despite the potential Aag involvement in the key brain processes, the impact of Aag loss on developing brain remains unknown. Here, by using Aag knockout (Aag-/-) mice, we show that Aag absence leads to reduced DNA break levels, evident in lowered number of γH2AX foci in postnatal day 5 (P5) hippocampi. This is accompanied by changes in 5hmC signal intensity in different hippocampal regions. Transcriptome analysis of hippocampi and prefrontal cortex, at different developmental stages, indicates that lack of Aag alters gene expression, primarily of genes involved in regulation of response to stress. Across all developmental stages tested aldehyde dehydrogenase 2 (Aldh2) emerged as one of the most prominent genes deregulated in Aag-dependent manner. In line with the changes in hippocampal DNA damage levels and the gene expression, adult Aag-/- mice exhibit altered behavior, evident in decreased anxiety levels determined in the Elevated Zero Maze and increased alternations in the Elevated T Maze tests. Taken together these results suggests that Aag has functions in modulation of genome dynamics during brain development, important for animal behavior.

Postnatal persistence of hippocampal Cajal-Retzius cells has a crucial role in the establishment of the hippocampal circuit.

Cajal-Retzius (CR) cells are a transient neuron type that populate the postnatal hippocampus. To understand how the persistence of CR cells influences the maturation of hippocampal circuits, we combined a specific transgenic mouse line with viral vector injection to selectively ablate CR cells from the postnatal hippocampus. We observed layer-specific changes in the dendritic complexity and spine density of CA1 pyramidal cells. In addition, transcriptomic analysis highlighted significant changes in the expression of synapse-related genes across development. Finally, we were able to identify significant changes in the expression levels of latrophilin 2, a postsynaptic guidance molecule known for its role in the entorhinal-hippocampal connectivity. These findings were supported by changes in the synaptic proteomic content in CA1 stratum lacunosum-moleculare. Our results reveal a crucial role for CR cells in the establishment of the hippocampal network.

Single-cell resolution of longitudinal blood transcriptome profiles in rheumatoid arthritis, systemic lupus erythematosus and healthy control pregnancies.

OBJECTIVES: Comparative longitudinal analyses of cellular composition and peripheral blood gene expression in Rheumatoid arthritis (RA), systemic lupus erythematosus (SLE) and healthy pregnancies. METHODS: In total, 335 whole blood samples from 84 RA, SLE and healthy controls before pregnancy, at each trimester, 6 weeks, 6 months and 12 months post partum were analysed. We combined bulk and single cell RNA analyses for cell-type estimation, validated by flow cytometry, before combining this in a cell-type adjusted analysis for an improved resolution of unrecognised gene expression changes associated with RA and SLE pregnancies. RESULTS: Patients were well regulated throughout pregnancy, and few had pregnancy complications. In SLE, the interferon signature was augmented during pregnancy, and the pregnancy signature was continued post partum. An altered cell type composition strongly influences the profile. In the pregnancy signature, transcripts involved in galactosylation potentially altering the effector functions of autoantibodies became more evident. Several genes in the adjusted RA signature are expressed in mucosal associated invariant T cells. CONCLUSION: We found distinct RA, SLE and pregnancy signatures, and no expression patterns could be attributed to medication or disease activity. Our results support the need for close postpartum follow-up of patients with SLE. Gene expression patterns in RA were closer to healthy controls than to SLE, and primarily became evident after cell-type adjustment. Adjusting for cell abundance unravelled gene expression signatures less associated with variation in cell-composition and highlighted genes with expression profiles associated with changes in specialised cell populations.

The Salmon Oil OmeGo Reduces Viability of Colorectal Cancer Cells and Potentiates the Anti-Cancer Effect of 5-FU.

Colorectal cancer (CRC) is one of the most common cancer types worldwide. Chemotherapy is toxic to normal cells, and combinatory treatment with natural well-tolerated products is being explored. Some omega-3 polyunsaturated fatty acids (n-3 PUFAs) and marine fish oils have anti-cancer effects on CRC cells. The salmon oil OmeGo (Hofseth BioCare) contains a spectrum of fatty acids, including the n-3 PUFAs docosahexaenoic acid (DHA) and eicosahexaenoic acid (EPA). We explored a potential anti-cancer effect of OmeGo on the four CRC cell lines DLD-1, HCT-8, LS411N, and LS513, alone and in combination with the chemotherapeutic agent 5-Fluorouracil (5-FU). Screening indicated a time- and dose-dependent effect of OmeGo on the viability of the DLD-1 and LS513 CRC cell lines. Treatment with 5-FU and OmeGo (IC20-IC30) alone indicated a significant reduction in viability. A combinatory treatment with OmeGo and 5-FU resulted in a further reduction in viability in DLD-1 and LS513 cells. Treatment of CRC cells with DHA + EPA in a concentration corresponding to the content in OmeGo alone or combined with 5-FU significantly reduced viability of all four CRC cell lines tested. The lowest concentration of OmeGo reduced viability to a higher degree both alone and in combination with 5-FU compared to the corresponding concentrations of DHA + EPA in three of the cell lines. Results suggest that a combination of OmeGo and 5-FU could have a potential as an alternative anti-cancer therapy for patients with CRC.

DNA hydroxymethylation differences underlie phenotypic divergence of somatic growth in Nile tilapia reared in common garden.

Phenotypic plasticity of metabolism and growth are essential for adaptation to new environmental conditions, such as those experienced during domestication. Epigenetic regulation plays a key role in this process but the underlying mechanisms are poorly understood, especially in the case of hydroxymethylation. Using reduced representation 5-hydroxymethylcytosine profiling, we compared the liver hydroxymethylomes in full-sib Nile tilapia with distinct growth rates (3.8-fold difference) and demonstrated that DNA hydroxymethylation is strongly associated with phenotypic divergence of somatic growth during the early stages of domestication. The 2677 differentially hydroxymethylated cytosines between fast- and slow-growing fish were enriched within gene bodies (79%), indicating a pertinent role in transcriptional regulation. Moreover, they were found in genes involved in biological processes related to skeletal system and muscle structure development, and there was a positive association between somatic growth and 5hmC levels in genes coding for growth factors, kinases and receptors linked to myogenesis. Single nucleotide polymorphism analysis revealed no genetic differentiation between fast- and slow-growing fish. In addition to unveiling a new link between DNA hydroxymethylation and epigenetic regulation of growth in fish during the initial stages of domestication, this study suggests that epimarkers may be applied in selective breeding programmes for superior phenotypes.

Increased expression of individual genes in whole blood is associated with late-stage lung cancer at and close to diagnosis.

Lung cancer (LC) mortality rates are still increasing globally. As survival is linked to stage, there is a need to identify markers for earlier LC diagnosis and individualized treatment. The whole blood transcriptome of LC patients represents a source of potential LC biomarkers. We compared expression of > 60,000 genes in whole blood specimens taken from LC cases at diagnosis (n = 128) and controls (n = 62) using genome-wide RNA sequencing, and identified 14 candidate genes associated with LC. High expression of ANXA3, ARG1 and HP was strongly associated with lower survival in late-stage LC cases (hazard ratios (HRs) = 2.81, 2.16 and 2.54, respectively). We validated these markers in two independent population-based studies with pre-diagnostic whole blood specimens taken up to eight years prior to LC diagnosis (n = 163 cases, 184 matched controls). ANXA3 and ARG1 expression was strongly associated with LC in these specimens, especially with late-stage LC within two years of diagnosis (odds ratios (ORs) = 3.47 and 5.00, respectively). Additionally, blood CD4 T cells, NK cells and neutrophils were associated with LC at diagnosis and improved LC discriminative ability beyond candidate genes. Our results indicate that in whole blood, increased expression levels of ANXA3, ARG1 and HP are diagnostic and prognostic markers of late-stage LC.

Randomised controlled trial of exercise training during lactation on breast milk composition in breastfeeding people with overweight/obesity: a study protocol for the MILKSHAKE trial.

Breast milk from people with overweight/obesity may differ in composition compared with that from normal-weight people. Exercise training can modify breast milk composition in rodent models, with a beneficial impact demonstrated on the offspring's metabolism, but whether these findings translate to humans is unclear. This trial aims to determine the effect of an exercise intervention on breast milk composition and whether an exercise-induced modification of breast milk impacts the infants' growth and body composition. Effect of Exercise Training on Breastmilk Composition is a randomised, controlled trial with two parallel groups, one exercise group and one control group, with a 1:1 allocation. We will include a minimum of 62 exclusively breastfeeding participants, 6 weeks postpartum. The exercise intervention lasts 8 weeks and comprises 25 supervised endurance exercise sessions with moderate or high intensity. The primary outcome measure is the change in the relative concentration of the human milk oligosaccharide 3'sialyllactose in breast milk from baseline at 6 weeks postpartum to the end of the intervention period. Secondary outcomes include breast milk concentrations of other metabolites, cytokines, hormones and microRNA, maternal health outcomes, infant growth, infant gut microbiome and infant circulating microRNA. Maternal and infant outcomes will be measured before, during and after the intervention period, with a follow-up of the infants until they are 24 months old. Trial registration number NCT05488964.

Associations of gene expression in blood with BMI and weight changes among women in the Norwegian Women and Cancer postgenome cohort.

OBJECTIVE: This study aimed to evaluate associations between blood gene expression profiles and (1) current BMI and (2) past weight changes (WCs) among women who had never been diagnosed with cancer in the Norwegian Women and Cancer (NOWAC) postgenome cohort. METHODS: This cross-sectional study (N = 1694) used gene expression profiles and information from three questionnaires: Q1 (baseline), Q2 (follow-up), and Q3 (blood collection). The authors performed gene-wise linear regression models to identify differentially expressed genes (DEGs) and functional enrichment analyses to identify their biological functions. RESULTS: When assessing BMIQ3 , the study observed 2394, 769, and 768 DEGs for the obesity-versus-normal weight, obesity-versus-overweight, and overweight-versus-normal weight comparisons, respectively. Up to 169 DEGs were observed when investigating WCQ3-Q1 (mean = 7 years, range = 5.5-14 years) and WCQ3-Q2 (mean = 1 year, range = <1 month-9 years) in interaction models with BMI categories, of which 1 to 169 genes were associated with WCs and 0 to 9 were associated with interaction effects of BMI and WCs. Biological functions of BMI-associated DEGs were linked to metabolism, erythrocytes, oxidative stress, and immune processes, whereas WC-associated DEGs were linked to signal transduction. CONCLUSIONS: Many BMI-associated but few WC-associated DEGs were identified in the blood of women in Norway. The biological functions of BMI-associated DEGs likely reflect systemic impacts of obesity, especially blood reticulocyte-erythrocyte ratio shifts.

Transcriptional development of the hippocampus and the dorsal-intermediate-ventral axis in rats.

Risk and resilience for neuropsychiatric illnesses are established during brain development, and transcriptional markers of risk may be identifiable in early development. The dorsal-ventral axis of the hippocampus has behavioral, electrophysiological, anatomical, and transcriptional gradients and abnormal hippocampus development is associated with autism, schizophrenia, epilepsy, and mood disorders. We previously showed that differential gene expression along the dorsoventral hippocampus in rats was present at birth (postnatal day 0, P0), and that a subset of differentially expressed genes (DEGs) was present at all postnatal ages examined (P0, P9, P18, and P60). Here, we extend the analysis of that gene expression data to understand the development of the hippocampus as a whole by examining DEGs that change with age. We additionally examine development of the dorsoventral axis by looking at DEGs along the axis at each age. Using both unsupervised and supervised analyses, we find that the majority of DEGs are present from P0 to P18, with many expression profiles presenting peaks or dips at P9/18. During development of the hippocampus, enriched pathways associated with learning, memory, and cognition increase with age, as do pathways associated with neurotransmission and synaptic function. Development of the dorsoventral axis is greatest at P9 and P18 and is marked by DEGs associated with metabolic functions. Our data indicate that neurodevelopmental disorders like epilepsy, schizophrenia and affective disorders are enriched with developmental DEGs in the hippocampus, regardless of dorsoventral location, with the greatest enrichment of these clinical disorders seen in genes whose expression changes from P0-9. When comparing DEGs from the ventral and dorsal poles, the greatest number of neurodevelopmental disorders is enriched with DEGs found at P18. Taken together, the developing hippocampus undergoes substantial transcriptional maturation during early postnatal development, with expression of genes involved in neurodevelopmental disorders also showing maximal expression changes within this developmental period.

Genome-wide hydroxymethylation profiles in liver of female Nile tilapia with distinct growth performance.

The mechanisms underlying the fast genome evolution that occurs during animal domestication are poorly understood. Here, we present a genome-wide epigenetic dataset that quantifies DNA hydroxymethylation at single nucleotide resolution among full-sib Nile tilapia (Oreochromis niloticus) with distinct growth performance. In total, we obtained 355 million, 75 bp reads from 5 large- and 5 small-sized fish on an Illumina NextSeq500 platform. We identified several growth-related genes to be differentially hydroxymethylated, especially within gene bodies and promoters. Previously, we proposed that DNA hydroxymethylation greatly affects the earliest responses to adaptation and potentially drives genome evolution through its targeted enrichment and elevated nucleotide transversion rates. This dataset can be analysed in various contexts (e.g., epigenetics, evolution and growth) and compared to other epigenomic datasets in the future, namely DNA methylation and histone modifications. With forthcoming advancements in genome research, this hydroxymethylation dataset will also contribute to better understand the epigenetic regulation of key genomic features, such as cis-regulatory and transposable elements.

Characterization of POLE c.1373A > T p.(Tyr458Phe), causing high cancer risk.

The cancer syndrome polymerase proofreading-associated polyposis results from germline mutations in the POLE and POLD1 genes. Mutations in the exonuclease domain of these genes are associated with hyper- and ultra-mutated tumors with a predominance of base substitutions resulting from faulty proofreading during DNA replication. When a new variant is identified by gene testing of POLE and POLD1, it is important to verify whether the variant is associated with PPAP or not, to guide genetic counseling of mutation carriers. In 2015, we reported the likely pathogenic (class 4) germline POLE c.1373A > T p.(Tyr458Phe) variant and we have now characterized this variant to verify that it is a class 5 pathogenic variant. For this purpose, we investigated (1) mutator phenotype in tumors from two carriers, (2) mutation frequency in cell-based mutagenesis assays, and (3) structural consequences based on protein modeling. Whole-exome sequencing of two tumors identified an ultra-mutator phenotype with a predominance of base substitutions, the majority of which are C > T. A SupF mutagenesis assay revealed increased mutation frequency in cells overexpressing the variant of interest as well as in isogenic cells encoding the variant. Moreover, exonuclease repair yeast-based assay supported defect in proofreading activity. Lastly, we present a homology model of human POLE to demonstrate structural consequences leading to pathogenic impact of the p.(Tyr458Phe) mutation. The three lines of evidence, taken together with updated co-segregation and previously published data, allow the germline variant POLE c.1373A > T p.(Tyr458Phe) to be reclassified as a class 5 variant. That means the variant is associated with PPAP.

TiSA: TimeSeriesAnalysis-a pipeline for the analysis of longitudinal transcriptomics data.

Improved transcriptomic sequencing technologies now make it possible to perform longitudinal experiments, thus generating a large amount of data. Currently, there are no dedicated or comprehensive methods for the analysis of these experiments. In this article, we describe our TimeSeries Analysis pipeline (TiSA) which combines differential gene expression, clustering based on recursive thresholding, and a functional enrichment analysis. Differential gene expression is performed for both the temporal and conditional axes. Clustering is performed on the identified differentially expressed genes, with each cluster being evaluated using a functional enrichment analysis. We show that TiSA can be used to analyse longitudinal transcriptomic data from both microarrays and RNA-seq, as well as small, large, and/or datasets with missing data points. The tested datasets ranged in complexity, some originating from cell lines while another was from a longitudinal experiment of severity in COVID-19 patients. We have also included custom figures to aid with the biological interpretation of the data, these plots include Principal Component Analyses, Multi Dimensional Scaling plots, functional enrichment dotplots, trajectory plots, and complex heatmaps showing the broad overview of results. To date, TiSA is the first pipeline to provide an easy solution to the analysis of longitudinal transcriptomics experiments.

Increased levels of microRNA-320 in blood serum and plasma is associated with imminent and advanced lung cancer.

Lung cancer (LC) incidence is increasing globally and altered levels of microRNAs (miRNAs) in blood may contribute to identification of individuals with LC. We identified miRNAs differentially expressed in peripheral blood at LC diagnosis and evaluated, in pre-diagnostic blood specimens, how long before diagnosis expression changes in such candidate miRNAs could be detected. We identified upregulated candidate miRNAs in plasma specimens from a hospital-based study sample of 128 patients with confirmed LC and 62 individuals with suspected but confirmed negative LC (FalsePos). We then evaluated the expression of candidate miRNAs in pre-diagnostic plasma or serum specimens of 360 future LC cases and 375 matched controls. There were 1663 miRNAs detected in diagnostic specimens, nine of which met our criteria for candidate miRNAs. Higher expression of three candidates, miR-320b, 320c, and 320d, was associated with poor survival, independent of LC stage and subtype. Moreover, miR-320c and miR-320d expression was higher in pre-diagnostic specimens collected within 2 years of LC diagnosis. Our results indicated that elevated levels of miR-320c and miR-320d may be early indications of imminent and advanced LC.

MicroRNA profiling of psoriatic skin identifies 11 miRNAs associated with disease severity.

MicroRNAs (miRNAs) are small non-coding RNAs that have emerged as central regulators of gene expression and powerful biomarkers of disease. Much is yet unknown about their role in psoriasis pathology. To globally characterize the miRNAome of psoriatic skin, skin biopsies were collected from psoriatic cases (n = 75) and non-psoriatic controls (n = 46) and RNA sequenced. Count data were meta-analysed with a previously published dataset (cases, n = 24, controls, n = 20), increasing the number of psoriatic cases fourfold from previously published studies. Differential gene expression analyses were performed comparing lesional psoriatic (PP), non-lesional psoriatic (PN) and control (NN) skin. Further, functional enrichment and cell-specific analyses were performed. Across all contrasts, we identified 439 significantly differentially expressed miRNAs (DEMs), of which 85 were novel for psoriasis and 11 were related to disease severity. Meta-analyses identified 20 DEMs between PN and NN, suggesting an inherent change in the constitution of all skin in psoriasis. By integrating the miRNA transcriptome with mRNA target interactions, we identified several functionally enriched terms, including "thyroid hormone signalling," "insulin resistance" and various infectious diseases. Cell-specific expression analyses revealed that the upregulated DEMs were enriched in epithelial and immune cells. This study provides the most comprehensive overview of the miRNAome in psoriatic skin to date and identifies a miRNA signature related to psoriasis severity. Our results may represent molecular links between psoriasis and related comorbidities and have outlined potential directions for future functional studies to identify biomarkers and treatment targets.

NEIL1 and NEIL2 DNA glycosylases modulate anxiety and learning in a cooperative manner in mice.

Oxidative DNA damage in the brain has been implicated in neurodegeneration and cognitive decline. DNA glycosylases initiate base excision repair (BER), the main pathway for oxidative DNA base lesion repair. NEIL1 and NEIL3 DNA glycosylases affect cognition in mice, while the role of NEIL2 remains unclear. Here, we investigate the impact of NEIL2 and its potential overlap with NEIL1 on behavior in knockout mouse models. Neil1-/-Neil2-/- mice display hyperactivity, reduced anxiety and improved learning. Hippocampal oxidative DNA base lesion levels are comparable between genotypes and no mutator phenotype is found. Thus, impaired canonical repair is not likely to explain the altered behavior. Electrophysiology suggests reduced axonal activation in the hippocampal CA1 region in Neil1-/-Neil2-/- mice and lack of NEIL1 and NEIL2 causes dysregulation of genes in CA1 relevant for synaptic function. We postulate a cooperative function of NEIL1 and NEIL2 in genome regulation, beyond canonical BER, modulating behavior in mice.

Non-Coding RNAs in Human Breast Milk: A Systematic Review.

Breast milk is the primary source of nutrition and hydration for the newborn infant but also plays an important role in the child's first immune defense. Additionally, several breast milk factors have been implicated in immune-related health outcomes later in life, including immunoglobulins, cytokines, chemokines, growth factors and, more recently, non-coding RNA (ncRNA) species. In this systematic review, we provide a comprehensive summary of the current literature on endogenous ncRNAs found in human breast milk. Thirty (30) relevant studies were identified and, whilst the majority studies focused on microRNAs (miRNAs), there is evidence that breast milk contains high quantities of RNA which also include long-coding RNAs, circular RNAs, as well as other short RNAs and fragmented tRNA and rRNAs. Among studies investigating miRNAs, miR-148a-3p, miR-30a/d-5p, miR-22-3p, miR-146b-5p, miR-200a/c-3p, and the 5p end of the let-7 miRNAs were commonly reported among the top 10 miRNAs in the cell, lipid, and skim milk fractions of breast milk. Methodological difference and small sample sizes limit the possibility of conclusively identifying which maternal and infant characteristics affect the miRNA profile. The highly expressed miRNAs were generally reported to be similar across lactational stage, milk fraction, maternal and infant characteristics, or infant growth and health. All the same, individual studies identify potential differences in miRNA expression levels which should be confirmed by future studies. Stability, uptake, and physiological functions of miRNAs were also considered in several studies. Breast milk miRNAs are relatively resistant to a range of harsh conditions and uptake experiments suggest that extracellular vesicles containing miRNAs and circular RNAs can be taken up by intestinal epithelial cells. Although the evidence regarding the functional effect of breast milk miRNAs is limited, the predicted functions range from metabolic and biosynthetic processes to signaling pathways, cellular adhesion, communication, growth, and differentiation. Finally, this systematic review highlights some of the methodological challenges and knowledge gaps which can help direct future research in this field. In particular, it is important to further investigate the bioavailability of miRNAs in different milk fractions, and to characterize other ncRNAs which are largely unstudied. Systematic Review Registration: PROSPERO https://www.crd.york.ac.uk/prospero/display_record.php?RecordID=138989, identifier CRD42020138989.

Joint changes in RNA, RNA polymerase II, and promoter activity through the cell cycle identify non-coding RNAs involved in proliferation.

Proper regulation of the cell cycle is necessary for normal growth and development of all organisms. Conversely, altered cell cycle regulation often underlies proliferative diseases such as cancer. Long non-coding RNAs (lncRNAs) are recognized as important regulators of gene expression and are often found dysregulated in diseases, including cancers. However, identifying lncRNAs with cell cycle functions is challenging due to their often low and cell-type specific expression. We present a highly effective method that analyses changes in promoter activity, transcription, and RNA levels for identifying genes enriched for cell cycle functions. Specifically, by combining RNA sequencing with ChIP sequencing through the cell cycle of synchronized human keratinocytes, we identified 1009 genes with cell cycle-dependent expression and correlated changes in RNA polymerase II occupancy or promoter activity as measured by histone 3 lysine 4 trimethylation (H3K4me3). These genes were highly enriched for genes with known cell cycle functions and included 57 lncRNAs. We selected four of these lncRNAs-SNHG26, EMSLR, ZFAS1, and EPB41L4A-AS1-for further experimental validation and found that knockdown of each of the four lncRNAs affected cell cycle phase distributions and reduced proliferation in multiple cell lines. These results show that many genes with cell cycle functions have concomitant cell-cycle dependent changes in promoter activity, transcription, and RNA levels and support that our multi-omics method is well suited for identifying lncRNAs involved in the cell cycle.