Evolutionary studies are often limited by missing data that are critical to understanding the history of selection. Selection experiments, which reproduce rapid evolution under controlled conditions, are excellent tools to study how genomes evolve under selection. Here we present a genomic dissection of the Longshanks selection experiment, in which mice were selectively bred over 20 generations for longer tibiae relative to body mass, resulting in 13% longer tibiae in two replicates. We synthesized evolutionary theory, genome sequences and molecular genetics to understand the selection response and found that it involved both polygenic adaptation and discrete loci of major effect, with the strongest loci tending to be selected in parallel between replicates. We show that selection may favor de-repression of bone growth through inactivating two limb enhancers of an inhibitor, Nkx3-2. Our integrative genomic analyses thus show that it is possible to connect individual base-pair changes to the overall selection response.

Genetic variation segregates as linked sets of variants or haplotypes. Haplotypes and linkage are central to genetics and underpin virtually all genetic and selection analysis. Yet, genomic data often omit haplotype information due to constraints in sequencing technologies. Here, we present “haplotagging,” a simple, low-cost linked-read sequencing technique that allows sequencing of hundreds of individuals while retaining linkage information. We apply haplotagging to construct megabase-size haplotypes for over 600 individual butterflies (Heliconius erato and H. melpomene), which form overlapping hybrid zones across an elevational gradient in Ecuador. Haplotagging identifies loci controlling distinctive high- and lowland wing color patterns. Divergent haplotypes are found at the same major loci in both species, while chromosome rearrangements show no parallelism. Remarkably, in both species, the geographic clines for the major wing-pattern loci are displaced by 18 km, leading to the rise of a novel hybrid morph in the center of the hybrid zone. We propose that shared warning signaling (Müllerian mimicry) may couple the cline shifts seen in both species and facilitate the parallel coemergence of a novel hybrid morph in both comimetic species. Our results show the power of efficient haplotyping methods when combined with large-scale sequencing data from natural populations.

One-sentence summary:Haplotagging, a novel linked-read sequencing technique that enables whole genome haplotyping in large populations, reveals the formation of a novel hybrid race in parallel hybrid zones of two co-mimicking Heliconius butterfly species through strikingly parallel divergences in their genomes. Short title:Haplotagging reveals parallel formation of hybrid races Abstract Genetic variation segregates as linked sets of variants, or haplotypes. Haplotypes and linkage are central to genetics and underpin virtually all genetic and selection analysis. And yet, genomic data often lack haplotype information, due to constraints in sequencing technologies. Here we present "haplotagging", a simple, low-cost linked-read sequencing technique that allows sequencing of hundreds of individuals while retaining linkage information. We apply haplotagging to construct megabasesize haplotypes for over 600 individual butterflies (Heliconius erato and H. melpomene), which form overlapping hybrid zones across an elevational gradient in Ecuador. Haplotagging identifies loci controlling distinctive high-and lowland wing color patterns. Divergent haplotypes are found at the same major loci in both species, while chromosome rearrangements show no parallelism. Remarkably, in both species the geographic clines for the major wing pattern loci are displaced by 18 km, leading to the rise of a novel hybrid morph in the centre of the hybrid zone.We propose that shared warning signalling (Müllerian mimicry) may couple the cline shifts seen in both species, and facilitate the parallel co-emergence of a novel hybrid morph in both co-mimetic species. Our results show the power of efficient haplotyping methods when combined with large-scale sequencing data from natural populations.

Pannexins are a newly discovered three-member family of proteins expressed in the brain and peripheral tissues that belong to the superfamily of gap junction proteins. However, in mammals pannexins do not form gap junctions, and their expression and function in the pituitary gland have not been studied. Here we show that the rat pituitary gland expresses mRNA and protein transcripts of pannexins 1 and 2 but not pannexin 3. Pannexin 1 was more abundantly expressed in the anterior lobe, whereas pannexin 2 was more abundantly expressed in the intermediate and posterior pituitary. Pannexin 1 was identified in corticotrophs and a fraction of somatotrophs, the S100-positive pituicytes of the posterior pituitary and AtT-20 (mouse pituitary adrenocorticotropin-secreting cells) and rat immortalized pituitary cells secreting prolactin, whereas pannexin 2 was detected in the S100-positive folliculostellate cells of the anterior pituitary, melanotrophs of the intermediate lobe, and vasopressin-containing axons and nerve endings in the posterior lobe. Overexpression of pannexins 1 and 2 in AtT-20 pituitary cells enhanced the release of ATP in the extracellular medium, which was blocked by the gap junction inhibitor carbenoxolone. Basal ATP release in At-T20 cells was also suppressed by down-regulating the expression of endogenous pannexin 1 but not pannexin 2 with their short interfering RNAs. These results indicate that pannexins may provide a pathway for delivery of ATP, which is a native agonist for numerous P2X cationic channels and G protein-coupled P2Y receptors endogenously expressed in the pituitary gland.

Atrazine, one of the most commonly used herbicides worldwide, acts as an endocrine disruptor, but the mechanism of its action has not been characterized. In this study, we show that atrazine rapidly increases cAMP levels in cultured rat pituitary and testicular Leydig cells in a concentration-dependent manner, but less effectively than 3-isobutyl-1-methylxanthine, a competitive non-specific inhibitor of phosphodiesterases (PDEs). In forskolin (an activator of adenylyl cyclase)- and probenecid (an inhibitor of cyclic nucleotide transporters)-treated cells, but not in 3-isobutyl-1-methylxanthine-treated cells, atrazine further increased cAMP levels, indicating that inhibition of PDEs accounts for accumulation of cAMP. In contrast to cAMP, atrazine did not alter cGMP levels, further indicating that it inhibits cAMP-specific PDEs. Atrazine-induced changes in cAMP levels were sufficient to stimulate prolactin release in pituitary cells and androgen production in Leydig cells, indicating that it acts as an endocrine disrupter both in cells that secrete by exocytosis of prestored hormones and in cells that secrete by de novo hormone synthesis. Rolipram abolished the stimulatory effect of atrazine on cAMP release in both cell types, suggesting that it acts as an inhibitor of PDE4s, isoforms whose mRNA transcripts dominate in pituitary and Leydig cells together with mRNA for PDE8A. In contrast, immortalized lacto-somatotrophs showed low expression of these mRNA transcripts and several fold higher cAMP levels compared to normal pituitary cells, and atrazine was unable to further increase cAMP levels. These results indicate that atrazine acts as a general endocrine disrupter by inhibiting cAMP-specific PDE4s.

Understanding the production, response, and genetics of signals used in mate choice can inform our understanding of the evolution of both intraspecific mate choice and reproductive isolation. Sex pheromones are important for courtship and mate choice in many insects, but we know relatively little of their role in butterflies. The butterfly Heliconius melpomene uses a complex blend of wing androconial compounds during courtship. Electroantennography in H. melpomene and its close relative Heliconius cydno showed that responses to androconial extracts were not species specific. Females of both species responded equally strongly to extracts of both species, suggesting conservation of peripheral nervous system elements across the two species. Individual blend components provoked little to no response, with the exception of octadecanal, a major component of the H. melpomene blend. Supplementing octadecanal on the wings of octadecanal‐rich H. melpomene males led to an increase in the time until mating, demonstrating the bioactivity of octadecanal in Heliconius. Using quantitative trait locus (QTL) mapping, we identified a single locus on chromosome 20 responsible for 41% of the parental species’ difference in octadecanal production. This QTL does not overlap with any of the major wing color or mate choice loci, nor does it overlap with known regions of elevated or reduced FST. A set of 16 candidate fatty acid biosynthesis genes lies underneath the QTL. Pheromones in Heliconius carry information relevant for mate choice and are under simple genetic control, suggesting they could be important during speciation.

Transgenic mice expressing the tdimer2(12) form of Discosoma red fluorescent protein under control of the proopiomelanocortin gene's regulatory elements are a useful model for studying corticotrophs. Using these mice, we studied the ion channels and mechanisms controlling corticotroph excitability. Corticotrophs were either quiescent or electrically active, with a 22-mV difference in the resting membrane potential (RMP) between the 2 groups. In quiescent cells, CRH depolarized the membrane, leading to initial single spiking and sustained bursting; in active cells, CRH further facilitated or inhibited electrical activity and calcium spiking, depending on the initial activity pattern and CRH concentration. The stimulatory but not inhibitory action of CRH on electrical activity was mimicked by cAMP independently of the presence or absence of arachidonic acid. Removal of bath sodium silenced spiking and hyperpolarized the majority of cells; in contrast, the removal of bath calcium did not affect RMP but reduced CRH-induced depolarization, which abolished bursting electrical activity and decreased the spiking frequency but not the amplitude of single spikes. Corticotrophs with inhibited voltage-gated sodium channels fired calcium-dependent action potentials, whereas cells with inhibited L-type calcium channels fired sodium-dependent spikes; blockade of both channels abolished spiking without affecting the RMP. These results indicate that the background voltage-insensitive sodium conductance influences RMP, the CRH-depolarization current is driven by a cationic conductance, and the interplay between voltage-gated sodium and calcium channels plays a critical role in determining the status and pattern of electrical activity and calcium signaling.

Strongyloidiasis is a much-neglected but sometimes fatal soil born helminthiasis. The causing agent, the small intestinal parasitic nematode Strongyloides stercoralis can reproduce sexually through the indirect/heterogonic life cycle, or asexually through the auto-infective or the direct/homogonic life cycles. Usually, among the progeny of the parasitic females both, parthenogenetic parasitic (females only) and sexual free-living (females and males) individuals, are present simultaneously. We isolated S . stercoralis from people living in a village with a high incidence of parasitic helminths, in particular liver flukes ( Clonorchis sinensis ) and hookworms, in the southern Chinese province Guangxi. We determined nuclear and mitochondrial DNA sequences of individual S . stercoralis isolated from this village and from close by hospitals and we compared these S . stercoralis among themselves and with selected published S . stercoralis from other geographic locations. For comparison, we also analyzed the hookworms present in the same location. We found that, compared to earlier studies of S . stercoralis populations in South East Asia, all S . stercoralis sampled in our study area were very closely related, suggesting a recent common source of infection for all patients. In contrast, the hookworms from the same location, while all belonging to the species Necator americanus , showed rather extensive genetic diversity even within host individuals. Different from earlier studies conducted in other geographic locations, almost all S . stercoralis in this study appeared heterozygous for different sequence variants of the 18S rDNA hypervariable regions (HVR) I and IV. In contrast to earlier investigations, except for three males, all S . stercoralis we isolated in this study were infective larvae, suggesting that the sampled population reproduces predominantly, if not exclusively through the clonal life cycles. Consistently, whole genome sequencing of individual worms revealed higher heterozygosity than reported earlier for likely sexual populations of S . stercoralis . Elevated heterozygosity is frequently associated with asexual clonal reproduction.