John H. Postlethwait

Stacks: Building and Genotyping Loci De Novo From Short-Read Sequences

Advances in sequencing technology provide special opportunities for genotyping individuals with speed and thrift, but the lack of software to automate the calling of tens of thousands of genotypes over hundreds of individuals has hindered progress. Stacks is a software system that uses short-read sequence data to identify and genotype loci in a set of individuals either de novo or by comparison to a reference genome. From reduced representation Illumina sequence data, such as RAD-tags, Stacks can recover thousands of single nucleotide polymorphism (SNP) markers useful for the genetic analysis of crosses or populations. Stacks can generate markers for ultra-dense genetic linkage maps, facilitate the examination of population phylogeography, and help in reference genome assembly. We report here the algorithms implemented in Stacks and demonstrate their efficacy by constructing loci from simulated RAD-tags taken from the stickleback reference genome and by recapitulating and improving a genetic map of the zebrafish, Danio rerio.

The Zebrafish Glypican Knypek Controls Cell Polarity during Gastrulation Movements of Convergent Extension

Mutations in the zebrafish knypek locus impair gastrulation movements of convergent extension that narrow embryonic body and elongate it from head to tail. We demonstrate that knypek regulates cellular movements but not cell fate specification. Convergent extension movement defects in knypek are associated with abnormal cell polarity, as mutant cells fail to elongate and align medio-laterally. Positional cloning reveals that knypek encodes a member of the glypican family of heparan sulfate proteoglycans. Double mutant and overexpression analyses show that Knypek potentiates Wnt11 signaling, mediating convergent extension. These studies provide experimental and genetic evidence that glypican Knypek acts during vertebrate gastrulation as a positive modulator of noncanonical Wnt signaling to establish polarized cell behaviors underlying convergent extension movements.

The African coelacanth genome provides insights into tetrapod evolution.

The discovery of a living coelacanth specimen in 1938 was remarkable, as this lineage of lobe-finned fish was thought to have become extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features. Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues show the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.The discovery of a living coelacanth specimen in 1938 was remarkable, as this lineage of lobe-finned fish was thought to have become extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features. Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues show the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.

Genome Evolution and Meiotic Maps by Massively Parallel DNA Sequencing: Spotted Gar, an Outgroup for the Teleost Genome Duplication

Genomic resources for hundreds of species of evolutionary, agricultural, economic, and medical importance are unavailable due to the expense of well-assembled genome sequences and difficulties with multigenerational studies. Teleost fish provide many models for human disease but possess anciently duplicated genomes that sometimes obfuscate connectivity. Genomic information representing a fish lineage that diverged before the teleost genome duplication (TGD) would provide an outgroup for exploring the mechanisms of evolution after whole-genome duplication. We exploited massively parallel DNA sequencing to develop meiotic maps with thrift and speed by genotyping F1 offspring of a single female and a single male spotted gar (Lepisosteus oculatus) collected directly from nature utilizing only polymorphisms existing in these two wild individuals. Using Stacks, software that automates the calling of genotypes from polymorphisms assayed by Illumina sequencing, we constructed a map containing 8406 markers. RNA-seq on two map-cross larvae provided a reference transcriptome that identified nearly 1000 mapped protein-coding markers and allowed genome-wide analysis of conserved synteny. Results showed that the gar lineage diverged from teleosts before the TGD and its genome is organized more similarly to that of humans than teleosts. Thus, spotted gar provides a critical link between medical models in teleost fish, to which gar is biologically similar, and humans, to which gar is genomically similar. Application of our F1 dense mapping strategy to species with no prior genome information promises to facilitate comparative genomics and provide a scaffold for ordering the numerous contigs arising from next generation genome sequencing.

A pair of Sox: distinct and overlapping functions of zebrafish sox9 co-orthologs in craniofacial and pectoral fin development

Understanding how developmental systems evolve after genome amplification is important for discerning the origins of vertebrate novelties, including neural crest, placodes, cartilage and bone. Sox9 is important for the development of these features, and zebrafish has two co-orthologs of tetrapod SOX9 stemming from an ancient genome duplication event in the lineage of ray-fin fish. We have used a genotype-driven screen to isolate a mutation deleting sox9b function, and investigated its phenotype and genetic interactions with a sox9a null mutation. Analysis of mutant phenotypes strongly supports the interpretation that ancestral gene functions partitioned spatially and temporally between Sox9 co-orthologs. Distinct subsets of the craniofacial skeleton, otic placode and pectoral appendage express each gene, and are defective in each single mutant. The double mutant phenotype is additive or synergistic. Ears are somewhat reduced in each single mutant but are mostly absent in the double mutant. Loss-of-function animals from mutations and morpholino injections, and gain-of-function animals injected with sox9a and sox9b mRNAs showed that sox9 helps regulate other early crest genes, including foxd3, sox10, snai1b and crestin, as well as the cartilage gene col2a1 and the bone gene runx2a; however, tfap2a was nearly unchanged in mutants. Chondrocytes failed to stack in sox9a mutants, failed to attain proper numbers in sox9b mutants and failed in both morphogenetic processes in double mutants. Pleiotropy can cause mutations in single copy tetrapod genes, such as Sox9, to block development early and obscure later gene functions. By contrast, subfunction partitioning between zebrafish co-orthologs of tetrapod genes, such as sox9a and sox9b, can relax pleiotropy and reveal both early and late developmental gene functions.

MicroRNA Mirn140 modulates Pdgf signaling during palatogenesis.

Disruption of signaling pathways such as those mediated by sonic hedgehog (Shh) or platelet-derived growth factor (Pdgf) causes craniofacial abnormalities, including cleft palate. The role that microRNAs play in modulating palatogenesis, however, is completely unknown. We show that, in zebrafish, the microRNA Mirn140 negatively regulates Pdgf signaling during palatal development, and we provide a mechanism for how disruption of Pdgf signaling causes palatal clefting. The pdgf receptor alpha (pdgfra) 3′ UTR contained a Mirn140 binding site functioning in the negative regulation of Pdgfra protein levels in vivo. pdgfra mutants and Mirn140-injected embryos shared a range of facial defects, including clefting of the crest-derived cartilages that develop in the roof of the larval mouth. Concomitantly, the oral ectoderm beneath where these cartilages develop lost pitx2 and shha expression. Mirn140 modulated Pdgf-mediated attraction of cranial neural crest cells to the oral ectoderm, where crest-derived signals were necessary for oral ectodermal gene expression. Mirn140 loss of function elevated Pdgfra protein levels, altered palatal shape and caused neural crest cells to accumulate around the optic stalk, a source of the ligand Pdgfaa. These results suggest that the conserved regulatory interactions of mirn140 and pdgfra define an ancient mechanism of palatogenesis, and they provide candidate genes for cleft palate.

Plasticity of Animal Genome Architecture Unmasked by Rapid Evolution of a Pelagic Tunicate

Ocean Dweller Sequenced The Tunicates, which include the solitary free-swimming larvaceans that are a major pelagic component of our oceans, are a basal lineage of the chordates. In order to investigate the major evolutionary transition represented by these organisms, Denoeud et al. (p. 1381, published online 18 November) sequenced the genome of Oikopleura dioica, a chordate placed by phylogeny between vertebrates and amphioxus. Surprisingly, the genome showed little conservation in genome architecture when compared to the genomes of other animals. Furthermore, this highly compacted genome contained intron gains and losses, as well as species-specific gene duplications and losses that may be associated with development. Thus, contrary to popular belief, global similarities of genome architecture from sponges to humans are not essential for the preservation of ancestral morphologies. A metazoan genome departs from the organization that appears rigidly established in other animal phyla. Genomes of animals as different as sponges and humans show conservation of global architecture. Here we show that multiple genomic features including transposon diversity, developmental gene repertoire, physical gene order, and intron-exon organization are shattered in the tunicate Oikopleura, belonging to the sister group of vertebrates and retaining chordate morphology. Ancestral architecture of animal genomes can be deeply modified and may therefore be largely nonadaptive. This rapidly evolving animal lineage thus offers unique perspectives on the level of genome plasticity. It also illuminates issues as fundamental as the mechanisms of intron gain.

Mapping of Mhc class I and class II regions to different linkage groups in the zebrafish, Danio rerio

Abstract The mammalian major histocompatibility complex (Mhc) consists of three closely linked regions, I, II, and III, occupying a single chromosomal segment. The class I loci in region I and the class II loci in region II are related in their structure, function, and evolution. Region III, which is intercalated between regions I and II, contains loci unrelated to the class I and II loci, and to one another. There are indications that a similar Mhc organization exists in birds and amphibians. Here, we demonstrate that in the zebrafish (Danio rerio), a representative of the teleost fishes, the class II loci are divided between two linkage groups which are distinct from the linkage group containing the class I loci. The β2-microglobulin-encoding gene is loosely linked to one of the class II loci. The gene coding for complement factor B, which is one of the region III genes in mammals, is linked neither to the class I nor to the class II loci in the zebrafish. These results, combined with preliminary data suggesting that the class I and class II regions in another order of teleost fish are also in different linkage groups, indicate that close linkage of the two regions is not necessary either for regulation of expression or for co-evolution of the class I and class II loci. They also raise the question of whether linkage of the class I and class II loci in tetrapods is a primitive or derived character.

Automated identification of conserved synteny after whole-genome duplication

An important objective for inferring the evolutionary history of gene families is the determination of orthologies and paralogies. Lineage-specific paralog loss following whole-genome duplication events can cause anciently related homologs to appear in some assays as orthologs. Conserved synteny-the tendency of neighboring genes to retain their relative positions and orders on chromosomes over evolutionary time-can help resolve such errors. Several previous studies examined genome-wide syntenic conservation to infer the contents of ancestral chromosomes and provided insights into the architecture of ancestral genomes, but did not provide methods or tools applicable to the study of the evolution of individual gene families. We developed an automated system to identify conserved syntenic regions in a primary genome using as outgroup a genome that diverged from the investigated lineage before a whole-genome duplication event. The product of this automated analysis, the Synteny Database, allows a user to examine fully or partially assembled genomes. The Synteny Database is optimized for the investigation of individual gene families in multiple lineages and can detect chromosomal inversions and translocations as well as ohnologs (paralogs derived by whole-genome duplication) gone missing. To demonstrate the utility of the system, we present a case study of gene family evolution, investigating the ARNTL gene family in the genomes of Ciona intestinalis, amphioxus, zebrafish, and human.

The spotted gar genome illuminates vertebrate evolution and facilitates human-teleost comparisons

To connect human biology to fish biomedical models, we sequenced the genome of spotted gar (Lepisosteus oculatus), whose lineage diverged from teleosts before teleost genome duplication (TGD). The slowly evolving gar genome has conserved in content and size many entire chromosomes from bony vertebrate ancestors. Gar bridges teleosts to tetrapods by illuminating the evolution of immunity, mineralization and development (mediated, for example, by Hox, ParaHox and microRNA genes). Numerous conserved noncoding elements (CNEs; often cis regulatory) undetectable in direct human-teleost comparisons become apparent using gar: functional studies uncovered conserved roles for such cryptic CNEs, facilitating annotation of sequences identified in human genome-wide association studies. Transcriptomic analyses showed that the sums of expression domains and expression levels for duplicated teleost genes often approximate the patterns and levels of expression for gar genes, consistent with subfunctionalization. The gar genome provides a resource for understanding evolution after genome duplication, the origin of vertebrate genomes and the function of human regulatory sequences.