Network Analysis of Genome-Wide Selective Constraint Reveals a Gene Network Active in Early Fetal Brain Intolerant of Mutation.
Ontology highlight
ABSTRACT:
SUBMITTER: Choi J
PROVIDER: S-EPMC4909280 | biostudies-literature | 2016 Jun
REPOSITORIES: biostudies-literature
Publications
Network Analysis of Genome-Wide Selective Constraint Reveals a Gene Network Active in Early Fetal Brain Intolerant of Mutation.
Choi Jinmyung J Shooshtari Parisa P Samocha Kaitlin E KE Daly Mark J MJ Cotsapas Chris C
PLoS genetics 20160615 6
Using robust, integrated analysis of multiple genomic datasets, we show that genes depleted for non-synonymous de novo mutations form a subnetwork of 72 members under strong selective constraint. We further show this subnetwork is preferentially expressed in the early development of the human hippocampus and is enriched for genes mutated in neurological Mendelian disorders. We thus conclude that carefully orchestrated developmental processes are under strong constraint in early brain development ...[more]
Similar Datasets
Project description:The mutation-selection model of coding sequence evolution has received renewed attention for its use in estimating site-specific amino acid propensities and selection coefficient distributions. Two computationally tractable mutation-selection inference frameworks have been introduced: One framework employs a fixed-effects, highly parameterized maximum likelihood approach, whereas the other employs a random-effects Bayesian Dirichlet Process approach. While both implementations follow the same model, they appear to make distinct predictions about the distribution of selection coefficients. The fixed-effects framework estimates a large proportion of highly deleterious substitutions, whereas the random-effects framework estimates that all substitutions are either nearly neutral or weakly deleterious. It remains unknown, however, how accurately each method infers evolutionary constraints at individual sites. Indeed, selection coefficient distributions pool all site-specific inferences, thereby obscuring a precise assessment of site-specific estimates. Therefore, in this study, we use a simulation-based strategy to determine how accurately each approach recapitulates the selective constraint at individual sites. We find that the fixed-effects approach, despite its extensive parameterization, consistently and accurately estimates site-specific evolutionary constraint. By contrast, the random-effects Bayesian approach systematically underestimates the strength of natural selection, particularly for slowly evolving sites. We also find that, despite the strong differences between their inferred selection coefficient distributions, the fixed- and random-effects approaches yield surprisingly similar inferences of site-specific selective constraint. We conclude that the fixed-effects mutation-selection framework provides the more reliable software platform for model application and future development.
| S-EPMC5062325 | biostudies-literature
Project description:While several studies have investigated general properties of the genetic architecture of natural variation in gene expression, few of these have considered natural, outbreeding populations. In parallel, systems biology has established that a general feature of biological networks is that they are scale-free, rendering them buffered against random mutations. To date, few studies have attempted to examine the relationship between the selective processes acting to maintain natural variation of gene expression and the associated co-expression network structure. Here we utilised RNA-Sequencing to assay gene expression in winter buds undergoing bud flush in a natural population of Populus tremula, an outbreeding forest tree species. We performed expression Quantitative Trait Locus (eQTL) mapping and identified 164,290 significant eQTLs associating 6,241 unique genes (eGenes) with 147,419 unique SNPs (eSNPs). We found approximately four times as many local as distant eQTLs, with local eQTLs having significantly higher effect sizes. eQTLs were primarily located in regulatory regions of genes (UTRs or flanking regions), regardless of whether they were local or distant. We used the gene expression data to infer a co-expression network and investigated the relationship between network topology, the genetic architecture of gene expression and signatures of selection. Within the co-expression network, eGenes were underrepresented in network module cores (hubs) and overrepresented in the periphery of the network, with a negative correlation between eQTL effect size and network connectivity. We additionally found that module core genes have experienced stronger selective constraint on coding and non-coding sequence, with connectivity associated with signatures of selection. Our integrated genetics and genomics results suggest that purifying selection is the primary mechanism underlying the genetic architecture of natural variation in gene expression assayed in flushing leaf buds of P. tremula and that connectivity within the co-expression network is linked to the strength of purifying selection.
| S-EPMC5407845 | biostudies-literature
Project description:BackgroundRegions of the genome that are under evolutionary constraint across multiple species have previously been used to identify functional sequences in the human genome. Furthermore, it is known that there is an inverse relationship between evolutionary constraint and the allele frequency of a mutation segregating in human populations, implying a direct relationship between interspecies divergence and fitness in humans. Here we utilise this relationship to test differences in the accumulation of putatively deleterious mutations both between populations and on the individual level.ResultsUsing whole genome and exome sequencing data from Phase 1 of the 1000 Genome Project for 1,092 individuals from 14 worldwide populations we show that minor allele frequency (MAF) varies as a function of constraint around both coding regions and non-coding sites genome-wide, implying that negative, rather than positive, selection primarily drives the distribution of alleles among individuals via background selection. We find a strong relationship between effective population size and the depth of depression in MAF around the most conserved genes, suggesting that populations with smaller effective size are carrying more deleterious mutations, which also translates into higher genetic load when considering the number of putatively deleterious alleles segregating within each population. Finally, given the extreme richness of the data, we are now able to classify individual genomes by the accumulation of mutations at functional sites using high coverage 1000 Genomes data. Using this approach we detect differences between 'healthy' individuals within populations for the distributions of putatively deleterious rare alleles they are carrying.ConclusionsThese findings demonstrate the extent of background selection in the human genome and highlight the role of population history in shaping patterns of diversity between human individuals. Furthermore, we provide a framework for the utility of personal genomic data for the study of genetic fitness and diseases.
| S-EPMC3727949 | biostudies-literature
Project description:As predicted by the nearly neutral model of evolution, numerous studies have shown that reduced N(e) accelerates the accumulation of slightly deleterious changes under genetic drift. While such studies have mostly focused on eukaryotes, bacteria also offer excellent models to explore the effects of N(e). Most notably, the genomes of host-dependent bacteria with small N(e) show signatures of genetic drift, including elevated K(a)/K(s). Here, I explore the utility of an alternative measure of selective constraint: the per-site rate of radical and conservative amino acid substitutions (D(r)/D(c)). I test the hypothesis that purifying selection against radical amino acid changes is less effective in two insect endosymbiont groups (Blochmannia of ants and Buchnera of aphids), compared to related gamma-Proteobacteria. Genome comparisons demonstrate a significant elevation in D(r)/D(c) in endosymbionts that affects the majority (66-79%) of shared orthologs examined. The elevation of D(r)/D(c) in endosymbionts affects all functional categories examined. Simulations indicate that D(r)/D(c) estimates are sensitive to codon frequencies and mutational parameters; however, estimation biases occur in the opposite direction as the patterns observed in genome comparisons, thereby making the inference of elevated D(r)/D(c) more conservative. Increased D(r)/D(c) and other signatures of genome degradation in endosymbionts are consistent with strong effects of genetic drift in their small populations, as well as linkage to selected sites in these asexual bacteria. While relaxed selection against radical substitutions may contribute, genome-wide processes such as genetic drift and linkage best explain the pervasive elevation in D(r)/D(c) across diverse functional categories that include basic cellular processes. Although the current study focuses on a few bacterial lineages, it suggests D(r)/D(c) is a useful gauge of selective constraint and may provide a valuable alternative to K(a)/K(s) when high sequence divergences preclude estimates of K(s). Broader application of D(r)/D(c) will benefit from approaches less prone to estimation biases.
| S-EPMC3237559 | biostudies-literature
Project description:Aim: To provide a comprehensive understanding of gene regulatory networks in the developing human brain and a foundation for interpreting pathogenic deregulation. Materials & methods: We generated reference epigenomes and transcriptomes of dissected brain regions and primary neural progenitor cells (NPCs) derived from cortical and ganglionic eminence tissues of four normal human fetuses. Results: Integration of these data across developmental stages revealed a directional increase in active regulatory states, transcription factor activities and gene transcription with developmental stage. Consistent with differences in their biology, NPCs derived from cortical and ganglionic eminence regions contained common, region specific, and gestational week specific regulatory states. Conclusion: We provide a high-resolution regulatory network for NPCs from different brain regions as a comprehensive reference for future studies.
| S-EPMC7857341 | biostudies-literature
Project description:Complicated molecular and cellular processes take place in a spatiotemporally heterogeneous and precisely regulated pattern in the human fetal brain, yielding not only dramatic morphological and microstructural changes, but also macroscale connectomic transitions. As the underlying substrate of the fetal brain structural network, both dynamic neuronal migration pathways and rapid developing fetal white matter (WM) fibers could fundamentally reshape early fetal brain connectome. Quantifying structural connectome development can not only shed light on the brain reconfiguration in this critical yet rarely studied developmental period, but also reveal alterations of the connectome under neuropathological conditions. However, transition of the structural connectome from the mid-fetal stage to birth is not yet known. The contribution of different types of neural fibers to the structural network in the mid-fetal brain is not known, either. In this study, diffusion tensor magnetic resonance imaging (DT-MRI or DTI) of 10 fetal brain specimens at the age of 20 postmenstrual weeks (PMW), 12 in vivo brains at 35 PMW, and 12 in vivo brains at term (40 PMW) were acquired. The structural connectome of each brain was established with evenly parcellated cortical regions as network nodes and traced fiber pathways based on DTI tractography as network edges. Two groups of fibers were categorized based on the fiber terminal locations in the cerebral wall in the 20 PMW fetal brains. We found that fetal brain networks become stronger and more efficient during 20-40 PMW. Furthermore, network strength and global efficiency increase more rapidly during 20-35 PMW than during 35-40 PMW. Visualization of the whole brain fiber distribution by the lengths suggested that the network reconfiguration in this developmental period could be associated with a significant increase of major long association WM fibers. In addition, non-WM neural fibers could be a major contributor to the structural network configuration at 20 PMW and small-world network organization could exist as early as 20 PMW. These findings offer a preliminary record of the fetal brain structural connectome maturation from the middle fetal stage to birth and reveal the critical role of non-WM neural fibers in structural network configuration in the middle fetal stage.
| S-EPMC5645529 | biostudies-literature
Project description:The comparative genomics revolution of the past decade has enabled the discovery of functional elements in the human genome via sequence comparison. While that is so, an important class of elements, those specific to humans, is entirely missed by searching for sequence conservation across species. Here we present an analysis based on variation data among human genomes that utilizes a supervised machine learning approach for the identification of human-specific purifying selection in the genome. Using only allele frequency information from the complete low-coverage 1000 Genomes Project data set in conjunction with a support vector machine trained from known functional and nonfunctional portions of the genome, we are able to accurately identify portions of the genome constrained by purifying selection. Our method identifies previously known human-specific gains or losses of function and uncovers many novel candidates. Candidate targets for gain and loss of function along the human lineage include numerous putative regulatory regions of genes essential for normal development of the central nervous system, including a significant enrichment of gain of function events near neurotransmitter receptor genes. These results are consistent with regulatory turnover being a key mechanism in the evolution of human-specific characteristics of brain development. Finally, we show that the majority of the genome is unconstrained by natural selection currently, in agreement with what has been estimated from phylogenetic methods but in sharp contrast to estimates based on transcriptomics or other high-throughput functional methods.
| S-EPMC4700959 | biostudies-literature
Project description:Cocaine use during pregnancy is deleterious to the newborn child, in part via its disruption of placental blood flow. However, the extent to which cocaine can affect the function of the fetal primate brain is still an unresolved question. Here we used PET and MRI and show that in third-trimester pregnant nonhuman primates, cocaine at doses typically used by drug abusers significantly increased brain glucose metabolism to the same extent in the mother as in the fetus (approximately 100%). Inasmuch as brain glucose metabolism is a sensitive marker of brain function, the current findings provide evidence that cocaine use by a pregnant mother will also affect the function of the fetal brain. We are also unique in showing that cocaine's effects in brain glucose metabolism differed in pregnant (increased) and nonpregnant (decreased) animals, which suggests that the psychoactive effects of cocaine are influenced by the state of pregnancy. Our findings have clinical implications because they imply that the adverse effects of prenatal cocaine exposure to the newborn child include not only cocaine's deleterious effects to the placental circulation, but also cocaine's direct pharmacological effect to the developing fetal brain.
| S-EPMC2824386 | biostudies-other
Project description:BackgroundA gene's position in regulatory, protein interaction or metabolic networks can be predictive of the strength of purifying selection acting on it, but these relationships are neither universal nor invariably strong. Following work in bacteria, fungi and invertebrate animals, we explore the relationship between selective constraint and metabolic function in mammals.ResultsWe measure the association between selective constraint, estimated by the ratio of nonsynonymous (Ka) to synonymous (Ks) substitutions, and several, primarily metabolic, measures of gene function. We find significant differences between the selective constraints acting on enzyme-coding genes from different cellular compartments, with the nucleus showing higher constraint than genes from either the cytoplasm or the mitochondria. Among metabolic genes, the centrality of an enzyme in the metabolic network is significantly correlated with Ka/Ks. In contrast to yeasts, gene expression magnitude does not appear to be the primary predictor of selective constraint in these organisms.ConclusionsOur results imply that the relationship between selective constraint and enzyme centrality is complex: the strength of selective constraint acting on mammalian genes is quite variable and does not appear to exclusively follow patterns seen in other organisms.
| S-EPMC3082228 | biostudies-literature
Project description:Glutaminase is a multifunctional enzyme encoded by gene Gls involved in energy metabolism, ammonia trafficking and regeneration of neurotransmitter glutamate. To address the proteomic basis for the neurophenotypes of glutaminase-deficient mice, brain proteins from late gestation wild type, Gls+/- and Gls-/- male mice were subjected to two-dimensional gel electrophoresis, with subsequent identification by mass spectrometry using nano-LC-ESI-MS/MS. Protein spots that showed differential genotypic variation were quantified by immunoblotting. Differentially expressed proteins unambiguously identified by MS/MS included neurocalcin delta, retinol binding protein-1, reticulocalbin-3, cytoskeleton proteins fascin and tropomyosin alpha-4-chain, dihydropyrimidinase-related protein-5, apolipoprotein IV and proteins from protein metabolism proteasome subunits alpha type 2, type 7, heterogeneous nuclear ribonucleoprotein C1/C2 and H, voltage-gated anion-selective channel proteins 1 and 2, ATP synthase subunit ? and transitional endoplasmic reticulum ATPase. An interaction network determined by Ingenuity Pathway Analysis revealed a link between glutaminase and calcium, Akt and retinol signaling, cytoskeletal elements, ATPases, ion channels, protein synthesis and the proteasome system, intermediary, nucleic acid and lipid metabolism, huntingtin, guidance cues, transforming growth factor beta-1 and hepatocyte nuclear factor 4-alpha. The network identified involves (a) cellular assembly and organization and (b) cell signaling and cell cycle, suggesting that Gls is crucial for neuronal maturation.
| S-EPMC3874211 | biostudies-literature