Top Banner
Asian Seabass Genome Project: A status report László Orbán, on behalf of the The Asian Seabass Genome Consortium Temasek Life Sciences Laboratory, Singapore; Dobzhansky Center, St. Petersburg, Russia; Pacific Biosystems, USA; SANBI, Capetown, South Africa; Max Planck Institute for Mol. Gen., Berlin, Germany; Chinese University of HK, Hong Kong; Centre for Comparative Genomics, Murdoch University, Australia; IGIB, New Delhi, India; Georgikon Faculty, University of Pannonnia, Keszthely, Hungary;
27

Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Jul 03, 2020

Download

Documents

dariahiddleston
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Asian Seabass Genome Project:

A status report

László Orbán, on behalf of the

The Asian Seabass Genome Consortium

Temasek Life Sciences Laboratory, Singapore; Dobzhansky Center, St.

Petersburg, Russia; Pacific Biosystems, USA; SANBI, Capetown,

South Africa; Max Planck Institute for Mol. Gen., Berlin, Germany;

Chinese University of HK, Hong Kong; Centre for Comparative

Genomics, Murdoch University, Australia; IGIB, New Delhi, India;

Georgikon Faculty, University of Pannonnia, Keszthely, Hungary;

Page 2: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Population 7 billion Depleted oceans

More people, less fish

Source: National Geographic

endoftheline.com

Page 3: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

• Introducing the Asian seabass and the selection

program;

• Status report on the Genome (and Transcriptome)

Project;

• Early applications: phylogeography and sex change;

• Summary.

Outline

Page 4: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Asian seabass

(Lates calcarifer)

• Barramundi (AUS);

• Euryhaline, catadromous predator;

• Protandrous (male-first) hermaphrodite;

• Distributed over a wide geographical area;

• Cultured mostly by smaller farms,a few

selection programs;

Page 5: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Asian seabass

(Lates calcarifer)

• Barramundi (AUS);

• Euryhaline, catadromous predator;

• Protandrous (male-first) hermaphrodite;

• Distributed over a wide geographical area;

• Cultured mostly by smaller farms, few

selection programs;

• TLL-AVA collaboration: Selection program

in its 9th year – first mature F2s are ready

for crosses;

• Increased growth rate.

Page 6: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Asian seabass selection

program (timeline)

• 2004-2011: marker-assisted selection

program with the Yue group

• Increased growth rate – F2 grows >20%

faster than unselected;

• 2011: received S$10M grant from NRF

for 5 yrs;

• Genomic selection – polygenic traits;

• Nutrigenomics and disease resistance;

• Start parallel program for salt tolerant

tilapia.

Page 7: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

• Introducing the Asian seabass and the selection

program;

• Status report on the Genome (and Transcriptome)

Project;

• Early applications: phylogeography and sex change;

• Summary.

Outline

Page 8: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Genome size: 700 Mb

2n = 48

Genetic linkage map (high density)

Physical map

Expected gene count: 26,000-27,000

Repeat inventory

Asian seabass genome – what was known

2um

Kuznetsova et al., Front. Genet. (2014)

Page 9: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Asian seabass

Genome Assembly Plan

BAC end

sequencing

1X Coverage

Illumina HiSeq 2500

2X100 nt PE

80X Coverage

PacBio SMRT

80 Smart Cells

30X Coverage

v1

v3

v4

Repeat Inventory

Transcriptome

QTL on LG2

Mate-pair libraries (?)

Pooled BAC sequencing

Improved High Quality Draft Genome

Gap filling (?)

v2

Lok Lab, CUHK, Hong Kong

Sivasubbu Lab, IGIB, New Delhi, India

Shubha Vij & team

Visit Poster 88 on the

Sequencing and assembly of a

10 Mb region of the Asian

seabass genome containing

growth-associated QTLs

Page 10: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Asian seabass - Genome Sequence Datasets

A partially inbred Asian seabass

individual

Illumina

HiSeq 2X100 nt PE

500&750 bp

80X

coverage

PacBio

SMRT ~4 kb avg read

length

~10 Kb library

30X

coverage

Sanger

seq Two BAC

libraries-

~120 kb insert

11.5K seq

1X coverage

PacBio

SMRT ~8 kb avg read

length

~20 Kb library

60X

coverage

Page 11: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Genome assembly - the first try

PacBio RS

Long reads

# contigs 3,807

N50 1.2 Mb

Maximum 18.9 Mb

de novo

assembly 90X

Siddarth Singh, Pacific Biosystems

Mike Schatz, CSH Labs

Page 12: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Final genome assembly

PacBio RS

Long reads

Illumina HiSeq

Short reads

# contigs 11,969,803

N50 1,001 bp

Maximum 26,868 bp

denovo

assembly de novo

assembly 90X 80X

Lok Lab, CUHK, Hong Kong

Sivasubbu Lab, IGIB, New Delhi, India

Siddarth Singh, Pacific Biosystems

Mike Schatz, CSH Labs

# contigs 3,807

N50 1.2 Mb

Maximum 18.9 Mb

Page 13: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

• Validation by 80X Illumina-reads mapping

− 96% of reads mapped as proper pairs

• Gene annotation underway

• Generated optical map data - de novo assembled and used to

place PacBio genome contigs

− 79.8% similarity between the optical map and sequence assembly

• Alternative Falcon assembly underway

Current status and improvement efforts

PacBio Optical Map

# contigs 3,807 3,333

Contig N50 1.2 Mb 6.2 Mb

Maximum 18.9 Mb 20.6 Mb

Page 14: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Transcriptome sequencing and assembly

• >1 billion reads from various organs of multiple individuals

and 3 NGS platforms were assembled in a step-wise manner

• >80% of the expected protein-coding loci obtained,

58% of these represented by a predicted FL-cDNA sequence

Total length (bp) 262,023,963

Number of contigs 267,616

Number of contigs ≥1 kb 70,588

Max length (bp) 31,251

Average length (bp) 979

Thevasagayam et al., J. Mar. Sci. Eng. (2015)

Page 15: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

• Introducing the Asian seabass and the selection

program;

• Status report on the Genome (and Transcriptome)

Project;

• Early applications: phylogeography and sex change;

• Summary.

Outline

Page 16: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Potential benefits of sequenced genomes

• Info on coding regions;

• Info on regulatory regions & pathways;

• Comparative/evolutionary genomics;

• New platforms (chips, RNAseq, GBS,

methylome, etc.);

• Re-sequencing of variants;

• Rapid identification of mutations.

Page 17: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Molecular evidence points to the existence

of two Asian seabass species

95/91/91

99/100/99

99/100/99

95/91/91

0.005

Indian

Subcontinent

(50)

S-E Asia

(35)

L. niloticus

16S rDNA 0 .0 2

99/99/99

99/99/99

98/97/97

99/99/99

98/97/97

90/99/99

0.02

Indian

Subcontinent

+

Myanmar

(52)

Australia

(9)

S-E Asia

(24)

L. niloticus

COI

Ward et al., J. Fish Biol. (2008);

Pethiyagoda & Gill, Zootaxa (2012)

Vij et al., Frontiers Marine Science (2014)

Page 18: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Andrey Yurchenko, St.

Petersburg, Russia; unpubl.

Resequencing 65 genomes reveals

clear signs of admixture n the SEA region

K=

2

3

4

7

6

8

India-Eastern Coast India-Western

Coast

PG AU-QLD AU-D ID-K ID-SJ KH TH-EC PH

AU/PG S-E Asia/PH Indian region

ID-SU VN SG

Page 19: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Asian seabass (Lates calcarifer)

• Problem: Protandrous (male-first) hermaphrodite;

• Zebrafish can be used as a model to understand sex

change.

Orban et al., Mol. Cell. Endocrinol. 312: 35-41 (2009)

Page 20: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Natural sex reversal in Asian seabass

Guiguen et al., Env. Biol. Fish. 39: 231-247 (1994)

Unusual individuals: ‘reluctant’ males and primary females

Page 21: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Array-based transcriptomics: The transforming gonad returns to

a near-undifferentiated stage before initiating feminization

Jiang et al., Ph.D. Thesis (2014)

Page 22: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

Ravi, Jiang et al., RB&E 12:5 (2014)

Page 23: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

• Genome is sequenced (170X);

• Assembly is based on PacBio data, Illumina

validates;

• Optical mapping showed improvements;

• International consortium for annotation;

• Multiple benefits;

• (Seq and assembly of Mozambique tilapia

genome is in progress.)

Summary – Asian seabass Genome

Page 24: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

• Different approach from the ASB

− 6 Short-insert PE (SIPE) – 60X

− 3 Long-insert MP (1, 2, 3 kb) – 30X

− Lucigen Long-insert MP (8kb) – 3X

− Lucigen Fosmid MP (40kb) – 3X

− Pacbio sequence data pending

• De novo assembly is underway

− V1: SIPE-only assembly by MaSuRCA

Genome of the Mozambique tilapia

Assembly size 964 Mb

# of contigs 173,146

Max length 257.3 kb

N50 length 15.3 kb

C-value (pg): 0.81-1.0 (Animal Genome Size Database)

Diploid Chrom number (n): 44

Genome Size: ca. 1 Gb

Liew Woei Chang, Shen Xueyan & team

Page 25: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

RGG:

Shubha Vij

Inna Kuznetcova

Woei Chang Liew

Xueyan Shen

Natascha May

Prakki Sri Datta

Jolly M. Saju

Purushothaman K.

Shawn Ngoh

Pranjali Bhandare

Farm teams

Genhua Yue & team

Huan Sein Lim &

team

TLL Facilities

:

Funding:

Collaborators:

Alan Christoffels (Capetown)

Si Lok (HK)

Matt Bellgard (Murdoch)

Mike Shatz (CSH)

PacBio

National

Research

Foundation

RGG Alumni:

Preethi Ravi

Rajini Sreenivasan

Laura Casas

Alex Chang

Mohd. Sorowar

Hossain

Xingang Wang

Richard Bartfai

Keh-Weei Tzung

Hsiao Yuen Kwan

Doreen Lau

Junhui Jiang

Aleks Komissarov & Andrey Yurchenko

(St. Petersburg)

Heiner Kuhl (Berlin)

Dean Jerry (Queensland)

Page 26: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass
Page 27: Asian Seabass Genome Project: A status reportisga2015.acuigen.es/uploads/presentations/S5_Orban.pdfVisit Poster 88 on the Sequencing and assembly of a 10 Mb region of the Asian seabass

An outlook for the future

• SPIRAL

Past: MAS on

seabass

Present:

Genomic selection on

A. seabass & tilapia

Near future:

Integrated systems applied to

several fish species

Photo: http://faithoncampus.com/