Clonal origin and genomic diversity in Lynch syndrome-associated endometrial cancer with multiple synchronous tumors: Identification of the pathogenicity of MLH1 p.L582H.

Lynch syndrome-associated endometrial cancer patients often present multiple synchronous tumors and this assessment can affect treatment strategies. We present a case of a 27-year-old woman with tumors in the uterine corpus, cervix, and ovaries who was diagnosed with endometrial cancer and exhibited cervical invasion and ovarian metastasis. Her family history suggested Lynch syndrome, and genetic testing identified a variant of uncertain significance, MLH1 p.L582H. We conducted immunohistochemical staining, microsatellite instability analysis, and Sanger sequencing for Lynch syndrome-associated cancers in three generations of the family and identified consistent MLH1 loss. Whole-exome sequencing for the corpus, cervical, and ovarian tumors of the proband identified a copy-neutral loss of heterozygosity (LOH) occurring at the MLH1 position in all tumors. This indicated that the germline variant and the copy-neutral LOH led to biallelic loss of MLH1 and was the cause of cancer initiation. All tumors shared a portion of somatic mutations with high mutant allele frequencies, suggesting a common clonal origin. There were no mutations shared only between the cervix and ovary samples. The profiles of mutant allele frequencies shared between the corpus and cervix or ovary indicated that two different subclones originating from the corpus independently metastasized to the cervix or ovary. Additionally, all tumors presented unique mutations in endometrial cancer-associated genes such as ARID1A and PIK3CA. In conclusion, we demonstrated clonal origin and genomic diversity in a Lynch syndrome-associated endometrial cancer, suggesting the importance of evaluating multiple sites in Lynch syndrome patients with synchronous tumors.

Metagenomic analyses of 7000 to 5500 years old coprolites excavated from the Torihama shell-mound site in the Japanese archipelago.

Coprolites contain various kinds of ancient DNAs derived from gut micro-organisms, viruses, and foods, which can help to determine the gut environment of ancient peoples. Their genomic information should be helpful in elucidating the interaction between hosts and microbes for thousands of years, as well as characterizing the dietary behaviors of ancient people. We performed shotgun metagenomic sequencing on four coprolites excavated from the Torihama shell-mound site in the Japanese archipelago. The coprolites were found in the layers of the Early Jomon period, corresponding stratigraphically to 7000 to 5500 years ago. After shotgun sequencing, we found that a significant number of reads showed homology with known gut microbe, viruses, and food genomes typically found in the feces of modern humans. We detected reads derived from several types of phages and their host bacteria simultaneously, suggesting the coexistence of viruses and their hosts. The food genomes provide biological evidence for the dietary behavior of the Jomon people, consistent with previous archaeological findings. These results indicate that ancient genomic analysis of coprolites is useful for understanding the gut environment and lifestyle of ancient peoples.

SLFN11 is a BRCA Independent Biomarker for the Response to Platinum-Based Chemotherapy in High-Grade Serous Ovarian Cancer and Clear Cell Ovarian Carcinoma.

BRCA1/2 mutations are robust biomarkers for platinum-based chemotherapy in epithelial ovarian cancers. However, BRCA1/2 mutations in clear cell ovarian carcinoma (CCC) are less frequent compared with high-grade serous ovarian cancer (HGSC). The discovery of biomarkers that can be applied to CCC is an unmet need in chemotherapy. Schlafen 11 (SLFN11) has attracted attention as a novel sensitizer for DNA-damaging agents including platinum. In this study, we investigated the utility of SLFN11 in HGSC and CCC for platinum-based chemotherapy. SLFN11 expression was analyzed retrospectively by IHC across 326 ovarian cancer samples. The clinicopathologic significance of SLFN11 expression was analyzed across 57 advanced HGSC as a discovery set, 96 advanced HGSC as a validation set, and 57 advanced CCC cases, all of whom received platinum-based chemotherapy. BRCA1/2 mutation was analyzed using targeted-gene sequencing. In the HGSC cohort, the SLFN11-positive and BRCA mutation group showed significantly longer whereas the SLFN11-negative and BRCA wild-type group showed significantly shorter progression-free survival and overall survival. Moreover, SLFN11-positive HGSC shrunk significantly better than SLFN11-negative HGSC after neoadjuvant chemotherapy. Comparable results were obtained with CCC but without consideration of BRCA1/2 mutation due to a small population. Multivariate analysis identified SLFN11 as an independent factor for better survival in HGSC and CCC. The SLFN11-dependent sensitivity to platinum and PARP inhibitors were validated with genetically modified non-HGSC ovarian cancer cell lines. Our study reveals that SLFN11 predicts platinum sensitivity in HGSC and CCC independently of BRCA1/2 mutation status, indicating that SLFN11 assessment can guide treatment selection in HGSC and CCC.

Detection of APP gene recombinant in human blood plasma.

The pathogenesis of Alzheimer's disease (AD) is believed to involve the accumulation of amyloid-ß in the brain, which is produced by the sequential cleavage of amyloid precursor protein (APP) by ß-secretase and γ-secretase. Recently, analysis of genomic DNA and mRNA from postmortem brain neurons has revealed intra-exonic recombinants of APP (gencDNA), which have been implicated in the accumulation of amyloid-ß. In this study, we computationally analyzed publicly available sequence data (SRA) using probe sequences we constructed to screen APP gencDNAs. APP gencDNAs were detected in SRAs constructed from both genomic DNA and RNA obtained from the postmortem brain and in the SRA constructed from plasma cell-free mRNA (cf-mRNA). The SRA constructed from plasma cf-mRNA showed a significant difference in the number of APP gencDNA reads between SAD and NCI: the p-value from the Mann-Whitney U test was 5.14 × 10-6. The transcripts were also found in circulating nucleic acids (CNA) from our plasma samples with NGS analysis. These data indicate that transcripts of APP gencDNA can be detected in blood plasma and suggest the possibility of using them as blood biomarkers for Alzheimer's disease.

Spatial genomic diversity associated with APOBEC mutagenesis in squamous cell carcinoma arising from ovarian teratoma.

Although the gross and microscopic features of squamous cell carcinoma arising from ovarian mature cystic teratoma (MCT-SCC) vary from case to case, the spatial spreading of genomic alterations within the tumor remains unclear. To clarify the spatial genomic diversity in MCT-SCCs, we performed whole-exome sequencing by collecting 16 samples from histologically different parts of two MCT-SCCs. Both cases showed histological diversity within the tumors (case 1: nonkeratinizing and keratinizing SCC and case 2: nonkeratinizing SCC and anaplastic carcinoma) and had different somatic mutation profiles by histological findings. Mutation signature analysis revealed a significantly enriched apolipoprotein B mRNA editing enzyme catalytic subunit (APOBEC) signature at all sites. Intriguingly, the spread of genomic alterations within the tumor and the clonal evolution patterns from nonmalignant epithelium to cancer sites differed between cases. TP53 mutation and copy number alterations were widespread at all sites, including the nonmalignant epithelium, in case 1. Keratinizing and nonkeratinizing SCCs were differentiated by the occurrence of unique somatic mutations from a common ancestral clone. In contrast, the nonmalignant epithelium showed almost no somatic mutations in case 2. TP53 mutation and the copy number alteration similarities were observed only in nonkeratinizing SCC samples. Nonkeratinizing SCC and anaplastic carcinoma shared almost no somatic mutations, suggesting that each locally and independently arose in the MCT. We demonstrated that two MCT-SCCs with different histologic findings were highly heterogeneous tumors with clearly different clones associated with APOBEC-mediated mutagenesis, suggesting the importance of evaluating intratumor histological and genetic heterogeneity among multiple sites of MCT-SCC.

AMBRA1 p.Gln30Arg Mutation, Identified in a Cowden Syndrome Family, Exhibits Hyperproliferative Potential in hTERT-RPE1 Cells.

Cowden syndrome (CS) is a rare autosomal dominant disorder associated with multiple hamartomatous and neoplastic lesions in various organs. Most CS patients have been found to have germline mutations in the PTEN tumor suppressor. In the present study, we investigated the causative gene of CS in a family of PTEN (phosphatase and tensin homolog deleted on chromosome 10) -negative CS patients. Whole exome sequencing analysis revealed AMBRA1 (Autophagy and Beclin 1 Regulator 1) as a novel candidate gene harboring two germline variants: p.Gln30Arg (Q30R) and p.Arg1195Ser (R1195S). AMBRA1 is a key regulator of the autophagy signaling network and a tumor suppressor. To functionally validate the role of AMBRA1 in the clinical manifestations of CS, we generated AMBRA1 depletion and Q30R mutation in hTERT-RPE1 (humanTelomerase Reverse Transcriptase-immortalized Retinal Pigmented Epithelial cells) using the CRISPR-Cas9 gene editing system. We observed that both AMBRA1-depleted and mutant cells showed accumulation in the S phase, leading to hyperproliferation, which is a characteristic of hamartomatous lesions. Specifically, the AMBRA1 Q30R mutation disturbed the G1/S transition of cells, leading to continuous mitotic entry of mutant cells, irrespective of the extracellular condition. From our analysis of primary ciliogenesis in these cells, we speculated that the mitotic entry of AMBRA1 Q30R mutants could be due to non-functional primary cilia that lead to impaired processing of extracellular sensory signals. Additionally, we observed a situs inversus phenotype in ambra1-depleted zebrafish, a developmental abnormality resulting from dysregulated primary ciliogenesis. Taken together, we established that the AMBRA1 Q30R mutation that we observed in CS patients might play an important role in inducing the hyperproliferative potential of cells through regulating primary ciliogenesis.

Spontaneous activity in whisker-innervating region of neonatal mouse trigeminal ganglion.

Spontaneous activity during the early postnatal period is thought to be crucial for the establishment of mature neural circuits. It remains unclear if the peripheral structure of the developing somatosensory system exhibits spontaneous activity, similar to that observed in the retina and cochlea of developing mammals. By establishing an ex vivo calcium imaging system, here we found that neurons in the whisker-innervating region of the trigeminal ganglion (TG) of neonatal mice generate spontaneous activity. A small percentage of neurons showed some obvious correlated activity, and these neurons were mostly located close to one another. TG spontaneous activity was majorly exhibited by medium-to-large diameter neurons, a characteristic of mechanosensory neurons, and was blocked by chelation of extracellular calcium. Moreover, this activity was diminished by the adult stage. Spontaneous activity in the TG during the first postnatal week could be a source of spontaneous activity observed in the neonatal mouse barrel cortex.

The landscape of genetic alterations of UVB-induced skin tumors in DNA repair-deficient mice.

Non-melanoma skin cancer (NMSC) is mainly caused by ultraviolet (UV)-induced somatic mutations and is characterized by UV signature modifications. Xeroderma pigmentosum group A (Xpa) knockout mice exhibit extreme UV-induced photo-skin carcinogenesis, along with a photosensitive phenotype. We performed whole-exome sequencing (WES) of squamous cell carcinoma (SCC) samples after repetitive ultraviolet B (UVB) exposure to investigate the differences in the landscape of somatic mutations between Xpa knockout and wild-type mice. Although the tumors that developed in mice harboured UV signature mutations in a similar set of cancer-related genes, the pattern of transcriptional strand asymmetry was largely different; UV signature mutations in Xpa knockout and wild-type mice preferentially occurred in transcribed and non-transcribed strands, respectively, reflecting a deficiency in transcription-coupled nucleotide excision repair in Xpa knockout mice. Serial time point analyses of WES for a tumor induced by only a single UVB exposure showed pathogenic mutations in Kras, Fat1, and Kmt2c, which may be driver genes for the initiation and promotion of SCC in Xpa knockout mice. Furthermore, the inhibitory effects on tumor production in Xpa knockout mice by the anti-inflammatory CXCL1 monoclonal antibody affected the pattern of somatic mutations, wherein the transcriptional strand asymmetry was attenuated and the activated signal transduction was shifted from the RAS/RAF/MAPK to the PIK3CA pathway.

High incidence of PI3K pathway gene mutations in South Indian cervical cancers.

INTRODUCTION: Cervical cancer is the second most common cancer in India. The phosphatidylinositol-3 kinase (PI3K) signaling is one of the most commonly activated pathways in cancer and comprises key molecules commonly targeted in cancer therapy. This study analyzed six PI3K pathway gene mutations. METHODS: We carried out targeted next-generation sequencing of six PI3K pathway genes (PIK3CA, PIK3R1, PTEN, AKT1, TSC2, and mTOR) in a total of 93 South Indian cervical cancer samples and confirmed them by sanger sequencing. RESULTS: The PI3K pathway gene mutations were observed in 54.8% (51/93) of the tumors and PIK3CA was the most mutated (34.4%, 32/93), followed by TSC2 (18.3%, 17/93), and PIK3R1 (14%, 13/93). The PIK3CA hotspot mutations E542K and E545K observed in this study were likely to disrupt the p110α-p85α interaction that could result in the PI3K pathway activation. We also found a few novel mutations in PIK3R1, PTEN, AKT1, TSC2, and mTOR genes while some of the tumors harbored multiple mutations in the genes of the PI3K pathway. The majority of the tumors were positive for high-risk HPV16/18 (60.7%). CONCLUSIONS: The high incidence of the PI3K pathway gene mutations observed in this study could be exploited for the therapeutic management of cervical cancers.

A sheet pocket to prevent cross-contamination of formalin-fixed paraffin-embedded block for application in next generation sequencing.

Formalin-fixed paraffin-embedded (FFPE) blocks are used as biomaterials for next-generation sequencing of cancer panels. Cross-contamination is detected in approximately 5% of the DNA extracted from FFPE samples, which reduces the detection rate of genetic abnormalities. There are no effective methods available for processing FFPE blocks that prevent cells from mixing with other specimens. The present study evaluated 897 sheets that could potentially prevent cell transmission but allow for the movement of various solvents used in FFPE blocks. According to the International Organization for Standardization and Japanese Industrial Standards, six requirements were established for the screening of packing sheets: 1) filter opening ≤5 µm, 2) thickness ≤100 µm, 3) chemical resistance, 4) permeability ≥1.0 × 10-3 cm/s, 5) water retention rate <200%, and 6) cell transit test (≤2 cells/10 high-power fields). Polyamide, polyethylene terephthalate, and polypropylene/polyethylene composite sheets met all criteria. A pocket, which was designed to wrap the tissue uniformly, was made of these sheets and was found to effectively block the entry of all cell types during FFPE block processing. Using a sheet pocket, no single cell from the cell pellet could pass through the outer layer. The presence or absence of the sheet pocket did not affect hematoxylin and eosin staining. When processing FFPE blocks as a biomaterial for next-generation sequencing, the sheet pocket was effective in preventing cross-contamination. This technology will in part support the precise translation of histopathological data into genome sequencing data in general pathology laboratories.

Spatiotemporal dynamics of clonal selection and diversification in normal endometrial epithelium.

It has become evident that somatic mutations in cancer-associated genes accumulate in the normal endometrium, but spatiotemporal understanding of the evolution and expansion of mutant clones is limited. To elucidate the timing and mechanism of the clonal expansion of somatic mutations in cancer-associated genes in the normal endometrium, we sequence 1311 endometrial glands from 37 women. By collecting endometrial glands from different parts of the endometrium, we show that multiple glands with the same somatic mutations occupy substantial areas of the endometrium. We demonstrate that "rhizome structures", in which the basal glands run horizontally along the muscular layer and multiple vertical glands rise from the basal gland, originate from the same ancestral clone. Moreover, mutant clones detected in the vertical glands diversify by acquiring additional mutations. These results suggest that clonal expansions through the rhizome structures are involved in the mechanism by which mutant clones extend their territories. Furthermore, we show clonal expansions and copy neutral loss-of-heterozygosity events occur early in life, suggesting such events can be tolerated many years in the normal endometrium. Our results of the evolutionary dynamics of mutant clones in the human endometrium will lead to a better understanding of the mechanisms of endometrial regeneration during the menstrual cycle and the development of therapies for the prevention and treatment of endometrium-related diseases.

Genome-wide meta-analysis between renal overload type and renal underexcretion type of clinically defined gout in Japanese populations.

Despite progress in understanding of the genetic basis of gout, the precise factors affecting differences in gout susceptibility among different gout subtypes remain unclear. Using clinically diagnosed gout patients, we conducted a genome-wide meta-analysis of two distinct gout subtypes: the renal overload type and the renal underexcretion type. We provide genetic evidence at a genome-wide level of significance that supports a positive association between ABCG2 dysfunction and acquisition of the renal overload type.

Biased expression of mutant alleles in cancer-related genes in esophageal squamous cell carcinoma.

BACKGROUND: Recent progress of large-scale international studies has provided comprehensive catalogs of somatic mutations in cancers. Additionally, it has become evident that allelic imbalance in the abundance of somatic mutations between DNA and RNA were pervasive in various types of cancer. However, the allelic imbalance of the abundance of somatic mutations in esophageal squamous cell carcinoma (ESCC) has not been fully analyzed. METHODS: We performed exome sequencing for 25 Japanese patients with ESCC to detect a comprehensive catalog of somatic mutations in ESCC. Additionally, we performed mRNA sequencing to evaluate the allelic imbalance of the identified somatic mutations at the transcriptional level by comparing the mutant allele frequencies between RNA and DNA. RESULTS: The exome sequencing showed that TP53 and ZNF750 were significantly mutated genes. The expression levels of TP53 and ZNF750 were different depending on the mutation status. In almost all the tumors with missense mutations in TP53 and ZNF750, the mutant allele frequencies were higher in the RNA sequencing than those in the exome sequencing, indicating that the mutant alleles were preferentially expressed. By examining the allelic imbalances for all the identified missense mutations, we demonstrated that genes showing preferential expressions of the mutant alleles were involved in the pathways including cell cycle, cell death, and chromatin modification. CONCLUSIONS: The results of this study suggest that the allelic imbalance of the abundance of somatic mutations plays important roles in the initiation and progression of ESCC by modulating cancer-related biological pathways.

APOBEC mediated mutagenesis drives genomic heterogeneity in endometriosis.

Endometriosis is a benign gynecologic condition, acting as a precursor of certain histological subtypes of ovarian cancers. The epithelial cells of endometriotic tissues and normal uterine endometrium accumulated somatic mutations in cancer-associated genes such as phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) and Kirsten rat sarcoma (KRAS) proto-oncogene. To determine the genomic characteristic of endometriotic epithelial cells and normal uterine endometrium and to identify the predominant mutational process acting on them, we studied the somatic mutation profiles obtained from whole exome sequencing of 14 endometriotic epithelium and 11 normal uterine endometrium tissues and classified them into mutational signatures. We observed that single base substitutions 2/13 (SBS), attributed to Apolipoprotein B mRNA Editing Enzyme Catalytic Subunit (APOBEC) induced mutagenesis, were significant in endometriotic tissues, but not in the normal uterine endometrium. Additionally, the larger number and wider allele frequency distribution of APOBEC signature mutations, compared to cancer-associated driver mutations in endometriotic epithelium suggested APOBEC mutagenesis as an important source of mutational burden and heterogeneity in endometriosis. Further, the relative risk of enriched APOBEC signature mutations was higher in endometriosis patients who were carriers of APOBEC3A/3B germline deletion, a common polymorphism in East Asians which involves the complete loss of APOBEC3B coding region. Our results illustrate the significance of APOBEC induced mutagenesis in driving the genomic heterogeneity of endometriosis.

Allelic and haplotypic HLA diversity in indigenous Malaysian populations explored using Next Generation Sequencing.

The heterogenous population of Malaysia includes more than 50 indigenous groups, and characterizing their HLA diversity would not only provide insights to their ancestry, but also on the effects of natural selection on their genome. We utilized hybridization-based sequence capture and short-read sequencing on the HLA region of 172 individuals representing seven indigenous groups in Malaysia (Jehai, Kintaq, Temiar, Mah Meri, Seletar, Temuan, Bidayuh). Allele and haplotype frequencies of HLA-A, -B, -C, -DRB1, -DQA1, -DQB1, -DPA1, and -DPB1 revealed several ancestry-informative markers. Using SNP-based heterozygosity and pairwise Fst, we observed signals of natural selection, particularly in HLA-A, -C and -DPB1 genes. Consequently, we showed the impact of natural selection on phylogenetic inference using HLA and non-HLA SNPs. We demonstrate the utility of Next Generation Sequencing for generating unambiguous, high-throughput, high-resolution HLA data that adds to our knowledge of HLA diversity and natural selection in indigenous minority groups.

Decoding the diversity of killer immunoglobulin-like receptors by deep sequencing and a high-resolution imputation method.

The killer cell immunoglobulin-like receptor (KIR) recognizes human leukocyte antigen (HLA) class I molecules and modulates the function of natural killer cells. Despite its role in immunity, the complex genomic structure has limited a deep understanding of the KIR genomic landscape. Here we conduct deep sequencing of 16 KIR genes in 1,173 individuals. We devise a bioinformatics pipeline incorporating copy number estimation and insertion or deletion (indel) calling for high-resolution KIR genotyping. We define 118 alleles in 13 genes and demonstrate a linkage disequilibrium structure within and across KIR centromeric and telomeric regions. We construct a KIR imputation reference panel (nreference = 689, imputation accuracy = 99.7%), apply it to biobank genotype (ntotal = 169,907), and perform phenome-wide association studies of 85 traits. We observe a dearth of genome-wide significant associations, even in immune traits implicated previously to be associated with KIR (the smallest p = 1.5 × 10-4). Our pipeline presents a broadly applicable framework to evaluate innate immunity in large-scale datasets.

Comprehensive discovery of CRISPR-targeted terminally redundant sequences in the human gut metagenome: Viruses, plasmids, and more.

Viruses are the most numerous biological entity, existing in all environments and infecting all cellular organisms. Compared with cellular life, the evolution and origin of viruses are poorly understood; viruses are enormously diverse, and most lack sequence similarity to cellular genes. To uncover viral sequences without relying on either reference viral sequences from databases or marker genes that characterize specific viral taxa, we developed an analysis pipeline for virus inference based on clustered regularly interspaced short palindromic repeats (CRISPR). CRISPR is a prokaryotic nucleic acid restriction system that stores the memory of previous exposure. Our protocol can infer CRISPR-targeted sequences, including viruses, plasmids, and previously uncharacterized elements, and predict their hosts using unassembled short-read metagenomic sequencing data. By analyzing human gut metagenomic data, we extracted 11,391 terminally redundant CRISPR-targeted sequences, which are likely complete circular genomes. The sequences included 2,154 tailed-phage genomes, together with 257 complete crAssphage genomes, 11 genomes larger than 200 kilobases, 766 genomes of Microviridae species, 56 genomes of Inoviridae species, and 95 previously uncharacterized circular small genomes that have no reliably predicted protein-coding gene. We predicted the host(s) of approximately 70% of the discovered genomes at the taxonomic level of phylum by linking protospacers to taxonomically assigned CRISPR direct repeats. These results demonstrate that our protocol is efficient for de novo inference of CRISPR-targeted sequences and their host prediction.

Proposing a molecular classification associated with hypercoagulation in ovarian clear cell carcinoma.

BACKGROUND: Although ovarian clear cell carcinoma (CCC) is associated with high incidence of thromboembolism, the clinicopathological and biological significance of hypercoagulable status in CCC remains unclear. MATERIALS AND METHODS: We retrospectively analyzed pretreatment D-dimer levels, thromboembolic status, and clinical outcome of 125 CCCs in the discovery set and 143 CCCs in two other independent validation sets. Next, we performed RNA sequencing of 93 CCCs and compared coagulation-related gene profiles with 2492 pan-cancer data. We investigated differences in molecular characteristics of CCC subclasses based on coagulation status. RESULTS: In the discovery dataset, D-dimer elevation above the normal range was significantly associated with shorter progression-free and overall survival, irrespective to thromboembolic status. Multivariate analysis identified D-dimer elevation and clinical stage as an independent prognostic factors. We confirmed the prognostic significance of D-dimer elevation in the validation sets. Tissue factor and IL6, which are considered key elements of cancer-induced hypercoagulation, were highly expressed in CCC than in other cancers regardless of D-dimer level. Higher activity of various oncogenic pathways was observed in CCC with compared to without D-dimer elevation. Moreover, hierarchical cluster analysis divided 57 CCCs with D-dimer elevation into immunologically hot and cold tumor subtypes. Hot tumors were characterized by enrichment of T-cell inflamed phenotype, inflammation, the epithelial-mesenchymal transition, and high serum levels of CRP, and cold tumors by enrichment of cell cycle and MYC pathways. CONCLUSIONS: CCC represents hypercoagulable disease and elevate D-dimer is a prognostic factor for decreased survival in CCC. D-dimer high CCC has distinct molecular characteristics into the inflammatory-driven pathway (hot tumor) and the immune-suppressive pathway (cold tumor). Treatment implication of our proposed molecular classification merits further investigation.

Three-dimensional understanding of the morphological complexity of the human uterine endometrium.

The fundamental morphology of the endometrial glands is not sufficiently understood by 2D observation because these glands have complicated winding and branching patterns. To construct a large picture of the endometrial gland structure, we performed tissue-clearing-based 3D imaging of human uterine endometrial tissue. Our 3D immunohistochemistry and layer analyses revealed that the endometrial glands form a plexus network in the stratum basalis and expand horizontally along the muscular layer, similar to the rhizome of grass. We then extended our method to assess the 3D morphology of tissue affected by adenomyosis, a representative "endometrium-related disease," and observed its 3D morphological features, including the direct invasion of endometrial glands into the myometrium and an ant colony-like network of ectopic endometrial glands within the myometrium. Thus, further understanding of the morphology of the human endometrium based on 3D analysis will lead to the identification of the pathogenesis of endometrium-related diseases.