Lizards are evolutionarily the most closely related vertebrates to humans that can lose and regrow an entire appendage. Regeneration in lizards involves differential expression of hundreds of genes that regulate wound healing, musculoskeletal development, hormonal response, and embryonic morphogenesis. While microRNAs are able to regulate large groups of genes, their role in lizard regeneration has not been investigated.

The molecular events leading to the development of the bat wing remain largely unknown, and are thought to be caused, in part, by changes in gene expression during limb development. These expression changes could be instigated by variations in gene regulatory enhancers. Here, we used a comparative genomics approach to identify regions that evolved rapidly in the bat ancestor, but are highly conserved in other vertebrates. We discovered 166 bat accelerated regions (BARs) that overlap H3K27ac and p300 ChIP-seq peaks in developing mouse limbs. Using a mouse enhancer assay, we show that five Myotis lucifugus BARs drive gene expression in the developing mouse limb, with the majority showing differential enhancer activity compared to the mouse orthologous BAR sequences. These include BAR116, which is located telomeric to the HoxD cluster and had robust forelimb expression for the M. lucifugus sequence and no activity for the mouse sequence at embryonic day 12.5. Developing limb expression analysis of Hoxd10-Hoxd13 in Miniopterus natalensis bats showed a high-forelimb weak-hindlimb expression for Hoxd10-Hoxd11, similar to the expression trend observed for M. lucifugus BAR116 in mice, suggesting that it could be involved in the regulation of the bat HoxD complex. Combined, our results highlight novel regulatory regions that could be instrumental for the morphological differences leading to the development of the bat wing.

Although inhaled corticosteroid (ICS) medication is considered the cornerstone treatment for patients with persistent asthma, few ICS pharmacogenomic studies have involved nonwhite populations.

Group 2 innate lymphoid cells (ILC2s) are distributed systemically and produce type 2 cytokines in response to a variety of stimuli, including the epithelial cytokines interleukin (IL)-25, IL-33, and thymic stromal lymphopoietin (TSLP). Transcriptional profiling of ILC2s from different tissues, however, grouped ILC2s according to their tissue of origin, even in the setting of combined IL-25-, IL-33-receptor-, and TSLP-receptor-deficiency. Single-cell profiling confirmed a tissue-organizing transcriptome and identified ILC2 subsets expressing distinct activating receptors, including the major subset of skin ILC2s, which were activated preferentially by IL-18. Tissue ILC2 subsets were unaltered in number and expression in germ-free mice, suggesting that endogenous, tissue-derived signals drive the maturation of ILC2 subsets by controlling expression of distinct patterns of activating receptors, thus anticipating tissue-specific perturbations occurring later in life.

Metformin is used as a first-line therapy for type 2 diabetes (T2D) and prescribed for numerous other diseases. However, its mechanism of action in the liver has yet to be characterized in a systematic manner. To comprehensively identify genes and regulatory elements associated with metformin treatment, we carried out RNA-seq and ChIP-seq (H3K27ac, H3K27me3) on primary human hepatocytes from the same donor treated with vehicle control, metformin or metformin and compound C, an AMP-activated protein kinase (AMPK) inhibitor (allowing to identify AMPK-independent pathways). We identified thousands of metformin responsive AMPK-dependent and AMPK-independent differentially expressed genes and regulatory elements. We functionally validated several elements for metformin-induced promoter and enhancer activity. These include an enhancer in an ataxia telangiectasia mutated (ATM) intron that has SNPs in linkage disequilibrium with a metformin treatment response GWAS lead SNP (rs11212617) that showed increased enhancer activity for the associated haplotype. Expression quantitative trait locus (eQTL) liver analysis and CRISPR activation suggest that this enhancer could be regulating ATM, which has a known role in AMPK activation, and potentially also EXPH5 and DDX10, its neighboring genes. Using ChIP-seq and siRNA knockdown, we further show that activating transcription factor 3 (ATF3), our top metformin upregulated AMPK-dependent gene, could have an important role in gluconeogenesis repression. Our findings provide a genome-wide representation of metformin hepatic response, highlight important sequences that could be associated with interindividual variability in glycemic response to metformin and identify novel T2D treatment candidates.

Standard cell culture guidelines often use media supplemented with antibiotics to prevent cell contamination. However, relatively little is known about the effect of antibiotic use in cell culture on gene expression and the extent to which this treatment could confound results. To comprehensively characterize the effect of antibiotic treatment on gene expression, we performed RNA-seq and ChIP-seq for H3K27ac on HepG2 cells, a human liver cell line commonly used for pharmacokinetic, metabolism and genomic studies, cultured in media supplemented with penicillin-streptomycin (PenStrep) or vehicle control. We identified 209 PenStrep-responsive genes, including transcription factors such as ATF3 that are likely to alter the regulation of other genes. Pathway analyses found a significant enrichment for "xenobiotic metabolism signaling" and "PXR/RXR activation" pathways. Our H3K27ac ChIP-seq identified 9,514 peaks that are PenStrep responsive. These peaks were enriched near genes that function in cell differentiation, tRNA modification, nuclease activity and protein dephosphorylation. Our results suggest that PenStrep treatment can significantly alter gene expression and regulation in a common liver cell type such as HepG2, advocating that antibiotic treatment should be taken into account when carrying out genetic, genomic or other biological assays in cultured cells.

Trauma to the spinal cord and brain can result in irreparable loss of function. This failure of recovery is in part due to inhibition of axon regeneration by myelin and chondroitin sulfate proteoglycans (CSPGs). Peripheral nervous system (PNS) neurons exhibit increased regenerative ability compared to central nervous system neurons, even in the presence of inhibitory environments. Previously, we identified over a thousand genes differentially expressed in PNS neurons relative to CNS neurons. These genes represent intrinsic differences that may account for the PNS's enhanced regenerative ability. Cerebellar neurons were transfected with cDNAs for each of these PNS genes to assess their ability to enhance neurite growth on inhibitory (CSPG) or permissive (laminin) substrates. Using high content analysis, we evaluated the phenotypic profile of each neuron to extract meaningful data for over 1100 genes. Several known growth associated proteins potentiated neurite growth on laminin. Most interestingly, novel genes were identified that promoted neurite growth on CSPGs (GPX3, EIF2B5, RBMX). Bioinformatic approaches also uncovered a number of novel gene families that altered neurite growth of CNS neurons.

Injury to the central nervous system (CNS) can result in lifelong loss of function due in part to the regenerative failure of CNS neurons. Inhibitory proteins derived from myelin and the astroglial scar are major barriers for the successful regeneration of injured CNS neurons. Previously, we described the identification of a novel compound, F05, which promotes neurite growth from neurons challenged with inhibitory substrates in vitro, and promotes axonal regeneration in vivo (Usher et al., 2010). To identify additional regeneration-promoting compounds, we used F05-induced gene expression profiles to query the Broad Institute Connectivity Map, a gene expression database of cells treated with >1300 compounds. Despite no shared chemical similarity, F05-induced changes in gene expression were remarkably similar to those seen with a group of piperazine phenothiazine antipsychotics (PhAPs). In contrast to antipsychotics of other structural classes, PhAPs promoted neurite growth of CNS neurons challenged with two different glial derived inhibitory substrates. Our pharmacological studies suggest a mechanism whereby PhAPs promote growth through antagonism of calmodulin signaling, independent of dopamine receptor antagonism. These findings shed light on mechanisms underlying neurite-inhibitory signaling, and suggest that clinically approved antipsychotic compounds may be repurposed for use in CNS injured patients.

The genetic control of gene expression is a core component of human physiology. For the past several years, transcriptome-wide association studies have leveraged large datasets of linked genotype and RNA sequencing information to create a powerful gene-based test of association that has been used in dozens of studies. While numerous discoveries have been made, the populations in the training data are overwhelmingly of European descent, and little is known about the generalizability of these models to other populations. Here, we test for cross-population generalizability of gene expression prediction models using a dataset of African American individuals with RNA-Seq data in whole blood. We find that the default models trained in large datasets such as GTEx and DGN fare poorly in African Americans, with a notable reduction in prediction accuracy when compared to European Americans. We replicate these limitations in cross-population generalizability using the five populations in the GEUVADIS dataset. Via realistic simulations of both populations and gene expression, we show that accurate cross-population generalizability of transcriptome prediction only arises when eQTL architecture is substantially shared across populations. In contrast, models with non-identical eQTLs showed patterns similar to real-world data. Therefore, generating RNA-Seq data in diverse populations is a critical step towards multi-ethnic utility of gene expression prediction.

Monitor lizards are unique among ectothermic reptiles in that they have high aerobic capacity and distinctive cardiovascular physiology resembling that of endothermic mammals. Here, we sequence the genome of the Komodo dragon Varanus komodoensis, the largest extant monitor lizard, and generate a high-resolution de novo chromosome-assigned genome assembly for V. komodoensis using a hybrid approach of long-range sequencing and single-molecule optical mapping. Comparing the genome of V. komodoensis with those of related species, we find evidence of positive selection in pathways related to energy metabolism, cardiovascular homoeostasis, and haemostasis. We also show species-specific expansions of a chemoreceptor gene family related to pheromone and kairomone sensing in V. komodoensis and other lizard lineages. Together, these evolutionary signatures of adaptation reveal the genetic underpinnings of the unique Komodo dragon sensory and cardiovascular systems, and suggest that selective pressure altered haemostasis genes to help Komodo dragons evade the anticoagulant effects of their own saliva. The Komodo dragon genome is an important resource for understanding the biology of monitor lizards and reptiles worldwide.

Cancer cells utilize multiple nutrient scavenging mechanisms to support growth and survival in nutrient-poor, hypoxic tumor microenvironments. Among these mechanisms, macropinocytosis has emerged as an important pathway of extracellular nutrient acquisition in cancer cells, particularly in tumors with activated RAS signaling, such as pancreatic cancer. However, the absence of a clinically available inhibitor, as well as the gap of knowledge in macropinocytosis regulation, remain a hurdle for its use for cancer therapy. Here, we use the Informer set library to identify novel regulators of macropinocytosis-dependent growth in pancreatic cancer cells. Understanding how these regulators function will allow us to provide novel opportunities for therapeutic intervention.

Asthma is associated with chronic changes in the airway epithelium, a key target of SARS-CoV-2. Many epithelial changes are driven by the type 2 cytokine IL-13, but the effects of IL-13 on SARS-CoV-2 infection are unknown.

Inguinal hernia repair is one of the most commonly performed operations in the world, yet little is known about the genetic mechanisms that predispose individuals to develop inguinal hernias. We perform a genome-wide association analysis of surgically confirmed inguinal hernias in 72,805 subjects (5,295 cases and 67,510 controls) and confirm top associations in an independent cohort of 92,444 subjects with self-reported hernia repair surgeries (9,701 cases and 82,743 controls). We identify four novel inguinal hernia susceptibility loci in the regions of EFEMP1, WT1, EBF2 and ADAMTS6. Moreover, we observe expression of all four genes in mouse connective tissue and network analyses show an important role for two of these genes (EFEMP1 and WT1) in connective tissue maintenance/homoeostasis. Our findings provide insight into the aetiology of hernia development and highlight genetic pathways for studies of hernia development and its treatment.

Lizards, which are amniote vertebrates like humans, are able to lose and regenerate a functional tail. Understanding the molecular basis of this process would advance regenerative approaches in amniotes, including humans. We have carried out the first transcriptomic analysis of tail regeneration in a lizard, the green anole Anolis carolinensis, which revealed 326 differentially expressed genes activating multiple developmental and repair mechanisms. Specifically, genes involved in wound response, hormonal regulation, musculoskeletal development, and the Wnt and MAPK/FGF pathways were differentially expressed along the regenerating tail axis. Furthermore, we identified 2 microRNA precursor families, 22 unclassified non-coding RNAs, and 3 novel protein-coding genes significantly enriched in the regenerating tail. However, high levels of progenitor/stem cell markers were not observed in any region of the regenerating tail. Furthermore, we observed multiple tissue-type specific clusters of proliferating cells along the regenerating tail, not localized to the tail tip. These findings predict a different mechanism of regeneration in the lizard than the blastema model described in the salamander and the zebrafish, which are anamniote vertebrates. Thus, lizard tail regrowth involves the activation of conserved developmental and wound response pathways, which are potential targets for regenerative medical therapies.

In addition to their protein coding function, exons can also serve as transcriptional enhancers. Mutations in these exonic-enhancers (eExons) could alter both protein function and transcription. However, the functional consequence of eExon mutations is not well known. Here, using massively parallel reporter assays, we dissect the enhancer activity of three liver eExons (SORL1 exon 17, TRAF3IP2 exon 2, PPARG exon 6) at single nucleotide resolution in the mouse liver. We find that both synonymous and non-synonymous mutations have similar effects on enhancer activity and many of the deleterious mutation clusters overlap known liver-associated transcription factor binding sites. Carrying a similar massively parallel reporter assay in HeLa cells with these three eExons found differences in their mutation profiles compared to the liver, suggesting that enhancers could have distinct operating profiles in different tissues. Our results demonstrate that eExon mutations could lead to multiple phenotypes by disrupting both the protein sequence and enhancer activity and that enhancers can have distinct mutation profiles in different cell types.

Neurons in the peripheral nervous system (PNS) display a higher capacity to regenerate after injury than those in the central nervous system, suggesting cell specific transcriptional modules underlying axon growth and inhibition. We report a systems biology based search for PNS specific transcription factors (TFs). Messenger RNAs enriched in dorsal root ganglion (DRG) neurons compared to cerebellar granule neurons (CGNs) were identified using subtractive hybridization and DNA microarray approaches. Network and transcription factor binding site enrichment analyses were used to further identify TFs that may be differentially active. Combining these techniques, we identified 32 TFs likely to be enriched and/or active in the PNS. Twenty-five of these TFs were then tested for an ability to promote CNS neurite outgrowth in an overexpression screen. Real-time PCR and immunohistochemical studies confirmed that one representative TF, STAT3, is intrinsic to PNS neurons, and that constitutively active STAT3 is sufficient to promote CGN neurite outgrowth.

T helper 17 cells (Th17) are critical for fighting infections at mucosal surfaces; however, they have also been found to contribute to the pathogenesis of multiple autoimmune diseases and have been targeted therapeutically. Due to the role of Th17 cells in autoimmune pathogenesis, it is important to understand the factors that control Th17 development. Here we identify the activating transcription factor 7 interacting protein (ATF7ip) as a critical regulator of Th17 differentiation. Mice with T cell-specific deletion of Atf7ip have impaired Th17 differentiation secondary to the aberrant overproduction of IL-2 with T cell receptor (TCR) stimulation and are resistant to colitis in vivo. ChIP-seq studies identified ATF7ip as an inhibitor of Il2 gene expression through the deposition of the repressive histone mark H3K9me3 in the Il2-Il21 intergenic region. These results demonstrate a new epigenetic pathway by which IL-2 production is constrained, and this may open up new avenues for modulating its production.

Synthetic lethal screens have the potential to identify new vulnerabilities incurred by specific cancer mutations but have been hindered by lack of agreement between studies. In the case of KRAS, we identify that published synthetic lethal screen hits significantly overlap at the pathway rather than gene level. Analysis of pathways encoded as protein networks could identify synthetic lethal candidates that are more reproducible than those previously reported. Lack of overlap likely stems from biological rather than technical limitations as most synthetic lethal phenotypes are strongly modulated by changes in cellular conditions or genetic context, the latter determined using a pairwise genetic interaction map that identifies numerous interactions that suppress synthetic lethal effects. Accounting for pathway, cellular and genetic context nominates a DNA repair dependency in KRAS-mutant cells, mediated by a network containing BRCA1. We provide evidence for why most reported synthetic lethals are not reproducible which is addressable using a multi-faceted testing framework.

Adolescent idiopathic scoliosis (AIS), a sideways curvature of the spine, is the most common pediatric musculoskeletal disorder, affecting ~3% of the population worldwide. However, its genetic bases and tissues of origin remain largely unknown. Several genome-wide association studies (GWAS) have implicated nucleotide variants in non-coding sequences that control genes with important roles in cartilage, muscle, bone, connective tissue and intervertebral disks (IVDs) as drivers of AIS susceptibility. Here, we set out to define the expression of AIS-associated genes and active regulatory elements by performing RNA-seq and chromatin immunoprecipitation-sequencing against H3 lysine 27 acetylation in these tissues in mouse and human. Our study highlights genetic pathways involving AIS-associated loci that regulate chondrogenesis, IVD development and connective tissue maintenance and homeostasis. In addition, we identify thousands of putative AIS-associated regulatory elements which may orchestrate tissue-specific expression in musculoskeletal tissues of the spine. Quantification of enhancer activity of several candidate regulatory elements from our study identifies three functional enhancers carrying AIS-associated GWAS SNPs at the ADGRG6 and BNC2 loci. Our findings provide a novel genome-wide catalog of AIS-relevant genes and regulatory elements and aid in the identification of novel targets for AIS causality and treatment.

Idiopathic scoliosis (IS) is a common paediatric musculoskeletal disease that displays a strong female bias. By performing a genome-wide association study (GWAS) of 3,102 individuals, we identify significant associations with 20p11.22 SNPs for females (P=6.89 × 10(-9)) but not males (P=0.71). This association with IS is also found in independent female cohorts from the United States of America and Japan (overall P=2.15 × 10(-10), OR=1.30 (rs6137473)). Unexpectedly, the 20p11.22 IS risk alleles were previously associated with protection from early-onset alopecia, another sexually dimorphic condition. The 174-kb associated locus is distal to PAX1, which encodes paired box 1, a transcription factor involved in spine development. We identify a sequence in the associated locus with enhancer activity in zebrafish somitic muscle and spinal cord, an activity that is abolished by IS-associated SNPs. We thus identify a sexually dimorphic IS susceptibility locus, and propose the first functionally defined candidate mutations in an enhancer that may regulate expression in specific spinal cells.