- Research article
- Open Access
- Open Peer Review
Epigenetic regulation of lncRNA connects ubiquitin-proteasome system with infection-inflammation in preterm births and preterm premature rupture of membranes
BMC Pregnancy and Childbirth volume 15, Article number: 35 (2015)
Preterm premature rupture of membranes (PPROM) is responsible for one third of all preterm births (PTBs). We have recently demonstrated that long noncoding RNAs (lncRNAs) are differentially expressed in human placentas derived from PPROM, PTB, premature rupture of the membranes (PROM), and full-term birth (FTB), and determined the major biological pathways involved in PPROM.
Here, we further investigated the relationship of lncRNAs, which are differentially expressed in spontaneous PTB (sPTB) and PPROM placentas and are found to overlap a coding locus, with the differential expression of transcribed mRNAs at the same locus. Ten lncRNAs (five up-regulated and five down-regulated) and the lncRNA-associated 10 mRNAs (six up- and four down-regulated), which were identified by microarray in comparing PPROM vs. sPTB, were then validated by real-time quantitative PCR.
A total of 62 (38 up- and 24 down-regulated) and 1,923 (790 up- and 1,133 down-regulated) lncRNAs were identified from placentas of premature labor (sPTB + PPROM), as compared to those from full-term labor (FTB + PROM) and from premature rupture of membranes (PPROM + PROM), as compared to those from non-rupture of membranes (sPTB + FTB), respectively. We found that a correlation existed between differentially expressed lncRNAs and their associated mRNAs, which could be grouped into four categories based on the gene strand (sense or antisense) of lncRNA and its paired transcript. These findings suggest that lncRNA regulates mRNA transcription through differential mechanisms. Differential expression of the transcripts PPP2R5C, STAM, TACC2, EML4, PAM, PDE4B, STAM, PPP2R5C, PDE4B, and EGFR indicated a co-expression among these mRNAs, which are involved in the ubiquitine-proteasome system (UPS), in addition to signaling transduction and beta adrenergic signaling, suggesting that imbalanced regulation of UPS may present an additional mechanism underlying the premature rupture of membrane in PPROM.
Differentially expressed lncRNAs that were identified from the human placentas of sPTB and PPROM may regulate their associated mRNAs through differential mechanisms and connect the ubiquitin-proteasome system with infection-inflammation pathways. Although the detailed mechanisms by which lncRNAs regulate their associated mRNAs in sPTB and PPROM are yet to be clarified, our findings open a new approach to explore the pathogenesis of sPTB and PPROM.
Preterm premature rupture of membranes (PPROM), which occurs in one-third of all preterm births (PTBs), is one of the major causes of prematurity . It accounts for a disproportionate amount of perinatal morbidity and mortality [2,3]. PPROM has been associated with several pathologic processes [4,5]. In addition, preterm contractions can lead to separation of the amnion and choriodecidua, with an overall reduction in membrane tensile strength, while cervical dilation can result in exposure of the membranes to vaginal microorganisms and reduce underlying tissue support . The molecular mechanisms, including epigenomic mechanisms, underlying sPTB and PPROM are not yet clearly understood.
Currently, high-throughput transcriptomic research has indicated that eukaryotic genomes transcribe up to 90% of the genomic DNA . Except for 1–2% of these transcripts encoding for proteins, the vast majority are transcribed as noncoding RNAs (ncRNAs) [7,8]. Although ncRNAs are defined by the lack of a protein-coding sequence, they can play important roles in a variety of biological processes [9,10]. ncRNAs are grouped into two major categories according to their lengths: the short noncoding RNAs, which include microRNAs (miRNAs) as well as other non-coding transcripts of less than 200 nucleotides (nts), and the more recently described long noncoding RNAs (lncRNAs) that are longer than 200 nt [11,12]. miRNAs, the most widely studied class of short ncRNAs, mediate post-transcriptional gene silencing by controlling the translation of mRNA into protein [13,14] and are involved in regulating proliferation, differentiation, apoptosis, and development . The disruption of expression of miRNAs has been found in many human diseases including cancers, neurological disorders, and cardiovascular disorders . For example, miR-15 and miR-16 are deregulated in B cell chronic lymphocytic leukemia ; miR-206 deficiency accelerates amyotrophic lateral sclerosis ; and miR-1, which is involved in heart development, has been linked with arrhythmias through down-regulating expression of ion channel genes [18,19]. In addition, the disruption of other classes of short ncRNAs, such as small nucleolar RNAs (snoRNAs) and piwi-interacting RNAs (piRNAs), can also lead to different human diseases . For example, the germline homozygous 2 bp (TT) deletion of the snoRNA U50 is associated with prostate cancer development , and overexpression of the piRNAs piwiL1 and piwiL2 is involved in somatic tumors [21-23]. Increasing evidence has revealed that lncRNAs can function as regulators of protein-coding gene expression and exert a variety of intrinsic functions in eukaryocytes . In genomic contexts, lncRNAs can be transcribed from enhancers, promoters, introns of genes, pseudogenes, and antisense to genes . They may influence almost every step in the life cycle of genes, and they carry out their biological roles through several different mechanisms, including regulating chromatin states and nuclear compartments [26-28], affecting the process of transcription [29-31], and mediating mRNA stability, splicing, and translation at the post-transcriptional level [32-34]. Recently, the disruption of lncRNAs was also found to be associated with different human diseases, as short ncRNAs were . ANRIL is the antisense lncRNA of the INK4 locus, and its altered activity could result in deregulated silencing of the INK4 locus, which contributes to the initiation of several cancers [36-39]. The lncRNA MALAT-1 is associated with early-stage non–small cell lung cancer , which depends on its ability in regulating the alternative splicing through interaction with nuclear phosphoproteins [41,42]. In addition, the antisense lncRNA BACE1-AS, which is encoded by the opposite strand of the gene BACE1, can increase BACE1 mRNA stability and protein abundance at the post-transcriptional level, as has been shown to be true in Alzheimer’s disease . Moreover, based on the recent expression analyses, multiple lines of evidence increasingly support the linkage of dysfunctions of lncRNAs to other human diseases, including neurodegenerative and psychiatric diseases , cardiovascular disease , and immune dysfunction and auto-immunity .
PPROM is a heterogeneous condition, although infection and inflammation are well documented as key etiological factors [38,39,46]. Earlier, we have hypothesized that the ncRNAs may play an epigenomic role in regulating the pathogenic development of PPROM. To test this, we applied microchip technology and identified more than a thousand placental lncRNAs. With these lncRNAs, we delineated more than 20 potential pathogenic pathways that are altered in PPROM . Among these, we showed that the pathways of infection and inflammatory response, extra cellular matrix (ECM)-receptor interactions, apoptosis, actin cytoskeleton, and smooth muscle contraction are the major pathogenic pathways involved in the development of PPROM. Here, we further investigated whether the differentially expressed lncRNAs, which were identified from placentas derived from PPROM, regulate related mRNA transcription.
The Ethics Committee of Lianyungang Maternal and Children’s Hospital reviewed and approved the research project. Written informed consent was obtained from the pregnant women who participated in this study. All material and data were previously coded and thus anonymous to the authors of this study.
A total of 40 placentas from age-matched pregnant women (25–30 years old) were divided into four groups of deliveries (10 placentas per group): sPTB, FTB, PPROM, and PROM, labeled as group A, B, C, and D, respectively. sPTB is defined as birth delivered at < 35 (= < 34+6) GW without premature rupture of amniochorionic membrane (PROM), FTB is defined as birth delivered at 39–40+6 GW without PROM; PPROM is defined as birth delivered at < 35 (= < 34+6) GW with PROM; and PROM is defined as birth delivered at 39–40+6 GW with premature rupture of the amniochorionic membrane before labor, which is the contraction of the uterus. Considering the epigenetic variation and high heterogeneity of PTB and PPROM, placentas were selected from pregnancies by using the criteria that there was no clinically recognized infection (no fever; no increase of white blood cell counts; no positive finding of amniotic fluid cultures; and no clinical intervention with antibiotics, steroid, or tocolytics during the pregnancy). Immediately following delivery, the placentas were flushed twice with 200 ml cold distilled water and dried with clean paper towels, sliced with a sterile scalpel into 1×1 cm2 cubes at a site 2 cm from the edge of the placenta, juxtaposed from the maternal side through the entire fetal membrane, quickly frozen in liquid nitrogen for 30 minutes, and stored at −80°C until use. The entire procedure was to be completed within 30 minutes after the placenta was delivered.
The Arraystar Human LncRNA Array v2.0 (www.arraystar.com) was employed in this study. This array covers 33,045 lncRNAs that were collected from the authoritative data sources RefSeq, UCSC Knowngenes, and Ensembl. RNA labeling and array hybridization were performed according to the Agilent One-Color Microarray-Based Gene Expression Analysis protocol (Agilent Technology, Santa Clara, CA) with minor modifications. Briefly, mRNA was purified from total RNA after removal of rRNA with mRNA-ONLY™ Eukaryotic mRNA Isolation Kit (Epicentre, Omaha, NE). Each sample was amplified and transcribed into fluorescent cRNA along the entire length of the transcripts without 3′ bias utilizing a random priming method. The labeled cRNAs were purified by RNeasy Mini Kit (Qiagen, Valencia, CA). The concentration and specific activity of the labeled cRNAs (pmol Cy3/μg cRNA) were measured by NanoDrop ND-1000. One μg of each labeled cRNA was fragmented by adding 5 μl 10 × blocking agent and 1 μl of 25 × fragmentation buffer, after which the mixture was heated at 60°C for 30 minutes, and 25 μl 2 × GE hybridization buffer was added to dilute the labeled cRNA. Fifty μl of hybridization solution was dispensed into the gasket slide and assembled to the lncRNA expression microarray slide. The slides were incubated for 17 hours at 65°C in an Agilent Hybridization Oven. The hybridized arrays were washed, fixed, and scanned by using the Agilent DNA Microarray Scanner (Agilent Technology, Santa Clara, CA). Agilent Feature Extraction software (version 22.214.171.124) was used to analyze the acquired array images. Quantile normalization and subsequent data processing were performed using the GeneSpring GX v12.1 software package (Agilent Technologies, Santa Clara, CA). After normalization of the raw data, lncRNAs and mRNAs that have flags (“All Targets Value”) were chosen for further data analysis. Differentially expressed lncRNAs and mRNAs with statistical significance between the two groups were identified through volcano plot filtering. Hierarchical clustering was performed using the Agilent GeneSpring GX software (Version 12.1). Both “GO analysis” and “Pathway analysis” were performed in the standard enrichment computation method. Results were also analyzed using the genetic and molecular interaction software GeneMANIA [48,49], an algorithm to indicate the co-relationship between these mRNA. The bio-functions and canonical pathways associated with our data were generated using the option of core-analysis in Ingenuity Pathway Analysis (IPA) (Ingenuity® Systems; http://www.ingenuity.com).
Quantitative real-time PCR analysis
Total RNA was extracted from placentas, and cDNA was synthesized. The expression level of lncRNAs, as well as of lncRNA-associated mRNAs, was determined by quantitative real-time PCR (RT-qPCR). Primer sequences used in qPCR were listed in Tables 1 and 2. qPCR reactions were performed by the ABI7900 system (Applied Biosystems, CA) and SYBR green dye PCR master mix (SuperArray, CA). Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used as an internal control, and lncRNAs’ or mRNAs’ values were normalized to GAPDH. Duplicated reactions were analyzed for each quantitative assay. For each lncRNA or mRNA, the result was finally reported as relative expression that was calculated relative to this control. All data were given in terms of relative expression of mean ± S.E. (N = 10). The data were subjected to one-way analysis of variance (one-way ANOVA) followed by an unpaired, two-tailed t-test. Differences were considered significant at P < 0.05 (labeled as “ * ” and extremely significant at P < 0.01 (labeled as “ ** ”).
Identification of differentially expressed lncRNAs in comparisons of group AC (PTB & PPROM) vs. group BD (FTB & PROM) and group CD (PPROM & PROM) vs. group AB (PTB & FTB)
Differential expression was both up- and down-regulated. A total of 62 (38 up- and 24 down-regulated) and 1,923 (790 up- and 1,133 down-regulated) lncRNAs were identified from group AC placentas compared to group BD, and from group CD compared to group AB, respectively (Figure 1A). To visualize differential expression between two different conditions, volcano plots (Figure 1B) were constructed by using fold-change (magnitude of change) values ≥ 2.0 and p-value < 0.05 as cut-offs, thus allowing visualization of the relationship between fold-change and statistical significance, which takes both magnitude of change and variability into consideration. The vertical lines correspond to 2.0-fold up and down, respectively, and the horizontal line represents a p-value of < 0.05. Therefore, the red point in the plot represents the differentially expressed lncRNAs with statistical significance. LncRNAs’ p-value < 0.0001, so those points were located in the axis Y = 4. Data of differentially expressed lncRNAs, generated by the microarray, has been deposited in Gene Expression Omnibus with an accession number GSE 50879 (http://www.ncbi.nlm.nih.gov/geo/info/linking.html).
Metabolic pathways involved in group AC (PTB & PPROM) vs. group BD (FTB & PROM) and group CD (PPROM & PROM) vs. group AB (PTB & FTB)
A functional analysis of mapping genes to KEGG (Kyoto Encyclopedia of Genes and Genomes, www.genome.jp/kegg/) pathways was performed with a p-value cut-off of 0.05. The p-value denotes the significance of the Pathway correlated to the conditions. The lower the p-value, the more significant is the Pathway. The up-regulated and down-regulated pathways with the Top 10 scores in group AC vs. group BD and group CD vs. group AB are shown in Figure 2, respectively.
Quantitative real-time PCR validation of differentially expressed lncRNAs and their associated mRNAs
In this study, 10 samples from each group were analyzed, which gave a total of 40 samples. To address the tissue/cell heterogeneity, we randomly selected five loci and performed an intragroup analysis. Our results showed that there was no significant difference (p > 0.05) among the 10 samples, but they did show intergroup differences that were statistically significant.
Ten paired lncRNAs and their associated mRNAs (PPP2R5C, STAM, TACC2, EML4, PAM, TATDN1, NBPF10, ADAMTS15, AQPEP, and PDE4B) with different expression patterns were selected from microarray data (Table 3). The selection was based on the following: 1) differentially expressed lncRNAs are natural antisense, 2) the sense strand of DNA has been transcribed into mRNA that is differentially expressed, and 3) the function of the gene is related to pathogenic pathway(s) that we have identified earlier from PTB and PPROM  or from this study as described above. RT-qPCR was used for validation of group PPROM compared to sPTB, PROM, and FTB; group sPTB compared to FTB and PROM; group PROM compared to FTB; the group of preterm labor [PPROM + sPTB] compared to full-term labor [PROM + FTB]; and the group of rupture of membrane [PPROM + PROM] compared to non-rupture of membrane [sPTB + FTB], respectively (Figure 3). The results of qPCR are summarized in Table 4. These comparisons revealed that, among 54 differentially expressed loci by microarray, one lncRNA (BC107431 in PPROM vs. sPTB) and seven mRNAs (ADAMTS15 in PPROM vs. sPTB; TATDN1 in PPROM vs. sPTB, PPROM vs. PROM, and [PPROM + PROM] vs. [sPTB + FTB]; AQPEP in PPROM vs. sPTB and sPTB vs. FTB; and TACC2 in PPROM vs. sPTB) by qPCR did not match the results of microarray, while 46 out of 54 (85.2%) were consistent with the results of microarray.
The target mRNAs of IncRNA interaction relationship analysis by GeneMANIA
Ten differentially transcribed mRNAs regulated by IncRNA, including PPP2R5C, STAM, TACC2, EML4, PAM, TATDN1, NBPF10, ADAMTS15, AQPEP, and PDE4B, were subjected to GeneMANIA analysis. Six—PPP2R5C, STAM, TACC2, EML4, PAM, and PDE4B—were found to be in a functional network in terms of co-expression (Figures 4 and 5). The core pathway analysis of the 10 loci of lncRNA-mRNA pairs showed that the top canonical pathway of the mRNAs involved was cardiac β-adrenergic signaling. Eight targeted mRNAs (ADAMTS15, EML4, NBPF10, PAM, PDE4B, PPP2R5C, STAM, TACC2, and TATDN1) and 27 ubiquitin protease pathway–associated genes (www.ncbi.nlm.nih.gov/geo/info/linking.html, accession # GSE 50879) could be composed together to construct a network (Figure 4). These genes can be subcategorized into different families, including complex, enzyme, kinase, peptidase, ion channel, phosphatase, and transporter, and their top function and disease was DNA replication, recombination, and repair; nucleic acid metabolism and small molecule biochemistry. The cellular locations of these molecules can be mapped by IPA program, as shown in Figure 5, and imply a possible signal transduction from extracellular to nucleus.
This current study is a continuing investigation to explore the molecular mechanisms of lncRNA’s regulatory functions in both PTB and PPROM, based on our findings that lncRNAs are differentially expressed in the placentas of PTB and PPROM . Histologically, the placenta contains both maternal (decidua) and fetal (chorionic villi and chorioamnionic membranes) tissues. Instead of using fetal membranes, which are essential for investigating the mechanisms directly related to the rupture of membrane, using the placenta may help us identify the lncRNAs that are differentially expressed in both maternal and fetal tissues and may provide a broad source of information related to the pathophysiological mechanisms underlying the premature contraction of the uterous, hemorrhage, and PROM, rather than only rupture of membrane. In addition, using the placenta as the study material may avoid the misidentifying of lncRNAs that maybe expressed in maternal tissues but not in fetal membranes. In our previous study, the lncRNAs from PPROM placentas and the pathogenic pathways they involved were identified based on microarray analysis . However, the regulatory pattern between lncRNAs and their associated mRNAs in PPROM is still unknown. In the present study, the relationship between lncRNAs and their associated mRNAs was detected, which may help us further understand the molecular mechanism underlying PPROM.
PPROM, defined as rupture of the membranes before 37 weeks of gestation, contains two major pathogenic factors, membrane rupture and preterm labor [50,51]. In order to study these two factors individually, groups with preterm labor or membrane rupture were analyzed by microarray. Compared to the full-term group with or without premature rupture of membrane (group BD), preterm labor is now thought to be a syndrome initiated by multiple mechanisms, including infection or inflammation, uteroplacental ischaemia or haemorrhage, uterine over-distension, stress, and other immunologically mediated processes . The up-regulated immune-related pathways in the preterm labor group (group AC) placentas further indicated that immunoreaction was one of the important risk factors in preterm labor. On the other hand, differentially expressed lncRNAs were identified from the premature membrane rupture group (group CD) placentas compared to the control group (group AB). Among the top ten up-regulated pathways, five were associated with microorganism infection, such as salmonella infection, influenza A, Chagas disease, measles, and legionellosis, and one (NOD-like receptor signaling pathway) was associated with the pathogen recognition process. One kind of etiologic factor linked to membrane rupture was microorganism infection and maternal or fetal host inflammatory response . Infection-induced activation of membrane matrix-specific enzymes (matrix metalloproteinases, MMPs) has recently been shown to be associated with excessive collagen turnover and membrane weakening, leading to rupture . The up-regulated pathways about microorganism infection and pathogen recognition further indicated that genital infections and inflammatory reactions were important pathological factors in membrane rupture. In addition, the up-regulated pathways in the membrane rupture group also contained the p53 signaling pathway, which was an important pathway leading to apoptosis and has been demonstrated to play an integral role in PPROM [4,53]. Up-regulated p53 signaling pathway might explain that microorganism infection may initiate membrane cell apoptosis through the p53 signaling pathway. The top ten down-regulated pathways in the premature membrane rupture group (group CD) were involved in amino acid, actin cytoskeleton regulation, and other signal transduction–related pathways. The down-regulation of these pathways might suggest that the placenta cells of the premature membrane rupture group have been in an apoptosis state and/or muscle contraction in a premature status, assuming the lncRNA(s) played a suppressing function for gene expression that resulted in increased transcription of mRNA(s).
It is possible for the lncRNAs to be located in the same strand or the antisense strand of their associated mRNA’s locus. Our microarray data showed that of the 10 paired lncRNAs and mRNAs studied, six pairs of locus TACC2, NBPF10, PPP2R5C, PAM, AQPEP, and PDE4B were located in the same strand (positive-strand), and the remaining four (ADAMTS15, EML4, STAM, and TATDN1) were in the antisense strands (Table 5). The differential expression pattern of lncRNAs between microarray data and RT-qPCR data among 27 comparisons (Table 3) showed 96% (26/27) agreement. However, the differential expression of the transcribed mRNA showed a lower rate (78%, 21/27) of agreement (Tables 3 and 4), suggesting that lncRNA is more stable than the coding mRNA. For example, both lncRNA and mRNA of ADAMTS15 were down-regulated in microarray assays; however, the lncRNA was down- but the mRNA was up-regulated in RT-qPCR. Nevertheless, Table 5 was constructed based on the qPCR results, in which we may see four categories of possible mechanisms by which lncRNA(s) regulate mRNA(s).
In the first category, lncRNAs and their associated mRNAs were located in the same strand, and both had the same expression pattern, which contains loci TACC2, NBPF10, and PPP2R5C. Recently, researchers have found that lncRNAs can compete for the miRNAs of mRNAs . Because both lncRNA and mRNA are transcribed as the positive strand (or leading strand, labeled as “+” in Table 5), the possible mechanism of lncRNA regulating mRNA is that lncRNAs may function to bind to miRNAs, which protects mRNA from the miRNA targeting and repressing. Consequently, the transcript of mRNA is up-regulated, along with the up-regulation of lncRNA. TACC2, which belongs to a conserved family of centrosome- and microtubule-interacting proteins and encodes a protein that is concentrated at centrosomes throughout the cell cycle , represented an example of the phenomenon that when lncRNA is up-regulated, the mRNA at the same locus would be up-regulated, too. PPP2R5C represents another example in which both lncRNA and mRNA were down-regulated. PPP2R5C is one of the four major Ser/Thr phosphatases . It was reported that a number of mammalian pseudogenes and lncRNAs have miRNA-binding sites in their 3′-UTRs and may therefore serve as “sponges” to sequester miRNAs away from miRNA-targets [56,57]. Considering that the lncRNAs and mRNAs of TACC2, NBPF10, and PPP2R5C were all in the sense strands and had the same expression pattern in regulating mRNA transcript in PPROM, the lncRNAs in the first category may regulate their associated mRNAs by competing for miRNA binding of mRNAs.
In the second category, lncRNAs and their associated mRNAs were also located in the same strand, but they had an opposite expression pattern, which contained loci of PAM, AQPEP, and PDE4B. In this group, lncRNA and mRNA might be transcribed bi-functionally. PDE4B is a member of the cyclic AMP (cAMP)-specific phosphodiesterase (PDE) family and encodes a protein that specifically hydrolyzes cAMP . Considering that cAMP is one of the important second messengers regulating and mediating a number of cellular responses to extracellular signals, such as hormones, light, and neurotransmitters [59-61], the down-regulated PDE4B mRNA in the PPROM group might imply that the cAMP-related signal transduction was involved in the formation of PPROM. Another gene, AQPEP, encoding a novel membrane-bound aminopeptidase termed aminopeptidase Q, was expressed on extravillous trophoblast (EVT) during placentation [62,63]. It can degrade several placenta-derived peptides including kisspeptin-10 , which suppresses trophoblast as well as cancer migration [64-66] and could also be involved in the regulation of EVT migration. The down-regulation of AQPEP mRNA in the PPROM group might indicate that the pathogenesis of PPROM was associated with the placental lesion. In addition, PAM, encoding a multifunctional protein with catalytic activities to catalyze neuroendocrine peptides to active alpha-amidated products [67,68], was up-regulated. Further studies are needed to investigate the role of this gene in PPROM. A recent interesting discovery about lncRNA indicated that SRA is a bifunctional transcript, with one isoform functioning as an ncRNA and another as an mRNA that is translated into the protein SRAP [69,70], which may in turn antagonize the function of its noncoding counterpart . Since lncRNAs and their associated mRNAs located in the same strand presented the opposite expression pattern, which is considered as bidirectional lncRNA in the second category, PAM, AQPEP, and PDE4B might possess similar transcripts as SRA in PPROM.
In the third category, lncRNAs and their associated mRNAs were located in different strands (lncRNAs in the antisense strand, and mRNAs in the sense strand) and had the opposite expression pattern, which contained genes of ADAMTS15, STAM, and EML4. ADAMTS15 encodes a member of the ADAMTS (a disintegrin and metalloproteinase with thrombospondin motifs) protein family. Proteins in the ADAMTS family share several distinct protein modules, including a propeptide region, a metalloproteinase domain, a disintegrin-like domain, and a thrombospondin type 1 motif [72,73]. Recently, researchers have found that degradation of extracellular matrix (ECM) by matrix metalloproteinases (MMPs) in fetal membrane plays significant roles in membrane rupture in PPROM [5,74]. Moreover, several other non-MMPs with a proteinase domain, such as serine proteases, cysteine proteases, and ADAMTS family members, can also break down amniochorion ECM substrates . The up-regulated mRNA expression of ADAMTS15 in the PPROM group might indicate that ADAMTS15 probably functioned as MMPs participating in the regulation of collagenolysis and ECM degradation in PPROM . Another interesting gene, STAM, encodes a member of the signal-transducing adaptor molecule family, which mediates downstream signaling of the cytokine-cytokine receptor interaction pathway . Host inflammatory responses induced by different pro-inflammatory cytokines (including IL-1β, IL-6, and TNF-α) were closely associated with PPROM [5,75]. The up-regulated expression of STAM mRNA might imply that STAM was probably involved in PPROM through regulation of pro-inflammatory cytokine-cytokine receptor interaction. To gain further understanding of the regulation mechanism of STAM in PPROM, KEGG analysis of this gene was processed. The results revealed that STAM was involved in the JAK-STAT signaling pathway (data not shown but available upon request), which could be activated by IL-6 and cytokine receptor interaction. Moreover, depending on the connection with the MAPK or PI3K-Akt signaling pathways, STAM could participate in the regulation of apoptosis. Apoptosis has been demonstrated to play an integral role in PPROM, and IL-6 also has been proved to promote MMP activation or apoptosis of fetal membranes in PPROM. Therefore, STAM might regulate apoptosis of fetal membranes in PPROM by mediating downstream signaling of the IL-6 and cytokine receptor interaction pathway. The last gene in this category was EML4, which encodes a novel microtubule-associated protein belonging to the conserved family of EMAP-like proteins. EML4 has been demonstrated to be essential for microtubule formation . Up-regulated expression of EML4 mRNA might imply that microtubule formation of the cytoskeleton in fetal membranes possibly was abnormal in PPROM. Because lncRNAs and mRNAs of ADAMTS15, STAM, and EML4 were located in opposite strands, lncRNAs in the third category could be complementary to their associated mRNAs. Although no evidence to date has demonstrated that lncRNAs can directly down-regulate their complementary mRNA through the RNAi-like pathway as miRNAs, the opposite expression pattern of antisense lncRNAs and sense mRNAs in this category could be explained by another new finding that lncRNAs could be the host genes for small RNAs . For example, lncRNA H19 is host to miR-675 , and the imprinted Gtl2, anti-Rtl1, and Mirg RNAs are hosts to almost 50 miRNAs and 40 snoRNAs . These lncRNAs could yield Dicer-dependent small RNAs and repress their complementary mRNAs through an RNAi-related pathway , which might be the mechanism whereby lncRNAs regulate mRNAs in this category.
In the fourth category, lncRNAs and their associated mRNAs were also located in opposite strands (lncRNAs in the antisense strand, and mRNAs in the sense strand), but they had the same expression pattern, which contained only one gene, TATDN1. TATDN1 encodes a protein known as a DNase domain-containing protein, which plays an important role in chromosomal segregation and cell cycle progression . TATDN1 is a conserved nuclease in both prokaryotes and eukaryotes and has been found to play a role in apoptotic DNA fragmentation in yeast and C. elegans [81,82]. Although it is unclear whether TATDN1 in vertebrates also functions in apoptotic DNA fragmentation, the down-regulation of TATDN1 mRNA in the PPROM group at least suggested that TATDN1 might be associated with PPROM by regulating DNA fragmentation. A recent study about Alzheimer’s disease found that the antisense transcript of the Alzheimer-associated b-secretase-1 (BACE), known as BACEAS, increases BACE mRNA stability , most likely by masking the binding sites for miR-485-5p , which suggested that lncRNAs could interfere with miRNA-mediated mRNA destabilization. Because lncRNAs and mRNA of TATDN1 were located in the opposite strand, and both were down-regulated in the PPROM group, lncRNAs of TATDN1 might regulate its mRNA as BACEAS by masking the miRNA-binding sites of its mRNA in PPROM.
Although the 10 lncRNA-targeting mRNAs may be regulated by lncRNA with different mechanisms, six of them—PPP2R5C, STAM, TACC2, EML4, PAM, and PDE4B—were found to be directly co-expressed [84-87] in one network by GeneMANIA analysis, which indicated that their expression levels are similar across conditions. Three of these mRNAs, STAM, PPP2R5C, and PDE4B, and another co-expressed gene, SMAD6, were found to be involved in the signal transduction pathway (REACTOME REACT_111102.4) . Signal transduction is a process in which extracellular signals elicit changes in cell state and activity. Depending on the cellular context, the signaling transduction pathway may impact cellular proliferation, differentiation, and survival. In these signaling pathways, PDE4B and PPP2R5C are involved in cardiac β-adrenergic signaling, which regulates the contraction of muscles. The earlier occurrence of uterine contraction was one of the main factors for PPROM. Therefore, the β-adrenergic signaling in uterine smooth muscles cells may be an important signaling pathway involved in PPROM. Further investigations should be pursued to explore how the mRNAs involved in this pathway are regulated by lncRNAs.
In addition to smooth muscle contraction, the ubiquitin C (UBC) network constructed by the eight lncRNA-targeted transcripts is shown in Figures 4 and 5. At the protein level, these targeted transcripts were found to interact with each other [89,90]. UBC is a polyubiquitin precursor. Conjugation of ubiquitin monomers or polymers can lead to various effects within a cell, depending on the residues to which ubiquitin is conjugated. Ubiquitination has been associated with protein degradation, DNA repair, cell cycle regulation, kinase modification, endocytosis, and regulation of other cell signaling pathways . Protein degradation was one of the pathophysiological changes for PPROM. Although MMP was one of the clarified factors associated with PPROM, the mechanism of PPROM was still unclear. Protein degradation through the ubiquitin-proteasome system (UPS) is the major pathway of non-lysosomal proteolysis of intracellular proteins and played a modulation role in the immune and inflammatory responses. Our analysis provided a possible explanation for the multiple pathways involved in PPROM under the regulation of UPS. In our analysis, STAM has two ubiquitin-binding domains, Vps27/Hrs/Stam (VHS) and ubiquitin-interacting motif (UIM), at its N-terminus, and its function was as a signal-transducing adaptor, which mediates downstream signaling of the cytokine-cytokine receptor interaction pathway. PPP2R5C is involved in multiple pathways, and its dephosphorization may also play a central role in regulation of target proteins’ degradation via the UPS. Our finding that transcription of mRNA of STAM may be regulated by lncRNA opens a new approach to study the pathogenesis of sPTB and PPROM through the connection between UPS and infection-inflammation pathways.
In conclusion, differentially expressed lncRNAs identified from the human placentas of PPROM might regulate their associated mRNAs through different mechanisms, as discussed above. The possible regulatory pattern between these lncRNAs and their associated mRNAs in PTB and PPROM might represent the potential molecular mechanism underlying PTB and PPROM. Although the detailed mechanisms through which lncRNAs regulate their associated mRNAs in sPTB and PPROM must be further clarified, our findings have identified a new way to explore the pathogenesis of sPTB and PPROM in the near future.
Lawn JE, Gravett MG, Nunes TM, Rubens CE, Stanton C. Global report on preterm birth and stillbirth (1 of 7): definitions, description of the burden and opportunities to improve data. BMC Pregnancy Childbirth. 2010;10 Suppl 1:S1.
Meis P, Ernest J, Moore M. Causes of low birth weight births in public and private patients. Am J Obstet Gynecol. 1987;156(5):1165–8.
Bott N, Nestor Demianczuk M. Influence of oligohydramnios on preterm premature rupture of the membranes at 30 to 36 weeks’ gestation. J Obstet Gynaecol Can. 2010;32(11):1030–4.
Srinivas S, Macones G. Preterm premature rupture of the fetal membranes: current concepts. Minerva Ginecol. 2005;57(4):389–96.
Romero R, Friel LA, Edwards DRV, Kusanovic JP, Hassan SS, Mazaki-Tovi S, et al. A genetic association study of maternal and fetal candidate genes that predispose to preterm prelabor rupture of membranes (PROM). Am J Obstet Gynecol. 2010;203:e1–361. e30. doi:101016/j.ajog.2010.05.026.
Strohl A, Kumar D, Novince R, Shaniuk P, Smith J, Bryant K, et al. Decreased adherence and spontaneous separation of fetal membrane layers–amnion and choriodecidua–a possible part of the normal weakening process. Placenta. 2010;31(1):18–24.
Feingold E, Good P, Guyer M, Kamholz S, Liefer L, Wetterstrand K, et al. The ENCODE (ENCyclopedia of DNA elements) project. Science. 2004;306(5696):636–40.
Kaikkonen MU, Lam MTY, Glass CK. Non-coding RNAs as regulators of gene expression and epigenetics. Cardiovasc Res. 2011;90(3):430–40.
Lipovich L, Johnson R, Lin C-Y. MacroRNA underdogs in a microRNA world: evolutionary, regulatory, and biomedical significance of mammalian long non-protein-coding RNA. Biochim Biophys Acta Gene Regulatory Mechanisms. 2010;1799(9):597–615.
Mattick JS. The genetic signatures of noncoding RNAs. PLoS Genet. 2009;5(4):e1000459.
Costa FF. Non-coding RNAs: Meet thy masters. Bioessays. 2010;32(7):599–608.
Jia H, Osak M, Bogu GK, Stanton LW, Johnson R, Lipovich L. Genome-wide computational identification and manual annotation of human long noncoding RNA genes. RNA. 2010;16(8):1478–87.
He L, Hannon GJ. MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet. 2004;5(7):522–31.
Mendell JT. MicroRNAs: critical regulators of development, cellular physiology and malignancy. Cell Cycle. 2005;4(9):1179–84.
Esteller M. Non-coding RNAs in human disease. Nat Rev Genet. 2011;12(12):861–74.
Calin GA, Dumitru CD, Shimizu M, Bichi R, Zupo S, Noch E, et al. Frequent deletions and down-regulation of micro-RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia. Proc Natl Acad Sci U S A. 2002;99(24):15524–9.
Williams AH, Valdez G, Moresi V, Qi X, McAnally J, Elliott JL, et al. MicroRNA-206 delays ALS progression and promotes regeneration of neuromuscular synapses in mice. Science. 2009;326(5959):1549–54.
Zhao Y, Ransom JF, Li A, Vedantham V, von Drehle M, Muth AN, et al. Dysregulation of cardiogenesis, cardiac conduction, and cell cycle in mice lacking miRNA-1-2. Cell. 2007;129(2):303–17.
Yang B, Lin H, Xiao J, Lu Y, Luo X, Li B, et al. The muscle-specific microRNA miR-1 regulates cardiac arrhythmogenic potential by targeting GJA1 and KCNJ2. Nat Med. 2007;13(4):486–91.
Dong X-Y, Guo P, Boyd J, Sun X, Li Q, Zhou W, et al. Implication of snoRNA U50 in human breast cancer. J Genet Genomics. 2009;36(8):447–54.
Taubert H, Greither T, Kaushal D, Würl P, Bache M, Bartel F, et al. Expression of the stem cell self-renewal gene Hiwi and risk of tumour-related death in patients with soft-tissue sarcoma. Oncogene. 2006;26(7):1098–100.
Liu JJ, Shen R, Chen L, Ye Y, He G, Hua K, et al. Piwil2 is expressed in various stages of breast cancers and has the potential to be used as a novel biomarker. Int J Clin Exp Pathol. 2010;3(4):328.
Sun G, Wang Y, Sun L, Luo H, Liu N, Fu Z, et al. Clinical significance of Hiwi gene expression in gliomas. Brain Res. 2011;1373:183–8.
Taft RJ, Pang KC, Mercer TR, Dinger M, Mattick JS. Non-coding RNAs: regulators of disease. J Pathol. 2010;220(2):126–39.
Kung JT, Colognori D, Lee JT. Long noncoding RNAs: past, present, and future. Genetics. 2013;193(3):651–69.
Campos EI, Reinberg D. Histones: annotating chromatin. Annu Rev Genet. 2009;43:559–99.
Pandey RR, Mondal T, Mohammad F, Enroth S, Redrup L, Komorowski J, et al. Kcnq1ot1 antisense noncoding RNA mediates lineage-specific transcriptional silencing through chromatin-level regulation. Mol Cell. 2008;32(2):232–46.
Zhang L-F, Huynh KD, Lee JT. Perinucleolar targeting of the inactive X during S phase: evidence for a role in the maintenance of silencing. Cell. 2007;129(4):693–706.
Hung T, Wang Y, Lin MF, Koegel AK, Kotake Y, Grant GD, et al. Extensive and coordinated transcription of noncoding RNAs within cell-cycle promoters. Nat Genet. 2011;43(7):621–9.
Kino T, Hurt DE, Ichijo T, Nader N, Chrousos GP. Noncoding RNA gas5 is a growth arrest-and starvation-associated repressor of the glucocorticoid receptor. Sci Signaling. 2010;3(107):ra8.
Sharma S, Findlay GM, Bandukwala HS, Oberdoerffer S, Baust B, Li Z, et al. Dephosphorylation of the nuclear factor of activated T cells (NFAT) transcription factor is regulated by an RNA-protein scaffold complex. Proc Natl Acad Sci U S A. 2011;108(28):11381–6.
Sone M, Hayashi T, Tarui H, Agata K, Takeichi M, Nakagawa S. The mRNA-like noncoding RNA Gomafu constitutes a novel nuclear domain in a subset of neurons. J Cell Sci. 2007;120(15):2498–506.
Matsui K, Nishizawa M, Ozaki T, Kimura T, Hashimoto I, Yamada M, et al. Natural antisense transcript stabilizes inducible nitric oxide synthase messenger RNA in rat hepatocytes. Hepatology. 2008;47(2):686–97.
Ebralidze AK, Guibal FC, Steidl U, Zhang P, Lee S, Bartholdy B, et al. PU. 1 expression is modulated by the balance of functional sense and antisense RNAs regulated by a shared cis-regulatory element. Genes Dev. 2008;22(15):2085–92.
Wapinski O, Chang HY. Long noncoding RNAs and human disease. Trends Cell Biol. 2011;21(6):354–61.
Yap K, Li S, Muñoz-Cabello A, Raguz S, Zeng L, Mujtaba S, et al. Molecular interplay of the noncoding RNA ANRIL and methylated histone H3 lysine 27 by polycomb CBX7 in transcriptional silencing of INK4a. Mol Cell. 2010;38(5):662.
Kotake Y, Nakagawa T, Kitagawa K, Suzuki S, Liu N, Kitagawa M, et al. Long non-coding RNA ANRIL is required for the PRC2 recruitment to and silencing of p15INK4B tumor suppressor gene. Oncogene. 2011;30(16):1956–62.
Iacobucci I, Sazzini M, Garagnani P, Ferrari A, Boattini A, Lonetti A, et al. A polymorphism in the chromosome 9p21 ANRIL locus is associated to Philadelphia positive acute lymphoblastic leukemia. Leuk Res. 2011;35(8):1052–9.
Pasmant E, Laurendeau I, Héron D, Vidaud M, Vidaud D, Bièche I. Characterization of a germ-line deletion, including the entire INK4/ARF locus, in a melanoma-neural system tumor family: identification of ANRIL, an antisense noncoding RNA whose expression coclusters with ARF. Cancer Res. 2007;67(8):3963–9.
Ji P, Diederichs S, Wang W, Böing S, Metzger R, Schneider PM, et al. MALAT-1, a novel noncoding RNA, and thymosin β4 predict metastasis and survival in early-stage non-small cell lung cancer. Oncogene. 2003;22(39):8031–41.
Bernard D, Prasanth KV, Tripathi V, Colasse S, Nakamura T, Xuan Z, et al. A long nuclear-retained non-coding RNA regulates synaptogenesis by modulating gene expression. EMBO J. 2010;29(18):3082–93.
Tripathi V, Ellis JD, Shen Z, Song DY, Pan Q, Watt AT, et al. The nuclear-retained noncoding RNA MALAT1 regulates alternative splicing by modulating SR splicing factor phosphorylation. Mol Cell. 2010;39(6):925–38.
Faghihi MA, Modarresi F, Khalil AM, Wood DE, Sahagan BG, Morgan TE, et al. Expression of a noncoding RNA is elevated in Alzheimer’s disease and drives rapid feed-forward regulation of β-secretase. Nat Med. 2008;14(7):723–30.
Faghihi MA, Wahlestedt C. Regulatory roles of natural antisense transcripts. Nat Rev Mol Cell Biol. 2009;10(9):637–43.
Annilo T, Kepp K, Laan M. Natural antisense transcript of natriuretic peptide precursor A (NPPA): structural organization and modulation of NPPA expression. BMC Mol Biol. 2009;10(1):81.
Cobo T, Jacobsson B, Kacerovsky M, Hougaard DM, Skogstrand K, Gratacos E, et al. Systemic and local inflammatory response in women with preterm prelabor rupture of membranes. PLoS One. 2014;9(1):e85277. doi:10.1371/journal.pone.0085277.
Luo X, Shi Q, Gu Y, Pan J, Hua M, Liu M, et al. LncRNA pathway involved in premature preterm rupture of membrane (PPROM): an epigenomic approach to study the pathogenesis of reproductive disorders. PLoS One. 2013;8(11):e79897.
Zuberi K, Franz M, Rodriguez H, Montojo J, Lopes CT, Bader GD, et al. GeneMANIA prediction server 2013 update. Nucleic Acids Res. 2013;41(W1):W115–22.
Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38 suppl 2:W214–20.
Menon R, Fortunato SJ. The role of matrix degrading enzymes and apoptosis in repture of membranes. J Soc Gynecol Investig. 2004;11(7):427–37.
Menon R, Fortunato SJ. Infection and the role of inflammation in preterm premature rupture of the membranes. Best Pract Res Clin Obstet Gynaecol. 2007;21(3):467–78.
Romero R, Espinoza J, Kusanovic JP, Gotsch F, Hassan S, Erez O, et al. The preterm parturition syndrome. BJOG. 2006;113(s3):17–42.
Kumar D, Fung W, Moore RM, Pandey V, Fox J, Stetzer B, et al. Proinflammatory cytokines found in amniotic fluid induce collagen remodeling, apoptosis, and biophysical weakening of cultured human fetal membranes. Biol Reprod. 2006;74(1):29–34.
K-i T, Horie-Inoue K, Suzuki T, Urano T, Ikeda K, Fujimura T, et al. TACC2 is an androgen-responsive cell cycle regulator promoting androgen-mediated and castration-resistant growth of prostate cancer. Mol Endocrinol. 2012;26(5):748–61.
Zheng H, Chen Y, Chen S, Niu Y, Yang L, Li B, et al. Expression and distribution of PPP2R5C gene in leukemia. J Hematol Oncol. 2011;4(21):1756–8722.
Cesana M, Cacchiarelli D, Legnini I, Santini T, Sthandier O, Chinappi M, et al. A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA. Cell. 2011;147(2):358–69.
Wang J, Liu X, Wu H, Ni P, Gu Z, Qiao Y, et al. CREB up-regulates long non-coding RNA, HULC expression through interaction with microRNA-372 in liver cancer. Nucleic Acids Res. 2010;38(16):5366–83.
Szpirer C, Szpirer J, Riviere M, Swinnen J, Vicini E, Conti M. Chromosomal localization of the human and rat genes (PDE4D and PDE4B) encoding the cAMP-specific phosphodiesterases 3 and 4. Cytogenet Genome Res. 1995;69(1–2):11–4.
Lauder JM. Neurotransmitters as growth regulatory signals: role of receptors and second messengers. Trends Neurosci. 1993;16(6):233–40.
Seger R, Krebs EG. The MAPK signaling cascade. FASEB J. 1995;9(9):726–35.
Meaney MJ, Diorio J, Francis D, Weaver S, Yau J, Chapman K, et al. Postnatal handling increases the expression of cAMP-inducible transcription factors in the rat hippocampus: the effects of thyroid hormones and serotonin. J Neurosci. 2000;20(10):3926–35.
Maruyama M, Hattori A, Goto Y, Ueda M, Maeda M, Fujiwara H, et al. Laeverin/aminopeptidase Q, a novel bestatin-sensitive leucine aminopeptidase belonging to the M1 family of aminopeptidases. J Biol Chem. 2007;282(28):20088–96.
Horie A, Fujiwara H, Sato Y, Suginami K, Matsumoto H, Maruyama M, et al. Laeverin/aminopeptidase Q induces trophoblast invasion during human early placentation. Hum Reprod. 2012;27(5):1267–76.
Bilban M, Ghaffari-Tabrizi N, Hintermann E, Bauer S, Molzer S, Zoratti C, et al. Kisspeptin-10, a KiSS-1/metastin-derived decapeptide, is a physiological invasion inhibitor of primary human trophoblasts. J Cell Sci. 2004;117(8):1319–28.
Navenot J-M, Wang Z, Chopin M, Fujii N, Peiper SC. Kisspeptin-10-induced signaling of GPR54 negatively regulates chemotactic responses mediated by CXCR4: a potential mechanism for the metastasis suppressor activity of kisspeptins. Cancer Res. 2005;65(22):10450–6.
Cho SG, Li D, Stafford LJ, Luo J, Rodriguez‐Villanueva M, Wang Y, et al. KiSS1 suppresses TNFα‐induced breast cancer cell invasion via an inhibition of RhoA‐Mediated NF‐κB activation. J Cell Biochem. 2009;107(6):1139–49.
Nieto L, Mascaraque A, Miller F, Glacial F, Ríos Martínez C, Kaiser M, et al. Synthesis and Antiprotozoal Activity of N-Alkoxy Analogues of the Trypanocidal Lead Compound 4, 4′-Bis (imidazolinylamino) diphenylamine with Improved Human Blood − Brain Barrier Permeability. J Med Chem. 2010;54(2):485–94.
Tsukamoto T, Noguchi M, Kayama H, Watanabe T, Asoh T, Yamamoto T. Increased peptidylglycine alpha-amidating monooxygenase activity in cerebrospinal fluid of patients with multiple sclerosis. Int Medicine (Tokyo, Japan). 1995;34(4):229–32.
Kawashima H, Takano H, Sugita S, Takahara Y, Sugimura K, Nakatani T. A novel steroid receptor co-activator protein (SRAP) as an alternative form of steroid receptor RNA-activator gene: expression in prostate cancer cells and enhancement of androgen receptor activity. Biochem J. 2003;369:163–71.
Chooniedass-Kothari S, Emberley E, Hamedani M, Troup S, Wang X, Czosnek A, et al. The steroid receptor RNA activator is the first functional RNA encoding a protein. FEBS Lett. 2004;566(1):43–7.
Hubé F, Velasco G, Rollin J, Furling D, Francastel C. Steroid receptor RNA activator protein binds to and counteracts SRA RNA-mediated activation of MyoD and muscle differentiation. Nucleic Acids Res. 2011;39(2):513–25.
Porter S, Clark I, Kevorkian L, Edwards D. The ADAMTS metalloproteinases. Biochem J. 2005;386:15–27.
Tanaka N, Kaneko K, Asao H, Kasai H, Endo Y, Fujita T, et al. Possible involvement of a novel STAM-associated molecule “AMSH” in intracellular signal transduction mediated by cytokines. J Biol Chem. 1999;274(27):19129–35.
Soydinc HE, Sak ME, Evliyaoglu O, Evsen MS, Turgut A, Ozler A, et al. Prolidase, matrix metalloproteinases 1 and 13 activity, oxidative-antioxidative status as a marker of preterm premature rupture of membranes and chorioamnionitis in maternal vaginal washing fluids. Int J Med Sci. 2013;10(10):1344–51. doi:10.7150/ijms.4802, PMCID: PMC3752722.
Hanada T, Yoshimura A. Regulation of cytokine signaling and inflammation. Cytokine Growth Factor Rev. 2002;13(4–5):413.
Pollmann M, Parwaresch R, Adam-Klages S, Kruse M-L, Buck F, Heidebrecht H-J. Human EML4, a novel member of the EMAP family, is essential for microtubule formation. Exp Cell Res. 2006;312(17):3241–51.
Jalali S, Bhartiya D, Lalwani MK, Sivasubbu S, Scaria V. Systematic transcriptome wide Analysis of lncRNA-miRNA interactions. PLoS ONE. 2013;8(2):e53823. doi:10.1371/journal.pone.0053823.
Keniry A, Oxley D, Monnier P, Kyba M, Dandolo L, Smits G, et al. The H19 lincRNA is a developmental reservoir of miR-675 that suppresses growth and Igf1r. Nat Cell Biol. 2012;14(7):659–65.
Rocha ST, Edwards CA, Ito M, Ogata T, Ferguson-Smith AC. Genomic imprinting at the mammalian Dlk1-Dio3 domain. Trends Genet. 2008;24(6):306–16.
Yang H, Liu C, Jamsen J, Wu Z, Wang Y, Chen J, et al. The DNase domain-containing protein TATDN1 plays an important role in chromosomal segregation and cell cycle progression during zebrafish eye development. Cell Cycle. 2012;11(24):4626–32.
Qiu J, Yoon J-H, Shen B. Search for apoptotic nucleases in yeast role of Tat-D nuclease in apoptotic DNA degradation. J Biol Chem. 2005;280(15):15370–9.
Parrish JZ, Xue D. Functional Genomic analysis of apoptotic DNA degradation in C. elegans. Mol Cell. 2003;11(4):987–96.
Caussy C, Charrière S, Marçais C, Di Filippo M, Sassolas A, Delay M, et al. An APOA5 3′UTR variant associated with plasma triglycerides triggers APOA5 downregulation by creating a functional miR-485-5p binding site. Am J Hum Genet. 2014;94(1):129.
Arijs I, De Hertogh G, Lemaire K, Quintens R, Van Lommel L, Van Steen K, et al. Mucosal gene expression of antimicrobial peptides in inflammatory bowel disease before and after first infliximab treatment. PLoS One. 2009;4(11):e7984.
Bild AH, Yao G, Chang JT, Wang Q, Potti A, Chasse D, et al. Oncogenic pathway signatures in human cancers as a guide to targeted therapies. Nature. 2006;439(7074):353–7.
Perou CM, Jeffrey SS, Van De Rijn M, Rees CA, Eisen MB, Ross DT, et al. Distinctive gene expression patterns in human mammary epithelial cells and breast cancers. Proc Natl Acad Sci U S A. 1999;96(16):9212–7.
Rieger KE, Hong W-J, Tusher VG, Tang J, Tibshirani R, Chu G. Toxicity from radiation therapy associated with abnormal transcriptional responses to DNA damage. Proc Natl Acad Sci U S A. 2004;101(17):6635–40.
Merico D, Isserlin R, Stueker O, Emili A, Bader GD. Enrichment map: a network-based method for gene-set enrichment visualization and interpretation. PLoS One. 2010;5(11):e13984.
Kim W, Bennett EJ, Huttlin EL, Guo A, Li J, Possemato A, et al. Systematic and quantitative assessment of the ubiquitin-modified proteome. Mol Cell. 2011;44(2):325–40.
Danielsen JM, Sylvestersen KB, Bekker-Jensen S, Szklarczyk D, Poulsen JW, Horn H, et al. Mass spectrometric analysis of lysine ubiquitylation reveals promiscuity at site level. Mol Cell Proteomics. 2011;10(3):M110. 003590.
Mukhopadhyay D, Riezman H. Proteasome-independent functions of ubiquitin in endocytosis and signaling. Science. 2007;315(5809):201–5.
The authors would like to thank the pregnant women and their families for participating in this study. We also thank KangChen Bio-tech (www.kangchen.com.cn) for their technical support. We appreciate everyone who helped in this research project.
This study was partly supported by the Jiangsu Provincial Natural Science Fund (BK20141244), the Lianyungang Maternal and Children’s Hospital, the New York State Institute for Basic Research in Developmental Disabilities, and the March of Dimes Foundation Global Network for Maternal and Infant Health.
The authors declare that they have no competing interests.
XL Participated in the design and coordination of the study, participated in lncRNA studies, and drafting the manuscript. JP: Participated in the design and coordination of the study, and participated in lncRNA studies. LW: Participated in preparing and drafting the manuscript and performed statistical analysis. PW: Participated in drafting the manuscript. MZ: Carried out specimen collection and processing. ML: Participated in data collection and processing, and drafting the manuscript. ZD: Participated in clinical case identifications, and drafting the manuscript. QM: Participated in clinical case identifications. XT: Participated in data processing. XZ and WJ: Participated in data analysis. JZ and EJ: Participated in drafting the manuscript. YG: Participated in coordinating clinical cases. WTB: Participated in data processing and drafting the manuscript. QS: Participated in coordinating clinical cases and seeking funding. NZ: Conceived, designed, and coordinated the study, consolidated funding, participated in drafting the manuscript, and final approved the publication. All authors read and approved the final manuscript.
Xiucui Luo and Jing Pan contributed equally to this work.