Integrative Genomic and Transcriptomic Profiling Reveals a Differential Molecular Signature in Uterine Leiomyoma versus Leiomyosarcoma.

The absence of standardized molecular profiling to differentiate uterine leiomyosarcomas versus leiomyomas represents a current diagnostic challenge. In this study, we aimed to search for a differential molecular signature for these myometrial tumors based on artificial intelligence. For this purpose, differential exome and transcriptome-wide research was performed on histologically confirmed leiomyomas (n = 52) and leiomyosarcomas (n = 44) to elucidate differences between and within these two entities. We identified a significantly higher tumor mutation burden in leiomyosarcomas vs. leiomyomas in terms of somatic single-nucleotide variants (171,863 vs. 81,152), indels (9491 vs. 4098), and copy number variants (8390 vs. 5376). Further, we discovered alterations in specific copy number variant regions that affect the expression of some tumor suppressor genes. A transcriptomic analysis revealed 489 differentially expressed genes between these two conditions, as well as structural rearrangements targeting ATRX and RAD51B. These results allowed us to develop a machine learning approach based on 19 differentially expressed genes that differentiate both tumor types with high sensitivity and specificity. Our findings provide a novel molecular signature for the diagnosis of leiomyoma and leiomyosarcoma, which could be helpful to complement the current morphological and immunohistochemical diagnosis and may lay the foundation for the future evaluation of malignancy risk.

Disrupted PGR-B and ESR1 signaling underlies defective decidualization linked to severe preeclampsia.

Background: Decidualization of the uterine mucosa drives the maternal adaptation to invasion by the placenta. Appropriate depth of placental invasion is needed to support a healthy pregnancy; shallow invasion is associated with the development of severe preeclampsia (sPE). Maternal contribution to sPE through failed decidualization is an important determinant of placental phenotype. However, the molecular mechanism underlying the in vivo defect linking decidualization to sPE is unknown. Methods: Global RNA sequencing was applied to obtain the transcriptomic profile of endometrial biopsies collected from nonpregnant women who suffer sPE in a previous pregnancy and women who did not develop this condition. Samples were randomized in two cohorts, the training and the test set, to identify the fingerprinting encoding defective decidualization in sPE and its subsequent validation. Gene Ontology enrichment and an interaction network were performed to deepen in pathways impaired by genetic dysregulation in sPE. Finally, the main modulators of decidualization, estrogen receptor 1 (ESR1) and progesterone receptor B (PGR-B), were assessed at the level of gene expression and protein abundance. Results: Here, we discover the footprint encoding this decidualization defect comprising 120 genes-using global gene expression profiling in decidua from women who developed sPE in a previous pregnancy. This signature allowed us to effectively segregate samples into sPE and control groups. ESR1 and PGR were highly interconnected with the dynamic network of the defective decidualization fingerprint. ESR1 and PGR-B gene expression and protein abundance were remarkably disrupted in sPE. Conclusions: Thus, the transcriptomic signature of impaired decidualization implicates dysregulated hormonal signaling in the decidual endometria in women who developed sPE. These findings reveal a potential footprint that could be leveraged for a preconception or early prenatal screening of sPE risk, thus improving prevention and early treatments. Funding: This work has been supported by the grant PI19/01659 (MCIU/AEI/FEDER, UE) from the Spanish Carlos III Institute awarded to TGG. NCM was supported by the PhD program FDGENT/2019/008 from the Spanish Generalitat Valenciana. IMB was supported by the PhD program PRE2019-090770 and funding was provided by the grant RTI2018-094946-B-100 (MCIU/AEI/FEDER, UE) from the Spanish Ministry of Science and Innovation with CS as principal investigator. This research was funded partially by Igenomix S.L.

Identification and Characterization of Extracellular Vesicles and Its DNA Cargo Secreted During Murine Embryo Development.

Extracellular vesicles (EVs) are known to transport DNA, but their implications in embryonic implantation are unknown. The aim of this study was to investigate EVs production and secretion by preimplantation embryos and assess their DNA cargo. Murine oocytes and embryos were obtained from six- to eight-week-old females, cultured until E4.5 and analyzed using transmission electron microscopy to examine EVs production. EVs were isolated from E4.5-day conditioned media and quantified by nanoparticle tracking analysis, characterized by immunogold, and their DNA cargo sequenced. Multivesicular bodies were observed in murine oocytes and preimplantation embryos together with the secretion of EVs to the blastocoel cavity and blastocyst spent medium. Embryo-derived EVs showed variable electron-densities and sizes (20-500 nm) and total concentrations of 1.74 × 107 ± 2.60 × 106 particles/mL. Embryo secreted EVs were positive for CD63 and ARF6. DNA cargo sequencing demonstrated no differences in DNA between apoptotic bodies or smaller EVs, although they showed significant gene enrichment compared to control medium. The analysis of sequences uniquely mapping the murine genome revealed that DNA contained in EVs showed higher representation of embryo genome than vesicle-free DNA. Murine blastocysts secrete EVs containing genome-wide sequences of DNA to the medium, reinforcing the relevance of studying these vesicles and their cargo in the preimplantation moment, where secreted DNA may help the assessment of the embryo previous to implantation.

Optimizing clinical exome design and parallel gene-testing for recessive genetic conditions in preconception carrier screening: Translational research genomic data from 14,125 exomes.

Limited translational genomic research data have been reported on the application of exome sequencing and parallel gene testing for preconception carrier screening (PCS). Here, we present individual-level data from a large PCS program in which exome sequencing was routinely performed on either gamete donors (5,845) or infertile patients (8,280) undergoing in vitro fertilization (IVF) treatment without any known family history of inheritable genetic conditions. Individual-level data on pathogenic variants were used to define conditions for PCS based on criteria for severity, penetrance, inheritance pattern, and age of onset. Fetal risk was defined based on actual carrier frequency data accounting for the specific inheritance pattern (fetal disease risk, FDR). In addition, large-scale application of exome sequencing for PCS allowed a deep investigation of the incidence of medically actionable secondary findings in this population. Exome sequencing achieved remarkable clinical sensitivity for reproductive risk of highly penetrant childhood-onset disorders (1/337 conceptions) through analysis of 114 selected gene-condition pairs. A significant contribution to fetal disease risk was observed for rare (carrier rate < 1:100) and X-linked conditions (16.7% and 41.2% of total FDR, respectively). Subgroup analysis of 776 IVF couples identified 37 at increased reproductive risk (4.8%; 95% CI = 3.4-6.5). Further, two additional couples had increased risk for very rare conditions when both members of a parental pair were treated as a unit and the search was extended to the entire exome. About 2.3% of participants showed at least one pathogenic variant for genes included in the updated American College of Medical Genetics and Genomics v2.0 list of secondary findings. Gamete donors and IVF couples showed similar carrier burden for both carrier screening and secondary findings, indicating no causal relationship to fertility. These translational research data will facilitate development of more effective PCS strategies that maximize clinical sensitivity with minimal counterproductive effects.

The differential diagnoses of uterine leiomyomas and leiomyosarcomas using DNA and RNA sequencing.

BACKGROUND: Although uterine leiomyomas and leiomyosarcomas are considered biologically unrelated tumors, they share morphologic and histologic characteristics that complicate their differential diagnosis. The long-term therapeutic option for leiomyoma is laparoscopic myomectomy with morcellation, particularly for patients who wish to preserve their fertility. However, because of the potential dissemination of undiagnosed or hidden leiomyosarcoma from morcellation, there is a need to develop a preoperative assessment of malignancy risk. OBJECTIVE: Through an integrated comparative genomic and transcriptomic analysis, we aim to identify differential genetic targets in leiomyomas vs leiomyosarcomas using next-generation sequencing as the first step toward preoperative differential diagnosis. STUDY DESIGN: Targeted sequencing of DNA and RNA coding regions for solid tumor-associated genes was performed on formalin-fixed paraffin-embedded samples from 13 leiomyomas and 13 leiomyosarcoma cases. DNA sequencing was used to identify copy number variations, single-nucleotide variants, and small insertions/deletions. RNA sequencing was used to identify gene fusions, splice variants, and/or differential gene expression profiles. RESULTS: In leiomyosarcomas, tumor mutation burden was higher in terms of copy number variations, single nucleotide variants, small insertions/deletions, and gene fusions compared with leiomyomas. For copy number variations, 20 genes were affected by deletions in leiomyosarcomas, compared with 6 observed losses in leiomyomas. Gains (duplications) were identified in 19 genes in leiomyosarcomas, but only 3 genes in leiomyomas. The most common mutations (single-nucleotide variants and insertions/deletions) for leiomyosarcomas were identified in 105 genes of all analyzed leiomyosarcomas; 82 genes were affected in leiomyomas. Of note, 1 tumor previously diagnosed as leiomyosarcoma was established as inflammatory myofibroblastic tumor along this study with a novel ALK-TNS1 fusion. Finally, a differential transcriptomic profile was observed for 11 of 55 genes analyzed in leiomyosarcomas; 8.5% of initially diagnosed leiomyosarcomas showed high-confidence, novel gene fusions that were associated with these tumors. CONCLUSION: Through integrated comparative genomic and transcriptomic analyses, we identified novel differential genetic targets that potentially differentiate leiomyosarcomas and leiomyomas. This provides a new insight into the differential diagnosis of these myometrial tumors.

Evidence that the endometrial microbiota has an effect on implantation success or failure.

BACKGROUND: Bacterial cells in the human body account for 1-3% of total body weight and are at least equal in number to human cells. Recent research has focused on understanding how the different bacterial communities in the body (eg, gut, respiratory, skin, and vaginal microbiomes) predispose to health and disease. The microbiota of the reproductive tract has been inferred from the vaginal bacterial communities, and the uterus has been classically considered a sterile cavity. However, while the vaginal microbiota has been investigated in depth, there is a paucity of consistent data regarding the existence of an endometrial microbiota and its possible impact in reproductive function. OBJECTIVE: This study sought to test the existence of an endometrial microbiota that differs from that in the vagina, assess its hormonal regulation, and analyze the impact of the endometrial microbial community on reproductive outcome in infertile patients undergoing in vitro fertilization. STUDY DESIGN: To identify the existence of an endometrial microbiota, paired samples of endometrial fluid and vaginal aspirates were obtained simultaneously from 13 fertile women in prereceptive and receptive phases within the same menstrual cycle (total samples analyzed n = 52). To investigate the hormonal regulation of the endometrial microbiota during the acquisition of endometrial receptivity, endometrial fluid was collected at prereceptive and receptive phases within the same cycle from 22 fertile women (n = 44). Finally, the reproductive impact of an altered endometrial microbiota in endometrial fluid was assessed by implantation, ongoing pregnancy, and live birth rates in 35 infertile patients undergoing in vitro fertilization (total samples n = 41) with a receptive endometrium diagnosed using the endometrial receptivity array. Genomic DNA was obtained either from endometrial fluid or vaginal aspirate and sequenced by 454 pyrosequencing of the V3-V5 region of the 16S ribosomal RNA (rRNA) gene; the resulting sequences were taxonomically assigned using QIIME. Data analysis was performed using R packages. The χ2 test, Student t test, and analysis of variance were used for statistical analyses. RESULTS: When bacterial communities from paired endometrial fluid and vaginal aspirate samples within the same subjects were interrogated, different bacterial communities were detected between the uterine cavity and the vagina of some subjects. Based on its composition, the microbiota in the endometrial fluid, comprising up to 191 operational taxonomic units, was defined as a Lactobacillus-dominated microbiota (>90% Lactobacillus spp.) or a non-Lactobacillus-dominated microbiota (<90% Lactobacillus spp. with >10% of other bacteria). Although the endometrial microbiota was not hormonally regulated during the acquisition of endometrial receptivity, the presence of a non-Lactobacillus-dominated microbiota in a receptive endometrium was associated with significant decreases in implantation [60.7% vs 23.1% (P = .02)], pregnancy [70.6% vs 33.3% (P = .03)], ongoing pregnancy [58.8% vs 13.3% (P = .02)], and live birth [58.8% vs 6.7% (P = .002)] rates. CONCLUSION: Our results demonstrate the existence of an endometrial microbiota that is highly stable during the acquisition of endometrial receptivity. However, pathological modification of its profile is associated with poor reproductive outcomes for in vitro fertilization patients. This finding adds a novel microbiological dimension to the reproductive process.

The Mutational Landscape of Acute Promyelocytic Leukemia Reveals an Interacting Network of Co-Occurrences and Recurrent Mutations.

Preliminary Acute Promyelocytic Leukemia (APL) whole exome sequencing (WES) studies have identified a huge number of somatic mutations affecting more than a hundred different genes mainly in a non-recurrent manner, suggesting that APL is a heterogeneous disease with secondary relevant changes not yet defined. To extend our knowledge of subtle genetic alterations involved in APL that might cooperate with PML/RARA in the leukemogenic process, we performed a comprehensive analysis of somatic mutations in APL combining WES with sequencing of a custom panel of targeted genes by next-generation sequencing. To select a reduced subset of high confidence candidate driver genes, further in silico analysis were carried out. After prioritization and network analysis we found recurrent deleterious mutations in 8 individual genes (STAG2, U2AF1, SMC1A, USP9X, IKZF1, LYN, MYCBP2 and PTPN11) with a strong potential of being involved in APL pathogenesis. Our network analysis of multiple mutations provides a reliable approach to prioritize genes for additional analysis, improving our knowledge of the leukemogenesis interactome. Additionally, we have defined a functional module in the interactome of APL. The hypothesis is that the number, or the specific combinations, of mutations harbored in each patient might not be as important as the disturbance caused in biological key functions, triggered by several not necessarily recurrent mutations.

Comprehensive carrier genetic test using next-generation deoxyribonucleic acid sequencing in infertile couples wishing to conceive through assisted reproductive technology.

OBJECTIVE: To develop an expanded pan-ethnic preconception carrier genetic screening test for use in assisted reproductive technology (ART) patients and donors. DESIGN: Retrospective analysis of results obtained from 2,570 analyses. SETTING: Reproductive genetic laboratory. PATIENT(S): The 2,570 samples comprised 1,170 individuals from the gamete donor programs; 1,124 individuals corresponding to the partner of the patient receiving the donated gamete; and 276 individuals from 138 couples seeking ART using their own gametes. INTERVENTION(S): None. MAIN OUTCOME MEASURE(S): Next-generation sequencing of 549 recessive and X-linked genes involved in severe childhood phenotypes reinforced with five complementary tests covering high prevalent mutations not detected by next-generation sequencing. RESULT(S): Preclinical validation included 48 DNA samples carrying known mutations for 27 genes, resulting in a sensitivity of 99%. In the clinical dataset, 2,161 samples (84%) tested positive, with an average carrier burden of 2.3 per sample. Five percent of the couples using their own gametes were found to have pathogenic variants conferring high risk for six different diseases. These high-risk couples and patients received genetic counseling and recommendations for preimplantation genetic diagnosis. For patients receiving gamete donation, we applied a genetic testing and blinded matching system to avoid high-risk combinations regardless of their carrier burden. For female donors, 1.94% were positive for X-linked conditions; they received genetic counselling and were discarded. CONCLUSION(S): We have developed a comprehensive carrier genetic screening test that, combined with our matching system and genetic counseling, constitutes a powerful tool to avoid more than 600 mendelian diseases in the offspring of patients undergoing ART.

The role of the interactome in the maintenance of deleterious variability in human populations.

Recent genomic projects have revealed the existence of an unexpectedly large amount of deleterious variability in the human genome. Several hypotheses have been proposed to explain such an apparently high mutational load. However, the mechanisms by which deleterious mutations in some genes cause a pathological effect but are apparently innocuous in other genes remain largely unknown. This study searched for deleterious variants in the 1,000 genomes populations, as well as in a newly sequenced population of 252 healthy Spanish individuals. In addition, variants causative of monogenic diseases and somatic variants from 41 chronic lymphocytic leukaemia patients were analysed. The deleterious variants found were analysed in the context of the interactome to understand the role of network topology in the maintenance of the observed mutational load. Our results suggest that one of the mechanisms whereby the effect of these deleterious variants on the phenotype is suppressed could be related to the configuration of the protein interaction network. Most of the deleterious variants observed in healthy individuals are concentrated in peripheral regions of the interactome, in combinations that preserve their connectivity, and have a marginal effect on interactome integrity. On the contrary, likely pathogenic cancer somatic deleterious variants tend to occur in internal regions of the interactome, often with associated structural consequences. Finally, variants causative of monogenic diseases seem to occupy an intermediate position. Our observations suggest that the real pathological potential of a variant might be more a systems property rather than an intrinsic property of individual proteins.

Exome sequencing identifies a new mutation in SERAC1 in a patient with 3-methylglutaconic aciduria.

3-Methylglutaconic aciduria (3-MGA-uria) is a heterogeneous group of syndromes characterized by an increased excretion of 3-methylglutaconic and 3-methylglutaric acids. Five types of 3-MGA-uria (I to V) with different clinical presentations have been described. Causative mutations in TAZ, OPA3, DNAJC19, ATP12, ATP5E, and TMEM70 have been identified. After excluding the known genetic causes of 3-MGA-uria we used exome sequencing to investigate a patient with Leigh syndrome and 3-MGA-uria. We identified a homozygous variant in SERAC1 (c.202C>T; p.Arg68*), that generates a premature stop codon at position 68 of SERAC1 protein. Western blot analysis in patient's fibroblasts showed a complete absence of SERAC1 that was consistent with the prediction of a truncated protein and supports the pathogenic role of the mutation. During the course of this project a parallel study identified mutations in SERAC1 as the genetic cause of the disease in 15 patients with MEGDEL syndrome, which was compatible with the clinical and biochemical phenotypes of the patient described here. In addition, our patient developed microcephaly and optic atrophy, two features not previously reported in MEGDEL syndrome. We highlight the usefulness of exome sequencing to reveal the genetic bases of human rare diseases even if only one affected individual is available.

A map of human microRNA variation uncovers unexpectedly high levels of variability.

BACKGROUND: MicroRNAs (miRNAs) are key components of the gene regulatory network in many species. During the past few years, these regulatory elements have been shown to be involved in an increasing number and range of diseases. Consequently, the compilation of a comprehensive map of natural variability in a healthy population seems an obvious requirement for future research on miRNA-related pathologies. METHODS: Data on 14 populations from the 1000 Genomes Project were analyzed, along with new data extracted from 60 exomes of healthy individuals from a population from southern Spain, sequenced in the context of the Medical Genome Project, to derive an accurate map of miRNA variability. RESULTS: Despite the common belief that miRNAs are highly conserved elements, analysis of the sequences of the 1,152 individuals indicated that the observed level of variability is double what was expected. A total of 527 variants were found. Among these, 45 variants affected the recognition region of the corresponding miRNA and were found in 43 different miRNAs, 26 of which are known to be involved in 57 diseases. Different parts of the mature structure of the miRNA were affected to different degrees by variants, which suggests the existence of a selective pressure related to the relative functional impact of the change. Moreover, 41 variants showed a significant deviation from the Hardy-Weinberg equilibrium, which supports the existence of a selective process against some alleles. The average number of variants per individual in miRNAs was 28. CONCLUSIONS: Despite an expectation that miRNAs would be highly conserved genomic elements, our study reports a level of variability comparable to that observed for coding genes.

Urinary 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxo-dG), a reliable oxidative stress marker in hypertension.

UNLABELLED: The potential use of oxidative stress products as disease markers and progression is an important aspect of biomedical research. In the present study, the quantification of urine 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxo-dG) concentration has been used to express the oxidation status of hypertensive subjects. 8-oxo-dG has been simultaneously isolated and assayed in nuclear (nDNA) and mitochondrial DNA (mtDNA). In addition, oxidative stress of mononuclear cells has been estimated by means of GSH and GSSG levels and GSSG/GSH ratio in hypertensive subjects before and after antihypertensive treatment. It is shown that oxidative stress decreases significantly in hypertensive patients after treatment the effect being accompanied by reduction of their blood pressure. A significant correlation is observed comparing the yield of urine 8-oxo-dG and that isolated from mitochondria DNA. Moreover, urinary excretion of 8-oxo-dG also correlates with the GSSG/GSH ratio of cells. CONCLUSION: urine 8-oxo-dG assay is a good marker for monitoring oxidative stress changes in hypertensives.