Sox Genes Research Articles

Anaerobic thiosulfate oxidation by the Roseobacter group is prevalent in marine biofilms

Thiosulfate oxidation by microbes has a major impact on global sulfur cycling. Here, we provide evidence that bacteria within various Roseobacter lineages are important for thiosulfate oxidation in marine biofilms. We isolate and sequence the genomes of 54 biofilm-associated Roseobacter strains, finding conserved sox gene clusters for thiosulfate oxidation and plasmids, pointing to a niche-specific lifestyle. Analysis of global ocean metagenomic data suggests that Roseobacter strains are abundant in biofilms and mats on various substrates, including stones, artificial surfaces, plant roots, and hydrothermal vent chimneys. Metatranscriptomic analysis indicates that the majority of active sox genes in biofilms belong to Roseobacter strains. Furthermore, we show that Roseobacter strains can grow and oxidize thiosulfate to sulfate under both aerobic and anaerobic conditions. Transcriptomic and membrane proteomic analyses of biofilms formed by a representative strain indicate that thiosulfate induces sox gene expression and alterations in cell membrane protein composition, and promotes biofilm formation and anaerobic respiration. We propose that bacteria of the Roseobacter group are major thiosulfate-oxidizers in marine biofilms, where anaerobic thiosulfate metabolism is preferred.

Coordinated microRNA/mRNA Expression Profiles Reveal Unique Skin Color Regulatory Mechanisms in Chinese Giant Salamander (Andrias davidianus).

The Chinese giant salamander (Andrias davidianus) has been increasingly popular in the aquaculture market in China in recent years. In the breeding process of Andrias davidianus, we found that some albino individuals were extremely rare and could not be inherited stably, which severely limits their commercialization in the aquaculture market. In this study, we performed transcriptome and small RNA (sRNA) sequencing analyses in the skin samples of wild-type (WT) and albino (AL) Andrias davidianus. In total, among 5517 differentially expressed genes (DEGs), 2911 DEGs were down-regulated in AL, including almost all the key genes involved in melanin formation. A total of 25 miRNAs were differentially expressed in AL compared to WT, of which 17 were up-regulated. Through the integrated analysis, no intersection was found between the target genes of the differentially expressed miRNAs and the key genes for melanin formation. Gene Ontology (GO) and KEGG pathway analyses on DEGs showed that these genes involved multiple processes relevant to melanin synthesis and the key signal pathway MAPK. Interestingly, the transcription factors SOX10 and PAX3 and the Wnt signaling pathway that play a key role in other species were not included, while the other two transcription factors in the SOX family, SOX21 and SOX7, were included. After analyzing the key genes for melanin formation, it was interesting to note an alternative splicing form of the MITF in WT and a critical mutation of the SLC24A5 gene in AL, which might be the main reason for the skin color change of Andrias davidianus. The results contributed to understanding the molecular mechanism of skin pigmentation in Andrias davidianus and accelerating the acquisition process of individuals with specific body colors by genetic means.

SOX8 is essential for male sexual differentiation in the Chinese soft-shelled turtle Pelodiscus sinensis‡.

SOX8, which belongs to SOXE transcription factor subfamily together with SOX9, participates in sex differentiation and testicular development by enhancing the function of SOX9 in mammals. However, the functional role of SOX8 in sex differentiation has not yet been identified in any non-mammalian vertebrates. Here, we found in the Chinese soft-shelled turtle Pelodiscus sinensis that SOX8 exhibited male-specific higher expression from stage 14 to 18, the critical period of sex determination, prior to the onset of gonadal differentiation. In addition, SOX8 was rapidly down-regulated during male-to-female sex reversal induced by estradiol. Moreover, knockdown of SOX8 led to complete feminization of ZZ P. sinensis, determined by gonadal morphology and distribution of germ cells, as well as the down-regulation of testicular marker DMRT1 and the up-regulation of ovarian regulator FOXL2. In conclusion, this study provides evidence that SOX8 is a key regulator of early male differentiation in P. sinensis, highlighting the significance of the SOX family in reptile sex determination.

Identification and characterization of a novel SoxB2 gene in Litopenaeus vannamei and its potential roles in response to ammonia-N stress and post-exposure recovery

Increased concentrations of ammonia-N can cause seriously physiological harm to crustaceans, especially farmed shrimp (Litopenaeus vannamei). Ammonia-N is generally removed by water change, and will gradually accumulate for unconsumed feed. However, the underlying adaptation mechanism of shrimp under such environmental conditions is still not fully understood. A novel Sox gene expressed differentially in previous comparative transcriptome studies was characterized in this study, designated as LvSoxB2–1. The LvSoxB2–1 encodes a protein of 420 amino acids, sharing the highest identity with the Chinese mitten crab (Eriocheir sinensis) Sox14b. The expression level of LvSoxB2–1 was higher in fertilized egg and cleavage stage of early embryonic development and declined in the blastula stage, which showed no significant difference compared with the subsequent stage. During larval development, LvSoxB2–1 increased significantly in nauplius stage III and zoea stage I and then restored to the initial level. LvSoxB2–1 was highly expressed in the gills, followed by eyestalks and brain ganglion. Results of subcellular localization have shown that LvSoxB2–1 is a nuclear-localized protein. In the ammonia-N stress and post-exposure recovery experiments, LvSoxB2–1 mRNA significantly increased in the gills, hepatopancreas, and hemocytes. The total hemocyte counts (THC) of shrimp was decreased through LvSoxB2–1 RNA interference (RNAi). Moreover, the inhibition of LvSoxB2–1 significantly exacerbated the hepatopancreas pathology and reduced the survival rates of shrimp exposed to ammonia-N. Dual-luciferase reporter assay confirmed the activation of the Wnt/β-catenin pathway by LvSoxB2–1. Under ammonia-N stress and recovery after exposure, LvSoxB2–1 silencing could decrease the expression of Wnt/β-catenin signaling molecules (β-catenin and Pangolin) and Rh protein in L. vannamei, with the higher malonyldialdehyde (MDA) and protein carbonylation (PCO) levels and hemolymph ammonia-N and urea contents. These results suggested that LvSoxB2–1 played positive roles in oxidative stress responses and ammonia-N excretion through the Wnt/β-catenin pathway.

Haemophilia B, severe childhood obesity and other extra-haematological features associated with similar 4Mb-deletions on Xq27: Clinical findings, molecular insights and literature update.

Haemophilia B (HB) is associated with pathogenic variants in F9. Hemizygous deletions encompassing the entire F9 and proximate genes may express extra-haematological clinical phenotypes. To analyse the genotype/phenotype correlations in two unrelated boys with severe early childhood obesity (SCO), global developmental delay (GDD) and similar bleeding phenotype associated with comparable Xq27 deletions spanning the entire F9 and proximate genes, and characterise the pathogenic events estimating the most likely mutational mechanism involved. Entire F9-deletions were detected in three hemizygous unrelated probands with HB: two cases, C#1/C#2, presented SCO and GDD and a control patient (Co), who only had severe bleeding symptoms. Dense SNP-array and case-specific STS walking scan allowed characterisation of the deletion breakpoints. Extensive use of bioinformatics, statistics and clinical databases allowed the investigation of genotype-phenotype associations. Patients C#1/C#2 and Co resulted in a complete F9 and additional gene deletions of variable extensions on Xq26.3-Xq27.2 (C#1/C#2/Co: 4.3Mb/3.9Mb/160Kb). C#1/C#2 common deleted gene SOX3 is directly associated with SCO, GDD and pituitary hypothyroidism (PH) whilst C#2 extra-deleted gene MAGEC2 indirectly relates to anal atresia (AA). Breakpoint analysis revealed the involvement of the mechanisms of Alu/Alu recombination for the first time in HB and non-homologous or alternative end-joining. Our results represent the first report of unrelated patients with HB, SCO and GDD. This study and the literature update expand the spectrum of clinical findings and molecular insights observed in patients with HB caused by complete F9 and nearby SOX3 and MAGEC2 gene deletions, which may configure a contiguous gene syndrome.

The NFkB activation domain is 14-amino-acid-long variant of the 9aaTAD.

The nine-amino-acid transactivation domains (9aaTAD) was identified in numerous transcription factors including Gal4, p53, E2A, MLL, c-Myc, N-Myc, and also in SP, KLF, and SOX families. Most of the 9aaTAD domains interact with the KIX domain of transcription mediators MED15 and CBP to activate transcription. The NFkB activation domain occupied the same position on the KIX domain as the 9aaTADs of MLL, E2A, and p53. Binding of the KIX domain is established by the two-point interaction involving 9aaTAD positions p3-4 and p6-7. The NFkB primary binding region (positions p3-4) is almost identical with MLL and E2A, but secondary NFkB binding region differs by the position and engages the distal NFkB region p10-11. Thus, the NFkB activation domain is five amino acids longer than the other 9aaTADs. The NFkB activation domain includes an additional region, which we called the Omichinski Insert extending activation domain length to 14 amino acids. By deletion, we demonstrated that Omichinski Insert is an entirely non-essential part of NFkB activation domain. In summary, we recognized the NFkB activation domain as prolonged 9aaTAD conserved in evolution from humans to amphibians.

A new insight into potential roles of Spfoxl-2 in the testicular development of Scylla paramamosain by RNAi and transcriptome analysis

In our previous study, we found that the Spfoxl-2 transcript was highly expressed in gonads and explored its potential target genes in the ovary of Scylla paramamosain. In the current study, we primally analyzed its potential target genes in the testis through RNAi and RNA-Seq technology and compared with that in the ovary. The results showed that a total of 7892 unigenes were differentially expressed after Spfoxl-2 silencing in the testis, including plenty of conserved genes involved in testicular development, such as Dmrt family genes, Sox family genes, Caspase family genes, Cdk family genes, Kinesin family genes, Fox family genes and other genes. Further analysis revealed that these differentially expressed genes (DEGs) were enriched in crucial pathways involved in spermatogenesis, such as DNA replication, Cell cycle, Spliceosome, Homologous recombination, Meiosis and Apoptosis. The comparison results of potential target genes in the ovary and testis reveal 135 common potential target genes, including some genes involved in the immune response. According to our knowledge, the present work was the first to disclose the functions of foxl-2 in the testis of crustacean species using transcriptome analysis. It not only identifies key genes and pathways involved in the testicular development of S. paramamosain, but also reveals a new molecular-level understanding of the function of foxl-2 in testicular development.

Impact of genome and epigenome on intratumor heterogeneity in colorectal cancer

Impact of genome and epigenome on intratumor heterogeneity in colorectal cancer

Transcription factor SOX15 regulates stem cell pluripotency and promotes neural fate during differentiation by activating the neurogenic gene Hes5

SOX2 and SOX15 are Sox family transcription factors enriched in embryonic stem cells (ESCs). The role of SOX2 in activating gene expression programs essential for stem cell self-renewal and acquisition of pluripotency during somatic cell reprogramming is well-documented. However, the contribution of SOX15 to these processes is unclear and often presumed redundant with SOX2 largely because overexpression of SOX15 can partially restore self-renewal in SOX2-deficient ESCs. Here, we show that SOX15 contributes to stem cell maintenance by cooperating with ESC-enriched transcriptional coactivators to ensure optimal expression of pluripotency-associated genes. We demonstrate that SOX15 depletion compromises reprogramming of fibroblasts to pluripotency which cannot be compensated by SOX2. Ectopic expression of SOX15 promotes the reversion of a post-implantation, epiblast stem cell state back to a pre-implantation, ESC-like identity even though SOX2 is expressed in both cell states. We also uncover a role of SOX15 in lineage specification, by showing that loss of SOX15 leads to defects in commitment of ESCs to neural fates. SOX15 promotes neural differentiation by binding to and activating a previously uncharacterized distal enhancer of a key neurogenic regulator, Hes5. Together, these findings identify a multifaceted role of SOX15 in induction and maintenance of pluripotency and neural differentiation.

Transcriptional Ontogeny Platform of Sex Determination/Sex Differentiation Genes in Almaco Jack Larvae

AbstractThe fish sex determination/sex differentiation (SD) system is still intriguing in many species, such as Seriola fish. Marine amberjacks have a conserved ZW/ZZ SD system. We aimed to identify how early SD genes were expressed in larval Almaco Jack Seriola rivoliana by RNA sequencing and to determine whether there was an early developmental differential expression up to the onset of exogenous feeding. Transcriptome analysis identified sox5 and sox7 (female) and sox3 and sox8 (male) genes; all of the known female SD genes, including hsd17b1 (master sex gene in Seriola), foxl2, and cyp19a, plus genes involved in female sexual development and maintenance, estrogen signaling, and osteoclast proliferation; and all of the known male SD genes, including dmrt1, nr5a1, gsdf, amhy, and amhr2. The expression of most SD genes occurred on each day during larval development. However, most of them showed low abundance levels and null differential expression except for the downregulation of the sox and foxl2 genes, which might have served an unknown function or were maternally derived rather than autonomously downregulated. We suggest that the concomitant expression of female and male SD genes in early larval stages is a biological transcriptional ontogeny sex platform that ensures posterior sex‐specific expression and sex differentiation. A single missense single‐nucleotide polymorphism (SNP) near the catalytic center of hsd17b1 has been suggested to differentiate phenotypic sex in Seriola by an endocrine regulation on estradiol synthesis. In that case, the ontogeny sex platform would be directed by the SNP but still critical to differentiate gonads.

Identification and functional analysis of Dmrt1 gene and the SoxE gene in the sexual development of sea cucumber, Apostichopus japonicus.

Members of the Doublesex and Mab-3-related transcription factor (Dmrt) gene family handle various vital functions in several biological processes, including sex determination/differentiation and gonad development. Dmrt1 and Sox9 (SoxE in invertebrates) exhibit a very conserved interaction function during testis formation in vertebrates. However, the dynamic expression pattern and functional roles of the Dmrt gene family and SoxE have not yet been identified in any echinoderm species. Herein, five members of the Dmrt gene family (Dmrt1, 2, 3a, 3b and 5) and the ancestor SoxE gene were identified from the genome of Apostichopus japonicus. Expression studies of Dmrt family genes and SoxE in different tissues of adult males and females revealed different expression patterns of each gene. Transcription of Dmrt2, Dmrt3a and Dmrt3b was higher expressed in the tube feet and coelomocytes instead of in gonadal tissues. The expression of Dmrt1 was found to be sustained throughout spermatogenesis. Knocking-down of Dmrt1 by means of RNA interference (RNAi) led to the downregulation of SoxE and upregulation of the ovarian regulator foxl2 in the testes. This indicates that Dmrt1 may be a positive regulator of SoxE and may play a role in the development of the testes in the sea cucumber. The expression level of SoxE was higher in the ovaries than in the testes, and knocking down of SoxE by RNAi reduced SoxE and Dmrt1 expression but conversely increased the expression of foxl2 in the testes. In summary, this study indicates that Dmrt1 and SoxE are indispensable for testicular differentiation, and SoxE might play a functional role during ovary differentiation in the sea cucumber.

Evolutionary Landscape of SOX Genes to Inform Genotype-to-Phenotype Relationships.

The SOX transcription factor family is pivotal in controlling aspects of development. To identify genotype-phenotype relationships of SOX proteins, we performed a non-biased study of SOX using 1890 open-reading frame and 6667 amino acid sequences in combination with structural dynamics to interpret 3999 gnomAD, 485 ClinVar, 1174 Geno2MP, and 4313 COSMIC human variants. We identified, within the HMG (High Mobility Group)- box, twenty-seven amino acids with changes in multiple SOX proteins annotated to clinical pathologies. These sites were screened through Geno2MP medical phenotypes, revealing novel SOX15 R104G associated with musculature abnormality and SOX8 R159G with intellectual disability. Within gnomAD, SOX18 E137K (rs201931544), found within the HMG box of ~0.8% of Latinx individuals, is associated with seizures and neurological complications, potentially through blood-brain barrier alterations. A total of 56 highly conserved variants were found at sites outside the HMG-box, including several within the SOX2 HMG-box-flanking region with neurological associations, several in the SOX9 dimerization region associated with Campomelic Dysplasia, SOX14 K88R (rs199932938) flanking the HMG box associated with cardiovascular complications within European populations, and SOX7 A379V (rs143587868) within an SOXF conserved far C-terminal domain heterozygous in 0.716% of African individuals with associated eye phenotypes. This SOX data compilation builds a robust genotype-to-phenotype association for a gene family through more robust ortholog data integration.

Additional evidence for the role of chromosomal imbalances and SOX8, ZNRF3 and HHAT gene variants in early human testis development

BackgroundForty-six ,XY Differences/Disorders of Sex Development (DSD) are characterized by a broad phenotypic spectrum ranging from typical female to male with undervirilized external genitalia, or more rarely testicular regression with a typical male phenotype. Despite progress in the genetic diagnosis of DSD, most 46,XY DSD cases remain idiopathic.MethodsTo determine the genetic causes of 46,XY DSD, we studied 165 patients of Tunisian ancestry, who presented a wide range of DSD phenotypes. Karyotyping, candidate gene sequencing, and whole-exome sequencing (WES) were performed.ResultsCytogenetic abnormalities, including a high frequency of sex chromosomal anomalies (85.4%), explained the phenotype in 30.9% (51/165) of the cohort. Sanger sequencing of candidate genes identified a novel pathogenic variant in the SRY gene in a patient with 46,XY gonadal dysgenesis. An exome screen of a sub-group of 44 patients with 46,XY DSD revealed pathogenic or likely pathogenic variants in 38.6% (17/44) of patients.ConclusionRare or novel pathogenic variants were identified in the AR, SRD5A2, ZNRF3, SOX8, SOX9 and HHAT genes. Overall our data indicate a genetic diagnosis rate of 41.2% (68/165) in the group of 46,XY DSD.

The Role of SOX Transcription Factors in Ageing and Age-Related Diseases.

The quest for eternal youth and immortality is as old as humankind. Ageing is an inevitable physiological process accompanied by many functional declines that are driving factors for age-related diseases. Stem cell exhaustion is one of the major hallmarks of ageing. The SOX transcription factors play well-known roles in self-renewal and differentiation of both embryonic and adult stem cells. As a consequence of ageing, the repertoire of adult stem cells present in various organs steadily declines, and their dysfunction/death could lead to reduced regenerative potential and development of age-related diseases. Thus, restoring the function of aged stem cells, inducing their regenerative potential, and slowing down the ageing process are critical for improving the health span and, consequently, the lifespan of humans. Reprograming factors, including SOX family members, emerge as crucial players in rejuvenation. This review focuses on the roles of SOX transcription factors in stem cell exhaustion and age-related diseases, including neurodegenerative diseases, visual deterioration, chronic obstructive pulmonary disease, osteoporosis, and age-related cancers. A better understanding of the molecular mechanisms of ageing and the roles of SOX transcription factors in this process could open new avenues for developing novel strategies that will delay ageing and prevent age-related diseases.

Single-cell transcriptome profiling reveals distinct expression patterns among Sox genes in the mouse incisor dental pulp.

SOX transcription factors play key roles in cell differentiation and cell fate determination during development. Using single-cell RNA-sequencing data, we examined the expression profiles of Sox genes in the mouse incisor dental pulp. Our analysis showed that Sox4, Sox5, Sox9, Sox11, and Sox12 are mainly expressed in mesenchymal stem/stromal cells (MSCs) representing osteogenic cells at different stages of differentiation. We found that in several MSCs, Sox genes co-expressed with regulatory genes such as Sp7, Satb2, Msx1, Snai2, Dlx1, Twist2, and Tfap2a. In addition, Sox family genes colocalized with Runx2 and Lef1, which are highly enriched in MSCs undergoing osteoblast differentiation. A protein interaction network analysis uncovered that CREBBP, CEBPB, TLE1, TWIST1, and members of the HDAC and SMAD families are interacting partners of RUNX2 and LEF1 during skeletal development. Collectively, the distinct expression patterns of the SOX transcription factors suggest that they play essential regulatory roles in directing lineage-specific gene expression during differentiation of MSCs.

The Initial Myelination in the Central Nervous System.

Myelination contributes not only to the rapid nerve conduction but also to axonal insulation and protection. In the central nervous system (CNS), the initial myelination features a multistep process where oligodendrocyte precursor cells undergo proliferation and migration before differentiating into mature oligodendrocytes. Mature oligodendrocytes then extend processes and wrap around axons to form the multilayered myelin sheath. These steps are tightly regulated by various cellular and molecular mechanisms, such as transcription factors (Olig family, Sox family), growth factors (PDGF, BDNF, FGF-2, IGF), chemokines/cytokines (TGF-β, IL-1β, TNFα, IL-6, IFN-γ), hormones (T3), axonal signals (PSA-NCAM, L1-CAM, LINGO-1, neural activity), and intracellular signaling pathways (Wnt/β-catenin, PI3 K/AKT/mTOR, ERK/MAPK). However, the fundamental mechanisms for initial myelination are yet to be fully elucidated. Identifying pivotal mechanisms for myelination onset, development, and repair will become the focus of future studies. This review focuses on the current understanding of how CNS myelination is initiated and also the regulatory mechanisms underlying the process.

SOX9 and HMGB3 co-operatively transactivate NANOG and promote prostate cancer progression.

The homeodomain-containing transcription factor NANOG is overexpressed in prostate adenocarcinoma (PCa) and predicts poor prognosis. The SOX family transcription factor SOX9, as well as the transcription co-activator HMGB3 of the HMGB family, are also overexpressed and may play pivotal roles in PCa. However, it is unknown whether SOX9 and HMGB3 interact with each other, or if they regulate NANOG gene transcription. We identified potential SOX9 responsive elements in NANOG promoter, and investigated if SOX9 regulated NANOG transcription in co-operation with HMGB3 by experimental analysis of potential SOX9 binding sites in NANOG promoter, reporter gene transcription assays with or without interference or artificial overexpression of SOX9 and/or HMGB3, and protein-binding assays of SOX9-HMGB3 interaction. Clinicopathologic and prognostic significance of SOX9-HMGB3 overexpression in PCa was analyzed. SOX9 activated NANOG gene transcription by preferentially binding to a highly conserved consensus cis-regulatory element (-573 to -568) in NANOG promoter, and promoted the expression of NANOG downstream oncogenic genes. Importantly, HMGB3 functioned as a partner of SOX9 to co-operatively enhance transactivation of NANOG by interacting with SOX9, predominantly via the HMG Box A domain of HMGB3. Overexpression of SOX9 and/or HMGB3 enhanced PCa cell survival and cell migration and were significantly associated with PCa progression. Notably, Cox proportional regression analysis showed that co-overexpression of both SOX9 and HMGB3 was an independent unfavorable prognosticator for both CRPC-free survival (relative risk [RR] = 3.779，95% confidence interval [CI]: 1.159-12.322, p = 0.028) and overall survival (RR = 3.615，95% CI: 1.101-11.876, p = 0.034). These findings showed a novel SOX9/HMGB3/NANOG regulatory mechanism, deregulation of which played important roles in PCa progression.

Cloning and Expression of Sox2 and Sox9 in Embryonic and Gonadal Development of Lutraria sieboldii

The Sox family plays essential roles as transcription factors in vertebrates; however, little is known about the Sox family in Lutraria sieboldii. L. sieboldii are pleasant to eat with a short growth cycle and have become one of the best bottom-seeded enrichment species in Guang Xi. In this study, Sox2 (named LsSox2) and Sox9 (named LsSox9) from L. sieboldii were cloned, and their expression patterns were analyzed. The length of the LsSox2 gene coding sequence was 1011 bp, encoding 336 amino acids, and LsSox9 was 1449 bp, encoding 482 amino acids. LsSox2 had its highest expression levels in the ovary, which were 356 times those in testis, whereas LsSox9 presented higher expression in testis, which was 6 times more highly expressed than in the ovary. LsSox2 exhibited the highest expression during the morula stage, which was 20 times that of the D-shaped larvae or zygote. LsSox9 exhibited two expression peaks, one at the four-cell stage and the other at the trochophore stage, while the lowest expression was in the zygote. LsSox9 was 73 times more highly expressed in the four-cell stage than in the zygote stage. During gonadal development, LsSox2 presented the highest expression in the mature ovary, which was 756 times more highly expressed than in mature testis. LsSox9 presented higher expression in testis at the emission stage which was 6 times more highly expressed than in the ovary. These results indicate that LsSox2 and LsSox9 may play important roles in embryonic and gonadal development.

Systematic Investigation of the Multifaceted Role of SOX11 in Cancer.

SRY-box transcription factor 11 (SOX11), as a member of the SOX family, is a transcription factor involved in the regulation of specific biological processes and has recently been found to be a prognostic marker for certain cancers. However, the roles of SOX11 in cancer remain controversial. Our study aimed to explore the various aspects of SOX11 in pan-cancer. The expression of SOX11 was investigated by the Genotype Tissue-Expression (GTEX) dataset and the Cancer Genome Atlas (TCGA) database. The protein level of SOX11 in tumor tissues and tumor-adjacent tissues was verified by human pan-cancer tissue microarray. Additionally, we used TCGA pan-cancer data to analyze the correlations among SOX11 expression and survival outcomes, clinical features, stemness, microsatellite instability (MSI), tumor mutation burden (TMB), mismatch repair (MMR) related genes and the tumor immune microenvironment. Furthermore, the cBioPortal database was applied to investigate the gene alterations of SOX11. The main biological processes of SOX11 in cancers were analyzed by Gene Set Enrichment Analysis (GSEA). As a result, aberrant expression of SOX11 has been implicated in 27 kinds of cancer types. Aberrant SOX11 expression was closely associated with survival outcomes, stage, tumor recurrence, MSI, TMB and MMR-related genes. In addition, the most frequent alteration of the SOX11 genome was mutation. Our study also showed the correlations of SOX11 with the level of immune infiltration in various cancers. In summary, our findings underline the multifaceted role and prognostic value of SOX11 in pan-cancer.

Human papillomavirus integration perspective in small cell cervical carcinoma

Small cell cervical carcinoma (SCCC) is a rare but aggressive malignancy. Here, we report human papillomavirus features and genomic landscape in SCCC via high-throughput HPV captured sequencing, whole-genome sequencing, whole-transcriptome sequencing, and OncoScan microarrays. HPV18 infections and integrations are commonly detected. Besides MYC family genes (37.9%), we identify SOX (8.4%), NR4A (6.3%), ANKRD (7.4%), and CEA (3.2%) family genes as HPV-integrated hotspots. We construct the genomic local haplotype around HPV-integrated sites, and find tandem duplications and amplified HPV long control regions (LCR). We propose three prominent HPV integration patterns: duplicating oncogenes (MYCN, MYC, and NR4A2), forming fusions (FGFR3–TACC3 and ANKRD12–NDUFV2), and activating genes (MYC) via the cis-regulations of viral LCRs. Moreover, focal CNA amplification peaks harbor canonical cancer genes including the HPV-integrated hotspots within MYC family, SOX2, and others. Our findings may provide potential molecular criteria for the accurate diagnosis and efficacious therapies for this lethal disease.