DOCSLIB.ORG

Within-Breed and Multi-Breed GWAS on Imputed Whole-Genome Sequence Variants Reveal Candidate Mutations Affecting Milk Protein Co

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Within-Breed and Multi-Breed GWAS on Imputed Whole-Genome Sequence Variants Reveal Candidate Mutations Affecting Milk Protein Co Sanchez et al. Genet Sel Evol (2017) 49:68 DOI 10.1186/s12711-017-0344-z Genetics Selection Evolution RESEARCH ARTICLE Open Access Within‑breed and multi‑breed GWAS on imputed whole‑genome sequence variants reveal candidate mutations afecting milk protein composition in dairy cattle Marie‑Pierre Sanchez1*, Armelle Govignon‑Gion1,2, Pascal Croiseau1, Sébastien Fritz1,3, Chris Hozé1,3, Guy Miranda1, Patrice Martin1, Anne Barbat‑Leterrier1, Rabia Letaïef1, Dominique Rocha1, Mickaël Brochard2, Mekki Boussaha1 and Didier Boichard1 Abstract Background: Genome-wide association studies (GWAS) were performed at the sequence level to identify candidate mutations that afect the expression of six major milk proteins in Montbéliarde (MON), Normande (NOR), and Holstein (HOL) dairy cattle. Whey protein (α-lactalbumin and β-lactoglobulin) and casein (αs1, αs2, β, and κ) contents were estimated by mid-infrared (MIR) spectrometry, with medium to high accuracy (0.59 R2 0.92), for 848,068 test-day milk samples from 156,660 cows in the frst three lactations. Milk composition was evaluated≤ ≤ as average test-day measurements adjusted for environmental efects. Next, we genotyped a subset of 8080 cows (2967 MON, 2737 NOR, and 2306 HOL) with the BovineSNP50 Beadchip. For each breed, genotypes were frst imputed to high-density (HD) using HD single nucleotide polymorphisms (SNPs) genotypes of 522 MON, 546 NOR, and 776 HOL bulls. The resulting HD SNP genotypes were subsequently imputed to the sequence level using 27 million high-quality sequence variants selected from Run4 of the 1000 Bull Genomes consortium (1147 bulls). Within-breed, multi-breed, and conditional GWAS were performed. Results: Thirty-four distinct genomic regions were identifed. Three regions on chromosomes 6, 11, and 20 had very signifcant efects on milk composition and were shared across the three breeds. Other signifcant efects, which partially overlapped across breeds, were found on almost all the autosomes. Multi-breed analyses provided a larger number of signifcant genomic regions with smaller confdence intervals than within-breed analyses. Combinations of within-breed, multi-breed, and conditional analyses led to the identifcation of putative causative variants in several candidate genes that presented signifcant protein–protein interactions enrichment, including those with previously described efects on milk composition (SLC37A1, MGST1, ABCG2, CSN1S1, CSN2, CSN1S2, CSN3, PAEP, DGAT1, AGPAT6) and those with efects reported for the frst time here (ALPL, ANKH, PICALM). Conclusions: GWAS applied to fne-scale phenotypes, multiple breeds, and whole-genome sequences seems to be efective to identify candidate gene variants. However, although we identifed functional links between some candi‑ date genes and milk phenotypes, the causality between candidate variants and milk protein composition remains to be demonstrated. Nevertheless, the identifcation of potential causative mutations that underlie milk protein compo‑ sition may have immediate applications for improvements in cheese-making. *Correspondence: marie‑[email protected] 1 GABI, INRA, AgroParisTech, Université Paris Saclay, 78350 Jouy‑en‑Josas, France Full list of author information is available at the end of the article © The Author(s) 2017. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/ publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated. Sanchez et al. Genet Sel Evol (2017) 49:68 Page 2 of 16 Background the animals and milk analyses are in Sanchez et al. [8]. In cattle, milk protein composition is mostly infuenced Briefy, MIR spectra were obtained for 848,068 milk sam- by genetic factors [1–4] and is of interest because it ples from 156,660 cows of the three main French dairy determines cheese-making properties [5]. Bovine milk breeds: Montbéliarde (MON), Normande (NOR), and protein composition can be predicted at a large scale by Holstein (HOL). Tese spectra were used to predict milk analyzing mid-infrared (MIR) spectra, which is routinely protein content (PC) and milk protein composition with performed [6, 7]. Combined with cow genotyping, this the equations derived as described by Ferrand et al. [7]. technique may open avenues to investigate the genomic More details about the method and the calibration popu- regions that infuence milk protein composition. In a pre- lation used are in Sanchez et al. [4]. Te contents of the vious genome-wide association study (GWAS) based on six main milk proteins (αs1-CN, αs2-CN, β-CN, κ-CN, the bovine 50 K single nucleotide polymorphism (SNP) α-LA, and β-LG) were predicted in g/100 g protein. array, we highlighted numerous genomic regions with Total casein content and total whey protein content were very signifcant efects on milk protein composition in the also analyzed (Σ-CN and Σ-WP, respectively). In order three main breeds of French dairy cattle: Holstein (HOL), to adjust phenotypes for non-genetic efects, a within- Montbéliarde (MON), and Normande (NOR) [8]. How- breed mixed model was applied to test-day data using the ever, because the 50 K SNP array contains only a small GENEKIT software [11]. Tis single-trait repeatability fraction of the total number of genomic variants, we were model included genetic, permanent environmental, and not able to directly pinpoint candidate mutations. residual random efects, as well as herd × test-day, par- In Run4 of the 1000 bull genome reference popula- ity × stage of lactation, year × month of calving, and tion, a database containing more than 56 million SNPs spectrometer × test month fxed efects. We applied this and small insertions/deletions (InDel) was constructed model to data from the frst two lactations that included by analyzing whole-genome sequences (WGS) from 1147 at least three test-day records across lactations during bulls representing 27 diferent breeds, including 288 the study period. Ten, test-day data were corrected for HOL, 28 MON and 24 NOR bulls. Tese data can then all non-genetic efects included in the model and aver- be used to impute WGS from experimentally or routinely aged per cow. Tus, for each trait and each cow, a sin- obtained 50 K SNP genotypes [9]. In this way, imputed gle phenotype was defned and subsequently used in WGS can be obtained for a large number of animals and GWAS analyses. In total, 293,780, 58,594, and 72,973 in particular, those with phenotypes. test-day records were analyzed, which corresponded to Since WGS contain almost all the genomic variants, they 44,959 MON, 12,428 NOR, and 14,530 HOL cows, i.e. should contain the causal mutations for a given trait and, an average of 6.5, 4.7, and 5.0 test-day records per cow, thus they provide a much higher GWAS resolution. How- respectively. ever, due to the long-range linkage disequilibrium that Among these cows, 8010 were genotyped with the exists within dairy cattle breeds, the resolution of within- Illumina BovineSNP50 BeadChip (Illumina Inc., San breed GWAS is often limited. For causal mutations that Diego). We applied the following quality control flters: are shared among breeds, a multi-breed model can be used the individual call rate had to be higher than 95%, the to refne regions that harbour quantitative trait loci (QTL). SNP call rate higher than 90%, the minor allele frequency Tis approach takes advantage of the historical recombi- (MAF) higher than 5%, and genotype frequencies had to nation events that have occurred in each breed, resulting be in Hardy–Weinberg equilibrium with P > 10−4. Te in linkage disequilibrium over shorter distances and better fnal dataset included between 37,332 and 41,028 SNPs resolution [10]. (Table 1), depending on the within-breed or multi-breed Here, we report the results of a GWAS at the sequence population considered, for 7907 cows (3032 MON, 2659 level for six major milk protein contents, namely NOR, and 2216 HOL) with phenotypes. α-lactalbumin and β-lactoglobulin and αs1, αs2, β, and κ caseins from HOL, MON, and NOR cows. Te results of Imputation to whole‑genome sequences within-breed, multi-breed, and conditional analyses, that Te 50 K SNP genotypes of the 7907 cows were imputed ft the most signifcant variant in addition to other tested to whole-genome sequence (WGS) using FImpute soft- variants, are examined together in order to pinpoint ware, which accurately and quickly analyzes large data- potential candidate variants in each genomic region. sets [12]. A two-step approach was applied in order to improve the accuracy of results: from 50 to 777 K high- Methods density (HD) SNPs, and then, from imputed HD SNPs Animals, phenotypes, and genotypes to WGS [13]. All imputations were performed separately For this study, we did not perform any animal experi- for each breed using either a breed-specifc (from 50 K ment, thus no ethical approval was required. Details on to HD SNPs) or a multi-breed (from HD SNPs to WGS) Sanchez et al. Genet Sel Evol (2017) 49:68 Page 3 of 16 Table 1 Features of the Montbéliarde (MON), Normande with 721 additional markers; and (2) a group of 2142 (NOR), Holstein (HOL), and multi-breed populations Montbéliarde cows that were genotyped with the fourth Number of MON NOR HOL Multi‑breed version (V4) of the same EuroG10k chip containing 3082 additional SNPs. Only SNPs with good technical qual- Phenotyped cows 44,959 12,428 14,530 71,917 ity (call rate > 95%, validation of the clusters by visual Total test‑day 293,780 58,594 72,973 425,347 records inspection, within-breed allelic frequency not signif- Test‑day records per 6.5 4.7 5 5.9 cantly diferent across chip versions) were used.
Load more

Recommended publications
  • Effects on Production Traits of Haplotypes Among Casein Genes in Norwegian Goats and Evidence for a Site of Preferential Recombination
    Copyright Ó 2006 by the Genetics Society of America DOI: 10.1534/genetics.106.058966 Effects on Production Traits of Haplotypes Among Casein Genes in Norwegian Goats and Evidence for a Site of Preferential Recombination Ben Hayes,1,2 Nina Hagesæther,1 Tormod A˚ dnøy, Grunde Pellerud, Paul R. Berg and Sigbjørn Lien Centre for Integrative Genetics and Department of Animal and Aquacultural Sciences, Norwegian University of Life Sciences, N-1432 A˚ s, Norway Manuscript received April 5, 2006 Accepted for publication July 2, 2006 ABSTRACT In goat milk the most abundant proteins are the casein genes, CSN1S1, CSN2, CSN1S2, and CSN3. Mutations have been identified within these genes affecting the level of gene expression, and effects on milk production traits have been reported. The aim of this study was to detect polymorphisms (SNPs) in the casein genes of Norwegian goats, resolve haplotype structures within the loci, and assess the effect of these haplotypes on milk production traits. Four hundred thirty-six Norwegian bucks were genotyped for 39 polymorphic sites across the four loci. The numbers of unique haplotypes present in each locus were 10, 6, 4, and 8 for CSN1S1, CSN2, CSN1S2, and CSN3, respectively. The effects of the CSN1S1 haplotypes on protein percentage and fat kilograms were significant, as were the effects of CSN3 haplotypes on fat percentage and protein percentage. A deletion in exon 12 of CSN1S1, unique to the Norwegian goat popu- lation, explained the effects of CSN1S1 haplotypes on fat kilograms, but not protein percentage. Inves- tigation of linkage disequilibrium between all possible pairs of SNPs revealed higher levels of linkage disequilbrium for SNP pairs within casein loci than for SNP pairs between casein loci, likely reflecting low levels of intragenic recombination.
    [Show full text]
  • Datasheet Blank Template
    SANTA CRUZ BIOTECHNOLOGY, INC. a-S1-casein (D-10): sc-373711 BACKGROUND PRODUCT a-S1-casein, also known as CSN1S1, CSN1 or CASA, is a 185 amino acid Each vial contains 200 µg IgG1 kappa light chain in 1.0 ml of PBS with < 0.1% secreted protein that is mammary gland-specific and belongs to the a-casein sodium azide and 0.1% gelatin. family. Existing as a disulfide-linked heterodimer with k-casein, a-S1-casein plays an important role in the ability of milk to transport calcium phosphate, APPLICATIONS a family of minerals that are key components of bone and teeth. a-S1-casein a-S1-casein (D-10) is recommended for detection of a-S1-casein of mouse exists as multiple alternatively spliced isoforms and is encoded by a gene and rat origin by Western Blotting (starting dilution 1:100, dilution range which maps to a region on human chromosome 4 that encodes other casein 1:100-1:1000), immunoprecipitation [1-2 µg per 100-500 µg of total protein family members. Chromosome 4 houses nearly 6% of the human genome (1 ml of cell lysate)], immunofluorescence (starting dilution 1:50, dilution and has the largest gene deserts (regions of the genome with no protein range 1:50-1:500) and solid phase ELISA (starting dilution 1:30, dilution encoding genes) of all of the human chromosomes. Defects in some of the range 1:30-1:3000). genes located on chromosome 4 are associated with Huntington’s disease, Ellis-van Creveld syndrome, methylmalonic acidemia and polycystic kidney Suitable for use as control antibody for a-S1-casein siRNA (m): sc-72411, disease.
    [Show full text]
  • Polymorphisms of the CSN1S1 Gene and Its Protein Variants in River and Swamp Buffalo (Bubalus Bubalis)
    Pakistan J. Zool., pp 1-10, 2021. DOI: https://dx.doi.org/10.17582/journal.pjz/20191202121250 Polymorphisms of the CSN1S1 Gene and its Protein Variants in River and Swamp Buffalo (Bubalus bubalis) Yongyun Zhang1, 2, Xinyang Fan1, Fangting Zhou1, Weizhen Li3, Yina Ouyang1,4 and Yongwang Miao1* 1Faculty of Animal Science and Technology, Yunnan Agricultural University, Kunming, Yunnan, 650201, China 2Teaching Demonstration Center of the Basic Experiments of Agricultural Majors, Yunnan Agricultural University, Kunming, Yunnan, 650201, China 3College of Veterinary Medicine, Yunnan Agricultural University, Kunming, Yunnan, 650201, China 4Yunnan Animal Science and Veterinary Institute, Kunming, Yunnan, 650224, China Article Information Yongyun Zhang and Xinyang Fan are the co-first authors. Received 02 December 2019 Revised 11 January 2020 Accepted 14 January 2020 ABSTRACT Available online 13 August 2020 Alpha S1-casein (α -CN) is a major casein in milk, which exerts a crucial role in casein transport and is Authors’ Contribution S1 YM and YZ conceived and designed related to individual milk components, nutritive value and production traits of milk. So far, α -CN coding S1 the research. YZ and XF performed gene (CSN1S1) has been widely studied in dairy cattle, but the polymorphisms of the CSN1S1 gene have the material preparation and not been fully understood in buffalo. In this study, the polymorphisms in coding sequence (CDS) of the experiments. FZ and WL performed CSN1S1 gene for river and swamp buffalo were detected using PCR product direct sequencing. The the data collection and analysis. YZ CDS for both types of buffalo was the same in length, which contained 645 nucleotides and encoded a and YM drafted the manuscript.
    [Show full text]
  • Rabbit Anti-Alpha S1 Casein/FITC Conjugated Antibody
    SunLong Biotech Co.,LTD Tel: 0086-571- 56623320 Fax:0086-571- 56623318 E-mail:[email protected] www.sunlongbiotech.com Rabbit Anti-Alpha s1 Casein/FITC Conjugated antibody SL10033R-FITC Product Name: Anti-Alpha s1 Casein/FITC Chinese Name: FITC标记的α-s1酪蛋白抗体 Alpha s1 Casein; alpha-S1-casein precursor; Alpha-S1-casein; CSN1S1; CSN1; Bos d Alias: 8; Antioxidant peptide; CASA1_BOVIN. Organism Species: Rabbit Clonality: Polyclonal React Species: Cow,Sheep, Applications: not yet tested in other applications. optimal dilutions/concentrations should be determined by the end user. Molecular weight: 23kDa Form: Lyophilized or Liquid Concentration: 1mg/ml immunogen: KLH conjugated synthetic peptide derived from cow Alpha s1 Casein Lsotype: IgG Purification: affinitywww.sunlongbiotech.com purified by Protein A Storage Buffer: 0.01M TBS(pH7.4) with 1% BSA, 0.03% Proclin300 and 50% Glycerol. Store at -20 ∑C for one year. Avoid repeated freeze/thaw cycles. The lyophilized antibody is stable at room temperature for at least one month and for greater than a year Storage: when kept at -20∑C. When reconstituted in sterile pH 7.4 0.01M PBS or diluent of antibody the antibody is stable for at least two weeks at 2-4 ∑C. background: alpha-S1-casein, also known as CSN1S1, CSN1 or CASA, is a 185 amino acid secreted protein that is mammary gland-specific and belongs to the alpha-casein family. Existing Product Detail: as a disulfide-linked heterodimer with ⅹ-casein, alpha-S1-casein plays an important role in the ability of milk to transport calcium phosphate, a family of minerals that are key components of bone and teeth.
    [Show full text]
  • Variant Discovery in the Sheep Milk Transcriptome Using RNA Sequencing
    Suárez-Vega et al. BMC Genomics (2017) 18:170 DOI 10.1186/s12864-017-3581-1 RESEARCH ARTICLE Open Access Variant discovery in the sheep milk transcriptome using RNA sequencing Aroa Suárez-Vega1, Beatriz Gutiérrez-Gil1, Christophe Klopp2, Gwenola Tosser-Klopp3 and Juan José Arranz1* Abstract Background: The identification of genetic variation underlying desired phenotypes is one of the main challenges of current livestock genetic research. High-throughput transcriptome sequencing (RNA-Seq) offers new opportunities for the detection of transcriptome variants (SNPs and short indels) in different tissues and species. In this study, we used RNA-Seq on Milk Sheep Somatic Cells (MSCs) with the goal of characterizing the genetic variation within the coding regions of the milk transcriptome in Churra and Assaf sheep, two common dairy sheep breeds farmed in Spain. Results: A total of 216,637 variants were detected in the MSCs transcriptome of the eight ewes analyzed. Among them, a total of 57,795 variants were detected in the regions harboring Quantitative Trait Loci (QTL) for milk yield, protein percentage and fat percentage, of which 21.44% were novel variants. Among the total variants detected, 561 (2.52%) and 1,649 (7.42%) were predicted to produce high or moderate impact changes in the corresponding transcriptional unit, respectively. In the functional enrichment analysis of the genes positioned within selected QTL regions harboring novel relevant functional variants (high and moderate impact), the KEGG pathway with the highest enrichment was “protein processing in endoplasmic reticulum”. Additionally, a total of 504 and 1,063 variants were identified in the genes encoding principal milk proteins and molecules involved in the lipid metabolism, respectively.
    [Show full text]
  • Novel Combination of Huanglianjiedu Decoction in Behavioral and Psychological Symptoms of Dementia in Alzheimer's Disease
    Novel Combination of HuanglianJiedu Decoction in Behavioral and Psychological Symptoms of Dementia in Alzheimer's Disease Su-Ya Ma Beijing University of Chinese Medicine Xu Wang Beijing University of Chinese Medicine Jing Shi ( [email protected] ) Department of Neurology, Dongzhimen Hospital, Beijing University of Chinese Medicine, Beijing, 100007, China Research Article Keywords: Combination of traditional Chinese medicine, Behavioral and psychological symptoms of dementia, HuanglianJiedu decoction, Network pharmacology, Molecular docking Posted Date: January 19th, 2021 DOI: https://doi.org/10.21203/rs.3.rs-124432/v2 License: This work is licensed under a Creative Commons Attribution 4.0 International License. Read Full License Page 1/15 Abstract Alzheimer's disease (AD) is characterized by progressive cognitive decline. Besides cognitive decit, AD is also characterized by behavioral and psychological symptoms in dementia (BPSD). However, therapeutic management of BPSD remains challenging. HuanglianJiedu decoction (HLJDD), a traditional Chinese prescription, consisting of four herbs, is applied to treat AD, especially AD with BPSD. Though HLJDD, has the traditional combination with the principal herb Coptidis rhizoma (Huang-lian), it might, however, not be suitable for treating BPSD. Elucidating the mechanism underlying each herb is critical to the disease-matched combination of HLJDD. In this study, network pharmacology was used to determine the targets and biological processes regulated by HLJDD in the treatment of BPSD. Moreover, molecular docking was utilized to evaluate the binding activity between the herbs' main active ingredients and neurotransmitter receptors. The results showed that Scutellariae radix (Huang-qin) and Phellodendri chinrnsis cortex (Huang-bai) exhibited better anti-BPSD effects when compared to Coptidis rhizoma and Gardeniae fructus (Zhi-zi).
    [Show full text]
  • Us 2018 / 0305689 A1
    US 20180305689A1 ( 19 ) United States (12 ) Patent Application Publication ( 10) Pub . No. : US 2018 /0305689 A1 Sætrom et al. ( 43 ) Pub . Date: Oct. 25 , 2018 ( 54 ) SARNA COMPOSITIONS AND METHODS OF plication No . 62 /150 , 895 , filed on Apr. 22 , 2015 , USE provisional application No . 62/ 150 ,904 , filed on Apr. 22 , 2015 , provisional application No. 62 / 150 , 908 , (71 ) Applicant: MINA THERAPEUTICS LIMITED , filed on Apr. 22 , 2015 , provisional application No. LONDON (GB ) 62 / 150 , 900 , filed on Apr. 22 , 2015 . (72 ) Inventors : Pål Sætrom , Trondheim (NO ) ; Endre Publication Classification Bakken Stovner , Trondheim (NO ) (51 ) Int . CI. C12N 15 / 113 (2006 .01 ) (21 ) Appl. No. : 15 /568 , 046 (52 ) U . S . CI. (22 ) PCT Filed : Apr. 21 , 2016 CPC .. .. .. C12N 15 / 113 ( 2013 .01 ) ; C12N 2310 / 34 ( 2013. 01 ) ; C12N 2310 /14 (2013 . 01 ) ; C12N ( 86 ) PCT No .: PCT/ GB2016 /051116 2310 / 11 (2013 .01 ) $ 371 ( c ) ( 1 ) , ( 2 ) Date : Oct . 20 , 2017 (57 ) ABSTRACT The invention relates to oligonucleotides , e . g . , saRNAS Related U . S . Application Data useful in upregulating the expression of a target gene and (60 ) Provisional application No . 62 / 150 ,892 , filed on Apr. therapeutic compositions comprising such oligonucleotides . 22 , 2015 , provisional application No . 62 / 150 ,893 , Methods of using the oligonucleotides and the therapeutic filed on Apr. 22 , 2015 , provisional application No . compositions are also provided . 62 / 150 ,897 , filed on Apr. 22 , 2015 , provisional ap Specification includes a Sequence Listing . SARNA sense strand (Fessenger 3 ' SARNA antisense strand (Guide ) Mathew, Si Target antisense RNA transcript, e . g . NAT Target Coding strand Gene Transcription start site ( T55 ) TY{ { ? ? Targeted Target transcript , e .
    [Show full text]
  • The Spotted Gar Genome Illuminates Vertebrate Evolution and Facilitates Humanteleost Comparisons
    ARTICLES OPEN The spotted gar genome illuminates vertebrate evolution and facilitates human-teleost comparisons Ingo Braasch1,34, Andrew R Gehrke2, Jeramiah J Smith3, Kazuhiko Kawasaki4, Tereza Manousaki5, Jeremy Pasquier6, Angel Amores1, Thomas Desvignes1, Peter Batzel1, Julian Catchen7, Aaron M Berlin8, Michael S Campbell9,34, Daniel Barrell10,11, Kyle J Martin12,34, John F Mulley13, Vydianathan Ravi14, Alison P Lee14, Tetsuya Nakamura2, Domitille Chalopin15,34, Shaohua Fan16,34, Dustin Wcisel17,18, Cristian Cañestro19,20, Jason Sydes1, Felix E G Beaudry21, Yi Sun22,23, Jana Hertel24,34, Michael J Beam1, Mario Fasold24,34, Mikio Ishiyama25, Jeremy Johnson8, Steffi Kehr24, Marcia Lara8, John H Letaw1, Gary W Litman26, Ronda T Litman26, Masato Mikami27, Tatsuya Ota28, Nil Ratan Saha29, Louise Williams8, Peter F Stadler24, Han Wang22,23, John S Taylor21, Quenton Fontenot30, Allyse Ferrara30, Stephen M J Searle10, Bronwen Aken10,11, Mark Yandell9, Igor Schneider31, Jeffrey A Yoder17,18, Jean-Nicolas Volff15, Axel Meyer16,32, Chris T Amemiya29, Byrappa Venkatesh14, Peter W H Holland12, Yann Guiguen6, Julien Bobe6, Neil H Shubin2, Federica Di Palma8,34, Jessica Alföldi8, Kerstin Lindblad-Toh8,33 & John H Postlethwait1 To connect human biology to fish biomedical models, we sequenced the genome of spotted gar (Lepisosteus oculatus), whose lineage diverged from teleosts before teleost genome duplication (TGD). The slowly evolving gar genome has conserved in content and size many entire chromosomes from bony vertebrate ancestors. Gar bridges teleosts to tetrapods by illuminating the evolution of immunity, mineralization and development (mediated, for example, by Hox, ParaHox and microRNA genes). Numerous conserved noncoding elements (CNEs; often cis regulatory) undetectable in direct human-teleost comparisons become apparent using gar: functional studies uncovered conserved roles for such cryptic CNEs, facilitating annotation of sequences Nature America, Inc.
    [Show full text]
  • UC Davis UC Davis Previously Published Works
    UC Davis UC Davis Previously Published Works Title G-NEST: a gene neighborhood scoring tool to identify co-conserved, co-expressed genes Permalink https://escholarship.org/uc/item/8991763q Journal BMC Bioinformatics, 13(1) ISSN 1471-2105 Authors Lemay, Danielle G Martin, William F Hinrichs, Angie S et al. Publication Date 2012-09-28 DOI http://dx.doi.org/10.1186/1471-2105-13-253 Supplemental Material https://escholarship.org/uc/item/8991763q#supplemental Peer reviewed eScholarship.org Powered by the California Digital Library University of California Lemay et al. BMC Bioinformatics 2012, 13:253 http://www.biomedcentral.com/1471-2105/13/253 RESEARCH ARTICLE Open Access G-NEST: a gene neighborhood scoring tool to identify co-conserved, co-expressed genes Danielle G Lemay1*, William F Martin2, Angie S Hinrichs3, Monique Rijnkels4, J Bruce German2, Ian Korf1 and Katherine S Pollard5 Abstract Background: In previous studies, gene neighborhoods—spatial clusters of co-expressed genes in the genome—have been defined using arbitrary rules such as requiring adjacency, a minimum number of genes, a fixed window size, or a minimum expression level. In the current study, we developed a Gene Neighborhood Scoring Tool (G-NEST) which combines genomic location, gene expression, and evolutionary sequence conservation data to score putative gene neighborhoods across all possible window sizes simultaneously. Results: Using G-NEST on atlases of mouse and human tissue expression data, we found that large neighborhoods of ten or more genes are extremely rare in mammalian genomes. When they do occur, neighborhoods are typically composed of families of related genes. Both the highest scoring and the largest neighborhoods in mammalian genomes are formed by tandem gene duplication.
    [Show full text]
  • Evolution of Cellular Diversity in Primary Motor Cortex of Human, Marmoset Monkey, and Mouse
    bioRxiv preprint doi: https://doi.org/10.1101/2020.03.31.016972; this version posted April 4, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. 1 Evolution of cellular diversity in primary motor cortex of human, marmoset monkey, and mouse 2 3 Trygve E. Bakken, Nikolas L. Jorstad, Qiwen Hu, Blue B. Lake, Wei Tian, Brian E. Kalmbach, Megan 4 Crow, Rebecca D. Hodge, Fenna M. Krienen, Staci A. Sorensen, Jeroen Eggermont, Zizhen Yao, Brian 5 D. Aevermann, Andrew I. Aldridge, Anna Bartlett, Darren Bertagnolli, Tamara Casper, Rosa G. 6 Castanon, Kirsten Crichton, Tanya L. Daigle, Rachel Dalley, Nick Dee, Nikolai Dembrow, Dinh Diep, 7 Song-Lin Ding, Weixiu Dong, Rongxin Fang, Stephan Fischer, Melissa Goldman, Jeff Goldy, Lucas T. 8 Graybuck, Brian R. Herb, Xiaomeng Hou, Jayaram Kancherla, Matthew Kroll, Kanan Lathia, Baldur van 9 Lew, Yang Eric Li, Christine S. Liu, Hanqing Liu, Jacinta D. Lucero, Anup Mahurkar, Delissa McMillen, 10 Jeremy A. Miller, Marmar Moussa, Joseph R. Nery, Philip R. Nicovich, Joshua Orvis, Julia K. Osteen, 11 Scott Owen, Carter R. Palmer, Thanh Pham, Nongluk Plongthongkum, Olivier Poirion, Nora M. Reed, 12 Christine Rimorin, Angeline Rivkin, William J. Romanow, Adriana E. Sedeño-Cortés, Kimberly Siletti, 13 Saroja Somasundaram, Josef Sulc, Michael Tieu, Amy Torkelson, Herman Tung, Xinxin Wang, 14 Fangming Xie, Anna Marie Yanny, Renee Zhang, Seth A. Ament, M. Margarita Behrens, Hector 15 Corrada Bravo, Jerold Chun, Alexander Dobin, Jesse Gillis, Ronna Hertzano, Patrick R.
    [Show full text]
  • Case Report Application of Microarray-Based Comparative Genomic Hybridization in Prenatal and Postnatal Settings: Three Case Reports
    SAGE-Hindawi Access to Research Genetics Research International Volume 2011, Article ID 976398, 9 pages doi:10.4061/2011/976398 Case Report Application of Microarray-Based Comparative Genomic Hybridization in Prenatal and Postnatal Settings: Three Case Reports Jing Liu, Francois Bernier, Julie Lauzon, R. Brian Lowry, and Judy Chernos Department of Medical Genetics, University of Calgary, 2888 Shaganappi Trail NW, Calgary, AB, Canada T3B 6A8 Correspondence should be addressed to Judy Chernos, [email protected] Received 16 February 2011; Revised 20 April 2011; Accepted 20 May 2011 Academic Editor: Reha Toydemir Copyright © 2011 Jing Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Microarray-based comparative genomic hybridization (array CGH) is a newly emerged molecular cytogenetic technique for rapid evaluation of the entire genome with sub-megabase resolution. It allows for the comprehensive investigation of thousands and millions of genomic loci at once and therefore enables the efficient detection of DNA copy number variations (a.k.a, cryptic genomic imbalances). The development and the clinical application of array CGH have revolutionized the diagnostic process in patients and has provided a clue to many unidentified or unexplained diseases which are suspected to have a genetic cause. In this paper, we present three clinical cases in both prenatal and postnatal settings. Among all, array CGH played a major discovery role to reveal the cryptic and/or complex nature of chromosome arrangements. By identifying the genetic causes responsible for the clinical observation in patients, array CGH has provided accurate diagnosis and appropriate clinical management in a timely and efficient manner.
    [Show full text]
  • UC Davis UC Davis Previously Published Works
    UC Davis UC Davis Previously Published Works Title From genes to milk: genomic organization and epigenetic regulation of the mammary transcriptome. Permalink https://escholarship.org/uc/item/4gz060ns Journal PloS one, 8(9) ISSN 1932-6203 Authors Lemay, Danielle G Pollard, Katherine S Martin, William F et al. Publication Date 2013 DOI 10.1371/journal.pone.0075030 Peer reviewed eScholarship.org Powered by the California Digital Library University of California From Genes to Milk: Genomic Organization and Epigenetic Regulation of the Mammary Transcriptome Danielle G. Lemay1*, Katherine S. Pollard3, William F. Martin2, Courtneay Freeman Zadrowski4, Joseph Hernandez4, Ian Korf1, J. Bruce German2, Monique Rijnkels4* 1 Genome Center, University of California Davis, Davis, California, United States of America, 2 Department of Food Science and Technology, University of California Davis, Davis, California, United States of America, 3 Gladstone Institutes, Institute for Human Genetics, and Department of Epidemiology and Biostatistics, University of California San Francisco, San Francisco, California, United States of America, 4 USDA/ARS Children’s Nutrition Research Center, Department of Pediatrics, Baylor College of Medicine, Houston, Texas, United States of America Abstract Even in genomes lacking operons, a gene’s position in the genome influences its potential for expression. The mechanisms by which adjacent genes are co-expressed are still not completely understood. Using lactation and the mammary gland as a model system, we explore the hypothesis that chromatin state contributes to the co-regulation of gene neighborhoods. The mammary gland represents a unique evolutionary model, due to its recent appearance, in the context of vertebrate genomes. An understanding of how the mammary gland is regulated to produce milk is also of biomedical and agricultural importance for human lactation and dairying.
    [Show full text]