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.

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.

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.

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.

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.

Acute wheeze-specific gene module shows correlation with vitamin D and asthma medication.

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 FEV1%/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.

Estrogen receptor ß, a regulator of androgen receptor signaling in the mouse ventral prostate.

As estrogen receptor ß-/- (ERß-/-) mice age, the ventral prostate (VP) develops increased numbers of hyperplastic, fibroplastic lesions and inflammatory cells. To identify genes involved in these changes, we used RNA sequencing and immunohistochemistry to compare gene expression profiles in the VP of young (2-mo-old) and aging (18-mo-old) ERß-/- mice and their WT littermates. We also treated young and old WT mice with an ERß-selective agonist and evaluated protein expression. The most significant findings were that ERß down-regulates androgen receptor (AR) signaling and up-regulates the tumor suppressor phosphatase and tensin homolog (PTEN). ERß agonist increased expression of the AR corepressor dachshund family (DACH1/2), T-cadherin, stromal caveolin-1, and nuclear PTEN and decreased expression of RAR-related orphan receptor c, Bcl2, inducible nitric oxide synthase, and IL-6. In the ERß-/- mouse VP, RNA sequencing revealed that the following genes were up-regulated more than fivefold: Bcl2, clusterin, the cytokines CXCL16 and -17, and a marker of basal/intermediate cells (prostate stem cell antigen) and cytokeratins 4, 5, and 17. The most down-regulated genes were the following: the antioxidant gene glutathione peroxidase 3; protease inhibitors WAP four-disulfide core domain 3 (WFDC3); the tumor-suppressive genes T-cadherin and caveolin-1; the regulator of transforming growth factor ß signaling SMAD7; and the PTEN ubiquitin ligase NEDD4. The role of ERß in opposing AR signaling, proliferation, and inflammation suggests that ERß-selective agonists may be used to prevent progression of prostate cancer, prevent fibrosis and development of benign prostatic hyperplasia, and treat prostatitis.

Gain-of-function CEBPE mutation causes noncanonical autoinflammatory inflammasomopathy.

BACKGROUND: CCAAT enhancer-binding protein epsilon (C/EBPε) is a transcription factor involved in late myeloid lineage differentiation and cellular function. The only previously known disorder linked to C/EBPε is autosomal recessive neutrophil-specific granule deficiency leading to severely impaired neutrophil function and early mortality. OBJECTIVE: The aim of this study was to molecularly characterize the effects of C/EBPε transcription factor Arg219His mutation identified in a Finnish family with previously genetically uncharacterized autoinflammatory and immunodeficiency syndrome. METHODS: Genetic analysis, proteomics, genome-wide transcriptional profiling by means of RNA-sequencing, chromatin immunoprecipitation (ChIP) sequencing, and assessment of the inflammasome function of primary macrophages were performed. RESULTS: Studies revealed a novel mechanism of genome-wide gain-of-function that dysregulated transcription of 464 genes. Mechanisms involved dysregulated noncanonical inflammasome activation caused by decreased association with transcriptional repressors, leading to increased chromatin occupancy and considerable changes in transcriptional activity, including increased expression of NLR family, pyrin domain-containing 3 protein (NLRP3) and constitutively expressed caspase-5 in macrophages. CONCLUSION: We describe a novel autoinflammatory disease with defective neutrophil function caused by a homozygous Arg219His mutation in the transcription factor C/EBPε. Mutated C/EBPε acts as a regulator of both the inflammasome and interferome, and the Arg219His mutation causes the first human monogenic neomorphic and noncanonical inflammasomopathy/immunodeficiency. The mechanism, including widely dysregulated transcription, is likely not unique for C/EBPε. Similar multiomics approaches should also be used in studying other transcription factor-associated diseases.

Guide for library design and bias correction for large-scale transcriptome studies using highly multiplexed RNAseq methods.

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.

The human PRD-like homeobox gene LEUTX has a central role in embryo genome activation.

Leucine twenty homeobox (LEUTX) is a paired (PRD)-like homeobox gene that is expressed almost exclusively in human embryos during preimplantation development. We previously identified a novel transcription start site for the predicted human LEUTX gene based on the transcriptional analysis of human preimplantation embryos. The novel variant encodes a protein with a complete homeodomain. Here, we provide a detailed description of the molecular cloning of the complete homeodomain-containing LEUTX Using a human embryonic stem cell overexpression model we show that the complete homeodomain isoform is functional and sufficient to activate the transcription of a large proportion of the genes that are upregulated in human embryo genome activation (EGA), whereas the previously predicted partial homeodomain isoform is largely inactive. Another PRD-like transcription factor, DPRX, is then upregulated as a powerful repressor of transcription. We propose a two-stage model of human EGA in which LEUTX acts as a transcriptional activator at the 4-cell stage, and DPRX as a balancing repressor at the 8-cell stage. We conclude that LEUTX is a candidate regulator of human EGA.

MUC20 expression marks the receptive phase of the human endometrium.

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.

Creating basis for introducing non-invasive prenatal testing in the Estonian public health setting.

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.

Impact of polyols on Oral microbiome of Estonian schoolchildren.

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.

MANF protects human pancreatic beta cells against stress-induced cell death.

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.

Mlh1 deficiency in normal mouse colon mucosa associates with chromosomally unstable colon cancer.

Colorectal cancer (CRC) genome is unstable and different types of instabilities, such as chromosomal instability (CIN) and microsatellite instability (MSI) are thought to reflect distinct cancer initiating mechanisms. Although 85% of sporadic CRC reveal CIN, 15% reveal mismatch repair (MMR) malfunction and MSI, the hallmarks of Lynch syndrome with inherited heterozygous germline mutations in MMR genes. Our study was designed to comprehensively follow genome-wide expression changes and their implications during colon tumorigenesis. We conducted a long-term feeding experiment in the mouse to address expression changes arising in histologically normal colonic mucosa as putative cancer preceding events, and the effect of inherited predisposition (Mlh1+/-) and Western-style diet (WD) on those. During the 21-month experiment, carcinomas developed mainly in WD-fed mice and were evenly distributed between genotypes. Unexpectedly, the heterozygote (B6.129-Mlh1tm1Rak) mice did not show MSI in their CRCs. Instead, both wildtype and heterozygote CRC mice showed a distinct mRNA expression profile and shortage of several chromosomal segregation gene-specific transcripts (Mlh1, Bub1, Mis18a, Tpx2, Rad9a, Pms2, Cenpe, Ncapd3, Odf2 and Dclre1b) in their colon mucosa, as well as an increased mitotic activity and abundant numbers of unbalanced/atypical mitoses in tumours. Our genome-wide expression profiling experiment demonstrates that cancer preceding changes are already seen in histologically normal colon mucosa and that decreased expressions of Mlh1 and other chromosomal segregation genes may form a field-defect in mucosa, which trigger MMR-proficient, chromosomally unstable CRC.

Endometrial receptivity revisited: endometrial transcriptome adjusted for tissue cellular heterogeneity.

STUDY QUESTION: Does cellular composition of the endometrial biopsy affect the gene expression profile of endometrial whole-tissue samples? SUMMARY ANSWER: The differences in epithelial and stromal cell proportions in endometrial biopsies modify the whole-tissue gene expression profiles and affect the results of differential expression analyses. WHAT IS ALREADY KNOWN: Each cell type has its unique gene expression profile. The proportions of epithelial and stromal cells vary in endometrial tissue during the menstrual cycle, along with individual and technical variation due to the method and tools used to obtain the tissue biopsy. STUDY DESIGN, SIZE, DURATION: Using cell-population specific transcriptome data and computational deconvolution approach, we estimated the epithelial and stromal cell proportions in whole-tissue biopsies taken during early secretory and mid-secretory phases. The estimated cellular proportions were used as covariates in whole-tissue differential gene expression analysis. Endometrial transcriptomes before and after deconvolution were compared and analysed in biological context. PARTICIPANTS/MATERIAL, SETTING, METHODS: Paired early- and mid-secretory endometrial biopsies were obtained from 35 healthy, regularly cycling, fertile volunteers, aged 23-36 years, and analysed by RNA sequencing. Differential gene expression analysis was performed using two approaches. In one of them, computational deconvolution was applied as an intermediate step to adjust for the proportions of epithelial and stromal cells in the endometrial biopsy. The results were then compared to conventional differential expression analysis. Ten paired endometrial samples were analysed with qPCR to validate the results. MAIN RESULTS AND THE ROLE OF CHANCE: The estimated average proportions of stromal and epithelial cells in early secretory phase were 65% and 35%, and during mid-secretory phase, 46% and 54%, respectively, correlating well with the results of histological evaluation (r = 0.88, P = 1.1 × 10-6). Endometrial tissue transcriptomic analysis showed that approximately 26% of transcripts (n = 946) differentially expressed in receptive endometrium in cell-type unadjusted analysis also remain differentially expressed after adjustment for biopsy cellular composition. However, the other 74% (n = 2645) become statistically non-significant after adjustment for biopsy cellular composition, underlining the impact of tissue heterogeneity on differential expression analysis. The results suggest new mechanisms involved in endometrial maturation, involving genes like LINC01320, SLC8A1 and GGTA1P, described for the first time in context of endometrial receptivity. LARGE-SCALE DATA: The RNA-seq data presented in this study is deposited in the Gene Expression Omnibus database with accession number GSE98386. LIMITATIONS REASONS FOR CAUTION: Only dominant endometrial cell types were considered in gene expression profile deconvolution; however, other less frequent endometrial cell types also contribute to the whole-tissue gene expression profile. WIDER IMPLICATIONS OF THE FINDINGS: The better understanding of molecular processes during transition from pre-receptive to receptive endometrium serves to improve the effectiveness and personalization of assisted reproduction protocols. Biopsy cellular composition should be taken into account in future endometrial 'omics' studies, where tissue heterogeneity could potentially influence the results. STUDY FUNDING/COMPETING INTEREST(S): This study was funded by: Estonian Ministry of Education and Research (grant IUT34-16); Enterprise Estonia (EU48695); the EU-FP7 Eurostars program (NOTED, EU41564); the EU-FP7 Marie Curie Industry-Academia Partnerships and Pathways (SARM, EU324509); Horizon 2020 innovation program (WIDENLIFE, EU692065); MSCA-RISE-2015 project MOMENDO (No 691058) and the Miguel Servet Program Type I of Instituto de Salud Carlos III (CP13/00038); Spanish Ministry of Economy, Industry and Competitiveness (MINECO) and European Regional Development Fund (FEDER): grants RYC-2016-21199 and ENDORE SAF2017-87526. Authors confirm no competing interests.

Damaging heterozygous mutations in NFKB1 lead to diverse immunologic phenotypes.

BACKGROUND: The nuclear factor κ light-chain enhancer of activated B cells (NF-κB) signaling pathway is a key regulator of immune responses. Accordingly, mutations in several NF-κB pathway genes cause immunodeficiency. OBJECTIVE: We sought to identify the cause of disease in 3 unrelated Finnish kindreds with variable symptoms of immunodeficiency and autoinflammation. METHODS: We applied genetic linkage analysis and next-generation sequencing and functional analyses of NFKB1 and its mutated alleles. RESULTS: In all affected subjects we detected novel heterozygous variants in NFKB1, encoding for p50/p105. Symptoms in variant carriers differed depending on the mutation. Patients harboring a p.I553M variant presented with antibody deficiency, infection susceptibility, and multiorgan autoimmunity. Patients with a p.H67R substitution had antibody deficiency and experienced autoinflammatory episodes, including aphthae, gastrointestinal disease, febrile attacks, and small-vessel vasculitis characteristic of Behçet disease. Patients with a p.R157X stop-gain experienced hyperinflammatory responses to surgery and showed enhanced inflammasome activation. In functional analyses the p.R157X variant caused proteasome-dependent degradation of both the truncated and wild-type proteins, leading to a dramatic loss of p50/p105. The p.H67R variant reduced nuclear entry of p50 and showed decreased transcriptional activity in luciferase reporter assays. The p.I553M mutation in turn showed no change in p50 function but exhibited reduced p105 phosphorylation and stability. Affinity purification mass spectrometry also demonstrated that both missense variants led to altered protein-protein interactions. CONCLUSION: Our findings broaden the scope of phenotypes caused by mutations in NFKB1 and suggest that a subset of autoinflammatory diseases, such as Behçet disease, can be caused by rare monogenic variants in genes of the NF-κB pathway.

Optimizing bone morphogenic protein 4-mediated human embryonic stem cell differentiation into trophoblast-like cells using fibroblast growth factor 2 and transforming growth factor-ß/activin/nodal signalling inhibition.

Several studies have demonstrated that human embryonic stem cells (hESC) can be differentiated into trophoblast-like cells if exposed to bone morphogenic protein 4 (BMP4) and/or inhibitors of fibroblast growth factor 2 (FGF2) and the transforming growth factor beta (TGF-ß)/activin/nodal signalling pathways. The goal of this study was to investigate how the inhibitors of these pathways improve the efficiency of hESC differentiation when compared with basic BMP4 treatment. RNA sequencing was used to analyse the effects of all possible inhibitor combinations on the differentiation of hESC into trophoblast-like cells over 12 days. Genes differentially expressed compared with untreated cells were identified at seven time points. Additionally, expression of total human chorionic gonadotrophin (HCG) and its hyperglycosylated form (HCG-H) were determined by immunoassay from cell culture media. We showed that FGF2 inhibition with BMP4 activation up-regulates syncytiotrophoblast-specific genes (CGA, CGB and LGALS16), induces several molecular pathways involved in embryo implantation and triggers HCG-H production. In contrast, inhibition of the TGF-ß/activin/nodal pathway decreases the ability of hESC to form trophoblast-like cells. Information about the conditions needed for hESC differentiation toward trophoblast-like cells helps us to find an optimal model for studying the early development of human trophoblasts in normal and in complicated pregnancy.

The pruritus- and TH2-associated cytokine IL-31 promotes growth of sensory nerves.

BACKGROUND: Pruritus is a cardinal symptom of atopic dermatitis, and an increased cutaneous sensory network is thought to contribute to pruritus. Although the immune cell-IL-31-neuron axis has been implicated in severe pruritus during atopic skin inflammation, IL-31's neuropoietic potential remains elusive. OBJECTIVE: We sought to analyze the IL-31-related transcriptome in sensory neurons and to investigate whether IL-31 promotes sensory nerve fiber outgrowth. METHODS: In vitro primary sensory neuron culture systems were subjected to whole-transcriptome sequencing, ingenuity pathway analysis, immunofluorescence, and nerve elongation, as well as branching assays after IL-31 stimulation. In vivo we investigated the cutaneous sensory neuronal network in wild-type, Il31-transgenic, and IL-31 pump-equipped mice. RESULTS: Transgenic Il31 overexpression and subcutaneously delivered IL-31 induced an increase in the cutaneous nerve fiber density in lesional skin in vivo. Transcriptional profiling of IL-31-activated dorsal root ganglia neurons revealed enrichment for genes promoting nervous system development and neuronal outgrowth and negatively regulating cell death. Moreover, the growth cones of primary small-diameter dorsal root ganglia neurons showed abundant IL-31 receptor α expression. Indeed, IL-31 selectively promoted nerve fiber extension only in small-diameter neurons. Signal transducer and activator of transcription 3 phosphorylation mediated IL-31-induced neuronal outgrowth, and pharmacologic inhibition of signal transducer and activator of transcription 3 completely abolished this effect. In contrast, transient receptor potential cation channel vanilloid subtype 1 channels were dispensable for IL-31-induced neuronal sprouting. CONCLUSIONS: The pruritus- and TH2-associated novel cytokine IL-31 induces a distinct transcriptional program in sensory neurons, leading to nerve elongation and branching both in vitro and in vivo. This finding might help us understand the clinical observation that patients with atopic dermatitis experience increased sensitivity to minimal stimuli inducing sustained itch.