Search

Abstract Background Although plasma fibrinogen levels are related to cardiovascular risk, data regarding the role of fibrinogen genetic variation in myocardial infarction (MI) or coronary artery disease (CAD) etiology remain inconsistent. The purpose of the present study was to investigate the effect of fibrinogen A (FGA), fibrinogen B (FGB) and fibrinogen G (FGG) gene SNPs and haplotypes on susceptibility to CAD in a homogeneous Greek population. Methods We genotyped for rs2070022, rs2070016, rs2070006 in FGA gene, the rs7673587, rs1800789, rs1800790, rs1800788, rs1800787, rs4681 and rs4220 in FGB gene and for the rs1118823, rs1800792 and rs2066865 SNPs in FGG gene applying an arrayed primer extension-based genotyping method (APEX-2) in a sample of CAD patients (n = 305) and controls (n = 305). Logistic regression analysis was used to calculate odds ratios (ORs) and 95% confidence intervals (CIs), before and after adjustment for potential confounders. Results None of the FGA and FGG SNPs and FGA, FGB, FGG and FGA-FGG haplotypes was associated with disease occurrence after adjustment. Nevertheless, rs1800787 and rs1800789 SNPs in FGB gene seem to decrease the risk of CAD, even after adjustment for potential confounders (OR = 0.42, 95%CI: 0.19-0.90, p = 0.026 and OR = 0.44, 95%CI:0.21-0.94, p = 0.039, respectively). Conclusions FGA and FGG SNPs as well as FGA, FGB, FGG and FGA-FGG haplotypes do not seem to be important contributors to CAD occurrence in our sample. On the contrary, FGB rs1800787 and rs1800789 SNPs seem to confer protection to disease onset lowering the risk by about 50% in homozygotes for the minor alleles.

Introduction: Genital tract microbial communities are tightly associated with reproductive function, and their imbalance may explain reproductive failure in some cases. It is highly likely that the male genital tract microbiota has an impact on the female microbiota. Previous studies have demonstrated an increase in vaginosis-type microbiota and coliforms in women, as well as adverse effects of male genital tract bacteria on in vitro fertilization and pregnancy outcome. We profiled the seminal and vaginal microbiome of infertile couples by applying Illumina sequencing. Methods: We investigated 23 infertile couples. Semen samples were collected during the menstruation of their partner. Two self-collected vaginal samples were taken on the 6th to 8th days of the menstrual cycle, before and 8-12 h after intercourse. Seminal and vaginal microbiomes were profiled by applying the method which uses combinatorial sequence tags attached to PCR primers that amplify the rRNA V6 region (967-985 and 1078-1061, Escherichia coli 16S rRNA segment). Amplified PCR products were sequenced using an Illumina paired-end protocol on the HiSeq 2000 platform. Overlapping segments of the forward and reverse reads were analyzed with SHERA software. Mothur software was applied to trim, screen, and align sequences, calculate distances, assign sequences to operational taxonomic units, and obtain distance matrices. Results: Semen samples were clearly divided into two groups based on community composition. One group was dominated by Gardnerella vaginalis, Corynebacterium spp., Bifidobacterium breve, Peptoniphilus spp., Porphyromonas spp., Prevotella spp., Lactobacillus spp., Campylobacter rectus, Staphylococcus hominis, Fam. Marinilabiaceae, Fam. Comamonadaceae, and some unclassified Proteobacteria and Firmicutes. The second group was dominated by Enterobacteriaceae, Prevotella spp., Lactobacillus coleohominis, Lactobacillus iners, Diaphorobacter spp., Stenotrophomonas spp., and some unclassified Firmicutes, Actinobacteria, and Proteobacteria. In men with inflammation in their semen (>1 M WBC/ml), the mean proportion of Proteobacteria was higher than in non-inflammation men (p=0.045). Most of vaginal samples were dominated by lactobacilli (mainly Lactobacillus iners and Lactobacillus crispatus), while a small proportion of women had Gardnerella vaginalis or other bacteria (Streptococcus spp., Veillonella spp., Atopobium vaginae, Pseudomonas spp., Fam. Enterobacteriaceae). Seminal communities were significantly more diverse than vaginal ones, but with lower bacterial concentrations. In the couples with more similar genital tract microbiota, the similarity decreased after intercourse, while in the couples with less similar microbiota, the similarity increased after intercourse. Discussion: According to our knowledge, this is the first study to compare the seminal and vaginal microbiome and the effects of intercourse on the latter by applying Illumina sequencing. At the same time, the weakness of the study is the absence of a fertile control group; however, this type of study is quite complicated and healthy couples are not sufficiently motivated to participate. We revealed that the microbiomes of semen samples of the male partner in infertile couples are less abundant but highly diverse, while those of vaginal samples are more homogenous. The impact of the seminal microbiome during intercourse caused surprising but significant shifts in the vaginal microbiome, with the similarity of couples being decreased when they are more different before intercourse, while similarity is increased in couples with lower initial similarity. Further studies should determine the functional meaning of these phenomena. [ABSTRACT FROM AUTHOR]

Study Question: What changes occur in the endometrium during aging, and do they impact fertility?, Summary Answer: Both the transcriptome and cellular composition of endometrial samples from women of advanced maternal age (AMA) are significantly different from that of samples from young women, suggesting specific changes in epithelial cells that may affect endometrial receptivity., What Is Known Already: Aging is associated with the accumulation of senescent cells in aging tissues. Reproductive aging is mostly attributed to the decline in ovarian reserve and oocyte quality, whereas the endometrium is a unique complex tissue that is monthly renewed under hormonal regulation. Several clinical studies have reported lower implantation and pregnancy rates in oocyte recipients of AMA during IVF. Molecular studies have indicated the presence of specific mutations within the epithelial cells of AMA endometrium, along with altered gene expression of bulk endometrial tissue., Study Design Size Duration: Endometrial transcriptome profiling was performed for 44 women undergoing HRT during the assessment of endometrial receptivity before IVF. Patients younger than 28 years were considered as the young maternal age (YMA) group (age 23-27 years) and women older than 45 years were considered as the AMA group (age 47-50 years). Endometrial biopsies were obtained on Day 5 of progesterone treatment and RNA was extracted. All endometrial samples were evaluated as being receptive based on the expression of 68 common endometrial receptivity markers. Endometrial samples from another 24 women classified into four age groups (YMA, intermediate age group 1 (IMA1, age 29-35), intermediate age group 2 (IMA2, age 36-44), and AMA) were obtained in the mid-secretory stage of a natural cycle (NC) and used for validation studies across the reproductive lifespan., Participants/materials Setting Methods: A total of 24 HRT samples (12 YMA and 12 AMA) were subject to RNA sequencing (RNA-seq) and differential gene expression analysis, 20 samples (10 YMA and 10 AMA) were used for qPCR validation, and 24 NC samples (6 YMA, 6 IMA1, 6 IMA2 and 6AMA) were used for RNA-seq validation of AMA genes across the woman's reproductive lifespan. Immunohistochemistry (IHC) was used to confirm some expression changes at the protein level. Computational deconvolution using six endometrial cell type-specific transcriptomic profiles was conducted to compare the cellular composition between the groups., Main Results and the Role of Chance: Comparisons between YMA and AMA samples identified a lower proportion of receptive endometria in the AMA group ( P  = 0.007). Gene expression profiling identified 491 differentially expressed age-sensitive genes ( P adj < 0.05) that revealed the effects of age on endometrial epithelial growth and receptivity, likely contributing to decreased reproductive performance. Our results indicate that changes in the expression of the cellular senescence marker p16 INK4a and genes associated with metabolism, inflammation, and hormone response are involved in endometrial aging. Importantly, we demonstrate that the proportion of multi-ciliated cells, as discovered based on RNA-seq data deconvolution and tissue IHC results, is affected by endometrial aging, and propose a putative onset of age-related changes. Furthermore, we propose that aging has an impact on the transcriptomic profile of endometrial tissue in the context of endometrial receptivity., Large Scale Data: The raw sequencing data reported in this article are deposited at the Gene Expression Omnibus under accession code GSE236128., Limitations Reasons for Caution: This retrospective study identified changes in the endometrium of patients undergoing hormonal replacement and validated these changes using samples obtained during a NC. However, future studies must clarify the importance of these findings on the clinical outcomes of assisted reproduction., Wider Implications of the Findings: The findings reported in this study have important implications for devising future strategies aimed at improving fertility management in women of advanced reproductive age., Study Funding/competing Interests: This research was funded by the Estonian Research Council (grant no. PRG1076), Horizon 2020 innovation grant (ERIN, grant no. EU952516), Enterprise Estonia (grant no. EU48695), MSCA-RISE-2020 project TRENDO (grant no. 101008193), EU 874867 project HUTER, the Horizon Europe NESTOR grant (grant no. 101120075) of the European Commission, the EVA specialty program (grant no. KP111513) of the Maastricht University Medical Center (MUMC+), MICIU/AEI/10.13039/501100011033 and FEDER, EU projects Endo-Map (grant no. PID2021-12728OB-100), ROSY (grant no. CNS2022-135999), and the National Science Fund of Bulgaria (grant no. KII-06 H31/2). The authors declare no competing interests., Competing Interests: The authors declare no competing interests., (© The Author(s) 2024. Published by Oxford University Press on behalf of European Society of Human Reproduction and Embryology.)

Study Question: What is the effect of the chemical in vitro activation (cIVA) protocol compared with fragmentation only (Frag, also known as mechanical IVA) on gene expression, follicle activation and growth in human ovarian tissue in vitro ?, Summary Answer: Although histological assessment shows that cIVA significantly increases follicle survival and growth compared to Frag, both protocols stimulate extensive and nearly identical transcriptomic changes in cultured tissue compared to freshly collected ovarian tissue, including marked changes in energy metabolism and inflammatory responses., What Is Known Already: Treatments based on cIVA of the phosphatase and tensin homolog (PTEN)-phosphatidylinositol 3-kinase (PI3K) pathway in ovarian tissue followed by auto-transplantation have been administered to patients with refractory premature ovarian insufficiency (POI) and resulted in live births. However, comparable effects with mere tissue fragmentation have been shown, questioning the added value of chemical stimulation that could potentially activate oncogenic responses., Study Design Size Duration: Fifty-nine ovarian cortical biopsies were obtained from consenting women undergoing elective caesarean section (C-section). The samples were fragmented for culture studies. Half of the fragments were exposed to bpV (HOpic)+740Y-P (Frag+cIVA group) during the first 24 h of culture, while the other half were cultured with medium only (Frag group). Subsequently, both groups were cultured with medium only for an additional 6 days. Tissue and media samples were collected for histological, transcriptomic, steroid hormone, and cytokine/chemokine analyses at various time points., Participants/materials Setting Methods: Effects on follicles were evaluated by counting and scoring serial sections stained with hematoxylin and eosin before and after the 7-day culture. Follicle function was assessed by quantification of steroids by ultra-performance liquid chromatography tandem-mass spectrometry at different time points. Cytokines and chemokines were measured by multiplex assay. Transcriptomic effects were measured by RNA-sequencing (RNA-seq) of the tissue after the initial 24-h culture. Selected differentially expressed genes (DEGs) were validated by quantitative PCR and immunofluorescence in cultured ovarian tissue as well as in KGN cell (human ovarian granulosa-like tumor cell line) culture experiments., Main Results and the Role of Chance: Compared to the Frag group, the Frag+cIVA group exhibited a significantly higher follicle survival rate, increased numbers of secondary follicles, and larger follicle sizes. Additionally, the tissue in the Frag+cIVA group produced less dehydroepiandrosterone compared to Frag. Cytokine measurement showed a strong inflammatory response at the start of the culture in both groups. The RNA-seq data revealed modest differences between the Frag+cIVA and Frag groups, with only 164 DEGs identified using a relaxed cut-off of false discovery rate (FDR) <0.1. Apart from the expected PI3K-protein kinase B (Akt) pathway, cIVA also regulated pathways related to hypoxia, cytokines, and inflammation. In comparison to freshly collected ovarian tissue, gene expression in general was markedly affected in both the Frag+cIVA and Frag groups, with a total of 3119 and 2900 DEGs identified (FDR < 0.001), respectively. The top enriched gene sets in both groups included several pathways known to modulate follicle growth such as mammalian target of rapamycin (mTOR)C1 signaling. Significant changes compared to fresh tissue were also observed in the expression of genes encoding for steroidogenesis enzymes and classical granulosa cell markers in both groups. Intriguingly, we discovered a profound upregulation of genes related to glycolysis and its upstream regulator in both Frag and Frag+cIVA groups, and these changes were further boosted by the cIVA treatment. Cell culture experiments confirmed glycolysis-related genes as direct targets of the cIVA drugs. In conclusion, cIVA enhances follicle growth, as expected, but the mechanisms may be more complex than PI3K-Akt-mTOR alone, and the impact on function and quality of the follicles after the culture period remains an open question., Large Scale Data: Data were deposited in the GEO data base, accession number GSE234765. The code for sequencing analysis can be found in https://github.com/tialiv/IVA&#95;project., Limitations Reasons for Caution: Similar to the published IVA protocols, the first steps in our study were performed in an in vitro culture model where the ovarian tissue was isolated from the regulation of hypothalamic-pituitary-ovarian axis. Further in vivo experiments will be needed, for example in xeno-transplantation models, to explore the long-term impacts of the discovered effects. The tissue collected from patients undergoing C-section may not be comparable to tissue of patients with POI., Wider Implications of the Findings: The general impact of fragmentation and short (24 h) in vitro culture on gene expression in ovarian tissue far exceeded the effects of cIVA. Yet, follicle growth was stimulated by cIVA, which may suggest effects on specific cell populations that may be diluted in bulk RNA-seq. Nevertheless, we confirmed the impact of cIVA on glycolysis using a cell culture model, suggesting impacts on cellular signaling beyond the PI3K pathway. The profound changes in inflammation and glycolysis following fragmentation and culture could contribute to follicle activation and loss in ovarian tissue culture, as well as in clinical applications, such as fertility preservation by ovarian tissue auto-transplantation., Study Funding/competing Interests: This study was funded by research grants from European Union's Horizon 2020 Research and Innovation Programme (Project ERIN No. 952516, FREIA No. 825100), Swedish Research Council VR (2020-02132), StratRegen funding from Karolinska Institutet, KI-China Scholarship Council (CSC) Programme and the Natural Science Foundation of Hunan (2022JJ40782). International Iberian Nanotechnology Laboratory Research was funded by the European Union's H2020 Project Sinfonia (857253) and SbDToolBox (NORTE-01-0145-FEDER-000047), supported by Norte Portugal Regional Operational Programme (NORTE 2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund. No competing interests are declared., Competing Interests: The authors have no relevant financial or non-financial interests to disclose., (© The Author(s) 2024. Published by Oxford University Press on behalf of European Society of Human Reproduction and Embryology.)

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., (© 2023 The Authors. American Journal of Human Biology published by Wiley Periodicals LLC.)

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 ., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Author(s). Published by Elsevier Inc. All rights reserved.)

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., (© 2023. Springer Nature Limited.)

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&#95;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., (© The Author(s) 2022. Published by Oxford University Press on behalf of European Society of Human Reproduction and Embryology.)

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., Competing Interests: Y.T. and A.K. are employees of K.K.DNAFORM. Y. T. and Y. M. are inventors on a patent related to NET-CAGE technology., (© 2022 The Authors.)

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., (© 2022. The Author(s).)

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., Competing Interests: The authors have declared that no competing interests exist. GIPseq authors (A.A., B.B. and J.V.) performed all their computational analyses ’blindly’, without having any information about the analysed samples. They also did not participate in the analysis and interpretation of corresponding results.

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of European Society of Human Reproduction and Embryology. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Hyperandrogenic women with PCOS show disrupted decidualization (DE) and placentation. Dihydrotestosterone (DHT) is reported to enhance DE in non-PCOS endometrial stromal cells (eSC Ctrl ); however, this has not been assessed in PCOS cells (eSC PCOS ). 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 eSC PCOS 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 eSC PCOS , 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., (© 2021. The Author(s).)

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., (© 2021 Paloviita et al.; Published by Cold Spring Harbor Laboratory Press.)

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.

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., Competing Interests: Conflict of interest: U. Käck reports a lecture fee from Thermo Fisher outside the submitted work. Conflict of interest: E. Einarsdottir has nothing to disclose. Conflict of interest: M. van Hage reports lecture fees from Thermo Fisher Scientific and ALK, and consultancy fees from Biomay AG, Vienna, Austria and Hycor Biomedical LLC, CA, US, outside the submitted work. Conflict of interest: A. Asarnoj has nothing to disclose. Conflict of interest: A. James has nothing to disclose. Conflict of interest: A. Nopp has nothing to disclose. Conflict of interest: K. Krjutškov has nothing to disclose. Conflict of interest: S. Katayama reports grants from the Jane and Aatos Erkko Foundation during the conduct of the study. Conflict of interest: J. Kere has nothing to disclose. Conflict of interest: G. Lilja has nothing to disclose. Conflict of interest: C. Söderhäll has nothing to disclose. Conflict of interest: J.R. Konradsen has received material from Thermo Fisher Scientific to perform IgE analysis in this project., (Copyright ©The authors 2021. This version is distributed under the terms of the Creative Commons Attribution Licence 4.0.)

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., (© 2020 The Author(s) Published by S. Karger AG, Basel.)

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.

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., (Copyright © 2020 The Authors. Published by Elsevier B.V. All rights reserved.)

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.

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., (Copyright © 2020 Raime, Krjutškov and Remm.)

Background: Airway obstruction and wheezing in preschool children with recurrent viral infections are a major clinical problem, and are recognised as a risk factor for the development of chronic asthma. We aimed to analyse whether gene expression profiling provides evidence for pathways that delineate distinct groups of children with wheeze, and in combination with clinical information could contribute to diagnosis and prognosis of disease development., Methods: We analysed leukocyte transcriptomes from preschool children (6 months-3 years) at acute wheeze (n=107), and at a revisit 2-3 months later, comparing them to age-matched healthy controls (n=66). RNA-sequencing applying GlobinLock was used. The cases were followed clinically until age 7 years. Differential expression tests, weighted correlation network analysis and logistic regression were applied and correlations to 76 clinical traits evaluated., Findings: Significant enrichment of genes involved in the innate immune responses was observed in children with wheeze. We identified a unique acute wheeze-specific gene-module, which was associated with vitamin D levels (p<0.005) in infancy, and asthma medication and FEV 1 %/FVC (forced expiratory volume in 1 s/forced vital capacity) ratio several years later, at age 7 years (p<0.005). A model that predicts leukotriene receptor antagonist medication at 7 years of age with high accuracy was developed (area under the curve 0.815, 95% CI 0.668-0.962)., Interpretation: Gene expression profiles in blood from preschool wheezers predict asthma symptoms at school age, and therefore serve as biomarkers. The acute wheeze-specific gene module suggests that molecular phenotyping in combination with clinical information already at an early episode of wheeze may help to distinguish children who will outgrow their wheeze from those who will develop chronic asthma., Competing Interests: Conflict of interest: S. Katayama has nothing to disclose. Conflict of interest: K. Stenberg Hammar has nothing to disclose. Conflict of interest: K. Krjutškov has a patent “GlobinLock – blocking oligonucleotides” pending. Conflict of interest: E. Einarsdottir has nothing to disclose. Conflict of interest: G. Hedlin has nothing to disclose. Conflict of interest: J. Kere has acted as scientific advisor for Blueprint Genetics, Inc., Helsinki, Finland, outside the submitted work; and has a patent GlobinLock method pending. Conflict of interest: C. Söderhäll has nothing to disclose., (Copyright ©ERS 2020.)

Objective: The study aimed to validate a whole-genome sequencing-based NIPT laboratory method and our recently developed NIPTmer aneuploidy detection software with the potential to integrate the pipeline into prenatal clinical care in Estonia., Method: In total, 424 maternal blood samples were included. Analysis pipeline involved cell-free DNA extraction, library preparation and massively parallel sequencing on Illumina platform. Aneuploidies were determined with NIPTmer software, which is based on counting pre-defined per-chromosome sets of unique k-mers from sequencing raw data. SeqFF was implemented to estimate cell-free fetal DNA (cffDNA) fraction., Results: NIPTmer identified correctly all samples of non-mosaic trisomy 21 (T21, 15/15), T18 (9/9), T13 (4/4) and monosomy X (4/4) cases, with the 100% sensitivity. However, one mosaic T18 remained undetected. Six false-positive (FP) results were observed (FP rate of 1.5%, 6/398), including three for T18 (specificity 99.3%) and three for T13 (specificity 99.3%). The level of cffDNA of <4% was estimated in eight samples, including one sample with T13 and T18. Despite low cffDNA level, these two samples were determined as aneuploid., Conclusion: We believe that the developed NIPT method can successfully be used as a universal primary screening test in combination with ultrasound scan for the first trimester fetal examination., (© 2019 John Wiley & Sons, Ltd.)

Research Question: How does mucin MUC20 expression change during the menstrual cycle in different cell types of human endometrium?, Design: Study involved examination of MUC20 expression in two previously published RNA-seq datasets in whole endometrial tissue (n = 10), sorted endometrial epithelial (n = 44) or stromal (n = 42) cell samples. RNA-Seq results were validated by quantitative reverse transcription polymerase chain reaction (qRT-PCR) in whole tissue (n = 10), sorted epithelial (n = 17) and stromal (n = 17) cell samples. MUC20 protein localization and expression were analysed in human endometrium by immunohistochemical analysis of intact endometrial tissue (n = 6) and also Western blot of cultured stromal and epithelial cells (n = 2)., Results: MUC20 is differentially expressed in the endometrium between the pre-receptive and receptive phases. We show that MUC20 is predominantly expressed by epithelial cells of the receptive endometrium, both at the mRNA (RNA-Seq, P = 0.005; qRT-PCR, P = 0.039) and protein levels (Western blot; immunohistochemistry, P = 0.029)., Conclusion: Our results indicate MUC20 as a novel marker of mid-secretory endometrial biology. We propose a model of MUC20 function in the hepatocyte growth factor (HGF)-activated mesenchymal-epithelial transition (MET) receptor signalling specifically in the receptive phase. Further investigations should reveal the precise function of MUC20 in human endometrium and the possible connection between MUC20 and HGF-activated MET receptor signalling. MUC20 could potentially be included in the list of endometrial receptivity markers after further clinical validation., (Copyright © 2019 Reproductive Healthcare Ltd. Published by Elsevier Ltd. All rights reserved.)

RMRP was the first non-coding nuclear RNA gene implicated in a disease. Its mutations cause cartilage-hair hypoplasia (CHH), an autosomal recessive skeletal dysplasia with growth failure, immunodeficiency, and a high risk for malignancies. This study aimed to gain further insight into the role of RNA Component of Mitochondrial RNA Processing Endoribonuclease (RMRP) in cellular physiology and disease pathogenesis. We combined transcriptome analysis with single-cell analysis using fibroblasts from CHH patients and healthy controls. To directly assess cell cycle progression, we followed CHH fibroblasts by pulse-labeling and time-lapse microscopy. Transcriptome analysis identified 35 significantly upregulated and 130 downregulated genes in CHH fibroblasts. The downregulated genes were significantly connected to the cell cycle. Multiple other pathways, involving regulation of apoptosis, bone and cartilage formation, and lymphocyte function, were also affected, as well as PI3K-Akt signaling. Cell-cycle studies indicated that the CHH cells were delayed specifically in the passage from G2 phase to mitosis. Our findings expand the mechanistic understanding of CHH, indicate possible pathways for therapeutic intervention and add to the limited understanding of the functions of RMRP.

Background: Standard RNAseq methods using bulk RNA and recent single-cell RNAseq methods use DNA barcodes to identify samples and cells, and the barcoded cDNAs are pooled into a library pool before high throughput sequencing. In cases of single-cell and low-input RNAseq methods, the library is further amplified by PCR after the pooling. Preparation of hundreds or more samples for a large study often requires multiple library pools. However, sometimes correlation between expression profiles among the libraries is low and batch effect biases make integration of data between library pools difficult., Results: We investigated 166 technical replicates in 14 RNAseq libraries made using the STRT method. The patterns of the library biases differed by genes, and uneven library yields were associated with library biases. The former bias was corrected using the NBGLM-LBC algorithm, which we present in the current study. The latter bias could not be corrected directly, but could be solved by omitting libraries with particularly low yields. A simulation experiment suggested that the library bias correction using NBGLM-LBC requires a consistent sample layout. The NBGLM-LBC correction method was applied to an expression profile for a cohort study of childhood acute respiratory illness, and the library biases were resolved., Conclusions: The R source code for the library bias correction named NBGLM-LBC is available at https://shka.github.io/NBGLM-LBC and https://shka.bitbucket.io/NBGLM-LBC . This method is applicable to correct the library biases in various studies that use highly multiplexed sequencing-based profiling methods with a consistent sample layout with samples to be compared (e.g., "cases" and "controls") equally distributed in each library.

Non-invasive prenatal testing (NIPT) enables accurate detection of fetal chromosomal trisomies. The majority of publicly available computational methods for sequencing-based NIPT analyses rely on low-coverage whole-genome sequencing (WGS) data and are not applicable for targeted high-coverage sequencing data from cell-free DNA samples. Here, we present a novel computational framework for a targeted high-coverage sequencing-based NIPT analysis. The developed framework uses a hidden Markov model (HMM) in conjunction with a supplemental machine learning model, such as decision tree (DT) or support vector machine (SVM), to detect fetal trisomy and parental origin of additional fetal chromosomes. These models were developed using simulated datasets covering a wide range of biologically relevant scenarios with various chromosomal quantities, parental origins of extra chromosomes, fetal DNA fractions, and sequencing read depths. Developed models were tested on simulated and experimental targeted sequencing datasets. Consequently, we determined the functional feasibility and limitations of each proposed approach and demonstrated that read count-based HMM achieved the best overall classification accuracy of 0.89 for detecting fetal euploidies and trisomies on simulated dataset. Furthermore, we show that by using the DT and SVM on the HMM classification results, it was possible to increase the final trisomy classification accuracy to 0.98 and 0.99, respectively. We demonstrate that read count and allelic ratio-based models can achieve a high accuracy (up to 0.98) for detecting fetal trisomy even if the fetal fraction is as low as 2%. Currently, existing commercial NIPT analysis requires at least 4% of fetal fraction, which can be possibly a challenge in case of early gestational age (<10 weeks) or high maternal body mass index (>35 kg/m2). More accurate detection can be achieved at higher sequencing depth using HMM in conjunction with supplemental models, which significantly improve the trisomy detection especially in borderline scenarios (e.g., very low fetal fraction) and enables to perform NIPT even earlier than 10 weeks of pregnancy., Competing Interests: The authors have declared that no competing interests exist.

Pleomorphic adenoma gene 1 (PLAG1) is a transcription factor involved in cancer and growth. We discovered a de novo DNA motif containing a PLAG1 binding site in the promoters of genes activated during zygotic genome activation (ZGA) in human embryos. This motif was located within an Alu element in a region that was conserved in the murine B1 element. We show that maternally provided Plag1 is needed for timely mouse preimplantation embryo development. Heterozygous mouse embryos lacking maternal Plag1 showed disrupted regulation of 1,089 genes, spent significantly longer time in the 2-cell stage, and started expressing Plag1 ectopically from the paternal allele. The de novo PLAG1 motif was enriched in the promoters of the genes whose activation was delayed in the absence of Plag1. Further, these mouse genes showed a significant overlap with genes upregulated during human ZGA that also contain the motif. By gene ontology, the mouse and human ZGA genes with de novo PLAG1 motifs were involved in ribosome biogenesis and protein synthesis. Collectively, our data suggest that PLAG1 affects embryo development in mice and humans through a conserved DNA motif within Alu/B1 elements located in the promoters of a subset of ZGA genes.

Objective: The effect of vitamin D at the transcriptome level is poorly understood, and furthermore, it is unclear if it differs between obese and normal-weight subjects. The objective of the study was to explore the transcriptome effects of vitamin D supplementation., Design and Methods: We analysed peripheral blood gene expression using GlobinLock oligonucleotides followed by RNA sequencing in individuals participating in a 12-week randomised double-blinded placebo-controlled vitamin D intervention study. The study involved 18 obese and 18 normal-weight subjects (of which 20 males) with mean (±s.d.) age 20.4 (±2.5) years and BMIs 36 (±10) and 23 (±4) kg/m2, respectively. The supplemental daily vitamin D dose was 50 µg (2000 IU). Data were available at baseline, 6- and 12-week time points and comparisons were performed between the vitamin D and placebo groups separately in obese and normal-weight subjects., Results: Significant transcriptomic changes were observed at 6 weeks, and only in the obese subjects: 1724 genes were significantly upregulated and 186 genes were downregulated in the vitamin D group compared with placebo. Further analyses showed several enriched gene categories connected to mitochondrial function and metabolism, and the most significantly enriched pathway was related to oxidative phosphorylation (adjusted P value 3.08 × 10-14). Taken together, our data suggest an effect of vitamin D supplementation on mitochondrial function in obese subjects., Conclusions: Vitamin D supplementation affects gene expression in obese, but not in normal-weight subjects. The altered genes are enriched in pathways related to mitochondrial function. The present study increases the understanding of the effects of vitamin D at the transcriptome level.

Background: Oral microbiome has significant impact on both oral and general health. Polyols have been promoted as sugar substitutes in prevention of oral diseases. We aimed to reveal the effect of candies containing erythritol, xylitol or control (sorbitol) on salivary microbiome., Methods: Ninety children (11.3 ± 0.6 years) consumed candies during 3 years. Microbial communities were profiled using Illumina HiSeq 2000 sequencing and real-time PCR., Results: The dominant phyla in saliva were Firmicutes (39.1%), Proteobacteria (26.1%), Bacteroidetes (14.7%), Actinobacteria (12%) and Fusobacteria (6%). The microbiome of erythritol group significantly differed from that of the other groups. Both erythritol and xylitol reduced the number of observed bacterial phylotypes in comparison to the control group. The relative abundance of the genera Veillonella, Streptococcus and Fusobacterium were higher while that of Bergeyella lower after erythritol intervention when comparing with control. The lowest prevalence of caries-related mutans streptococci corresponded with the lowest clinical caries markers in the erythritol group., Conclusions: Daily consumption of erythritol, xylitol or control candies has a specific influence on the salivary microbiome composition in schoolchildren. Erythritol is associated with the lowest prevalence of caries-related mutans streptococci and the lowest levels of clinical caries experience., Trial Registration: ClinicalTrials.gov Identifier NCT01062633.

Targeted next-generation sequencing (NGS) methods have become essential in medical research and diagnostics. In addition to NGS sensitivity and high-throughput capacity, precise biomolecule counting based on unique molecular identifier (UMI) has potential to increase biomolecule detection accuracy. Although UMIs are widely used in basic research its introduction to clinical assays is still in progress. Here, we present a robust and cost-effective TAC-seq (Targeted Allele Counting by sequencing) method that uses UMIs to estimate the original molecule counts of mRNAs, microRNAs, and cell-free DNA. We applied TAC-seq in three different clinical applications and compared the results with standard NGS. RNA samples extracted from human endometrial biopsies were analyzed using previously described 57 mRNA-based receptivity biomarkers and 49 selected microRNAs at different expression levels. Cell-free DNA aneuploidy testing was based on cell line (47,XX, +21) genomic DNA. TAC-seq mRNA profiling showed identical clustering results to transcriptome RNA sequencing, and microRNA detection demonstrated significant reduction in amplification bias, allowing to determine minor expression changes between different samples that remained undetermined by standard NGS. The mimicking experiment for cell-free DNA fetal aneuploidy analysis showed that TAC-seq can be applied to count highly fragmented DNA, detecting significant ( p  = 7.6 × 10 -4 ) excess of chromosome 21 molecules at 10% fetal fraction level. Based on three proof-of-principle applications we demonstrate that TAC-seq is an accurate and highly potential biomarker profiling method for advanced medical research and diagnostics., Competing Interests: Kaarel Krjutškov, Mariann Koel, Juha Kere, and Andres Salumets declare that they have submitted an international patent application on the TAC-seq method. The other authors declare no competing interests.

Extracellular small RNAs (sRNAs), including microRNAs (miRNAs), are promising biomarkers for diseases such as Duchenne muscular dystrophy (DMD), although their biological relevance is largely unknown. To investigate the relationship between intracellular and extracellular sRNA levels on a global scale, we performed sRNA sequencing in four muscle types and serum from wild-type, dystrophic mdx, and mdx mice in which dystrophin protein expression was restored by exon skipping. Differentially abundant sRNAs were identified in serum (mapping to miRNA, small nuclear RNA [snRNA], and PIWI-interacting RNA [piRNA] loci). One novel candidate biomarker, miR-483, was increased in both mdx serum and muscle, and also elevated in DMD patient sera. Dystrophin restoration induced global shifts in miRNA (including miR-483) and snRNA-fragment abundance toward wild-type levels. Specific serum piRNA-like sRNAs also responded to exon skipping therapy. Absolute miRNA expression in muscle was positively correlated with abundance in the circulation, although multiple highly expressed miRNAs in muscle were not elevated in mdx serum, suggesting that both passive and selective release mechanisms contribute to serum miRNA levels. In conclusion, this study has revealed new insights into the sRNA biology of dystrophin deficiency and identified novel DMD biomarkers., (Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.)

Aims/hypothesis: There is a great need to identify factors that could protect pancreatic beta cells against apoptosis or stimulate their replication and thus prevent or reverse the development of diabetes. One potential candidate is mesencephalic astrocyte-derived neurotrophic factor (MANF), an endoplasmic reticulum (ER) stress inducible protein. Manf knockout mice used as a model of diabetes develop the condition because of increased apoptosis and reduced proliferation of beta cells, apparently related to ER stress. Given this novel association between MANF and beta cell death, we studied the potential of MANF to protect human beta cells against experimentally induced ER stress., Methods: Primary human islets were challenged with proinflammatory cytokines, with or without MANF. Cell viability was analysed and global transcriptomic analysis performed. Results were further validated using the human beta cell line EndoC-βH1., Results: There was increased expression and secretion of MANF in human beta cells in response to cytokines. Addition of recombinant human MANF reduced cytokine-induced cell death by 38% in human islets (p < 0.05). MANF knockdown in EndoC-βH1 cells led to increased ER stress after cytokine challenge. Mechanistic studies showed that the protective effect of MANF was associated with repression of the NF-κB signalling pathway and amelioration of ER stress. MANF also increased the proliferation of primary human beta cells twofold when TGF-β signalling was inhibited (p < 0.01)., Conclusions/interpretation: Our studies show that exogenous MANF protein can provide protection to human beta cells against death induced by inflammatory stress. The antiapoptotic and mitogenic properties of MANF make it a potential therapeutic agent for beta cell protection.

CRISPR-Cas9-based gene activation (CRISPRa) is an attractive tool for cellular reprogramming applications due to its high multiplexing capacity and direct targeting of endogenous loci. Here we present the reprogramming of primary human skin fibroblasts into induced pluripotent stem cells (iPSCs) using CRISPRa, targeting endogenous OCT4, SOX2, KLF4, MYC, and LIN28A promoters. The low basal reprogramming efficiency can be improved by an order of magnitude by additionally targeting a conserved Alu-motif enriched near genes involved in embryo genome activation (EEA-motif). This effect is mediated in part by more efficient activation of NANOG and REX1. These data demonstrate that human somatic cells can be reprogrammed into iPSCs using only CRISPRa. Furthermore, the results unravel the involvement of EEA-motif-associated mechanisms in cellular reprogramming.