A recombineering pipeline to clone large and complex genes in Chlamydomonas.

The ability to clone genes has greatly advanced cell and molecular biology research, enabling researchers to generate fluorescent protein fusions for localization and confirm genetic causation by mutant complementation. Most gene cloning is polymerase chain reaction (PCR)�or DNA synthesis-dependent, which can become costly and technically challenging as genes increase in size, particularly if they contain complex regions. This has been a long-standing challenge for the Chlamydomonas reinhardtii research community, as this alga has a high percentage of genes containing complex sequence structures. Here we overcame these challenges by developing a recombineering pipeline for the rapid parallel cloning of genes from a Chlamydomonas bacterial artificial chromosome collection. To generate fluorescent protein fusions for localization, we applied the pipeline at both batch and high-throughput scales to 203 genes related to the Chlamydomonas CO2 concentrating mechanism (CCM), with an overall cloning success rate of 77%. Cloning success was independent of gene size and complexity, with cloned genes as large as 23 kb. Localization of a subset of CCM targets confirmed previous mass spectrometry data, identified new pyrenoid components, and enabled complementation of mutants. We provide vectors and detailed protocols to facilitate easy adoption of this technology, which we envision will open up new possibilities in algal and plant research.

A novel terpene synthase controls differences in anti-aphrodisiac pheromone production between closely related Heliconius butterflies.

Plants and insects often use the same compounds for chemical communication, but not much is known about the genetics of convergent evolution of chemical signals. The terpene (E)-ß-ocimene is a common component of floral scent and is also used by the butterfly Heliconius melpomene as an anti-aphrodisiac pheromone. While the biosynthesis of terpenes has been described in plants and microorganisms, few terpene synthases (TPSs) have been identified in insects. Here, we study the recent divergence of 2 species, H. melpomene and Heliconius cydno, which differ in the presence of (E)-ß-ocimene; combining linkage mapping, gene expression, and functional analyses, we identify 2 novel TPSs. Furthermore, we demonstrate that one, HmelOS, is able to synthesise (E)-ß-ocimene in vitro. We find no evidence for TPS activity in HcydOS (HmelOS ortholog of H. cydno), suggesting that the loss of (E)-ß-ocimene in this species is the result of coding, not regulatory, differences. The TPS enzymes we discovered are unrelated to previously described plant and insect TPSs, demonstrating that chemical convergence has independent evolutionary origins.

Superoxide is promoted by sucrose and affects amplitude of circadian rhythms in the evening.

Plants must coordinate photosynthetic metabolism with the daily environment and adapt rhythmic physiology and development to match carbon availability. Circadian clocks drive biological rhythms which adjust to environmental cues. Products of photosynthetic metabolism, including sugars and reactive oxygen species (ROS), are closely associated with the plant circadian clock, and sugars have been shown to provide metabolic feedback to the circadian oscillator. Here, we report a comprehensive sugar-regulated transcriptome of Arabidopsis and identify genes associated with redox and ROS processes as a prominent feature of the transcriptional response. We show that sucrose increases levels of superoxide (O2-), which is required for transcriptional and growth responses to sugar. We identify circadian rhythms of O2--regulated transcripts which are phased around dusk and find that O2- is required for sucrose to promote expression of TIMING OF CAB1 (TOC1) in the evening. Our data reveal a role for O2- as a metabolic signal affecting transcriptional control of the circadian oscillator in Arabidopsis.

Recombination rate variation shapes barriers to introgression across butterfly genomes.

Hybridisation and introgression can dramatically alter the relationships among groups of species, leading to phylogenetic discordance across the genome and between populations. Introgression can also erode species differences over time, but selection against introgression at certain loci acts to maintain postmating species barriers. Theory predicts that species barriers made up of many loci throughout the genome should lead to a broad correlation between introgression and recombination rate, which determines the extent to which selection on deleterious foreign alleles will affect neutral alleles at physically linked loci. Here, we describe the variation in genealogical relationships across the genome among three species of Heliconius butterflies: H. melpomene (mel), H. cydno (cyd), and H. timareta (tim), using whole genomes of 92 individuals, and ask whether this variation can be explained by heterogeneous barriers to introgression. We find that species relationships vary predictably at the chromosomal scale. By quantifying recombination rate and admixture proportions, we then show that rates of introgression are predicted by variation in recombination rate. This implies that species barriers are highly polygenic, with selection acting against introgressed alleles across most of the genome. In addition, long chromosomes, which have lower recombination rates, produce stronger barriers on average than short chromosomes. Finally, we find a consistent difference between two species pairs on either side of the Andes, which suggests differences in the architecture of the species barriers. Our findings illustrate how the combined effects of hybridisation, recombination, and natural selection, acting at multitudes of loci over long periods, can dramatically sculpt the phylogenetic relationships among species.

Genetic dissection of assortative mating behavior.

The evolution of new species is made easier when traits under divergent ecological selection are also mating cues. Such ecological mating cues are now considered more common than previously thought, but we still know little about the genetic changes underlying their evolution or more generally about the genetic basis for assortative mating behaviors. Both tight physical linkage and the existence of large-effect preference loci will strengthen genetic associations between behavioral and ecological barriers, promoting the evolution of assortative mating. The warning patterns of Heliconius melpomene and H. cydno are under disruptive selection due to increased predation of nonmimetic hybrids and are used during mate recognition. We carried out a genome-wide quantitative trait locus (QTL) analysis of preference behaviors between these species and showed that divergent male preference has a simple genetic basis. We identify three QTLs that together explain a large proportion (approximately 60%) of the difference in preference behavior observed between the parental species. One of these QTLs is just 1.2 (0-4.8) centiMorgans (cM) from the major color pattern gene optix, and, individually, all three have a large effect on the preference phenotype. Genomic divergence between H. cydno and H. melpomene is high but broadly heterogenous, and admixture is reduced at the preference-optix color pattern locus but not the other preference QTLs. The simple genetic architecture we reveal will facilitate the evolution and maintenance of new species despite ongoing gene flow by coupling behavioral and ecological aspects of reproductive isolation.

Ancestral Admixture Is the Main Determinant of Global Biodiversity in Fission Yeast.

Mutation and recombination are key evolutionary processes governing phenotypic variation and reproductive isolation. We here demonstrate that biodiversity within all globally known strains of Schizosaccharomyces pombe arose through admixture between two divergent ancestral lineages. Initial hybridization was inferred to have occurred ∼20-60 sexual outcrossing generations ago consistent with recent, human-induced migration at the onset of intensified transcontinental trade. Species-wide heritable phenotypic variation was explained near-exclusively by strain-specific arrangements of alternating ancestry components with evidence for transgressive segregation. Reproductive compatibility between strains was likewise predicted by the degree of shared ancestry. To assess the genetic determinants of ancestry block distribution across the genome, we characterized the type, frequency, and position of structural genomic variation using nanopore and single-molecule real-time sequencing. Despite being associated with double-strand break initiation points, over 800 segregating structural variants exerted overall little influence on the introgression landscape or on reproductive compatibility between strains. In contrast, we found strong ancestry disequilibrium consistent with negative epistatic selection shaping genomic ancestry combinations during the course of hybridization. This study provides a detailed, experimentally tractable example that genomes of natural populations are mosaics reflecting different evolutionary histories. Exploiting genome-wide heterogeneity in the history of ancestral recombination and lineage-specific mutations sheds new light on the population history of S. pombe and highlights the importance of hybridization as a creative force in generating biodiversity.

How do predators generalize warning signals in simple and complex prey communities? Insights from a videogame.

The persistence of distinct warning signals within and between sympatric mimetic communities is a puzzling evolutionary question because selection favours convergence of colour patterns among toxic species. Such convergence is partly shaped by predators' reaction to similar but not identical stimulus (i.e. generalization behaviour), and generalization by predators is likely to be shaped by the diversity of local prey. However, studying generalization behaviour is generally limited to simple variations of prey colour patterns. Here, we used a computer game played by humans as surrogate predators to investigate generalization behaviours in simple (4 morphs) and complex (10 morphs) communities of unprofitable (associated with a penalty) and profitable butterflies. Colour patterns used in the game are observed in the natural populations of unprofitable butterfly species such as Heliconius numata. Analyses of 449 game participants' behaviours show that players avoided unprofitable prey more readily in simple than in complex communities. However, generalization was observed only in players that faced complex communities, enhancing the protection of profitable prey that looked similar to at least one unprofitable morph. Additionally, similarity among unprofitable prey also reduced attack rates only in complex communities. These results are consistent with previous studies using avian predators but artificial colour patterns and suggest that mimicry is more likely to evolve in complex communities where increases in similarity are more likely to be advantageous.

Sexually dimorphic gene expression and transcriptome evolution provide mixed evidence for a fast-Z effect in Heliconius.

Sex chromosomes have different evolutionary properties compared to autosomes due to their hemizygous nature. In particular, recessive mutations are more readily exposed to selection, which can lead to faster rates of molecular evolution. Here, we report patterns of gene expression and molecular evolution for a group of butterflies. First, we improve the completeness of the Heliconius melpomene reference annotation, a neotropical butterfly with a ZW sex determination system. Then, we analyse RNA from male and female whole abdomens and sequence female ovary and gut tissue to identify sex- and tissue-specific gene expression profiles in H. melpomene. Using these expression profiles, we compare (a) sequence divergence and polymorphism; (b) the strength of positive and negative selection; and (c) rates of adaptive evolution, for Z and autosomal genes between two species of Heliconius butterflies, H. melpomene and H. erato. We show that the rate of adaptive substitutions is higher for Z than autosomal genes, but contrary to expectation, it is also higher for male-biased than female-biased genes. Additionally, we find no significant increase in the rate of adaptive evolution or purifying selection on genes expressed in ovary tissue, a heterogametic-specific tissue. Our results contribute to a growing body of literature from other ZW systems that also provide mixed evidence for a fast-Z effect where hemizygosity influences the rate of adaptive substitutions.

Genome-wide genetic marker discovery and genotyping using next-generation sequencing.

The advent of next-generation sequencing (NGS) has revolutionized genomic and transcriptomic approaches to biology. These new sequencing tools are also valuable for the discovery, validation and assessment of genetic markers in populations. Here we review and discuss best practices for several NGS methods for genome-wide genetic marker development and genotyping that use restriction enzyme digestion of target genomes to reduce the complexity of the target. These new methods -- which include reduced-representation sequencing using reduced-representation libraries (RRLs) or complexity reduction of polymorphic sequences (CRoPS), restriction-site-associated DNA sequencing (RAD-seq) and low coverage genotyping -- are applicable to both model organisms with high-quality reference genome sequences and, excitingly, to non-model species with no existing genomic data.

Evaluating the use of ABBA-BABA statistics to locate introgressed loci.

Several methods have been proposed to test for introgression across genomes. One method tests for a genome-wide excess of shared derived alleles between taxa using Patterson's D statistic, but does not establish which loci show such an excess or whether the excess is due to introgression or ancestral population structure. Several recent studies have extended the use of D by applying the statistic to small genomic regions, rather than genome-wide. Here, we use simulations and whole-genome data from Heliconius butterflies to investigate the behavior of D in small genomic regions. We find that D is unreliable in this situation as it gives inflated values when effective population size is low, causing D outliers to cluster in genomic regions of reduced diversity. As an alternative, we propose a related statistic ƒ(d), a modified version of a statistic originally developed to estimate the genome-wide fraction of admixture. ƒ(d) is not subject to the same biases as D, and is better at identifying introgressed loci. Finally, we show that both D and ƒ(d) outliers tend to cluster in regions of low absolute divergence (d(XY)), which can confound a recently proposed test for differentiating introgression from shared ancestral variation at individual loci.

Estimation of the spontaneous mutation rate in Heliconius melpomene.

We estimated the spontaneous mutation rate in Heliconius melpomene by genome sequencing of a pair of parents and 30 of their offspring, based on the ratio of number of de novo heterozygotes to the number of callable site-individuals. We detected nine new mutations, each one affecting a single site in a single offspring. This yields an estimated mutation rate of 2.9 × 10(-9) (95% confidence interval, 1.3 × 10(-9)-5.5 × 10(-9)), which is similar to recent estimates in Drosophila melanogaster, the only other insect species in which the mutation rate has been directly estimated. We infer that recent effective population size of H. melpomene is about 2 million, a substantially lower value than its census size, suggesting a role for natural selection reducing diversity. We estimate that H. melpomene diverged from its Müllerian comimic H. erato about 6 Ma, a somewhat later date than estimates based on a local molecular clock.

A physiologically required G protein-coupled receptor (GPCR)-regulator of G protein signaling (RGS) interaction that compartmentalizes RGS activity.

G protein-coupled receptors (GPCRs) can interact with regulator of G protein signaling (RGS) proteins. However, the effects of such interactions on signal transduction and their physiological relevance have been largely undetermined. Ligand-bound GPCRs initiate by promoting exchange of GDP for GTP on the Gα subunit of heterotrimeric G proteins. Signaling is terminated by hydrolysis of GTP to GDP through intrinsic GTPase activity of the Gα subunit, a reaction catalyzed by RGS proteins. Using yeast as a tool to study GPCR signaling in isolation, we define an interaction between the cognate GPCR (Mam2) and RGS (Rgs1), mapping the interaction domains. This reaction tethers Rgs1 at the plasma membrane and is essential for physiological signaling response. In vivo quantitative data inform the development of a kinetic model of the GTPase cycle, which extends previous attempts by including GPCR-RGS interactions. In vivo and in silico data confirm that GPCR-RGS interactions can impose an additional layer of regulation through mediating RGS subcellular localization to compartmentalize RGS activity within a cell, thus highlighting their importance as potential targets to modulate GPCR signaling pathways.

Quantitative analysis of human ras localization and function in the fission yeast Schizosaccharomyces pombe.

Ras signalling is central to fundamental and diverse cellular processes. In higher eukaryotes ras signalling is highly complex, involving multiple isoforms, regulatory proteins and effectors. As a consequence, the study of ras activity in mammalian systems presents a number of technical challenges. The model organism Schizosaccharomyces pombe has previously proved a key system for the study of human signalling components and provides an ideal model for the study of ras, as it contains just one ras protein (Ras1p), which is non-essential and controls a number of downstream processes. Here we present data demonstrating the quantitative analysis of three distinct Ras1-related signalling outputs, utilizing the three most abundant human ras isoforms, H-Ras, N-Ras and K-Ras4B, in Sz. pombe. Further, we have characterized the localization of these three human ras isoforms in Sz. pombe, utilizing quantitative image analysis techniques. These data indicate that all three human ras isoforms are functional in fission yeast, displaying differing localization patterns which correlate strongly with function in the regulation of pheromone response and cell shape. These data demonstrate that such yeast strains could provide powerful tools for the investigation of ras biology, and potentially in the development of cancer therapies.

Special features of RAD Sequencing data: implications for genotyping.

Restriction site-associated DNA Sequencing (RAD-Seq) is an economical and efficient method for SNP discovery and genotyping. As with other sequencing-by-synthesis methods, RAD-Seq produces stochastic count data and requires sensitive analysis to develop or genotype markers accurately. We show that there are several sources of bias specific to RAD-Seq that are not explicitly addressed by current genotyping tools, namely restriction fragment bias, restriction site heterozygosity and PCR GC content bias. We explore the performance of existing analysis tools given these biases and discuss approaches to limiting or handling biases in RAD-Seq data. While these biases need to be taken seriously, we believe RAD loci affected by them can be excluded or processed with relative ease in most cases and that most RAD loci will be accurately genotyped by existing tools.

RAD-Seq derived markers flank the shell colour and banding loci of the Cepaea nemoralis supergene.

Studies on the classic shell colour and banding polymorphism of the land snail Cepaea played a crucial role in establishing the importance of natural selection in maintaining morphological variation. Cepaea is also a pre-eminent model for ecological genetics because the outward colour and banding phenotype is entirely genetically determined, primarily by a 'supergene' of at least five loci. Unfortunately, progress in understanding the evolution and maintenance of the Cepaea polymorphism stalled, partly because of a lack of genetic markers. With a view to re-establish Cepaea as a prominent model of molecular ecology, we made six laboratory crosses of Cepaea nemoralis, five of which segregated for shell ground colour (C) and the presence or absence of bands (B). First, scoring of colour and banding in 323 individuals found no recombination between the C and B loci of the supergene. Second, using restriction site-associated DNA sequencing (RAD-Seq) of two parents and 22 offspring, we identified 44 anonymous markers putatively linked to the colour (C) and banding (B) loci. The genotype of eleven of the most promising RAD-Seq markers was independently validated in the same 22 offspring, then up to a further 146 offspring were genotyped. The closest RAD-Seq markers scored are within ~0.6 centimorgan (cM) of the C-B supergene linkage group, with the combined loci together forming a 35.8 cM linkage map of markers that flank both sides of the Cepaea C-B supergene.

Genome-wide patterns of divergence and gene flow across a butterfly radiation.

The Heliconius butterflies are a diverse recent radiation comprising multiple levels of divergence with ongoing gene flow between species. The recently sequenced genome of Heliconius melpomene allowed us to investigate the genomic evolution of this group using dense RAD marker sequencing. Phylogenetic analysis of 54 individuals robustly supported reciprocal monophyly of H. melpomene and Heliconius cydno and refuted previous phylogenetic hypotheses that H. melpomene may be paraphylectic with respect to H. cydno. Heliconius timareta also formed a monophyletic clade closely related but distinct from H. cydno with Heliconius heurippa falling within this clade. We find evidence for genetic admixture between sympatric populations of the sister clades H. melpomene and H. cydno/timareta, particularly between H. cydno and H. melpomene from Central America and between H. timareta and H. melpomene from the eastern slopes of the Andes. Between races, divergence is primarily explained by isolation by distance and there is no detectable genetic population structure between parapatric races, suggesting that hybrid zones between races are not zones of secondary contact. Our results also support previous findings that colour pattern loci are shared between populations and species with similar colour pattern elements. Furthermore, this pattern is almost unique to these genomic regions, with only a very small number of other loci showing significant similarity between populations and species with similar colour patterns.

Characterisation of QTL-linked and genome-wide restriction site-associated DNA (RAD) markers in farmed Atlantic salmon.

BACKGROUND: Restriction site-associated DNA sequencing (RAD-Seq) is a genome complexity reduction technique that facilitates large-scale marker discovery and genotyping by sequencing. Recent applications of RAD-Seq have included linkage and QTL mapping with a particular focus on non-model species. In the current study, we have applied RAD-Seq to two Atlantic salmon families from a commercial breeding program. The offspring from these families were classified into resistant or susceptible based on survival/mortality in an Infectious Pancreatic Necrosis (IPN) challenge experiment, and putative homozygous resistant or susceptible genotype at a major IPN-resistance QTL. From each family, the genomic DNA of the two heterozygous parents and seven offspring of each IPN phenotype and genotype was digested with the SbfI enzyme and sequenced in multiplexed pools. RESULTS: Sequence was obtained from approximately 70,000 RAD loci in both families and a filtered set of 6,712 segregating SNPs were identified. Analyses of genome-wide RAD marker segregation patterns in the two families suggested SNP discovery on all 29 Atlantic salmon chromosome pairs, and highlighted the dearth of male recombination. The use of pedigreed samples allowed us to distinguish segregating SNPs from putative paralogous sequence variants resulting from the relatively recent genome duplication of salmonid species. Of the segregating SNPs, 50 were linked to the QTL. A subset of these QTL-linked SNPs were converted to a high-throughput assay and genotyped across large commercial populations of IPNV-challenged salmon fry. Several SNPs showed highly significant linkage and association with resistance to IPN, and population linkage-disequilibrium-based SNP tests for resistance were identified. CONCLUSIONS: We used RAD-Seq to successfully identify and characterise high-density genetic markers in pedigreed aquaculture Atlantic salmon. These results underline the effectiveness of RAD-Seq as a tool for rapid and efficient generation of QTL-targeted and genome-wide marker data in a large complex genome, and its possible utility in farmed animal selection programs.

The expression and activity of ß-catenin in the thalamus and its projections to the cerebral cortex in the mouse embryo.

BACKGROUND: The mammalian thalamus relays sensory information from the periphery to the cerebral cortex for cognitive processing via the thalamocortical tract. The thalamocortical tract forms during embryonic development controlled by mechanisms that are not fully understood. ß-catenin is a nuclear and cytosolic protein that transduces signals from secreted signaling molecules to regulate both cell motility via the cytoskeleton and gene expression in the nucleus. In this study we tested whether ß-catenin is likely to play a role in thalamocortical connectivity by examining its expression and activity in developing thalamic neurons and their axons. RESULTS: At embryonic day (E)15.5, the time when thalamocortical axonal projections are forming, we found that the thalamus is a site of particularly high ß-catenin mRNA and protein expression. As well as being expressed at high levels in thalamic cell bodies, ß-catenin protein is enriched in the axons and growth cones of thalamic axons and its growth cone concentration is sensitive to Netrin-1. Using mice carrying the ß-catenin reporter BAT-gal we find high levels of reporter activity in the thalamus. Further, Netrin-1 induces BAT-gal reporter expression and upregulates levels of endogenous transcripts encoding ß-actin and L1 proteins in cultured thalamic cells. We found that ß-catenin mRNA is enriched in thalamic axons and its 3'UTR is phylogenetically conserved and is able to direct heterologous mRNAs along the thalamic axon, where they can be translated. CONCLUSION: We provide evidence that ß-catenin protein is likely to be an important player in thalamocortcial development. It is abundant both in the nucleus and in the growth cones of post-mitotic thalamic cells during the development of thalamocortical connectivity and ß-catenin mRNA is targeted to thalamic axons and growth cones where it could potentially be translated. ß-catenin is involved in transducing the Netrin-1 signal to thalamic cells suggesting a mechanism by which Netrin-1 guides thalamocortical development.

Does a prefabricated gentamicin-impregnated, load-bearing spacer control periprosthetic hip infection?

INTRODUCTION: Treating deep infection following THA has been a challenge. While the standard treatment has remained a two-stage revision, spacer designs, incorporated antibiotics, and concentrations have varied. Since control of infection may relate to choice and concentration of antibiotics, it is important to report rates of control from various spacers. QUESTIONS/PURPOSES: We therefore determined (1) the rate of infection control and (2) complications associated with a prefabricated, load-bearing, gentamicin-impregnated hip spacer in treating periprosthetic infections of the hip. METHODS: We retrospectively reviewed 33 patients with periprosthetic THA infections treated with a prefabricated, partial load-bearing, gentamicin-impregnated hemiarthroplasty spacer. Thirty of the 33 patients underwent second stage reimplantation after a mean 15 weeks. We collected patient demographic data, laboratory values, infecting organism, size of spacer mold, antibiotic selection, complications, and infection control rates from two academic centers. Recurrent infection at last followup was determined by the presence of physical symptoms or signs or elevated serologic tests. The minimum followup was 24 months (mean, 43 months; range, 24-70 months). RESULTS: Twenty-eight of the 30 patients who underwent reimplantation remained infection-free at last followup: one patient became reinfected with a different organism secondary to wound problems; one became reinfected with the same organism, but was restaged with the mold used in this study, reimplanted, and subsequently remained free of infection. Two of the 33 patients had persistently elevated inflammatory markers at the completion of their first stage and were restaged with this mold; both underwent reimplantation and remained free of infection at latest followup. One of the 33 patients was satisfied and ambulatory with their spacer mold. There were no major complications. CONCLUSION: Our data supported the use of a partial load-bearing, gentamicin-impregnated hemiarthroplasty spacer in treating deep periprosthetic THA infections. LEVEL OF EVIDENCE: Level IV, therapeutic study. See Guidelines for Authors for a complete description of levels of evidence.

Chromosomal assembly of the nuclear genome of the endosymbiont-bearing trypanosomatid Angomonas deanei.

Angomonas deanei is an endosymbiont-bearing trypanosomatid with several highly fragmented genome assemblies and unknown chromosome number. We present an assembly of the A. deanei nuclear genome based on Oxford Nanopore sequence that resolves into 29 complete or close-to-complete chromosomes. The assembly has several previously unknown special features; it has a supernumerary chromosome, a chromosome with a 340-kb inversion, and there is a translocation between two chromosomes. We also present an updated annotation of the chromosomal genome with 10,365 protein-coding genes, 59 transfer RNAs, 26 ribosomal RNAs, and 62 noncoding RNAs.