Effects of RNA integrity on transcript quantification by total RNA sequencing of clinically collected human placental samples

Revealing the molecular landscape of human placenta: a systematic review and meta-analysis of single-cell RNA sequencing studies

BACKGROUND

With increasing significance of developmental programming effects associated with placental dysfunction, more investigations are devoted to improving the characterization and understanding of placental signatures in health and disease. The placenta is a transitory but dynamic organ adapting to the shifting demands of fetal development and available resources of the maternal supply throughout pregnancy. Trophoblasts (cytotrophoblasts, syncytiotrophoblasts, and extravillous trophoblasts) are placental-specific cell types responsible for the main placental exchanges and adaptations. Transcriptomic studies with single-cell resolution have led to advances in understanding the placenta's role in health and disease. These studies, however, often show discrepancies in characterization of the different placental cell types.

OBJECTIVE AND RATIONALE

We aim to review the knowledge regarding placental structure and function gained from the use of single-cell RNA sequencing (scRNAseq), followed by comparing cell-type-specific genes, highlighting their similarities and differences. Moreover, we intend to identify consensus marker genes for the various trophoblast cell types across studies. Finally, we will discuss the contributions and potential applications of scRNAseq in studying pregnancy-related diseases.

SEARCH METHODS

We conducted a comprehensive systematic literature review to identify different cell types and their functions at the human maternal-fetal interface, focusing on all original scRNAseq studies on placentas published before March 2023 and published reviews (total of 28 studies identified) using PubMed search. Our approach involved curating cell types and subtypes that had previously been defined using scRNAseq and comparing the genes used as markers or identified as potential new markers. Next, we reanalyzed expression matrices from the six available scRNAseq raw datasets with cell annotations (four from first trimester and two at term), using Wilcoxon rank-sum tests to compare gene expression among studies and annotate trophoblast cell markers in both first trimester and term placentas. Furthermore, we integrated scRNAseq raw data available from 18 healthy first trimester and nine term placentas, and performed clustering and differential gene expression analysis. We further compared markers obtained with the analysis of annotated and raw datasets with the literature to obtain a common signature gene list for major placental cell types.

OUTCOMES

Variations in the sampling site, gestational age, fetal sex, and subsequent sequencing and analysis methods were observed between the studies. Although their proportions varied, the three trophoblast types were consistently identified across all scRNAseq studies, unlike other non-trophoblast cell types. Notably, no marker genes were shared by all studies for any of the investigated cell types. Moreover, most of the newly defined markers in one study were not observed in other studies. These discrepancies were confirmed by our analysis on trophoblast cell types, where hundreds of potential marker genes were identified in each study but with little overlap across studies. From 35 461 and 23 378 cells of high quality in the first trimester and term placentas, respectively, we obtained major placental cell types, including perivascular cells that previously had not been identified in the first trimester. Importantly, our meta-analysis provides marker genes for major placental cell types based on our extensive curation.

WIDER IMPLICATIONS

This review and meta-analysis emphasizes the need for establishing a consensus for annotating placental cell types from scRNAseq data. The marker genes identified here can be deployed for defining human placental cell types, thereby facilitating and improving the reproducibility of trophoblast cell annotation.

Measuring Up: A Comparison of TapeStation 4200 and Bioanalyzer 2100 as Measurement Tools for RNA Quality in Postmortem Human Brain Samples

The determination of RNA integrity is a critical quality assessment tool for gene expression studies where the experiment's success is highly dependent on the sample quality. Since its introduction in 1999, the gold standard in the scientific community has been the Agilent 2100 Bioanalyzer's RNA integrity number (RIN), which uses a 1-10 value system, from 1 being the most degraded, to 10 being the most intact. In 2015, Agilent launched 4200 TapeStation's RIN equivalent, and reported a strong correlation of r2 of 0.936 and a median error < ±0.4 RIN units. To evaluate this claim, we compared the Agilent 4200 TapeStation's RIN equivalent (RINe) and DV200 to the Agilent 2100 Bioanalyzer's RIN for 183 parallel RNA samples. In our study, using RNA from a total of 183 human postmortem brain samples, we found that the RIN and RINe values only weakly correlate, with an r2 of 0.393 and an average difference of 3.2 RIN units. DV200 also only weakly correlated with RIN (r2 of 0.182) and RINe (r2 of 0.347). Finally, when applying a cut-off value of 6.5 for both metrics, we found that 95.6% of samples passed with RIN, while only 23.5% passed with RINe. Our results suggest that even though RIN (Bioanalyzer) and RINe (TapeStation) use the same 1-10 value system, they should not be used interchangeably, and cut-off values should be calculated independently.

The rat Achilles and patellar tendons have similar increases in mechanical properties but become transcriptionally divergent during postnatal development

The molecular events that drive post-natal tendon development are poorly characterized. In this study, we profiled morphological, mechanical, and transcriptional changes in the rat Achilles and patellar tendon before walking (P7), shortly after onset of walking (P14), and at motor maturity (P28). The Achilles and patellar tendons increased collagen content and mechanical strength similarly throughout post-natal development. However, at P28 the patellar tendon tended to display a higher maximal tensile load (MTL) (P = 0.0524) than the Achilles tendon, but a similar ultimate tensile strength (UTS; load relative to cross-sectional area) probably due to its increased cross-sectional area during development. The tendons started transcriptionally similar, with overlapping PCA clusters at P7 and P14, before becoming transcriptionally distinct at P28. In both tendons, there was an increase in extracellular matrix (ECM) gene expression and a concomitant decrease in cell cycle and mitochondrial gene expression. The transcriptional divergence at P28 suggested that STAT signalling was lower in the patellar tendon where MTL increased the most. Treating engineered human ligaments with the STAT inhibitor itacitinib increased collagen content and MTL. Our results suggest that during post-natal development, cellular resources are initially allocated towards cell proliferation before shifting towards extracellular matrix development following the onset of mechanical load and provide potential targets for improving tendon function. KEY POINTS: Little is known about mechanisms of post-natal tendon growth. We characterized morphological, mechanical, and transcriptional changes that occur before (P7), and early (P14) and late after (P28) rats begin to walk. From P7 to P28, the Achilles tendon increased in length, whereas the patellar tendon increased in cross-sectional area. Mechanical and material properties of the Achilles and patellar tendon increased from P7 to P28. From P7 to P28, the Achilles and patellar tendons increased expression of ECM genes and decreased mitochondrial and cell cycle gene expression. Ribosomal gene expression also significantly decreased in the Achilles and tended to decrease in the patellar tendon. At P28, STAT1 signalling tended to be lower in the patellar tendon which had grown by increasing cross-sectional area and inhibiting STAT activation in vitro improved mechanical properties in engineered human ligaments.

High purity and integrity RNA from human cell lines stored in liquid nitrogen for over 20 years

Biobanks are a key resource for obtaining human cell lines for biomedical research, including for drug development projects. Such projects often include comparative RNA-sequencing of large panels of human cell lines from individuals affected by certain disorders and healthy controls, or from individuals with different drug response phenotypes. RNA extractions are typically done from growing cell cultures, a process that may take several weeks. However, maintaining large numbers of cell lines in parallel increases the project workload. Here, we show that extracting RNAs directly from frozen vials of human cell lines stored for over 20 years in a liquid nitrogen freezer yields RNAs with the high purity and integrity parameters that conform to those required for optimal RNA-sequencing and are closely similar to those obtained for RNAs extracted from growing human cell lines.

Pervasive effects of RNA degradation on Nanopore direct RNA sequencing

Oxford Nanopore direct RNA sequencing (DRS) is capable of sequencing complete RNA molecules and accurately measuring gene and isoform expression. However, as DRS is designed to profile intact RNA, expression quantification may be more heavily dependent upon RNA integrity than alternative RNA sequencing methodologies. It is currently unclear how RNA degradation impacts DRS or whether it can be corrected for. To assess the impact of RNA integrity on DRS, we performed a degradation time series using SH-SY5Y neuroblastoma cells. Our results demonstrate that degradation is a significant and pervasive factor that can bias DRS measurements, including a reduction in library complexity resulting in an overrepresentation of short genes and isoforms. Degradation also biases differential expression analyses; however, we find that explicit correction can almost fully recover meaningful biological signal. In addition, DRS provided less biased profiling of partially degraded samples than Nanopore PCR-cDNA sequencing. Overall, we find that samples with RNA integrity number (RIN) > 9.5 can be treated as undegraded and samples with RIN > 7 can be utilized for DRS with appropriate correction. These results establish the suitability of DRS for a wide range of samples, including partially degraded in vivo clinical and post-mortem samples, while limiting the confounding effect of degradation on expression quantification.

Nuclei on the Rise: When Nuclei-Based Methods Meet Next-Generation Sequencing

In the last decade, we have witnessed an upsurge in nuclei-based studies, particularly coupled with next-generation sequencing. Such studies aim at understanding the molecular states that exist in heterogeneous cell populations by applying increasingly more affordable sequencing approaches, in addition to optimized methodologies developed to isolate and select nuclei. Although these powerful new methods promise unprecedented insights, it is important to understand and critically consider the associated challenges. Here, we provide a comprehensive overview of the rise of nuclei-based studies and elaborate on their advantages and disadvantages, with a specific focus on their utility for transcriptomic sequencing analyses. Improved designs and appropriate use of the various experimental strategies will result in acquiring biologically accurate and meaningful information.

Optimization of RNA extraction for bacterial whole transcriptome studies of low-biomass samples

We developed a procedure for extracting maximal amounts of high-quality RNA from low-biomass producing (autotrophic) bacteria for experiments where sample volume is limited. Large amounts of high-quality RNA for downstream analyses cannot be obtained using larger quantities of culture volume. The performance of standard commercial silica-column based kit protocols and these procedures amended by ultrasonication or enzymatic lysis were assessed. The ammonium-oxidizing Nitrosomonas europaea and nitrite-oxidizing Nitrobacter winogradskyi were used as model organisms for optimization of the RNA isolation protocol. Enzymatic lysis through lysozyme digestion generated high-quality, high-yield RNA samples. Subsequent RNA-seq analysis resulted in qualitative data for both strains. The RNA extraction procedure is suitable for experiments with volume and/or biomass limitations, e.g., as encountered during space flight experiments. Furthermore, it will also result in higher RNA yields for whole transcriptome experiments where sample volume and/or biomass was increased to compensate the low-biomass characteristic of autotrophs.

•

Lysis through ultrasonication results in a high RNA yield but in a low RNA quality

•

Enzymatic lysis using lysozyme provides high-yield, high-quality RNA samples

•

RNA from low-biomass bacterial samples is suitable for downstream RNA-seq

Laser Capture Microdissection Optimization for High-Quality RNA in Mouse Brain Tissue

Laser Capture Microdissection (LCM) is a method that allows one to select and dissect well-defined structures, specific cell subpopulations, or even single cells from different types of tissue for subsequent extraction of DNA, RNA, or proteins. Its precision allows the dissection of specific groups of cells, avoiding unwanted cells. However, despite its efficiency, several steps can affect the sample RNA integrity. RNA instability represents a challenge in the LCM method, and low RNA integrity can introduce biases, as different transcripts often have different degradation rates. Here we describe an optimized protocol to provide good-concentration and high-quality RNA from specific structures: dentate gyrus and CA1 in the hippocampus, basolateral amygdala, and anterior cingulate cortex of mouse brain tissue. However, the protocol is applicable to other areas of interest. © 2022 Wiley Periodicals LLC. Basic Protocol: Laser capture microdissection of mouse brain tissue.

The transcript integrity index (TII) provides a standard measure of sperm RNA

Standardizing RNA quality is key to interpreting RNA-seq data as a compromised sample can mask the underlying biology. The challenge remains when evaluating RNA quality in samples with high RNA fragmentation. For example, programmed fragmentation and cytoplasmic expulsion, integral to sperm maturation, is a prime example of the complexities of interpreting RNA-seq data, given that fragmentation can be random and\or targeted. To meet this challenge, we developed an algorithm that accurately measures RNA quality in samples with high fragmentation, such as spermatozoa. The integrity of 1,000 previously identified abundant sperm transcripts were independently visualized and evaluated using the Transcript Integrity Index (TII) algorithm to identify intact transcripts. Full-length transcripts from visual and the TII algorithm were evaluated for testis preference in humans using the GTEx tissues database. Samples were then filtered by the Interquartile Range (IQR), identifying those in which the greatest number of transcripts failed to pass the visual or TII thresholds. Transcript lists were overlapped, forming the set of intact transcripts used as TII standards. Each sample was re-evaluated as a function of this TII set of intact transcripts, with poor quality samples identified as those failing in the largest number of transcripts. While ontologically enriched in roles related to spermatogenesis and/or fertilization, samples did not segregate based on birth outcome. The TII algorithm proved an effective means to identify samples of similar quality from sperm, a cell type enriched in biologically fragmented RNAs. The algorithm should facilitate other studies using samples compromised by high levels of RNA fragmentation, such as Formalin-Fixed Paraffin-Embedded samples. Requisite to assessing male health, TII provides a solution to the long-sought-after standard that identifies samples of similar quality.

Assessment of RNA extraction protocols from cladocerans

The usage of cladocerans as non-model organisms in ecotoxicological and risk assessment studies has intensified in recent years due to their ecological importance in aquatic ecosystems. The molecular assessment such as gene expression analysis has been introduced in ecotoxicological and risk assessment to link the expression of specific genes to a biological process in the cladocerans. The validity and accuracy of gene expression analysis depends on the quantity, quality and integrity of extracted ribonucleic acid (RNA) of the sample. However, the standard methods of RNA extraction from the cladocerans are still lacking. This study evaluates the extraction of RNA from tropical freshwater cladocerans Moina micrura using two methods: the phenol-chloroform extraction method (QIAzol) and a column-based kit (Qiagen Micro Kit). Glycogen was introduced in both approaches to enhance the recovery of extracted RNA and the extracted RNA was characterised using spectrophotometric analysis (NanoDrop), capillary electrophoresis (Bioanalyzer). Then, the extracted RNA was analysed with reverse transcription polymerase chain reaction (RT-PCR) to validate the RNA extraction method towards downstream gene expression analysis. The results indicate that the column-based kit is most suitable for the extraction of RNA from M. micrura, with the quantity (RNA concentration = 26.90 ± 6.89 ng/μl), quality (A260:230 = 1.95 ± 0.15, A280:230 = 1.85 ± 0.09) and integrity (RNA integrity number, RIN = 7.20 ± 0.16). The RT-PCR analysis shows that the method successfully amplified both alpha tubulin and actin gene at 33-35 cycles (i.e. Ct = 32.64 to 33.48). The results demonstrate that the addition of glycogen is only suitable for the phenol-chloroform extraction method. RNA extraction with high and comprehensive quality control assessment will increase the accuracy and reliability of downstream gene expression, thus providing more ecotoxicological data at the molecular biological level on other freshwater zooplankton species.

Dissecting the transcriptome in cardiovascular disease

The human transcriptome comprises a complex network of coding and non-coding RNAs implicated in a myriad of biological functions. Non-coding RNAs exhibit highly organized spatial and temporal expression patterns and are emerging as critical regulators of differentiation, homeostasis, and pathological states, including in the cardiovascular system. This review defines the current knowledge gaps, unmet methodological needs, and describes the challenges in dissecting and understanding the role and regulation of the non-coding transcriptome in cardiovascular disease. These challenges include poor annotation of the non-coding genome, determination of the cellular distribution of transcripts, assessment of the role of RNA processing and identification of cell-type specific changes in cardiovascular physiology and disease. We highlight similarities and differences in the hurdles associated with the analysis of the non-coding and protein-coding transcriptomes. In addition, we discuss how the lack of consensus and absence of standardized methods affect reproducibility of data. These shortcomings should be defeated in order to make significant scientific progress and foster the development of clinically applicable non-coding RNA-based therapeutic strategies to lessen the burden of cardiovascular disease.

Altered human alveolar bone gene expression in type 2 diabetes-A cross-sectional study

OBJECTIVE

The objective of this cross-sectional study is to investigate alveolar bone gene expression in health and diabetes through ribonucleic acid (RNA) sequencing and bioinformatics analysis.

BACKGROUND

It is relatively unknown how type 2 diabetes modulates gene expression in alveolar bone in humans. Clinical concern regarding increased implant failure rate in patients with diabetes has been discussed in the literature. Previous studies in animal models and humans have suggested an imbalance between the genes regulating bone formation with data suggesting bone resorption in diabetes. However, there is lack of data regarding a comprehensive gene expression from human alveolar bone in diabetes.

METHODS

Alveolar bone was collected from healthy and type 2 diabetic subjects undergoing periodontal and implant surgeries. The homogenized RNA sample was then extracted and analyzed for quantity and quality. RNA samples were further purified using ribosomal RNA depletion technique and processed for RNA sequencing and analysis. Expression levels for mRNAs were performed by calculating FPKM ([total_exon_fragments/mapped reads (millions) × exon length (kB)]), and differentially expressed mRNAs were selected with log₂ (fold change) >1 or log₂ (fold change) ≤1 and with a parametric F test comparing nested linear models.

RESULTS

Eighteen bone samples (10 healthy, 8 patients with diabetes) were analyzed for gene expression. The mean age and HbA1c% of healthy versus diabetic subjects were as follows: age (55.3 ± 17.5 vs 63.9 ± 8.7 years) and HbA1c% (5.6 ± 0.29 vs 7.3 ± 2.4), respectively. Sequencing analysis showed that expression of genes that regulate bone turnover like TGFB1, LTBP4, IGF1, BMP2, BMP4, BMP6, SMAD1, RUNX2, MCSF, and THRA was significantly downregulated in diabetes samples compared with healthy controls with overall reduced expression of genes in the bone regulation pathway in patients with diabetes. Bioinformatics analysis for the altered genes highlighted several pathways related to bone homeostasis and inflammation in diabetes. Periodontitis did not affect the gene expression pattern based on diabetes status.

CONCLUSIONS

Altered expression of genes due to downregulation of certain pathways that are involved in bone turnover and inflammation suggests that overall wound healing and bone homeostasis may be compromised in type 2 diabetes.

Optimization of RNA extraction methods from human metabolic tissue samples of the COMET biobank

Constitution of biobank of human tissues requires careful handling and storage of biological material, to guarantee the quality of samples. Tissue preparation is also critical for further applications such as transcriptomic profiling. In this study, our aim was to evaluate the impact of different disruption techniques (FastPrep-24 instrument, GentleMACS dissociator, and syringe/needle) and homogenizing buffers (RLT versus QIAzol) on RNA purity and quality of metabolic tissues (adipose tissues, liver and skeletal muscle) present in the COMET Biobank. For all homogenization methods used and tissue types, the A260/280 ratios reached values ≥ 1.8, which are in the range of what is found in human tissues and cell lines, while the A260/230 ratios were however ≤ 1.8, with the lowest value obtained with GentleMACS Dissociator. In addition, GentleMACS Dissociator combined with QIAzol reagent gave the highest RIN value and 28S/18S ratio for all tissues tested, except for muscle. Performing RT-qPCR, Ct values for different housekeeping genes can be influenced by extraction methods and RNA quality of samples. In conclusion, we have demonstrated that different disruption techniques and homogenizing buffers impact the purity and some quality markers of RNA, and can also impact quantification of mRNAs by RT-qPCR in human metabolic tissues.

A new pipeline for clinico-pathological and molecular placental research utilizing FFPE tissues

The placenta is at the core of many pregnancy pathologies, but we have limited knowledge about placental function because of two key research barriers: 1) lack of guidelines for sample collection and pathologic diagnosis; and 2) limited tools are available for molecular analysis of stored placental samples. We aimed to create a searchable, population-based placental database of pathologic diagnoses, and to validate molecular methods for gene expression studies of matching formalin fixed paraffin embedded (FFPE) placental blocks. Our database has over 1000 pregnancies coded for clinical diagnosis with corresponding FFPE blocks that are available for gene expression studies. RNA harvested from FFPE tissues is of sufficient quality for downstream applications. We successfully used this pipeline to identify FFPE placenta from term and preterm pregnancies, and compared their gene expression. The establishment of this platform, which links clinicopathological data and molecular gene expression, will increase our understanding of obstetrical diseases.

Coordinated Expressional Landscape of the Human Placental miRNome and Transcriptome

Placenta is a unique organ that serves its own function, and contributes to maternal gestational adaptation and fetal development. Coordination of its transcriptome to satisfy all the maternal-fetal needs across gestation is not fully understood. MicroRNAs are powerful transcriptome modulators capable to adjust rapidly the expression level and dynamics of large gene sets. This MiR-Seq based study presents a multi-omics investigation of the human placental miRNome and its synergy with the transcriptome. The analysis included 52 placentas representing three trimesters of normal pregnancy, and term cases of late-onset preeclampsia (LO-PE), gestational diabetes and affected fetal growth. Gestational-age dependent differential expression (FDR < 0.05) was detected for 319 of 417 tested miRNAs (76.5%). A shared list of target genes of dynamic miRNAs suggested their coordinated action. The most abundant miR-143-3p revealed as a marker for pregnancy progression. The data suggested critical, but distinct roles of placenta-specific imprinted C19MC and C14MC miRNA clusters. Paternally encoded primate-specific C19MC was highly transcribed during first trimester, potentially fine-tuning the early placental transcriptome in dosage-sensitive manner. Maternally encoded eutherian C14MC showed high expression until term, underlining its key contribution across gestation. A major shift in placental miRNome (16% miRNAs) was observed in LO-PE, but not in other term pregnancy complications. Notably, 13/38 upregulated miRNAs were transcribed from C19MC and only one from C14MC, whereas 11/28 downregulated miRNAs represented C14MC and none C19MC. miR-210-3p, miR-512-5p, miR-32-5p, miR-19a-3p, miR-590-3p, miR-379-5p were differentially expressed in LO-PE and cases of small-for-gestational-age newborns, supporting a shared etiology. Expression correlation analysis with the RNA-Seq data (16,567 genes) of the same samples clustered PE-linked miRNAs into five groups. Large notable clusters of miRNA–gene pairs showing directly and inversely correlated expression dynamics suggested potential functional relationships in both scenarios. The first genome-wide study of placental miR-eQTLs identified 66 placental SNVs associated with the expression of neighboring miRNAs, including PE-linked miRNAs miR-30a-5p, miR-210-3p, miR-490-3p and miR-518-5p. This study provided a rich catalog of miRNAs for further in-depth investigations of their individual and joint effect on placental transcriptome. Several highlighted miRNAs may serve as potential biomarkers for pregnancy monitoring and targets to prevent or treat gestational disorders.

Ischemia-induced upregulation of autophagy preludes dysfunctional lysosomal storage and associated synaptic impairments in neurons

Macroautophagy/autophagy is vital for neuronal homeostasis and functions. Accumulating evidence suggest that autophagy is impaired during cerebral ischemia, contributing to neuronal dysfunction and neurodegeneration. However, the outcomes after transient modification in autophagy machinery are not fully understood. This study investigated the effects of ischemic stress on autophagy and synaptic structures using a rat model of oxygen-glucose deprivation (OGD) in hippocampal neurons and a mouse model of middle cerebral artery occlusion (MCAO). Upon acute ischemia, an initial autophagy modification occurred in an upregulation manner. Following, the number of lysosomes increased, as well as lysosomal volume, indicating dysfunctional lysosomal storage. These changes were prevented by inhibiting autophagy via 3-methyladenine (3-MA) treatment or ATG7 (autophagy related 7) knockdown, or were mimicked by rapamycin (RAPA), a known activator of autophagy. This suggests that dysfunctional lysosomal storage is associated with the early burst of autophagy. Dysfunctional lysosomal storage contributed to autophagy dysfunction because the basal level of MTOR-dependent lysosomal biogenesis in the reperfusion was not sufficient to clear undegraded cargoes after transient autophagy upregulation. Further investigation revealed that impairment of synaptic ultra-structures, accompanied by dysfunctional lysosomal storage, may result from a failure in dynamic turnover of synaptic proteins. This indicates a vital role of autophagy-lysosomal machinery in the maintenance of synaptic structures. This study supports previous evidence that dysfunctional lysosomal storage may occur following the upregulation of autophagy in neurons. Appropriate autophagosome-lysosomal functioning is vital for maintenance of neuronal synaptic function and impacts more than the few known synaptic proteins.Abbreviations: 3-MA: 3-methyladenine; ACTB: actin beta; AD: Alzheimer disease; ALR: autophagic lysosome reformation; ATG7: autophagy related 7; CTSB: cathepsin B; CTSD: cathepsin D; DAPI: 4',6-diamidino-2-phenylindole; DEGs: differentially expressed genes; DMEM: Dulbecco's modified Eagle's medium; DMSO: dimethyl sulfoxide; GO: Gene Ontology; HBSS: Hanks' balanced salt solution; HPCA: hippocalcin; i.c.v: intracerebroventricular; KEGG: kyoto encyclopedia of genes and genomes; LAMP1: lysosomal-associated membrane protein 1; MAP1LC3B/LC3: microtubule-associated protein 1 light chain 3 beta; LSDs: lysosomal storage disorders; MAP2: microtubule-associated protein 2; MCAO: middle cerebral artery occlusion; mCTSB: mature CTSB; mCTSD: mature CTSD; MOI: multiplicity of infection; MTOR: mechanistic target of rapamycin kinase; OGD/R: oxygen-glucose deprivation/reoxygenation; PBS: phosphate-buffered saline; PRKAA/AMPKα: protein kinase AMP-activated catalytic subunit alpha; proCTSD: pro-cathepsin D; RAPA: rapamycin; RNA-seq: RNA sequencing; RPS6KB/p70S6K: ribosomal protein S6 kinase; SDS-PAGE: sodium dodecyl sulfate-polyacrylamide gel electrophoresis; SIM: Structured Illumination Microscopy; SNAP25: synaptosomal-associated protein 25; SQSTM1/p62: sequestosome 1; SYN1: synapsin I; SYT1: synaptotagmin I; TBST: tris-buffered saline Tween-20; TEM: transmission electron microscopy; TFEB: transcription factor EB; tMCAO: transient middle cerebral artery occlusion; TTC: 2,3,5-triphenyltetrazolium chloride; TUBB3: tubulin, beta 3 class III.

Current approaches and developments in transcript profiling of the human placenta

BACKGROUND

The placenta is the active interface between mother and foetus, bearing the molecular marks of rapid development and exposures in utero. The placenta is routinely discarded at delivery, providing a valuable resource to explore maternal-offspring health and disease in pregnancy. Genome-wide profiling of the human placental transcriptome provides an unbiased approach to study normal maternal–placental–foetal physiology and pathologies.

OBJECTIVE AND RATIONALE

To date, many studies have examined the human placental transcriptome, but often within a narrow focus. This review aims to provide a comprehensive overview of human placental transcriptome studies, encompassing those from the cellular to tissue levels and contextualize current findings from a broader perspective. We have consolidated studies into overarching themes, summarized key research findings and addressed important considerations in study design, as a means to promote wider data sharing and support larger meta-analysis of already available data and greater collaboration between researchers in order to fully capitalize on the potential of transcript profiling in future studies.

SEARCH METHODS

The PubMed database, National Center for Biotechnology Information and European Bioinformatics Institute dataset repositories were searched, to identify all relevant human studies using ‘placenta’, ‘decidua’, ‘trophoblast’, ‘transcriptome’, ‘microarray’ and ‘RNA sequencing’ as search terms until May 2019. Additional studies were found from bibliographies of identified studies.

OUTCOMES

The 179 identified studies were classifiable into four broad themes: healthy placental development, pregnancy complications, exposures during pregnancy and in vitro placental cultures. The median sample size was 13 (interquartile range 8–29). Transcriptome studies prior to 2015 were predominantly performed using microarrays, while RNA sequencing became the preferred choice in more recent studies. Development of fluidics technology, combined with RNA sequencing, has enabled transcript profiles to be generated of single cells throughout pregnancy, in contrast to previous studies relying on isolated cells. There are several key study aspects, such as sample selection criteria, sample processing and data analysis methods that may represent pitfalls and limitations, which need to be carefully considered as they influence interpretation of findings and conclusions. Furthermore, several areas of growing importance, such as maternal mental health and maternal obesity are understudied and the profiling of placentas from these conditions should be prioritized.

WIDER IMPLICATIONS

Integrative analysis of placental transcriptomics with other ‘omics’ (methylome, proteome and metabolome) and linkage with future outcomes from longitudinal studies is crucial in enhancing knowledge of healthy placental development and function, and in enabling the underlying causal mechanisms of pregnancy complications to be identified. Such understanding could help in predicting risk of future adversity and in designing interventions that can improve the health outcomes of both mothers and their offspring. Wider collaboration and sharing of placental transcriptome data, overcoming the challenges in obtaining sufficient numbers of quality samples with well-defined clinical characteristics, and dedication of resources to understudied areas of pregnancy will undoubtedly help drive the field forward.

Time-dependent changes in placental mRNA expression after delivery due to delayed specimen collection

PROBLEM

Current scientific guidelines recommend collecting placental specimens within two hours of delivery for gene expression analysis. However, collecting samples in a narrow time window is a challenge in the dynamic and unpredictable clinical setting, so delays in placental specimen collection are possible. The purpose of our analysis was to investigate temporal changes in placental gene expression by longitudinally sampling placentas over a 24 h period.

METHOD OF STUDY

Eight placentas from individuals with uncomplicated, term pregnancies delivered by scheduled cesarean section were collected and sampled following the placental delivery and again at 1, 2, 4, 6, and 24 h post-delivery. At each time point, biopsies of chorionic villous tissue were taken from 3 cotyledons to account for intra-placental heterogeneity. The 3 biopsies from each time point were pooled prior to RNA extraction. Expression of 382 mRNA transcripts was quantified using the NanoString nCounter System. Fold change values were calculated for each time point relative to delivery, and a fold change threshold of 1.25 was used to determine a meaningful change from delivery.

RESULTS

Based on a fold change threshold of 1.25, 84.3% of transcripts were stable for at least 1 h, 80.2% were stable for at least two hours, and 20.6% of transcripts were stable through the collection at 24 h.

CONCLUSION

Our results suggest that for some mRNA transcripts, expression changes as time to sample collection increases. We have developed a Web application to allow investigators to explore transcripts relevant to their research interests and to set appropriate thresholds to aid in determining whether placentas with delayed sample collection can be included in analyses (https://placentaexpression.foundationsofhealth.org/).

INTACT vs. FANS for Cell-Type-Specific Nuclei Sorting: A Comprehensive Qualitative and Quantitative Comparison

Increasing numbers of studies seek to characterize the different cellular sub-populations present in mammalian tissues. The techniques “Isolation of Nuclei Tagged in Specific Cell Types” (INTACT) or “Fluorescence-Activated Nuclei Sorting” (FANS) are frequently used for isolating nuclei of specific cellular subtypes. These nuclei are then used for molecular characterization of the cellular sub-populations. Despite the increasing popularity of both techniques, little is known about their isolation efficiency, advantages, and disadvantages or downstream molecular effects. In our study, we compared the physical and molecular attributes of sfGFP+ nuclei isolated by the two methods—INTACT and FANS—from the neocortices of Arc-CreERT2 × CAG-Sun1/sfGFP animals. We identified differences in efficiency of sfGFP+ nuclei isolation, nuclear size as well as transcriptional (RNA-seq) and chromatin accessibility (ATAC-seq) states. Therefore, our study presents a comprehensive comparison between the two widely used nuclei sorting techniques, identifying the advantages and disadvantages for both INTACT and FANS. Our conclusions are summarized in a table to guide researchers in selecting the most suitable methodology for their individual experimental design.

Detection of H3K4me3 Identifies NeuroHIV Signatures, Genomic Effects of Methamphetamine and Addiction Pathways in Postmortem HIV+ Brain Specimens that Are Not Amenable to Transcriptome Analysis

Human postmortem specimens are extremely valuable resources for investigating translational hypotheses. Tissue repositories collect clinically assessed specimens from people with and without HIV, including age, viral load, treatments, substance use patterns and cognitive functions. One challenge is the limited number of specimens suitable for transcriptional studies, mainly due to poor RNA quality resulting from long postmortem intervals. We hypothesized that epigenomic signatures would be more stable than RNA for assessing global changes associated with outcomes of interest. We found that H3K27Ac or RNA Polymerase (Pol) were not consistently detected by Chromatin Immunoprecipitation (ChIP), while the enhancer H3K4me3 histone modification was abundant and stable up to the 72 h postmortem. We tested our ability to use HeK4me3 in human prefrontal cortex from HIV+ individuals meeting criteria for methamphetamine use disorder or not (Meth +/-) which exhibited poor RNA quality and were not suitable for transcriptional profiling. Systems strategies that are typically used in transcriptional metadata were applied to H3K4me3 peaks revealing consistent genomic activity differences in regions where addiction and neuronal synapses pathway genes are represented, including genes of the dopaminergic system, as well as inflammatory pathways. The resulting comparisons mirrored previously observed effects of Meth on suppressing gene expression and provided insights on neurological processes affected by Meth. The results suggested that H3K4me3 detection in chromatin may reflect transcriptional patterns, thus providing opportunities for analysis of larger numbers of specimens from cases with substance use and neurological deficits. In conclusion, the detection of H3K4me3 in isolated chromatin can be an alternative to transcriptome strategies to increase the power of association using specimens with long postmortem intervals and low RNA quality.

Multiple freeze-thaw cycles lead to a loss of consistency in poly(A)-enriched RNA sequencing

Background

Both RNA-Seq and sample freeze-thaw are ubiquitous. However, knowledge about the impact of freeze-thaw on downstream analyses is limited. The lack of common quality metrics that are sufficiently sensitive to freeze-thaw and RNA degradation, e.g. the RNA Integrity Score, makes such assessments challenging.

Results

Here we quantify the impact of repeated freeze-thaw cycles on the reliability of RNA-Seq by examining poly(A)-enriched and ribosomal RNA depleted RNA-seq from frozen leukocytes drawn from a toddler Autism cohort. To do so, we estimate the relative noise, or percentage of random counts, separating technical replicates. Using this approach we measured noise associated with RIN and freeze-thaw cycles. As expected, RIN does not fully capture sample degradation due to freeze-thaw. We further examined differential expression results and found that three freeze-thaws should extinguish the differential expression reproducibility of similar experiments. Freeze-thaw also resulted in a 3′ shift in the read coverage distribution along the gene body of poly(A)-enriched samples compared to ribosomal RNA depleted samples, suggesting that library preparation may exacerbate freeze-thaw-induced sample degradation.

Conclusion

The use of poly(A)-enrichment for RNA sequencing is pervasive in library preparation of frozen tissue, and thus, it is important during experimental design and data analysis to consider the impact of repeated freeze-thaw cycles on reproducibility.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07381-z.

Low RIN Value for RNA-Seq Library Construction from Long-Term Stored Seeds: A Case Study of Barley Seeds

Seed aging is a complex biological process and its fundamentals and mechanisms have not yet been fully recognized. This is a key issue faced by research teams involved in the collection and storage of plant genetic resources in gene banks every day. Transcriptomic changes associated with seed aging in the dry state have barely been studied. The aim of the study was to develop an efficient protocol for construction of RNA-Seq libraries from long-term stored seeds with very low viability and low RNA integrity number (RIN). Here, barley seeds that have almost completely lost their viability as a result of long-term storage were used. As a control, fully viable seeds obtained in the course of field regeneration were used. The effectiveness of protocols dedicated to RNA samples with high and low RIN values was compared. The experiment concluded that library construction from low viable or long-term stored seeds with degraded RNA (RIN < 3) should be carried out with extraordinary attention due to the possibility of uneven degradation of different RNA fractions.

Characterization of left ventricular myocardial sodium-glucose cotransporter 1 expression in patients with end-stage heart failure

Background

Whereas selective sodium-glucose cotransporter 2 (SGLT2) inhibitors consistently showed cardiovascular protective effects in large outcome trials independent of the presence of type 2 diabetes mellitus (T2DM), the cardiovascular effects of dual SGLT1/2 inhibitors remain to be elucidated. Despite its clinical relevance, data are scarce regarding left ventricular (LV) SGLT1 expression in distinct heart failure (HF) pathologies. We aimed to characterize LV SGLT1 expression in human patients with end-stage HF, in context of the other two major glucose transporters: GLUT1 and GLUT4.

Methods

Control LV samples (Control, n = 9) were harvested from patients with preserved LV systolic function who went through mitral valve replacement. LV samples from HF patients undergoing heart transplantation (n = 71) were obtained according to the following etiological subgroups: hypertrophic cardiomyopathy (HCM, n = 7); idiopathic dilated cardiomyopathy (DCM, n = 12); ischemic heart disease without T2DM (IHD, n = 14), IHD with T2DM (IHD + T2DM, n = 11); and HF patients with cardiac resynchronization therapy (DCM:CRT, n = 9, IHD:CRT, n = 9 and IHD-T2DM:CRT, n = 9). We measured LV SGLT1, GLUT1 and GLUT4 gene expressions with qRT-PCR. The protein expression of SGLT1, and activating phosphorylation of AMP-activated protein kinase (AMPKα) and extracellular signal-regulated kinase 1/2 (ERK1/2) were quantified by western blotting. Immunohistochemical staining of SGLT1 was performed.

Results

Compared with controls, LV SGLT1 mRNA and protein expressions were significantly and comparably upregulated in HF patients with DCM, IHD and IHD + T2DM (all P < 0.05), but not in HCM. LV SGLT1 mRNA and protein expressions positively correlated with LVEDD and negatively correlated with EF (all P < 0.01). Whereas AMPKα phosphorylation was positively associated with SGLT1 protein expression, ERK1/2 phosphorylation showed a negative correlation (both P < 0.01). Immunohistochemical staining revealed that SGLT1 expression was predominantly confined to cardiomyocytes, and not fibrotic tissue. Overall, CRT was associated with reduction of LV SGLT1 expression, especially in patients with DCM.

Conclusions

Myocardial LV SGLT1 is upregulated in patients with HF (except in those with HCM), correlates significantly with parameters of cardiac remodeling (LVEDD) and systolic function (EF), and is downregulated in DCM patients with CRT. The possible role of SGLT1 in LV remodeling needs to be elucidated.

C-allele of rs4769613 Near FLT1 Represents a High-Confidence Placental Risk Factor for Preeclampsia

The variant rs4769613 T/C within the enhancer element near FLT1, an acknowledged gene in preeclampsia, was previously identified as a risk factor for preeclampsia in the genome-wide association study (GWAS) targeting placental genotypes. We aimed to test the robustness of this association in 2 Estonian cohorts. Both placental sample sets HAPPY PREGNANCY (Development of novel non-invasive biomarkers for fertility and healthy pregnancy; preeclampsia, n=44 versus nonpreeclampsia, n=1724) and REPROMETA (REPROgrammed fetal and/or maternal METAbolism; 52/277) exhibited suggestive association between rs4769613[C] variant and preeclampsia (logistic regression adjusted for gestational age and fetal sex, nominal P<0.05). Meta-analysis across 2 samples (96/2001) replicated the genome-wide association study outcome (Bonferroni corrected P=4×10^-3; odds ratio, 1.75 [95% CI, 1.23-2.49]). No association was detected with gestational diabetes mellitus, preterm birth, and newborn parameters. Also, neither maternal nor paternal rs4769613 genotypes predisposed to preeclampsia. The exact role of placental rs4769613 genotype in the preeclampsia pathogenesis is to be clarified as no effect was detected on maternal baseline serum sFlt-1 (soluble fms-related receptor tyrosine kinase 1) levels. However, when placental FLT1 gene expression and maternal serum sFlt-1 measurements were stratified by placental rs4769613 genotypes, significantly higher transcript and biomarker levels were detected in preeclampsia versus nonpreeclampsia cases in the CC- and CT- (Student t test, P≤0.02), but not in the TT-genotype subgroup. We suggest that rs4769613 represents a conditional expression Quantitative Trait Locus, whereby only the enhancer with the C-allele reacts to promote the FLT1 expression in unfavorable placental conditions. The study highlighted that the placental FLT1 rs4769613 C-allele is a preeclampsia-specific risk factor. It may contribute to early identification of high-risk women, for example, when genotyped in the cffDNA available in maternal blood plasma.

A Guide to Understanding "State-of-the-Art" Basic Research Techniques in Anesthesiology

Perioperative medicine is changing from a "protocol-based" approach to a progressively personalized care model. New molecular techniques and comprehensive perioperative medical records allow for detection of patient-specific phenotypes that may better explain, or even predict, a patient's response to perioperative stress and anesthetic care. Basic science technology has significantly evolved in recent years with the advent of powerful approaches that have translational relevance. It is incumbent on us as a primarily clinical specialty to have an in-depth understanding of rapidly evolving underlying basic science techniques to incorporate such approaches into our own research, critically interpret the literature, and improve future anesthesia patient care. This review focuses on 3 important and most likely practice-changing basic science techniques: next-generation sequencing (NGS), clustered regularly interspaced short palindromic repeat (CRISPR) modulations, and inducible pluripotent stem cells (iPSCs). Each technique will be described, potential advantages and limitations discussed, open questions and challenges addressed, and future developments outlined. We hope to provide insight for practicing physicians when confronted with basic science articles and encourage investigators to apply "state-of-the-art" technology to their future experiments.

Congenital cystic adenomatoid malformations of the lung: an epithelial transcriptomic approach

Background

The pathophysiology of congenital cystic adenomatoid malformations (CCAM) of the lung remains poorly understood.

Aim

This study aimed to identify more precisely the molecular mechanisms limited to a compartment of lung tissue, through a transcriptomic analysis of the epithelium of macrocystic forms.

Methods

Tissue fragments displaying CCAM were obtained during planned surgical resections. Epithelial mRNA was obtained from cystic and normal areas after laser capture microdissection (LCM). Transcriptomic analyses were performed and the results were confirmed by RT-PCR and immunohistochemistry in independent samples.

Results

After controlling for RNA quality, we analysed the transcriptomes of six cystic areas and five control areas. In total, 393 transcripts were differentially expressed in the epithelium, between CCAM and control areas. The most highly redundant genes involved in biological functions and signalling pathways differentially expressed between CCAM and control epithelium included TGFB2, TGFBR1, and MAP 2 K1. These genes were considered particularly relevant as they have been implicated in branching morphogenesis. RT-qPCR analysis confirmed in independent samples that TGFBR1 was more strongly expressed in CCAM than in control tissues (p < 0.03). Immunohistochemistry analysis showed TGFBR1 (p = 0.0007) and TGFB2 (p < 0.02) levels to be significantly higher in the epithelium of CCAM than in that of control tissues.

Conclusions

This compartmentalised transcriptomic analysis of the epithelium of macrocystic lung malformations identified a dysregulation of TGFB signalling at the mRNA and protein levels, suggesting a possible role of this pathway in CCAM pathogenesis.

Trial registration

ClinicalTrials.gov Identifier: NCT01732185.

In search of an evidence-based strategy for quality assessment of human tissue samples: report of the tissue Biospecimen Research Working Group of the Spanish Biobank Network

The purpose of the present work is to underline the importance of obtaining a standardized procedure to ensure and evaluate both clinical and research usability of human tissue samples. The study, which was carried out by the Biospecimen Science Working Group of the Spanish Biobank Network, is based on a general overview of the current situation about quality assurance in human tissue biospecimens. It was conducted an exhaustive review of the analytical techniques used to evaluate the quality of human tissue samples over the past 30 years, as well as their reference values if they were published, and classified them according to the biomolecules evaluated: (i) DNA, (ii) RNA, and (iii) soluble or/and fixed proteins for immunochemistry. More than 130 publications released between 1989 and 2019 were analysed, most of them reporting results focused on the analysis of tumour and biopsy samples. A quality assessment proposal with an algorithm has been developed for both frozen tissue samples and formalin-fixed paraffin-embedded (FFPE) samples, according to the expected quality of sample based on the available pre-analytical information and the experience of the participants in the Working Group. The high heterogeneity of human tissue samples and the wide number of pre-analytic factors associated to quality of samples makes it very difficult to harmonize the quality criteria. However, the proposed method to assess human tissue sample integrity and antigenicity will not only help to evaluate whether stored human tissue samples fit for the purpose of biomarker development, but will also allow to perform further studies, such as assessing the impact of different pre-analytical factors on very well characterized samples or evaluating the readjustment of tissue sample collection, processing and storing procedures. By ensuring the quality of the samples used on research, the reproducibility of scientific results will be guaranteed.

Considerations when processing and interpreting genomics data of the placenta

The application of genomic approaches to placental research has opened exciting new avenues to help us understand basic biological properties of the placenta, improve prenatal screening/diagnosis, and measure effects of in utero exposures on child health outcomes. In the last decade, such large-scale genomic data (including epigenomics and transcriptomics) have become more easily accessible to researchers from many disciplines due to the increasing ease of obtaining such data and the rapidly evolving computational tools available for analysis. While the potential of large-scale studies has been widely promoted, less attention has been given to some of the challenges associated with processing and interpreting such data. We hereby share some of our experiences in assessing data quality, reproducibility, and interpretation in the context of genome-wide studies of the placenta, with the aim to improve future studies. There is rarely a single "best" approach, as that can depend on the study question and sample cohort. However, being consistent, thoroughly assessing potential confounders in the data, and communicating key variables in the methods section of the manuscript are critically important to help researchers to collaborate and build on each other's work.

An Optimized Method of RNA Isolation from Goat Milk Somatic Cells for Transcriptomic Analysis

The goat (Capra hircus) is a perfect animal model for analyzing the transcriptome of milk somatic cells (MSCs), as sufficient numbers of somatic cells in goat milk, i.e., exfoliated epithelial cells, can be obtained using noninvasive methods. RNA integrity and purity are the first and most important parameters qualifying samples for transcriptomic tests and next-generation sequencing, as RNA quality influences experimental results. The aim of this study was to optimize a method for obtaining high-quality RNA from goat MSCs, irrespective of effects like breed, lactation stage, health status (e.g., with or without small ruminant lentivirus [SRLV] infection), or number of somatic cells. Milk samples were obtained from goats of two Polish breeds in various lactation stages and in different parities, and from goats infected and not infected with SRLV. Altogether, 412 MSC samples were examined: 206 using method A with fenozol and 206 using method B with QIAzol. Though the overall purity (measured as absorbance ratios at 260 nm/280 nm and 260 nm/230 nm) of the RNA material was comparable, the average yield of RNA isolated using method A was 11.9 µg, while method B’s average yield was 29.9 µg. Moreover, method B resulted in good quality RNA suitable for transcriptome analysis. Results were confirmed by RT-qPCR, using 18S rRNA and RPLP0 as the reference genes. The application of our modified treatment method was successful in obtaining high-integrity samples for transcriptomic or next-generation sequencing analysis. Using a 400 mL milk sample cooled in ice directly after milking, securing the cooling chain process from milking to MSC isolation, and applying method B to isolate RNA, we obtained good RNA quality irrespective of the goats’ breed, lactation stage, parity, milk yield, SRLV infection, and even milk yield and number of somatic cells in milk.

Parent-of-origin-specific allelic expression in the human placenta is limited to established imprinted loci and it is stably maintained across pregnancy

Background

Genomic imprinting, mediated by parent-of-origin-specific epigenetic silencing, adjusts the gene expression dosage in mammals. We aimed to clarify parental allelic expression in the human placenta for 396 claimed candidate imprinted genes and to assess the evidence for the proposed enrichment of imprinted expression in the placenta. The study utilized RNA-Seq-based transcriptome and genotyping data from 54 parental-placental samples representing the three trimesters of gestation, and term cases of preeclampsia, gestational diabetes, and fetal growth disturbances.

Results

Almost half of the targeted genes (n = 179; 45%) were either not transcribed or showed limited expression in the human placenta. After filtering for the presence of common exonic SNPs, adequacy of sequencing reads and informative families, 91 genes were retained (43 loci form Geneimprint database; 48 recently proposed genes). Only 11/91 genes (12.1%) showed confident signals of imprinting (binomial test, Bonferroni corrected P < 0.05; > 90% transcripts originating from one parental allele). The confirmed imprinted genes exhibit enriched placental expression (PHLDA2, H19, IGF2, MEST, ZFAT, PLAGL1, AIM1) or are transcribed additionally only in the adrenal gland (MEG3, RTL1, PEG10, DLK1). Parental monoallelic expression showed extreme stability across gestation and in term pregnancy complications. A distinct group of additional 14 genes exhibited a statistically significant bias in parental allelic proportions defined as having 65–90% of reads from one parental allele (e.g., KLHDC10, NLRP2, RHOBTB3, DNMT1). Molecular mechanisms behind biased parental expression are still to be clarified. However, 66 of 91 (72.5%) analyzed candidate imprinted genes showed no signals of deviation from biallelic expression.

Conclusions

As placental tissue is not included in The Genotype-Tissue Expression (GTEx) project, the study contributed to fill the gap in the knowledge concerning parental allelic expression. A catalog of parental allelic proportions and gene expression of analyzed loci across human gestation and in term pregnancy complications is provided as additional files. The study outcome suggested that true imprinting in the human placenta is restricted to well-characterized loci. High expression of imprinted genes during mid-pregnancy supports their critical role in developmental programming. Consistent with the data on other GTEx tissues, the number of human imprinted genes appears to be overestimated.

Electronic supplementary material

The online version of this article (10.1186/s13148-019-0692-3) contains supplementary material, which is available to authorized users.

The Effect of Genetic Variation on the Placental Transcriptome in Humans

The knowledge of genetic variants shaping human placental transcriptome is limited and they are not cataloged in the Genotype-Tissue Expression project. So far, only one whole genome analysis of placental expression quantitative trait loci (eQTLs) has been published by Peng et al. (2017) with no external independent validation. We report the second study on the landscape of placental eQTLs. The study aimed to generate a high-confidence list of placental cis-eQTLs and to investigate their potential functional implications. Analysis of cis-eQTLs (±100 kbp from the gene) utilized 40 placental RNA sequencing and respective whole genome genotyping datasets. The identified 199 placental cis-eSNPs represented 88 independent eQTL signals (FDR < 5%). The most significant placental eQTLs (FDR < 10-5) modulated the expression of ribosomal protein RPL9, transcription factor ZSCAN9 and aminopeptidase ERAP2. The analysis confirmed 50 eSNP-eGenes pairs reported by Peng et al. (2017) and thus, can be claimed as robust placental eQTL signals. The study identified also 13 novel placental eGenes. Among these, ZSCAN9 is modulated by several eSNPs (experimentally validated: rs1150707) that have been also shown to affect the methylation level of genes variably escaping X-chromosomal inactivation. The identified 63 placental eGenes exhibited mostly mixed or ubiquitous expression. Functional enrichment analysis highlighted 35 Gene Ontology categories with the top ranking pathways "ruffle membrane" (FDR = 1.81 × 10-2) contributing to the formation of motile cell surface and "ATPase activity, coupled" (FDR = 2.88 × 10-2), critical for the membrane transport. Placental eGenes were also significantly enriched in pathways implicated in development, signaling and immune function. However, this study was not able to confirm a significant overrepresentation of genome-wide association studies top hits among the placental eSNP and eGenes, reported by Peng et al. (2017). The identified eSNPs were further analyzed in association with newborn and pregnancy traits. In the discovery step, a suggestive association was detected between the eQTL of ALPG (rs11678251) and reduced placental, newborn's and infant's weight. Meta-analysis across REPROMETA, HAPPY PREGNANCY, ALSPAC cohorts (n = 6830) did not replicate these findings. In summary, the study emphasizes the role of genetic variation in driving the transcriptome profile of the human placenta and the importance to explore further its functional implications.

Efficient and High-Quality RNA Isolation from Metabolite-Rich Tissues of Stevia rebaudiana, an Important Commercial Crop

Stevia rebaudiana, a perennial herb native to northeastern Paraguay, has gained immense attention globally over the recent decades due to the natural sweetness of its leaves. Like in most plants, this particular species contains high amount of secondary metabolites, thus rendering the isolation of high quality and quantity RNA extract for molecular applications rather challenging. An effective, high-yield and high-quality RNA isolation protocol for this economically important plant species was devised here based on the cetyltrimethylammonium bromide (CTAB) extraction method, with an additional genomic DNA (gDNA) removal step. DNA and other contaminants that may affect downstream applications were effectively removed. Our results exhibited that RNA samples isolated from the leaves and stems of Stevia rebaudiana using this improvised method are high in integrity and quality with RNA integrity number (RIN) of more than 8 and low in contaminants.

Quality Control System in an Obstetrics and Gynecology Disease Biobank

AIM

To ensure that sample quality meets the requirements of experimental research, the gynecology and obstetrics biobank of the Nanjing Drum Tower hospital designed different quality control methods for relevant types of samples. A range of quality control procedures has been formulated.

METHODS

The sample types were frozen tissue, paraffin-embedded tissue, optimal cutting temperature (OCT)-embedded tissue, plasma, buffy coat, serum, blood clots, and urine. Different categories of samples from a random selection of 1% of cases were analyzed for quality control experiments: (i) frozen tissue, buffy coat, and blood clots: RNA and DNA were extracted and the concentration, purity, and integrity were determined; (ii) paraffin-embedded tissue: morphological observations were made after hematoxylin-eosin staining and immunohistochemical detection of β-actin or CD10; (iii) OCT-embedded tissue: hematoxylin-eosin staining and immunofluorescence detection of β-actin; and (iv) frozen tissue samples derived from different organs of 18 fetal autopsy specimens with different cold ischemia times (CITs), 0-12 hours, 12-18 hours, 18-24 hours, and 24-48 hours, were chosen to study RNA quality. There is no universally recognized quality control index for plasma, serum, and urine, so the quality of samples was evaluated from feedback from the research projects in which the samples were used.

RESULTS

Currently, there are ∼2000 cases and 360,000 sample vials in the biobank. According to the experiments, (i) the concentration and purity of all nucleic acids of selected samples were qualified; (ii) for frozen tissues with a CIT ≤1 hour, using a qualified standard RNA quality number (RQN) ≥7, the qualification rate was 90%; (iii) frozen tissues with CIT between 1 and 18 hours, using a qualified standard RQN ≥5, the qualification rate was 61.1%; (iv) all of the paraffin-embedded tissues qualified for morphological observation; (v) the qualification rate of OCT-embedded tissue was 89%; and (vi) CIT had a great influence on the integrity of frozen tissue RNA. As the tissue CIT lengthened, the integrity of the RNA decreased. The RNA integrity parameters of different tissue types in the same specimen were significantly different.

CONCLUSIONS

A quality control system was constructed in an obstetrics and gynecology disease biobank with various types of diseases and abundant samples. Using specific quality control experiments for different types of samples was a reliable operating strategy that can be beneficial for providing qualified research resources. For birth defect autopsy specimens, the samples used for RNA research should have a CIT of at least <12 hours.

Application of ribonucleoside vanadyl complex (RVC) for developing a multifunctional tissue preservative solution

The quality of biological samples greatly affects the accuracy of scientific results. However, RNA in cryopreserved tissues gradually degrades during storage, leading to errors in the results of subsequent experiments. A suitable sample preservative solution can prolong storage and enhance the research value of samples. Here, we developed a sample preservative solution using the properties of the ribonucleoside vanadyl complex (RVC) and compared its effects on RNA and DNA quality, protein activity, and tissue morphology with the commercially available and widely used RNAlater® Stabilization Solution. The results showed that both the RVC-based preservative solution and RNAlater can effectively delay RNA degradation in tissue samples stored at 4°C or −80°C compared with samples stored without any preservative solution. In contrast to RNAlater, the RVC-based preservative solution did not result in damage to the tissue morphology or a loss of protein activity. Additionally, the RVC-based preservative solution did not affect the RNA and genomic DNA contents of the tissue samples or the results of subsequent experimental analyses. An RVC-based reagent can be used as a multifunctional yet relatively inexpensive tissue preservative solution to provide a comprehensive and cost-effective method for preserving samples for tissue banks.