SlideShare a Scribd company logo
Santiago Montero-Mendieta
Qualifications: Biologist
Specializations: Genomics, systematics, …
Current position: PhD student
2009-2013: Degree in Biology (University of Girona, Spain)
2013-2014: MSc in Biodiversity, focusing on Evolutionary
Biology (University of Barcelona, Spain)
2015-Present: PhD Student, Estación Biológica de
Doñana (CSIC) Seville, Spain (1 year and a half)
# ABOUT ME
1/25
A GENOMIC VIEW ON THE
DIVERSIFICATION OF NEOTROPICAL FROGS
(provisional title)
Main advisor: Carles Vilà
Collaborators: Jennifer Leonard, Matthew Webster,
José Manuel Padial & Ignacio De la Riva
The history of life
2/25
The theory of evolution is based on the idea that all species are
related and gradually change over time.
Speciation
Mutation
Recombination
Natural selection
Genetic driftHybridization
Polyploidy
Isolation
What is speciation?
3/25
The formation of new and distinct species in the course of evolution
Model & non-model organisms
4/25
High-throughput sequencing (e.g. Illumina) makes non-model
organisms increasingly accessible for speciation studies, mainly
through proteomics
Speciation genomics
5/25
LIMITED
POWER
# Pre-selected
markers are used
# Need of sampling
multiple individuals
Large
amounts of
orthologous
loci can be
obtained,
allowing the
use of less
individuals
1. It allows finding genes
involved in speciation
2. It allows finding genes
homogenized by gene-flow (or
those that resist introgression)
3. It allows finding genes related
to adaptation
Humans & Neanderthals mated in the past
Green et al. (2010)
Neotropical amphibians
~ 50% of world’s amphibians
6/25
Our study model
 Frogs of the genus Oreobates
24 species up to date
Distributed
across a wide
range of
habitats and
altitudes
7/25
Cloud forests
(2500-3500m)
Montane forests
(500-1500m)
Lowland
forests
(0-500m)
Dry forests (500-1500m)
Still little is known
Kohler & Padial
(unpublished data) 8/25
?47 Million Years Ago
46.7 Mya
The best non-model
Extremely difficult to sample
– Difficult to find
– Few museums have specimens
– Logistic problems (permits)
Few genomic data available in frogs
– Xenopus tropicalis (206.6 Mya)
– Nanorana parkeri (156.0 Mya)
Oreobates zongoensis
Some species have been only
found once (by our collaborators)
>>> Challenge effect <<<
9/25
Why did we chose Oreobates?
10/25
 Oreobates amarakaeri (Padial et al., 2012)
 Oreobates ayacucho (Lehr, 2007)
 Oreobates barituensis (Vaira & Ferrari, 2008)
 Oreobates berdemenos (Pereyra et al., 2014)
 Oreobates choristolemma (Harvey & Sheehy, 2005)
 Oreobates crepitans (Bokermann, 1965)
 Oreobates cruralis (Boulenger, 1902)
 Oreobates discoidalis (Peracca, 1895)
 Oreobates gemcare (Padial et al., 2012)
 Oreobates granulosus (Boulenger, 1902)
 Oreobates heterodactylus (Miranda-Ribeiro, 1937)
 Oreobates ibischi (Reichle, et al. 2001)
 Oreobates lehri (Padial et al., 2007)
 Oreobates lundbergi (Lehr, 2005)
 Oreobates machiguenga (Padial et al., 2012)
 Oreobates madidi (Padial et al., 2005)
 Oreobates pereger (Lynch, 1975)
 Oreobates quixensis (Jiménez de la Espada, 1872)
 Oreobates remotus (Teixeira et al., 2012)
 Oreobates sanctaecrucis (Harvey & Keck, 1995)
 Oreobates sanderi (Padial, et al., 2005)
 Oreobates saxatilis (Duellman, 1990)
 Oreobates yanucu (Kohler & Padial 2016)
 Oreobates zongoensis (Reichle & Köhler, 1997)
We have access to
(almost) all the
Oreobates species
José Manuel
Padial
Ignacio De la Riva
Our collaborators
are experts on
these frogs
Research goals
1. Phylogenomics: genetic relationship among Oreobates
2. Evolutionary history: study variation in evolution rates
3. Demographic history: track demographic changes
through time and correspondence with habitat changes
4. Adaptation: identify genes that have been
differentiated between populations (adaptation)
To study evolution rates, demographic history and
adaptation patterns on the frogs of the genus
Oreobates
11/25
1st stage : phylogenomics
12/25
Wet puna
Amazonian lowland forest
Lower montane forest +
dry forest
Upper montane forest
Cloud forest
4000 m
3000 m
2000 m
1000 m
0
Andes mountain diagram
How many colonization
events to lowland rainforest
have occurred?
H-1: Oreobates
emerged on the
Andes highlands
Fact: Now there are
species living on both
highlands and lowlands
When the genus Oreobates
was originated? How many
Oreobates species?
Goal: to build a highly supported
tree for downstream analysis
2nd stage: evolutionary history
13/25
Is the evolution rate lower in
the highland species?
Goal: study the variation in the
evolution rate of Oreobates
H-2: Ectotherm
metabolism slows
down at low
temperatures
Fact: Previous studies in glass-
frogs proved a reduction in the
rate of evolution in highland
environments
Wet puna
Amazonian lowland forest
Lower montane forest +
dry forest
Upper montane forest
Cloud forest
4000 m
3000 m
2000 m
1000 m
0
Andes mountain diagram
3rd stage: demographic history
14/25
Goal: study the effect of the
past environmental conditions on
the Oreobates demography
H-3: Species with similar habitat
requirements will show parallel
demographic changes during Pleistocene
climate changes
Do highland species show
different demographic
trends compared to lowland?
Is there any hybridization
between diverging linages
living on the lowlands?
?
O. quixensis
O. saxatilis
Wet puna
Amazonian lowland forest
Lower montane forest +
dry forest
Upper montane forest
Cloud forest
4000 m
3000 m
2000 m
1000 m
0
Andes mountain diagram
4th stage: study of adaptation
15/25
Are there any genes related
to adaptation to dry forest?
Fact: Some
Oreobates inhabit
a wide diversity
of habitats
Goal: study the genomic
signatures of speciation
H-4: Genomic regions associated
with adaptation to environment
should show larger genetic
divergences
Are there any genes related
to adaptation to highlands?
O. cruralis
BARRIER
Wet puna
Amazonian lowland forest
Lower montane forest +
dry forest
Upper montane forest
Cloud forest
4000 m
3000 m
2000 m
1000 m
0
Andes mountain diagram
Our initial idea
16/25
 1st : Transcriptome
sequencing (as a reference)
 2nd : Whole genome
sequencing (for the others)
 3rd : Exome assembly
 4th : SNP detection and
analyzing data
DRAWBACK: big waste
Adapted from Lamichhaney et al. (2012)
NEED OF: genome size
Amphibians have big genomes
17/25
# The C-value Enigma
Previous
work
Reduced representation of the
genome: transcriptome
How to do it?
Total
RNA
Transcriptome sequencing
Bioinformatics
analysis Liver
 Intestine
 Spleen
De-novo
transcriptome
assembly
Uppsala University
Oreobates
cruralis
19/25
Transcriptome-
based exon
capture
(on-going)
Phylogenomics & more NOW
Transcriptome workflow
20/25
# Quick guide to build de-novo assemblies
1. Get raw reads (RNAseq data)
2. Quality control [FastQC]
3. Ribosomal RNA removal [SortMeRNA]
4. Quality trimming & adapter removal
[Trimmomatic]
5. Quality control (again) [FastQC]
6. In silico normalization [Trinity]
7. Merge data (when multiple tissues per sample)
8. In silico normalization (again) [Trinity]
9. De-novo transcriptome assembly [Trinity]
9.1. Assembly validation [Bowtie]
10. Functional annotation [Trinotate]
Transcriptome results (I)
523M reads
520M reads
426M reads
36M reads
Raw data (intestine, liver, spleen)
After ribosomal RNA removal
After quality trimming & adapter removal
After normalization
550,871 transcripts De-novo Transcriptome Assembly
45,885 unigenes Unigenes with a Gene Ontology (GO) match
40,362 unigenes Unigenes that have 40-60% GC content
21,369 unigenes Unigenes that are & 500-850bp in length
17,879 unigenes Removing repeats and low complexity regions
Pre-
processing
21/25
422,999 unigenes Selecting longest isoform for each trinity ‘gene’
Post-
processing
Exon capture
0 10000 20000 30000 40000 50000 60000
Xenopus_BLASTX
SwissProt_BLASTX
GO_SwissProt_Pfam
GO_SwisProt
Kegg
eggNOG
SwissProt_BLASTP
Pfam
GO_Pfam
# Annotated unigenes
Transcriptome results (II)
22/25
10,463
orthologous
genes
O.cruralis
X.tropicalis
Transcriptome results (III)
23/25
Cellular process........................ 78%
Single-organism process........ 64%
Metabolic process………………….63%
Binding………………………………………70%
Catalytic activity……………………39%
Cell…………………………………………….82%
Cell part……………………………………81%
Organelle………………………………... 67%
# Gene Ontology (GO) categories
Biological process
Molecular function
Cellular component
Percentage of
unigenes
0 5000 10000 15000 20000 25000
Homo sapiens
Mus musculus
Arabidopsis thaliana
Rattus norvegicus
Saccharomyces cerevisiae
Bos taurus
Schizosaccharomyces pombe
Escherichia coli
Bacillus subtilis
Dictyostelium discoideum
Caenorhabditis elegans
Oryza sativa
Xenopus laevis
Drosophila melanogaster
Danio rerio
Gallus gallus
Pongo abelii
Mycobacterium tuberculosis
Escherichia coli
Mycobacterium tuberculosis
Top-20 most represented species in
UniprotKB/SwissProt
Transcriptome results (IV)
24/25
Homo sapiens
Mus musculus
Caenorhabditis
elegans
Xenopus laevis
Rattus norvegicus
Bos taurus
Xenopus tropicalis
Caenorhabditis briggsae
Gallus gallus
Others
TOP BLASTX-HIT SPECIES DISTRIBUTION
IN OREOBATES CURALIS
19%
11%
10%
6%
5%
5%
4%
3% 5%
32% 32%
# Transcriptome-based exon capture
So, what is next?
We are using the
17,879 unigene
sequences from
O. cruralis to design
capture probes for all
other Oreobates
species.
“ We gotta catch ‘em all ! ” 25/25
Orthologous genes
will be identified
and used to test
initial hypothesis
Thanks for your attention!
QUESTIONS?

More Related Content

PPTX
Job Talk Iowa State University Ag Bio Engineering
PDF
Open Tree of Life Phyloseminar 2014
PDF
ISU ENVSCI690 Graduate Seminar Slides
PDF
J. Mollus. Stud.-2015-Carvalho-mollus-eyv014
PPTX
Global patterns of insect diiversity, distribution and evolutionary distinctness
PDF
BM405 Lecture Slides 21/11/2014 University of Strathclyde
PPTX
Rapid Impact Assessment of Climatic and Physio-graphic Changes on Flagship G...
PDF
Towards identifying novel phenotypes in climate adapted livestock production
Job Talk Iowa State University Ag Bio Engineering
Open Tree of Life Phyloseminar 2014
ISU ENVSCI690 Graduate Seminar Slides
J. Mollus. Stud.-2015-Carvalho-mollus-eyv014
Global patterns of insect diiversity, distribution and evolutionary distinctness
BM405 Lecture Slides 21/11/2014 University of Strathclyde
Rapid Impact Assessment of Climatic and Physio-graphic Changes on Flagship G...
Towards identifying novel phenotypes in climate adapted livestock production

Similar to A genomic view on the diversification of Neotropical frogs (20)

PDF
2017 - Environmental Ordination of Filamentous Bacteria in Activated Sludge
PDF
Morphology, molecular, and nutritional value of Amphora sp. from coastal wate...
PDF
2016 10-27 timbers
PDF
So many different kinds of mistakes: Or why systematic error is the 21st cent...
PDF
Diversity Diversity Diversity Diversity ....
PPTX
Post genomic microbiology rodriguez valera
PDF
Identification of fish species using dna barcode from visakhapatnam, east coa...
PDF
Ecological Monitoring Techniques
PPTX
CCBC tutorial beiko
PPTX
Big data nebraska
PPTX
Big data nebraska
PPTX
Ecosytem services for food and nutritional security
PDF
5 mohammad chamani
PDF
Microbial Agrogenomics 4/2/2015, UK-MX Workshop
PDF
Determinació per factors ambientals de la proporció de sexes a les poblacions...
PDF
Paper to Upload, MOLECULAR PHYLOGENY OF CATFISHES.pdf
PPT
13128931.ppt
PPTX
Open Data and Ecological and Evolutionary synthesis
PDF
Maquia et al., 2013
2017 - Environmental Ordination of Filamentous Bacteria in Activated Sludge
Morphology, molecular, and nutritional value of Amphora sp. from coastal wate...
2016 10-27 timbers
So many different kinds of mistakes: Or why systematic error is the 21st cent...
Diversity Diversity Diversity Diversity ....
Post genomic microbiology rodriguez valera
Identification of fish species using dna barcode from visakhapatnam, east coa...
Ecological Monitoring Techniques
CCBC tutorial beiko
Big data nebraska
Big data nebraska
Ecosytem services for food and nutritional security
5 mohammad chamani
Microbial Agrogenomics 4/2/2015, UK-MX Workshop
Determinació per factors ambientals de la proporció de sexes a les poblacions...
Paper to Upload, MOLECULAR PHYLOGENY OF CATFISHES.pdf
13128931.ppt
Open Data and Ecological and Evolutionary synthesis
Maquia et al., 2013
Ad

Recently uploaded (20)

PPTX
Classification Systems_TAXONOMY_SCIENCE8.pptx
PPTX
Introduction to Fisheries Biotechnology_Lesson 1.pptx
PPTX
neck nodes and dissection types and lymph nodes levels
PPTX
TOTAL hIP ARTHROPLASTY Presentation.pptx
PDF
SEHH2274 Organic Chemistry Notes 1 Structure and Bonding.pdf
PPTX
Protein & Amino Acid Structures Levels of protein structure (primary, seconda...
PPTX
ognitive-behavioral therapy, mindfulness-based approaches, coping skills trai...
PPTX
The KM-GBF monitoring framework – status & key messages.pptx
PPTX
Vitamins & Minerals: Complete Guide to Functions, Food Sources, Deficiency Si...
PPTX
Comparative Structure of Integument in Vertebrates.pptx
PPTX
Cell Membrane: Structure, Composition & Functions
PPTX
microscope-Lecturecjchchchchcuvuvhc.pptx
PPTX
GEN. BIO 1 - CELL TYPES & CELL MODIFICATIONS
PDF
CAPERS-LRD-z9:AGas-enshroudedLittleRedDotHostingaBroad-lineActive GalacticNuc...
PPTX
2. Earth - The Living Planet earth and life
PPTX
G5Q1W8 PPT SCIENCE.pptx 2025-2026 GRADE 5
PDF
VARICELLA VACCINATION: A POTENTIAL STRATEGY FOR PREVENTING MULTIPLE SCLEROSIS
PDF
HPLC-PPT.docx high performance liquid chromatography
PDF
Sciences of Europe No 170 (2025)
PPTX
EPIDURAL ANESTHESIA ANATOMY AND PHYSIOLOGY.pptx
Classification Systems_TAXONOMY_SCIENCE8.pptx
Introduction to Fisheries Biotechnology_Lesson 1.pptx
neck nodes and dissection types and lymph nodes levels
TOTAL hIP ARTHROPLASTY Presentation.pptx
SEHH2274 Organic Chemistry Notes 1 Structure and Bonding.pdf
Protein & Amino Acid Structures Levels of protein structure (primary, seconda...
ognitive-behavioral therapy, mindfulness-based approaches, coping skills trai...
The KM-GBF monitoring framework – status & key messages.pptx
Vitamins & Minerals: Complete Guide to Functions, Food Sources, Deficiency Si...
Comparative Structure of Integument in Vertebrates.pptx
Cell Membrane: Structure, Composition & Functions
microscope-Lecturecjchchchchcuvuvhc.pptx
GEN. BIO 1 - CELL TYPES & CELL MODIFICATIONS
CAPERS-LRD-z9:AGas-enshroudedLittleRedDotHostingaBroad-lineActive GalacticNuc...
2. Earth - The Living Planet earth and life
G5Q1W8 PPT SCIENCE.pptx 2025-2026 GRADE 5
VARICELLA VACCINATION: A POTENTIAL STRATEGY FOR PREVENTING MULTIPLE SCLEROSIS
HPLC-PPT.docx high performance liquid chromatography
Sciences of Europe No 170 (2025)
EPIDURAL ANESTHESIA ANATOMY AND PHYSIOLOGY.pptx
Ad

A genomic view on the diversification of Neotropical frogs

  • 1. Santiago Montero-Mendieta Qualifications: Biologist Specializations: Genomics, systematics, … Current position: PhD student 2009-2013: Degree in Biology (University of Girona, Spain) 2013-2014: MSc in Biodiversity, focusing on Evolutionary Biology (University of Barcelona, Spain) 2015-Present: PhD Student, Estación Biológica de Doñana (CSIC) Seville, Spain (1 year and a half) # ABOUT ME
  • 2. 1/25 A GENOMIC VIEW ON THE DIVERSIFICATION OF NEOTROPICAL FROGS (provisional title) Main advisor: Carles Vilà Collaborators: Jennifer Leonard, Matthew Webster, José Manuel Padial & Ignacio De la Riva
  • 3. The history of life 2/25 The theory of evolution is based on the idea that all species are related and gradually change over time.
  • 4. Speciation Mutation Recombination Natural selection Genetic driftHybridization Polyploidy Isolation What is speciation? 3/25 The formation of new and distinct species in the course of evolution
  • 5. Model & non-model organisms 4/25 High-throughput sequencing (e.g. Illumina) makes non-model organisms increasingly accessible for speciation studies, mainly through proteomics
  • 6. Speciation genomics 5/25 LIMITED POWER # Pre-selected markers are used # Need of sampling multiple individuals Large amounts of orthologous loci can be obtained, allowing the use of less individuals 1. It allows finding genes involved in speciation 2. It allows finding genes homogenized by gene-flow (or those that resist introgression) 3. It allows finding genes related to adaptation Humans & Neanderthals mated in the past Green et al. (2010)
  • 7. Neotropical amphibians ~ 50% of world’s amphibians 6/25
  • 8. Our study model  Frogs of the genus Oreobates 24 species up to date Distributed across a wide range of habitats and altitudes 7/25 Cloud forests (2500-3500m) Montane forests (500-1500m) Lowland forests (0-500m) Dry forests (500-1500m)
  • 9. Still little is known Kohler & Padial (unpublished data) 8/25 ?47 Million Years Ago 46.7 Mya
  • 10. The best non-model Extremely difficult to sample – Difficult to find – Few museums have specimens – Logistic problems (permits) Few genomic data available in frogs – Xenopus tropicalis (206.6 Mya) – Nanorana parkeri (156.0 Mya) Oreobates zongoensis Some species have been only found once (by our collaborators) >>> Challenge effect <<< 9/25
  • 11. Why did we chose Oreobates? 10/25  Oreobates amarakaeri (Padial et al., 2012)  Oreobates ayacucho (Lehr, 2007)  Oreobates barituensis (Vaira & Ferrari, 2008)  Oreobates berdemenos (Pereyra et al., 2014)  Oreobates choristolemma (Harvey & Sheehy, 2005)  Oreobates crepitans (Bokermann, 1965)  Oreobates cruralis (Boulenger, 1902)  Oreobates discoidalis (Peracca, 1895)  Oreobates gemcare (Padial et al., 2012)  Oreobates granulosus (Boulenger, 1902)  Oreobates heterodactylus (Miranda-Ribeiro, 1937)  Oreobates ibischi (Reichle, et al. 2001)  Oreobates lehri (Padial et al., 2007)  Oreobates lundbergi (Lehr, 2005)  Oreobates machiguenga (Padial et al., 2012)  Oreobates madidi (Padial et al., 2005)  Oreobates pereger (Lynch, 1975)  Oreobates quixensis (Jiménez de la Espada, 1872)  Oreobates remotus (Teixeira et al., 2012)  Oreobates sanctaecrucis (Harvey & Keck, 1995)  Oreobates sanderi (Padial, et al., 2005)  Oreobates saxatilis (Duellman, 1990)  Oreobates yanucu (Kohler & Padial 2016)  Oreobates zongoensis (Reichle & Köhler, 1997) We have access to (almost) all the Oreobates species José Manuel Padial Ignacio De la Riva Our collaborators are experts on these frogs
  • 12. Research goals 1. Phylogenomics: genetic relationship among Oreobates 2. Evolutionary history: study variation in evolution rates 3. Demographic history: track demographic changes through time and correspondence with habitat changes 4. Adaptation: identify genes that have been differentiated between populations (adaptation) To study evolution rates, demographic history and adaptation patterns on the frogs of the genus Oreobates 11/25
  • 13. 1st stage : phylogenomics 12/25 Wet puna Amazonian lowland forest Lower montane forest + dry forest Upper montane forest Cloud forest 4000 m 3000 m 2000 m 1000 m 0 Andes mountain diagram How many colonization events to lowland rainforest have occurred? H-1: Oreobates emerged on the Andes highlands Fact: Now there are species living on both highlands and lowlands When the genus Oreobates was originated? How many Oreobates species? Goal: to build a highly supported tree for downstream analysis
  • 14. 2nd stage: evolutionary history 13/25 Is the evolution rate lower in the highland species? Goal: study the variation in the evolution rate of Oreobates H-2: Ectotherm metabolism slows down at low temperatures Fact: Previous studies in glass- frogs proved a reduction in the rate of evolution in highland environments Wet puna Amazonian lowland forest Lower montane forest + dry forest Upper montane forest Cloud forest 4000 m 3000 m 2000 m 1000 m 0 Andes mountain diagram
  • 15. 3rd stage: demographic history 14/25 Goal: study the effect of the past environmental conditions on the Oreobates demography H-3: Species with similar habitat requirements will show parallel demographic changes during Pleistocene climate changes Do highland species show different demographic trends compared to lowland? Is there any hybridization between diverging linages living on the lowlands? ? O. quixensis O. saxatilis Wet puna Amazonian lowland forest Lower montane forest + dry forest Upper montane forest Cloud forest 4000 m 3000 m 2000 m 1000 m 0 Andes mountain diagram
  • 16. 4th stage: study of adaptation 15/25 Are there any genes related to adaptation to dry forest? Fact: Some Oreobates inhabit a wide diversity of habitats Goal: study the genomic signatures of speciation H-4: Genomic regions associated with adaptation to environment should show larger genetic divergences Are there any genes related to adaptation to highlands? O. cruralis BARRIER Wet puna Amazonian lowland forest Lower montane forest + dry forest Upper montane forest Cloud forest 4000 m 3000 m 2000 m 1000 m 0 Andes mountain diagram
  • 17. Our initial idea 16/25  1st : Transcriptome sequencing (as a reference)  2nd : Whole genome sequencing (for the others)  3rd : Exome assembly  4th : SNP detection and analyzing data DRAWBACK: big waste Adapted from Lamichhaney et al. (2012) NEED OF: genome size
  • 18. Amphibians have big genomes 17/25 # The C-value Enigma Previous work
  • 19. Reduced representation of the genome: transcriptome
  • 20. How to do it? Total RNA Transcriptome sequencing Bioinformatics analysis Liver  Intestine  Spleen De-novo transcriptome assembly Uppsala University Oreobates cruralis 19/25 Transcriptome- based exon capture (on-going) Phylogenomics & more NOW
  • 21. Transcriptome workflow 20/25 # Quick guide to build de-novo assemblies 1. Get raw reads (RNAseq data) 2. Quality control [FastQC] 3. Ribosomal RNA removal [SortMeRNA] 4. Quality trimming & adapter removal [Trimmomatic] 5. Quality control (again) [FastQC] 6. In silico normalization [Trinity] 7. Merge data (when multiple tissues per sample) 8. In silico normalization (again) [Trinity] 9. De-novo transcriptome assembly [Trinity] 9.1. Assembly validation [Bowtie] 10. Functional annotation [Trinotate]
  • 22. Transcriptome results (I) 523M reads 520M reads 426M reads 36M reads Raw data (intestine, liver, spleen) After ribosomal RNA removal After quality trimming & adapter removal After normalization 550,871 transcripts De-novo Transcriptome Assembly 45,885 unigenes Unigenes with a Gene Ontology (GO) match 40,362 unigenes Unigenes that have 40-60% GC content 21,369 unigenes Unigenes that are & 500-850bp in length 17,879 unigenes Removing repeats and low complexity regions Pre- processing 21/25 422,999 unigenes Selecting longest isoform for each trinity ‘gene’ Post- processing Exon capture
  • 23. 0 10000 20000 30000 40000 50000 60000 Xenopus_BLASTX SwissProt_BLASTX GO_SwissProt_Pfam GO_SwisProt Kegg eggNOG SwissProt_BLASTP Pfam GO_Pfam # Annotated unigenes Transcriptome results (II) 22/25 10,463 orthologous genes O.cruralis X.tropicalis
  • 24. Transcriptome results (III) 23/25 Cellular process........................ 78% Single-organism process........ 64% Metabolic process………………….63% Binding………………………………………70% Catalytic activity……………………39% Cell…………………………………………….82% Cell part……………………………………81% Organelle………………………………... 67% # Gene Ontology (GO) categories Biological process Molecular function Cellular component Percentage of unigenes
  • 25. 0 5000 10000 15000 20000 25000 Homo sapiens Mus musculus Arabidopsis thaliana Rattus norvegicus Saccharomyces cerevisiae Bos taurus Schizosaccharomyces pombe Escherichia coli Bacillus subtilis Dictyostelium discoideum Caenorhabditis elegans Oryza sativa Xenopus laevis Drosophila melanogaster Danio rerio Gallus gallus Pongo abelii Mycobacterium tuberculosis Escherichia coli Mycobacterium tuberculosis Top-20 most represented species in UniprotKB/SwissProt Transcriptome results (IV) 24/25 Homo sapiens Mus musculus Caenorhabditis elegans Xenopus laevis Rattus norvegicus Bos taurus Xenopus tropicalis Caenorhabditis briggsae Gallus gallus Others TOP BLASTX-HIT SPECIES DISTRIBUTION IN OREOBATES CURALIS 19% 11% 10% 6% 5% 5% 4% 3% 5% 32% 32%
  • 26. # Transcriptome-based exon capture So, what is next? We are using the 17,879 unigene sequences from O. cruralis to design capture probes for all other Oreobates species. “ We gotta catch ‘em all ! ” 25/25 Orthologous genes will be identified and used to test initial hypothesis
  • 27. Thanks for your attention! QUESTIONS?

Editor's Notes

  • #3: A genomic view on the diversification of neotropical frogs
  • #4: The theory of evolution is based on the idea that all species are related and gradually change over the time. The accumulation of small changes over the time can produce that some populations of the same species get enough differences that they will be “reproductively isolated” and so they will no be able to mate between them anymore. We call this “speciation event”.
  • #5: What is speciation? The formation of new and distinct species in the course of evolution. There are several ways in which a species can evolve to become two (or more) distinct species. These ways are: isolation, mutation, recombination, natural selection, genetic drift, hybridization and polyploidy.
  • #6: The study of speciation by using genetic resources has been mainly applied on model organisms, such as human & mouse. However, nowadays we can use high-throughput sequencing techniques (such as illumina) to study the speciation in non-model organisms too. We call this “genomics speciation”.
  • #7: Before, with Sanger sequencing (and other similar approaches) we had to pre-select some markers and we had to “believe” that this markers were showing us “the truth”. For this reason it was necessary to use multiple individuals from several populations in order to know if ”the truth that we saw” was actually “the right truth”. Nowadays, by using high-throughput sequencing techniques, we can obtain thousands of orthologous loci, and so this allow us to use less individuals. So, “speciation genomics” mainly refers to: 1) finding genes involved in speciation, 2) finding genes that are homogenized by gene-flow, and 3) finding genes that are related to adaptation processes. But genomics is not only about finding new genes, we can also use it to reconstruct evolutionary and demographic histories.
  • #8: Neotropical amphibians are the 50% of amphibians in the world. Many new species are being described every year. At the same time that many others become extinct too. However, not many studies have applied genomics to study the speciation in Neotropical frogs. This means that we still know little about the patterns and processes that produced this huge diversity on the tropical regions.
  • #9: Our study model is a genus of Neotropical frogs named Oreobates. To the present date, this genus includes 24 species that are distributed across a wide range of habitats and altitudes. All of them are found in South America (as you can see in the map). It’s quite interesting that Oreobates species can be mainly found up to four types of habitat: cloud forests, montane forests, lowland forests, and even dry forests.
  • #10: Oreobates is a very interesting genus of frogs not only because they live in a wide variety of habitats but also because little is know about them. More than half of the 24-named species have been described in the last ten years. And, the phylogenetic relationships among them are not clear. Some ”trees” have been published (using few molecular markers) but still the evolutionary history is not clear either (as you can see).
  • #11: We think that Oreobates are one of the best non-model species in the world. They are extremely difficult to find. Also, there are very few museums that have Oreobates samples. And if you want to get samples from them, there is always problems with permits. Moreover, there are some species that have been only found once. Luckily our collaborators were the ones collecting it. And that’s probably another big reason why we choose Oreobates.
  • #12: Our collaborators (José Manuel Padial and Ignacio de la Riva, from Madrid) they probably are the ones who know the most about these frogs. As you can see “in blue” most of the Oreobates species have been described by them. And “in red” you can see which species are available in the MNCN (Madrid). As you can see, we have access to almost all the species. And this quite good, because it means that nobody can steal us this work.
  • #13: So, the main goal of my PhD project is to study the evolution rates, demographic history and adaptation patterns on the frogs of the genus Oreobates (using genomic data). As you can see, we have four different goals (phylogenomics, evolutionary history, demographic history and adaptation). I’m going to explain them in more detail in the following slides.
  • #14: First stage: phylogenomics. Our goal is to reconstruct a highly supported tree among Oreobates species (using genomic data). This will give us a very good idea about the relationships among Oreobates. We will use this tree not only for downstream analysis but also to address questions related to the phylogeny. We know that there are Oreobates species living on both highlands and lowlands. It has been hypothesized that Oreobates emerged on the Andes highlands. This is quite weird, since usually it’s the contrary (species first emerge on the lowlands and the colonize the highlands). We want to know WHEN the genus Oreobates was originated and HOW many colonization events to lowland rainforest have occurred. We will use a time-calibrated phylogenomics tree for looking at this.
  • #15: Second stage: evolutionary history. Here our goal is to study the variation in the evolution rate of the different Oreobates species. Previous studies (by Alvaro Dugo-Cota) have shown that glass-frogs living in the highlands are subjected to a reduction in the rate of evolution. In this study, the people from our group studied over 100 species of frogs, but using few molecular markers. The results were crystal clear specially for mitochondrial markers but not so clear for nuclear. They hypothesized that is because ectotherm metabolism is reduced at low temperatures. In our case, using thousands of molecular markers, we want to test if Oroebates species from the highlands also show a lower rate of evolution compared to the species living on the lowlands.
  • #16: Third stage: demographic history. Our goal is to study the effect of the past environmental conditions on the Oreobates demography. The Andes formation was about 20 to 30 MYA and so the origin of Oreobates too. As you may know, about 2.5 MYA (on the Pleistocene) the Earth faced repeated glacial cycles. And so species living in the highlands are expected to suffer populations increasing and reductions following climatic changes. We hypothesized that species with similar habitat requirements will show parallel demographic changes during the Pleistocene glaciations. And so we want to know if highland species show different demographic trends compared to lowland. But also, we want to use genomics to study hybridization as a speciation mechanisms between diverging lineages that inhabit the lowlands.
  • #17: Fourth stage: study of adaptation. Our goal is to study the genomics regions that could be associated with adaptation. We know that some Oreobates species inhabit a wide diversity of habitats. For example, Oreobates cruralis (which is the species on the picture here) is usually found on the Amazonian lowland forests, but the same species can also be found in montane forests and even on dry forests. There are some other species that have found living on both highlands and lowlands environments. So we want to study if there are any genes related to adaptation to altitude or genes related to adaptation to dry environmental conditions.
  • #18: Regarding the methods, since we don’t have a closely related genome of reference, our initial idea was to (first) sequence the transcriptome for one species and (then) do whole genome sequencing for the other species. Our idea was to use (then) an exome assembly to address our hypothesis. However, we realized that this was quite expensive and also a big waste of resources (since –depending on the genome size– we would be throwing away about 95% of the data). For this reason, we decided to first try to find out the genome size for Oreobates species.
  • #19: Actually I spent a lot of time try to develop a method for estimating genome size based on real-time PCR. However, I’m not going into details since this work was not conclusive enough. Anyway, amphibians have a very large variation in the genome size (as you can see) compared to other vertebrates such as birds, reptiles or mammals. For this reason we decided to change our initial approach, and use a small representation of the genome.
  • #20: In my PhD project, we will use the transcriptome as a reduced representation of the genome
  • #21: How to do it? The main problem is that we don’t have RNA for all the Oreobates species. As you may know, RNA is very instable at room temperature and so it’s very difficult to preserve it. In fact, we only had RNA available for one single Oreobates species (O. cruralis). So, we extracted the RNA for the three tissues that we had (liver, intestine and spleen) and did transcriptome sequencing. Then, I spent 4 months in Sweden, doing fancy bioinformatics analysis to reconstruct a de-novo transcriptome assembly for one species. Don’t’ worry I’m going to explain you this a little bit more in detail in the next slide. Using this transcriptome as a reference, now we are working on a transcriptome-based exon capture (which, by the way, is a very novel technique) that will allow us to obtain the transcriptome for all the others Oreobates species. I will also explain this on the last slide. And so finally, once we get all the data, we will be able to test the initial hyphotesis.
  • #24: Pfam: database of protein families SwissProt: curated protein sequence database Xenopus: protein database of Xenopus tropicalis GO_Pfam: Gene Ontology (using only Pfam results) GO_SwissProt: Gene Ontology (using only SwissProt results) GO_SwissProt_Pfam: Gene Ontology (using SwissProt + Pfam results) Kegg: Kyoto Encyclopedia of Genes and Genomes eggNOG: Evolutionary Genealogy of Genes: Non-supervised Orthologous Groups
  • #27: In brief, using the transcriptome of one species, we will synthesize thousands of short nucleic acid probes that will specifically hybridize to genomic DNA of other closely related species (Maricic et al. 2010). Thereby, we will be able to retain and enrich homologous sequence-based loci across our target species (that is, transcriptome enrichment). Our main goal is to use this transcriptomic information to assess divergence patterns between Oreobates species. As this genus includes species adapted to highland and lowland environments, it could be interesting to know if highlands frogs have differences in gene expression compared to lowland frogs.
  • #29: My thesis is about diversification so…
  • #30: My thesis is about diversification so…
  • #31: My thesis is about diversification so…