The pale spear-nosed bat: A neuromolecular and transgenic model for vocal learning.

Vocal learning, the ability to produce modified vocalizations via learning from acoustic signals, is a key trait in the evolution of speech. While extensively studied in songbirds, mammalian models for vocal learning are rare. Bats present a promising study system given their gregarious natures, small size, and the ability of some species to be maintained in captive colonies. We utilize the pale spear-nosed bat (Phyllostomus discolor) and report advances in establishing this species as a tractable model for understanding vocal learning. We have taken an interdisciplinary approach, aiming to provide an integrated understanding across genomics (Part I), neurobiology (Part II), and transgenics (Part III). In Part I, we generated new, high-quality genome annotations of coding genes and noncoding microRNAs to facilitate functional and evolutionary studies. In Part II, we traced connections between auditory-related brain regions and reported neuroimaging to explore the structure of the brain and gene expression patterns to highlight brain regions. In Part III, we created the first successful transgenic bats by manipulating the expression of FoxP2, a speech-related gene. These interdisciplinary approaches are facilitating a mechanistic and evolutionary understanding of mammalian vocal learning and can also contribute to other areas of investigation that utilize P. discolor or bats as study species.

ARHGEF39, a Gene Implicated in Developmental Language Disorder, Activates RHOA and Is Involved in Cell De-Adhesion and Neural Progenitor Cell Proliferation.

ARHGEF39 was previously implicated in developmental language disorder (DLD) via a functional polymorphism that can disrupt post-transcriptional regulation by microRNAs. ARHGEF39 is part of the family of Rho guanine nucleotide exchange factors (RhoGEFs) that activate small Rho GTPases to regulate a wide variety of cellular processes. However, little is known about the function of ARHGEF39, or how its function might contribute to neurodevelopment or related disorders. Here, we explore the molecular function of ARHGEF39 and show that it activates the Rho GTPase RHOA and that high ARHGEF39 expression in cell cultures leads to an increase of detached cells. To explore its role in neurodevelopment, we analyse published single cell RNA-sequencing data and demonstrate that ARHGEF39 is a marker gene for proliferating neural progenitor cells and that it is co-expressed with genes involved in cell division. This suggests a role for ARHGEF39 in neurogenesis in the developing brain. The co-expression of ARHGEF39 with other RHOA-regulating genes supports RHOA as substrate of ARHGEF39 in neural cells, and the involvement of RHOA in neuropsychiatric disorders highlights a potential link between ARHGEF39 and neurodevelopment and disorder. Understanding the GTPase substrate, co-expression network, and processes downstream of ARHGEF39 provide new avenues for exploring the mechanisms by which altered expression levels of ARHGEF39 may contribute to neurodevelopment and associated disorders.

DNA methylation predicts age and provides insight into exceptional longevity of bats.

Exceptionally long-lived species, including many bats, rarely show overt signs of aging, making it difficult to determine why species differ in lifespan. Here, we use DNA methylation (DNAm) profiles from 712 known-age bats, representing 26 species, to identify epigenetic changes associated with age and longevity. We demonstrate that DNAm accurately predicts chronological age. Across species, longevity is negatively associated with the rate of DNAm change at age-associated sites. Furthermore, analysis of several bat genomes reveals that hypermethylated age- and longevity-associated sites are disproportionately located in promoter regions of key transcription factors (TF) and enriched for histone and chromatin features associated with transcriptional regulation. Predicted TF binding site motifs and enrichment analyses indicate that age-related methylation change is influenced by developmental processes, while longevity-related DNAm change is associated with innate immunity or tumorigenesis genes, suggesting that bat longevity results from augmented immune response and cancer suppression.

Six reference-quality genomes reveal evolution of bat adaptations.

Bats possess extraordinary adaptations, including flight, echolocation, extreme longevity and unique immunity. High-quality genomes are crucial for understanding the molecular basis and evolution of these traits. Here we incorporated long-read sequencing and state-of-the-art scaffolding protocols1 to generate, to our knowledge, the first reference-quality genomes of six bat species (Rhinolophus ferrumequinum, Rousettus aegyptiacus, Phyllostomus discolor, Myotis myotis, Pipistrellus kuhlii and Molossus molossus). We integrated gene projections from our 'Tool to infer Orthologs from Genome Alignments' (TOGA) software with de novo and homology gene predictions as well as short- and long-read transcriptomics to generate highly complete gene annotations. To resolve the phylogenetic position of bats within Laurasiatheria, we applied several phylogenetic methods to comprehensive sets of orthologous protein-coding and noncoding regions of the genome, and identified a basal origin for bats within Scrotifera. Our genome-wide screens revealed positive selection on hearing-related genes in the ancestral branch of bats, which is indicative of laryngeal echolocation being an ancestral trait in this clade. We found selection and loss of immunity-related genes (including pro-inflammatory NF-κB regulators) and expansions of anti-viral APOBEC3 genes, which highlights molecular mechanisms that may contribute to the exceptional immunity of bats. Genomic integrations of diverse viruses provide a genomic record of historical tolerance to viral infection in bats. Finally, we found and experimentally validated bat-specific variation in microRNAs, which may regulate bat-specific gene-expression programs. Our reference-quality bat genomes provide the resources required to uncover and validate the genomic basis of adaptations of bats, and stimulate new avenues of research that are directly relevant to human health and disease1.

Tissue Collection of Bats for -Omics Analyses and Primary Cell Culture.

As high-throughput sequencing technologies advance, standardized methods for high quality tissue acquisition and preservation allow for the extension of these methods to non-model organisms. A series of protocols to optimize tissue collection from bats has been developed for a series of high-throughput sequencing approaches. Outlined here are protocols for the capture of bats, desired demographics to be collected for each bat, and optimized methods to minimize stress on a bat during tissue collection. Specifically outlined are methods for collecting and treating tissue to obtain (i) DNA for high molecular weight genomic analyses, (ii) RNA for tissue-specific transcriptomes, and (iii) proteins for proteomic-level analyses. Lastly, also outlined is a method to avoid lethal sampling by creating viable primary cell cultures from wing clips. A central motivation of these methods is to maximize the amount of potential molecular and morphological data for each bat and suggest optimal ways to preserve tissues so they retain their value as new methods develop in the future. This standardization has become particularly important as initiatives to sequence chromosome-level, error-free genomes of species across the world have emerged, in which multiple scientific parties are spearheading the sequencing of different taxonomic groups. The protocols outlined herein define the ideal tissue collection and tissue preservation methods for Bat1K, the consortium that is sequencing the genomes of every species of bat.

P300/CBP-associated factor regulates transcription and function of isocitrate dehydrogenase 2 during muscle differentiation.

The epigenetic enzyme p300/CBP-associated factor (PCAF) belongs to the GCN5-related N-acetyltransferase (GNAT) family together with GCN5. Although its transcriptional and post-translational function is well characterized, little is known about its properties as regulator of cell metabolism. Here, we report the mitochondrial localization of PCAF conferred by an 85 aa mitochondrial targeting sequence (MTS) at the N-terminal region of the protein. In mitochondria, one of the PCAF targets is the isocitrate dehydrogenase 2 (IDH2) acetylated at lysine 180. This PCAF-regulated post-translational modification might reduce IDH2 affinity for isocitrate as a result of a conformational shift involving predictively the tyrosine at position 179. Site-directed mutagenesis and functional studies indicate that PCAF regulates IDH2, acting at dual level during myoblast differentiation: at a transcriptional level together with MyoD, and at a post-translational level by direct modification of lysine acetylation in mitochondria. The latter event determines a decrease in IDH2 function with negative consequences on muscle fiber formation in C2C12 cells. Indeed, a MTS-deprived PCAF does not localize into mitochondria, remains enriched into the nucleus, and contributes to a significant increase of muscle-specific gene expression enhancing muscle differentiation. The role of PCAF in mitochondria is a novel finding shedding light on metabolic processes relevant to early muscle precursor differentiation.-Savoia, M., Cencioni, C., Mori, M., Atlante, S., Zaccagnini, G., Devanna, P., Di Marcotullio, L., Botta, B., Martelli, F., Zeiher, A. M., Pontecorvi, A., Farsetti, A., Spallotta, F., Gaetano, C. P300/CBP-associated factor regulates transcription and function of isocitrate dehydrogenase 2 during muscle differentiation.

Genome-wide investigation of an ID cohort reveals de novo 3'UTR variants affecting gene expression.

Intellectual disability (ID) is a severe neurodevelopmental disorder with genetically heterogeneous causes. Large-scale sequencing has led to the identification of many gene-disrupting mutations; however, a substantial proportion of cases lack a molecular diagnosis. As such, there remains much to uncover for a complete understanding of the genetic underpinnings of ID. Genetic variants present in non-coding regions of the genome have been highlighted as potential contributors to neurodevelopmental disorders given their role in regulating gene expression. Nevertheless the functional characterization of non-coding variants remains challenging. We describe the identification and characterization of de novo non-coding variation in 3'UTR regulatory regions within an ID cohort of 50 patients. This cohort was previously screened for structural and coding pathogenic variants via CNV, whole exome and whole genome analysis. We identified 44 high-confidence single nucleotide non-coding variants within the 3'UTR regions of these 50 genomes. Four of these variants were located within predicted miRNA binding sites and were thus hypothesised to have regulatory consequences. Functional testing showed that two of the variants interfered with miRNA-mediated regulation of their target genes, AMD1 and FAIM. Both these variants were found in the same individual and their functional consequences may point to a potential role for such variants in intellectual disability.

Mapping of Human FOXP2 Enhancers Reveals Complex Regulation.

Mutations of the FOXP2 gene cause a severe speech and language disorder, providing a molecular window into the neurobiology of language. Individuals with FOXP2 mutations have structural and functional alterations affecting brain circuits that overlap with sites of FOXP2 expression, including regions of the cortex, striatum, and cerebellum. FOXP2 displays complex patterns of expression in the brain, as well as in non-neuronal tissues, suggesting that sophisticated regulatory mechanisms control its spatio-temporal expression. However, to date, little is known about the regulation of FOXP2 or the genomic elements that control its expression. Using chromatin conformation capture (3C), we mapped the human FOXP2 locus to identify putative enhancer regions that engage in long-range interactions with the promoter of this gene. We demonstrate the ability of the identified enhancer regions to drive gene expression. We also show regulation of the FOXP2 promoter and enhancer regions by candidate regulators - FOXP family and TBR1 transcription factors. These data point to regulatory elements that may contribute to the temporal- or tissue-specific expression patterns of human FOXP2. Understanding the upstream regulatory pathways controlling FOXP2 expression will bring new insight into the molecular networks contributing to human language and related disorders.

Mapping the distribution of language related genes FoxP1, FoxP2, and CntnaP2 in the brains of vocal learning bat species.

Genes including FOXP2, FOXP1, and CNTNAP2, have been implicated in human speech and language phenotypes, pointing to a role in the development of normal language-related circuitry in the brain. Although speech and language are unique to humans a comparative approach is possible by addressing language-relevant traits in animal systems. One such trait, vocal learning, represents an essential component of human spoken language, and is shared by cetaceans, pinnipeds, elephants, some birds and bats. Given their vocal learning abilities, gregarious nature, and reliance on vocalizations for social communication and navigation, bats represent an intriguing mammalian system in which to explore language-relevant genes. We used immunohistochemistry to detail the distribution of FoxP2, FoxP1, and Cntnap2 proteins, accompanied by detailed cytoarchitectural histology in the brains of two vocal learning bat species; Phyllostomus discolor and Rousettus aegyptiacus. We show widespread expression of these genes, similar to what has been previously observed in other species, including humans. A striking difference was observed in the adult P. discolor bat, which showed low levels of FoxP2 expression in the cortex that contrasted with patterns found in rodents and nonhuman primates. We created an online, open-access database within which all data can be browsed, searched, and high resolution images viewed to single cell resolution. The data presented herein reveal regions of interest in the bat brain and provide new opportunities to address the role of these language-related genes in complex vocal-motor and vocal learning behaviors in a mammalian model system.

Understanding Neurodevelopmental Disorders: The Promise of Regulatory Variation in the 3'UTRome.

Neurodevelopmental disorders have a strong genetic component, but despite widespread efforts, the specific genetic factors underlying these disorders remain undefined for a large proportion of affected individuals. Given the accessibility of exome sequencing, this problem has thus far been addressed from a protein-centric standpoint; however, protein-coding regions only make up ∼1% to 2% of the human genome. With the advent of whole genome sequencing we are in the midst of a paradigm shift as it is now possible to interrogate the entire sequence of the human genome (coding and noncoding) to fill in the missing heritability of complex disorders. These new technologies bring new challenges, as the number of noncoding variants identified per individual can be overwhelming, making it prudent to focus on noncoding regions of known function, for which the effects of variation can be predicted and directly tested to assess pathogenicity. The 3'UTRome is a region of the noncoding genome that perfectly fulfills these criteria and is of high interest when searching for pathogenic variation related to complex neurodevelopmental disorders. Herein, we review the regulatory roles of the 3'UTRome as binding sites for microRNAs or RNA binding proteins, or during alternative polyadenylation. We detail existing evidence that these regions contribute to neurodevelopmental disorders and outline strategies for identification and validation of novel putatively pathogenic variation in these regions. This evidence suggests that studying the 3'UTRome will lead to the identification of new risk factors, new candidate disease genes, and a better understanding of the molecular mechanisms contributing to neurodevelopmental disorders.

A chromosomal rearrangement in a child with severe speech and language disorder separates FOXP2 from a functional enhancer.

Mutations of FOXP2 in 7q31 cause a rare disorder involving speech apraxia, accompanied by expressive and receptive language impairments. A recent report described a child with speech and language deficits, and a genomic rearrangement affecting chromosomes 7 and 11. One breakpoint mapped to 7q31 and, although outside its coding region, was hypothesised to disrupt FOXP2 expression. We identified an element 2 kb downstream of this breakpoint with epigenetic characteristics of an enhancer. We show that this element drives reporter gene expression in human cell-lines. Thus, displacement of this element by translocation may disturb gene expression, contributing to the observed language phenotype.

FOXP2 drives neuronal differentiation by interacting with retinoic acid signaling pathways.

FOXP2 was the first gene shown to cause a Mendelian form of speech and language disorder. Although developmentally expressed in many organs, loss of a single copy of FOXP2 leads to a phenotype that is largely restricted to orofacial impairment during articulation and linguistic processing deficits. Why perturbed FOXP2 function affects specific aspects of the developing brain remains elusive. We investigated the role of FOXP2 in neuronal differentiation and found that FOXP2 drives molecular changes consistent with neuronal differentiation in a human model system. We identified a network of FOXP2 regulated genes related to retinoic acid signaling and neuronal differentiation. FOXP2 also produced phenotypic changes associated with neuronal differentiation including increased neurite outgrowth and reduced migration. Crucially, cells expressing FOXP2 displayed increased sensitivity to retinoic acid exposure. This suggests a mechanism by which FOXP2 may be able to increase the cellular differentiation response to environmental retinoic acid cues for specific subsets of neurons in the brain. These data demonstrate that FOXP2 promotes neuronal differentiation by interacting with the retinoic acid signaling pathway and regulates key processes required for normal circuit formation such as neuronal migration and neurite outgrowth. In this way, FOXP2, which is found only in specific subpopulations of neurons in the brain, may drive precise neuronal differentiation patterns and/or control localization and connectivity of these FOXP2 positive cells.

A direct molecular link between the autism candidate gene RORa and the schizophrenia candidate MIR137.

Retinoic acid-related orphan receptor alpha gene (RORa) and the microRNA MIR137 have both recently been identified as novel candidate genes for neuropsychiatric disorders. RORa encodes a ligand-dependent orphan nuclear receptor that acts as a transcriptional regulator and miR-137 is a brain enriched small non-coding RNA that interacts with gene transcripts to control protein levels. Given the mounting evidence for RORa in autism spectrum disorders (ASD) and MIR137 in schizophrenia and ASD, we investigated if there was a functional biological relationship between these two genes. Herein, we demonstrate that miR-137 targets the 3'UTR of RORa in a site specific manner. We also provide further support for MIR137 as an autism candidate by showing that a large number of previously implicated autism genes are also putatively targeted by miR-137. This work supports the role of MIR137 as an ASD candidate and demonstrates a direct biological link between these previously unrelated autism candidate genes.

Diagnostic potential of plasmatic MicroRNA signatures in stable and unstable angina.

PURPOSE: We examined circulating miRNA expression profiles in plasma of patients with coronary artery disease (CAD) vs. matched controls, with the aim of identifying novel discriminating biomarkers of Stable (SA) and Unstable (UA) angina. METHODS: An exploratory analysis of plasmatic expression profile of 367 miRNAs was conducted in a group of SA and UA patients and control donors, using TaqMan microRNA Arrays. Screening confirmation and expression analysis were performed by qRT-PCR: all miRNAs found dysregulated were examined in the plasma of troponin-negative UA (n=19) and SA (n=34) patients and control subjects (n=20), matched for sex, age, and cardiovascular risk factors. In addition, the expression of 14 known CAD-associated miRNAs was also investigated. RESULTS: Out of 178 miRNAs consistently detected in plasma samples, 3 showed positive modulation by CAD when compared to controls: miR-337-5p, miR-433, and miR-485-3p. Further, miR-1, -122, -126, -133a, -133b, and miR-199a were positively modulated in both UA and SA patients, while miR-337-5p and miR-145 showed a positive modulation only in SA or UA patients, respectively. ROC curve analyses showed a good diagnostic potential (AUC ≥ 0.85) for miR-1, -126, and -483-5p in SA and for miR-1, -126, and -133a in UA patients vs. controls, respectively. No discriminating AUC values were observed comparing SA vs. UA patients. Hierarchical cluster analysis showed that the combination of miR-1, -133a, and -126 in UA and of miR-1, -126, and -485-3p in SA correctly classified patients vs. controls with an efficiency ≥ 87%. No combination of miRNAs was able to reliably discriminate patients with UA from patients with SA. CONCLUSIONS: This work showed that specific plasmatic miRNA signatures have the potential to accurately discriminate patients with angiographically documented CAD from matched controls. We failed to identify a plasmatic miRNA expression pattern capable to differentiate SA from UA patients.

Histone deacetylase inhibition enhances self renewal and cardioprotection by human cord blood-derived CD34 cells.

BACKGROUND: Use of peripheral blood- or bone marrow-derived progenitors for ischemic heart repair is a feasible option to induce neo-vascularization in ischemic tissues. These cells, named Endothelial Progenitors Cells (EPCs), have been extensively characterized phenotypically and functionally. The clinical efficacy of cardiac repair by EPCs cells remains, however, limited, due to cell autonomous defects as a consequence of risk factors. The devise of "enhancement" strategies has been therefore sought to improve repair ability of these cells and increase the clinical benefit. PRINCIPAL FINDINGS: Pharmacologic inhibition of histone deacetylases (HDACs) is known to enhance hematopoietic stem cells engraftment by improvement of self renewal and inhibition of differentiation in the presence of mitogenic stimuli in vitro. In the present study cord blood-derived CD34(+) were pre-conditioned with the HDAC inhibitor Valproic Acid. This treatment affected stem cell growth and gene expression, and improved ischemic myocardium protection in an immunodeficient mouse model of myocardial infarction. CONCLUSIONS: Our results show that HDAC blockade leads to phenotype changes in CD34(+) cells with enhanced self renewal and cardioprotection.

Endothelial fate and angiogenic properties of human CD34+ progenitor cells in zebrafish.

OBJECTIVE: The vascular competence of human-derived hematopoietic progenitors for postnatal vascularization is still poorly characterized. It is unclear whether, in the absence of ischemia, hematopoietic progenitors participate in neovascularization and whether they play a role in new blood vessel formation by incorporating into developing vessels or by a paracrine action. METHODS AND RESULTS: In the present study, human cord blood-derived CD34(+) (hCD34(+)) cells were transplanted into pre- and postgastrulation zebrafish embryos and in an adult vascular regeneration model induced by caudal fin amputation. When injected before gastrulation, hCD34(+) cells cosegregated with the presumptive zebrafish hemangioblasts, characterized by Scl and Gata2 expression, in the anterior and posterior lateral mesoderm and were involved in early development of the embryonic vasculature. These morphogenetic events occurred without apparent lineage reprogramming, as shown by CD45 expression. When transplanted postgastrulation, hCD34(+) cells were recruited into developing vessels, where they exhibited a potent paracrine proangiogenic action. Finally, hCD34(+) cells rescued vascular defects induced by Vegf-c in vivo targeting and enhanced vascular repair in the zebrafish fin amputation model. CONCLUSIONS: These results indicate an unexpected developmental ability of human-derived hematopoietic progenitors and support the hypothesis of an evolutionary conservation of molecular pathways involved in endothelial progenitor differentiation in vivo.

Human cord blood CD34+ progenitor cells acquire functional cardiac properties through a cell fusion process.

The efficacy of cardiac repair by stem cell administration relies on a successful functional integration of injected cells into the host myocardium. Safety concerns have been raised about the possibility that stem cells may induce foci of arrhythmia in the ischemic myocardium. In a previous work (36), we showed that human cord blood CD34(+) cells, when cocultured on neonatal mouse cardiomyocytes, exhibit excitation-contraction coupling features similar to those of cardiomyocytes, even though no human genes were upregulated. The aims of the present work are to investigate whether human CD34(+) cells, isolated after 1 wk of coculture with neonatal ventricular myocytes, possess molecular and functional properties of cardiomyocytes and to discriminate, using a reporter gene system, whether cardiac differentiation derives from a (trans)differentiation or a cell fusion process. Umbilical cord blood CD34(+) cells were isolated by a magnetic cell sorting method, transduced with a lentiviral vector carrying the enhanced green fluorescent protein (EGFP) gene, and seeded onto primary cultures of spontaneously beating rat neonatal cardiomyocytes. Cocultured EGFP(+)/CD34(+)-derived cells were analyzed for their electrophysiological features at different time points. After 1 wk in coculture, EGFP(+) cells, in contact with cardiomyocytes, were spontaneously contracting and had a maximum diastolic potential (MDP) of -53.1 mV, while those that remained isolated from the surrounding myocytes did not contract and had a depolarized resting potential of -11.4 mV. Cells were then resuspended and cultured at low density to identify EGFP(+) progenitor cell derivatives. Under these conditions, we observed single EGFP(+) beating cells that had acquired an hyperpolarization-activated current typical of neonatal cardiomyocytes (EGFP(+) cells, -2.24 ± 0.89 pA/pF; myocytes, -1.99 ± 0.63 pA/pF, at -125 mV). To discriminate between cell autonomous differentiation and fusion, EGFP(+)/CD34(+) cells were cocultured with cardiac myocytes infected with a red fluorescence protein-lentiviral vector; under these conditions we found that 100% of EGFP(+) cells were also red fluorescent protein positive, suggesting cell fusion as the mechanism by which cardiac functional features are acquired.

Circulating microRNAs are new and sensitive biomarkers of myocardial infarction.

AIMS: Circulating microRNAs (miRNAs) may represent a novel class of biomarkers; therefore, we examined whether acute myocardial infarction (MI) modulates miRNAs plasma levels in humans and mice. METHODS AND RESULTS: Healthy donors (n = 17) and patients (n = 33) with acute ST-segment elevation MI (STEMI) were evaluated. In one cohort (n = 25), the first plasma sample was obtained 517 ± 309 min after the onset of MI symptoms and after coronary reperfusion with percutaneous coronary intervention (PCI); miR-1, -133a, -133b, and -499-5p were ~15- to 140-fold control, whereas miR-122 and -375 were ~87-90% lower than control; 5 days later, miR-1, -133a, -133b, -499-5p, and -375 were back to baseline, whereas miR-122 remained lower than control through Day 30. In additional patients (n = 8; four treated with thrombolysis and four with PCI), miRNAs and troponin I (TnI) were quantified simultaneously starting 156 ± 72 min after the onset of symptoms and at different times thereafter. Peak miR-1, -133a, and -133b expression and TnI level occurred at a similar time, whereas miR-499-5p exhibited a slower time course. In mice, miRNAs plasma levels and TnI were measured 15 min after coronary ligation and at different times thereafter. The behaviour of miR-1, -133a, -133b, and -499-5p was similar to STEMI patients; further, reciprocal changes in the expression levels of these miRNAs were found in cardiac tissue 3-6 h after coronary ligation. In contrast, miR-122 and -375 exhibited minor changes and no significant modulation. In mice with acute hind-limb ischaemia, there was no increase in the plasma level of the above miRNAs. CONCLUSION: Acute MI up-regulated miR-1, -133a, -133b, and -499-5p plasma levels, both in humans and mice, whereas miR-122 and -375 were lower than control only in STEMI patients. These miRNAs represent novel biomarkers of cardiac damage.