. 2014 Sep 18;513(7518):375-381.

doi: 10.1038/nature13726. Epub 2014 Sep 3.

The genomic substrate for adaptive radiation in African cichlid fish

David Brawand^#^{1

2}, Catherine E Wagner^#^{3

4}, Yang I Li^#², Milan Malinsky^{5

6}, Irene Keller⁴, Shaohua Fan⁷, Oleg Simakov^{7

8}, Alvin Y Ng⁹, Zhi Wei Lim⁹, Etienne Bezault¹⁰, Jason Turner-Maier¹, Jeremy Johnson¹, Rosa Alcazar¹¹, Hyun Ji Noh¹, Pamela Russell¹², Bronwen Aken⁶, Jessica Alföldi¹, Chris Amemiya¹³, Naoual Azzouzi¹⁴, Jean-François Baroiller¹⁵, Frederique Barloy-Hubler¹⁴, Aaron Berlin¹, Ryan Bloomquist¹⁶, Karen L Carleton¹⁷, Matthew A Conte¹⁷, Helena D'Cotta¹⁵, Orly Eshel¹⁸, Leslie Gaffney¹, Francis Galibert¹⁴, Hugo F Gante¹⁹, Sante Gnerre¹, Lucie Greuter^{3

4}, Richard Guyon¹⁴, Natalie S Haddad¹⁶, Wilfried Haerty², Rayna M Harris²⁰, Hans A Hofmann²⁰, Thibaut Hourlier⁶, Gideon Hulata¹⁸, David B Jaffe¹, Marcia Lara¹, Alison P Lee⁹, Iain MacCallum¹, Salome Mwaiko³, Masato Nikaido²¹, Hidenori Nishihara²¹, Catherine Ozouf-Costaz²², David J Penman²³, Dariusz Przybylski¹, Michaelle Rakotomanga¹⁴, Suzy C P Renn¹⁰, Filipe J Ribeiro¹, Micha Ron¹⁸, Walter Salzburger¹⁹, Luis Sanchez-Pulido², M Emilia Santos¹⁹, Steve Searle⁶, Ted Sharpe¹, Ross Swofford¹, Frederick J Tan²⁴, Louise Williams¹, Sarah Young¹, Shuangye Yin¹, Norihiro Okada^{21

25}, Thomas D Kocher¹⁷, Eric A Miska⁵, Eric S Lander¹, Byrappa Venkatesh⁹, Russell D Fernald¹¹, Axel Meyer⁷, Chris P Ponting², J Todd Streelman¹⁶, Kerstin Lindblad-Toh^{1

26}, Ole Seehausen^{3

4}, Federica Di Palma^{1

27}

Affiliations

¹ Broad Institute of MIT and Harvard, Cambridge, Massachusetts 02142, USA.
² MRC Functional Genomics Unit, University of Oxford, Oxford OX1 3QX, UK.
³ Department of Fish Ecology and Evolution, Eawag Swiss Federal Institute of Aquatic Science and Technology, Center for Ecology, Evolution & Biogeochemistry, CH-6047 Kastanienbaum, Switzerland.
⁴ Division of Aquatic Ecology, Institute of Ecology & Evolution, University of Bern, CH-3012 Bern, Switzerland.
⁵ Gurdon Institute, Cambridge CB2 1QN, UK.
⁶ Wellcome Trust Sanger Institute, Hinxton CB10 1SA, UK.
⁷ Department of Biology, University of Konstanz, D-78457 Konstanz, Germany.
⁸ European Molecular Biology Laboratory, 69117 Heidelberg, Germany.
⁹ Institute of Molecular and Cell Biology, ASTAR, 138673 Singapore.
¹⁰ Department of Biology, Reed College, Portland, Oregon 97202, USA.
¹¹ Biology Department, Stanford University, Stanford, California 94305-5020, USA.
¹² Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, California 91125, USA.
¹³ Benaroya Research Institute at Virginia Mason, Seattle, Washington 98101, USA.
¹⁴ Institut Génétique et Développement, CNRS/University of Rennes, 35043 Rennes, France.
¹⁵ CIRAD, Campus International de Baillarguet, TA B-110/A, 34398 Montpellier cedex 5, France.
¹⁶ School of Biology, Georgia Institute of Technology, Atlanta, Georgia 30332-0230, USA.
¹⁷ Department of Biology, University of Maryland, College Park, Maryland 20742, USA.
¹⁸ Animal Genetics, Institute of Animal Science, ARO, The Volcani Center, Bet-Dagan, 50250 Israel.
¹⁹ Zoological Institute, University of Basel, CH-4051 Basel, Switzerland.
²⁰ Department of Integrative Biology, Center for Computational Biology and Bioinformatics; The University of Texas at Austin, Austin, Texas 78712, USA.
²¹ Department of Biological Sciences, Tokyo Institute of Technology, Tokyo, 226-8501 Yokohama, Japan.
²² Systématique, Adaptation, Evolution, National Museum of Natural History, 75005 Paris, France.
²³ Institute of Aquaculture, University of Stirling, Stirling FK9 4LA, UK.
²⁴ Carnegie Institution of Washington, Department of Embryology, 3520 San Martin Drive Baltimore, Maryland 21218, USA.
²⁵ National Cheng Kung University, Tainan City, 704 Taiwan.
²⁶ Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, 751 23 Uppsala, Sweden.
²⁷ Vertebrate and Health Genomics, The Genome Analysis Centre, Norwich NR18 7UH, UK.

^# Contributed equally.

PMID: 25186727
PMCID: PMC4353498
DOI: 10.1038/nature13726

The genomic substrate for adaptive radiation in African cichlid fish

David Brawand et al. Nature. 2014.

. 2014 Sep 18;513(7518):375-381.

doi: 10.1038/nature13726. Epub 2014 Sep 3.

VSports - Authors

Affiliations

¹ Broad Institute of MIT and Harvard, Cambridge, Massachusetts 02142, USA.
² MRC Functional Genomics Unit, University of Oxford, Oxford OX1 3QX, UK.
³ Department of Fish Ecology and Evolution, Eawag Swiss Federal Institute of Aquatic Science and Technology, Center for Ecology, Evolution & Biogeochemistry, CH-6047 Kastanienbaum, Switzerland.
⁴ Division of Aquatic Ecology, Institute of Ecology & Evolution, University of Bern, CH-3012 Bern, Switzerland.
⁵ Gurdon Institute, Cambridge CB2 1QN, UK.
⁶ Wellcome Trust Sanger Institute, Hinxton CB10 1SA, UK.
⁷ Department of Biology, University of Konstanz, D-78457 Konstanz, Germany.
⁸ European Molecular Biology Laboratory, 69117 Heidelberg, Germany.
⁹ Institute of Molecular and Cell Biology, ASTAR, 138673 Singapore.
¹⁰ Department of Biology, Reed College, Portland, Oregon 97202, USA.
¹¹ Biology Department, Stanford University, Stanford, California 94305-5020, USA.
¹² Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, California 91125, USA.
¹³ Benaroya Research Institute at Virginia Mason, Seattle, Washington 98101, USA.
¹⁴ Institut Génétique et Développement, CNRS/University of Rennes, 35043 Rennes, France.
¹⁵ CIRAD, Campus International de Baillarguet, TA B-110/A, 34398 Montpellier cedex 5, France.
¹⁶ School of Biology, Georgia Institute of Technology, Atlanta, Georgia 30332-0230, USA.
¹⁷ Department of Biology, University of Maryland, College Park, Maryland 20742, USA.
¹⁸ Animal Genetics, Institute of Animal Science, ARO, The Volcani Center, Bet-Dagan, 50250 Israel.
¹⁹ Zoological Institute, University of Basel, CH-4051 Basel, Switzerland.
²⁰ Department of Integrative Biology, Center for Computational Biology and Bioinformatics; The University of Texas at Austin, Austin, Texas 78712, USA.
²¹ Department of Biological Sciences, Tokyo Institute of Technology, Tokyo, 226-8501 Yokohama, Japan.
²² Systématique, Adaptation, Evolution, National Museum of Natural History, 75005 Paris, France.
²³ Institute of Aquaculture, University of Stirling, Stirling FK9 4LA, UK.
²⁴ Carnegie Institution of Washington, Department of Embryology, 3520 San Martin Drive Baltimore, Maryland 21218, USA.
²⁵ National Cheng Kung University, Tainan City, 704 Taiwan.
²⁶ Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, 751 23 Uppsala, Sweden.
²⁷ Vertebrate and Health Genomics, The Genome Analysis Centre, Norwich NR18 7UH, UK.

^# Contributed equally.

PMID: 25186727
PMCID: PMC4353498 (V体育2025版)
DOI: 10.1038/nature13726

"VSports" Abstract

Cichlid fishes are famous for large, diverse and replicated adaptive radiations in the Great Lakes of East Africa. To understand the molecular mechanisms underlying cichlid phenotypic diversity, we sequenced the genomes and transcriptomes of five lineages of African cichlids: the Nile tilapia (Oreochromis niloticus), an ancestral lineage with low diversity; and four members of the East African lineage: Neolamprologus brichardi/pulcher (older radiation, Lake Tanganyika), Metriaclima zebra (recent radiation, Lake Malawi), Pundamilia nyererei (very recent radiation, Lake Victoria), and Astatotilapia burtoni (riverine species around Lake Tanganyika). We found an excess of gene duplications in the East African lineage compared to tilapia and other teleosts, an abundance of non-coding element divergence, accelerated coding sequence evolution, expression divergence associated with transposable element insertions, and regulation by novel microRNAs. In addition, we analysed sequence data from sixty individuals representing six closely related species from Lake Victoria, and show genome-wide diversifying selection on coding and regulatory variants, some of which were recruited from ancient polymorphisms. We conclude that a number of molecular mechanisms shaped East African cichlid genomes, and that amassing of standing variation during periods of relaxed purifying selection may have been important in facilitating subsequent evolutionary diversification. VSports手机版.

PubMed Disclaimer

"V体育2025版" Figures

**Extended Data Figure 1. Genome assembly and evolutionary rates**
a, Genome assembly and annotation. b, Genome-wide dN/dS. Rates are calculated from 20 resampled sets of 200 orthologous genes. Gene annotations from interspecies projections (see Methods in Supplementary Information) were excluded from the data set.

**Extended Data Figure 2. Rapid Evolution of EDNRB1**
a, Alignments of EDNRB1 in cichlids with human (HS), zebrafish (DAR) and medaka (ORL). Black star denotes site shown to be required to activate SRF in human by interacting with the G protein G13 (ref. 40). Red star denotes site that may affect the anchoring of the C terminus of EDNRB1to the transmembrane domain. Highlighted are amino acid substitution in the ancestor of haplochromine and lamprologini (blue) and in the ancestor of haplochromine (red). b, Location of substitutions on 7 transmembrane domain representation (Adapted from ref. Science 318, 1453–1455. Reprinted with permission from AAAS.). c, Sites (spheres) on the structure of the human kappa opioid receptor in complex (4DJH). Only the right homodimer is annotated.

**Extended Data Figure 3. Duplication in the cichlid genomes**
a, The number of the recently duplicated genomic regions identified by the read depth method in the five East African cichlid genomes. The numbers in red is the number of duplicated genes and the numbers in black is the corresponding branch length. b, Summary of duplication regions in the five African cichlid genomes. c, Venn diagram of the duplicated genes detected by aCGH across cichlid species relative to *O. niloticus*. d, Expression patterns of duplicate genes. Matrix represents the expression level of retained duplicate genes from the cichlid common ancestor in the specified tissues. Expression is showed as an inverse logit function of log₂-transformed, relative sequence fragment numbers. Uncoloured fields designate missing expression data or absence of either of the paralogue copies in the annotation set.

**Extended Data Figure 4. Cichlid OR and *TAAR* genes**
a, Cichlid OR and *TAAR* genes identified in this study. b, PHYML tree based on the fish TAAR and OR amino acid sequences. A phylogeny tree was constructed with all OR and TAAR cichlid proteins identified in this study (n = 503 + 119) plus 229 OR and 173 *TAAR* genes identified in zebrafish, fugu, tetraodon, medaka and stickleback. The amino-acid sequences were aligned with MAFFT version 7 and a tree constructed with PHYML and visualized with Fig Tree (version 1.3.1). The TAAR branches are in pink and the OR branches in blue or yellow. Colours indicate the composition of the branches. Dark blue branches are made of cichlid OR only, light blue indicates the presence of cichlid and model fish OR proteins. Yellow is for branches made of OR model fish proteins, only. Light pink branches correspond to cichlid and model fish TAAR proteins and dark pink to model fish TAAR proteins only. Letters correspond to family names.

**Extended Data Figure 5. Comparison of TEs among cichlids and other vertebrate genomes**
a, Repeat content of selected vertebrate genomes. Table, no legend. b, Proportions of TEs in the genomes. c, Proportions of each TE class among all TEs. The TE proportions are much lower in cichlid genomes than that in zebrafish.

**Extended Data Figure 6. A comparison of TEs in the African cichlids and Medaka genomes**
a–f, The x-axis indicates a specific TE family at a given divergence from the consensus sequence and y-axis indicates its percentage of the genome.

**Extended Data Figure 7. Cichlid transposable elements**
a, Association between TE insertions and gene expression levels of orthologous genes. All 5 cichlids are merged into the same data set. Groupings are based on whether one gene copy lies within 20 kb up or downstream of a TE. b, Orientation bias of transposable elements within or near non-duplicated genes. TE orientation bias in intron sequence of 5 cichlid species. Bias is shown as log₂(sense/antisense) of TE counts. c, Orientation bias of transposable elements in introns of protein coding genes. The x-axis denotes the maximum age of the TEs as divergence from the consensus sequence. The y-axis shows the proportion of TE insertions in the sense of transcription. Data points with large confidence intervals (exceeding the display range) are omitted. d, Orientation bias of LINE insertions in introns in 4% divergence wide windows in *O. niloticus*, *N. brichardi* and combined haplochromines. Proportion of sense oriented LINEs in introns is shown on the y-axis. Age is shown on the x-axis as percent divergence from the TE consensus.

**Extended Data Figure 8. Reporter gene expression of a selectedO. niloticus hCNE–*P. nyererei* aCNE pair in transgenic zebrafish**
a, *O. niloticus* *Pbx1a* locus showing the conservation track and alignment of an hCNE (LG18.20714) in *O. latipes* and East Africa cichlids. b, Reporter gene expression in 72 hours post-fertilization (hpf) G1 transgenic zebrafish. Expression is shown for the hCNE in *O. niloticus* and the corresponding aCNE in *P. nyererei*. The *P. nyererei* aCNE also shows expression in circulating blood cells.

**Extended Data Figure 9. Reporter gene expression of a selected hCNE–aCNE pairs in transgenic zebrafish (G0)**
a, b, Comparison of expression pattern driven by *O. niloticus* and *N. brichardi* aCNE #911 (*UNCX* locus) in 72 hours post-fertilization (hpf) zebrafish embryos. c, d, Comparison of expression pattern driven by tilapia and *N. brichardi* aCNE #7012 (*SERPINH1* locus) in 72 hpf zebrafish embryos. e, f, Comparison of expression pattern driven by tilapia and *N. brichardi* aCNE #1649 (*TBX2* locus) in 72 hpf zebrafish embryos. g, h, Comparison of expression pattern driven by tilapia and *M. zebra* aCNE #26432 (*FOXP4* locus) in 72 hpf zebrafish embryos. i, j, Comparison of expression pattern driven by tilapia and *A. burtoni* aCNE #5509 (*PROX1* locus) in 72 hpf zebrafish embryos.

**Extended Data Figure 10. Cichlid microRNAs**
a, Novelty in microRNAs mapped on the phylogenetic tree of the five cichlid species. Complementary expression of novel cichlid miRNA *mir-10032* (c, e, g) and predicted target gene *neurod2* (b, d, f) in stage 23 (9–10 days post-fertilization) *Metriaclima zebra* embryos. d, e are 18-μm sagittal sections. In d and e, arrows point to expression in the medulla (left), cerebellum (middle) and optic tectum (right). *neurod2* is expressed in the neural tube (b, f), while *mir-10032* is expressed in the surrounding somites (c, g). In all panels, anterior is to the right.

**Figure 1. The adaptive radiation of African cichlid fish**
Top left, map of Africa showing lakes in which cichlid fish have radiated. Right, the five sequenced species: *Pundamilia nyererei* (endemic of Lake Victoria); *Neolamprologus brichardi* (endemic of Lake Tanganyika); *Metriaclima zebra* (endemic of Lake Malawi); *Oreochromis niloticus* (from rivers across northern Africa); *Astatotilapia burtoni* (from rivers connected to Lake Tanganyika). Major ecotypes are shown from each lake: a, pelagic zooplanktivore; b, rock-dwelling algaescraper; c, paedophage (absent from Lake Tanganyika); d, scale eater; e, snail crusher; f, reef-dwelling planktivore; g, lobe-lipped insect eater; h, pelagic piscivore; i, ancestral river-dweller also found in lakes (absent from Lake Tanganyika). Bottom left, phylogenetic tree illustrating relationships between the five sequenced species (red), major adaptive radiations and major river lineages. The tree is from ref. , pruned to the major lineages. Upper timescale (4), lower timescale (32). Photos by Ad Konings (Tanganyika a, b, d, e, g, h; Malawi a, c, d, e, f, g, h, i), O.S. (Victoria a–g, i; Malawi b), Frans Witte (Victoria h), W.S. (Tanganyika f), Oliver Selz (Victoria f, *A. burtoni*), Marcel Haesler (*O. niloticus*).

**Figure 2. Gene duplication in the ancestry of East African lake cichlids**
Black numbers represents species divergence calculated as neutral genomic divergence between the sequenced species using ~2.7 million fourfold degenerate sites from the alignment of 9 teleost genomes. This neutral substitution model suggests ~2% pairwise divergence between the three haplochromines and a ~6% divergence to *N. brichardi*. Red numbers represent duplicated genes. Asterisks indicate excluded branches owing to incomplete lineage sorting in haplochromines or weak support of consensus species tree.

**Figure 3. Novel cichlid microRNAs**
a–f, Complementary expression of *mir*-10029 (b, d, f) and its predicted target gene *bmpr1b* (a, c, e) in stage 18 (6 days post-fertilization) *Metriaclima zebra* embryos. c–f are 18-μm sagittal sections. In c and d arrows point to expression (black) or lack of expression (white) in the somites, presumptive cerebellum, and optic tectum (from left to right). In e and f, arrows point to expression and lack of expression in the somites (dorsal) and the gut (ventral). In all panels, anterior is to the right.

**Figure 4. Genomic divergence stems from incomplete lineage sorting (ILS) and both old and novel coding and noncoding variation**
a, Coalescence times and trees supporting ILS among the genomes of allopatric East African cichlid lineages were inferred by coalHMM. The most common genealogy matches the known species tree and represents a *M. zebra*–*P. nyererei* coalescence that falls between the two speciation times, T_zn (speciation *M. zebra*–*P. nyererei*) and T_znb (speciation *M. zebra*–*P. nyererei*–*A. burtoni*). In genealogies 1 (dashed line), 2, and 3, all coalescence events are ancient and occur before time T_znb. b, Phylogenetic analysis of RAD-sequence data showing well-supported differentiation among young Victoria species. The complete data set (top) renders the genus *Mbipia* non-monophyletic, exclusion of the top 1% divergent loci (bottom) supports monophyly of each genus. c, Genomic divergence in paired comparisons of Lake Victoria cichlids (per-site F_ST; black/grey are chromosomes). Sister species from top: *Pundamilia nyererei*/*P. pundamilia* and *Mbipia lutea*/*M. mbipi* differ in male breeding coloration but have conserved morphology; *Neochromis omnicaeruleus*/N. sp. “unicuspid scraper” and distant relatives *P. pundamilia*/*M. mbipi* and *P. nyererei*/M. lutea have similar coloration but differ in morphology. Red-highlighted SNPs indicate significantly divergent sites between colour-contrasting species, but not between same-colour species. Bar plots show the proportion of SNPs in four annotation categories: exons (orange), introns (dark blue), 25-kb flanking genes (turquoise), or none of the above (grey), for thresholds of increasing F_ST. In “All sites” and “Ancient variant sites” analyses, symbols indicate an excess of SNPs in a given annotation category compared to expectations from the full data set or from all non-ancient variant sites, respectively (FDR q-values: *q < 0.05; †q = 0.05), (Supplementary Information, Data Portals, Supplementary Population Genomics FTP files).

See this image and copyright information in PMC

V体育平台登录 - Comment in

Evolutionary biology: Radiating genomes.
Jiggins CD. Jiggins CD. Nature. 2014 Sep 18;513(7518):318-9. doi: 10.1038/nature13742. Epub 2014 Sep 3. Nature. 2014. PMID: 25186726 No abstract available.
"VSports注册入口" Evolution: a fishy tale of adaptive radiation.
Jones B. Jones B. Nat Rev Genet. 2014 Nov;15(11):704. doi: 10.1038/nrg3835. Epub 2014 Sep 23. Nat Rev Genet. 2014. PMID: 25246195 No abstract available.

References

1. Darwin C. On the Origin of Species. 6th edn John Murray; 1859.
1. Simpson GG. Tempo and Mode in Evolution. Columbia Univ. Press; 1944.
1. Kocher TD. Adaptive evolution and explosive speciation: the cichlid fish model. Nature Rev. Genet. 2004;5:288–298. - "VSports手机版" PubMed
1. Wagner CE, Harmon LJ, Seehausen O. Ecological opportunity and sexual selection together predict adaptive radiation. Nature. 2012;487:366–369. - V体育平台登录 - PubMed
1. Meyer A. Morphometrics and allometry in the trophically polymorphic cichlid fish, Cichlasoma citrinellum: alternative adaptations and ontogenic changes in shape. J. Zool. 1990;221:237–260.

"V体育官网入口" Publication types

Actions
"V体育平台登录" Actions

MeSH terms

Actions
"VSports手机版" Actions
Actions
Actions
V体育ios版 - Actions
Actions
V体育官网 - Actions
VSports - Actions
Actions
VSports注册入口 - Actions
Actions
Actions
"V体育官网" Actions
Actions
Actions

Substances (VSports最新版本)

Actions
Actions (VSports手机版)

Associated data

Grants and funding (V体育安卓版)

LinkOut - more resources (V体育官网)

Full Text Sources
"VSports手机版" Other Literature Sources
- The Lens - Patent Citations Database
- scite Smart Citations
Molecular Biology Databases
- SILVA

"V体育安卓版" Save citation to file

VSports app下载 - Email citation

Add to Collections

Add to My Bibliography

Your saved search

Create a file for external citation management software

VSports注册入口 - Your RSS Feed

The genomic substrate for adaptive radiation in African cichlid fish

Affiliations

The genomic substrate for adaptive radiation in African cichlid fish

VSports - Authors

Affiliations

"VSports" Abstract

"V体育2025版" Figures

V体育平台登录 - Comment in

References

"V体育官网入口" Publication types

MeSH terms

Substances (VSports最新版本)

Associated data

Grants and funding (V体育安卓版)

LinkOut - more resources (V体育官网)

Full Text Sources

"VSports手机版" Other Literature Sources

Molecular Biology Databases