The cardiovascular risk factor LPA has a puzzling distribution among mammals, its presence being limited to a subset of primates and a member of the insectivore lineage, the hedgehog. To explore the evolutionary history of LPA, we performed extensive genomic sequence comparisons of multiple species with and without an LPA gene product, such as human, baboon, hedgehog, lemur, and mouse. This analysis indicated that LPA arose independently in a subset of primates, including baboon and human, and an insectivore, the hedgehog, and was not simply lost by species lacking it. The similar structural domains shared by the hedgehog and primate LPA indicate that they were formed by a unique molecular mechanism involving the convergent evolution of paralogous genes in these distant species.

Members of the apolipoprotein gene cluster (APOA1/C3/A4/A5) on human chromosome 11q23 play an important role in lipid metabolism. Polymorphisms in both APOA5 and APOC3 are strongly associated with plasma triglyceride concentrations. The close genomic locations of these two genes as well as their functional similarity have hindered efforts to define whether each gene independently influences human triglyceride concentrations. In this study, we examined the linkage disequilibrium and haplotype structure of 49 SNPs in a 150-kb region spanning the gene cluster. We identified a total of five common APOA5 haplotypes with a frequency of greater than 8% in samples of northern European origin. The APOA5 haplotype block did not extend past the 7 SNPs in the gene and was separated from the other apolipoprotein gene in the cluster by a region of significantly increased recombination. Furthermore, one previously identified triglyceride risk haplotype of APOA5 (APOA5*3) showed no association with three APOC3 SNPs previously associated with triglyceride concentrations, in contrast to the other risk haplotype (APOA5*2), which was associated with all three minor APOC3 SNP alleles. These results highlight the complex genetic relationship between APOA5 and APOC3 and support the notion that APOA5 represents an independent risk gene affecting plasma triglyceride concentrations in humans.

Massively parallel sequencing of barcoded DNA samples significantly increases screening efficiency for clinically important genes. Short read aligners are well suited to single nucleotide and indel detection. However, methods for CNV detection from targeted enrichment are lacking. We present a method combining coverage with map information for the identification of deletions and duplications in targeted sequence data.

Little is known about genetic mechanisms that regulate the ratio of cortical excitatory and inhibitory neurons. We show that NPAS1 and NPAS3 transcription factors (TFs) are expressed in progenitor domains of the mouse basal ganglia (subpallium, MGE, and CGE). NPAS1(-/-) mutants had increased proliferation, ERK signaling, and expression of Arx in the MGE and CGE. NPAS1(-/-) mutants also had increased neocortical inhibition (sIPSC and mIPSC) and generated an excess of somatostatin(+) (SST) (MGE-derived) and vasoactive intestinal polypeptide(+) (VIP) (CGE-derived) neocortical interneurons, but had a normal density of parvalbumin(+) (PV) (MGE-derived) interneurons. In contrast, NPAS3(-/-) mutants showed decreased proliferation and ERK signaling in progenitors of the ganglionic eminences and had fewer SST(+) and VIP(+) interneurons. NPAS1 repressed activity of an Arx enhancer, and Arx overexpression resulted in increased proliferation of CGE progenitors. These results provide insights into genetic regulation of cortical interneuron numbers and cortical inhibitory tone.

Harnessing the biotechnological potential of the large number of proteins available in sequence databases requires scalable methods for functional characterization. Here we propose a workflow to address this challenge by combining phylogenomic guided DNA synthesis with high-throughput mass spectrometry and apply it to the systematic characterization of GH1 β-glucosidases, a family of enzymes necessary for biomass hydrolysis, an important step in the conversion of lignocellulosic feedstocks to fuels and chemicals. We synthesized and expressed 175 GH1s, selected from over 2000 candidate sequences to cover maximum sequence diversity. These enzymes were functionally characterized over a range of temperatures and pHs using nanostructure-initiator mass spectrometry (NIMS), generating over 10,000 data points. When combined with HPLC-based sugar profiling, we observed GH1 enzymes active over a broad temperature range and toward many different β-linked disaccharides. For some GH1s we also observed activity toward laminarin, a more complex oligosaccharide present as a major component of macroalgae. An area of particular interest was the identification of GH1 enzymes compatible with the ionic liquid 1-ethyl-3-methylimidazolium acetate ([C2mim][OAc]), a next-generation biomass pretreatment technology. We thus searched for GH1 enzymes active at 70 °C and 20% (v/v) [C2mim][OAc] over the course of a 24-h saccharification reaction. Using our unbiased approach, we identified multiple enzymes of different phylogentic origin with such activities. Our approach of characterizing sequence diversity through targeted gene synthesis coupled to high-throughput screening technologies is a broadly applicable paradigm for a wide range of biological problems.

Short non-coding transcripts can be transcribed from distant-acting transcriptional enhancer loci, but the prevalence of such enhancer RNAs (eRNAs) within the transcriptome, and the association of eRNA expression with tissue-specific enhancer activity in vivo remain poorly understood. Here, we investigated the expression dynamics of tissue-specific non-coding RNAs in embryonic mouse tissues via deep RNA sequencing. Overall, approximately 80% of validated in vivo enhancers show tissue-specific RNA expression that correlates with tissue-specific enhancer activity. Globally, we identified thousands of tissue-specifically transcribed non-coding regions (TSTRs) displaying various genomic hallmarks of bona fide enhancers. In transgenic mouse reporter assays, over half of tested TSTRs functioned as enhancers with reproducible activity in the predicted tissue. Together, our results demonstrate that tissue-specific eRNA expression is a common feature of in vivo enhancers, as well as a major source of extragenic transcription, and that eRNA expression signatures can be used to predict tissue-specific enhancers independent of known epigenomic enhancer marks.

The accurate and comprehensive identification of functional regulatory sequences in mammalian genomes remains a major challenge. Here we describe site-specific integration fluorescence-activated cell sorting followed by sequencing (SIF-seq), an unbiased, medium-throughput functional assay for the discovery of distant-acting enhancers. Targeted single-copy genomic integration into pluripotent cells, reporter assays and flow cytometry are coupled with high-throughput DNA sequencing to enable parallel screening of large numbers of DNA sequences. By functionally interrogating >500 kilobases (kb) of mouse and human sequence in mouse embryonic stem cells for enhancer activity we identified enhancers at pluripotency loci including NANOG. In in vitro-differentiated cardiomyocytes and neural progenitor cells, we identified cardiac enhancers and neuronal enhancers, respectively. SIF-seq is a powerful and flexible method for de novo functional identification of mammalian enhancers in a potentially wide variety of cell types.

The embryonic basal ganglia generates multiple projection neurons and interneuron subtypes from distinct progenitor domains. Combinatorial interactions of transcription factors and chromatin are thought to regulate gene expression. In the medial ganglionic eminence, the NKX2-1 transcription factor controls regional identity and, with LHX6, is necessary to specify pallidal projection neurons and forebrain interneurons. Here, we dissected the molecular functions of NKX2-1 by defining its chromosomal binding, regulation of gene expression, and epigenetic state. NKX2-1 binding at distal regulatory elements led to a repressed epigenetic state and transcriptional repression in the ventricular zone. Conversely, NKX2-1 is required to establish a permissive chromatin state and transcriptional activation in the sub-ventricular and mantle zones. Moreover, combinatorial binding of NKX2-1 and LHX6 promotes transcriptionally permissive chromatin and activates genes expressed in cortical migrating interneurons. Our integrated approach provides a foundation for elucidating transcriptional networks guiding the development of the MGE and its descendants.

Enteric fermentation by farmed ruminant animals is a major source of methane and constitutes the second largest anthropogenic contributor to global warming. Reducing methane emissions from ruminants is needed to ensure sustainable animal production in the future. Methane yield varies naturally in sheep and is a heritable trait that can be used to select animals that yield less methane per unit of feed eaten. We previously demonstrated elevated expression of hydrogenotrophic methanogenesis pathway genes of methanogenic archaea in the rumens of high methane yield (HMY) sheep compared to their low methane yield (LMY) counterparts. Methane production in the rumen is strongly connected to microbial hydrogen production through fermentation processes. In this study, we investigate the contribution that rumen bacteria make to methane yield phenotypes in sheep.

Analysis of the increasing wealth of metagenomic data collected from diverse environments can lead to the discovery of novel branches on the tree of life. Here we analyse 5.2 Tb of metagenomic data collected globally to discover a novel bacterial phylum ('Candidatus Kryptonia') found exclusively in high-temperature pH-neutral geothermal springs. This lineage had remained hidden as a taxonomic 'blind spot' because of mismatches in the primers commonly used for ribosomal gene surveys. Genome reconstruction from metagenomic data combined with single-cell genomics results in several high-quality genomes representing four genera from the new phylum. Metabolic reconstruction indicates a heterotrophic lifestyle with conspicuous nutritional deficiencies, suggesting the need for metabolic complementarity with other microbes. Co-occurrence patterns identifies a number of putative partners, including an uncultured Armatimonadetes lineage. The discovery of Kryptonia within previously studied geothermal springs underscores the importance of globally sampled metagenomic data in detection of microbial novelty, and highlights the extraordinary diversity of microbial life still awaiting discovery.

The mammalian telencephalon plays critical roles in cognition, motor function, and emotion. Though many of the genes required for its development have been identified, the distant-acting regulatory sequences orchestrating their in vivo expression are mostly unknown. Here, we describe a digital atlas of in vivo enhancers active in subregions of the developing telencephalon. We identified more than 4,600 candidate embryonic forebrain enhancers and studied the in vivo activity of 329 of these sequences in transgenic mouse embryos. We generated serial sets of histological brain sections for 145 reproducible forebrain enhancers, resulting in a publicly accessible web-based data collection comprising more than 32,000 sections. We also used epigenomic analysis of human and mouse cortex tissue to directly compare the genome-wide enhancer architecture in these species. These data provide a primary resource for investigating gene regulatory mechanisms of telencephalon development and enable studies of the role of distant-acting enhancers in neurodevelopmental disorders.

Development and function of the human heart depend on the dynamic control of tissue-specific gene expression by distant-acting transcriptional enhancers. To generate an accurate genome-wide map of human heart enhancers, we used an epigenomic enhancer discovery approach and identified ∼6,200 candidate enhancer sequences directly from fetal and adult human heart tissue. Consistent with their predicted function, these elements were markedly enriched near genes implicated in heart development, function and disease. To further validate their in vivo enhancer activity, we tested 65 of these human sequences in a transgenic mouse enhancer assay and observed that 43 (66%) drove reproducible reporter gene expression in the heart. These results support the discovery of a genome-wide set of noncoding sequences highly enriched in human heart enhancers that is likely to facilitate downstream studies of the role of enhancers in development and pathological conditions of the heart.

Several methods to identify tagging single-nucleotide polymorphisms (SNPs) are in common use for genetic epidemiologic studies; however, there may be loss of information when using only a subset of SNPs. We sought to compare the ability of commonly used pairwise, multimarker, and haplotype-based tagging SNP selection methods to detect known associations with quantitative expression phenotypes. Using data from HapMap release 21 on unrelated Utah residents with ancestors from northern and western Europe (CEPH-Utah, CEU), we selected tagging SNPs in five chromosomal regions using ldSelect, Tagger, and TagSNPs. We found that SNP subsets did not substantially overlap, and that the use of trio data did not greatly impact SNP selection. We then tested associations between HapMap genotypes and expression phenotypes on 28 CEU individuals as part of Genetic Analysis Workshop 15. Relative to the use of all SNPs (n = 210 SNPs across all regions), most subset methods were able to detect single-SNP and haplotype associations. Generally, pairwise selection approaches worked extremely well, relative to use of all SNPs, with marked reductions in the number of SNPs required. Haplotype-based approaches, which had identified smaller SNP subsets, missed associations in some regions. We conclude that the optimal tagging SNP method depends on the true model of the genetic association (i.e., whether a SNP or haplotype is responsible); unfortunately, this is often unknown at the time of SNP selection. Additional evaluations using empirical and simulated data are needed.

Large-scale genome sequencing is poised to provide a substantial increase in the rate of discovery of disease-associated mutations, but the functional interpretation of such mutations remains challenging. Here we show that deletions of a sequence on human chromosome 16 that we term the intestine-critical region (ICR) cause intractable congenital diarrhoea in infants1,2. Reporter assays in transgenic mice show that the ICR contains a regulatory sequence that activates transcription during the development of the gastrointestinal system. Targeted deletion of the ICR in mice caused symptoms that recapitulated the human condition. Transcriptome analysis revealed that an unannotated open reading frame (Percc1) flanks the regulatory sequence, and the expression of this gene was lost in the developing gut of mice that lacked the ICR. Percc1-knockout mice displayed phenotypes similar to those observed upon ICR deletion in mice and patients, whereas an ICR-driven Percc1 transgene was sufficient to rescue the phenotypes found in mice that lacked the ICR. Together, our results identify a gene that is critical for intestinal function and underscore the need for targeted in vivo studies to interpret the growing number of clinical genetic findings that do not affect known protein-coding genes.

Medial ganglionic eminence (MGE)-derived somatostatin (SST)+ and parvalbumin (PV)+ cortical interneurons (CINs), have characteristic molecular, anatomical and physiological properties. However, mechanisms regulating their diversity remain poorly understood. Here, we show that conditional loss of the Tuberous Sclerosis Complex (TSC) gene, Tsc1, which inhibits the mammalian target of rapamycin (MTOR), causes a subset of SST+ CINs, to express PV and adopt fast-spiking (FS) properties, characteristic of PV+ CINs. Milder intermediate phenotypes also occur when only one allele of Tsc1 is deleted. Notably, treatment of adult mice with rapamycin, which inhibits MTOR, reverses the phenotypes. These data reveal novel functions of MTOR signaling in regulating PV expression and FS properties, which may contribute to TSC neuropsychiatric symptoms. Moreover, they suggest that CINs can exhibit properties intermediate between those classically associated with PV+ or SST+ CINs, which may be dynamically regulated by the MTOR signaling.

In utero exposure to maternal immune activation (MIA) is an environmental risk factor for neurodevelopmental and neuropsychiatric disorders. Animal models provide an opportunity to identify mechanisms driving neuropathology associated with MIA. We performed time-course transcriptional profiling of mouse cortical development following induced MIA via poly(I:C) injection at E12.5. MIA-driven transcriptional changes were validated via protein analysis, and parallel perturbations to cortical neuroanatomy were identified via imaging. MIA-induced acute upregulation of genes associated with hypoxia, immune signaling, and angiogenesis, by 6 hr following exposure. This acute response was followed by changes in proliferation, neuronal and glial specification, and cortical lamination that emerged at E14.5 and peaked at E17.5. Decreased numbers of proliferative cells in germinal zones and alterations in neuronal and glial populations were identified in the MIA-exposed cortex. Overall, paired transcriptomic and neuroanatomical characterization revealed a sequence of perturbations to corticogenesis driven by mid-gestational MIA.

The TSC1 and TSC2 genes are connected to multiple syndromes from Tuberous Sclerosis Complex (TSC) to autism spectrum disorder (ASD), with uncertainty if genetic variants cause all or subsets of phenotypes based on the location and type of change. For TSC1, few have addressed if non-TSC associated genetic variants have direct contributions to changes in neurological genotype-to-phenotype impacts, including elevated rates of ASD and seizures. Dominant variants cause TSC, yet TSC1 has many heritable variants not dominant for TSC that are poorly understood in neurological function, with some associated with ASD. Herein, we examined how missense variants in TSC1, R336W, T360N, T393I, S403L, and H732Y, impacted the development of cortical inhibitory interneurons, cell-types whose molecular, cellular, and physiological properties are altered after the loss of mouse TSC1. We found these variants complemented a known phenotype caused by loss of TSC1, increased cell size. However, distinct variants, particularly S403L showed deficits in complementing an increase in parvalbumin levels and exhibited smaller amplitude after hyperpolarizations. Overall, these data show that subtle phenotypes can be induced by some TSC1 missense variants and provide an in vivo system to assess TSC1 variants' neurological impact better.

Genetic studies of autism have revealed causal roles for chromatin remodeling gene mutations. Chromodomain helicase DNA binding protein 8 (CHD8) encodes a chromatin remodeler with significant de novo mutation rates in sporadic autism. However, relationships between CHD8 genomic function and autism-relevant biology remain poorly elucidated. Published studies utilizing ChIP-seq to map CHD8 protein-DNA interactions have high variability, consistent with technical challenges and limitations associated with this method. Thus, complementary approaches are needed to establish CHD8 genomic targets and regulatory functions in developing brain. We used in utero CHD8 Targeted DamID followed by sequencing (TaDa-seq) to characterize CHD8 binding in embryonic mouse cortex. CHD8 TaDa-seq reproduced interaction patterns observed from ChIP-seq and further highlighted CHD8 distal interactions associated with neuronal loci. This study establishes TaDa-seq as a useful alternative for mapping protein-DNA interactions in vivo and provides insights into the regulatory targets of CHD8 and autism-relevant pathophysiology associated with CHD8 mutations.

RAS/MAPK gene dysfunction underlies various cancers and neurocognitive disorders. Although the roles of RAS/MAPK genes have been well studied in cancer, less is known about their function during neurodevelopment. There are many genes that work in concert to regulate RAS/MAPK signaling, suggesting that if common brain phenotypes could be discovered they could have a broad impact on the many other disorders caused by distinct RAS/MAPK genes. We assessed the cellular and molecular consequences of hyperactivating the RAS/MAPK pathway using two distinct genes in a cell type previously implicated in RAS/MAPK-mediated cognitive changes, cortical GABAergic interneurons. We uncovered some GABAergic core programs that are commonly altered in each of the mutants. Notably, hyperactive RAS/MAPK mutants bias developing cortical interneurons towards those that are somatostatin positive. The increase in somatostatin-positive interneurons could also be prevented by pharmacological inhibition of the core RAS/MAPK signaling pathway. Overall, these findings present new insights into how different RAS/MAPK mutations can converge on GABAergic interneurons, which may be important for other RAS/MAPK genes and related disorders.

Epidemiological and clinical evidence points to cancer as a comorbidity in people with autism spectrum disorders (ASD). A significant overlap of genes and biological processes between both diseases has also been reported.