Single-Cell Transcriptome Atlas of Murine Endothelial Cells.

The heterogeneity of endothelial cells (ECs) across tissues remains incompletely inventoried. We constructed an atlas of >32,000 single-EC transcriptomes from 11 mouse tissues and identified 78 EC subclusters, including Aqp7+ intestinal capillaries and angiogenic ECs in healthy tissues. ECs from brain/testis, liver/spleen, small intestine/colon, and skeletal muscle/heart pairwise expressed partially overlapping marker genes. Arterial, venous, and lymphatic ECs shared more markers in more tissues than did heterogeneous capillary ECs. ECs from different vascular beds (arteries, capillaries, veins, lymphatics) exhibited transcriptome similarity across tissues, but the tissue (rather than the vessel) type contributed to the EC heterogeneity. Metabolic transcriptome analysis revealed a similar tissue-grouping phenomenon of ECs and heterogeneous metabolic gene signatures in ECs between tissues and between vascular beds within a single tissue in a tissue-type-dependent pattern. The EC atlas taxonomy enabled identification of EC subclusters in public scRNA-seq datasets and provides a powerful discovery tool and resource value.

Single-cell exome sequencing and monoclonal evolution of a JAK2-negative myeloproliferative neoplasm.

Tumor heterogeneity presents a challenge for inferring clonal evolution and driver gene identification. Here, we describe a method for analyzing the cancer genome at a single-cell nucleotide level. To perform our analyses, we first devised and validated a high-throughput whole-genome single-cell sequencing method using two lymphoblastoid cell line single cells. We then carried out whole-exome single-cell sequencing of 90 cells from a JAK2-negative myeloproliferative neoplasm patient. The sequencing data from 58 cells passed our quality control criteria, and these data indicated that this neoplasm represented a monoclonal evolution. We further identified essential thrombocythemia (ET)-related candidate mutations such as SESN2 and NTRK1, which may be involved in neoplasm progression. This pilot study allowed the initial characterization of the disease-related genetic architecture at the single-cell nucleotide level. Further, we established a single-cell sequencing method that opens the way for detailed analyses of a variety of tumor types, including those with high genetic complex between patients.

Single-cell exome sequencing reveals single-nucleotide mutation characteristics of a kidney tumor.

Clear cell renal cell carcinoma (ccRCC) is the most common kidney cancer and has very few mutations that are shared between different patients. To better understand the intratumoral genetics underlying mutations of ccRCC, we carried out single-cell exome sequencing on a ccRCC tumor and its adjacent kidney tissue. Our data indicate that this tumor was unlikely to have resulted from mutations in VHL and PBRM1. Quantitative population genetic analysis indicates that the tumor did not contain any significant clonal subpopulations and also showed that mutations that had different allele frequencies within the population also had different mutation spectrums. Analyses of these data allowed us to delineate a detailed intratumoral genetic landscape at a single-cell level. Our pilot study demonstrates that ccRCC may be more genetically complex than previously thought and provides information that can lead to new ways to investigate individual tumors, with the aim of developing more effective cellular targeted therapies.

Sequencing and de novo assembly of 150 genomes from Denmark as a population reference.

Hundreds of thousands of human genomes are now being sequenced to characterize genetic variation and use this information to augment association mapping studies of complex disorders and other phenotypic traits. Genetic variation is identified mainly by mapping short reads to the reference genome or by performing local assembly. However, these approaches are biased against discovery of structural variants and variation in the more complex parts of the genome. Hence, large-scale de novo assembly is needed. Here we show that it is possible to construct excellent de novo assemblies from high-coverage sequencing with mate-pair libraries extending up to 20 kilobases. We report de novo assemblies of 150 individuals (50 trios) from the GenomeDenmark project. The quality of these assemblies is similar to those obtained using the more expensive long-read technology. We use the assemblies to identify a rich set of structural variants including many novel insertions and demonstrate how this variant catalogue enables further deciphering of known association mapping signals. We leverage the assemblies to provide 100 completely resolved major histocompatibility complex haplotypes and to resolve major parts of the Y chromosome. Our study provides a regional reference genome that we expect will improve the power of future association mapping studies and hence pave the way for precision medicine initiatives, which now are being launched in many countries including Denmark.

Genome-wide annotation of protein-coding genes in pig.

BACKGROUND: There is a need for functional genome-wide annotation of the protein-coding genes to get a deeper understanding of mammalian biology. Here, a new annotation strategy is introduced based on dimensionality reduction and density-based clustering of whole-body co-expression patterns. This strategy has been used to explore the gene expression landscape in pig, and we present a whole-body map of all protein-coding genes in all major pig tissues and organs. RESULTS: An open-access pig expression map ( www.rnaatlas.org ) is presented based on the expression of 350 samples across 98 well-defined pig tissues divided into 44 tissue groups. A new UMAP-based classification scheme is introduced, in which all protein-coding genes are stratified into tissue expression clusters based on body-wide expression profiles. The distribution and tissue specificity of all 22,342 protein-coding pig genes are presented. CONCLUSIONS: Here, we present a new genome-wide annotation strategy based on dimensionality reduction and density-based clustering. A genome-wide resource of the transcriptome map across all major tissues and organs in pig is presented, and the data is available as an open-access resource ( www.rnaatlas.org ), including a comparison to the expression of human orthologs.

Haplotyping by CRISPR-mediated DNA circularization (CRISPR-hapC) broadens allele-specific gene editing.

Allele-specific protospacer adjacent motif (asPAM)-positioning SNPs and CRISPRs are valuable resources for gene therapy of dominant disorders. However, one technical hurdle is to identify the haplotype comprising the disease-causing allele and the distal asPAM SNPs. Here, we describe a novel CRISPR-based method (CRISPR-hapC) for haplotyping. Based on the generation (with a pair of CRISPRs) of extrachromosomal circular DNA in cells, the CRISPR-hapC can map haplotypes from a few hundred bases to over 200 Mb. To streamline and demonstrate the applicability of the CRISPR-hapC and asPAM CRISPR for allele-specific gene editing, we reanalyzed the 1000 human pan-genome and generated a high frequency asPAM SNP and CRISPR database (www.crispratlas.com/knockout) for four CRISPR systems (SaCas9, SpCas9, xCas9 and Cas12a). Using the huntingtin (HTT) CAG expansion and transthyretin (TTR) exon 2 mutation as examples, we showed that the asPAM CRISPRs can specifically discriminate active and dead PAMs for all 23 loci tested. Combination of the CRISPR-hapC and asPAM CRISPRs further demonstrated the capability for achieving highly accurate and haplotype-specific deletion of the HTT CAG expansion allele and TTR exon 2 mutation in human cells. Taken together, our study provides a new approach and an important resource for genome research and allele-specific (haplotype-specific) gene therapy.

Single-Cell RNA Sequencing Reveals Renal Endothelium Heterogeneity and Metabolic Adaptation to Water Deprivation.

BACKGROUND: Renal endothelial cells from glomerular, cortical, and medullary kidney compartments are exposed to different microenvironmental conditions and support specific kidney processes. However, the heterogeneous phenotypes of these cells remain incompletely inventoried. Osmotic homeostasis is vitally important for regulating cell volume and function, and in mammals, osmotic equilibrium is regulated through the countercurrent system in the renal medulla, where water exchange through endothelium occurs against an osmotic pressure gradient. Dehydration exposes medullary renal endothelial cells to extreme hyperosmolarity, and how these cells adapt to and survive in this hypertonic milieu is unknown. METHODS: We inventoried renal endothelial cell heterogeneity by single-cell RNA sequencing >40,000 mouse renal endothelial cells, and studied transcriptome changes during osmotic adaptation upon water deprivation. We validated our findings by immunostaining and functionally by targeting oxidative phosphorylation in a hyperosmolarity model in vitro and in dehydrated mice in vivo. RESULTS: We identified 24 renal endothelial cell phenotypes (of which eight were novel), highlighting extensive heterogeneity of these cells between and within the cortex, glomeruli, and medulla. In response to dehydration and hypertonicity, medullary renal endothelial cells upregulated the expression of genes involved in the hypoxia response, glycolysis, and-surprisingly-oxidative phosphorylation. Endothelial cells increased oxygen consumption when exposed to hyperosmolarity, whereas blocking oxidative phosphorylation compromised endothelial cell viability during hyperosmotic stress and impaired urine concentration during dehydration. CONCLUSIONS: This study provides a high-resolution atlas of the renal endothelium and highlights extensive renal endothelial cell phenotypic heterogeneity, as well as a previously unrecognized role of oxidative phosphorylation in the metabolic adaptation of medullary renal endothelial cells to water deprivation.

LION: a simple and rapid method to achieve CRISPR gene editing.

The RNA-guided CRISPR-Cas9 technology has paved the way for rapid and cost-effective gene editing. However, there is still a great need for effective methods for rapid generation and validation of CRISPR/Cas9 gRNAs. Previously, we have demonstrated that highly efficient generation of multiplexed CRISPR guide RNA (gRNA) expression array can be achieved with Golden Gate Assembly (GGA). Here, we present an optimized and rapid method for generation and validation in less than 1 day of CRISPR gene targeting vectors. The method (LION) is based on ligation of double-stranded gRNA oligos into CRISPR vectors with GGA followed by nucleic acid purification. Using a dual-fluorescent reporter vector (C-Check), T7E1 assay, TIDE assay and a traffic light reporter assay, we proved that the LION-based generation of CRISPR vectors are functionally active, and equivalent to CRISPR plasmids generated by traditional methods. We also tested the activity of LION CRISPR vectors in different human cell types. The LION method presented here advances the rapid functional validation and application of CRISPR system for gene editing and simplified the CRISPR gene-editing procedures.

CRISPR-C: circularization of genes and chromosome by CRISPR in human cells.

Extrachromosomal circular DNA (eccDNA) and ring chromosomes are genetic alterations found in humans with genetic disorders. However, there is a lack of genetic engineering tools to recapitulate and study the biogenesis of eccDNAs. Here, we created a dual-fluorescence biosensor cassette, which upon the delivery of pairs of CRISPR/Cas9 guide RNAs, CRISPR-C, allows us to study the biogenesis of a specific fluorophore expressing eccDNA in human cells. We show that CRISPR-C can generate functional eccDNA, using the novel eccDNA biosensor system. We further reveal that CRISPR-C also can generate eccDNAs from intergenic and genic loci in human embryonic kidney 293T cells and human mammary fibroblasts. EccDNAs mainly forms by end-joining mediated DNA-repair and we show that CRISPR-C is able to generate endogenous eccDNAs in sizes from a few hundred base pairs and ranging up to 207 kb. Even a 47.4 megabase-sized ring chromosome 18 can be created by CRISPR-C. Our study creates a new territory for CRISPR gene editing and highlights CRISPR-C as a useful tool for studying the cellular impact, persistence and function of eccDNAs.

WEGO 2.0: a web tool for analyzing and plotting GO annotations, 2018 update.

WEGO (Web Gene Ontology Annotation Plot), created in 2006, is a simple but useful tool for visualizing, comparing and plotting GO (Gene Ontology) annotation results. Owing largely to the rapid development of high-throughput sequencing and the increasing acceptance of GO, WEGO has benefitted from outstanding performance regarding the number of users and citations in recent years, which motivated us to update to version 2.0. WEGO uses the GO annotation results as input. Based on GO's standardized DAG (Directed Acyclic Graph) structured vocabulary system, the number of genes corresponding to each GO ID is calculated and shown in a graphical format. WEGO 2.0 updates have targeted four aspects, aiming to provide a more efficient and up-to-date approach for comparative genomic analyses. First, the number of input files, previously limited to three, is now unlimited, allowing WEGO to analyze multiple datasets. Also added in this version are the reference datasets of nine model species that can be adopted as baselines in genomic comparative analyses. Furthermore, in the analyzing processes each Chi-square test is carried out for multiple datasets instead of every two samples. At last, WEGO 2.0 provides an additional output graph along with the traditional WEGO histogram, displaying the sorted P-values of GO terms and indicating their significant differences. At the same time, WEGO 2.0 features an entirely new user interface. WEGO is available for free at http://wego.genomics.org.cn.

Analyses of pig genomes provide insight into porcine demography and evolution.

For 10,000 years pigs and humans have shared a close and complex relationship. From domestication to modern breeding practices, humans have shaped the genomes of domestic pigs. Here we present the assembly and analysis of the genome sequence of a female domestic Duroc pig (Sus scrofa) and a comparison with the genomes of wild and domestic pigs from Europe and Asia. Wild pigs emerged in South East Asia and subsequently spread across Eurasia. Our results reveal a deep phylogenetic split between European and Asian wild boars ∼1 million years ago, and a selective sweep analysis indicates selection on genes involved in RNA processing and regulation. Genes associated with immune response and olfaction exhibit fast evolution. Pigs have the largest repertoire of functional olfactory receptor genes, reflecting the importance of smell in this scavenging animal. The pig genome sequence provides an important resource for further improvements of this important livestock species, and our identification of many putative disease-causing variants extends the potential of the pig as a biomedical model.

Systematic in vitro and in vivo characterization of Leukemia-inhibiting factor- and Fibroblast growth factor-derived porcine induced pluripotent stem cells.

Derivation and stable maintenance of porcine induced pluripotent stem cells (piPSCs) is challenging. We herein systematically analyzed two piPSC lines, derived by lentiviral transduction and cultured under either leukemia inhibitory factor (LIF) or fibroblast growth factor (FGF) conditions, to shed more light on the underlying biological mechanisms of porcine pluripotency. LIF-derived piPSCs were more successful than their FGF-derived counterparts in the generation of in vitro chimeras and in teratoma formation. When LIF piPSCs chimeras were transferred into surrogate sows and allowed to develop, only their prescence within the embryonic membranes could be detected. Whole-transcriptome analysis of the piPSCs and porcine neonatal fibroblasts showed that they clustered together, but apart from the two pluripotent cell populations of early porcine embryos, indicating incomplete reprogramming. Indeed, bioinformatic analysis of the pluripotency-related gene network of the LIF- versus FGF-derived piPSCs revealed that ZFP42 (REX1) expression was absent in both piPSC-like cells, whereas it was expressed in the porcine inner cell mass at Day 7/8. A second striking difference was the expression of ATOH1 in piPSC-like cells, which was absent in the inner cell mass. Moreover, our gene expression analyses plus correlation analyses of known pluripotency genes identified unique relationships between pluripotency genes in the inner cell mass, which are to some extent, in the piPSC-like cells. This deficiency in downstream gene activation and divergent gene expression may be underlie the inability to derive germ line-transmitting piPSCs, and provides unique insight into which genes are necessary to achieve fully reprogrammed piPSCs. 84: 229-245, 2017. © 2016 Wiley Periodicals, Inc.

Isolation and whole genome sequencing of fetal cells from maternal blood towards the ultimate non-invasive prenatal testing.

OBJECTIVE: The purpose of this study were to develop a methodology of isolating fetal cells from maternal blood and use deep sequence demonstrating the promise for complete and accurate genetic screening compared to other non-invasive prenatal testing. METHODS: Here in this study, we developed a double negative selection (DNS) procedure to unbiasedly enrich fetal cells. After validated by short tandem repeat (STR), the isolated circulating fetal cells (CFCs) were subjected to deep whole genome sequencing analysis. RESULTS: Our DNS protocol significantly increasing the purity of the mimic fetal cells from 1 in 1 million nucleated cells in whole blood to 1:8 to 1:30 (12.5%-3.33%) after 2 steps of enrichment. Isolated single fetal cell obtained a coverage rate (86.8%) and allelic dropout rate (24.90%) comparative to the reported results of human cell line. Several disease-associated variants were identified in the whole genome sequencing data of isolated CFCs and further confirmed in the sequencing data of unamplified gDNA. CONCLUSION: In conclusion, the robustness of DNS and STR to collect CFCs from peripheral maternal blood for the first time coupled with deep sequencing technique demonstrates the possibility of comprehensive non-invasive prenatal testing of genetic disorders using isolated CFCs.

Golden Gate Assembly of CRISPR gRNA expression array for simultaneously targeting multiple genes.

The engineered CRISPR/Cas9 technology has developed as the most efficient and broadly used genome editing tool. However, simultaneously targeting multiple genes (or genomic loci) in the same individual cells using CRISPR/Cas9 remain one technical challenge. In this article, we have developed a Golden Gate Assembly method for the generation of CRISPR gRNA expression arrays, thus enabling simultaneous gene targeting. Using this method, the generation of CRISPR gRNA expression array can be accomplished in 2 weeks, and contains up to 30 gRNA expression cassettes. We demonstrated in the study that simultaneously targeting 10 genomic loci or simultaneously inhibition of multiple endogenous genes could be achieved using the multiplexed gRNA expression array vector in human cells. The complete set of plasmids is available through the non-profit plasmid repository Addgene.

Enhanced genome editing in mammalian cells with a modified dual-fluorescent surrogate system.

Programmable DNA nucleases such as TALENs and CRISPR/Cas9 are emerging as powerful tools for genome editing. Dual-fluorescent surrogate systems have been demonstrated by several studies to recapitulate DNA nuclease activity and enrich for genetically edited cells. In this study, we created a single-strand annealing-directed, dual-fluorescent surrogate reporter system, referred to as C-Check. We opted for the Golden Gate Cloning strategy to simplify C-Check construction. To demonstrate the utility of the C-Check system, we used the C-Check in combination with TALENs or CRISPR/Cas9 in different scenarios of gene editing experiments. First, we disrupted the endogenous pIAPP gene (3.0 % efficiency) by C-Check-validated TALENs in primary porcine fibroblasts (PPFs). Next, we achieved gene-editing efficiencies of 9.0-20.3 and 4.9 % when performing single- and double-gene targeting (MAPT and SORL1), respectively, in PPFs using C-Check-validated CRISPR/Cas9 vectors. Third, fluorescent tagging of endogenous genes (MYH6 and COL2A1, up to 10.0 % frequency) was achieved in human fibroblasts with C-Check-validated CRISPR/Cas9 vectors. We further demonstrated that the C-Check system could be applied to enrich for IGF1R null HEK293T cells and CBX5 null MCF-7 cells with frequencies of nearly 100.0 and 86.9 %, respectively. Most importantly, we further showed that the C-Check system is compatible with multiplexing and for studying CRISPR/Cas9 sgRNA specificity. The C-Check system may serve as an alternative dual-fluorescent surrogate tool for measuring DNA nuclease activity and enrichment of gene-edited cells, and may thereby aid in streamlining programmable DNA nuclease-mediated genome editing and biological research.

The stepwise evolution of the exome during acquisition of docetaxel resistance in breast cancer cells.

BACKGROUND: Resistance to taxane-based therapy in breast cancer patients is a major clinical problem that may be addressed through insight of the genomic alterations leading to taxane resistance in breast cancer cells. In the current study we used whole exome sequencing to discover somatic genomic alterations, evolving across evolutionary stages during the acquisition of docetaxel resistance in breast cancer cell lines. RESULTS: Two human breast cancer in vitro models (MCF-7 and MDA-MB-231) of the step-wise acquisition of docetaxel resistance were developed by exposing cells to 18 gradually increasing concentrations of docetaxel. Whole exome sequencing performed at five successive stages during this process was used to identify single point mutational events, insertions/deletions and copy number alterations associated with the acquisition of docetaxel resistance. Acquired coding variation undergoing positive selection and harboring characteristics likely to be functional were further prioritized using network-based approaches. A number of genomic changes were found to be undergoing evolutionary selection, some of which were likely to be functional. Of the five stages of progression toward resistance, most resistance relevant genomic variation appeared to arise midway towards fully resistant cells corresponding to passage 31 (5 nM docetaxel) for MDA-MB-231 and passage 16 (1.2 nM docetaxel) for MCF-7, and where the cells also exhibited a period of reduced growth rate or arrest, respectively. MCF-7 cell acquired several copy number gains on chromosome 7, including ABC transporter genes, including ABCB1 and ABCB4, as well as DMTF1, CLDN12, CROT, and SRI. For MDA-MB-231 numerous copy number losses on chromosome X involving more than 30 genes was observed. Of these genes, CASK, POLA1, PRDX4, MED14 and PIGA were highly prioritized by the applied network-based gene ranking approach. At higher docetaxel concentration MCF-7 subclones exhibited a copy number loss in E2F4, and the gene encoding this important transcription factor was down-regulated in MCF-7 resistant cells. CONCLUSIONS: Our study of the evolution of acquired docetaxel resistance identified several genomic changes that might explain development of docetaxel resistance. Interestingly, the most relevant resistance-associated changes appeared to originate midway through the evolution towards fully resistant cell lines. Our data suggest that no single genomic event sufficiently predicts resistance to docetaxel, but require genomic alterations affecting multiple pathways that in concert establish the final resistance stage.

SOAPindel: efficient identification of indels from short paired reads.

We present a new approach to indel calling that explicitly exploits that indel differences between a reference and a sequenced sample make the mapping of reads less efficient. We assign all unmapped reads with a mapped partner to their expected genomic positions and then perform extensive de novo assembly on the regions with many unmapped reads to resolve homozygous, heterozygous, and complex indels by exhaustive traversal of the de Bruijn graph. The method is implemented in the software SOAPindel and provides a list of candidate indels with quality scores. We compare SOAPindel to Dindel, Pindel, and GATK on simulated data and find similar or better performance for short indels (<10 bp) and higher sensitivity and specificity for long indels. A validation experiment suggests that SOAPindel has a false-positive rate of ∼10% for long indels (>5 bp), while still providing many more candidate indels than other approaches.

BS-SNPer: SNP calling in bisulfite-seq data.

UNLABELLED: Sodium bisulfite conversion followed by sequencing (BS-Seq, such as whole genome bisulfite sequencing or reduced representation bisulfite sequencing) has become popular for studying human epigenetic profiles. Identifying single nucleotide polymorphisms (SNPs) is important for quantification of methylation levels and for study of allele-specific epigenetic events such as imprinting. However, SNP calling in such data is complex and time consuming. Here, we present an ultrafast and memory-efficient package named BS-SNPer for the exploration of SNP sites from BS-Seq data. Compared with Bis-SNP, a popular BS-Seq specific SNP caller, BS-SNPer is over 100 times faster and uses less memory. BS-SNPer also offers higher sensitivity and specificity compared with existing methods. AVAILABILITY AND IMPLEMENTATION: BS-SNPer is written in C++ and Perl, and is freely available at https://github.com/hellbelly/BS-Snper.

EWS and FUS bind a subset of transcribed genes encoding proteins enriched in RNA regulatory functions.

BACKGROUND: FUS (TLS) and EWS (EWSR1) belong to the FET-protein family of RNA and DNA binding proteins. FUS and EWS are structurally and functionally related and participate in transcriptional regulation and RNA processing. FUS and EWS are identified in translocation generated cancer fusion proteins and involved in the human neurological diseases amyotrophic lateral sclerosis and fronto-temporal lobar degeneration. RESULTS: To determine the gene regulatory functions of FUS and EWS at the level of chromatin, we have performed chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq). Our results show that FUS and EWS bind to a subset of actively transcribed genes, that binding often is downstream the poly(A)-signal, and that binding overlaps with RNA polymerase II. Functional examinations of selected target genes identified that FUS and EWS can regulate gene expression at different levels. Gene Ontology analyses showed that FUS and EWS target genes preferentially encode proteins involved in regulatory processes at the RNA level. CONCLUSIONS: The presented results yield new insights into gene interactions of EWS and FUS and have identified a set of FUS and EWS target genes involved in pathways at the RNA regulatory level with potential to mediate normal and disease-associated functions of the FUS and EWS proteins.

The potential role of Alu Y in the development of resistance to SN38 (Irinotecan) or oxaliplatin in colorectal cancer.

BACKGROUND: Irinotecan (SN38) and oxaliplatin are chemotherapeutic agents used in the treatment of colorectal cancer. However, the frequent development of resistance to these drugs represents a considerable challenge in the clinic. Alus as retrotransposons comprise 11% of the human genome. Genomic toxicity induced by carcinogens or drugs can reactivate Alus by altering DNA methylation. Whether or not reactivation of Alus occurs in SN38 and oxaliplatin resistance remains unknown. RESULTS: We applied reduced representation bisulfite sequencing (RRBS) to investigate the DNA methylome in SN38 or oxaliplatin resistant colorectal cancer cell line models. Moreover, we extended the RRBS analysis to tumor tissue from 14 patients with colorectal cancer who either did or did not benefit from capecitabine + oxaliplatin treatment. For the clinical samples, we applied a concept of 'DNA methylation entropy' to estimate the diversity of DNA methylation states of the identified resistance phenotype-associated methylation loci observed in the cell line models. We identified different loci being characteristic for the different resistant cell lines. Interestingly, 53% of the identified loci were Alu sequences- especially the Alu Y subfamily. Furthermore, we identified an enrichment of Alu Y sequences that likely results from increased integration of new copies of Alu Y sequence in the drug-resistant cell lines. In the clinical samples, SOX1 and other SOX gene family members were shown to display variable DNA methylation states in their gene regions. The Alu Y sequences showed remarkable variation in DNA methylation states across the clinical samples. CONCLUSION: Our findings imply a crucial role of Alu Y in colorectal cancer drug resistance. Our study underscores the complexity of colorectal cancer aggravated by mobility of Alu elements and stresses the importance of personalized strategies, using a systematic and dynamic view, for effective cancer therapy.