Comparative genomics provides insights into molecular adaptation to hypermetamorphosis and cantharidin metabolism in blister beetles (Coleoptera: Meloidae)

Blister beetles (Coleoptera: Meloidae) are currently subdivided into three subfamilies: Eleticinae (a basal group), Nemognathinae, and Meloinae. These are all characterized by the endogenous production of the defensive terpene cantharidin (CA), whereas the two most derived subfamilies show a hypermetamorphic larval development. Here, we provide novel draft genome assemblies of five species sampled across the three blister beetle subfamilies (Iselma pallidipennis, Stenodera caucasica, Zonitis immaculata, Lydus trimaculatus, and Mylabris variabilis) and performed a comparative analysis with other available Meloidae genomes and the closely-related canthariphilous species (Pyrochroa serraticornis) to disclose adaptations at a molecular level. Our results highlighted the expansion and selection of genes potentially responsible for CA production and metabolism, as well as its mobilization and vesicular compartmentalization. Furthermore, we observed adaptive selection patterns and gain of genes devoted to epigenetic regulation, development, and morphogenesis, possibly related to hypermetamorphosis. We hypothesize that most genetic adaptations occurred to support both CA biosynthesis and hypermetamorphosis, two crucial aspects of Meloidae biology that likely contributed to their evolutionary success.

Genome-wide diversity, population structure and signatures of inbreeding in the African buffalo in Mozambique

The African buffalo, Syncerus caffer, is a key species in African ecosystems. Like other large herbivores, it plays a fundamental role in its habitat acting as an ecosystem engineer. Over the last few centuries, African buffalo populations have declined because of range contraction and demographic decline caused by direct or indirect human activities. In Mozambique, historically home to large buffalo herds, the combined effect of colonialism and subsequent civil wars has created a critical situation that urgently needs to be addressed. In this study, we focused on the analysis of genetic diversity of Syncerus caffer caffer populations from six areas of Mozambique. Using genome-wide SNPs obtained from ddRAD sequencing, we examined the population structure across the country, estimated gene flow between areas under conservation management, including national reserves, and assessed the inbreeding coefficients. Our results indicate that all studied populations of Syncerus caffer caffer are genetically depauperate, with a high level of inbreeding. Moreover, buffaloes in Mozambique present a significant population differentiation between southern and central areas. We found an unexpected genotype in the Gorongosa National Park, where buffaloes experienced a dramatic population size reduction, that shares a common ancestry with southern populations of Catuane and Namaacha. This could suggest the past occurrence of a connection between southern and central Mozambique and that the observed population structuring could reflect recent events of anthropogenic origin. All the populations analysed showed high levels of homozygosity, likely due to extensive inbreeding over the last few decades, which could have increased the frequency of recessive deleterious alleles. Improving the resilience of Syncerus caffer caffer in Mozambique is essential for preserving the ecosystem integrity. The most viable approach appears to be facilitating translocations and re-establishing connectivity between isolated herds. However, our results also highlight the importance of assessing intraspecific genetic diversity when considering interventions aimed at enhancing population viability such as selecting suitable source populations.

Limited introgression from non-native commercial strains and signatures of adaptation in the key pollinator Bombus terrestris

Insect pollination is fundamental for natural ecosystems and agricultural crops. The bumblebee species Bombus terrestris has become a popular choice for commercial crop pollination worldwide due to its effectiveness and ease of mass rearing. Bumblebee colonies are mass produced for the pollination of more than 20 crops and imported into over 50 countries including countries outside their native ranges, and the risk of invasion by commercial non-native bumblebees is considered an emerging issue for global conservation and biological diversity. Here, we use genome-wide data from seven wild populations close to and far from farms using commercial colonies, as well as commercial populations, to investigate the implications of utilizing commercial bumblebee subspecies in the UK. We find evidence for generally low levels of introgression between commercial and wild bees, with higher admixture proportions in the bees occurring close to farms. We identify genomic regions putatively involved in local and global adaptation, and genes in locally adaptive regions were found to be enriched for functions related to taste receptor activity, oxidoreductase activity, fatty acid and lipid biosynthetic processes. Despite more than 30 years of bumblebee colony importation into the UK, we observe low impact on the genetic integrity of local B. terrestris populations, but we highlight that even limited introgression might negatively affect locally adapted populations.

Comparative ontogenetic and transcriptomic analyses shed light on color pattern divergence in cichlid fishes

Stripe patterns are a striking example for a repeatedly evolved color pattern. In the African adaptive radiations of cichlid fishes, stripes evolved several times independently. Previously, it has been suggested that regulatory evolution of a single gene, agouti-related-peptide 2 (agrp2), explains the evolutionary lability of this trait. Here, using a comparative transcriptomic approach, we performed comparisons between (adult) striped and nonstriped cichlid fishes of representatives of Lake Victoria and the two major clades of Lake Malawi (mbuna and non-mbuna lineage). We identify agrp2 to be differentially expressed across all pairwise comparisons, reaffirming its association with stripe pattern divergence. We therefore also provide evidence that agrp2 is associated with the loss of the nonstereotypic oblique stripe of Mylochromis mola. Complementary ontogenetic data give insights into the development of stripe patterns as well as vertical bar patterns that both develop postembryonically. Lastly, using the Lake Victoria species pair Haplochromis sauvagei and Pundamilia nyererei, we investigated the differences between melanic and non-melanic regions to identify additional genes that contribute to the formation of stripes. Expression differences-that most importantly also do not include agrp2-are surprisingly small. This suggests, at least in this species pair, that the stripe phenotype might be caused by a combination of more subtle transcriptomic differences or cellular changes without transcriptional correlates. In summary, our comprehensive analysis highlights the ontogenetic and adult transcriptomic differences between cichlids with different color patterns and serves as a basis for further investigation of the mechanistic underpinnings of their diversification.

Comparative genomics of the waterfowl innate immune system

Animal species differ considerably in their ability to fight off infections. Finding the genetic basis of these differences is not easy, as the immune response is comprised of a complex network of proteins that interact with one another to defend the body against infection. Here, we used population- and comparative genomics to study the evolutionary forces acting on the innate immune system in natural hosts of the avian influenza virus (AIV). For this purpose, we used a combination of hybrid capture, next- generation sequencing and published genomes to examine genetic diversity, divergence, and signatures of selection in 127 innate immune genes at a micro- and macroevolutionary time scale in 26 species of waterfowl. We show across multiple immune pathways (AIV-, toll-like-, and RIG-I -like receptors signalling pathways) that genes involved genes in pathogen detection (i.e., toll-like receptors) and direct pathogen inhibition (i.e., antimicrobial peptides and interferon-stimulated genes), as well as host proteins targeted by viral antagonist proteins (i.e., mitochondrial antiviral-signaling protein, [MAVS]) are more likely to be polymorphic, genetically divergent, and under positive selection than other innate immune genes. Our results demonstrate that selective forces vary across innate immune signaling signalling pathways in waterfowl, and we present candidate genes that may contribute to differences in susceptibility and resistance to infectious diseases in wild birds, and that may be manipulated by viruses. Our findings improve our understanding of the interplay between host genetics and pathogens, and offer the opportunity for new insights into pathogenesis and potential drug targets.

Bacterial microcompartments for isethionate desulfonation in the taurine-degrading human-gut bacterium Bilophila wadsworthia

Background

Bilophila wadsworthia, a strictly anaerobic, sulfite-reducing bacterium and common member of the human gut microbiota, has been associated with diseases such as appendicitis and colitis. It is specialized on organosulfonate respiration for energy conservation, i.e., utilization of dietary and host-derived organosulfonates, such as taurine (2-aminoethansulfonate), as sulfite donors for sulfite respiration, producing hydrogen sulfide (H₂S), an important intestinal metabolite that may have beneficial as well as detrimental effects on the colonic environment. Its taurine desulfonation pathway involves the glycyl radical enzyme (GRE) isethionate sulfite-lyase (IslAB), which cleaves isethionate (2-hydroxyethanesulfonate) into acetaldehyde and sulfite.

Results

We demonstrate that taurine metabolism in B. wadsworthia 3.1.6 involves bacterial microcompartments (BMCs). First, we confirmed taurine-inducible production of BMCs by proteomic, transcriptomic and ultra-thin sectioning and electron-microscopical analyses. Then, we isolated BMCs from taurine-grown cells by density-gradient ultracentrifugation and analyzed their composition by proteomics as well as by enzyme assays, which suggested that the GRE IslAB and acetaldehyde dehydrogenase are located inside of the BMCs. Finally, we are discussing the recycling of cofactors in the IslAB-BMCs and a potential shuttling of electrons across the BMC shell by a potential iron-sulfur (FeS) cluster-containing shell protein identified by sequence analysis.

Conclusions

We characterized a novel subclass of BMCs and broadened the spectrum of reactions known to take place enclosed in BMCs, which is of biotechnological interest. We also provided more details on the energy metabolism of the opportunistic pathobiont B. wadsworthia and on microbial H₂S production in the human gut.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02386-w.

Phosphitispora fastidiosa gen. nov. sp. nov., a new dissimilatory phosphite-oxidizing anaerobic bacterium isolated from anaerobic sewage sludge

A new strictly anaerobic bacterium, strain DYL19^T, was enriched and isolated with phosphite as the sole electron donor and CO₂ as a single carbon source and electron acceptor from anaerobic sewage sludge sampled at a sewage treatment plant in Constance, Germany. It is a Gram-positive, spore-forming, slightly curved, rod-shaped bacterium which oxidizes phosphite to phosphate while reducing CO₂ to biomass and small amounts of acetate. Optimal growth is observed at 30 °C, pH 7.2, with a doubling time of 3 days. Beyond phosphite, no further inorganic or organic electron donor can be used, and no other electron acceptor than CO₂ is reduced. Sulphate inhibits growth with phosphite and CO₂. The G+C content is 45.95 mol%, and dimethylmenaquinone-7 is the only quinone detectable in the cells. On the basis of 16S rRNA gene sequence analysis and other chemotaxonomic properties, strain DYL19^T is described as the type strain of a new genus and species, Phosphitispora fastidiosa gen. nov., sp. nov.

Rapid adaptive radiation in a hillstream cyprinid fish in the East African White Nile River basin

Adaptive radiation of freshwater fishes was long thought to be possible only in lacustrine environments. Recently, several studies have shown that riverine and stream environments also provide the ecological opportunity for adaptive radiation. In this study, we report on a riverine adaptive radiation of six ecomorphs of cyprinid hillstream fishes of the genus Garra in a river located in the Ethiopian Highlands in East Africa. Garra are predominantly highly specialized algae-scrapers with a wide distribution ranging from Southeast Asia to West Africa. However, adaptive phenotypic diversification in mouth type, sucking disc morphology, gut length and body shape have probably been found among these ecomorphs in a single Ethiopian river. Moreover, we found two novel phenotypes of Garra ("thick-lipped" and "predatory") that had not been discovered before in this species-rich genus (>160 species). Mitochondrial and genome-wide data suggest monophyletic, intrabasin evolution of Garra phenotypic diversity with signatures of gene flow from other local populations. Although sympatric ecomorphs are genetically distinct and can be considered to being young species as suggested by genome-wide single nucleotide polymorphism data, mitochondrial DNA was unable to identify any genetic structure suggesting recent and rapid speciation events. Some data suggest a hybrid origin of the novel "thick-lipped" ecomorph. Here we highlight how, driven by ecological opportunity, an ancestral trophically highly specialized lineage is likely to have rapidly radiated in a riverine environment promoted by the evolution of novel feeding strategies.

Capturing the rapidly evolving study of adaptation

Research on the genomics of adaptation is rapidly changing. In the last few decades, progress in this area has been driven by methodological advances, not only in the way increasingly large amounts of molecular data are generated (e.g. with high-throughput sequencing), but also in the way these data are analysed. This includes a growing appreciation and quantitative treatment of covariation among units within the same data type (e.g. genes) or across data types (e.g. genes and phenotypes). The development and adoption of more and more integrative tools have resulted in richer and more interesting empirical work. This special issue - comprising methodological, empirical, and review papers - aims to capture a 'snapshot' of this rapidly evolving field. We discuss in particular three important themes in the study of adaptation: the genetic architecture of adaptive variation, protein-coding and regulatory changes, and parallel evolution. We highlight how more traditional key themes in the study of genetic architecture (e.g. the number of loci underlying adaptive traits and the distribution of their effects) are now being complemented by other factors (e.g. how patterns of linkage and number of loci interact to affect the ability to adapt). Similarly, apart from addressing the relative importance of protein-coding and regulatory changes, we now have the tools to look in-depth at specific types of regulatory variation to gain a clearer picture of regulatory networks. Finally, parallel evolution has always been central to the study of adaptation, but now we are often able to address the question of whether - and to what extent - parallelism at the organismal or phenotypic level is matched by parallelism at the genetic level. Perhaps most importantly, we can now determine what mechanisms are driving parallelism (or lack thereof) across levels of biological organization. All these recent methodological developments open up new directions for future studies of adaptive changes across traits, levels of biological organization, demographic contexts and time scales.

Pelorhabdus rhamnosifermentans gen. nov., sp. nov., a strictly anaerobic rhamnose degrader from freshwater lake sediment

A rhamnose-degrading bacterium, strain BoRhaA^T, was isolated from profundal sediment of Lake Constance in agar dilution series with l-rhamnose as substrate and with a background lawn of Methanospirillum hungatei. The isolated strain was a motile rod that stained Gram positive. Growth was observed within a pH range of 4.0-7.5 and a temperature range of 15-30°C. Fermentation products of rhamnose or glucose were acetate, propionate, ethanol, butyrate, and 1-propanol. The G+C content was 40.6% G+C. The dominant fatty acids are C_16:1ω9c, i-C_13:03OH, C_16:0 and C_17:1ω8c with 8-21% relative abundance. Polar lipids were glycolipids, phosphatidylethanolamine, phosphoaminolipid and other lipids, of which phosphatidylethanolamine was most abundant. The sequence of the 16S rRNA gene of the new isolate matches the sequence of its closest relative Anaerosporomusa subterranea to 92.4%. A comparison of the genome with this strain showed 60.2% genome-wide average amino acid identity (AAI), comparisons with other type strains showed a maximum of 62.7% AAI. Thus, the definition of a new genus is justified for which we propose the name Pelorhabdus. For strain BoRhaA^T, we propose the name Pelorhabdus rhamnosifermentans gen. nov., sp. nov., with strain BoRhaA^T (DSM 111565^T = JCM 39158^T) as the type strain.

Towards complete and error-free genome assemblies of all vertebrate species

High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species1-4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.

The comparative genomic landscape of adaptive radiation in crater lake cichlid fishes

Factors ranging from ecological opportunity to genome composition might explain why only some lineages form adaptive radiations. While being rare, particular systems can provide natural experiments within an identical ecological setting where species numbers and phenotypic divergence in two closely related lineages are notably different. We investigated one such natural experiment using two de novo assembled and 40 resequenced genomes and asked why two closely related Neotropical cichlid fish lineages, the Amphilophus citrinellus species complex (Midas cichlids; radiating) and Archocentrus centrarchus (Flyer cichlid; nonradiating), have resulted in such disparate evolutionary outcomes. Although both lineages inhabit many of the same Nicaraguan lakes, whole-genome inferred demography suggests that priority effects are not likely to be the cause of the dissimilarities. Also, genome-wide levels of selection, transposable element dynamics, gene family expansion, major chromosomal rearrangements and the number of genes under positive selection were not markedly different between the two lineages. To more finely investigate particular subsets of the genome that have undergone adaptive divergence in Midas cichlids, we also examined if there was evidence for 'molecular pre-adaptation' in regions identified by QTL mapping of repeatedly diverging adaptive traits. Although most of our analyses failed to pinpoint substantial genomic differences, we did identify functional categories containing many genes under positive selection that provide candidates for future studies on the propensity of Midas cichlids to radiate. Our results point to a disproportionate role of local, rather than genome-wide factors underlying the propensity for these cichlid fishes to adaptively radiate.

Giant lungfish genome elucidates the conquest of land by vertebrates

Lungfishes belong to lobe-fined fish (Sarcopterygii) that, in the Devonian period, 'conquered' the land and ultimately gave rise to all land vertebrates, including humans1-3. Here we determine the chromosome-quality genome of the Australian lungfish (Neoceratodus forsteri), which is known to have the largest genome of any animal. The vast size of this genome, which is about 14× larger than that of humans, is attributable mostly to huge intergenic regions and introns with high repeat content (around 90%), the components of which resemble those of tetrapods (comprising mainly long interspersed nuclear elements) more than they do those of ray-finned fish. The lungfish genome continues to expand independently (its transposable elements are still active), through mechanisms different to those of the enormous genomes of salamanders. The 17 fully assembled lungfish macrochromosomes maintain synteny to other vertebrate chromosomes, and all microchromosomes maintain conserved ancient homology with the ancestral vertebrate karyotype. Our phylogenomic analyses confirm previous reports that lungfish occupy a key evolutionary position as the closest living relatives to tetrapods4,5, underscoring the importance of lungfish for understanding innovations associated with terrestrialization. Lungfish preadaptations to living on land include the gain of limb-like expression in developmental genes such as hoxc13 and sall1 in their lobed fins. Increased rates of evolution and the duplication of genes associated with obligate air-breathing, such as lung surfactants and the expansion of odorant receptor gene families (which encode proteins involved in detecting airborne odours), contribute to the tetrapod-like biology of lungfishes. These findings advance our understanding of this major transition during vertebrate evolution.

Does divergence from normal patterns of integration increase as chromosomal fusions increase in number? A test on a house mouse hybrid zone

Chromosomal evolution is widely considered an important driver of speciation because it can promote the establishment of reproductive barriers. Karyotypic reorganization is also expected to affect the mean phenotype, as well as its development and patterns of phenotypic integration, through processes such as variation in genetic linkage between quantitative trait loci or between regulatory regions and their targets. Here we explore the relationship between chromosomal evolution and phenotypic integration by analyzing a well-known house mouse parapatric contact zone between a highly derived Robertsonian (Rb) race (2n = 22) and populations with standard karyotype (2n = 40). Populations with hybrid karyotypes are scattered throughout the hybrid zone connecting the two parental races. Using mandible shape data and geometric morphometrics, we test the hypothesis that patterns of integration progressively diverge from the “normal” integration pattern observed in the standard race as they accumulate Rb fusions. We find that the main pattern of integration observed between the posterior and anterior part of the mandible can be largely attributed to allometry. We find no support for a gradual increase in divergence from normal patterns of integration as fusions accumulate. Surprisingly, however, we find that the derived Rb race (2n = 22) has a distinct allometric trajectory compared with the standard race. Our results suggest that either individual fusions disproportionately affect patterns of integration or that there are mechanisms which “purge” extreme variants in hybrids (e.g. reduced fitness of hybrid shape).

Contrasting signatures of genomic divergence during sympatric speciation

The transition from 'well-marked varieties' of a single species into 'well-defined species'-especially in the absence of geographic barriers to gene flow (sympatric speciation)-has puzzled evolutionary biologists ever since Darwin1,2. Gene flow counteracts the buildup of genome-wide differentiation, which is a hallmark of speciation and increases the likelihood of the evolution of irreversible reproductive barriers (incompatibilities) that complete the speciation process3. Theory predicts that the genetic architecture of divergently selected traits can influence whether sympatric speciation occurs4, but empirical tests of this theory are scant because comprehensive data are difficult to collect and synthesize across species, owing to their unique biologies and evolutionary histories5. Here, within a young species complex of neotropical cichlid fishes (Amphilophus spp.), we analysed genomic divergence among populations and species. By generating a new genome assembly and re-sequencing 453 genomes, we uncovered the genetic architecture of traits that have been suggested to be important for divergence. Species that differ in monogenic or oligogenic traits that affect ecological performance and/or mate choice show remarkably localized genomic differentiation. By contrast, differentiation among species that have diverged in polygenic traits is genomically widespread and much higher overall, consistent with the evolution of effective and stable genome-wide barriers to gene flow. Thus, we conclude that simple trait architectures are not always as conducive to speciation with gene flow as previously suggested, whereas polygenic architectures can promote rapid and stable speciation in sympatry.

The mole genome reveals regulatory rearrangements associated with adaptive intersexuality

Linking genomic variation to phenotypical traits remains a major challenge in evolutionary genetics. In this study, we use phylogenomic strategies to investigate a distinctive trait among mammals: the development of masculinizing ovotestes in female moles. By combining a chromosome-scale genome assembly of the Iberian mole, Talpa occidentalis, with transcriptomic, epigenetic, and chromatin interaction datasets, we identify rearrangements altering the regulatory landscape of genes with distinct gonadal expression patterns. These include a tandem triplication involving CYP17A1, a gene controlling androgen synthesis, and an intrachromosomal inversion involving the pro-testicular growth factor gene FGF9, which is heterochronically expressed in mole ovotestes. Transgenic mice with a knock-in mole CYP17A1 enhancer or overexpressing FGF9 showed phenotypes recapitulating mole sexual features. Our results highlight how integrative genomic approaches can reveal the phenotypic impact of noncoding sequence changes.

Environmental and Intestinal Phylum Firmicutes Bacteria Metabolize the Plant Sugar Sulfoquinovose via a 6-Deoxy-6-sulfofructose Transaldolase Pathway

Bacterial degradation of the sugar sulfoquinovose (SQ, 6-deoxy-6-sulfoglucose) produced by plants, algae, and cyanobacteria, is an important component of the biogeochemical carbon and sulfur cycles. Here, we reveal a third biochemical pathway for primary SQ degradation in an aerobic Bacillus aryabhattai strain. An isomerase converts SQ to 6-deoxy-6-sulfofructose (SF). A novel transaldolase enzyme cleaves the SF to 3-sulfolactaldehyde (SLA), while the non-sulfonated C3-(glycerone)-moiety is transferred to an acceptor molecule, glyceraldehyde phosphate (GAP), yielding fructose-6-phosphate (F6P). Intestinal anaerobic bacteria such as Enterococcus gilvus, Clostridium symbiosum, and Eubacterium rectale strains also express transaldolase pathway gene clusters during fermentative growth with SQ. The now three known biochemical strategies for SQ catabolism reflect adaptations to the aerobic or anaerobic lifestyle of the different bacteria. The occurrence of these pathways in intestinal (family) Enterobacteriaceae and (phylum) Firmicutes strains further highlights a potential importance of metabolism of green-diet SQ by gut microbial communities to, ultimately, hydrogen sulfide.

Reconstructing the Evolutionary History of Chromosomal Races on Islands: A Genome-Wide Analysis of Natural House Mouse Populations

Chromosomal evolution is widely considered to be an important driver of speciation, as karyotypic reorganization can bring about the establishment of reproductive barriers between incipient species. One textbook example for genetic mechanisms of speciation are large-scale chromosomal rearrangements such as Robertsonian (Rb) fusions, a common class of structural variants that can drastically change the recombination landscape by suppressing crossing-over and influence gene expression by altering regulatory networks. Here, we explore the population structure and demographic patterns of a well-known house mouse Rb system in the Aeolian archipelago in Southern Italy using genome-wide data. By analyzing chromosomal regions characterized by different levels of recombination, we trace the evolutionary history of a set of Rb chromosomes occurring in different geographical locations and test whether chromosomal fusions have a single shared origin or occurred multiple times. Using a combination of phylogenetic and population genetic approaches, we find support for multiple, independent origins of three focal Rb chromosomes. The elucidation of the demographic patterns of the mouse populations within the Aeolian archipelago shows that an interplay between fixation of newly formed Rb chromosomes and hybridization events has contributed to shaping their current karyotypic distribution. Overall, our results illustrate that chromosome structure is much more dynamic than anticipated and emphasize the importance of large-scale chromosomal translocations in speciation.

Divergent Allometric Trajectories in Gene Expression and Coexpression Produce Species Differences in Sympatrically Speciating Midas Cichlid Fish

The mechanisms of speciation without geographic isolation (i.e., sympatric speciation) remain debated. This is due in part to the fact that the genomic landscape that could promote or hinder species divergence in the presence of gene flow is still largely unknown. However, intensive research is now centered on understanding the genetic architecture of adaptive traits associated with this process as well as how gene expression might affect these traits. Here, using RNA-Seq data, we investigated gene expression of sympatrically speciating benthic and limnetic Neotropical cichlid fishes at two developmental stages. First, we identified groups of coexpressed genes (modules) at each stage. Although there are a few large and well-preserved modules, most of the other modules are not preserved across life stages. Second, we show that later in development more and larger coexpression modules are associated with divergence between benthic and limnetic fish compared with the earlier life stage. This divergence between benthic and limnetic fish in coexpression mirrors divergence in overall expression between benthic and limnetic fish, which is more pronounced later in life. Our results reveal that already at 1-day posthatch benthic and limnetic fish diverge in (co)expression, and that this divergence becomes more substantial when fish are free-swimming but still unlikely to have divergent swimming and feeding habits. More importantly, our study describes how the coexpression of several genes through development, as opposed to individual genes, is associated with benthic–limnetic species differences, and how two morphogenetic trajectories diverge as fish grow older.

Conservation and novelty in the microRNA genomic landscape of hyperdiverse cichlid fishes

MicroRNAs (miRNAs) play crucial roles in the post-transcriptional control of messenger RNA (mRNA). These miRNA-mRNA regulatory networks are present in nearly all organisms and contribute to development, phenotypic divergence, and speciation. To examine the miRNA landscape of cichlid fishes, one of the most species-rich families of vertebrates, we profiled the expression of both miRNA and mRNA in a diverse set of cichlid lineages. Among these, we found that conserved miRNAs differ from recently arisen miRNAs (i.e. lineage specific) in average expression levels, number of target sites, sequence variability, and physical clustering patterns in the genome. Furthermore, conserved miRNA target sites tend to be enriched at the 5' end of protein-coding gene 3' UTRs. Consistent with the presumed regulatory role of miRNAs, we detected more negative correlations between the expression of miRNA-mRNA functional pairs than in random pairings. Finally, we provide evidence that novel miRNA targets sites are enriched in genes involved in protein synthesis pathways. Our results show how conserved and evolutionarily novel miRNAs differ in their contribution to the genomic landscape and highlight their particular evolutionary roles in the adaptive diversification of cichlids.

MicroRNA Gene Regulation in Extremely Young and Parallel Adaptive Radiations of Crater Lake Cichlid Fish

Cichlid fishes provide textbook examples of explosive phenotypic diversification and sympatric speciation, thereby making them ideal systems for studying the molecular mechanisms underlying rapid lineage divergence. Despite the fact that gene regulation provides a critical link between diversification in gene function and speciation, many genomic regulatory mechanisms such as microRNAs (miRNAs) have received little attention in these rapidly diversifying groups. Therefore, we investigated the posttranscriptional regulatory role of miRNAs in the repeated sympatric divergence of Midas cichlids (Amphilophus spp.) from Nicaraguan crater lakes. Using miRNA and mRNA sequencing of embryos from five Midas species, we first identified miRNA binding sites in mRNAs and highlighted the presences of a surprising number of novel miRNAs in these adaptively radiating species. Then, through analyses of expression levels, we identified putative miRNA/gene target pairs with negatively correlated expression level that were consistent with the role of miRNA in downregulating mRNA. Furthermore, we determined that several miRNA/gene pairs show convergent expression patterns associated with the repeated benthic/limnetic sympatric species divergence implicating these miRNAs as potential molecular mechanisms underlying replicated sympatric divergence. Finally, as these candidate miRNA/gene pairs may play a central role in phenotypic diversification in these cichlids, we characterized the expression domains of selected miRNAs and their target genes via in situ hybridization, providing further evidence that miRNA regulation likely plays a role in the Midas cichlid adaptive radiation. These results provide support for the hypothesis that extremely quickly evolving miRNA regulation can contribute to rapid evolutionary divergence even in the presence of gene flow.

Long-term experimental hybridisation results in the evolution of a new sex chromosome in swordtail fish

The remarkable diversity of sex determination mechanisms known in fish may be fuelled by exceptionally high rates of sex chromosome turnovers or transitions. However, the evolutionary causes and genomic mechanisms underlying this variation and instability are yet to be understood. Here we report on an over 30-year evolutionary experiment in which we tested the genomic consequences of hybridisation and selection between two Xiphophorus fish species with different sex chromosome systems. We find that introgression and imposing selection for pigmentation phenotypes results in the retention of an unexpectedly large maternally derived genomic region. During the hybridisation process, the sex-determining region of the X chromosome from one parental species was translocated to an autosome in the hybrids leading to the evolution of a new sex chromosome. Our results highlight the complexity of factors contributing to patterns observed in hybrid genomes, and we experimentally demonstrate that hybridisation can catalyze rapid evolution of a new sex chromosome.

Fish have a high diversity of sex-determining systems, but the mechanisms responsible for this are not well understood. Here, Franchini et al. show how hybridization and backcrossing have led to the evolution of a new sex chromosome in swordtail fish during 30 years of experimental evolution.

Anaerobic Degradation of the Plant Sugar Sulfoquinovose Concomitant With H2S Production: Escherichia coli K-12 and Desulfovibrio sp. Strain DF1 as Co-culture Model

Sulfoquinovose (SQ, 6-deoxy-6-sulfoglucose) is produced by plants and other phototrophs and its biodegradation is a relevant component of the biogeochemical carbon and sulfur cycles. SQ is known to be degraded by aerobic bacterial consortia in two tiers via C₃-organosulfonates as transient intermediates to CO₂, water and sulfate. In this study, we present a first laboratory model for anaerobic degradation of SQ by bacterial consortia in two tiers to acetate and hydrogen sulfide (H₂S). For the first tier, SQ-degrading Escherichia coli K-12 was used. It catalyzes the fermentation of SQ to 2,3-dihydroxypropane-1-sulfonate (DHPS), succinate, acetate and formate, thus, a novel type of mixed-acid fermentation. It employs the characterized SQ Embden-Meyerhof-Parnas pathway, as confirmed by mutational and proteomic analyses. For the second tier, a DHPS-degrading Desulfovibrio sp. isolate from anaerobic sewage sludge was used, strain DF1. It catalyzes another novel fermentation, of the DHPS to acetate and H₂S. Its DHPS desulfonation pathway was identified by differential proteomics and demonstrated by heterologously produced enzymes: DHPS is oxidized via 3-sulfolactaldehyde to 3-sulfolactate (SL) by two NAD⁺-dependent dehydrogenases (DhpA, SlaB); the SL is cleaved by an SL sulfite-lyase known from aerobic bacteria (SuyAB) to pyruvate and sulfite. The pyruvate is oxidized to acetate, while the sulfite is used as electron acceptor in respiration and reduced to H₂S. In conclusion, anaerobic sulfidogenic SQ degradation was demonstrated as a novel link in the biogeochemical sulfur cycle. SQ is also a constituent of the green-vegetable diet of herbivores and omnivores and H₂S production in the intestinal microbiome has many recognized and potential contributions to human health and disease. Hence, it is important to examine bacterial SQ degradation also in the human intestinal microbiome, in relation to H₂S production, dietary conditions and human health.

Phylogenomics uncovers early hybridization and adaptive loci shaping the radiation of Lake Tanganyika cichlid fishes

Lake Tanganyika is the oldest and phenotypically most diverse of the three East African cichlid fish adaptive radiations. It is also the cradle for the younger parallel haplochromine cichlid radiations in Lakes Malawi and Victoria. Despite its evolutionary significance, the relationships among the main Lake Tanganyika lineages remained unresolved, as did the general timescale of cichlid evolution. Here, we disentangle the deep phylogenetic structure of the Lake Tanganyika radiation using anchored phylogenomics and uncover hybridization at its base, as well as early in the haplochromine radiation. This suggests that hybridization might have facilitated these speciation bursts. Time-calibrated trees support that the radiation of Tanganyika cichlids coincided with lake formation and that Gondwanan vicariance concurred with the earliest splits in the cichlid family tree. Genes linked to key innovations show signals of introgression or positive selection following colonization of lake habitats and species' dietary adaptations are revealed as major drivers of colour vision evolution. These findings shed light onto the processes shaping the evolution of adaptive radiations.

Evolutionary divergence of 3' UTRs in cichlid fishes

Background

Post-transcriptional regulation is crucial for the control of eukaryotic gene expression and might contribute to adaptive divergence. The three prime untranslated regions (3’ UTRs), that are located downstream of protein-coding sequences, play important roles in post-transcriptional regulation. These regions contain functional elements that influence the fate of mRNAs and could be exceptionally important in groups such as rapidly evolving cichlid fishes.

Results

To examine cichlid 3’ UTR evolution, we 1) identified gene features in nine teleost genomes and 2) performed comparative analyses to assess evolutionary variation in length, functional motifs, and evolutionary rates of 3’ UTRs. In all nine teleost genomes, we found a smaller proportion of repetitive elements in 3’ UTRs than in the whole genome. We found that the 3’ UTRs in cichlids tend to be longer than those in non-cichlids, and this was associated, on average, with one more miRNA target per gene in cichlids. Moreover, we provided evidence that 3’ UTRs on average have evolved faster in cichlids than in non-cichlids. Finally, analyses of gene function suggested that both the top 5% longest and 5% most rapidly evolving 3’ UTRs in cichlids tended to be involved in ribosome-associated pathways and translation.

Conclusions

Our results reveal novel patterns of evolution in the 3’ UTRs of teleosts in general and cichlids in particular. The data suggest that 3’ UTRs might serve as important meta-regulators, regulators of other mechanisms governing post-transcriptional regulation, especially in groups like cichlids that have undergone extremely fast rates of phenotypic diversification and speciation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4821-8) contains supplementary material, which is available to authorized users.

Lateralized Feeding Behavior is Associated with Asymmetrical Neuroanatomy and Lateralized Gene Expressions in the Brain in Scale-Eating Cichlid Fish

Lateralized behavior ("handedness") is unusual, but consistently found across diverse animal lineages, including humans. It is thought to reflect brain anatomical and/or functional asymmetries, but its neuro-molecular mechanisms remain largely unknown. Lake Tanganyika scale-eating cichlid fish, Perissodus microlepis show pronounced asymmetry in their jaw morphology as well as handedness in feeding behavior-biting scales preferentially only from one or the other side of their victims. This makes them an ideal model in which to investigate potential laterality in neuroanatomy and transcription in the brain in relation to behavioral handedness. After determining behavioral handedness in P. microlepis (preferred attack side), we estimated the volume of the hemispheres of brain regions and captured their gene expression profiles. Our analyses revealed that the degree of behavioral handedness is mirrored at the level of neuroanatomical asymmetry, particularly in the tectum opticum. Transcriptome analyses showed that different brain regions (tectum opticum, telencephalon, hypothalamus, and cerebellum) display distinct expression patterns, potentially reflecting their developmental interrelationships. For numerous genes in each brain region, their extent of expression differences between hemispheres was found to be correlated with the degree of behavioral lateralization. Interestingly, the tectum opticum and telencephalon showed divergent biases on the direction of up- or down-regulation of the laterality candidate genes (e.g., grm2) in the hemispheres, highlighting the connection of handedness with gene expression profiles and the different roles of these brain regions. Hence, handedness in predation behavior may be caused by asymmetric size of brain hemispheres and also by lateralized gene expressions in the brain.

Animal tracking meets migration genomics: transcriptomic analysis of a partially migratory bird species

Seasonal migration is a widespread phenomenon, which is found in many different lineages of animals. This spectacular behaviour allows animals to avoid seasonally adverse environmental conditions to exploit more favourable habitats. Migration has been intensively studied in birds, which display astonishing variation in migration strategies, thus providing a powerful system for studying the ecological and evolutionary processes that shape migratory behaviour. Despite intensive research, the genetic basis of migration remains largely unknown. Here, we used state-of-the-art radio-tracking technology to characterize the migratory behaviour of a partially migratory population of European blackbirds (Turdus merula) in southern Germany. We compared gene expression of resident and migrant individuals using high-throughput transcriptomics in blood samples. Analyses of sequence variation revealed a nonsignificant genetic structure between blackbirds differing by their migratory phenotype. We detected only four differentially expressed genes between migrants and residents, which might be associated with hyperphagia, moulting and enhanced DNA replication and transcription. The most pronounced changes in gene expression occurred between migratory birds depending on when, in relation to their date of departure, blood was collected. Overall, the differentially expressed genes detected in this analysis may play crucial roles in determining the decision to migrate, or in controlling the physiological processes required for the onset of migration. These results provide new insights into, and testable hypotheses for, the molecular mechanisms controlling the migratory phenotype and its underlying physiological mechanisms in blackbirds and other migratory bird species.

quaddRAD: a new high-multiplexing and PCR duplicate removal ddRAD protocol produces novel evolutionary insights in a nonradiating cichlid lineage

The identification of thousands of variants across the genomes and their accurate genotyping are crucial for estimating the genetic parameters needed to address a host of molecular ecological and evolutionary questions. With rapid advances of massively parallel high-throughput sequencing technologies, several methods have recently been developed to access genomewide data on population variation. One of the most successful and widely used techniques relies on the combination of restriction enzymes and sequencing-by-synthesis: restriction-site-associated DNA sequencing (RADSeq). We developed a new, more time- and cost-efficient double-digest RAD paired-end protocol (quaddRAD) that simplifies and speeds up the identification of PCR duplicates and permits large-scale multiplexing. Assessing its performance on a technical data set, we also applied the quaddRAD method on population samples of a Neotropical cichlid fish lineage (Archocentrus centrarchus) to assess its genetic structure and demographic history. While we identified allopatric interlake genetic divergence, most likely driven by drift, no signature of sympatric divergence was detected. This differs from what has been observed in the clade of Midas cichlids (Amphilophus citrinellus spp.), another cichlid lineage that inhabits the same lakes and shares a similar demographic history, but has evolved into small-scale adaptive radiations via sympatric speciation. We demonstrate that quaddRAD is a robust and efficient method for genotyping a massive number and widely overlapping set of loci with high accuracy. Furthermore, the results on A. centrarchus open new research avenues providing an ideal system to investigate genome-level mechanisms that could alter the speciation potential of different but closely related cichlid lineages.

Genetic linkage of distinct adaptive traits in sympatrically speciating crater lake cichlid fish

Our understanding of how biological diversity arises is limited, especially in the case of speciation in the face of gene flow. Here we investigate the genomic basis of adaptive traits, focusing on a sympatrically diverging species pair of crater lake cichlid fishes. We identify the main quantitative trait loci (QTL) for two eco-morphological traits: body shape and pharyngeal jaw morphology. These traits diverge in parallel between benthic and limnetic species in the repeated adaptive radiations of this and other fish lineages. Remarkably, a single chromosomal region contains the highest effect size QTL for both traits. Transcriptomic data show that the QTL regions contain genes putatively under selection. Independent population genomic data corroborate QTL regions as areas of high differentiation between the sympatric sister species. Our results provide empirical support for current theoretical models that emphasize the importance of genetic linkage and pleiotropy in facilitating rapid divergence in sympatry.

The Role of microRNAs in the Repeated Parallel Diversification of Lineages of Midas Cichlid Fish from Nicaragua

Cichlid fishes are an ideal model system for studying biological diversification because they provide textbook examples of rapid speciation. To date, there has been little focus on the role of gene regulation during cichlid speciation. However, in recent years, gene regulation has been recognized as a powerful force linking diversification in gene function to speciation. Here, we investigated the potential role of miRNA regulation in the diversification of six cichlid species of the Midas cichlid lineage (Amphilophus spp.) inhabiting the Nicaraguan crater lakes. Using several genomic resources, we inferred 236 Midas miRNA genes that were used to predict the miRNA target sites on 8,232 Midas 3′-UTRs. Using population genomic calculations of SNP diversity, we found the miRNA genes to be more conserved than protein coding genes. In contrast to what has been observed in other cichlid fish, but similar to what has been typically found in other groups, we observed genomic signatures of purifying selection on the miRNA targets by comparing these sites with the less conserved nontarget portion of the 3′-UTRs. However, in one species pair that has putatively speciated sympatrically in crater Lake Apoyo, we recovered a different pattern of relaxed purifying selection and high genetic divergence at miRNA targets. Our results suggest that sequence evolution at miRNA binding sites could be a critical genomic mechanism contributing to the rapid phenotypic evolution of Midas cichlids.

Are sympatrically speciating Midas cichlid fish special? Patterns of morphological and genetic variation in the closely related species Archocentrus centrarchus

Established empirical cases of sympatric speciation are scarce, although there is an increasing consensus that sympatric speciation might be more common than previously thought. Midas cichlid fish are one of the few substantiated cases of sympatric speciation, and they formed repeated radiations in crater lakes. In contrast, in the same environment, such radiation patterns have not been observed in other species of cichlids and other families of fish. We analyze morphological and genetic variation in a cichlid species (Archocentrus centrarchus) that co-inhabits several crater lakes with the Midas species complex. In particular, we analyze variation in body and pharyngeal jaw shape (two ecologically important traits in sympatrically divergent Midas cichlids) and relate that to genetic variation in mitochondrial control region and microsatellites. Using these four datasets, we analyze variation between and within two Nicaraguan lakes: a crater lake where multiple Midas cichlids have been described and a lake where the source population lives. We do not observe any within-lake clustering consistent across morphological traits and genetic markers, suggesting the absence of sympatric divergence in A. centrarchus. Genetic differentiation between lakes was low and morphological divergence absent. Such morphological similarity between lakes is found not only in average morphology, but also when analyzing covariation between traits and degree of morphospace occupation. A combined analysis of the mitochondrial control region in A. centrarchus and Midas cichlids suggests that a difference between lineages in the timing of crater lake colonization cannot be invoked as an explanation for the difference in their levels of diversification. In light of our results, A. centrarchus represents the ideal candidate to study the genomic differences between these two lineages that might explain why some lineages are more likely to speciate and diverge in sympatry than others.

Chromosomal rearrangements, phenotypic variation and modularity: a case study from a contact zone between house mouse Robertsonian races in Central Italy

The Western European house mouse, Mus musculus domesticus, is well‐known for the high frequency of Robertsonian fusions that have rapidly produced more than 50 karyotipic races, making it an ideal model for studying the mechanisms of chromosomal speciation. The mouse mandible is one of the traits studied most intensively to investigate the effect of Robertsonian fusions on phenotypic variation within and between populations. This complex bone structure has also been widely used to study the level of integration between different morphogenetic units. Here, with the aim of testing the effect of different karyotypic assets on the morphology of the mouse mandible and on its level of modularity, we performed morphometric analyses of mice from a contact area between two highly metacentric races in Central Italy. We found no difference in size, while the mandible shape was found to be different between the two Robertsonian races, even after accounting for the genetic relationships among individuals and geographic proximity. Our results support the existence of two modules that indicate a certain degree of evolutionary independence, but no difference in the strength of modularity between chromosomal races. Moreover, the ascending ramus showed more pronounced interpopulation/race phenotypic differences than the alveolar region, an effect that could be associated to their different polygenic architecture. This study suggests that chromosomal rearrangements play a role in the house mouse phenotypic divergence, and that the two modules of the mouse mandible are differentially affected by environmental factors and genetic makeup.

Permanent draft genome sequence of sulfoquinovose-degrading Pseudomonas putida strain SQ1

Pseudomonas putida SQ1 was isolated for its ability to utilize the plant sugar sulfoquinovose (6-deoxy-6-sulfoglucose) for growth, in order to define its SQ-degradation pathway and the enzymes and genes involved. Here we describe the features of the organism, together with its draft genome sequence and annotation. The draft genome comprises 5,328,888 bp and is predicted to encode 5,824 protein-coding genes; the overall G + C content is 61.58 %. The genome annotation is being used for identification of proteins that might be involved in SQ degradation by peptide fingerprinting-mass spectrometry.

The venomous cocktail of the vampire snail Colubraria reticulata (Mollusca, Gastropoda)

Background

Hematophagy arose independently multiple times during metazoan evolution, with several lineages of vampire animals particularly diversified in invertebrates. However, the biochemistry of hematophagy has been studied in a few species of direct medical interest and is still underdeveloped in most invertebrates, as in general is the study of venom toxins. In cone snails, leeches, arthropods and snakes, the strong target specificity of venom toxins uniquely aligns them to industrial and academic pursuits (pharmacological applications, pest control etc.) and provides a biochemical tool for studying biological activities including cell signalling and immunological response. Neogastropod snails (cones, oyster drills etc.) are carnivorous and include active predators, scavengers, grazers on sessile invertebrates and hematophagous parasites; most of them use venoms to efficiently feed. It has been hypothesized that trophic innovations were the main drivers of rapid radiation of Neogastropoda in the late Cretaceous.

We present here the first molecular characterization of the alimentary secretion of a non-conoidean neogastropod, Colubraria reticulata. Colubrariids successfully feed on the blood of fishes, throughout the secretion into the host of a complex mixture of anaesthetics and anticoagulants. We used a NGS RNA-Seq approach, integrated with differential expression analyses and custom searches for putative secreted feeding-related proteins, to describe in detail the salivary and mid-oesophageal transcriptomes of this Mediterranean vampire snail, with functional and evolutionary insights on major families of bioactive molecules.

Results

A remarkably low level of overlap was observed between the gene expression in the two target tissues, which also contained a high percentage of putatively secreted proteins when compared to the whole body. At least 12 families of feeding-related proteins were identified, including: 1) anaesthetics, such as ShK Toxin-containing proteins and turripeptides (ion-channel blockers), Cysteine-rich secretory proteins (CRISPs), Adenosine Deaminase (ADA); 2) inhibitors of primary haemostasis, such as novel vWFA domain-containing proteins, the Ectonucleotide pyrophosphatase/phosphodiesterase family member 5 (ENPP5) and the wasp Antigen-5; 3) anticoagulants, such as TFPI-like multiple Kunitz-type protease inhibitors, Peptidases S1 (PS1), CAP/ShKT domain-containing proteins, Astacin metalloproteases and Astacin/ShKT domain-containing proteins; 4) additional proteins, such the Angiotensin-Converting Enzyme (ACE: vasopressive) and the cytolytic Porins.

Conclusions

Colubraria feeding physiology seems to involve inhibitors of both primary and secondary haemostasis, anaesthetics, a vasoconstrictive enzyme to reduce feeding time and tissue-degrading proteins such as Porins and Astacins. The complexity of Colubraria venomous cocktail and the divergence from the arsenal of the few neogastropods studied to date (mostly conoideans) suggest that biochemical diversification of neogastropods might be largely underestimated and worth of extensive investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1648-4) contains supplementary material, which is available to authorized users.

Transcriptomics of two evolutionary novelties: how to make a sperm-transfer organ out of an anal fin and a sexually selected "sword" out of a caudal fin

Swords are exaggerated male ornaments of swordtail fishes that have been of great interest to evolutionary biologists ever since Darwin described them in the Descent of Man (1871). They are a novel sexually selected trait derived from modified ventral caudal fin rays and are only found in the genus Xiphophorus. Another phylogenetically more widespread and older male trait is the gonopodium, an intromittent organ found in all poeciliid fishes, that is derived from a modified anal fin. Despite many evolutionary and behavioral studies on both traits, little is known so far about the molecular mechanisms underlying their development. By investigating transcriptomic changes (utilizing a RNA-Seq approach) in response to testosterone treatment in the swordtail fish, Xiphophorus hellerii, we aimed to better understand the architecture of the gene regulatory networks underpinning the development of these two evolutionary novelties. Large numbers of genes with tissue-specific expression patterns were identified. Among the "sword genes" those involved in embryonic organ development, sexual character development and coloration were highly expressed, while in the gonopodium rather more morphogenesis-related genes were found. Interestingly, many genes and genetic pathways are shared between both developing novel traits derived from median fins: the sword and the gonopodium. Our analyses show that a larger set of gene networks was co-opted during the development and evolution of the "older" gonopodium than in the "younger," and morphologically less complex trait, the sword. We provide a catalog of candidate genes for future efforts to dissect the development of those sexually selected exaggerated male traits in swordtails.

Genetic mapping of horizontal stripes in Lake Victoria cichlid fishes: benefits and pitfalls of using RAD markers for dense linkage mapping

The genetic dissection of naturally occurring phenotypes sheds light on many fundamental and longstanding questions in speciation and adaptation and is a central research topic in evolutionary biology. Until recently, forward-genetic approaches were virtually impossible to apply to nonmodel organisms, but the development of next-generation sequencing techniques eases this difficulty. Here, we use the ddRAD-seq method to map a colour trait with a known adaptive function in cichlid fishes, well-known textbook examples for rapid rates of speciation and astonishing phenotypic diversification. A suite of phenotypic key innovations is related to speciation and adaptation in cichlids, among which body coloration features prominently. The focal trait of this study, horizontal stripes, evolved in parallel in several cichlid radiations and is associated with piscivorous foraging behaviour. We conducted interspecific crosses between Haplochromis sauvagei and H. nyererei and constructed a linkage map with 867 SNP markers distributed on 22 linkage groups and total size of 1130.63 cM. Lateral stripes are inherited as a Mendelian trait and map to a single genomic interval that harbours a paralog of a gene with known function in stripe patterning. Dorsolateral and mid-lateral stripes were always coinherited and are thus under the same genetic control. Additionally, we directly quantify the genotyping error rates in RAD markers and offer guidelines for identifying and dealing with errors. Uncritical marker selection was found to severely impact linkage map construction. Fortunately, by applying appropriate quality control steps, a genotyping accuracy of >99.9% can be reached, thus allowing for efficient linkage mapping of evolutionarily relevant traits.

Genomics and the origin of species

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics.

Genomic architecture of ecologically divergent body shape in a pair of sympatric crater lake cichlid fishes

Determining the genetic bases of adaptations and their roles in speciation is a prominent issue in evolutionary biology. Cichlid fish species flocks are a prime example of recent rapid radiations, often associated with adaptive phenotypic divergence from a common ancestor within a short period of time. In several radiations of freshwater fishes, divergence in ecomorphological traits - including body shape, colour, lips and jaws - is thought to underlie their ecological differentiation, specialization and, ultimately, speciation. The Midas cichlid species complex (Amphilophus spp.) of Nicaragua provides one of the few known examples of sympatric speciation where species have rapidly evolved different but parallel morphologies in young crater lakes. This study identified significant QTL for body shape using SNPs generated via ddRAD sequencing and geometric morphometric analyses of a cross between two ecologically and morphologically divergent, sympatric cichlid species endemic to crater Lake Apoyo: an elongated limnetic species (Amphilophus zaliosus) and a high-bodied benthic species (Amphilophus astorquii). A total of 453 genome-wide informative SNPs were identified in 240 F2 hybrids. These markers were used to construct a genetic map in which 25 linkage groups were resolved. Seventy-two segregating SNPs were linked to 11 QTL. By annotating the two most highly supported QTL-linked genomic regions, genes that might contribute to divergence in body shape along the benthic-limnetic axis in Midas cichlid sympatric adaptive radiations were identified. These results suggest that few genomic regions of large effect contribute to early stage divergence in Midas cichlids.

Shaping development through mechanical strain: the transcriptional basis of diet-induced phenotypic plasticity in a cichlid fish

Adaptive phenotypic plasticity, the ability of an organism to change its phenotype to match local environments, is increasingly recognized for its contribution to evolution. However, few empirical studies have explored the molecular basis of plastic traits. The East African cichlid fish Astatoreochromis alluaudi displays adaptive phenotypic plasticity in its pharyngeal jaw apparatus, a structure that is widely seen as an evolutionary key innovation that has contributed to the remarkable diversity of cichlid fishes. It has previously been shown that in response to different diets, the pharyngeal jaws change their size, shape and dentition: hard diets induce an adaptive robust molariform tooth phenotype with short jaws and strong internal bone structures, while soft diets induce a gracile papilliform tooth phenotype with elongated jaws and slender internal bone structures. To gain insight into the molecular underpinnings of these adaptations and enable future investigations of the role that phenotypic plasticity plays during the formation of adaptive radiations, the transcriptomes of the two divergent jaw phenotypes were examined. Our study identified a total of 187 genes whose expression differs in response to hard and soft diets, including immediate early genes, extracellular matrix genes and inflammatory factors. Transcriptome results are interpreted in light of expression of candidate genes-markers for tooth size and shape, bone cells and mechanically sensitive pathways. This study opens up new avenues of research at new levels of biological organization into the roles of phenotypic plasticity during speciation and radiation of cichlid fishes.

The evolutionary history of Xiphophorus fish and their sexually selected sword: a genome-wide approach using restriction site-associated DNA sequencing

Next-generation sequencing (NGS) techniques are now key tools in the detection of population genomic and gene expression differences in a large array of organisms. However, so far few studies have utilized such data for phylogenetic estimations. Here, we use NGS data obtained from genome-wide restriction site-associated DNA (RAD) (∼66000 SNPs) to estimate the phylogenetic relationships among all 26 species of swordtail and platyfish (genus Xiphophorus) from Central America. Past studies, both sequence and morphology-based, have differed in their inferences of the evolutionary relationships within this genus, particularly at the species-level and among monophyletic groupings. We show that using a large number of markers throughout the genome, we are able to infer the phylogenetic relationships with unparalleled resolution for this genus. The relationships among all three major clades and species within each of them are highly resolved and consistent under maximum likelihood, Bayesian inference and maximum parsimony. However, we also highlight the current cautions with this data type and analyses. This genus exhibits a particularly interesting evolutionary history where at least two species may have arisen through hybridization events. Here, we are able to infer the paternal lineages of these putative hybrid species. Using the RAD-marker-based tree we reconstruct the evolutionary history of the sexually selected sword trait and show that it may have been present in the common ancestor of the genus. Together our results highlight the outstanding capacity that RAD sequencing data has for resolving previously problematic phylogenetic relationships, particularly among relatively closely related species.

Transcriptomics of morphological color change in polychromatic Midas cichlids

Background

Animal pigmentation has received much attention in evolutionary biology research due to its strong implications for adaptation and speciation. However, apart from a few cases the genetic changes associated with these evolutionary processes remain largely unknown. The Midas cichlid fish from Central America are an ideal model system for investigating pigmentation traits that may also play a role in speciation. Most Midas cichlids maintain their melanophores and exhibit a grayish (normal) color pattern throughout their lives. A minority of individuals, however, undergo color change and exhibit a distinctive gold or even white coloration in adulthood. The ontogenetic color change in the Midas cichlids may also shed light on the molecular mechanisms underlying pigmentation disorders in humans.

Results

Here we use next-generation sequencing (Illumina) RNAseq analyses to compare skin transcriptome-wide expression levels in three distinct stages of color transformation in Midas cichlids. cDNA libraries of scale tissue, for six biological replicates of each group, were generated and sequenced using Illumina technology. Using a combination of three differential expression (DE) analyses we identified 46 candidate genes that showed DE between the color morphs. We find evidence for two key DE patterns: a) genes involved in melanosomal pathways are up-regulated in normally pigmented fish; and b) immediate early and inflammatory response genes were up-regulated in transitional fish, a response that parallels some human skin disorders such as melanoma formation and psoriasis. One of the DE genes segregates with the gold phenotype in a genetic cross and might be associated with incipient speciation in this highly “species-rich” lineage of cichlids.

Conclusions

Using transcriptomic analyses we successfully identified key expression differences between different color morphs of Midas cichlid fish. These differentially expressed genes have important implications for our understanding of the molecular mechanisms underlying speciation in this lineage of extremely young species since they mate strongly assortatively, and new species may arise by sexual selection due to this color polymorphism. Some of the human orthologues of the genes identified here may also be involved in pigmentation differences and diseases and therefore provide genetic markers for the detection of human pigmentation disorders.

Differential growth-related gene expression in abalone (Haliotis midae)

The slow growth rate of Haliotis midae impedes the optimal commercial production of this most profitable South African aquaculture species. To date, no comprehensive effort has been made to identify genes associated with growth variation in farmed H. midae. The aim of this study was therefore to investigate growth variation in H. midae and to identify and quantify the expression of selected growth-related genes. Towards this aim, molecular methodologies and cell cultures were combined as a time-efficient and economical way of studying abalone transcriptomics and cell biology. Modern Illumina sequencing-by-synthesis technology and subsequent sequence annotation were used to elucidate differential gene expression between two sibling groups of abalone demonstrating significant growth variation. The expression of selected target genes involved in growth was subsequently analysed by quantitative real-time PCR (qPCR). Fast- and slow-growing abalone and in vitro primary haemocyte cultures treated with different growth-stimulating factors were used. The results obtained from transcriptome analysis and qPCR revealed significant differences in gene expression between large and small abalone, and between treated and untreated haemocyte cell cultures. Throughout in vivo and in vitro qPCR experiments, the up-regulation of genes involved in the insulin signalling pathway suggests that insulin may be involved in enhanced growth rate for various H. midae tissues.

Observation of the decay B- → D(s)((*)+) K- ℓ- ν(ℓ)

We report the observation of the decay B- → D(s)((*)+) K- ℓ- ν(ℓ) based on 342  fb(-1) of data collected at the Υ(4S) resonance with the BABAR detector at the PEP-II e+ e- storage rings at SLAC. A simultaneous fit to three D(s)(+) decay chains is performed to extract the signal yield from measurements of the squared missing mass in the B meson decay. We observe the decay B- → D(s)((*)+) K- ℓ- ν(ℓ) with a significance greater than 5 standard deviations (including systematic uncertainties) and measure its branching fraction to be B(B- → D(s)((*)+) K- ℓ- ν(ℓ)) = [6.13(-1.03)(+1.04)(stat)±0.43(syst)±0.51(B(D(s)))]×10(-4), where the last error reflects the limited knowledge of the D(s) branching fractions.

Search for production of invisible final states in single-photon decays of Υ(1S)

We search for single-photon decays of the Υ(1S) resonance, Υ → γ + invisible, where the invisible state is either a particle of definite mass, such as a light Higgs boson A⁰, or a pair of dark matter particles, χχ. Both A⁰ and χ are assumed to have zero spin. We tag Υ(1S) decays with a dipion transition Υ(2S) → π⁺π⁻Υ(1S) and look for events with a single energetic photon and significant missing energy. We find no evidence for such processes in the mass range m(A⁰) ≤ 9.2 GeV and m(χ) ≤ 4.5 GeV in the sample of 98 × 10⁶ Υ(2S) decays collected with the BABAR detector and set stringent limits on new physics models that contain light dark matter states.