Leafy and weedy seadragon genomes connect genic and repetitive DNA features to the extravagant biology of syngnathid fishes.

Seadragons are a remarkable lineage of teleost fishes in the family Syngnathidae, renowned for having evolved male pregnancy. Comprising three known species, seadragons are widely recognized and admired for their fantastical body forms and coloration, and their specific habitat requirements have made them flagship representatives for marine conservation and natural history interests. Until recently, a gap has been the lack of significant genomic resources for seadragons. We have produced gene-annotated, chromosome-scale genome models for the leafy and weedy seadragon to advance investigations of evolutionary innovation and elaboration of morphological traits in seadragons as well as their pipefish and seahorse relatives. We identified several interesting features specific to seadragon genomes, including divergent noncoding regions near a developmental gene important for integumentary outgrowth, a high genome-wide density of repetitive DNA, and recent expansions of transposable elements and a vesicular trafficking gene family. Surprisingly, comparative analyses leveraging the seadragon genomes and additional syngnathid and outgroup genomes revealed striking, syngnathid-specific losses in the family of fibroblast growth factors (FGFs), which likely involve reorganization of highly conserved gene regulatory networks in ways that have not previously been documented in natural populations. The resources presented here serve as important tools for future evolutionary studies of developmental processes in syngnathids and hold value for conservation of the extravagant seadragons and their relatives.

Genome-wide analysis facilitates estimation of the amount of male contribution in meiotic gynogenetic three-spined stickleback (Gasterosteus aculeatus).

Gynogenetic embryos - those inheriting only maternal DNA - can be experimentally created by fertilizing eggs with radiation-treated sperm containing inactivated paternal chromosomes. Diploidy in the zygotes can be maintained through prevention of the second meiosis or restored by preventing the first mitosis after the maternal chromosome complement has been replicated. These gynogenetic organisms are useful in many fields including aquaculture, evolutionary biology and genomics. Although gynogenetic organisms have been created in numerous species, the completeness of uni-parental inheritance has often been assumed rather than thoroughly quantified across the genome. Instead, when tests of uni-parental inheritance occur, they typically rely on well-studied genetically determined phenotypes that represent a very small sub-set of the genome. Only assessing small genomic regions for paternal inheritance leaves the question of whether some paternal contributions to offspring might still have occurred. In this study, the authors quantify the efficacy of creating gynogenetic diploid three-spined stickleback fish (Gasterosteus aculeatus). To this end, the authors mirrored previous assessments of paternal contribution using well-studied genetically determined phenotypes including sex and genetically dominant morphological traits but expanded on previous studies using dense restriction site-associated DNA sequencing (RAD-seq) markers in parents and offspring to assess paternal inheritance genome-wide. In the gynogenetic diploids, the authors found no male genotypes underlying their phenotypes of interest - sex and dominant phenotypic traits. Using genome-wide assessments of paternal contribution, nevertheless, the authors found evidence of a small, yet potentially important, amount of paternally "leaked" genetic material. The application of this genome-wide approach identifies the need for more widespread assessment of paternal contributions to gynogenetic animals and promises benefits for many aspects of aquaculture, evolutionary biology and genomics.

Genome-wide single nucleotide polymorphism markers reveal population structure and dispersal direction of an expanding nuisance algal bloom species.

Species invasion and range expansion are currently under scrutiny due to increasing anthropogenic impact on the natural environment. This is also true for harmful algal blooms, which have been reported to have increased in frequency. However, this research is challenging due to the ephemeral nature, small size and mostly low concentrations of microalgae in the environment. One such species is the nuisance microalga Gonyostomum semen (Raphidophyceae), which has increased in occurrence in northern Europe in recent decades. The question of whether the species has expanded its habitat range or if it was already present in the lakes but was too rare to be detected remains unanswered. The aim of the present study was to determine the genetic structure and dispersal pathways of G. semen using RAD (restriction-site-associated DNA) tag sequencing. For G. semen, which has a huge genome (32 Gbp), we faced particular challenges, but were nevertheless able to recover over 1000 single nucleotide polymorphisms at high coverage. Our data revealed a distinct population genetic structure, demonstrating a divide of western and eastern populations that probably represent different lineages. Despite significant genetic differentiation among lakes, we found only limited isolation-by-distance. While we had expected a pattern of recent expansion northwards, the data demonstrated gene flow from the northeast/east towards the southwest/west. This genetic signature suggests that the observed gene flow may be due to dispersal by autumn migratory birds, which act as dispersal vectors of resistant resting propagules that form at the end of the G. semen blooms.

Evolution of stickleback in 50 years on earthquake-uplifted islands.

How rapidly can animal populations in the wild evolve when faced with sudden environmental shifts? Uplift during the 1964 Great Alaska Earthquake abruptly created freshwater ponds on multiple islands in Prince William Sound and the Gulf of Alaska. In the short time since the earthquake, the phenotypes of resident freshwater threespine stickleback fish on at least three of these islands have changed dramatically from their oceanic ancestors. To test the hypothesis that these freshwater populations were derived from oceanic ancestors only 50 y ago, we generated over 130,000 single-nucleotide polymorphism genotypes from more than 1,000 individuals using restriction site-associated DNA sequencing (RAD-seq). Population genomic analyses of these data support the hypothesis of recent and repeated, independent colonization of freshwater habitats by oceanic ancestors. We find evidence of recurrent gene flow between oceanic and freshwater ecotypes where they co-occur. Our data implicate natural selection in phenotypic diversification and support the hypothesis that the metapopulation organization of this species helps maintain a large pool of genetic variation that can be redeployed rapidly when oceanic stickleback colonize freshwater environments. We find that the freshwater populations, despite population genetic analyses clearly supporting their young age, have diverged phenotypically from oceanic ancestors to nearly the same extent as populations that were likely founded thousands of years ago. Our results support the intriguing hypothesis that most stickleback evolution in fresh water occurs within the first few decades after invasion of a novel environment.

Developmental timing differences underlie armor loss across threespine stickleback populations.

Comparing ontogenetic patterns within a well-described evolutionary context aids in inferring mechanisms of change, including heterochronies or deletion of developmental pathways. Because selection acts on phenotypes throughout ontogeny, any within-taxon developmental variation has implications for evolvability. We compare ontogenetic order and timing of locomotion and defensive traits in three populations of threespine stickleback that have evolutionarily divergent adult forms. This analysis adds to the growing understanding of developmental genetic mechanisms of adaptive change in this evolutionary model species by delineating when chondrogenesis and osteogenesis in two derived populations begin to deviate from the developmental pattern in their immediate ancestors. We found that differences in adult defensive morphologies arise through abolished or delayed initiation of these traits rather than via an overall heterochronic shift, that intra-population ontogenetic variation is increased for some derived traits, and that altered armor developmental timing differentiates the derived populations from each other despite parallels in adult lateral plate armor phenotypes. We found that changes in ossified elements of the pelvic armor are linked to delayed and incomplete development of an early-forming pelvic cartilage, and that this disruption likely presages the variable pelvic vestiges documented in many derived populations.

Remembering James Alan Bassham (1922-2012).

James Alan Bassham, known to many as Al, was born on November 26, 1922, in Sacramento, California (CA), USA. He died on November 19, 2012, in El Cerrito, CA. To celebrate his life at his 3rd death anniversary, we present here a brief biography, comments on his discoveries, but most importantly, remembrances from family and friends; we remember this wonderful and modest person who had played a major pivotal role in the discoveries that led to what he would like to call the P(hotosynthetic) C(arbon) R(eduction) cycle, known to many as the Calvin Cycle, the Calvin-Benson Cycle, or the Calvin-Benson-Bassham Cycle. Based on a personal request by Bassham himself to one of us (Govindjee), we refrain from including his name in the cycle-in recognition of his many students and associates he would have liked to honor.

Effects of crossovers between homeologs on inheritance and population genomics in polyploid-derived salmonid fishes.

A whole genome duplication occurred in the ancestor of all salmonid fishes some 50-100 million years ago. Early inheritance studies with allozymes indicated that loci in the salmonid genome are inherited disomically in females. However, some pairs of duplicated loci showed patterns of inheritance in males indicating pairing and recombination between homeologous chromosomes. Nearly 20% of loci in the salmonid genome are duplicated and share the same alleles (isoloci), apparently due to homeologous recombination. Half-tetrad analysis revealed that isoloci tend to be telomeric. These results suggested that residual tetrasomic inheritance of isoloci results from homeologous recombination near chromosome ends and that continued disomic inheritance resulted from homologous pairing of centromeric regions. Many current genetic maps of salmonids are based on single nucleotide polymorphisms and microsatellites that are no longer duplicated. Therefore, long sections of chromosomes on these maps are poorly represented, especially telomeric regions. In addition, preferential multivalent pairing of homeologs from the same species in F1 hybrids results in an excess of nonparental gametes (so-called pseudolinkage). We consider how not including duplicated loci has affected our understanding of population and evolutionary genetics of salmonids, and we discuss how incorporating these loci will benefit our understanding of population genomics.

Single Cell Sequencing Provides Clues about the Developmental Genetic Basis of Evolutionary Adaptations in Syngnathid Fishes.

Seahorses, pipefishes, and seadragons are fishes from the family Syngnathidae that have evolved extraordinary traits including male pregnancy, elongated snouts, loss of teeth, and dermal bony armor. The developmental genetic and cellular changes that led to the evolution of these traits are largely unknown. Recent syngnathid genome assemblies revealed suggestive gene content differences and provide the opportunity for detailed genetic analyses. We created a single cell RNA sequencing atlas of Gulf pipefish embryos to understand the developmental basis of four traits: derived head shape, toothlessness, dermal armor, and male pregnancy. We completed marker gene analyses, built genetic networks, and examined spatial expression of select genes. We identified osteochondrogenic mesenchymal cells in the elongating face that express regulatory genes bmp4, sfrp1a, and prdm16. We found no evidence for tooth primordia cells, and we observed re-deployment of osteoblast genetic networks in developing dermal armor. Finally, we found that epidermal cells expressed nutrient processing and environmental sensing genes, potentially relevant for the brooding environment. The examined pipefish evolutionary innovations are composed of recognizable cell types, suggesting derived features originate from changes within existing gene networks. Future work addressing syngnathid gene networks across multiple stages and species is essential for understanding how their novelties evolved.

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.

The population structure and recent colonization history of Oregon threespine stickleback determined using restriction-site associated DNA-sequencing.

Understanding how genetic variation is partitioned across genomes within and among populations is a fundamental problem in ecological and evolutionary genetics. To address this problem, we studied the threespine stickleback fish, which has repeatedly undergone parallel phenotypic and genetic differentiation when oceanic fish have invaded freshwater habitats. While significant evolutionary genetic research has been performed using stickleback from geographic regions that have been deglaciated in the last 20 000 years, less research has focused on freshwater populations that predate the last glacial maximum. We performed restriction-site associated DNA-sequencing (RAD-seq) based population genomic analyses on stickleback from across Oregon, which was not glaciated during the last maximum. We sampled stickleback from coastal, Willamette Basin and central Oregon sites, analysed their genetic diversity using RAD-seq, performed structure analyses, reconstructed their phylogeographic history and tested the hypothesis of recent stickleback introduction into central Oregon, where incidence of this species was only recently documented. Our results showed a clear phylogeographic break between coastal and inland populations, with oceanic populations exhibiting the lowest levels of divergence from one another. Willamette Basin and central Oregon populations formed a clade of closely related populations, a finding consistent with a recent introduction of stickleback into central Oregon. Finally, genome-wide analysis of genetic diversity (π) and correlations of alleles within individuals in subpopulations (FIS) supported a role for introgressive hybridization in coastal populations and a recent expansion in central Oregon. Our results exhibit the power of next-generation sequencing genomic approaches such as RAD-seq to identify both historical population structure and recent colonization history.

Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags.

Next-generation sequencing technology provides novel opportunities for gathering genome-scale sequence data in natural populations, laying the empirical foundation for the evolving field of population genomics. Here we conducted a genome scan of nucleotide diversity and differentiation in natural populations of threespine stickleback (Gasterosteus aculeatus). We used Illumina-sequenced RAD tags to identify and type over 45,000 single nucleotide polymorphisms (SNPs) in each of 100 individuals from two oceanic and three freshwater populations. Overall estimates of genetic diversity and differentiation among populations confirm the biogeographic hypothesis that large panmictic oceanic populations have repeatedly given rise to phenotypically divergent freshwater populations. Genomic regions exhibiting signatures of both balancing and divergent selection were remarkably consistent across multiple, independently derived populations, indicating that replicate parallel phenotypic evolution in stickleback may be occurring through extensive, parallel genetic evolution at a genome-wide scale. Some of these genomic regions co-localize with previously identified QTL for stickleback phenotypic variation identified using laboratory mapping crosses. In addition, we have identified several novel regions showing parallel differentiation across independent populations. Annotation of these regions revealed numerous genes that are candidates for stickleback phenotypic evolution and will form the basis of future genetic analyses in this and other organisms. This study represents the first high-density SNP-based genome scan of genetic diversity and differentiation for populations of threespine stickleback in the wild. These data illustrate the complementary nature of laboratory crosses and population genomic scans by confirming the adaptive significance of previously identified genomic regions, elucidating the particular evolutionary and demographic history of such regions in natural populations, and identifying new genomic regions and candidate genes of evolutionary significance.

Improvements to the Gulf pipefish Syngnathus scovelli genome.

The Gulf pipefish Syngnathus scovelli has emerged as an important species for studying sexual selection, development, and physiology. Comparative evolutionary genomics research involving fishes from Syngnathidae depends on having a high-quality genome assembly and annotation. However, the first S. scovelli genome assembled using short-read sequences and a smaller RNA-sequence dataset has limited contiguity and a relatively poor annotation. Here, using PacBio long-read high-fidelity sequences and a proximity ligation library, we generate an improved assembly to obtain 22 chromosome-level scaffolds. Compared to the first assembly, the gaps in the improved assembly are smaller, the N75 is larger, and our genome is ~95% BUSCO complete. Using a large body of RNA-Seq reads from different tissue types and NCBI's Eukaryotic Annotation Pipeline, we discovered 28,162 genes, of which 8,061 are non-coding genes. Our new genome assembly and annotation are tagged as a RefSeq genome by NCBI and provide enhanced resources for research work involving S. scovelli..

Host genomic variation shapes gut microbiome diversity in threespine stickleback fish.

IMPORTANCE: A major focus of host-microbe research is to understand how genetic differences, of various magnitudes, among hosts translate to differences in their microbiomes. This has been challenging for animal hosts, including humans, because it is difficult to control environmental variables tightly enough to isolate direct genetic effects on the microbiome. Our work in stickleback fish is a significant contribution because our experimental approach allowed strict control over environmental factors, including standardization of the microbiome from the earliest stage of development and unrestricted co-housing of fish in a truly common environment. Furthermore, we measured host genetic variation over 2,000 regions of the stickleback genome, comparing this information and microbiome composition data among fish from very similar and very different genetic backgrounds. Our findings highlight how differences in the host genome influence microbiome diversity and make a case for future manipulative microbiome experiments that use host systems with naturally occurring genetic variation.

Extreme intraspecific divergence in mitochondrial haplotypes makes the threespine stickleback fish an emerging evolutionary mutant model for mito-nuclear interactions.

Mitochondrial DNA is primarily maternally inherited in most animals and evolves about 10 times faster than biparentally inherited nuclear DNA. Mitochondrial dysfunction (mt-dys) arises when interactions between the co-evolving mitochondrial and nuclear genomes are perturbed in essential processes like oxidative phosphorylation (OXPHOS). Over time mt-dys can lead to mitochondrial diseases (mt-diseases), which are surprisingly prevalent and include common diseases such as Alzheimer's, Parkinson's, and diabetes. Unfortunately, the strong impact that intraspecific mitochondrial and nuclear genetic variation has on mt-disease complicates its study and the development of effective treatments. Animal models have advanced our understanding of mt-disease but their relevance to human conditions is often limited by their relatively low nuclear genetic diversity. Many traditional laboratory models also typically have a single mitochondrial haplotype (mitotype), in stark contrast to over 5,000 mitotypes in humans worldwide. The threespine stickleback fish has an evolutionary history that has made it a favorable evolutionary mutant model (EMM) for studying mito-nuclear interactions and possibly mt-diseases. EMMs are species with naturally evolved states that mimic maladaptive human diseases. In threespine stickleback, a period of isolation followed by introgression of the mitochondrial genome from a sister species resulted in the maintenance of two distinct mitochondrial haplotypes which continue to segregate within many populations of wild stickleback. The existence of two mitogenomes segregating in numerous genetically diverse populations provides a unique system for exploring complex mito-nuclear dynamics. Here we provide the first complete coding region analysis of the two threespine stickleback mitotypes, whose mitogenomic divergence exceeds that of other mammalian models for mitochondrial disease and even that between ancient and modern humans. We find that divergence is not uniform across the mitogenome, but primarily impacts protein coding genes, and significantly impacts proteins in Complex I of OXPHOS. The full characterization of these highly divergent intraspecific mitotypes provides a foundation for the development of threespine stickleback as an EMM for mito-nuclear interactions.

Single-cell Iso-Sequencing enables rapid genome annotation for scRNAseq analysis.

Single-cell RNA sequencing is a powerful technique that continues to expand across various biological applications. However, incomplete 3'-UTR annotations can impede single-cell analysis resulting in genes that are partially or completely uncounted. Performing single-cell RNA sequencing with incomplete 3'-UTR annotations can hinder the identification of cell identities and gene expression patterns and lead to erroneous biological inferences. We demonstrate that performing single-cell isoform sequencing in tandem with single-cell RNA sequencing can rapidly improve 3'-UTR annotations. Using threespine stickleback fish (Gasterosteus aculeatus), we show that gene models resulting from a minimal embryonic single-cell isoform sequencing dataset retained 26.1% greater single-cell RNA sequencing reads than gene models from Ensembl alone. Furthermore, pooling our single-cell sequencing isoforms with a previously published adult bulk Iso-Seq dataset from stickleback, and merging the annotation with the Ensembl gene models, resulted in a marginal improvement (+0.8%) over the single-cell isoform sequencing only dataset. In addition, isoforms identified by single-cell isoform sequencing included thousands of new splicing variants. The improved gene models obtained using single-cell isoform sequencing led to successful identification of cell types and increased the reads identified of many genes in our single-cell RNA sequencing stickleback dataset. Our work illuminates single-cell isoform sequencing as a cost-effective and efficient mechanism to rapidly annotate genomes for single-cell RNA sequencing.

Chromonomer: A Tool Set for Repairing and Enhancing Assembled Genomes Through Integration of Genetic Maps and Conserved Synteny.

The pace of the sequencing and computational assembly of novel reference genomes is accelerating. Though DNA sequencing technologies and assembly software tools continue to improve, biological features of genomes such as repetitive sequence as well as molecular artifacts that often accompany sequencing library preparation can lead to fragmented or chimeric assemblies. If left uncorrected, defects like these trammel progress on understanding genome structure and function, or worse, positively mislead this research. Fortunately, integration of additional, independent streams of information, such as a marker-dense genetic map and conserved orthologous gene order from related taxa, can be used to scaffold together unlinked, disordered fragments and to restructure a reference genome where it is incorrectly joined. We present a tool set for automating these processes, one that additionally tracks any changes to the assembly and to the genetic map, and which allows the user to scrutinize these changes with the help of web-based, graphical visualizations. Chromonomer takes a user-defined reference genome, a map of genetic markers, and, optionally, conserved synteny information to construct an improved reference genome of chromosome models: a "chromonome". We demonstrate Chromonomer's performance on genome assemblies and genetic maps that have disparate characteristics and levels of quality.

Evolution of developmental roles of Pax2/5/8 paralogs after independent duplication in urochordate and vertebrate lineages.

BACKGROUND: Gene duplication provides opportunities for lineage diversification and evolution of developmental novelties. Duplicated genes generally either disappear by accumulation of mutations (nonfunctionalization), or are preserved either by the origin of positively selected functions in one or both duplicates (neofunctionalization), or by the partitioning of original gene subfunctions between the duplicates (subfunctionalization). The Pax2/5/8 family of important developmental regulators has undergone parallel expansion among chordate groups. After the divergence of urochordate and vertebrate lineages, two rounds of independent gene duplications resulted in the Pax2, Pax5, and Pax8 genes of most vertebrates (the sister group of the urochordates), and an additional duplication provided the pax2a and pax2b duplicates in teleost fish. Separate from the vertebrate genome expansions, a duplication also created two Pax2/5/8 genes in the common ancestor of ascidian and larvacean urochordates. RESULTS: To better understand mechanisms underlying the evolution of duplicated genes, we investigated, in the larvacean urochordate Oikopleura dioica, the embryonic gene expression patterns of Pax2/5/8 paralogs. We compared the larvacean and ascidian expression patterns to infer modular subfunctions present in the single pre-duplication Pax2/5/8 gene of stem urochordates, and we compared vertebrate and urochordate expression to infer the suite of Pax2/5/8 gene subfunctions in the common ancestor of olfactores (vertebrates + urochordates). Expression pattern differences of larvacean and ascidian Pax2/5/8 orthologs in the endostyle, pharynx and hindgut suggest that some ancestral gene functions have been partitioned differently to the duplicates in the two urochordate lineages. Novel expression in the larvacean heart may have resulted from the neofunctionalization of a Pax2/5/8 gene in the urochordates. Expression of larvacean Pax2/5/8 in the endostyle, in sites of epithelial remodeling, and in sensory tissues evokes like functions of Pax2, Pax5 and Pax8 in vertebrate embryos, and may indicate ancient origins for these functions in the chordate common ancestor. CONCLUSION: Comparative analysis of expression patterns of chordate Pax2/5/8 duplicates, rooted on the single-copy Pax2/5/8 gene of amphioxus, whose lineage diverged basally among chordates, provides new insights into the evolution and development of the heart, thyroid, pharynx, stomodeum and placodes in chordates; supports the controversial conclusion that the atrial siphon of ascidians and the otic placode in vertebrates are homologous; and backs the notion that Pax2/5/8 functioned in ancestral chordates to engineer epithelial fusions and perforations, including gill slit openings.

Genetic divergence outpaces phenotypic evolution among threespine stickleback populations in old freshwater habitats.

Species such as threespine stickleback (Gasterosteus aculeatus) that inhabit divergent selective environments and that have diversified on different time scales can be of value for understanding evolutionary processes. Here we synthesize high-resolution genotypic and phenotypic data to explore a largely unstudied distribution of threespine stickleback populations living in oceanic and freshwater habitats along coastal and inland regions of Oregon. Many inland aquatic habitats of Oregon remained unglaciated during the last ice age, meaning that some extant Oregon lake and river stickleback may have descended from freshwater populations established long before more well-studied, post-glacial freshwater populations. To address the degree of congruence between genetic and phenotypic divergence, we directly compared Oregon stickleback to much younger (post-glacial) Alaskan populations. We found phenotypic variation in Oregon stickleback to be primarily partitioned between oceanic and freshwater habitats, as has been documented in other stickleback systems. However, the main axis of genetic divergence was between coastal and inland regions regardless of habitat type. Furthermore, when comparing patterns between Oregon and Alaska we found similar levels of phenotypic divergence, but much greater genetic divergence among Oregon's populations. The Oregon stickleback system therefore appears well suited for future studies linking genotypic and phenotypic change, further extending the utility of this small fish to provide general insights into evolutionary processes.

Highly Reproducible 16S Sequencing Facilitates Measurement of Host Genetic Influences on the Stickleback Gut Microbiome.

Multicellular organisms interact with resident microbes in important ways, and a better understanding of host-microbe interactions is aided by tools such as high-throughput 16S sequencing. However, rigorous evaluation of the veracity of these tools in a different context from which they were developed has often lagged behind. Our goal was to perform one such critical test by examining how variation in tissue preparation and DNA isolation could affect inferences about gut microbiome variation between two genetically divergent lines of threespine stickleback fish maintained in the same laboratory environment. Using careful experimental design and intensive sampling of individuals, we addressed technical and biological sources of variation in 16S-based estimates of microbial diversity. After employing a two-tiered bead beating approach that comprised tissue homogenization followed by microbial lysis in subsamples, we found an extremely minor effect of DNA isolation protocol relative to among-host microbial diversity differences. Abundance estimates for rare operational taxonomic units (OTUs), however, showed much lower reproducibility. Gut microbiome composition was highly variable across fish-even among cohoused siblings-relative to technical replicates, but a subtle effect of host genotype (stickleback line) was nevertheless detected for some microbial taxa.IMPORTANCE Our findings demonstrate the importance of appropriately quantifying biological and technical variance components when attempting to understand major influences on high-throughput microbiome data. Our focus was on understanding among-host (biological) variance in community metrics and its magnitude in relation to within-host (technical) variance, because meaningful comparisons among individuals are necessary in addressing major questions in host-microbe ecology and evolution, such as heritability of the microbiome. Our study design and insights should provide a useful example for others desiring to quantify microbiome variation at biological levels in the face of various technical factors in a variety of systems.