A pan-transcriptome analysis shows that disease resistance genes have undergone more selection pressure during barley domestication

BackgroundIt has become clear in recent years that many genes in a given species may not be found in a single genotype thus using sequences from a single genotype as reference may not be adequate for various applications.ResultsIn this study we constructed a pan-transcriptome for barley by de novo assembling 288 sets of RNA-seq data from 32 cultivated barley genotypes and 31 wild barley genotypes. The pan-transcriptome consists of 756,632 transcripts with an average N50 length of 1240 bp. Of these, 289,697 (38.2%) were not found in the genome of the international reference genotype Morex. The novel transcripts are enriched with genes associated with responses to different stresses and stimuli. At the pan-transcriptome level, genotypes of wild barley have a higher proportion of disease resistance genes than cultivated ones.ConclusionsWe demonstrate that the use of the pan-transcriptome dramatically improved the efficiency in detecting variation in barley. Analysing the pan-transcriptome also found that, compared with those in other categories, disease resistance genes have gone through stronger selective pressures during domestication.

[1]  Patrick J. Biggs,et al.  SolexaQA: At-a-glance quality assessment of Illumina second-generation sequencing data , 2010, BMC Bioinformatics.

[2]  Robert D. Finn,et al.  Rfam 12.0: updates to the RNA families database , 2014, Nucleic Acids Res..

[3]  Monika S. Doblin,et al.  De novo transcriptome assembly and analysis of differentially expressed genes of two barley genotypes reveal root-zone-specific responses to salt exposure , 2016, Scientific Reports.

[4]  Sean R. Eddy,et al.  Infernal 1.1: 100-fold faster RNA homology searches , 2013, Bioinform..

[5]  A. El-Osta,et al.  Analysis of the barley leaf transcriptome under salinity stress using mRNA-Seq , 2013, Acta Physiologiae Plantarum.

[6]  M. A. Pedraza,et al.  Insights into the Maize Pan-Genome and Pan-Transcriptome[W][OPEN] , 2014, Plant Cell.

[7]  M. Parker THE PLEIOTROPY THEORY FOR POLYMORPHISM OF DISEASE RESISTANCE GENES IN PLANTS , 1990, Evolution; international journal of organic evolution.

[8]  Jing Zhao,et al.  A maize wall-associated kinase confers quantitative resistance to head smut , 2014, Nature Genetics.

[9]  Ning Wang,et al.  Transcriptome profiling reveals mosaic genomic origins of modern cultivated barley , 2014, Proceedings of the National Academy of Sciences.

[10]  Hui Wang,et al.  AgBase: supporting functional modeling in agricultural organisms , 2010, Nucleic Acids Res..

[11]  RNA-Seq analysis identifies genes associated with differential reproductive success under drought-stress in accessions of wild barley Hordeum spontaneum , 2015, BMC Plant Biology.

[12]  A. Pankin,et al.  Global Transcriptome Profiling of Developing Leaf and Shoot Apices Reveals Distinct Genetic and Environmental Control of Floral Transition and Inflorescence Development in Barley[OPEN] , 2015, The Plant Cell.

[13]  Jaideep P. Sundaram,et al.  Genome analysis of multiple pathogenic isolates of Streptococcus agalactiae: implications for the microbial "pan-genome". , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[14]  Guo‐Liang Wang,et al.  Balancing Immunity and Yield in Crop Plants. , 2017, Trends in plant science.

[15]  Robert K Jansen,et al.  RNA-Seq analysis of the wild barley (H. spontaneum) leaf transcriptome under salt stress. , 2015, Comptes rendus biologies.

[16]  P. Gao,et al.  Transcriptome Assembly and Analysis of Tibetan Hulless Barley (Hordeum vulgare L. var. nudum) Developing Grains, with Emphasis on Quality Properties , 2014, PloS one.

[17]  Chengdao Li,et al.  Involvement of Alternative Splicing in Barley Seed Germination , 2016, PloS one.

[18]  Pablo Cingolani,et al.  © 2012 Landes Bioscience. Do not distribute. , 2022 .

[19]  Ruiqiang Li,et al.  De novo assembly of soybean wild relatives for pan-genome analysis of diversity and agronomic traits , 2014, Nature Biotechnology.

[20]  T. Unver,et al.  High-throughput transcriptome analysis of barley (Hordeum vulgare) exposed to excessive boron. , 2015, Gene.

[21]  P. Schulze-Lefert,et al.  Diversity at the Mla powdery mildew resistance locus from cultivated barley reveals sites of positive selection. , 2010, Molecular plant-microbe interactions : MPMI.

[22]  R. Wing Harvesting rice’s dispensable genome , 2015, Genome Biology.

[23]  N. Friedman,et al.  Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data , 2011, Nature Biotechnology.

[24]  M. Gerstein,et al.  RNA-Seq: a revolutionary tool for transcriptomics , 2009, Nature Reviews Genetics.

[25]  Michele Morgante,et al.  Structural variation and genome complexity: is dispensable really dispensable? , 2014, Current opinion in plant biology.

[26]  M. Rees,et al.  How did the domestication of Fertile Crescent grain crops increase their yields? , 2016, Functional ecology.

[27]  Wei Wu,et al.  NONCODE 2016: an informative and valuable data source of long non-coding RNAs , 2015, Nucleic Acids Res..

[28]  Thomas D. Wu,et al.  GMAP: a genomic mapping and alignment program for mRNA and EST sequence , 2005, Bioinform..

[29]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[30]  I. Jakobsen,et al.  Heat Stress Affects Pi-related Genes Expression and Inorganic Phosphate Deposition/Accumulation in Barley , 2016, Front. Plant Sci..

[31]  Robert D. Finn,et al.  The Pfam protein families database: towards a more sustainable future , 2015, Nucleic Acids Res..

[32]  Chun-Chao Wang,et al.  RPAN: rice pan-genome browser for ∼3000 rice genomes , 2016, Nucleic acids research.

[33]  T. Vogt Phenylpropanoid biosynthesis. , 2010, Molecular plant.

[34]  C. K. Chan,et al.  The pangenome of an agronomically important crop plant Brassica oleracea , 2016, Nature Communications.

[35]  Domestication affected the basal and induced disease resistance in common bean (Phaseolus vulgaris) , 2012, European Journal of Plant Pathology.

[36]  Cheng He,et al.  Maize pan-transcriptome provides novel insights into genome complexity and quantitative trait variation , 2015, Scientific Reports.

[37]  T. Kroj,et al.  Several wall-associated kinases participate positively and negatively in basal defense against rice blast fungus , 2016, BMC Plant Biology.

[38]  Hiroaki Sakai,et al.  Comprehensive Sequence Analysis of 24,783 Barley Full-Length cDNAs Derived from 12 Clone Libraries1[W][OA] , 2011, Plant Physiology.

[39]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[40]  M. Hood,et al.  The evolution of disease resistance and tolerance in spatially structured populations , 2012, Ecology and evolution.

[41]  Matko Bosnjak,et al.  REVIGO Summarizes and Visualizes Long Lists of Gene Ontology Terms , 2011, PloS one.

[42]  John K. McCooke,et al.  A chromosome conformation capture ordered sequence of the barley genome , 2017, Nature.

[43]  Jonathan D. G. Jones,et al.  NLR-parser: rapid annotation of plant NLR complements , 2015, Bioinform..

[44]  M. Kreitman,et al.  Fitness costs of R-gene-mediated resistance in Arabidopsis thaliana , 2003, Nature.

[45]  C. K. Chan,et al.  The pangenome of hexaploid bread wheat , 2017, The Plant journal : for cell and molecular biology.

[46]  J. Doležel,et al.  Chromosome-specific sequencing reveals an extensive dispensable genome component in wheat , 2016, Scientific Reports.

[47]  Adam Godzik,et al.  Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences , 2006, Bioinform..

[48]  M. Albrecht,et al.  Structure-function analysis of the NB-ARC domain of plant disease resistance proteins. , 2008, Journal of experimental botany.

[49]  Gonçalo R. Abecasis,et al.  The variant call format and VCFtools , 2011, Bioinform..

[50]  Paul D. Shaw,et al.  Natural variation in a homolog of Antirrhinum CENTRORADIALIS contributed to spring growth habit and environmental adaptation in cultivated barley , 2012, Nature Genetics.

[51]  Robert D. Finn,et al.  HMMER web server: 2015 update , 2015, Nucleic Acids Res..

[52]  P. Visscher,et al.  GCTA: a tool for genome-wide complex trait analysis. , 2011, American journal of human genetics.

[53]  Mihaela M. Martis,et al.  A physical, genetic and functional sequence assembly of the barley genome. , 2022 .

[54]  R. Edwards,et al.  Fast Identification and Removal of Sequence Contamination from Genomic and Metagenomic Datasets , 2011, PloS one.

[55]  Bo Yu,et al.  CDD/SPARCLE: functional classification of proteins via subfamily domain architectures , 2016, Nucleic Acids Res..

[56]  T. R. Hvidsten,et al.  Comparative transcriptomics provides insight into the evolution of cold response in Pooideae , 2017, bioRxiv.

[57]  Martin Wolfe,et al.  Genetics of Powdery Mildew Resistance in Barley , 1994 .

[58]  Dong-Lei Yang,et al.  Roles of plant hormones and their interplay in rice immunity. , 2013, Molecular plant.

[59]  Heng Li,et al.  A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data , 2011, Bioinform..

[60]  M. Bayer,et al.  The low-recombining pericentromeric region of barley restricts gene diversity and evolution but not gene expression , 2014, The Plant journal : for cell and molecular biology.

[61]  A. Diévart,et al.  LRR-containing receptors regulating plant development and defense , 2003, Development.

[62]  Zhou Du,et al.  agriGO v2.0: a GO analysis toolkit for the agricultural community, 2017 update , 2017, Nucleic Acids Res..