Extensive linkage disequilibrium and parallel adaptive divergence across threespine stickleback genomes

Population genomic studies are beginning to provide a more comprehensive view of dynamic genome-scale processes in evolution. Patterns of genomic architecture, such as genomic islands of increased divergence, may be important for adaptive population differentiation and speciation. We used next-generation sequencing data to examine the patterns of local and long-distance linkage disequilibrium (LD) across oceanic and freshwater populations of threespine stickleback, a useful model for studies of evolution and speciation. We looked for associations between LD and signatures of divergent selection, and assessed the role of recombination rate variation in generating LD patterns. As predicted under the traditional biogeographic model of unidirectional gene flow from ancestral oceanic to derived freshwater stickleback populations, we found extensive local and long-distance LD in fresh water. Surprisingly, oceanic populations showed similar patterns of elevated LD, notably between large genomic regions previously implicated in adaptation to fresh water. These results support an alternative biogeographic model for the stickleback radiation, one of a metapopulation with appreciable bi-directional gene flow combined with strong divergent selection between oceanic and freshwater populations. As predicted by theory, these processes can maintain LD within and among genomic islands of divergence. These findings suggest that the genomic architecture in oceanic stickleback populations may provide a mechanism for the rapid re-assembly and evolution of multi-locus genotypes in newly colonized freshwater habitats, and may help explain genetic mapping of parallel phenotypic variation to similar loci across independent freshwater populations.     

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 F(2) 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.     

A publish-subscribe model of genetic networks

We present a simple model of genetic regulatory networks in which regulatory connections among genes are mediated by a limited number of signaling molecules. Each gene in our model produces (publishes) a single gene product, which regulates the expression of other genes by binding to regulatory regions that correspond (subscribe) to that product. We explore the consequences of this publish-subscribe model of regulation for the properties of single networks and for the evolution of populations of networks. Degree distributions of randomly constructed networks, particularly multimodal in-degree distributions, which depend on the length of the regulatory sequences and the number of possible gene products, differed from simpler Boolean NK models. In simulated evolution of populations of networks, single mutations in regulatory or coding regions resulted in multiple changes in regulatory connections among genes, or alternatively in neutral change that had no effect on phenotype. This resulted in remarkable evolvability in both number and length of attractors, leading to evolved networks far beyond the expectation of these measures based on random distributions. Surprisingly, this rapid evolution was not accompanied by changes in degree distribution; degree distribution in the evolved networks was not substantially different from that of randomly generated networks. The publish-subscribe model also allows exogenous gene products to create an environment, which may be noisy or stable, in which dynamic behavior occurs. In simulations, networks were able to evolve moderate levels of both mutational and environmental robustness.     

The ecology of an adaptive radiation of three‐spined stickleback from North Uist,Scotland

There has been a large focus on the genetics of traits involved in adaptation, but knowledge of the environmental variables leading to adaptive changes is surprisingly poor. Combined use of environmental data with morphological and genomic data should allow us to understand the extent to which patterns of phenotypic and genetic diversity within a species can be explained by the structure of the environment. Here, we analyse the variation of populations of three‐spined stickleback from 27 freshwater lakes on North Uist, Scotland, that vary greatly in their environment, to understand how environmental and genetic constraints contribute to phenotypic divergence. We collected 35 individuals per population and 30 abiotic and biotic environmental parameters to characterize variation across lakes and analyse phenotype–environment associations. Additionally, we used RAD sequencing to estimate the genetic relationships among a subset of these populations. We found a large amount of phenotypic variation among populations, most prominently in armour and spine traits. Despite large variation in the abiotic environment, namely in ion composition, depth and dissolved organic Carbon, more phenotypic variation was explained by the biotic variables (presence of predators and density of predator and competitors), than by associated abiotic variables. Genetic structure among populations was partly geographic, with closer populations being more similar. Altogether, our results suggest that differences in body shape among stickleback populations are the result of both canalized genetic and plastic responses to environmental factors, which shape fish morphology in a predictable direction regardless of their genetic starting point.     

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction‐site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small N_e (range = 2.1–89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome‐wide divergence. Nonetheless, outlier tests identified 3.6–6.6% of loci as high F_ST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high F_ST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small N_e in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.     

An essential saccharide binding domain for the mAb 2C7 established for Neisseria gonorrhoeae LOS by ES-MS and MSn

A study of bacterial surface oligosaccharides were investigated among different strains of Neisseria gonorrhoeae to correlate structural features essential for binding to the MAb 2C7. This epitope is widely expressed and conserved in gonococcal isolates, characteristics essential to an effective candidate vaccine antigen. Sample lipooligosaccharides (LOS), was prepared by a modification of the hot phenol-water method from which de-O-acetylated LOS and oligosaccharide (OS) components were analyzed by ES-MS-CID-MS and ES-MSnin a triple quadrupole and an ion trap mass spectrometer, respectively. Previously documented natural heterogeneity was apparent from both LOS and OS preparations which was admixed with fragments induced by hydrazine and mild acid treatment. Natural heterogeneity was limited to phosphorylation and antenni extensions to the alpha-chain. Mild acid hydrolysis to release OS also hydrolyzed the beta(1-->6) glycosidic linkage of lipid A. OS structures were determined by collisional and resonance excitation combined with MS and multistep MSn which provided sequence information from both neutral loss, and nonreducing terminal fragments. A comparison of OS structures, with earlier knowledge of MAb binding, enzyme treatment, and partial acid hydrolysis indicates a generic overlapping domain for 2C7 binding. Reoccurring structural features include a Hepalpha(1-->3)Hepbeta(1-->5)KDO trisaccharide core branched on the nonreducing terminus (Hep-2) with an alpha(1-->2) linked GlcNAc (gamma-chain), and an alpha-linked lactose (beta-chain) residue. From the central heptose (Hep-1), a beta(1-->4) linked lactose (alpha-chain), moiety is required although extensions to this residue appear unnecessary.      

Ecological genomics in full colour

Colour patterns in animals have long offered an opportunity to observe adaptive traits in natural populations. Colour plays myriad roles in interactions within and among species, from reproductive signalling to predator avoidance, leading to multiple targets of natural and sexual selection and opportunities for diversification. Understanding the genetic and developmental underpinnings of variation in colour promises a fuller understanding of these evolutionary processes, but the path to unravelling these connections can be arduous. The advent of genomic techniques suitable for nonmodel organisms is now beginning to light the way. Two new studies in this issue of Molecular Ecology use genomic sequencing of laboratory crosses to map colour traits in cichlid fishes, a remarkably diverse group in which coloration has played a major role in diversification. They illustrate how genomic approaches, specifically RAD sequencing, can rapidly identify both simple and more complex genetic variation underlying ecologically important traits. In the first, Henning et al. ( 2014 ) detect a single locus that appears to control in a Mendelian fashion the presence of horizontal stripes, a trait that has evolved in numerous cichlid lineages. In the second, Albertson et al. ( 2014 ) identify several genes and epistatic interactions affecting multiple colour traits, as well as a novel metric describing integration across colour traits. Albertson et al. ( 2014 ) go further, by quantifying differential expression of parental alleles at a candidate locus and by relating differentiation among natural populations at mapped loci to trait divergence. Herein lies the promise of ecological genomics – efficiently integrating genetic mapping of phenotypes with population genomic data to both identify functional genes and unravel their evolutionary history. These studies offer guidance on how genomic techniques can be tailored to a research question or study system, and they also add to the growing body of empirical examples addressing basic questions about how ecologically important traits evolve in natural populations.     

Conservation genomics of natural and managed populations: building a conceptual and practical framework

The boom of massive parallel sequencing (MPS) technology and its applications in conservation of natural and managed populations brings new opportunities and challenges to meet the scientific questions that can be addressed. Genomic conservation offers a wide range of approaches and analytical techniques, with their respective strengths and weaknesses that rely on several implicit assumptions. However, finding the most suitable approaches and analysis regarding our scientific question are often difficult and time‐consuming. To address this gap, a recent workshop entitled ‘ConGen 2015’ was held at Montana University in order to bring together the knowledge accumulated in this field and to provide training in conceptual and practical aspects of data analysis applied to the field of conservation and evolutionary genomics. Here, we summarize the expertise yield by each instructor that has led us to consider the importance of keeping in mind the scientific question from sampling to management practices along with the selection of appropriate genomics tools and bioinformatics challenges.     

RAD sequencing and genomic simulations resolve hybrid origins within North American Canis

Top predators are disappearing worldwide, significantly changing ecosystems that depend on top-down regulation. Conflict with humans remains the primary roadblock for large carnivore conservation, but for the eastern wolf (Canis lycaon), disagreement over its evolutionary origins presents a significant barrier to conservation in Canada and has impeded protection for grey wolves (Canis lupus) in the USA. Here, we use 127 235 single-nucleotide polymorphisms (SNPs) identified from restriction-site associated DNA sequencing (RAD-seq) of wolves and coyotes, in combination with genomic simulations, to test hypotheses of hybrid origins of Canis types in eastern North America. A principal components analysis revealed no evidence to support eastern wolves, or any other Canis type, as the product of grey wolf × western coyote hybridization. In contrast, simulations that included eastern wolves as a distinct taxon clarified the hybrid origins of Great Lakes-boreal wolves and eastern coyotes. Our results support the eastern wolf as a distinct genomic cluster in North America and help resolve hybrid origins of Great Lakes wolves and eastern coyotes. The data provide timely information that will shed new light on the debate over wolf conservation in eastern North America.     

Harnessing genomics for delineating conservation units

Genomic data have the potential to revolutionize the delineation of conservation units (CUs) by allowing the detection of adaptive genetic variation, which is otherwise difficult for rare, endangered species. In contrast to previous recommendations, we propose that the use of neutral versus adaptive markers should not be viewed as alternatives. Rather, neutral and adaptive markers provide different types of information that should be combined to make optimal management decisions. Genetic patterns at neutral markers reflect the interaction of gene flow and genetic drift that affects genome-wide variation within and among populations. This population genetic structure is what natural selection operates on to cause adaptive divergence. Here, we provide a new framework to integrate data on neutral and adaptive markers to protect biodiversity.