Evolutionary and developmental analysis reveals KANK genes were co-opted for vertebrate vascular development

Evolution of two-pore domain potassium channels and their gene expression in zebrafish embryos

BACKGROUND

The two-pore domain potassium (K2P) channels are a major type of potassium channels that maintain the cell membrane potential by conducting passive potassium leak currents independent of voltage change. They play prominent roles in multiple physiological processes, including neuromodulation, perception of pain, breathing and mood control, and response to volatile anesthetics. Mutations in K2P channels have been linked to many human diseases, such as neuronal and cardiovascular disorders and cancers. Significant progress has been made to understand their protein structures, physiological functions, and pharmacological modifiers. However, their expression and function during embryonic development remain largely unknown.

RESULTS

We employed the zebrafish model and identified 23 k2p genes using BLAST search and gene cloning. We first analyzed vertebrate K2P channel evolution by phylogenetic and syntenic analyses. Our data revealed that the six subtypes of the K2P genes have already evolved in invertebrates long before the emergence of vertebrates. Moreover, the vertebrate K2P gene number increased, most likely due to two whole-genome duplications. Furthermore, we examined zebrafish k2p gene expression during early embryogenesis by in situ hybridization. Each subgroup's genes showed similar but distinct gene expression domains with some exceptions. Most of them were expressed in neural tissues consistent with their known function of neural excitability regulation. However, a few k2p genes were expressed temporarily in specific tissues or organs, suggesting that these K2P channels may be needed for embryonic development.

CONCLUSIONS

Our phylogenetic and developmental analyses of K2P channels shed light on their evolutionary history and potential roles during embryogenesis related to their physiological functions and human channelopathies.

A germline chimeric KANK1-DMRT1 transcript derived from a complex structural variant is associated with a congenital heart defect segregating across five generations

Structural variants (SVs) pose a challenge to detect and interpret, but their study provides novel biological insights and molecular diagnosis underlying rare diseases. The aim of this study was to resolve a 9p24 rearrangement segregating in a family through five generations with a congenital heart defect (congenital pulmonary and aortic valvular stenosis and pulmonary artery stenosis), by applying a combined genomic analysis. The analysis involved multiple techniques, including karyotype, chromosomal microarray analysis (CMA), FISH, genome sequencing (GS), RNA-seq, and optical genome mapping (OGM). A complex 9p24 SV was hinted at by CMA results, showing three interspersed duplicated segments. Combined GS and OGM analyses revealed that the 9p24 duplications constitute a complex SV, on which a set of breakpoints matches the boundaries of the CMA duplicated sequences. The proposed structure for this complex rearrangement implies three duplications associated with an inversion of ~ 2 Mb region on chromosome 9 and a SINE element insertion at the more distal breakpoint. Interestingly, this genomic structure of rearrangement forms a chimeric transcript of the KANK1/DMRT1 loci, which was confirmed by both RNA-seq and Sanger sequencing on blood samples from 9p24 rearrangement carriers. Altogether with breakpoint amplification and FISH analysis, this combined approach allowed a deep characterization of this complex rearrangement. Although the genotype-phenotype correlation remains elusive from the molecular mechanism point of view, this study identified a large genomic rearrangement at 9p24 segregating with a familial congenital heart defect, revealing a genetic biomarker that was successfully applied for embryo selection, changing the reproductive perspective of affected individuals.

A germline chimeric KANK1-DMRT1 transcript derived from a complex structural variant is associated with a congenital heart defect segregating across five generations

Structural variants (SVs) pose a challenge to detect and interpret, but their study provides novel biological insights and molecular diagnosis underlying rare diseases. The aim of this study was to resolve a 9p24 rearrangement segregating in a family through five generations with a congenital heart defect (congenital pulmonary and aortic valvular stenosis, and pulmonary artery stenosis), by applying a combined genomic analysis. The analysis involved multiple techniques, including karyotype, chromosomal microarray analysis (CMA), FISH, whole-genome sequencing (WGS), RNA-seq and optical genome mapping (OGM). A complex 9p24 SV was hinted at by CMA results, showing three interspersed duplicated segments. Combined WGS and OGM analyses revealed that the 9p24 duplications constitute a complex SV, on which a set of breakpoints match the boundaries of the CMA duplicated sequences. The proposed structure for this complex rearrangement implies three duplications associated with an inversion of ~ 2Mb region on chromosome 9 with a SINE element insertion at the more distal breakpoint. Interestingly, this hypothesized genomic structure of rearrangement forms a chimeric transcript of the KANK1/DMRT1 loci, which was confirmed by RNA-seq on blood from 9p24 rearrangement carriers. Altogether with breakpoint amplification and FISH analysis, this combined approach allowed a deep characterization of this complex rearrangement. Although the genotype-phenotype correlation remains elusive from the molecular mechanism point of view, this study identified a large genomic rearrangement at 9p segregating with a familial congenital clinical trait, revealing a genetic biomarker that was successfully applied for embryo selection, changing the reproductive perspective of affected individuals.

KANK4 Promotes Arteriogenesis by Potentiating VEGFR2 Signaling in a TALIN-1-Dependent Manner

BACKGROUND

Arteriogenesis plays a critical role in maintaining adequate tissue blood supply and is related to a favorable prognosis in arterial occlusive diseases. Strategies aimed at promoting arteriogenesis have thus far not been successful because the factors involved in arteriogenesis remain incompletely understood. Previous studies suggest that evolutionarily conserved KANK4 (KN motif and ankyrin repeat domain-containing proteins 4) might involve in vertebrate vessel development. However, how the KANK4 regulates vessel function remains unknown. We aim to determine the role of endothelial cell-specifically expressed KANK4 in arteriogenesis.

METHODS

The role of KANK4 in regulating arteriogenesis was evaluated using Kank4^-/- and KANK4^iECOE mice. Molecular mechanisms underlying KANK4-potentiated arteriogenesis were investigated by employing RNA transcriptomic profiling and mass spectrometry analysis.

RESULTS

By analyzing Kank4-EGFP reporter mice, we showed that KANK4 was specifically expressed in endothelial cells. In particular, KANK4 displayed a dynamic expression pattern from being ubiquitously expressed in all endothelial cells of the developing vasculature to being explicitly expressed in the endothelial cells of arterioles and arteries in matured vessels. In vitro microfluidic chip-based vascular morphology analysis and in vivo hindlimb ischemia assays using Kank4^-/- and KANK4^iECOE mice demonstrated that deletion of KANK4 impaired collateral artery growth and the recovery of blood perfusion, whereas KANK4 overexpression leads to increased vessel caliber and blood perfusion. Bulk RNA sequencing and Co-immunoprecipitation/mass spectrometry (Co-IP/MS) analysis identified that KANK4 promoted EC proliferation and collateral artery remodeling through coupling VEGFR2 (vascular endothelial growth factor receptor 2) to TALIN-1, which augmented the activation of the VEGFR2 signaling cascade.

CONCLUSIONS

This study reveals a novel role for KANK4 in arteriogenesis in response to ischemia. KANK4 links VEGFR2 to TALIN-1, resulting in enhanced VEGFR2 activation and increased EC proliferation, highlighting that KANK4 is a potential therapeutic target for promoting arteriogenesis for arterial occlusive diseases.

Proteomic profiling of concurrently isolated primary microvascular endothelial cells, pericytes, and vascular smooth muscle cells from adult mouse heart

The microcirculation serves crucial functions in adult heart, distinct from those carried out by epicardial vessels. Microvessels are governed by unique regulatory mechanisms, impairment of which leads to microvessel-specific pathology. There are few treatment options for patients with microvascular heart disease, primarily due to limited understanding of underlying pathology. High throughput mRNA sequencing and protein expression profiling in specific cells can improve our understanding of microvessel biology and disease at the molecular level. Understanding responses of individual microvascular cells to the same physiological or pathophysiological stimuli requires the ability to isolate the specific cell types that comprise the functional units of the microcirculation in the heart, preferably from the same heart, to ensure that different cells have been exposed to the same in-vivo conditions. We developed an integrated process for simultaneous isolation and culture of the main cell types comprising the microcirculation in adult mouse heart: endothelial cells, pericytes, and vascular smooth muscle cells. These cell types were characterized with isobaric labeling quantitative proteomics and mRNA sequencing. We defined microvascular cell proteomes, identified novel protein markers, and confirmed established cell-specific markers. Our results allow identification of unique markers and regulatory proteins that govern microvascular physiology and pathology.

Ancient Origins of Cytoskeletal Crosstalk: Spectraplakin-like Proteins Precede the Emergence of Cortical Microtubule Stabilization Complexes as Crosslinkers

Adhesion between cells and the extracellular matrix (ECM) is one of the prerequisites for multicellularity, motility, and tissue specialization. Focal adhesions (FAs) are defined as protein complexes that mediate signals from the ECM to major components of the cytoskeleton (microtubules, actin, and intermediate filaments), and their mutual communication determines a variety of cellular processes. In this study, human cytoskeletal crosstalk proteins were identified by comparing datasets with experimentally determined cytoskeletal proteins. The spectraplakin dystonin was the only protein found in all datasets. Other proteins (FAK, RAC1, septin 9, MISP, and ezrin) were detected at the intersections of FAs, microtubules, and actin cytoskeleton. Homology searches for human crosstalk proteins as queries were performed against a predefined dataset of proteomes. This analysis highlighted the importance of FA communication with the actin and microtubule cytoskeleton, as these crosstalk proteins exhibit the highest degree of evolutionary conservation. Finally, phylogenetic analyses elucidated the early evolutionary history of spectraplakins and cortical microtubule stabilization complexes (CMSCs) as model representatives of the human cytoskeletal crosstalk. While spectraplakins probably arose at the onset of opisthokont evolution, the crosstalk between FAs and microtubules is associated with the emergence of metazoans. The multiprotein complexes contributing to cytoskeletal crosstalk in animals gradually gained in complexity from the onset of metazoan evolution.

The functional genetic architecture of egg-laying and live-bearing reproduction in common lizards

All amniotes reproduce either by egg-laying (oviparity), which is ancestral to vertebrates or by live-bearing (viviparity), which has evolved many times independently. However, the genetic basis of these parity modes has never been resolved and, consequently, its convergence across evolutionary scales is currently unknown. Here, we leveraged natural hybridizations between oviparous and viviparous common lizards (Zootoca vivipara) to describe the functional genes and genetic architecture of parity mode and its key traits, eggshell and gestation length, and compared our findings across vertebrates. In these lizards, parity trait genes were associated with progesterone-binding functions and enriched for tissue remodelling and immune system pathways. Viviparity involved more genes and complex gene networks than did oviparity. Angiogenesis, vascular endothelial growth and adrenoreceptor pathways were enriched in the viviparous female reproductive tissue, while pathways for transforming growth factor were enriched in the oviparous. Natural selection on these parity mode genes was evident genome-wide. Our comparison to seven independent origins of viviparity in mammals, squamates and fish showed that genes active in pregnancy were related to immunity, tissue remodelling and blood vessel generation. Therefore, our results suggest that pre-established regulatory networks are repeatedly recruited for viviparity and that these are shared at deep evolutionary scales.

Loss of smarcad1a accelerates tumorigenesis of malignant peripheral nerve sheath tumors in zebrafish

Malignant peripheral nerve sheath tumors (MPNSTs) are a type of sarcoma that generally originates from Schwann cells. The prognosis for this type of malignancy is relatively poor due to complicated genetic alterations and the lack of specific targeted therapy. Chromosome fragment 4q22-23 is frequently deleted in MPNSTs and other human tumors, suggesting tumor suppressor genes may reside in this region. Here, we provide evidence that SMARCAD1, a known chromatin remodeler, is a novel tumor suppressor gene located in 4q22-23. We identified two human homologous smarcad1 genes (smarcad1a and smarcad1b) in zebrafish, and both genes share overlapping expression patterns during embryonic development. We demonstrated that two smarcad1a loss-of-function mutants, sa1299 and p403, can accelerate MPNST tumorigenesis in the tp53 mutant background, suggesting smarcad1a is a bona fide tumor suppressor gene for MPNSTs. Moreover, we found that DNA double-strand break (DSB) repair might be compromised in both mutants compared to wildtype zebrafish, as indicated by pH2AX, a DNA DSB marker. In addition, both SMARCAD1 gene knockdown and overexpression in human cells were able to inhibit tumor growth and displayed similar DSB repair responses, suggesting proper SMARCAD1 gene expression level or gene dosage is critical for cell growth. Given that mutations of SMARCAD1 sensitize cells to poly ADP ribose polymerase inhibitors in yeast and the human U2OS osteosarcoma cell line, the identification of SMARCAD1 as a novel tumor suppressor gene might contribute to the development of new cancer therapies for MPNSTs.

Evolution of inwardly rectifying potassium channels and their gene expression in zebrafish embryos

BACKGROUND

Inwardly rectifying potassium channels are essential for normal potassium homeostasis, maintaining the cellular resting membrane potential, and regulating electrolyte transportation. Mutations in Kir channels have been known to cause debilitating diseases ranging from neurological abnormalities to renal and cardiac failures. Many efforts have been made to understand their protein structures, physiological functions, and pharmacological modifiers. However, their expression and functions during embryonic development remain largely unknown.

RESULTS

Using zebrafish as a model, we identified and renamed 31 kir genes. We also analyzed Kir gene evolution by phylogenetic and syntenic analyses. Our data indicated that the four subtypes of the Kir genes might have already evolved out in chordates. These vertebrate Kir genes most likely resulted from both whole-genome duplications and tandem duplications. In addition, we examined zebrafish kir gene expression during early embryogenesis. Each subgroup's genes showed similar but distinct gene expression domains. The gene expression of ohnologous genes from teleost-specific whole-genome duplication indicated subfunctionalization. Varied temporal gene expression domains suggest that Kir channels may be needed for embryonic patterning or regulation.

CONCLUSIONS

Our phylogenetic and developmental analyses of Kir channels shed light on their evolutionary history and potential functions during embryogenesis related to congenital diseases and human channelopathies.

Phylogenetic and developmental analyses indicate complex functions of calcium-activated potassium channels in zebrafish embryonic development

BACKGROUND

Calcium-activated potassium channels (KCa) are a specific type of potassium channel activated by intracellular calcium concentration changes. This group of potassium channels plays fundamental roles ranging from regulating neuronal excitability to immune cell activation. Many human diseases such as schizophrenia, hypertension, epilepsy, and cancers have been linked to mutations in this group of potassium channels. Although the KCa channels have been extensively studied electrophysiologically and pharmacologically, their spatiotemporal gene expression during embryogenesis remains mostly unknown.

RESULTS

Using zebrafish as a model, we identified and renamed 14 KCa genes. We further performed phylogenetic and syntenic analyses on vertebrate KCa genes. Our data revealed that the number of KCa genes in zebrafish was increased, most likely due to teleost-specific whole-genome duplication. Moreover, we examined zebrafish KCa gene expression during early embryogenesis. The duplicated ohnologous genes show distinct and overlapped gene expression. Furthermore, we found that zebrafish KCa genes are expressed in various tissues and organs (somites, fins, olfactory regions, eye, kidney, and so on) and neuronal tissues, suggesting that they may play important roles during zebrafish embryogenesis.

CONCLUSIONS

Our phylogenetic and developmental analyses shed light on the potential functions of the KCa genes during embryogenesis related to congenital diseases and human channelopathies.

Tissue distribution and subcellular localization of the family of Kidney Ankyrin Repeat Domain (KANK) proteins

Kidney Ankyrin Repeat-containing Proteins (KANKs) comprise a family of four evolutionary conserved proteins (KANK1 to 4) that localize to the belt of mature focal adhesions (FAs) where they regulate integrin-mediated adhesion, actomyosin contractility, and link FAs to the cortical microtubule stabilization complex (CMSC). The human KANK proteins were first identified in kidney and have been associated with kidney cancer and nephrotic syndrome. Here, we report the distributions and subcellular localizations of the four Kank mRNAs and proteins in mouse tissues. We found that the KANK family members display distinct and rarely overlapping expression patterns. Whereas KANK1 is expressed at the basal side of epithelial cells of all tissues tested, KANK2 expression is mainly observed at the plasma membrane and/or cytoplasm of mesenchymal cells and KANK3 exclusively in vascular and lymphatic endothelial cells. KANK4 shows the least widespread expression pattern and when present, overlaps with KANK2 in contractile cells, such as smooth muscle cells and pericytes. Our findings show that KANKs are widely expressed in a cell type-specific manner, which suggests that they have cell- and tissue-specific functions.

Potassium Channel-Associated Bioelectricity of the Dermomyotome Determines Fin Patterning in Zebrafish

The roles of bioelectric signaling in developmental patterning remain largely unknown, although recent work has implicated bioelectric signals in cellular processes such as proliferation and migration. Here, we report a mutation in the inwardly rectifying potassium channel (kir) gene, kcnj13/kir7.1, that causes elongation of the fins in the zebrafish insertional mutant Dhi2059. A viral DNA insertion into the noncoding region of kcnj13 results in transient activation and ectopic expression of kcnj13 in the somite and dermomyotome, from which the fin ray progenitors originate. We made an allele-specific loss-of-function kcnj13 mutant by CRISPR (clustered regularly interspaced short palindromic repeats) and showed that it could reverse the long-finned phenotype, but only when located on the same chromosome as the Dhi2059 viral insertion. Also, we showed that ectopic expression of kcnj13 in the dermomyotome of transgenic zebrafish produces phenocopies of the Dhi2059 mutant in a gene dosage-sensitive manner. Finally, to determine whether this developmental function is specific to kcnj13, we ectopically expressed three additional potassium channel genes: kcnj1b, kcnj10a, and kcnk9 We found that all induce the long-finned phenotype, indicating that this function is conserved among potassium channel genes. Taken together, our results suggest that dermomyotome bioelectricity is a new fin-patterning mechanism, and we propose a two-stage bioelectricity model for zebrafish fin patterning. This ion channel-regulated bioelectric developmental patterning mechanism may provide with us new insight into vertebrate morphological evolution and human congenital malformations.

Identification of RECK as an evolutionarily conserved tumor suppressor gene for zebrafish malignant peripheral nerve sheath tumors

Malignant peripheral nerve sheath tumors (MPNSTs) are a type of sarcoma with poor prognosis due to their complex genetic changes, invasive growth, and insensitivity to chemo- and radiotherapies. One of the most frequently lost chromosome arms in human MPNSTs is chromosome 9p. However, the cancer driver genes located on it remain largely unknown, except the tumor suppressor gene, p16 (INK4)/CDKN2A. Previously, we identified RECK as a tumor suppressor gene candidate on chromosome 9p using zebrafish-human comparative oncogenomics. In this study, we investigated the tumorigenesis of the reck gene using zebrafish genetic models in both tp53 and ribosomal protein gene mutation background. We also examined the biological effects of RECK gene restoration in human MPNST cell lines. These results provide the first genetic evidence that reck is a bona fide tumor suppressor gene for MPNSTs in zebrafish. In addition, restoration of the RECK gene in human MPNST cells leads to growth inhibition suggesting that the reactivation of RECK could serve as a potential therapeutic strategy for MPNSTs.

Evolutionary and Expression Analyses Show Co-option of khdrbs Genes for Origin of Vertebrate Brain

Genes generated by whole genome duplications (WGD) can be co-opted by changing their regulation process or altering their coding proteins, which has been shown contributable to the emergence of vertebrate morphological novelties such as vertebrate cartilage. Mouse khdrbs genes, differing from its invertebrate orthologs, were mainly expressed in brain, hinting that khdrbs gene family as a member of genetic toolkit may be linked to vertebrate brain development. However, the evolutionary relationship between khdrbs gene family and vertebrate brain development is unclear. First, we analyzed the evolutionary history of khdrbs gene family in metazoans, and then investigated their expression patterns during early development and in adulthood of zebrafish. We found that the duplication of khdrbs gene family by WGD took place in zebrafish, and all zebrafish khdrbs genes were predominantly expressed in the substructures of brain during early development. Given the expression of invertebrate khdrbs gene in germ line, the distinct expression domains of zebrafish khdrbs genes in brain suggested that the duplicated khdrbs genes are co-opted for promoting the evolutionary origin of vertebrate brain.

KANK1 inhibits cell growth by inducing apoptosis through regulating CXXC5 in human malignant peripheral nerve sheath tumors

Malignant peripheral nerve sheath tumors (MPNSTs) are a type of rare sarcomas with a poor prognosis due to its highly invasive nature and limited treatment options. Currently there is no targeted-cancer therapy for this type of malignancy. Thus, it is important to identify more cancer driver genes that may serve as targets of cancer therapy. Through comparative oncogenomics, we have found that KANK1 was a candidate tumor suppressor gene (TSG) for human MPNSTs. Although KANK1 is known as a cytoskeleton regulator, its tumorigenic function in MPNSTs remains largely unknown. In this study, we report that restoration of KANK1 in human MPNST cells inhibits cell growth both in human cell culture and xenograft mice by increasing apoptosis. Consistently, knockdown of KANK1 in neurofibroma cells promoted cell growth. Using RNA-seq analysis, we identified CXXC5 and other apoptosis-related genes, and demonstrated that CXXC5 is regulated by KANK1. Knockdown of CXXC5 was found to diminish KANK1-induced apoptosis in MPNST cells. Thus, KANK1 inhibits MPNST cell growth though CXXC5 mediated apoptosis. Our results suggest that KANK1 may function as a tumor suppressor in human MPNSTs, and thus it may be useful for targeted therapy.

Molecular Evolution of MDM1, a "Duplication-Resistant" Gene in Vertebrates

BACKGROUND

The mouse double minute 1 (Mdm1) gene was first reported and cloned in mouse tumor cell lines as an oncogene candidate. Later, it was found that mutation of Mdm1 might cause age-related retinal degeneration 2 in mice by genetic linkage analysis. Additionally, the MDM1 protein was found to be expressed in the centrosomes, cilia, and the nucleus of multiciliated tracheal epithelial cells in mice. These observations suggest that MDM1 may have some basal functions in cell physiology. However, the evolutionary history of this gene and its expression during embryonic development remain largely unexplored.

RESULTS

Using molecular phylogenetic analysis, we found that the MDM1 gene encoded an evolutionarily conserved protein across all metazoans. We also found that the MDM1 gene was in a conserved synteny in vertebrates. In almost all the species that were analyzed, there was only one MDM1 gene based on current genome annotations. Since vertebrate genomes underwent two to three rounds of whole-genome duplications around the origin of the vertebrates, it is interesting that only one MDM1 ohnolog was retained. This observation implies that other MDM1 ohnologs were lost after the whole-genome duplications. Furthermore, using whole-mount in situ hybridization, we found that mdm1 was expressed in the forebrain, nephric ducts, and tail buds during zebrafish early embryonic development.

CONCLUSION

MDM1 is an evolutionary conserved gene, and its homologous genes can be traced back to basal metazoan lineages. In vertebrates, the MDM1 gene is in a conserved synteny and there is only one MDM1 ohnolog suggesting it is a "duplication-resistant" gene. Its expression patterns in early zebrafish embryos indicate that mdm1 may play important roles in the development of the central nervous system, kidneys, and hematopoietic system.