Age prediction from human blood plasma using proteomic and small RNA data: a comparative analysis

Aging clocks, built from comprehensive molecular data, have emerged as promising tools in medicine, forensics, and ecological research. However, few studies have compared the suitability of different molecular data types to predict age in the same cohort and whether combining them would improve predictions. Here, we explored this at the level of proteins and small RNAs in 103 human blood plasma samples. First, we used a two-step mass spectrometry approach measuring 612 proteins to select and quantify 21 proteins that changed in abundance with age. Notably, proteins increasing with age were enriched for components of the complement system. Next, we used small RNA sequencing to select and quantify a set of 315 small RNAs that changed in abundance with age. Most of these were microRNAs (miRNAs), downregulated with age, and predicted to target genes related to growth, cancer, and senescence. Finally, we used the collected data to build age-predictive models. Among the different types of molecules, proteins yielded the most accurate model (R² = 0.59 ± 0.02), followed by miRNAs as the best-performing class of small RNAs (R² = 0.54 ± 0.02). Interestingly, the use of protein and miRNA data together improved predictions (R2 = 0.70 ± 0.01). Future work using larger sample sizes and a validation dataset will be necessary to confirm these results. Nevertheless, our study suggests that combining proteomic and miRNA data yields superior age predictions, possibly by capturing a broader range of age-related physiological changes. It will be interesting to determine if combining different molecular data types works as a general strategy to improve future aging clocks.

Corrigendum: Identification of extremely GC-rich micro RNAs for RT-qPCR data normalization in human plasma

[This corrects the article DOI: 10.3389/fgene.2022.1058668.].

Contribution of mutant HSC clones to immature and mature cells in MDS and CMML, and variations with AZA therapy

Myelodysplastic neoplasms (MDSs) and chronic myelomonocytic leukemia (CMML) are clonal disorders driven by progressively acquired somatic mutations in hematopoietic stem cells (HSCs). Hypomethylating agents (HMAs) can modify the clinical course of MDS and CMML. Clinical improvement does not require eradication of mutated cells and may be related to improved differentiation capacity of mutated HSCs. However, in patients with established disease it is unclear whether (1) HSCs with multiple mutations progress through differentiation with comparable frequency to their less mutated counterparts or (2) improvements in peripheral blood counts following HMA therapy are driven by residual wild-type HSCs or by clones with particular combinations of mutations. To address these questions, the somatic mutations of individual stem cells, progenitors (common myeloid progenitors, granulocyte monocyte progenitors, and megakaryocyte erythroid progenitors), and matched circulating hematopoietic cells (monocytes, neutrophils, and naïve B cells) in MDS and CMML were characterized via high-throughput single-cell genotyping, followed by bulk analysis in immature and mature cells before and after AZA treatment. The mutational burden was similar throughout differentiation, with even the most mutated stem and progenitor clones maintaining their capacity to differentiate to mature cell types in vivo. Increased contributions from productive mutant progenitors appear to underlie improved hematopoiesis in MDS following HMA therapy.

Identification of extremely GC-rich micro RNAs for RT-qPCR data normalization in human plasma

We aimed at extending the repertoire of high-quality miRNA normalizers for reverse transcription-quantitative PCR (RT-qPCR) of human plasma with special emphasis on the extremely guanine-cytosine-rich portion of the miRNome. For high-throughput selection of stable candidates, microarray technology was preferred over small-RNA sequencing (sRNA-seq) since the latter underrepresented miRNAs with a guanine-cytosine (GC) content of at least 75% (p = 0.0002, n = 2). miRNA abundances measured on the microarray were ranked for consistency and uniformity using nine normalization approaches. The eleven most stable sequences included miRNAs of moderate, but also extreme GC content (45%-65%: miR-320d, miR-425-5p, miR-185-5p, miR-486-5p; 80%-95%: miR-1915-3p, miR-3656-5p, miR-3665-5p, miR-3960-5p, miR-4488-5p, miR-4497 and miR-4787-5p). In contrast, the seven extremely GC-rich miRNAs were not found in the two plasma miRNomes screened by sRNA-seq. Stem-loop RT-qPCR was employed for stability verification in 32 plasma samples of healthy male Caucasians (age range: 18-55 years). In general, inter-individual variance of miRNA abundance was low or very low as indicated by coefficient of variation (CV) values of 0.6%-8.2%. miR-3665 and miR-1915-3p outperformed in this analysis (CVs: 0.6 and 2.4%, respectively). The eight most stable sequences included four extremely GC-rich miRNAs (miR-1915-3p, miR-3665, miR-4787-5p and miR-4497). The best-performing duo normalization factor (NF) for the condition of human plasma, miR-320d and miR-4787-5p, also included a GC-extreme miRNA. In summary, the identification of extremely guanine-cytosine-rich plasma normalizers will help to increase accuracy of PCR-based miRNA quantification, thus raise the potential that miRNAs become markers for psychological stress reactions or early and precise diagnosis of clinical phenotypes. The novel miRNAs might also be useful for orthologous contexts considering their conservation in related animal genomes.

OUP accepted manuscript

STUDY QUESTION

Is the composition of microRNAs (miRNAs) in uterine fluid (UF) of women with recurrent implantation failure (RIF) different from that of healthy fertile women?

SUMMARY ANSWER

The composition of miRNAs in UF of women with RIF is different from that of healthy fertile women and the dysregulated miRNAs are associated with impaired endometrial receptivity and embryo implantation.

WHAT IS KNOWN ALREADY

It has previously been demonstrated that the miRNAs secreted from endometrial cells into the UF contribute to the achievement of endometrial receptivity. Endometrial miRNAs are dysregulated in women with RIF.

STUDY DESIGN, SIZE, DURATION

In this descriptive laboratory case–control study, miRNA abundancy was compared between UF collected during implantation phase from healthy fertile women (n = 17) and women with RIF (n = 34), which was defined as three failed IVF cycles with high-quality embryos.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Recruitment of study subjects and sampling of UF were performed at two university clinics in Stockholm, Sweden and Tartu, Estonia. The study participants monitored their menstrual cycles using an LH test kit. The UF samples were collected on Day LH + 7–9 by flushing with saline. Samples were processed for small RNA sequencing and mapped for miRNAs. The differential abundance of miRNAs in UF was compared between the two groups using differential expression analysis (DESeq2). Further downstream analyses, including miRNA target gene prediction (miRTarBase), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis (g:Profiler) and external validation using relevant published data, were performed on the dysregulated miRNAs. Two miRNAs were technically validated with quantitative real-time PCR (RT-PCR).

MAIN RESULTS AND THE ROLE OF CHANCE

After processing of the sequencing data, there were 15 samples in the healthy fertile group and 33 samples in the RIF group. We found 61 differentially abundant UF miRNAs (34 upregulated and 27 downregulated) in RIF compared to healthy women with a false discovery rate of <0.05 and a fold change (FC) of ≤−2 or ≥2. When analyzed with published literature, we found that several of the differentially abundant miRNAs are expressed in endometrial epithelial cells and have been reported in endometrial extracellular vesicles and in association with endometrial receptivity and RIF. Their predicted target genes were further expressed both in the trophectodermal cells of blastocyst-stage embryos and endometrial mid-secretory epithelial cells, as assessed by publicly available single-cell transcriptome-sequencing studies. Pathway analysis further revealed that 25 pathways, having key roles in endometrial receptivity and implantation, were significantly enriched. Hsa-miR-486-5p (FC −20.32; P-value = 0.004) and hsa-miR-92b-3p (FC −9.72; P-value = 0.004) were successfully technically validated with RT-PCR.

LARGE SCALE DATA

The data are available in Gene Expression Omnibus (GEO) at https://www.ncbi.nlm.nih.gov/geo/ with GEO accession number: GSE173289.

LIMITATIONS, REASONS FOR CAUTION

This is a descriptive study with a limited number of study participants. Moreover, the identified differentially abundant miRNAs should be validated in a larger study cohort, and the predicted miRNA target genes and enriched pathways in RIF need to be confirmed and further explored in vitro.

WIDER IMPLICATIONS OF THE FINDINGS

RIF is a major challenge in the current IVF setting with no diagnostic markers nor effective treatment options at hand. For the first time, total miRNAs have been extensively mapped in receptive phase UF of both healthy women with proven fertility and women diagnosed with RIF. Our observations shed further light on the molecular mechanisms behind RIF, with possible implications in future biomarker and clinical treatment studies.

STUDY FUNDING/COMPETING INTEREST(S)

This work was financially supported by the Swedish Research Council (2017-00932), a joint grant from Region Stockholm and Karolinska Institutet (ALF Medicine 2020, FoUI-954072), Estonian Research Council (PRG1076), Horizon 2020 innovation (ERIN, EU952516) and European Commission and Enterprise Estonia (EU48695). The authors have no competing interests to declare for the current study.

Trophoblast derived extracellular vesicles specifically alter the transcriptome of endometrial cells and may constitute a critical component of embryo-maternal communication

BACKGROUND

The period of time when the embryo and the endometrium undergo significant morphological alterations to facilitate a successful implantation-known as "window of implantation"-is a critical moment in human reproduction. Embryo and the endometrium communicate extensively during this period, and lipid bilayer bound nanoscale extracellular vesicles (EVs) are purported to be integral to this communication.

METHODS

To investigate the nature of the EV-mediated embryo-maternal communication, we have supplemented trophoblast analogue spheroid (JAr) derived EVs to an endometrial analogue (RL 95-2) cell layer and characterized the transcriptomic alterations using RNA sequencing. EVs derived from non-trophoblast cells (HEK293) were used as a negative control. The cargo of the EVs were also investigated through mRNA and miRNA sequencing.

RESULTS

Trophoblast spheroid derived EVs induced drastic transcriptomic alterations in the endometrial cells while the non-trophoblast cell derived EVs failed to induce such changes demonstrating functional specificity in terms of EV origin. Through gene set enrichment analysis (GSEA), we found that the response in endometrial cells was focused on extracellular matrix remodelling and G protein-coupled receptors' signalling, both of which are of known functional relevance to endometrial receptivity. Approximately 9% of genes downregulated in endometrial cells were high-confidence predicted targets of miRNAs detected exclusively in trophoblast analogue-derived EVs, suggesting that only a small proportion of reduced expression in endometrial cells can be attributed directly to gene silencing by miRNAs carried as cargo in the EVs.

CONCLUSION

Our study reveals that trophoblast derived EVs have the ability to modify the endometrial gene expression, potentially with functional importance for embryo-maternal communication during implantation, although the exact underlying signalling mechanisms remain to be elucidated.

In Vivo Effects of Mesenchymal Stromal Cells in Two Patients With Severe Acute Respiratory Distress Syndrome

Mesenchymal stromal cells (MSCs) have been investigated as a treatment for various inflammatory diseases because of their immunomodulatory and reparative properties. However, many basic questions concerning their mechanisms of action after systemic infusion remain unanswered. We performed a detailed analysis of the immunomodulatory properties and proteomic profile of MSCs systemically administered to two patients with severe refractory acute respiratory distress syndrome (ARDS) on a compassionate use basis and attempted to correlate these with in vivo anti-inflammatory actions. Both patients received 2×10(6) cells per kilogram, and each subsequently improved with resolution of respiratory, hemodynamic, and multiorgan failure. In parallel, a decrease was seen in multiple pulmonary and systemic markers of inflammation, including epithelial apoptosis, alveolar-capillary fluid leakage, and proinflammatory cytokines, microRNAs, and chemokines. In vitro studies of the MSCs demonstrated a broad anti-inflammatory capacity, including suppression of T-cell responses and induction of regulatory phenotypes in T cells, monocytes, and neutrophils. Some of these in vitro potency assessments correlated with, and were relevant to, the observed in vivo actions. These experiences highlight both the mechanistic information that can be gained from clinical experience and the value of correlating in vitro potency assessments with clinical effects. The findings also suggest, but do not prove, a beneficial effect of lung protective strategies using adoptively transferred MSCs in ARDS. Appropriate randomized clinical trials are required to further assess any potential clinical efficacy and investigate the effects on in vivo inflammation.

SIGNIFICANCE

This article describes the cases of two patients with severe refractory adult respiratory syndrome (ARDS) who failed to improve after both standard life support measures, including mechanical ventilation, and additional measures, including extracorporeal ventilation (i.e., in a heart-lung machine). Unlike acute forms of ARDS (such in the current NIH-sponsored study of mesenchymal stromal cells in ARDS), recovery does not generally occur in such patients.

Full-length RNA-seq from single cells using Smart-seq2

Emerging methods for the accurate quantification of gene expression in individual cells hold promise for revealing the extent, function and origins of cell-to-cell variability. Different high-throughput methods for single-cell RNA-seq have been introduced that vary in coverage, sensitivity and multiplexing ability. We recently introduced Smart-seq for transcriptome analysis from single cells, and we subsequently optimized the method for improved sensitivity, accuracy and full-length coverage across transcripts. Here we present a detailed protocol for Smart-seq2 that allows the generation of full-length cDNA and sequencing libraries by using standard reagents. The entire protocol takes ∼2 d from cell picking to having a final library ready for sequencing; sequencing will require an additional 1-3 d depending on the strategy and sequencer. The current limitations are the lack of strand specificity and the inability to detect nonpolyadenylated (polyA(-)) RNA.

Smart-seq2 for sensitive full-length transcriptome profiling in single cells

Single-cell gene expression analyses hold promise for characterizing cellular heterogeneity, but current methods compromise on either the coverage, the sensitivity or the throughput. Here, we introduce Smart-seq2 with improved reverse transcription, template switching and preamplification to increase both yield and length of cDNA libraries generated from individual cells. Smart-seq2 transcriptome libraries have improved detection, coverage, bias and accuracy compared to Smart-seq libraries and are generated with off-the-shelf reagents at lower cost.

Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells

In the last decade, genome-wide transcriptome analyses have been routinely used to monitor tissue-, disease- and cell type-specific gene expression, but it has been technically challenging to generate expression profiles from single cells. Here we describe a novel and robust mRNA-Seq protocol (Smart-Seq) that is applicable down to single cell levels. Compared with existing methods, Smart-Seq has improved read coverage across transcripts, which significantly enhances detailed analyses of alternative transcript isoforms and identification of SNPs. We have determined the sensitivity and quantitative accuracy of Smart-Seq for single-cell transcriptomics by evaluating it on total RNA dilution series. Applying Smart-Seq to circulating tumor cells from melanomas, we identified distinct gene expression patterns, including new candidate biomarkers for melanoma circulating tumor cells. Importantly, our protocol can easily be utilized for addressing fundamental biological problems requiring genome-wide transcriptome profiling in rare cells.

Analysis of microRNA signatures using size-coded ligation-mediated PCR

The expression pattern and regulatory functions of microRNAs (miRNAs) are intensively investigated in various tissues, cell types and disorders. Differential miRNA expression signatures have been revealed in healthy and unhealthy tissues using high-throughput profiling methods. For further analyses of miRNA signatures in biological samples, we describe here a simple and efficient method to detect multiple miRNAs simultaneously in total RNA. The size-coded ligation-mediated polymerase chain reaction (SL-PCR) method is based on size-coded DNA probe hybridization in solution, followed-by ligation, PCR amplification and gel fractionation. The new method shows quantitative and specific detection of miRNAs. We profiled miRNAs of the let-7 family in a number of organisms, tissues and cell types and the results correspond with their incidence in the genome and reported expression levels. Finally, SL-PCR detected let-7 expression changes in human embryonic stem cells as they differentiate to neuron and also in young and aged mice brain and bone marrow. We conclude that the method can efficiently reveal miRNA signatures in a range of biological samples.

Specific ligation to double-stranded RNA for analysis of cellular RNA::RNA interactions

Double-stranded RNA (dsRNA) is formed in cells as intra- and intermolecular RNA interactions and is involved in a range of biological processes including RNA metabolism, RNA interference and translation control mediated by natural antisense RNA and microRNA. Despite this breadth of activities, few molecular tools are available to analyse dsRNA as native hybrids. We describe a two-step ligation method for enzymatic joining of dsRNA adaptors to any dsRNA molecule in its duplex form without a need for prior sequence or termini information. The method is specific for dsRNA and can ligate various adaptors to label, map or amplify dsRNA sequences. When combined with reverse transcription–polymerase chain reaction, the method is sensitive and can detect low nanomolar concentrations of dsRNA in total RNA. As examples, we mapped dsRNA/single-stranded RNA junctions within Escherichia coli hok mRNA and the human immunodeficiency virus TAR element using RNA from bacteria and mammalian cells.

Antisense PNA accumulates in Escherichia coli and mediates a long post-antibiotic effect

Antisense agents that target growth-essential genes display surprisingly potent bactericidal properties. In particular, peptide nucleic acid (PNA) and phosphorodiamidate morpholino oligomers linked to cationic carrier peptides are effective in time kill assays and as inhibitors of bacterial peritonitis in mice. It is unclear how these relatively large antimicrobials overcome stringent bacterial barriers and mediate killing. Here we determined the transit kinetics of peptide-PNAs and observed an accumulation of cell-associated PNA in Escherichia coli and slow efflux. An inhibitor of drug efflux pumps did not alter peptide-PNA potency, indicating a lack of active efflux from cells. Consistent with cell retention, the post-antibiotic effect (PAE) of the anti-acyl carrier protein (acpP) peptide-PNA was greater than 11 hours. Bacterial cell accumulation and a long PAE are properties of significant interest for antimicrobial development.

Competitive inhibition of natural antisense Sok-RNA interactions activates Hok-mediated cell killing in Escherichia coli

Short regulatory RNAs are widespread in bacteria, and many function through antisense recognition of mRNA. Among the best studied antisense transcripts are RNA antitoxins that repress toxin mRNA translation. The hok/sok locus of plasmid R1 from Escherichia coli is an established model for RNA antitoxin action. Base-pairing between hok mRNA and Sok-antisense-RNA increases plasmid maintenance through post-segregational-killing of plasmid-free progeny cells. To test the model and the idea that sequestration of Sok-RNA activity could provide a novel antimicrobial strategy, we designed anti Sok peptide nucleic acid (PNA) oligomers that, according to the model, would act as competitive inhibitors of hok mRNA::Sok-RNA interactions. In hok/sok-carrying cells, anti Sok PNAs were more bactericidal than rifampicin. Also, anti Sok PNAs induced ghost cell morphology and an accumulation of mature hok mRNA, consistent with cell killing through synthesis of Hok protein. The results support the sense/antisense model for hok mRNA repression by Sok-RNA and demonstrate that antisense agents can be used to out-compete RNA::RNA interactions in bacteria. Finally, BLAST analyses of approximately 200 prokaryotic genomes revealed that many enteric bacteria have multiple hok/sok homologous and analogous RNA-regulated toxin-antitoxin loci. Therefore, it is possible to activate suicide in bacteria by targeting antitoxins.