RNA-Seq and expression microarray highlight different aspects of the fetal amniotic fluid transcriptome

RNA-Seq of amniotic fluid cell-free RNA: a discovery phase study of the pathophysiology of congenital cytomegalovirus infection

BACKGROUND

Congenital cytomegalovirus infection is the most common perinatal infection and a significant cause of sensorineural hearing loss, cerebral palsy, and neurodevelopmental disability. There is a paucity of human gene expression studies examining the pathophysiology of cytomegalovirus infection.

OBJECTIVE

This study aimed to perform a whole transcriptomic assessment of amniotic fluid from pregnancies with live fetuses to identify differentially expressed genes and enriched Gene Ontology categories associated with congenital cytomegalovirus infection.

STUDY DESIGN

Amniotic fluid supernatant was prospectively collected from pregnant women undergoing amniocentesis for suspected congenital cytomegalovirus infection because of first-trimester maternal primary infection or ultrasound features suggestive of fetal infection. Women who had received therapy to prevent fetal infection were excluded. Congenital cytomegalovirus infection was diagnosed via viral polymerase chain reaction of amniotic fluid; cytomegalovirus-infected fetuses were paired with noninfected controls, matched for gestational age and fetal sex. Paired-end RNA sequencing was performed on amniotic fluid cell-free RNA with the Novaseq 6000 at a depth of 30 million reads per sample. Following quality control and filtering, reads were mapped to the human genome and counts summarized across genes. Differentially expressed genes were identified using 2 approaches: voomWithQualityWeights in conjunction with limma and RUVSeq with edgeR. Genes with a false discovery rate <0.05 were considered statistically significant. Differential exon use was analyzed using DEXSeq. Functional analysis was performed using gene set enrichment analysis and Ingenuity Pathway Analysis. Manual curation of differentially regulated genes was also performed.

RESULTS

Amniotic fluid samples were collected from 50 women; 16 (32%) had congenital cytomegalovirus infection confirmed by polymerase chain reaction. After excluding 3 samples without matched controls, 13 cytomegalovirus-infected samples collected at 18 to 23 weeks and 13 cytomegalovirus-negative gestation-matched controls were submitted for RNA sequencing and analysis (N=26). Ten of the 13 pregnancies with cytomegalovirus-infected fetuses had amniocentesis because of serologic evidence of maternal primary infection with normal fetal ultrasound, and 3 had amniocentesis because of ultrasound abnormality suggestive of cytomegalovirus infection. Four cytomegalovirus-infected pregnancies ended in termination (n=3) or fetal death (n=1), and 9 resulted in live births. Pregnancy outcomes were available for 11 of the 13 cytomegalovirus-negative controls; all resulted in live births of clinically-well infants. Differential gene expression analysis revealed 309 up-regulated and 32 down-regulated genes in the cytomegalovirus-infected group compared with the cytomegalovirus-negative group. Gene set enrichment analysis showed significant enrichment of multiple Gene Ontology categories involving the innate immune response to viral infection and interferon signaling. Of the 32 significantly down-regulated genes, 8 were known to be involved in neurodevelopment and preferentially expressed by the brain. Six specific cellular restriction factors involved in host defense to cytomegalovirus infection were up-regulated in the cytomegalovirus-infected group. Ingenuity Pathway Analysis predicted the activation of pathways involved in progressive neurologic disease and inflammatory neurologic disease.

CONCLUSION

In this next-generation sequencing study, we revealed new insights into the pathophysiology of congenital cytomegalovirus infection. These data on the up-regulation of the intraamniotic innate immune response to cytomegalovirus infection and the dysregulation of neurodevelopmental genes may inform future approaches to developing prognostic markers and assessing fetal responses to in utero therapy.

Emerging technologies for prenatal diagnosis: The application of whole genome and RNA sequencing

DNA sequencing technologies for clinical genetic testing have been rapidly evolving in recent years, and steadily become more important within the field of prenatal diagnostics. This review aims to give an overview of recent developments and to describe how they have the potential to fill the gaps of the currently clinically implemented methods for prenatal diagnosis of various genetic disorders. It has been shown for postnatal testing that whole genome sequencing provides a set of added benefits compared to exome sequencing, and it is to be expected that this will be the case for prenatal testing as well. RNA-sequencing, already used postnatally, can provide valuable complementary data to DNA-based testing, and aid in variant interpretation. While not ready for clinical implementation, emerging technologies such as long-read and Hi-C sequencing analyses might add to the toolbox for interpreting the expanding genetic data sets generated by genome-wide sequencing. Lastly, we also discuss some more practical implications of introducing these emerging technologies, which generate larger and larger genomic data sets, in the prenatal field.

The Amniotic Fluid Cell-Free Transcriptome Provides Novel Information about Fetal Development and Placental Cellular Dynamics

The amniotic fluid (AF) is a complex biofluid that reflects fetal well-being during development. AF con be divided into two fractions, the supernatant and amniocytes. The supernatant contains cell-free components, including placenta-derived microparticles, protein, cell-free fetal DNA, and cell-free fetal RNA from the fetus. Cell-free mRNA (cfRNA) analysis holds a special position among high-throughput analyses, such as transcriptomics, proteomics, and metabolomics, owing to its ease of profiling. The AF cell-free transcriptome differs from the amniocyte transcriptome and alters with the progression of pregnancy and is often associated with the development of various organ systems including the fetal lung, skin, brain, pancreas, adrenal gland, gastrointestinal system, etc. The AF cell-free transcriptome is affected not only by normal physiologies, such as fetal sex, gestational age, and fetal maturity, but also by pathologic mechanisms such as maternal obesity, and genetic syndromes (Down, Edward, Turner, etc.), as well as pregnancy complications (preeclampsia, intrauterine growth restriction, preterm birth, etc.). cfRNA in the amniotic fluid originates from the placenta and fetal organs directly contacting the amniotic fluid as well as from the fetal plasma across the placenta. The AF transcriptome may reflect the fetal and placental development and therefore aid in the monitoring of normal and abnormal development.

Transcriptomic Analysis of Cell-free Fetal RNA in the Amniotic Fluid of Vervet Monkeys (Chlorocebus sabaeus)

NHP are important translational models for understanding the genomic underpinnings of growth, development, fetal programming, and predisposition to disease, with potential for the development of early health biomarkers. Understanding how prenatal gene expression is linked to pre- and postnatal health and development requires methods for assessing the fetal transcriptome. Here we used RNAseq methodology to analyze the expression of cell-free fetal RNA in the amniotic fluid supernatant (AFS) of vervet monkeys. Despite the naturally high level of degradation of free-floating RNA, we detected more than 10,000 gene transcripts in vervet AFS. The most highly expressed genes were H19, IGF2, and TPT1, which are involved in embryonic growth and glycemic health. We noted global similarities in expression profiles between vervets and humans, with genes involved in embryonic growth and glycemic health among the genes most highly expressed in AFS. Our study demonstrates both the feasibility and usefulness of prenatal transcriptomic profiles, by using amniocentesis procedures to obtain AFS and cell-free fetal RNA from pregnant vervets.

Utility of Circulating Cell-Free RNA Analysis for the Characterization of Global Transcriptome Profiles of Multiple Myeloma Patients

In this study, we evaluated the utility of extracellular RNA (exRNA) derived from the plasma of multiple myeloma (MM) patients for whole transcriptome characterization. exRNA from 10 healthy controls (HC), five newly diagnosed (NDMM), and 12 relapsed and refractory (RRMM) MM patients were analyzed and compared. We showed that ~45% of the exRNA genes were protein-coding genes and ~85% of the identified genes were covered >70%. Compared to HC, we identified 632 differentially expressed genes (DEGs) in MM patients, of which 26 were common to NDMM and RRMM. We further identified 54 and 191 genes specific to NDMM and RRMM, respectively, and these included potential biomarkers such as LINC00863, MIR6754, CHRNE, ITPKA, and RGS18 in NDMM, and LINC00462, PPBP, RPL5, IER3, and MIR425 in RRMM, that were subsequently validated using droplet digital PCR. Moreover, single nucleotide polymorphisms and small indels were identified in the exRNA, including mucin family genes that demonstrated different rates of mutations between NDMM and RRMM. This is the first whole transcriptome study of exRNA in hematological malignancy and has provided the basis for the utilization of exRNA to enhance our understanding of the MM biology and to identify potential biomarkers relevant to the diagnosis and prognosis of MM patients.

Targeted expression profiling by RNA-Seq improves detection of cellular dynamics during pregnancy and identifies a role for T cells in term parturition

Development of maternal blood transcriptomic markers to monitor placental function and risk of obstetrical complications throughout pregnancy requires accurate quantification of gene expression. Herein, we benchmark three state-of-the-art expression profiling techniques to assess in maternal circulation the expression of cell type-specific gene sets previously discovered by single-cell genomics studies of the placenta. We compared Affymetrix Human Transcriptome Arrays, Illumina RNA-Seq, and sequencing-based targeted expression profiling (DriverMap^TM) to assess transcriptomic changes with gestational age and labor status at term, and tested 86 candidate genes by qRT-PCR. DriverMap identified twice as many significant genes (q < 0.1) than RNA-Seq and five times more than microarrays. The gap in the number of significant genes remained when testing only protein-coding genes detected by all platforms. qRT-PCR validation statistics (PPV and AUC) were high and similar among platforms, yet dynamic ranges were higher for sequencing based platforms than microarrays. DriverMap provided the strongest evidence for the association of B-cell and T-cell gene signatures with gestational age, while the T-cell expression was increased with spontaneous labor at term according to all three platforms. We concluded that sequencing-based techniques are more suitable to quantify whole-blood gene expression compared to microarrays, as they have an expanded dynamic range and identify more true positives. Targeted expression profiling achieved higher coverage of protein-coding genes with fewer total sequenced reads, and it is especially suited to track cell type-specific signatures discovered in the placenta. The T-cell gene expression signature was increased in women who underwent spontaneous labor at term, mimicking immunological processes at the maternal-fetal interface and placenta.

Next-generation sequencing and prenatal 'omics: advanced diagnostics and new insights into human development

Prenatal genetics has evolved over the last decade to include application of new 'omics technologies to improve perinatal care. The clinical utility of these technologies when applied to direct fetal specimens from amniocentesis or chorionic villus sampling is being explored. In this review, we provide an overview of use of prenatal exome sequencing and role in evaluation of the structurally abnormal fetus, potential applications of genome sequencing, and finally, use of transcriptomics to assess placental and fetal well-being.

RNA-Seq and Expression Arrays: Selection Guidelines for Genome-Wide Expression Profiling

The development of genome-wide gene expression profiling technologies over the past two decades has produced great opportunity for researchers to explore the transcriptome and to better understand biological systems and their perturbation. In this chapter we provide an overview of microarray and massively parallel sequencing technologies and their application to gene expression analysis. We discuss factors that impact expression data generation and analysis that which should be considered in the application of these technology platforms. We further present the results of a simple illustration study to highlight performance similarities and differences in expression profiling of protein-coding mRNAs with each platform. Based on technical and analytical differences between the two platforms, reports in the literature comparing arrays and RNA-Seq for gene expression, and our own example study and experience, we provide recommendations for platform selection for gene expression studies.

Amniotic fluid transcriptomics reflects novel disease mechanisms in fetuses with myelomeningocele

BACKGROUND

Cell-free RNA in amniotic fluid supernatant reflects developmental changes in gene expression in the living fetus, which includes genes that are specific to the central nervous system. Although it has been previously shown that central nervous system-specific transcripts are present in amniotic fluid supernatant, it is not known whether changes in the amniotic fluid supernatant transcriptome reflect the specific pathophysiologic condition of fetal central nervous system disorders. In myelomeningocele, there is open communication between the central nervous system and amniotic fluid.

OBJECTIVES

The purpose of this study was to identify molecular pathophysiologic changes and novel disease mechanisms that are specific to myelomeningocele by the analysis of amniotic fluid supernatant cell-free RNA in fetuses with open myelomeningocele.

STUDY DESIGN

Amniotic fluid supernatant was collected from 10 pregnant women at the time of the open myelomeningocele repair in the second trimester (24.5±1.0 weeks); 10 archived amniotic fluid supernatant from sex and gestational age-matched euploid fetuses without myelomeningocele were used as controls (20.9±0.9 weeks). Differentially regulated gene expression patterns were analyzed with the use of human genome expression arrays.

RESULTS

Fetuses with myelomeningocele had 284 differentially regulated genes (176 up- and 108 down-regulated) in amniotic fluid supernatant. Known genes that were associated with myelomeningocele (PRICKLE2, GLI3, RAB23, HES1, FOLR1) and novel dysregulated genes were identified in association with neurodevelopment and neuronal regeneration (up-regulated, GAP43 and ZEB1) or axonal growth and guidance (down-regulated, ACAP1). Pathway analysis demonstrated a significant contribution of inflammation to disease and a broad influence of Wnt signaling pathways (Wnt1, Wnt5A, ITPR1).

CONCLUSION

Transcriptomic analyses of living fetuses with myelomeningocele with the use of amniotic fluid supernatant cell-free RNA demonstrated differential regulation of specific genes and molecular pathways relevant to this central nervous system disorder, which resulted in a new understanding of pathophysiologic changes. The data also suggested the importance of pathways that involve secondary disease, such as inflammation, in myelomeningocele. These newly identified pathways may lead to hypotheses that can test novel therapeutic targets as adjuncts to fetal surgical repair.

Extracellular RNA in aging

Since the discovery of extracellular RNA (exRNA) in circulation and other bodily fluids, there has been considerable effort to catalog and assess whether exRNAs can be used as markers for health and disease. A variety of exRNA species have been identified including messenger RNA and noncoding RNA such as microRNA (miRNA), small nucleolar RNA, transfer RNA, and long noncoding RNA. Age-related changes in exRNA abundance have been observed, and it is likely that some of these transcripts play a role in aging. In this review, we summarize the current state of exRNA profiling in various body fluids and discuss age-related changes in exRNA abundance that have been identified in humans and other model organisms. miRNAs, in particular, are a major focus of current research and we will highlight and discuss the potential role that specific miRNAs might play in age-related phenotypes and disease. We will also review challenges facing this emerging field and various strategies that can be used for the validation and future use of exRNAs as markers of aging and age-related disease. WIREs RNA 2017, 8:e1385. doi: 10.1002/wrna.1385 For further resources related to this article, please visit the WIREs website.

Systems biology evaluation of cell-free amniotic fluid transcriptome of term and preterm infants to detect fetal maturity

Background

Amniotic fluid (AF) is a proximal fluid to the fetus containing higher amounts of cell-free fetal RNA/DNA than maternal serum, thereby making it a promising source for identifying novel biomarkers that predict fetal development and organ maturation. Our aim was to compare AF transcriptomic profiles at different time points in pregnancy to demonstrate unique genetic signatures that would serve as potential biomarkers indicative of fetal maturation.

Methods

We isolated AF RNA from 16 women at different time points in pregnancy: 4 from 18 to 24 weeks, 6 from 34 to 36 weeks, and 6 from 39 to 40 weeks. RNA-sequencing was performed on cell-free RNA. Gene expression and splicing analyses were performed in conjunction with cell-type and pathway predictions.

Results

Sample-level analysis at different time points in pregnancy demonstrated a strong correlation with cell types found in the intrauterine environment and fetal respiratory, digestive and external barrier tissues of the fetus, using high-confidence cellular molecular markers. While some RNAs and splice variants were present throughout pregnancy, many transcripts were uniquely expressed at different time points in pregnancy and associated with distinct neonatal co-morbidities (respiratory distress and gavage feeding), indicating fetal immaturity.

Conclusion

The AF transcriptome exhibits unique cell/organ-selective expression patterns at different time points in pregnancy that can potentially identify fetal organ maturity and predict neonatal morbidity. Developing novel biomarkers indicative of the maturation of multiple organ systems can improve upon our current methods of fetal maturity testing which focus solely on the lung, and will better inform obstetrical decisions regarding delivery timing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-015-0138-5) contains supplementary material, which is available to authorized users.

A RNA-Seq Analysis of the Rat Supraoptic Nucleus Transcriptome: Effects of Salt Loading on Gene Expression

Magnocellular neurons (MCNs) in the hypothalamo-neurohypophysial system (HNS) are highly specialized to release large amounts of arginine vasopressin (Avp) or oxytocin (Oxt) into the blood stream and play critical roles in the regulation of body fluid homeostasis. The MCNs are osmosensory neurons and are excited by exposure to hypertonic solutions and inhibited by hypotonic solutions. The MCNs respond to systemic hypertonic and hypotonic stimulation with large changes in the expression of their Avp and Oxt genes, and microarray studies have shown that these osmotic perturbations also cause large changes in global gene expression in the HNS. In this paper, we examine gene expression in the rat supraoptic nucleus (SON) under normosmotic and chronic salt-loading SL) conditions by the first time using "new-generation", RNA sequencing (RNA-Seq) methods. We reliably detect 9,709 genes as present in the SON by RNA-Seq, and 552 of these genes were changed in expression as a result of chronic SL. These genes reflect diverse functions, and 42 of these are involved in either transcriptional or translational processes. In addition, we compare the SON transcriptomes resolved by RNA-Seq methods with the SON transcriptomes determined by Affymetrix microarray methods in rats under the same osmotic conditions, and find that there are 6,466 genes present in the SON that are represented in both data sets, although 1,040 of the expressed genes were found only in the microarray data, and 2,762 of the expressed genes are selectively found in the RNA-Seq data and not the microarray data. These data provide the research community a comprehensive view of the transcriptome in the SON under normosmotic conditions and the changes in specific gene expression evoked by salt loading.

The amniotic fluid transcriptome as a guide to understanding fetal disease

Numerous recent studies have shown the power of cell-free fetal RNA, obtained from amniotic fluid supernatant, to report on the development of the living fetus in real time. Examination of these transcripts on a genome-wide basis has led to new insights into the prenatal pathophysiology of multiple genetic, developmental, and environmental diseases. Each studied condition presents a unique, characteristic fetal transcriptome, which points to specific disrupted molecular pathways. These studies have also improved our knowledge of the normal development of the human fetus, revealing gestational age-related dynamic gene expression from a variety of organs. Analysis of the fetal transcriptome in normal and abnormal development has led to novel approaches for in utero prenatal treatment.