Identification of differentially expressed placental transcripts during multiple gestations in the Eurasian beaver (Castor fiber L.)

Identification and Characterization of Extrachromosomal Circular DNA in Human Placentas With Fetal Growth Restriction

Many studies have confirmed that extrachromosomal circular DNAs (eccDNAs/ecDNAs) exist in tumor and normal cells independently of the chromosome and are essential for oncogene plasticity and drug resistance. Studies have confirmed that there are many eccDNAs/ecDNAs in maternal plasma derived from the fetus. Fetal growth restriction (FGR) is a pregnancy-related disease associated with high newborn morbidity and mortality. However, the characteristics and nature of eccDNAs/ecDNAs in FGR are poorly understood. This study aims to deconstruct the properties and potential functions of eccDNAs/ecDNAs in FGR. We performed circle-seq to identify the expression profile of eccDNAs/ecDNAs, analyzed by bioinformatics, and verified by real-time Polymerase Chain Reaction (PCR) combined with southern blot in FGR compared with the normal groups. A total of 45,131 eccDNAs/ecDNAs (including 2,118 unique ones) were identified, which had significantly higher abundance in FRG group than in normal group, and was bimodal in length, peaking at ~146bp and ~340bp, respectively. Gestational age may be one independent factor affecting the production of eccDNAs/ecDNAs, most of which come from genomic regions with high gene density, with a 4~12bp repeat around the junction, and their origin had a certain genetic preference. In addition, some of the host-genes overlapped with non-coding RNAs (ncRNAs) partially or even completely. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis revealed that host-genes on the differentially expressed eccDNAs/ecDNAs (DEEECs/DEECs) were mainly enriched in immune-related functions and pathways. The presence of some ecDNAs were verified, and whose variability were consistent with the circle-seq results. We identified and characterized eccDNAs/ecDNAs in placentas with FGR, and elucidated the formation mechanisms and the networks with ncRNAs, which provide a new vision for the screening of new biomarkers and therapeutic targets for FGR.

Sex-Biased lncRNA Signature in Fetal Growth Restriction (FGR)

Impaired fetal growth is one of the most important causes of prematurity, stillbirth and infant mortality. The pathogenesis of idiopathic fetal growth restriction (FGR) is poorly understood but is thought to be multifactorial and comprise a range of genetic causes. This research aimed to investigate non-coding RNAs (lncRNAs) in the placentas of male and female fetuses affected by FGR. RNA-Seq data were analyzed to detect lncRNAs, their potential target genes and circular RNAs (circRNAs); a differential analysis was also performed. The multilevel bioinformatic analysis enabled the detection of 23,137 placental lncRNAs and 4263 of them were classified as novel. In FGR-affected female fetuses’ placentas (ff-FGR), among 19 transcriptionally active regions (TARs), five differentially expressed lncRNAs (DELs) and 12 differentially expressed protein-coding genes (DEGs) were identified. Within 232 differentially expressed TARs identified in male fetuses (mf-FGR), 33 encompassed novel and 176 known lncRNAs, and 52 DEGs were upregulated, while 180 revealed decreased expression. In ff-FGR ACTA2-AS1, lncRNA expression was significantly correlated with five DEGs, and in mf-FGR, 25 TARs were associated with DELs correlated with 157 unique DEGs. Backsplicing circRNA processes were detected in the range of H19 lncRNA, in both ff- and mf-FGR placentas. The performed global lncRNAs characteristics in terms of fetal sex showed dysregulation of DELs, DEGs and circRNAs that may affect fetus growth and pregnancy outcomes. In female placentas, DELs and DEGs were associated mainly with the vasculature, while in male placentas, disturbed expression predominantly affected immune processes.

Placenta Transcriptome Profiling in Intrauterine Growth Restriction (IUGR)

Intrauterine growth restriction (IUGR) is a serious pathological complication associated with compromised fetal development during pregnancy. The aim of the study was to broaden knowledge about the transcriptomic complexity of the human placenta by identifying genes potentially involved in IUGR pathophysiology. RNA-Seq data were used to profile protein-coding genes, detect alternative splicing events (AS), single nucleotide variant (SNV) calling, and RNA editing sites prediction in IUGR-affected placental transcriptome. The applied methodology enabled detection of 37,501 transcriptionally active regions and the selection of 28 differentially-expressed genes (DEGs), among them 10 were upregulated and 18 downregulated in IUGR-affected placentas. Functional enrichment annotation indicated that most of the DEGs were implicated in the processes of inflammation and immune disorders related to IUGR and preeclampsia. Additionally, we revealed that some genes (S100A13, GPR126, CTRP1, and TFPI) involved in the alternation of splicing events were mainly implicated in angiogenic-related processes. Significant SNVs were overlapped with 6533 transcripts and assigned to 2386 coding sequence (CDS), 1528 introns, 345 5' untranslated region (UTR), 1260 3'UTR, 918 non-coding RNA (ncRNA), and 10 intergenic regions. Within CDS regions, 543 missense substitutions with functional effects were recognized. Two known mutations (rs4575, synonymous; rs3817, on the downstream region) were detected within the range of AS and DEG candidates: PA28β and PINLYP, respectively. Novel genes that are dysregulated in IUGR were detected in the current research. Investigating genes underlying the IUGR is crucial for identification of mechanisms regulating placental development during a complicated pregnancy.

De novo characterization of placental transcriptome in the Eurasian beaver (Castor fiber L.)

Our pioneering data provide the first comprehensive view of placental transcriptome of the beaver during single and multiple gestation. RNA-Seq and a de novo approach allowed global pattern identification of C. fiber placental transcriptome. Non-redundant beaver transcriptome comprised 211,802,336 nt of placental transcripts, grouped into 128,459 contigs and clustered into 83,951 unigenes. An Ensembl database search revealed 14,487, 14,994, 15,004, 15,267 and 15,892 non-redundant homologs for Ictidomys tridecemlineatus, Rattus norvegicus, Mus musculus, Homo sapiens and Castor canadensis, respectively. Due to expression levels, the identified transcripts were divided into two sets: non-redundant and highly expressed (FPKM > 2 in at least three examined samples), analysed simultaneously. Among 17,009 highly expressed transcripts, 12,147 had BLASTx hits. GO annotations (175,882) were found for 4301 transcripts that were assigned to biological process (16,386), cellular component (9149) and molecular function (8338) categories; 666 unigenes were also classified into 122 KEGG pathways. Comprehensive analyses were performed for 411 and 3078 highly expressed transcripts annotated with a list of processes linked to ‘placenta’ (31 GO terms) or ‘embryo’ (324 GO terms), respectively. Among transcripts with entire CDS annotation, 281 (placenta) and 34 (embryo) alternative splicing events were identified. A total of 8499 putative SNVs (~ 6.2 SNV/transcript and 1.7 SNV/1 kb) were predicted with 0.1 minimum frequency and maximum variant quality (p value 10e⁻⁹). Our results provide a broad-based characterization of the global expression pattern of the beaver placental transcriptome. Enhancement of transcriptomic resources for C. fiber should improve understanding of crucial pathways relevant to proper placenta development and successful reproduction.

Preliminary RNA-Seq Analysis of Long Non-Coding RNAs Expressed in Human Term Placenta

Development of particular structures and proper functioning of the placenta are under the influence of sophisticated pathways, controlled by the expression of substantial genes that are additionally regulated by long non-coding RNAs (lncRNAs). To date, the expression profile of lncRNA in human term placenta has not been fully established. This study was conducted to characterize the lncRNA expression profile in human term placenta and to verify whether there are differences in the transcriptomic profile between the sex of the fetus and pregnancy multiplicity. RNA-Seq data were used to profile, quantify, and classify lncRNAs in human term placenta. The applied methodology enabled detection of the expression of 4463 isoforms from 2899 annotated lncRNA loci, plus 990 putative lncRNA transcripts from 607 intergenic regions. Those placentally expressed lncRNAs displayed features such as shorter transcript length, longer exon length, fewer exons, and lower expression levels compared to messenger RNAs (mRNAs). Among all placental transcripts, 175,268 were classified as mRNAs and 15,819 as lncRNAs, and 56,727 variants were discovered within unannotated regions. Five differentially expressed lncRNAs (HAND2-AS1, XIST, RP1-97J1.2, AC010084.1, TTTY15) were identified by a sex-bias comparison. Splicing events were detected within 37 genes and 4 lncRNA loci. Functional analysis of cis-related potential targets for lncRNAs identified 2021 enriched genes. It is presumed that the obtained data will expand the current knowledge of lncRNAs in placenta and human non-coding catalogs, making them more contemporary and specific.