Prioritizing cardiovascular disease-associated variants altering NKX2-5 and TBX5 binding through an integrative computational approach.

Cardiovascular diseases (CVDs) are the leading cause of death worldwide and are heavily influenced by genetic factors. Genome-wide association studies have mapped >90% of CVD-associated variants within the noncoding genome, which can alter the function of regulatory proteins, such as transcription factors (TFs). However, due to the overwhelming number of single-nucleotide polymorphisms (SNPs) (>500,000) in genome-wide association studies, prioritizing variants for in vitro analysis remains challenging. In this work, we implemented a computational approach that considers support vector machine (SVM)-based TF binding site classification and cardiac expression quantitative trait loci (eQTL) analysis to identify and prioritize potential CVD-causing SNPs. We identified 1535 CVD-associated SNPs within TF footprints and putative cardiac enhancers plus 14,218 variants in linkage disequilibrium with genotype-dependent gene expression in cardiac tissues. Using ChIP-seq data from two cardiac TFs (NKX2-5 and TBX5) in human-induced pluripotent stem cell-derived cardiomyocytes, we trained a large-scale gapped k-mer SVM model to identify CVD-associated SNPs that altered NKX2-5 and TBX5 binding. The model was tested by scoring human heart TF genomic footprints within putative enhancers and measuring in vitro binding through electrophoretic mobility shift assay. Five variants predicted to alter NKX2-5 (rs59310144, rs6715570, and rs61872084) and TBX5 (rs7612445 and rs7790964) binding were prioritized for in vitro validation based on the magnitude of the predicted change in binding and are in cardiac tissue eQTLs. All five variants altered NKX2-5 and TBX5 DNA binding. We present a bioinformatic approach that considers tissue-specific eQTL analysis and SVM-based TF binding site classification to prioritize CVD-associated variants for in vitro analysis.

Prioritizing Cardiovascular Disease-Associated Variants Altering NKX2-5 Binding through an Integrative Computational Approach.

Cardiovascular diseases (CVDs) are the leading cause of death worldwide and are heavily influenced by genetic factors. Genome-wide association studies (GWAS) have mapped > 90% of CVD-associated variants within the non-coding genome, which can alter the function of regulatory proteins, like transcription factors (TFs). However, due to the overwhelming number of GWAS single nucleotide polymorphisms (SNPs) (>500,000), prioritizing variants for in vitro analysis remains challenging. In this work, we implemented a computational approach that considers support vector machine (SVM)-based TF binding site classification and cardiac expression quantitative trait loci (eQTL) analysis to identify and prioritize potential CVD-causing SNPs. We identified 1,535 CVD-associated SNPs that occur within human heart footprints/enhancers and 9,309 variants in linkage disequilibrium (LD) with differential gene expression profiles in cardiac tissue. Using hiPSC-CM ChIP-seq data from NKX2-5 and TBX5, two cardiac TFs essential for proper heart development, we trained a large-scale gapped k-mer SVM (LS-GKM-SVM) predictive model that can identify binding sites altered by CVD-associated SNPs. The computational predictive model was tested by scoring human heart footprints and enhancers in vitro through electrophoretic mobility shift assay (EMSA). Three variants (rs59310144, rs6715570, and rs61872084) were prioritized for in vitro validation based on their eQTL in cardiac tissue and LS-GKM-SVM prediction to alter NKX2-5 DNA binding. All three variants altered NKX2-5 DNA binding. In summary, we present a bioinformatic approach that considers tissue-specific eQTL analysis and SVM-based TF binding site classification to prioritize CVD-associated variants for in vitro experimental analysis.

Characterization and Expression of Holothurian Wnt Signaling Genes during Adult Intestinal Organogenesis.

Wnt signaling has been shown to play multiple roles in regenerative processes, one of the most widely studied of which is the regeneration of the intestinal luminal epithelia. Most studies in this area have focused on self-renewal of the luminal stem cells; however, Wnt signaling may also have more dynamic functions, such as facilitating intestinal organogenesis. To explore this possibility, we employed the sea cucumber Holothuria glaberrima that can regenerate a full intestine over the course of 21 days after evisceration. We collected RNA-seq data from various intestinal tissues and regeneration stages and used these data to define the Wnt genes present in H. glaberrima and the differential gene expression (DGE) patterns during the regenerative process. Twelve Wnt genes were found, and their presence was confirmed in the draft genome of H. glaberrima. The expressions of additional Wnt-associated genes, such as Frizzled and Disheveled, as well as genes from the Wnt/ß-catenin and Wnt/Planar Cell Polarity (PCP) pathways, were also analyzed. DGE showed unique distributions of Wnt in early- and late-stage intestinal regenerates, consistent with the Wnt/ß-catenin pathway being upregulated during early-stages and the Wnt/PCP pathway being upregulated during late-stages. Our results demonstrate the diversity of Wnt signaling during intestinal regeneration, highlighting possible roles in adult organogenesis.

Characterization of Two Novel EF-Hand Proteins Identifies a Clade of Putative Ca2+-Binding Protein Specific to the Ambulacraria.

In recent years, transcriptomic databases have become one of the main sources for protein discovery. In our studies of nervous system and digestive tract regeneration in echinoderms, we have identified several transcripts that have attracted our attention. One of these molecules corresponds to a previously unidentified transcript (Orpin) from the sea cucumber Holothuria glaberrima that appeared to be upregulated during intestinal regeneration. We have now identified a second highly similar sequence and analyzed the predicted proteins using bioinformatics tools. Both sequences have EF-hand motifs characteristic of calcium-binding proteins (CaBPs) and N-terminal signal peptides. Sequence comparison analyses such as multiple sequence alignments and phylogenetic analyses only showed significant similarity to sequences from other echinoderms or from hemichordates. Semi-quantitative RT-PCR analyses revealed that transcripts from these sequences are expressed in various tissues including muscle, haemal system, gonads, and mesentery. However, contrary to previous reports, there was no significant differential expression in regenerating tissues. Nonetheless, the identification of unique features in the predicted proteins and their presence in the holothurian draft genome suggest that these might comprise a novel subfamily of EF-hand containing proteins specific to the Ambulacraria clade.

Transcriptomic analysis of early stages of intestinal regeneration in Holothuria glaberrima.

Echinoderms comprise a group of animals with impressive regenerative capabilities. They can replace complex internal organs following injury or autotomy. In holothurians or sea cucumbers, cellular processes of intestinal regeneration have been extensively studied. The molecular machinery behind this faculty, however, remains to be understood. Here we assembled and annotated a de novo transcriptome using RNA-seq data consisting of regenerating and non-regenerating intestinal tissues from the sea cucumber Holothuria glaberrima. Comparisons of differential expression were made using the mesentery as a reference against 24 h and 3 days regenerating intestine, revealing a large number of differentially expressed transcripts. Gene ontology and pathway enrichment analysis showed evidence of increasing transcriptional activity. Further analysis of transcripts associated with transcription factors revealed diverse expression patterns with mechanisms involving developmental and cancer-related activity that could be related to the regenerative process. Our study demonstrates the broad and diversified gene expression profile during the early stages of the process using the mesentery as the focal point of intestinal regeneration. It also establishes the genes that are the most important candidates in the cellular processes that underlie regenerative responses.

Regeneration in Echinoderms: Molecular Advancements.

Which genes and gene signaling pathways mediate regenerative processes? In recent years, multiple studies, using a variety of animal models, have aimed to answer this question. Some answers have been obtained from transcriptomic and genomic studies where possible gene and gene pathway candidates thought to be involved in tissue and organ regeneration have been identified. Several of these studies have been done in echinoderms, an animal group that forms part of the deuterostomes along with vertebrates. Echinoderms, with their outstanding regenerative abilities, can provide important insights into the molecular basis of regeneration. Here we review the available data to determine the genes and signaling pathways that have been proposed to be involved in regenerative processes. Our analyses provide a curated list of genes and gene signaling pathways and match them with the different cellular processes of the regenerative response. In this way, the molecular basis of echinoderm regenerative potential is revealed, and is available for comparisons with other animal taxa.

Draft Genome of the Sea Cucumber Holothuria glaberrima, a Model for the Study of Regeneration.

Regeneration is one of the most fascinating and yet least understood biological processes. Echinoderms, one of the closest related invertebrate groups to humans, can contribute to our understanding of the genetic basis of regenerative processes. Among echinoderms, sea cucumbers have the ability to grow back most of their body parts following injury, including the intestine and nervous tissue. The cellular and molecular events underlying these abilities in sea cucumbers have been most extensively studied in the species Holothuria glaberrima. However, research into the regenerative abilities of this species has been impeded due to the lack of adequate genomic resources. Here, we report the first draft genome assembly of H. glaberrima and demonstrate its value for future genetic studies. Using only short sequencing reads, we assembled the genome into 89,105 scaffolds totaling 1.1 gigabases with an N50 of 25 kilobases. Our BUSCO assessment of the genome resulted in 894 (91.4%) complete and partial genes from 978 genes queried. We incorporated transcriptomic data from several different life history stages to annotate 51,415 genes in our final assembly. To demonstrate the usefulness of the genome, we fully annotated the melanotransferrin (Mtf) gene family, which have a potential role in the regeneration of the sea cucumber intestine. Using these same data, we extracted the mitochondrial genome, which showed high conservation to that of other holothuroids. Thus, these data will be a critical resource for ongoing studies of regeneration and other studies in sea cucumbers.