Divergent evolution in the genomes of closely related lacertids, Lacerta viridis and L. bilineata, and implications for speciation

Abstract Background Lacerta viridis and Lacerta bilineata are sister species of European green lizards (eastern and western clades, respectively) that, until recently, were grouped together as the L. viridis complex. Genetic incompatibilities were observed between lacertid populations through crossing experiments, which led to the delineation of two separate species within the L. viridis complex. The population history of these sister species and processes driving divergence are unknown. We constructed the first high-quality de novo genome assemblies for both L. viridis and L. bilineata through Illumina and PacBio sequencing, with annotation support provided from transcriptome sequencing of several tissues. To estimate gene flow between the two species and identify factors involved in reproductive isolation, we studied their evolutionary history, identified genomic rearrangements, detected signatures of selection on non-coding RNA, and on protein-coding genes. Findings Here we show that gene flow was primarily unidirectional from L. bilineata to L. viridis after their split at least 1.15 million years ago. We detected positive selection of the non-coding repertoire; mutations in transcription factors; accumulation of divergence through inversions; selection on genes involved in neural development, reproduction, and behavior, as well as in ultraviolet-response, possibly driven by sexual selection, whose contribution to reproductive isolation between these lacertid species needs to be further evaluated. Conclusion The combination of short and long sequence reads resulted in one of the most complete lizard genome assemblies. The characterization of a diverse array of genomic features provided valuable insights into the demographic history of divergence among European green lizards, as well as key species differences, some of which are candidates that could have played a role in speciation. In addition, our study generated valuable genomic resources that can be used to address conservation-related issues in lacertids.

[1]  L. Bernatchez,et al.  Eco-Evolutionary Genomics of Chromosomal Inversions. , 2018, Trends in ecology & evolution.

[2]  Katja Nowick,et al.  Temporal ordering of substitutions in RNA evolution: Uncovering the structural evolution of the Human Accelerated Region 1. , 2018, Journal of theoretical biology.

[3]  Stuart R. Dennis,et al.  Comparative Genomics Reveals Accelerated Evolution in Conserved Pathways during the Diversification of Anole Lizards , 2018, Genome biology and evolution.

[4]  M. Kirkpatrick,et al.  Chromosome Inversions, Local Adaptation and Speciation , 2017, Genetics.

[5]  C. Villagra,et al.  The Importance of ncRNAs as Epigenetic Mechanisms in Phenotypic Variation and Organic Evolution , 2017, Front. Microbiol..

[6]  R. Faria,et al.  Interpreting the genomic landscape of speciation: a road map for finding barriers to gene flow , 2017, Journal of evolutionary biology.

[7]  Michael C. Schatz,et al.  Accurate detection of complex structural variations using single molecule sequencing , 2017, Nature Methods.

[8]  Richard Durbin,et al.  No evidence for maintenance of a sympatric Heliconius species barrier by chromosomal inversions , 2017, Evolution letters.

[9]  Huanming Yang,et al.  Sequencing, de novo assembling, and annotating the genome of the endangered Chinese crocodile lizard Shinisaurus crocodilurus , 2017, GigaScience.

[10]  K. Cadigan Faculty of 1000 evaluation for Stem cell signaling. An integral program for tissue renewal and regeneration: Wnt signaling and stem cell control. , 2017 .

[11]  Son K. Pham,et al.  Improved genome assembly of American alligator genome reveals conserved architecture of estrogen signaling. , 2017, Genome research.

[12]  B. Ni,et al.  Regulatory effects of autophagy on spermatogenesis , 2017, Biology of Reproduction.

[13]  Sree Rohit Raj Kolora,et al.  The complete mitochondrial genome of Lacerta bilineata and comparison with its closely related congener L. Viridis , 2015, Mitochondrial DNA. Part A, DNA mapping, sequencing, and analysis.

[14]  P. Stadler,et al.  SSS-test : A Novel Test for Detecting Selection on the Secondary Structures of non-coding RNAs , 2017 .

[15]  Q. Xi,et al.  In low protein diets, microRNA-19b regulates urea synthesis by targeting SIRT5 , 2016, Scientific Reports.

[16]  J. H. Chae,et al.  Reciprocal localization of transcription factors YY1 and CP2c in spermatogonial stem cells and their putative roles during spermatogenesis. , 2016, Acta histochemica.

[17]  H. Ellegren,et al.  Direct estimate of the rate of germline mutation in a bird , 2016, Genome research.

[18]  M. Blaxter,et al.  Para‐allopatry in hybridizing fire‐bellied toads (Bombina bombina and B. variegata): Inference from transcriptome‐wide coalescence analyses , 2016, Evolution; international journal of organic evolution.

[19]  Chao Liu,et al.  Autophagy regulates spermatid differentiation via degradation of PDLIM1 , 2016, Autophagy.

[20]  Qiye Li,et al.  Anchoring genome sequence to chromosomes of the central bearded dragon (Pogona vitticeps) enables reconstruction of ancestral squamate macrochromosomes and identifies sequence content of the Z chromosome , 2016, BMC Genomics.

[21]  P. Stadler,et al.  An updated human snoRNAome , 2016, Nucleic acids research.

[22]  Masashi Takahashi,et al.  Requirement for nuclear autoantigenic sperm protein mRNA expression in bovine preimplantation development. , 2016, Animal science journal = Nihon chikusan Gakkaiho.

[23]  V. Scaria,et al.  Genomic variations in non-coding RNAs: Structure, function and regulation. , 2016, Genomics.

[24]  A. Regev,et al.  Evolutionary analysis across mammals reveals distinct classes of long non-coding RNAs , 2015, Genome Biology.

[25]  M. Nachman,et al.  Gene regulation and speciation in house mice , 2016, Genome research.

[26]  U. Joger,et al.  Phylogeography of the Lacerta viridis complex: mitochondrial and nuclear markers provide taxonomic insights , 2016 .

[27]  Huimin Kang,et al.  Structural Variant Detection by Large-scale Sequencing Reveals New Evolutionary Evidence on Breed Divergence between Chinese and European Pigs , 2016, Scientific Reports.

[28]  Samuel S. Hunter,et al.  The population genomics of rapid adaptation: disentangling signatures of selection and demography in white sands lizards , 2016, Molecular ecology.

[29]  Daniel R. Zerbino,et al.  Ensembl 2016 , 2015, Nucleic Acids Res..

[30]  Huanming Yang,et al.  Gekko japonicus genome reveals evolution of adhesive toe pads and tail regeneration , 2015, Nature Communications.

[31]  Gabor T. Marth,et al.  An integrated map of structural variation in 2,504 human genomes , 2015, Nature.

[32]  Guojie Zhang,et al.  High-coverage sequencing and annotated assembly of the genome of the Australian dragon lizard Pogona vitticeps , 2015, GigaScience.

[33]  Laurent A. F. Frantz,et al.  Inferring Bottlenecks from Genome-Wide Samples of Short Sequence Blocks , 2015, Genetics.

[34]  T. Laisk-Podar,et al.  Genetic variants associated with female reproductive ageing--potential markers for assessing ovarian function and ovarian stimulation outcome. , 2015, Reproductive biomedicine online.

[35]  Rob W. Ness,et al.  Assessing Recent Selection and Functionality at Long Noncoding RNA Loci in the Mouse Genome , 2015, Genome biology and evolution.

[36]  M. Ritchie,et al.  Genome‐wide tests for introgression between cactophilic Drosophila implicate a role of inversions during speciation , 2015, Evolution; international journal of organic evolution.

[37]  Mark Gerstein,et al.  MetaSV: an accurate and integrative structural-variant caller for next generation sequencing , 2015, Bioinform..

[38]  Sarah W. Burge,et al.  Conservation and Losses of Non-Coding RNAs in Avian Genomes , 2015, PloS one.

[39]  Martin Chmelik,et al.  Efficient Strategies for Calculating Blockwise Likelihoods Under the Coalescent , 2015, Genetics.

[40]  Peter F. Stadler,et al.  The Expansion of Animal MicroRNA Families Revisited , 2015, Life.

[41]  Adhemar Zerlotini,et al.  POTION: an end-to-end pipeline for positive Darwinian selection detection in genome-scale data through phylogenetic comparison of protein-coding genes , 2015, BMC Genomics.

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

[43]  B. Lu,et al.  RNA metabolism in the pathogenesis of Parkinson׳s disease , 2014, Brain Research.

[44]  R. Nusse,et al.  An integral program for tissue renewal and regeneration: Wnt signaling and stem cell control , 2014, Science.

[45]  E. Font,et al.  Ultraviolet vision in lacertid lizards: evidence from retinal structure, eye transmittance, SWS1 visual pigment genes and behaviour , 2014, Journal of Experimental Biology.

[46]  B. Nickel,et al.  Illuminating the base of the annelid tree using transcriptomics. , 2014, Molecular biology and evolution.

[47]  B. Sipos,et al.  Establishment of robust controls for the normalization of miRNA expression in neuroendocrine tumors of the ileum and pancreas , 2014, Endocrine.

[48]  Åke Brännström,et al.  Genomics and the origin of species , 2014, Nature Reviews Genetics.

[49]  L. Stubbs,et al.  Deep Vertebrate Roots for Mammalian Zinc Finger Transcription Factor Subfamilies , 2014, Genome biology and evolution.

[50]  Peter F. Stadler,et al.  snoStrip: a snoRNA annotation pipeline , 2014, Bioinform..

[51]  Drew R. Schield,et al.  The Burmese python genome reveals the molecular basis for extreme adaptation in snakes , 2013, Proceedings of the National Academy of Sciences.

[52]  J. Logan,et al.  The king cobra genome reveals dynamic gene evolution and adaptation in the snake venom system , 2013, Proceedings of the National Academy of Sciences.

[53]  J. Endler,et al.  Functional characterization of spectral tuning mechanisms in the great bowerbird short-wavelength sensitive visual pigment (SWS1), and the origins of UV/violet vision in passerines and parrots , 2013, BMC Evolutionary Biology.

[54]  H. Burger,et al.  Genomic basis of aromatase excess syndrome: recombination- and replication-mediated rearrangements leading to CYP19A1 overexpression. , 2013, The Journal of clinical endocrinology and metabolism.

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

[56]  Colin N. Dewey,et al.  De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis , 2013, Nature Protocols.

[57]  Michael Ott,et al.  De novo transcript sequence reconstruction from RNA-Seq: reference generation and analysis with Trinity , 2013 .

[58]  Jie Cao,et al.  The Expression of FOXJ1 in Neurogenesis after Transient Focal Cerebral Ischemia , 2013, Canadian Journal of Neurological Sciences / Journal Canadien des Sciences Neurologiques.

[59]  Bronwen L. Aken,et al.  The draft genomes of soft–shell turtle and green sea turtle yield insights into the development and evolution of the turtle–specific body plan , 2013, Nature Genetics.

[60]  David Haussler,et al.  HAL: a hierarchical format for storing and analyzing multiple genome alignments , 2013, Bioinform..

[61]  Heng Li Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM , 2013, 1303.3997.

[62]  R. Faria,et al.  A prominent role of KRAB-ZNF transcription factors in mammalian speciation? , 2013, Trends in genetics : TIG.

[63]  S. Gabriel,et al.  Somatic rearrangements across cancer reveal classes of samples with distinct patterns of DNA breakage and rearrangement-induced hypermutability , 2012, Genome research.

[64]  Daniel E. Warren,et al.  The western painted turtle genome, a model for the evolution of extreme physiological adaptations in a slowly evolving lineage , 2013, Genome Biology.

[65]  J. Ausió,et al.  Vertebrate nucleoplasmin and NASP: egg histone storage proteins with multiple chaperone activities , 2012, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[66]  P. Marynen,et al.  Copy-number gains of HUWE1 due to replication- and recombination-based rearrangements. , 2012, American journal of human genetics.

[67]  Gabor T. Marth,et al.  Haplotype-based variant detection from short-read sequencing , 2012, 1207.3907.

[68]  S. Boissinot,et al.  Multi-Locus Phylogeographic and Population Genetic Analysis of Anolis carolinensis: Historical Demography of a Genomic Model Species , 2012, PloS one.

[69]  K. Fukui,et al.  RBMX: a regulator for maintenance and centromeric protection of sister chromatid cohesion. , 2012, Cell reports.

[70]  Rolf Backofen,et al.  Global or local? Predicting secondary structure and accessibility in mRNAs , 2012, Nucleic acids research.

[71]  M. Noor,et al.  Genomic impacts of chromosomal inversions in parapatric Drosophila species , 2012, Philosophical Transactions of the Royal Society B: Biological Sciences.

[72]  Martin Kircher,et al.  Double indexing overcomes inaccuracies in multiplex sequencing on the Illumina platform , 2011, Nucleic acids research.

[73]  P. Andolfatto,et al.  Methods to detect selection on noncoding DNA. , 2012, Methods in molecular biology.

[74]  O. Molnár,et al.  Ultraviolet nuptial colour determines fight success in male European green lizards (Lacerta viridis) , 2011, Biology Letters.

[75]  Peter F. Stadler,et al.  ViennaRNA Package 2.0 , 2011, Algorithms for Molecular Biology.

[76]  R. J. Harrison,et al.  A General Method for Calculating Likelihoods Under the Coalescent Process , 2011, Genetics.

[77]  Caitlin A. Kuczynski,et al.  Phylogeny of iguanian lizards inferred from 29 nuclear loci, and a comparison of concatenated and species-tree approaches for an ancient, rapid radiation. , 2011, Molecular phylogenetics and evolution.

[78]  O. Ishikawa,et al.  Changes in epidermal hyaluronan metabolism following UVB irradiation. , 2011, Journal of dermatological science.

[79]  F. Maquart,et al.  Ultraviolet‐B Irradiation Induces Differential Regulations of Hyaluronidase Expression and Activity in Normal Human Keratinocytes , 2011, Photochemistry and photobiology.

[80]  Jacob D. Jaffe,et al.  The genome of the green anole lizard and a comparative analysis with birds and mammals , 2011, Nature.

[81]  Katja Nowick,et al.  Gain, Loss and Divergence in Primate Zinc-Finger Genes: A Rich Resource for Evolution of Gene Regulatory Differences between Species , 2011, PloS one.

[82]  H. T. Ghashghaei,et al.  Specification of a Foxj1-Dependent Lineage in the Forebrain Is Required for Embryonic-to-Postnatal Transition of Neurogenesis in the Olfactory Bulb , 2011, The Journal of Neuroscience.

[83]  M. Gerstein,et al.  CNVnator: an approach to discover, genotype, and characterize typical and atypical CNVs from family and population genome sequencing. , 2011, Genome research.

[84]  Bradley P. Coe,et al.  Genome structural variation discovery and genotyping , 2011, Nature Reviews Genetics.

[85]  Tao Wang,et al.  miR-27 promotes osteoblast differentiation by modulating Wnt signaling. , 2010, Biochemical and biophysical research communications.

[86]  R. Faria,et al.  Chromosomal speciation revisited: rearranging theory with pieces of evidence. , 2010, Trends in ecology & evolution.

[87]  M. Kirkpatrick How and Why Chromosome Inversions Evolve , 2010, PLoS biology.

[88]  B. Sinervo,et al.  Lizards as model organisms for linking phylogeographic and speciation studies , 2010, Molecular ecology.

[89]  T. Houtani,et al.  GPR155: Gene organization, multiple mRNA splice variants and expression in mouse central nervous system. , 2010, Biochemical and biophysical research communications.

[90]  O. Molnár,et al.  Female European green lizards (Lacerta viridis) prefer males with high ultraviolet throat reflectance , 2010, Behavioral Ecology and Sociobiology.

[91]  C. Ponting,et al.  which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Molecular evolution of genes in avian genomes , 2010 .

[92]  Matthias Meyer,et al.  Illumina sequencing library preparation for highly multiplexed target capture and sequencing. , 2010, Cold Spring Harbor protocols.

[93]  Toralf Kirsten,et al.  Genomic organization of eukaryotic tRNAs , 2010, BMC Genomics.

[94]  T. Kiss,et al.  Box H/ACA small ribonucleoproteins. , 2010, Molecular cell.

[95]  D. Presgraves,et al.  The molecular evolutionary basis of species formation , 2010, Nature Reviews Genetics.

[96]  R. Wilson,et al.  BreakDancer: An algorithm for high resolution mapping of genomic structural variation , 2009, Nature Methods.

[97]  J. Mattick,et al.  Small RNAs derived from snoRNAs. , 2009, RNA.

[98]  Kai Ye,et al.  Pindel: a pattern growth approach to detect break points of large deletions and medium sized insertions from paired-end short reads , 2009, Bioinform..

[99]  B. Charlesworth Effective population size and patterns of molecular evolution and variation , 2009, Nature Reviews Genetics.

[100]  R. Emerson,et al.  Adaptive Evolution in Zinc Finger Transcription Factors , 2009, PLoS genetics.

[101]  Chao Xie,et al.  CNV-seq, a new method to detect copy number variation using high-throughput sequencing , 2009, BMC Bioinformatics.

[102]  G. Odierna,et al.  Karyotype, chromosome structure, reproductive modalities of three Southern Eurasian populations of the common lacertid lizard, Zootoca vivipara (Jacquin, 1787) , 2008 .

[103]  G. Meister,et al.  Tdrd3 is a novel stress granule-associated protein interacting with the Fragile-X syndrome protein FMRP. , 2008, Human molecular genetics.

[104]  Y. Surget-Groba,et al.  The relative importance of ecology and geographic isolation for speciation in anoles , 2008, Philosophical Transactions of the Royal Society B: Biological Sciences.

[105]  P. Marynen,et al.  Nonrecurrent MECP2 duplications mediated by genomic architecture-driven DNA breaks and break-induced replication repair. , 2008, Genome research.

[106]  R. Nusse,et al.  Wnt signaling and stem cell control , 2008, Cell Research.

[107]  Z. Nagy,et al.  Phylogeography of western Palaearctic reptiles - Spatial and temporal speciation patterns $ , 2007 .

[108]  JaneR . Taylor,et al.  Inhibition of Rho via Arg and p190RhoGAP in the Postnatal Mouse Hippocampus Regulates Dendritic Spine Maturation, Synapse and Dendrite Stability, and Behavior , 2007, The Journal of Neuroscience.

[109]  R. Ricklefs,et al.  Evolutionary diversification of clades of squamate reptiles , 2007, Journal of evolutionary biology.

[110]  D. Presgraves,et al.  High-Resolution Genome-Wide Dissection of the Two Rules of Speciation in Drosophila , 2007, PLoS biology.

[111]  Ziheng Yang PAML 4: phylogenetic analysis by maximum likelihood. , 2007, Molecular biology and evolution.

[112]  N. Kane,et al.  Meeting review: American Genetics Association Symposium on the genetics of speciation , 2007, Molecular ecology.

[113]  M. Böhme,et al.  Phylogeography and cryptic variation within the Lacerta viridis complex (Lacertidae, Reptilia) , 2007 .

[114]  J. Sleeman,et al.  Differential regulation of hyaluronan metabolism in the epidermal and dermal compartments of human skin by UVB irradiation. , 2007, The Journal of investigative dermatology.

[115]  M. A. U. B. Öhme Phylogeography and cryptic variation within the Lacerta viridis complex ( Lacertidae , Reptilia ) , 2007 .

[116]  Robert S. Harris,et al.  Improved pairwise alignment of genomic dna , 2007 .

[117]  M. Böhme,et al.  Small edge populations at risk: genetic diversity of the green lizard (Lacerta viridis viridis) in Germany and implications for conservation management , 2007, Conservation Genetics.

[118]  H. Vaucheret,et al.  Arabidopsis RPA2: A Genetic Link among Transcriptional Gene Silencing, DNA Repair, and DNA Replication , 2005, Current Biology.

[119]  E. Olmo Rate of Chromosome changes and Speciation in Reptiles , 2005, Genetica.

[120]  R. Bleiweiss Ultraviolet plumage reflectance distinguishes sibling bird species. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[121]  Heiko Schoof,et al.  Conservation, diversification and expansion of C2H2 zinc finger proteins in the Arabidopsis thaliana genome , 2004, BMC Genomics.

[122]  D. Haussler,et al.  Aligning multiple genomic sequences with the threaded blockset aligner. , 2004, Genome research.

[123]  Chau-Ti Ting,et al.  Genes and speciation , 2001, Nature Reviews Genetics.

[124]  E. Olmo,et al.  C-band variability and phylogeny of Lacertidae , 1986, Genetica.

[125]  Raul Urrutia,et al.  KRAB-containing zinc-finger repressor proteins , 2003, Genome Biology.

[126]  D. Haussler,et al.  Evolution's cauldron: Duplication, deletion, and rearrangement in the mouse and human genomes , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[127]  N. Barton,et al.  Chromosomal Speciation and Molecular Divergence--Accelerated Evolution in Rearranged Chromosomes , 2003, Science.

[128]  N. Barton,et al.  ACCUMULATING POSTZYGOTIC ISOLATION GENES IN PARAPATRY: A NEW TWIST ON CHROMOSOMAL SPECIATION , 2003, Evolution; international journal of organic evolution.

[129]  Daniel J. Gaffney,et al.  Quantifying the slightly deleterious mutation model of molecular evolution. , 2002, Molecular biology and evolution.

[130]  W. Paschen,et al.  Mechanisms Underlying Suppression of Protein Synthesis Induced by Transient Focal Cerebral Ischemia in Mouse Brain , 2002, Experimental Neurology.

[131]  E. Mayr Wu’s genic view of speciation , 2001 .

[132]  L. Rieseberg,et al.  A genic view of species integration , 2001 .

[133]  M. Noor,et al.  Chromosomal inversions and the reproductive isolation of species , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[134]  L H. Rieseberg,et al.  Chromosomal rearrangements and speciation. , 2001, Trends in ecology & evolution.

[135]  C. Pabo,et al.  DNA recognition by Cys2His2 zinc finger proteins. , 2000, Annual review of biophysics and biomolecular structure.

[136]  N. Tvrtković,et al.  THE WESTERN GREEN LIZARD Lacerta (viridis) bilineata Daudin, 1804 (Sauria: Lacertidae) IN SLOVENIA AND CROATIA , 1998 .

[137]  S. Eddy,et al.  tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. , 1997, Nucleic acids research.

[138]  H. A. Orr,et al.  Dobzhansky, Bateson, and the genetics of speciation. , 1996, Genetics.

[139]  S. Rykena Kreuzungsexperimente zur Prüfung der Artgrenzen im Genus Lacerta sensu stricto , 1991 .

[140]  J. Baron,et al.  Démographie comparée de deux populations françaises de Lacerta viridis (Laurenti, 1768) , 1989, Revue d'Écologie (La Terre et La Vie).

[141]  Nicholas H. Barton,et al.  The Relative Rates of Evolution of Sex Chromosomes and Autosomes , 1987, The American Naturalist.

[142]  M. Kimura,et al.  Genetic variability and effective population size when local extinction and recolonization of subpopulations are frequent. , 1980, Proceedings of the National Academy of Sciences of the United States of America.

[143]  G. Caughley Parameters for Seasonally Breeding Populations. , 1967, Ecology.