Genetic drivers and cellular selection of female mosaic X chromosome loss.

Mosaic loss of the X chromosome (mLOX) is the most common clonal somatic alteration in leukocytes of female individuals1,2, but little is known about its genetic determinants or phenotypic consequences. Here, to address this, we used data from 883,574 female participants across 8 biobanks; 12% of participants exhibited detectable mLOX in approximately 2% of leukocytes. Female participants with mLOX had an increased risk of myeloid and lymphoid leukaemias. Genetic analyses identified 56 common variants associated with mLOX, implicating genes with roles in chromosomal missegregation, cancer predisposition and autoimmune diseases. Exome-sequence analyses identified rare missense variants in FBXO10 that confer a twofold increased risk of mLOX. Only a small fraction of associations was shared with mosaic Y chromosome loss, suggesting that distinct biological processes drive formation and clonal expansion of sex chromosome missegregation. Allelic shift analyses identified X chromosome alleles that are preferentially retained in mLOX, demonstrating variation at many loci under cellular selection. A polygenic score including 44 allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Our results support a model in which germline variants predispose female individuals to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of clonal expansion.

Population analyses of mosaic X chromosome loss identify genetic drivers and widespread signatures of cellular selection.

Mosaic loss of the X chromosome (mLOX) is the most commonly occurring clonal somatic alteration detected in the leukocytes of women, yet little is known about its genetic determinants or phenotypic consequences. To address this, we estimated mLOX in >900,000 women across eight biobanks, identifying 10% of women with detectable X loss in approximately 2% of their leukocytes. Out of 1,253 diseases examined, women with mLOX had an elevated risk of myeloid and lymphoid leukemias and pneumonia. Genetic analyses identified 49 common variants influencing mLOX, implicating genes with established roles in chromosomal missegregation, cancer predisposition, and autoimmune diseases. Complementary exome-sequence analyses identified rare missense variants in FBXO10 which confer a two-fold increased risk of mLOX. A small fraction of these associations were shared with mosaic Y chromosome loss in men, suggesting different biological processes drive the formation and clonal expansion of sex chromosome missegregation events. Allelic shift analyses identified alleles on the X chromosome which are preferentially retained, demonstrating that variation at many loci across the X chromosome is under cellular selection. A novel polygenic score including 44 independent X chromosome allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Collectively our results support a model where germline variants predispose women to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of subsequent clonal expansion.

MicroRNA-27a-3p targets FoxO signalling to induce tumour-like phenotypes in bile duct cells.

BACKGROUND & AIMS: Cholangiocarcinoma (CCA) is a heterogeneous and lethal malignancy, the molecular origins of which remain poorly understood. MicroRNAs (miRs) target diverse signalling pathways, functioning as potent epigenetic regulators of transcriptional output. We aimed to characterise miRNome dysregulation in CCA, including its impact on transcriptome homeostasis and cell behaviour. METHODS: Small RNA sequencing was performed on 119 resected CCAs, 63 surrounding liver tissues, and 22 normal livers. High-throughput miR mimic screens were performed in three primary human cholangiocyte cultures. Integration of patient transcriptomes and miRseq together with miR screening data identified an oncogenic miR for characterization. MiR-mRNA interactions were investigated by a luciferase assay. MiR-CRISPR knockout cells were generated and phenotypically characterized in vitro (proliferation, migration, colony, mitochondrial function, glycolysis) and in vivo using subcutaneous xenografts. RESULTS: In total, 13% (140/1,049) of detected miRs were differentially expressed between CCA and surrounding liver tissues, including 135 that were upregulated in tumours. CCA tissues were characterised by higher miRNome heterogeneity and miR biogenesis pathway expression. Unsupervised hierarchical clustering of tumour miRNomes identified three subgroups, including distal CCA-enriched and IDH1 mutant-enriched subgroups. High-throughput screening of miR mimics uncovered 71 miRs that consistently increased proliferation of three primary cholangiocyte models and were upregulated in CCA tissues regardless of anatomical location, among which only miR-27a-3p had consistently increased expression and activity in several cohorts. FoxO signalling was predominantly downregulated by miR-27a-3p in CCA, partially through targeting of FOXO1. MiR-27a knockout increased FOXO1 levels in vitro and in vivo, impeding tumour behaviour and growth. CONCLUSIONS: The miRNomes of CCA tissues are highly remodelled, impacting transcriptome homeostasis in part through regulation of transcription factors like FOXO1. MiR-27a-3p arises as an oncogenic vulnerability in CCA. IMPACT AND IMPLICATIONS: Cholangiocarcinogenesis entails extensive cellular reprogramming driven by genetic and non-genetic alterations, but the functional roles of these non-genetic events remain poorly understood. By unveiling global miRNA upregulation in patient tumours and their functional ability to increase proliferation of cholangiocytes, these small non-coding RNAs are implicated as critical non-genetic alterations promoting biliary tumour initiation. These findings identify possible mechanisms for transcriptome rewiring during transformation, with potential implications for patient stratification.

Whole blood microRNAs capture systemic reprogramming and have diagnostic potential in patients with biliary tract cancer.

BACKGROUND & AIMS: Late diagnosis is a critical factor undermining clinical management of patients with biliary tract cancer (BTC). While biliary tumours display extensive inter-patient heterogeneity, the host immune response may be comparatively homogenous, providing diagnostic opportunities. Herein, we investigated whether cancer-associated systemic reprogramming could be detected non-invasively to improve diagnosis of BTC. METHODS: In this prospective Danish study, whole blood (WB) microRNA (miRNA) profiling was performed in samples from 218 patients with BTC, 99 healthy participants, and 69 patients with differential diagnoses split into discovery (small RNA-sequencing) and validation (RT-qPCR) cohorts. miRNA expression and activity were further investigated in 119 and 660 BTC tissues, respectively. RESULTS: Four WB miRNAs (let-7a-3p, miR-92b-5p, miR-145-3p, miR-582-3p) were identified and validated as diagnostic of BTC on univariable analysis. Two diagnostic miRNA indexes were subsequently identified that were elevated in patients with BTC and in patients with differential diagnoses, compared to healthy participants. The combination of these miRNA indexes with serum CA 19-9 significantly improved the diagnostic performance of CA 19-9 alone, consistently achieving superior AUC values irrespective of clinical setting (minimum AUC >0.84) or tumour location (minimum AUC >0.87). The diagnostic information captured by miRNA indexes was not recapitulated by routine clinical measurements. Index miRNA expression in BTC tissues was associated with distinct pathobiological and immune features. CONCLUSIONS: WB miRNA profiles are altered in patients with BTC. Quantification of miRNA indexes in combination with serum CA 19-9 has the potential to improve early diagnosis of BTC, pending further validation. LAY SUMMARY: Surgery is currently the only curative intervention for patients with biliary tract cancer (BTC). However, resection is not possible for most patients who are diagnosed with late-stage disease. With the aim of identifying new early diagnostic opportunities, we analysed circulating microRNAs (small non-coding RNAs whose role in cancer is being increasingly recognised) in whole blood samples. We identified a microRNA signature that could distinguish patients with BTC from healthy participants. These miRNAs significantly improved the diagnostic potential of the routinely measured biomarker, CA 19-9, and were implicated in distinct immune processes in tumour tissues.

ZBED6 modulates the transcription of myogenic genes in mouse myoblast cells.

ZBED6 is a recently discovered transcription factor, unique to placental mammals, that has evolved from a domesticated DNA transposon. It acts as a repressor at the IGF2 locus. Here we show that ZBED6 acts as a transcriptional modulator in mouse myoblast cells, where more than 700 genes were differentially expressed after Zbed6-silencing. The most significantly enriched GO term was muscle protein and contractile fiber, which was consistent with increased myotube formation. Twenty small nucleolar RNAs all showed increased expression after Zbed6-silencing. The co-localization of histone marks and ZBED6 binding sites and the effect of Zbed6-silencing on distribution of histone marks was evaluated by ChIP-seq analysis. There was a strong association between ZBED6 binding sites and the H3K4me3, H3K4me2 and H3K27ac modifications, which are usually found at active promoters, but no association with the repressive mark H3K27me3. Zbed6-silencing led to increased enrichment of active marks at myogenic genes, in agreement with the RNA-seq findings. We propose that ZBED6 preferentially binds to active promoters and modulates transcriptional activity without recruiting repressive histone modifications.

ZBED evolution: repeated utilization of DNA transposons as regulators of diverse host functions.

ZBED genes originate from domesticated hAT DNA transposons and encode regulatory proteins of diverse function in vertebrates. Here we reveal the evolutionary relationship between ZBED genes and demonstrate that they are derived from at least two independent domestication events in jawed vertebrate ancestors. We show that ZBEDs form two monophyletic clades, one of which has expanded through several independent duplications in host lineages. Subsequent diversification of ZBED genes has facilitated regulation of multiple diverse fundamental functions. In contrast to known examples of transposable element exaptation, our results demonstrate a novel unprecedented capacity for the repeated utilization of a family of transposable element-derived protein domains sequestered as regulators during the evolution of diverse host gene functions in vertebrates. Specifically, ZBEDs have contributed to vertebrate regulatory innovation through the donation of modular DNA and protein interacting domains. We identify that C7ORF29, ZBED2, 3, 4, and ZBEDX form a monophyletic group together with ZBED6, that is distinct from ZBED1 genes. Furthermore, we show that ZBED5 is related to Buster DNA transposons and is phylogenetically separate from other ZBEDs. Our results offer new insights into the evolution of regulatory pathways, and suggest that DNA transposons have contributed to regulatory complexity during genome evolution in vertebrates.

Phylogenetic and structural analysis of HCV nonstructural protein 4A from Pakistani patients.

Hepatitis C virus nonstructural protein, NS4A, is a small protein comprising of about 54 amino acids. Despite its small size, it plays key role in many viral and cellular functions. The most important of which is its role as the co-factor of viral serine protease and helicase (NS3). Our study examines the phylogenetic and structural analysis of this coding region after isolation from Pakistani HCV patient samples. Phylogenetic analysis of the gene revealed that Pakistani 3a HCV strains do not show significant divergence from those reported from the rest of the world. The findings of this study also depict that NS4A sequence is conserved within genotypes, whereas it shows variations among different genotypes. While predicting the tertiary structure of the protein two important mutations (H28Y & E32G) were observed when comparing the Pakistani sequences with that of a reference HCV (genotype 3a) strain NZL (D17763). These mutations were observed in the central domain of NS4A which is responsible for interaction with NS3. Taken together, these mutations within the NS4A coding region can play an important role in the binding capacity of NS4A with HCV serine protease NS3.