Initial Characterization of Canine T-cell Receptor Repertoires Using RNA-seq Data from Different Diseases and Tissues.

The dog serves as a key translational model in cancer immunotherapy. Understanding the T cell receptor (TCR) repertoire is needed for various cancer immunotherapies. Compared to humans where >300 million TCRs have been identified, <100 canine TCRs are reported. To address this deficiency, we assembled >200,000 complete TCR complementarity-determining region 3 (CDR3) sequences from RNA-seq data published for ~2,000 canine samples of blood, lymph node, and other tissues, of which 613 are tumors. We collected 1,324 human RNA-seq samples to compare the similarities and differences in T-cell repertoires between humans and dogs. Notably, our analysis revealed distinct variable gene usage patterns between blood samples and solid tissues in both canine and human samples for TRA and TRB loci. Moreover, our investigation led to the discovery of novel V gene and allele candidates in the canine genome. Our findings also revealed that the canine CDR3 resembled human CDR3 in terms of length and motifs. Additionally, our study unveiled shared traits in cancer TCRs between dogs and humans, including longer lengths and higher hydrophobicity of private CDR3s. Our results indicated the diversity of canine to be more comparable to that of humans than mice. Our study provides an initial landscape of the canine TCR repertoire, highlighting both its similarities and differences with the human counterpart, thus laying the groundwork for future research in comparative immunology and vaccine development.

Canine major histocompatibility complex class I (MHC-I) diversity landscape.

The genes of the Major Histocompatibility Complex class I (MHC-I) are among the most diverse in the mammalian genome, playing a crucial role in immunology. Understanding the diversity landscape of MHC-I is therefore of paramount importance. The dog is a key translational model in various biomedical fields. However, our understanding of the canine MHC-I diversity landscape lags significantly behind that of humans. To address this deficiency, we used our newly developed software, KPR de novo assembler and genotyper, to genotype 1,325 samples from 1,025 dogs with paired-end RNA-seq data from 43 BioProjects, after extensive quality control. Among 926 dogs that pass the QC, 591 dogs (64%) have at least one allele genotyped, and a total of 97 known alleles and 52 putative new alleles were identified. Further analysis reveals that DLA-I gene expression levels vary among the tissues, with lowest for testis and brain tissues and highest for blood, corpus luteum, and spleen. We identified dominant alleles in each of the 17 canine breeds, as well as among the entire canine population. Furthermore, our analysis also identifies breed-specific alleles and mutually co-occurred/exclusive alleles. Our study indicates that canine DLA-88 is as diversified as human HLA-A/B/C genes within the entire population, but less diversified within a breed than with HLA-A/B/C within an ethnic group. Lastly, we examined the hypervariable regions (HVR) within or across human/canine MHC-I alleles and found that 80% of the HVRs overlap between the two species. We further noted that 80% of the HVRs are within 4A contact with the peptides, and that the dog-human difference overlaps with only 20% HVRs. Our research offers valuable insights for immunological studies involving dogs.

Leading the pack: Best practices in comparative canine cancer genomics to inform human oncology.

Pet dogs develop spontaneous cancers at a rate estimated to be five times higher than that of humans, providing a unique opportunity to study disease biology and evaluate novel therapeutic strategies in a model system that possesses an intact immune system and mirrors key aspects of human cancer biology. Despite decades of interest, effective utilization of pet dog cancers has been hindered by a limited repertoire of necessary cellular and molecular reagents for both in vitro and in vivo studies, as well as a dearth of information regarding the genomic landscape of these cancers. Recently, many of these critical gaps have been addressed through the generation of a highly annotated canine reference genome, the creation of several tools necessary for multi-omic analysis of canine tumours, and the development of a centralized repository for key genomic and associated clinical information from canine cancer patients, the Integrated Canine Data Commons. Together, these advances have catalysed multidisciplinary efforts designed to integrate the study of pet dog cancers more effectively into the translational continuum, with the ultimate goal of improving human outcomes. The current review summarizes this recent progress and provides a guide to resources and tools available for comparative study of pet dog cancers.

Human basal-like breast cancer is represented by one of the two mammary tumor subtypes in dogs.

BACKGROUND: About 20% of breast cancers in humans are basal-like, a subtype that is often triple-negative and difficult to treat. An effective translational model for basal-like breast cancer is currently lacking and urgently needed. To determine whether spontaneous mammary tumors in pet dogs could meet this need, we subtyped canine mammary tumors and evaluated the dog-human molecular homology at the subtype level. METHODS: We subtyped 236 canine mammary tumors from 3 studies by applying various subtyping strategies on their RNA-seq data. We then performed PAM50 classification with canine tumors alone, as well as with canine tumors combined with human breast tumors. We identified feature genes for human BLBC and luminal A subtypes via machine learning and used these genes to repeat canine-alone and cross-species tumor classifications. We investigated differential gene expression, signature gene set enrichment, expression association, mutational landscape, and other features for dog-human subtype comparison. RESULTS: Our independent genome-wide subtyping consistently identified two molecularly distinct subtypes among the canine tumors. One subtype is mostly basal-like and clusters with human BLBC in cross-species PAM50 and feature gene classifications, while the other subtype does not cluster with any human breast cancer subtype. Furthermore, the canine basal-like subtype recaptures key molecular features (e.g., cell cycle gene upregulation, TP53 mutation) and gene expression patterns that characterize human BLBC. It is enriched in histological subtypes that match human breast cancer, unlike the other canine subtype. However, about 33% of canine basal-like tumors are estrogen receptor negative (ER-) and progesterone receptor positive (PR+), which is rare in human breast cancer. Further analysis reveals that these ER-PR+ canine tumors harbor additional basal-like features, including upregulation of genes of interferon-γ response and of the Wnt-pluripotency pathway. Interestingly, we observed an association of PGR expression with gene silencing in all canine tumors and with the expression of T cell exhaustion markers (e.g., PDCD1) in ER-PR+ canine tumors. CONCLUSIONS: We identify a canine mammary tumor subtype that molecularly resembles human BLBC overall and thus could serve as a vital translational model of this devastating breast cancer subtype. Our study also sheds light on the dog-human difference in the mammary tumor histology and the hormonal cycle.

Shared hotspot mutations in oncogenes position dogs as an unparalleled comparative model for precision therapeutics.

Naturally occurring canine cancers have remarkable similarities to their human counterparts. To better understand these similarities, we investigated 671 client-owned dogs from 96 breeds with 23 common tumor types, including those whose mutation profile are unknown (anal sac carcinoma and neuroendocrine carcinoma) or understudied (thyroid carcinoma, soft tissue sarcoma and hepatocellular carcinoma). We discovered mutations in 50 well-established oncogenes and tumor suppressors, and compared them to those reported in human cancers. As in human cancer, TP53 is the most commonly mutated gene, detected in 22.5% of canine tumors overall. Canine tumors share mutational hotspots with human tumors in oncogenes including PIK3CA, KRAS, NRAS, BRAF, KIT and EGFR. Hotspot mutations with significant association to tumor type include NRAS G61R and PIK3CA H1047R in hemangiosarcoma, ERBB2 V659E in pulmonary carcinoma, and BRAF V588E (equivalent of V600E in humans) in urothelial carcinoma. Our findings better position canines as a translational model of human cancer to investigate a wide spectrum of targeted therapies.

Human basal-like breast cancer is represented by one of the two mammary tumor subtypes in dogs.

Background: About 20% of breast cancers in humans are basal-like, a subtype that is often triple negative and difficult to treat. An effective translational model for basal-like breast cancer (BLBC) is currently lacking and urgently needed. To determine if spontaneous mammary tumors in pet dogs could meet this need, we subtyped canine mammary tumors and evaluated the dog-human molecular homology at the subtype level. Methods: We subtyped 236 canine mammary tumors from 3 studies by applying various subtyping strategies on their RNA-seq data. We then performed PAM50 classification with canine tumors alone, as well as with canine tumors combined with human breast tumors. We investigated differential gene expression, signature gene set enrichment, expression association, mutational landscape, and other features for dog-human subtype comparison. Results: Our independent genome-wide subtyping consistently identified two molecularly distinct subtypes among the canine tumors. One subtype is mostly basal-like and clusters with human BLBC in cross-species PAM50 classification, while the other subtype does not cluster with any human breast cancer subtype. Furthermore, the canine basal-like subtype recaptures key molecular features (e.g., cell cycle gene upregulation, TP53 mutation) and gene expression patterns that characterize human BLBC. It is enriched histological subtypes that match human breast cancer, unlike the other canine subtype. However, about 33% of canine basal-like tumors are estrogen receptor negative (ER-) and progesterone receptor positive (PR+), which is rare in human breast cancer. Further analysis reveals that these ER-PR+ canine tumors harbor additional basal-like features, including upregulation of genes of interferon-γ response and of the Wnt-pluripotency pathway. Interestingly, we observed an association of PGR expression with gene silencing in all canine tumors, and with the expression of T cell exhaustion markers (e.g., PDCD1 ) in ER-PR+ canine tumors. Conclusions: We identify a canine mammary tumor subtype that molecularly resembles human BLBC overall, and thus could serve as a vital spontaneous animal model of this devastating breast cancer subtype. Our study also sheds light on the dog-human difference in the mammary tumor histology and the hormonal cycle.

A Kmer-based paired-end read de novo assembler and genotyper for canine MHC class I genotyping.

The major histocompatibility complex class I (MHC-I) genes are highly polymorphic. MHC-I genotyping is required for determining the peptide epitopes available to an individual's T-cell repertoire. Current genotyping software tools do not work for the dog, due to very limited known canine alleles. To address this, we developed a Kmer-based paired-end read (KPR) de novo assembler and genotyper, which assemble paired-end RNA-seq reads from MHC-I regions into contigs, and then genotype each contig and estimate its expression level. KPR tools outperform other popular software examined in typing new alleles. We used KPR tools to successfully genotype152 dogs from a published dataset. The study discovers 33 putative new alleles, finds dominant alleles in 4 dog breeds, and builds allele diversity and expression landscapes among the 152 dogs. Our software meets a significant need in biomedical research.

Canine tumor mutational burden is correlated with TP53 mutation across tumor types and breeds.

Spontaneous canine cancers are valuable but relatively understudied and underutilized models. To enhance their usage, we reanalyze whole exome and genome sequencing data published for 684 cases of >7 common tumor types and >35 breeds, with rigorous quality control and breed validation. Our results indicate that canine tumor alteration landscape is tumor type-dependent, but likely breed-independent. Each tumor type harbors major pathway alterations also found in its human counterpart (e.g., PI3K in mammary tumor and p53 in osteosarcoma). Mammary tumor and glioma have lower tumor mutational burden (TMB) (median < 0.5 mutations per Mb), whereas oral melanoma, osteosarcoma and hemangiosarcoma have higher TMB (median ≥ 1 mutations per Mb). Across tumor types and breeds, TMB is associated with mutation of TP53 but not PIK3CA, the most mutated genes. Golden Retrievers harbor a TMB-associated and osteosarcoma-enriched mutation signature. Here, we provide a snapshot of canine mutations across major tumor types and breeds.

A Qualitative Change in the Transcriptome Occurs after the First Cell Cycle and Coincides with Lumen Establishment during MDCKII Cystogenesis.

Madin-Darby canine kidney II (MDCKII) cells are widely used to study epithelial morphogenesis. To better understand this process, we performed time course RNA-seq analysis of MDCKII 3D cystogenesis, along with polarized 2D cells for comparison. Our study reveals a biphasic change in the transcriptome that occurs after the first cell cycle and coincides with lumen establishment. This change appears to be linked to translocation of ß-catenin, supported by analyses with AVL9- and DENND5A-knockdown clones, and regulation by HNF1B, supported by ATAC-seq study. These findings indicate a qualitative change model for transcriptome remodeling during epithelial morphogenesis, leading to cell proliferation decrease and cell polarity establishment. Furthermore, our study reveals that active mitochondria are retained and chromatin accessibility decreases in 3D cysts but not in 2D polarized cells. This indicates that 3D culture is a better model than 2D culture for studying epithelial morphogenesis.

Proteins with Evolutionarily Hypervariable Domains are Associated with Immune Response and Better Survival of Basal-like Breast Cancer Patients.

Maltase-glucoamylase (MGAM) and MGAM2 both belong to the glycoside hydrolase family 31. MGAM, a therapeutic target for type 2 diabetes, is α-1,4-glucosidase and expressed in the intestine to catalyze starch digestion. MGAM2, however, is largely uncharacterized. By investigating The Cancer Genome Atlas data, we found that among breast cancer subtypes, MGAM2 expression is nearly exclusive to basal-like breast cancers (BLBCs), whereas MGAM tends to express in luminal A breast cancers. Moreover, MGAM2 expression is associated with better patient survival and correlated with immune genes/signatures, unlike MGAM. Both genes have emerged in mammals, but diverged after the placental-marsupial split. In placentals, MGAM2 has likely lost its α-1,4-glucosidase activity due to mutations in key catalytic sites, and has acquired a large domain that is extracellular, threonine-rich and evolutionarily hypervariable (EHV). Guided by MGAM2 findings, our genome-wide search identified >1000 human proteins with EHV regions. These proteins are enriched in immune functions and molecules, including major histocompatibility complex proteins. Their genes are expressed higher in BLBCs and are associated with better patient survival, like MGAM2. Their EHV-coding sequences are rich in simple repeats and harbor more cancer passenger mutations. In conclusion, MGAM2 diverges from MGAM structurally and likely functionally in placentals. MGAM2 is among >1000 human proteins with EHV regions and associated with immune response. We propose that these EHV molecules may have significant implication in cancer immunotherapy and BLBC treatment.

Proliferative and Invasive Colorectal Tumors in Pet Dogs Provide Unique Insights into Human Colorectal Cancer.

Spontaneous tumors in pet dogs represent a valuable but undercharacterized cancer model. To better use this resource, we performed an initial global comparison between proliferative and invasive colorectal tumors from 20 canine cases, and evaluated their molecular homology to human colorectal cancer (CRC). First, proliferative canine tumors harbor overactivated WNT/ß-catenin pathways and recurrent CTNNB1 (ß-catenin) mutations S45F/P, D32Y and G34E. Invasive canine tumors harbor prominent fibroblast proliferation and overactivated stroma. Both groups have recurrent TP53 mutations. We observed three invasion patterns in canine tumors: collective, crypt-like and epithelial⁻mesenchymal transition (EMT). We detected enriched Helicobacter bilis and Alistipes finegoldii in proliferative and crypt-like tumors, but depleted mucosa-microbes in the EMT tumor. Second, guided by our canine findings, we classified 79% of 478 human colon cancers from The Cancer Genome Atlas into four subtypes: primarily proliferative, or with collective, crypt-like or EMT invasion features. Their molecular characteristics match those of canine tumors. We showed that consensus molecular subtype 4 (mesenchymal) of human CRC should be further divided into EMT and crypt-like subtypes, which differ in TGF-ß activation and mucosa-microbe content. Our canine tumors share the same pathogenic pathway as human CRCs. Dog-human integration identifies three CRC invasion patterns and improves CRC subtyping.

SEG - A Software Program for Finding Somatic Copy Number Alterations in Whole Genome Sequencing Data of Cancer.

As next-generation sequencing technology advances and the cost decreases, whole genome sequencing (WGS) has become the preferred platform for the identification of somatic copy number alteration (CNA) events in cancer genomes. To more effectively decipher these massive sequencing data, we developed a software program named SEG, shortened from the word "segment". SEG utilizes mapped read or fragment density for CNA discovery. To reduce CNA artifacts arisen from sequencing and mapping biases, SEG first normalizes the data by taking the log2-ratio of each tumor density against its matching normal density. SEG then uses dynamic programming to find change-points among a contiguous log2-ratio data series along a chromosome, dividing the chromosome into different segments. SEG finally identifies those segments having CNA. Our analyses with both simulated and real sequencing data indicate that SEG finds more small CNAs than other published software tools.

Collaborating genomic, transcriptomic and microbiomic alterations lead to canine extreme intestinal polyposis.

Extreme intestinal polyposis in pet dogs has not yet been reported in literature. We identified a dog patient who developed numerous intestinal polyps, with the severity resembling human classic familial adenomatous polyposis (FAP), except the jejunum-ileum junction being the most polyp-dense. We investigated this dog, in comparison with 22 other dogs with spontaneous intestinal tumors but no severe polyposis, and with numerous published human cancers. We found, not APC mutation, but three other alteration pathways as likely reasons of this canine extreme polyposis. First, somatic truncation mutation W411X of FBXW7, a component of an E3 ubiquitin ligase, over-activates MYC and cell cycle-promoting network, accelerating crypt cell proliferation. Second, genes of protein trafficking and localization are downregulated, likely associated with germline mutation G406D of STAMBPL1, a K63-deubiquitinase, and MYC network activation. This inhibits epithelial apical-basolateral polarity establishment, preventing crypt cell differentiation. Third, Bacteroides uniformis, a commensal gut anaerobe, thrives and expresses abundantly thioredoxin and nitroreductase. These bacterial products could reduce oxidative stress linked to host germline mutation R51X of CYB5RL, a cytochrome b5 reductase homologue, decreasing cell death. Our work emphasizes the close collaboration of alterations across the genome, transcriptome and microbiome in promoting tumorigenesis.

Pharmacologically targeting the myristoylation of the scaffold protein FRS2α inhibits FGF/FGFR-mediated oncogenic signaling and tumor progression.

Fibroblast growth factor (FGF)/FGF receptor (FGFR) signaling facilitates tumor initiation and progression. Although currently approved inhibitors of FGFR kinase have shown therapeutic benefit in clinical trials, overexpression or mutations of FGFRs eventually confer drug resistance and thereby abrogate the desired activity of kinase inhibitors in many cancer types. In this study, we report that loss of myristoylation of fibroblast growth factor receptor substrate 2 (FRS2α), a scaffold protein essential for FGFR signaling, inhibits FGF/FGFR-mediated oncogenic signaling and FGF10-induced tumorigenesis. Moreover, a previously synthesized myristoyl-CoA analog, B13, which targets the activity of N-myristoyltransferases, suppressed FRS2α myristoylation and decreased the phosphorylation with mild alteration of FRS2α localization at the cell membrane. B13 inhibited oncogenic signaling induced by WT FGFRs or their drug-resistant mutants (FGFRsDRM). B13 alone or in combination with an FGFR inhibitor suppressed FGF-induced WT FGFR- or FGFRDRM-initiated phosphoinositide 3-kinase (PI3K) activity or MAPK signaling, inducing cell cycle arrest and thereby inhibiting cell proliferation and migration in several cancer cell types. Finally, B13 significantly inhibited the growth of xenograft tumors without pathological toxicity to the liver, kidney, or lung in vivo In summary, our study suggests a possible therapeutic approach for inhibiting FGF/FGFR-mediated cancer progression and drug-resistant FGF/FGFR mutants.

Nucleosome positioning changes during human embryonic stem cell differentiation.

Nucleosomes are the basic unit of chromatin. Nucleosome positioning (NP) plays a key role in transcriptional regulation and other biological processes. To better understand NP we used MNase-seq to investigate changes that occur as human embryonic stem cells (hESCs) transition to nascent mesoderm and then to smooth muscle cells (SMCs). Compared to differentiated cell derivatives, nucleosome occupancy at promoters and other notable genic sites, such as exon/intron junctions and adjacent regions, in hESCs shows a stronger correlation with transcript abundance and is less influenced by sequence content. Upon hESC differentiation, genes being silenced, but not genes being activated, display a substantial change in nucleosome occupancy at their promoters. Genome-wide, we detected a shift of NP to regions of higher G+C content as hESCs differentiate to SMCs. Notably, genomic regions with higher nucleosome occupancy harbor twice as many G↔C changes but fewer than half A↔T changes, compared to regions with lower nucleosome occupancy. Finally, our analysis indicates that the hESC genome is not rearranged and has a sequence mutation rate resembling normal human genomes. Our study reveals another unique feature of hESC chromatin, and sheds light on the relationship between nucleosome occupancy and sequence G+C content.

Cell-Cycle Control of Bivalent Epigenetic Domains Regulates the Exit from Pluripotency.

Here we show that bivalent domains and chromosome architecture for bivalent genes are dynamically regulated during the cell cycle in human pluripotent cells. Central to this is the transient increase in H3K4-trimethylation at developmental genes during G1, thereby creating a "window of opportunity" for cell-fate specification. This mechanism is controlled by CDK2-dependent phosphorylation of the MLL2 (KMT2B) histone methyl-transferase, which facilitates its recruitment to developmental genes in G1. MLL2 binding is required for changes in chromosome architecture around developmental genes and establishes promoter-enhancer looping interactions in a cell-cycle-dependent manner. These cell-cycle-regulated loops are shown to be essential for activation of bivalent genes and pluripotency exit. These findings demonstrate that bivalent domains are established to control the cell-cycle-dependent activation of developmental genes so that differentiation initiates from the G1 phase.

Canine spontaneous head and neck squamous cell carcinomas represent their human counterparts at the molecular level.

Spontaneous canine head and neck squamous cell carcinoma (HNSCC) represents an excellent model of human HNSCC but is greatly understudied. To better understand and utilize this valuable resource, we performed a pilot study that represents its first genome-wide characterization by investigating 12 canine HNSCC cases, of which 9 are oral, via high density array comparative genomic hybridization and RNA-seq. The analyses reveal that these canine cancers recapitulate many molecular features of human HNSCC. These include analogous genomic copy number abnormality landscapes and sequence mutation patterns, recurrent alteration of known HNSCC genes and pathways (e.g., cell cycle, PI3K/AKT signaling), and comparably extensive heterogeneity. Amplification or overexpression of protein kinase genes, matrix metalloproteinase genes, and epithelial-mesenchymal transition genes TWIST1 and SNAI1 are also prominent in these canine tumors. This pilot study, along with a rapidly growing body of literature on canine cancer, reemphasizes the potential value of spontaneous canine cancers in HNSCC basic and translational research.

Comprehensive characterization of the genomic alterations in human gastric cancer.

Gastric cancer is one of the most prevalent and aggressive cancers worldwide, and its molecular mechanism remains largely elusive. Here we report the genomic landscape in primary gastric adenocarcinoma of human, based on the complete genome sequences of five pairs of cancer and matching normal samples. In total, 103,464 somatic point mutations, including 407 nonsynonymous ones, were identified and the most recurrent mutations were harbored by Mucins (MUC3A and MUC12) and transcription factors (ZNF717, ZNF595 and TP53). 679 genomic rearrangements were detected, which affect 355 protein-coding genes; and 76 genes show copy number changes. Through mapping the boundaries of the rearranged regions to the folded three-dimensional structure of human chromosomes, we determined that 79.6% of the chromosomal rearrangements happen among DNA fragments in close spatial proximity, especially when two endpoints stay in a similar replication phase. We demonstrated evidences that microhomology-mediated break-induced replication was utilized as a mechanism in inducing ∼40.9% of the identified genomic changes in gastric tumor. Our data analyses revealed potential integrations of Helicobacter pylori DNA into the gastric cancer genomes. Overall a large set of novel genomic variations were detected in these gastric cancer genomes, which may be essential to the study of the genetic basis and molecular mechanism of the gastric tumorigenesis.

Molecular homology and difference between spontaneous canine mammary cancer and human breast cancer.

Spontaneously occurring canine mammary cancer represents an excellent model of human breast cancer, but is greatly understudied. To better use this valuable resource, we performed whole-genome sequencing, whole-exome sequencing, RNA-seq, and/or high-density arrays on twelve canine mammary cancer cases, including seven simple carcinomas and four complex carcinomas. Canine simple carcinomas, which histologically match human breast carcinomas, harbor extensive genomic aberrations, many of which faithfully recapitulate key features of human breast cancer. Canine complex carcinomas, which are characterized by proliferation of both luminal and myoepithelial cells and are rare in human breast cancer, seem to lack genomic abnormalities. Instead, these tumors have about 35 chromatin-modification genes downregulated and are abnormally enriched with active histone modification H4-acetylation, whereas aberrantly depleted with repressive histone modification H3K9me3. Our findings indicate the likelihood that canine simple carcinomas arise from genomic aberrations, whereas complex carcinomas originate from epigenomic alterations, reinforcing their unique value. Canine complex carcinomas offer an ideal system to study myoepithelial cells, the second major cell lineage of the mammary gland. Canine simple carcinomas, which faithfully represent human breast carcinomas at the molecular level, provide indispensable models for basic and translational breast cancer research.

Cancer driver candidate genes AVL9, DENND5A and NUPL1 contribute to MDCK cystogenesis.

AVL9, DENND5A and NUPL1 are among the cancer driver candidate genes previously identified via dog-human comparison, and may function in epithelial cell polarity as indicated by bioinformatics analysis. To better understand their cellular functions and roles in cancer, we knocked down each gene in MDCKII cells through shRNA and performed three-dimensional culture. Compared to the control, the knockdown clones developed significantly more abnormal cysts, e.g., cysts with the lumen harboring dead and/or live cells, or cysts having multiple lumens. Further analysis revealed that abnormalities initiated at the first cell division and persisted throughout the entire cystogenesis process. For NUPL1-knockdown cells, abnormal cytogenesis largely arose from faulty cell divisions, notably monopolar spindles or spindles with poorly separated poles. For AVL9- or DENND5A-knockdown cells, abnormalities originated from both aberrant intracellular trafficking and defective mitosis. Moreover, while all knockdown clones displayed an accelerated rate of both cell proliferation and death, only AVL9- and DENND5A-knockdowns, but not NUPL1-knockdown, promoted cell migration. These observations indicate that NUPL1 contributes to bipolar spindle formation, whereas AVL9 and DENND5A participate in both intracellular trafficking and cell cycle progression. Our study shed lights on these genes' normal cellular functions and on how their alteration contributes to carcinogenesis.