Angel Amores

Stacks: an analysis tool set for population genomics

Massively parallel short‐read sequencing technologies, coupled with powerful software platforms, are enabling investigators to analyse tens of thousands of genetic markers. This wealth of data is rapidly expanding and allowing biological questions to be addressed with unprecedented scope and precision. The sizes of the data sets are now posing significant data processing and analysis challenges. Here we describe an extension of the Stacks software package to efficiently use genotype‐by‐sequencing data for studies of populations of organisms. Stacks now produces core population genomic summary statistics and SNP‐by‐SNP statistical tests. These statistics can be analysed across a reference genome using a smoothed sliding window. Stacks also now provides several output formats for several commonly used downstream analysis packages. The expanded population genomics functions in Stacks will make it a useful tool to harness the newest generation of massively parallel genotyping data for ecological and evolutionary genetics.

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.

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.

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.

The genome of the platyfish, Xiphophorus maculatus , provides insights into evolutionary adaptation and several complex traits

Several attributes intuitively considered to be typical mammalian features, such as complex behavior, live birth and malignant disease such as cancer, also appeared several times independently in lower vertebrates. The genetic mechanisms underlying the evolution of these elaborate traits are poorly understood. The platyfish, X. maculatus, offers a unique model to better understand the molecular biology of such traits. We report here the sequencing of the platyfish genome. Integrating genome assembly with extensive genetic maps identified an unexpected evolutionary stability of chromosomes in fish, in contrast to in mammals. Genes associated with viviparity show signatures of positive selection, identifying new putative functional domains and rare cases of parallel evolution. We also find that genes implicated in cognition show an unexpectedly high rate of duplicate gene retention after the teleost genome duplication event, suggesting a hypothesis for the evolution of the behavioral complexity in fish, which exceeds that found in amphibians and reptiles.

Multiple Sex-Associated Regions and a Putative Sex Chromosome in Zebrafish Revealed by RAD Mapping and Population Genomics

Within vertebrates, major sex determining genes can differ among taxa and even within species. In zebrafish (Danio rerio), neither heteromorphic sex chromosomes nor single sex determination genes of large effect, like Sry in mammals, have yet been identified. Furthermore, environmental factors can influence zebrafish sex determination. Although progress has been made in understanding zebrafish gonad differentiation (e.g. the influence of germ cells on gonad fate), the primary genetic basis of zebrafish sex determination remains poorly understood. To identify genetic loci associated with sex, we analyzed F2 offspring of reciprocal crosses between Oregon *AB and Nadia (NA) wild-type zebrafish stocks. Genome-wide linkage analysis, using more than 5,000 sequence-based polymorphic restriction site associated (RAD-tag) markers and population genomic analysis of more than 30,000 single nucleotide polymorphisms in our *ABxNA crosses revealed a sex-associated locus on the end of the long arm of chr-4 for both cross families, and an additional locus in the middle of chr-3 in one cross family. Additional sequencing showed that two SNPs in dmrt1 previously suggested to be functional candidates for sex determination in a cross of ABxIndia wild-type zebrafish, are not associated with sex in our AB fish. Our data show that sex determination in zebrafish is polygenic and that different genes may influence sex determination in different strains or that different genes become more important under different environmental conditions. The association of the end of chr-4 with sex is remarkable because, unique in the karyotype, this chromosome arm shares features with known sex chromosomes: it is highly heterochromatic, repetitive, late replicating, and has reduced recombination. Our results reveal that chr-4 has functional and structural properties expected of a sex chromosome.

Genome duplication, subfunction partitioning, and lineage divergence: Sox9 in stickleback and zebrafish

Teleosts are the most species‐rich group of vertebrates, and a genome duplication (tetraploidization) event in ray‐fin fish appears to have preceded this remarkable explosion of biodiversity. What is the relationship of the ray‐fin genome duplication to the teleost radiation? Genome duplication may have facilitated lineage divergence by partitioning different ancestral gene subfunctions among co‐orthologs of tetrapod genes in different teleost lineages. To test this hypothesis, we investigated gene expression patterns for Sox9 gene duplicates in stickleback and zebrafish, teleosts whose lineages diverged early in Euteleost evolution. Most expression domains appear to have been partitioned between Sox9a and Sox9b before the divergence of stickleback and zebrafish lineages, but some ancestral expression domains were distributed differentially in each lineage. We conclude that some gene subfunctions, as represented by lineage‐specific expression domains, may have assorted differently in separate lineages and that these may have contributed to lineage diversification during teleost evolution. Developmental Dynamics, 2003.

The zebrafish genome.

Publisher Summary Analysis of the zebrafish gene map suggests that during the course of vertebrate evolution, the zebrafish and human lineages have shared two rounds of whole genome duplication, and a third whole genome duplication event probably occurred before the teleost radiation. For many genes in humans, there may be two copies in zebrafish, and for many human chromosome segments, zebrafish is likely to have two such segments. The haploid zebrafish genome has 25 chromosomes, most of which are difficult to distinguish. These chromosomes contain about 1.7 X l0 9 base pairs of DNA, about half the mammalian genome size. The first genetic map of the zebrafish included about 400 genetic markers, mostly random amplified polymorphic DNAs (RAPDs), along with a few genes and mutations. The conservation of genome organization between humans and zebrafish can facilitate the transferal of information from one species to the other. The inference that zebrafish lineage appeared to have suffered an extra polyploidization step not experienced by mammals simply means that both zebrafish copies of a mammalian gene should be studied for better connectivity of molecular genetic mechanisms of development in zebrafish to the tetrapod—especially human—condition.

Wild Sex in Zebrafish: Loss of the Natural Sex Determinant in Domesticated Strains

Sex determination can be robustly genetic, strongly environmental, or genetic subject to environmental perturbation. The genetic basis of sex determination is unknown for zebrafish (Danio rerio), a model for development and human health. We used RAD-tag population genomics to identify sex-linked polymorphisms. After verifying this “RAD-sex” method on medaka (Oryzias latipes), we studied two domesticated zebrafish strains (AB and TU), two natural laboratory strains (WIK and EKW), and two recent isolates from nature (NA and CB). All four natural strains had a single sex-linked region at the right tip of chromosome 4, enabling sex genotyping by PCR. Genotypes for the single nucleotide polymorphism (SNP) with the strongest statistical association to sex suggested that wild zebrafish have WZ/ZZ sex chromosomes. In natural strains, “male genotypes” became males and some “female genotypes” also became males, suggesting that the environment or genetic background can cause female-to-male sex reversal. Surprisingly, TU and AB lacked detectable sex-linked loci. Phylogenomics rooted on D. nigrofasciatus verified that all strains are monophyletic. Because AB and TU branched as a monophyletic clade, we could not rule out shared loss of the wild sex locus in a common ancestor despite their independent domestication. Mitochondrial DNA sequences showed that investigated strains represent only one of the three identified zebrafish haplogroups. Results suggest that zebrafish in nature possess a WZ/ZZ sex-determination mechanism with a major determinant lying near the right telomere of chromosome 4 that was modified during domestication. Strains providing the zebrafish reference genome lack key components of the natural sex-determination system but may have evolved variant sex-determining mechanisms during two decades in laboratory culture.