Effects of exercise on whole-blood transcriptome profile in children with overweight/obesity.

BACKGROUND: The current knowledge about the molecular mechanisms underlying the health benefits of exercise is still limited, especially in childhood. We set out to investigate the effects of a 20-week exercise intervention on whole-blood transcriptome profile (RNA-seq) in children with overweight/obesity. METHODS: Twenty-four children (10.21 ± 1.33 years, 46% girls) with overweight/obesity, were randomized to either a 20-week exercise program (intervention group; n = 10), or to a no-exercise control group (n = 14). Whole-blood transcriptome profile was analyzed using RNA-seq by STRT technique with GlobinLock technology. RESULTS: Following the 20-week exercise intervention program, 161 genes were differentially expressed between the exercise and the control groups among boys, and 121 genes among girls (p-value <0.05), while after multiple correction, no significant difference between exercise and control groups persisted in gene expression profiles (FDR >0.05). Genes enriched in GO processes and molecular pathways showed different immune response in boys (antigen processing and presentation, infections, and T cell receptor complex) and in girls (Fc epsilon RI signaling pathway) (FDR <0.05). CONCLUSION: These results suggest that 20-week exercise intervention program alters the molecular pathways involved in immune processes in children with overweight/obesity.

Targeted gene expression profiling for accurate endometrial receptivity testing.

Expressional profiling of the endometrium enables the personalised timing of the window of implantation (WOI). This study presents and evaluates a novel analytical pipeline based on a TAC-seq (Targeted Allele Counting by sequencing) method for endometrial dating. The expressional profiles were clustered, and differential expression analysis was performed on the model development group, using 63 endometrial biopsies spanning over proliferative (PE, n = 18), early-secretory (ESE, n = 18), mid-secretory (MSE, n = 17) and late-secretory (LSE, n = 10) endometrial phases of the natural cycle. A quantitative predictor model was trained on the development group and validated on sequenced samples from healthy women, consisting of 52 paired samples taken from ESE and MSE phases and five LSE phase samples from 31 individuals. Finally, the developed test was applied to 44 MSE phase samples from a study group of patients diagnosed with recurrent implantation failure (RIF). In validation samples (n = 57), we detected displaced WOI in 1.8% of the samples from fertile women. In the RIF study group, we detected a significantly higher proportion of the samples with shifted WOI than in the validation set of samples from fertile women, 15.9% and 1.8% (p = 0.012), respectively. The developed model was evaluated with an average cross-validation accuracy of 98.8% and an accuracy of 98.2% in the validation group. The developed beREADY screening model enables sensitive and dynamic detection of selected transcriptome biomarkers, providing a quantitative and accurate prediction of endometrial receptivity status.

Optimized single-cell RNA sequencing protocol to study early genome activation in mammalian preimplantation development.

Here, we present a modification of single-cell tagged reverse transcription protocol to study gene expression on a single-cell level or with limited RNA input. We describe different enzymes for reverse transcription and cDNA amplification, modified lysis buffer, and additional clean-up steps before cDNA amplification. We also detail an optimized single-cell RNA sequencing method for handpicked single cells, or tens to hundreds of cells, as input material to study mammalian preimplantation development. For complete details on the use and execution of this protocol, please refer to Ezer et al.1.

Human endometrial cell-type-specific RNA sequencing provides new insights into the embryo-endometrium interplay.

STUDY QUESTION: Which genes regulate receptivity in the epithelial and stromal cellular compartments of the human endometrium, and which molecules are interacting in the implantation process between the blastocyst and the endometrial cells? SUMMARY ANSWER: A set of receptivity-specific genes in the endometrial epithelial and stromal cells was identified, and the role of galectins (LGALS1 and LGALS3), integrin ß1 (ITGB1), basigin (BSG) and osteopontin (SPP1) in embryo-endometrium dialogue among many other protein-protein interactions were highlighted. WHAT IS KNOWN ALREADY: The molecular dialogue taking place between the human embryo and the endometrium is poorly understood due to ethical and technical reasons, leaving human embryo implantation mostly uncharted. STUDY DESIGN SIZE DURATION: Paired pre-receptive and receptive phase endometrial tissue samples from 16 healthy women were used for RNA sequencing. Trophectoderm RNA sequences were from blastocysts. PARTICIPANTS/MATERIALS SETTING METHODS: Cell-type-specific RNA-seq analysis of freshly isolated endometrial epithelial and stromal cells using fluorescence-activated cell sorting (FACS) from 16 paired pre-receptive and receptive tissue samples was performed. Endometrial transcriptome data were further combined in silico with trophectodermal gene expression data from 466 single cells originating from 17 blastocysts to characterize the first steps of embryo implantation. We constructed a protein-protein interaction network between endometrial epithelial and embryonal trophectodermal cells, and between endometrial stromal and trophectodermal cells, thereby focusing on the very first phases of embryo implantation, and highlighting the molecules likely to be involved in the embryo apposition, attachment and invasion. MAIN RESULTS AND THE ROLE OF CHANCE: In total, 499 epithelial and 581 stromal genes were up-regulated in the receptive phase endometria when compared to pre-receptive samples. The constructed protein-protein interactions identified a complex network of 558 prioritized protein-protein interactions between trophectodermal, epithelial and stromal cells, which were grouped into clusters based on the function of the involved molecules. The role of galectins (LGALS1 and LGALS3), integrin ß1 (ITGB1), basigin (BSG) and osteopontin (SPP1) in the embryo implantation process were highlighted. LARGE SCALE DATA: RNA-seq data are available at www.ncbi.nlm.nih.gov/geo under accession number GSE97929. LIMITATIONS REASONS FOR CAUTION: Providing a static snap-shot of a dynamic process and the nature of prediction analysis is limited to the known interactions available in databases. Furthermore, the cell sorting technique used separated enriched epithelial cells and stromal cells but did not separate luminal from glandular epithelium. Also, the use of biopsies taken from non-pregnant women and using spare IVF embryos (due to ethical considerations) might miss some of the critical interactions characteristic of natural conception only. WIDER IMPLICATIONS OF THE FINDINGS: The findings of our study provide new insights into the molecular embryo-endometrium interplay in the first steps of implantation process in humans. Knowledge about the endometrial cell-type-specific molecules that coordinate successful implantation is vital for understanding human reproduction and the underlying causes of implantation failure and infertility. Our study results provide a useful resource for future reproductive research, allowing the exploration of unknown mechanisms of implantation. We envision that those studies will help to improve the understanding of the complex embryo implantation process, and hopefully generate new prognostic and diagnostic biomarkers and therapeutic approaches to target both infertility and fertility, in the form of new contraceptives. STUDY FUNDING/COMPETING INTERESTS: This research was funded by the Estonian Research Council (grant PRG1076); Horizon 2020 innovation grant (ERIN, grant no. EU952516); Enterprise Estonia (grant EU48695); the EU-FP7 Marie Curie Industry-Academia Partnerships and Pathways (IAPP, grant SARM, EU324509); Spanish Ministry of Economy, Industry and Competitiveness (MINECO) and European Regional Development Fund (FEDER) (grants RYC-2016-21199, ENDORE SAF2017-87526-R, and Endo-Map PID2021-127280OB-100); Programa Operativo FEDER Andalucía (B-CTS-500-UGR18; A-CTS-614-UGR20), Junta de Andalucía (PAIDI P20_00158); Margarita Salas program for the Requalification of the Spanish University system (UJAR01MS); the Knut and Alice Wallenberg Foundation (KAW 2015.0096); Swedish Research Council (2012-2844); and Sigrid Jusélius Foundation; Academy of Finland. A.S.-L. is funded by the Spanish Ministry of Science, Innovation and Universities (PRE2018-085440). K.G.-D. has received consulting fees and/or honoraria from RemovAid AS, Norway Bayer, MSD, Gedeon Richter, Mithra, Exeltis, MedinCell, Natural cycles, Exelgyn, Vifor, Organon, Campus Pharma and HRA-Pharma and NIH support to the institution; D.B. is an employee of IGENOMIX. The rest of the authors declare no conflict of interest.

DUX4 is a multifunctional factor priming human embryonic genome activation.

Double homeobox 4 (DUX4) is expressed at the early pre-implantation stage in human embryos. Here we show that induced human DUX4 expression substantially alters the chromatin accessibility of non-coding DNA and activates thousands of newly identified transcribed enhancer-like regions, preferentially located within ERVL-MaLR repeat elements. CRISPR activation of transcribed enhancers by C-terminal DUX4 motifs results in the increased expression of target embryonic genome activation (EGA) genes ZSCAN4 and KHDC1P1. We show that DUX4 is markedly enriched in human zygotes, followed by intense nuclear DUX4 localization preceding and coinciding with minor EGA. DUX4 knockdown in human zygotes led to changes in the EGA transcriptome but did not terminate the embryos. We also show that the DUX4 protein interacts with the Mediator complex via the C-terminal KIX binding motif. Our findings contribute to the understanding of DUX4 as a regulator of the non-coding genome.

Prenatal diagnosis of a 46,XY karyotype female fetus with an SRY-associated gonadal dysgenesis, conceived through an intracytoplasmic sperm injection: a case report.

BACKGROUND: Permanent progression of paternal age and development of reproductive medicine lead to increase in number of children conceived with assisted reproductive techniques (ART). Although it is uncertain if ARTs have direct influence on offspring health, advanced paternal age, associated comorbidities and reduced fertility possess significant risks of genetic disorders to the offspring. With a broad implementation of a non-invasive prenatal testing (NIPT), more cases of genetic disorders, including sex discordance are revealed. Among biological causes of sex discordance are disorders of sexual development, majority of which are associated with the SRY gene. CASE PRESENTATION: We report a case of a non-invasive prenatal testing and ultrasound sex discordance in a 46,XY karyotype female fetus with an SRY pathogenic variant, who was conceived through an intracytoplasmic sperm injection (ICSI) due to severe oligozoospermia of the father. Advanced mean age of ICSI patients is associated with risk of de novo mutations and monogenic disorders in the offspring. Additionally, ICSI patients have higher risk to harbour infertility-predisposing mutations, including mutations in the SRY gene. These familial and de novo genetic factors predispose ICSI-conceived children to congenital malformations and might negatively affect reproductive health of ICSI-patients' offspring. CONCLUSIONS: Oligozoospermic patients planning assisted reproduction are warranted to undergo genetic counselling and testing for possible inherited and mosaic mutations, and risk factors for de novo mutations.

Generation of RNA sequencing libraries for transcriptome analysis of globin-rich tissues of the domestic dog.

We have developed a protocol for barcoded cDNA libraries of 48 samples to study gene expression across tissues in the domestic dog, Canis familiaris, by modifying the Single-Cell Tagged Reverse Transcription (STRT) protocol (Islam et al., 2012, 2014). The cDNA reads represent mRNA 5' ends, enabling the study of transcription start sites (TSS). Our modifications include longer UMIs for molecular counting and Globin-Lock® to deplete globin mRNAs that are abundant in blood and blood-rich tissues dominating all reads.

Systematic evaluation of NIPT aneuploidy detection software tools with clinically validated NIPT samples.

Non-invasive prenatal testing (NIPT) is a powerful screening method for fetal aneuploidy detection, relying on laboratory and computational analysis of cell-free DNA. Although several published computational NIPT analysis tools are available, no prior comprehensive, head-to-head accuracy comparison of the various tools has been published. Here, we compared the outcome accuracies obtained for clinically validated samples with five commonly used computational NIPT aneuploidy analysis tools (WisecondorX, NIPTeR, NIPTmer, RAPIDR, and GIPseq) across various sequencing depths (coverage) and fetal DNA fractions. The sample set included cases of fetal trisomy 21 (Down syndrome), trisomy 18 (Edwards syndrome), and trisomy 13 (Patau syndrome). We determined that all of the compared tools were considerably affected by lower sequencing depths, such that increasing proportions of undetected trisomy cases (false negatives) were observed as the sequencing depth decreased. We summarised our benchmarking results and highlighted the advantages and disadvantages of each computational NIPT software. To conclude, trisomy detection for lower coverage NIPT samples (e.g. 2.5M reads per sample) is technically possible but can, with some NIPT tools, produce troubling rates of inaccurate trisomy detection, especially in low-FF samples.

Small RNA expression and miRNA modification dynamics in human oocytes and early embryos.

Small noncoding RNAs (sRNAs) play important roles during the oocyte-to-embryo transition (OET), when the maternal phenotype is reprogrammed and the embryo genome is gradually activated. The transcriptional program driving early human development has been studied with the focus mainly on protein-coding RNAs, and expression dynamics of sRNAs remain largely unexplored. We profiled sRNAs in human oocytes and early embryos using an RNA-sequencing (RNA-seq) method suitable for low inputs of material. We show that OET in humans is temporally coupled with the transition from predominant expression of oocyte short piRNAs (os-piRNAs) in oocytes, to activation of microRNA (miRNA) expression in cleavage stage embryos. Additionally, 3' mono- and oligoadenylation of miRNAs is markedly increased in zygotes. We hypothesize that this may modulate the function or stability of maternal miRNAs, some of which are retained throughout the first cell divisions in embryos. This study is the first of its kind elucidating the dynamics of sRNA expression and miRNA modification along a continuous trajectory of early human development and provides a valuable data set for in-depth interpretative analyses.

Decidualized endometrial stromal cells present with altered androgen response in PCOS.

Hyperandrogenic women with PCOS show disrupted decidualization (DE) and placentation. Dihydrotestosterone (DHT) is reported to enhance DE in non-PCOS endometrial stromal cells (eSCCtrl); however, this has not been assessed in PCOS cells (eSCPCOS). Therefore, we studied the transcriptome profile of non-decidualized (non-DE) and DE eSCs from women with PCOS and Ctrl in response to short-term estradiol (E2) and/or progesterone (P4) exposure with/without (±) DHT. The non-DE eSCs were subjected to E2 ± DHT treatment, whereas the DE (0.5 mM 8-Br-cAMP, 96 h) eSCs were post-treated with E2 and P4 ± DHT, and RNA-sequenced. Validation was performed by immunofluorescence and immunohistochemistry. The results showed that, regardless of treatment, the PCOS and Ctrl samples clustered separately. The comparison of DE vs. non-DE eSCPCOS without DHT revealed PCOS-specific differentially expressed genes (DEGs) involved in mitochondrial function and progesterone signaling. When further adding DHT, we detected altered responses for lysophosphatidic acid (LPA), inflammation, and androgen signaling. Overall, the results highlight an underlying defect in decidualized eSCPCOS, present with or without DHT exposure, and possibly linked to the altered pregnancy outcomes. We also report novel factors which elucidate the mechanisms of endometrial dysfunction in PCOS.

Profiling of Extracellular Small RNAs Highlights a Strong Bias towards Non-Vesicular Secretion.

The extracellular environment consists of a plethora of molecules, including extracellular miRNA that can be secreted in association with extracellular vesicles (EVs) or soluble protein complexes (non-EVs). Yet, interest in therapeutic short RNA carriers lies mainly in EVs, the vehicles conveying the great majority of the biological activity. Here, by overexpressing miRNA and shRNA sequences in parent cells and using size exclusion liquid chromatography (SEC) to separate the secretome into EV and non-EV fractions, we saw that >98% of overexpressed miRNA was secreted within the non-EV fraction. Furthermore, small RNA sequencing studies of native miRNA transcripts revealed that although the abundance of miRNAs in EVs, non-EVs and parent cells correlated well (R2 = 0.69-0.87), quantitatively an outstanding 96.2-99.9% of total miRNA was secreted in the non-EV fraction. Nevertheless, though EVs contained only a fraction of secreted miRNAs, these molecules were stable at 37 °C in a serum-containing environment, indicating that if sufficient miRNA loading is achieved, EVs can remain delivery-competent for a prolonged period of time. This study suggests that the passive endogenous EV loading strategy might be a relatively wasteful way of loading miRNA to EVs, and active miRNA loading approaches are needed for developing advanced EV miRNA therapies in the future.

From late fatherhood to prenatal screening of monogenic disorders: evidence and ethical concerns.

BACKGROUND: With the help of ART, an advanced parental age is not considered to be a serious obstacle for reproduction anymore. However, significant health risks for future offspring hide behind the success of reproductive medicine for the treatment of reduced fertility associated with late parenthood. Although an advanced maternal age is a well-known risk factor for poor reproductive outcomes, understanding the impact of an advanced paternal age on offspring is yet to be elucidated. De novo monogenic disorders (MDs) are highly associated with late fatherhood. MDs are one of the major sources of paediatric morbidity and mortality, causing significant socioeconomic and psychological burdens to society. Although individually rare, the combined prevalence of these disorders is as high as that of chromosomal aneuploidies, indicating the increasing need for prenatal screening. With the help of advanced reproductive technologies, families with late paternity have the option of non-invasive prenatal testing (NIPT) for multiple MDs (MD-NIPT), which has a sensitivity and specificity of almost 100%. OBJECTIVE AND RATIONALE: The main aims of the current review were to examine the effect of late paternity on the origin and nature of MDs, to highlight the role of NIPT for the detection of a variety of paternal age-associated MDs, to describe clinical experiences and to reflect on the ethical concerns surrounding the topic of late paternity and MD-NIPT. SEARCH METHODS: An extensive search of peer-reviewed publications (1980-2021) in English from the PubMed and Google Scholar databases was based on key words in different combinations: late paternity, paternal age, spermatogenesis, selfish spermatogonial selection, paternal age effect, de novo mutations (DNMs), MDs, NIPT, ethics of late fatherhood, prenatal testing and paternal rights. OUTCOMES: An advanced paternal age provokes the accumulation of DNMs, which arise in continuously dividing germline cells. A subset of DNMs, owing to their effect on the rat sarcoma virus protein-mitogen-activated protein kinase signalling pathway, becomes beneficial for spermatogonia, causing selfish spermatogonial selection and outgrowth, and in some rare cases may lead to spermatocytic seminoma later in life. In the offspring, these selfish DNMs cause paternal age effect (PAE) disorders with a severe and even life-threatening phenotype. The increasing tendency for late paternity and the subsequent high risk of PAE disorders indicate an increased need for a safe and reliable detection procedure, such as MD-NIPT. The MD-NIPT approach has the capacity to provide safe screening for pregnancies at risk of PAE disorders and MDs, which constitute up to 20% of all pregnancies. The primary risks include pregnancies with a paternal age over 40 years, a previous history of an affected pregnancy/child, and/or congenital anomalies detected by routine ultrasonography. The implementation of NIPT-based screening would support the early diagnosis and management needed in cases of affected pregnancy. However, the benefits of MD-NIPT need to be balanced with the ethical challenges associated with the introduction of such an approach into routine clinical practice, namely concerns regarding reproductive autonomy, informed consent, potential disability discrimination, paternal rights and PAE-associated issues, equity and justice in accessing services, and counselling. WIDER IMPLICATIONS: Considering the increasing parental age and risks of MDs, combined NIPT for chromosomal aneuploidies and microdeletion syndromes as well as tests for MDs might become a part of routine pregnancy management in the near future. Moreover, the ethical challenges associated with the introduction of MD-NIPT into routine clinical practice need to be carefully evaluated. Furthermore, more focus and attention should be directed towards the ethics of late paternity, paternal rights and paternal genetic guilt associated with pregnancies affected with PAE MDs.

Nasal upregulation of CST1 in dog-sensitised children with severe allergic airway disease.

BACKGROUND: The clinical presentation of children sensitised to dog dander varies from asymptomatic to severe allergic airway disease, but the genetic mechanisms underlying these differences are not clear. The objective of the present study was to investigate nasal transcriptomic profiles associated with dog dander sensitisation in school children and to reveal clinical symptoms related with these profiles. METHODS: RNA was extracted from nasal epithelial cell brushings of children sensitised to dog dander and healthy controls. Blood sample analyses included IgE against dog dander, dog allergen molecules, other airborne and food allergens, basophil activation and white blood cell counts. Clinical history of asthma and rhinitis was recorded, and lung function was assessed (spirometry, methacholine provocation and exhaled nitric oxide fraction). RESULTS: The most overexpressed gene in children sensitised to dog dander compared to healthy controls was CST1, coding for Cystatin 1. A cluster of these children with enhanced CST1 expression showed lower forced expiratory volume in 1 s, increased bronchial hyperreactivity, pronounced eosinophilia and higher basophil allergen threshold sensitivity compared with other children sensitised to dog dander. In addition, multi-sensitisation to lipocalins was more common in this group. CONCLUSIONS: Overexpression of CST1 is associated with more severe allergic airway disease in children sensitised to dog dander. CST1 is thus a possible biomarker of the severity of allergic airway disease and a possible therapeutic target for the future treatment of airborne allergy.

Motor Function Deficits in the Estrogen Receptor Beta Knockout Mouse: Role on Excitatory Neurotransmission and Myelination in the Motor Cortex.

BACKGROUND: Male estrogen receptor beta (ERß) knockout (BERKO) mice display anxiety and aggression linked to, among others, altered serotonergic signaling in the basolateral amygdala and dorsal raphe, impaired cortical radial glia migration, and reduced GABAergic signaling. The effects on primary motor cortex (M1 cortex) and locomotor activity as a consequence of ERß loss have not been investigated. OBJECTIVE: The aim of this study was to determine whether locomotor activity is altered as a consequence of the changes in the M1 cortex. METHODS: The locomotor activity of male wild-type (WT) and BERKO mice was evaluated using the open-field and rotarod tests. Molecular changes in the M1 cortex were analyzed by RNA sequencing, electron microscopy, electrophysiology, and immunohistological techniques. In addition, we established oligodendrocyte (OL) cultures from WT and BERKO mouse embryonic stem cells to evaluate OL function. RESULTS: Locomotor profiling revealed that BERKO mice were more active than WT mice but had impaired motor coordination. Analysis of the M1 cortex pointed out differences in synapse function and myelination. There was a reduction in GABAergic signaling resulting in imbalanced excitatory and inhibitory neurotransmission as well as a defective OL differentiation accompanied by myelin defects. The effects of ERß loss on OL differentiation were confirmed in vitro. CONCLUSION: ERß is an important regulator of GABAergic interneurons and OL differentiation, which impacts on adult M1 cortex function and may be linked to increased locomotor activity and decreased motor coordination in BERKO mice.

Distinct whole-blood transcriptome profile of children with metabolic healthy overweight/obesity compared to metabolic unhealthy overweight/obesity.

BACKGROUND: Youth populations with overweight/obesity (OW/OB) exhibit heterogeneity in cardiometabolic health phenotypes. The underlying mechanisms for those differences are still unclear. This study aimed to analyze the whole-blood transcriptome profile (RNA-seq) of children with metabolic healthy overweight/obesity (MHO) and metabolic unhealthy overweight/obesity (MUO) phenotypes. METHODS: Twenty-seven children with OW/OB (10.1 ± 1.3 years, 59% boys) from the ActiveBrains project were included. MHO was defined as having none of the following criteria for metabolic syndrome: elevated fasting glucose, high serum triglycerides, low high-density lipoprotein-cholesterol, and high systolic or diastolic blood pressure, while MUO was defined as presenting one or more of these criteria. Inflammatory markers were additionally determined. Total blood RNA was analyzed by 5'-end RNA-sequencing. RESULTS: Whole-blood transcriptome analysis revealed a distinct pattern of gene expression in children with MHO compared to MUO children. Thirty-two genes differentially expressed were linked to metabolism, mitochondrial, and immune functions. CONCLUSIONS: The identified gene expression patterns related to metabolism, mitochondrial, and immune functions contribute to a better understanding of why a subset of the population remains metabolically healthy despite having overweight/obesity. IMPACT: A distinct pattern of whole-blood transcriptome profile (RNA-seq) was identified in children with metabolic healthy overweight/obesity (MHO) compared to metabolic unhealthy overweight/obesity (MUO) phenotype. The most relevant genes in understanding the molecular basis underlying the MHO/MUO phenotypes in children could be: RREB1, FAM83E, SLC44A1, NRG1, TMC5, CYP3A5, TRIM11, and ADAMTSL2. The identified whole-blood transcriptome profile related to metabolism, mitochondrial, and immune functions contribute to a better understanding of why a subset of the population remains metabolically healthy despite having overweight/obesity.

Identification of fetal unmodified and 5-hydroxymethylated CG sites in maternal cell-free DNA for non-invasive prenatal testing.

BACKGROUND: Massively parallel sequencing of maternal cell-free DNA (cfDNA) is widely used to test fetal genetic abnormalities in non-invasive prenatal testing (NIPT). However, sequencing-based approaches are still of high cost. Building upon previous knowledge that placenta, the main source of fetal circulating DNA, is hypomethylated in comparison to maternal tissue counterparts of cfDNA, we propose that targeting either unmodified or 5-hydroxymethylated CG sites specifically enriches fetal genetic material and reduces numbers of required analytical sequencing reads thereby decreasing cost of a test. METHODS: We employed uTOPseq and hmTOP-seq approaches which combine covalent derivatization of unmodified or hydroxymethylated CG sites, respectively, with next generation sequencing, or quantitative real-time PCR. RESULTS: We detected increased 5-hydroxymethylcytosine (5hmC) levels in fetal chorionic villi (CV) tissue samples as compared with peripheral blood. Using our previously developed uTOP-seq and hmTOP-seq approaches we obtained whole-genome uCG and 5hmCG maps of 10 CV tissue and 38 cfDNA samples in total. Our results indicated that, in contrast to conventional whole genome sequencing, such epigenomic analysis highly specifically enriches fetal DNA fragments from maternal cfDNA. While both our approaches yielded 100% accuracy in detecting Down syndrome in fetuses, hmTOP-seq maintained such accuracy at ultra-low sequencing depths using only one million reads. We identified 2164 and 1589 placenta-specific differentially modified and 5-hydroxymethylated regions, respectively, in chromosome 21, as well as 3490 and 2002 Down syndrome-specific differentially modified and 5-hydroxymethylated regions, respectively, that can be used as biomarkers for identification of Down syndrome or other epigenetic diseases of a fetus. CONCLUSIONS: uTOP-seq and hmTOP-seq approaches provide a cost-efficient and sensitive epigenetic analysis of fetal abnormalities in maternal cfDNA. The results demonstrated that T21 fetuses contain a perturbed epigenome and also indicated that fetal cfDNA might originate from fetal tissues other than placental chorionic villi. Robust covalent derivatization followed by targeted analysis of fetal DNA by sequencing or qPCR presents an attractive strategy that could help achieve superior sensitivity and specificity in prenatal diagnostics.

Fetal HLA-G mediated immune tolerance and interferon response in preeclampsia.

BACKGROUND: Fetal immune tolerance is crucial for pregnancy success. We studied the link between preeclampsia, a severe pregnancy disorder with uncertain pathogenesis, and fetal human leukocyte antigen G (HLA-G) and other genes regulating maternal immune responses. METHODS: We assessed sex ratios and regulatory HLA-G haplotypes in population cohorts and series of preeclampsia and stillbirth. We studied placental mRNA expression of 136 genes by sequencing and HLA-G and interferon alpha (IFNα) protein expression by immunohistochemistry. FINDINGS: We found underrepresentation of males in preeclamptic births, especially those delivered preterm or small for gestational age. Balancing selection at HLA-G associated with the sex ratio, stillbirth, and preeclampsia. We observed downregulation of HLA-G, its receptors, and many other tolerogenic genes, and marked upregulation of IFNA1 in preeclamptic placentas. INTERPRETATION: These findings indicate that an evolutionary trade-off between immune tolerance and protection against infections at the maternal-fetal interface promotes genetic diversity in fetal HLA-G, thereby affecting survival, preeclampsia, and sex ratio. We highlight IFNA1 as a potential mediator of preeclampsia and a target for therapeutic trials. FUNDING: Finnish Medical Foundation, Päivikki and Sakari Sohlberg Foundation, Karolinska Institutet Research Foundation, Scandinavia-Japan Sasakawa Foundation, Japan Eye Bank Association, Astellas Foundation for Research on Metabolic Disorders, Japan Society for the Promotion of Science, Knut and Alice Wallenberg Foundation, Swedish Research Council, Medical Society Liv och Hälsa, Sigrid Jusélius Foundation, Helsinki University Hospital and University of Helsinki, Jane and Aatos Erkko Foundation, Academy of Finland, Finska Läkaresällskapet, Novo Nordisk Foundation, Finnish Foundation for Pediatric Research, and Emil Aaltonen Foundation.

Method for the Identification of Plant DNA in Food Using Alignment-Free Analysis of Sequencing Reads: A Case Study on Lupin.

Fast and reliable analytical methods for the identification of plants from metagenomic samples play an important role in identifying the components of complex mixtures of processed biological materials, including food, herbal products, gut contents or environmental samples. Different PCR-based methods that are commonly used for plant identification from metagenomic samples are often inapplicable due to DNA degradation, a low level of successful amplification or a lack of detection power. We introduce a method that combines metagenomic sequencing and an alignment-free k-mer based approach for the identification of plant DNA in processed metagenomic samples. Our method identifies plant DNA directly from metagenomic sequencing reads and does not require mapping or assembly of the reads. We identified more than 31,000 Lupinus-specific 32-mers from assembled chloroplast genome sequences. We demonstrate that lupin DNA can be detected from controlled mixtures of sequences from target species (different Lupinus species) and closely related non-target species (Arachis hypogaea, Glycine max, Pisum sativum, Vicia faba, Phaseolus vulgaris, Lens culinaris, and Cicer arietinum). Moreover, these 32-mers are detectable in the following processed samples: lupin flour, conserved seeds and baked cookies containing different amounts of lupin flour. Under controlled conditions, lupin-specific components are detectable in baked cookies containing a minimum of 0.05% of lupin flour in wheat flour.

Dyslexia Candidate Gene and Ciliary Gene Expression Dynamics During Human Neuronal Differentiation.

Developmental dyslexia (DD) is a neurodevelopmental condition with complex genetic mechanisms. A number of candidate genes have been identified, some of which are linked to neuronal development and migration and to ciliary functions. However, expression and regulation of these genes in human brain development and neuronal differentiation remain uncharted. Here, we used human long-term self-renewing neuroepithelial stem (lt-NES, here termed NES) cells derived from human induced pluripotent stem cells to study neuronal differentiation in vitro. We characterized gene expression changes during differentiation by using RNA sequencing and validated dynamics for selected genes by qRT-PCR. Interestingly, we found that genes related to cilia were significantly enriched among upregulated genes during differentiation, including genes linked to ciliopathies with neurodevelopmental phenotypes. We confirmed the presence of primary cilia throughout neuronal differentiation. Focusing on dyslexia candidate genes, 33 out of 50 DD candidate genes were detected in NES cells by RNA sequencing, and seven candidate genes were upregulated during differentiation to neurons, including DYX1C1 (DNAAF4), a highly replicated DD candidate gene. Our results suggest a role of ciliary genes in differentiating neuronal cells and show that NES cells provide a relevant human neuronal model to study ciliary and DD candidate genes.