Search

It is still difficult to construct the genomes of higher organisms as their genome sequences must be extended to the length of the chromosome by linkage analysis. In this study, we attempted to provide an innovative alternative to conventional linkage analysis by devising a method to genotype sperm using 10× Genomics single-cell genome sequencing libraries to generate a linkage map without interbreeding individuals. A genome was assembled using sperm from the Japanese stickleback Gasterosteus nipponicus, with single-cell genotyping yielding 1 864 430 very dense hetero-SNPs and an average coverage per sperm cell of 0.13×. In total, 1665 sperm were used, which is an order of magnitude higher than the number of recombinations used for conventional linkage analysis. We then improved the linkage analysis tool scaffold extender with low depth linkage analysis (SELDLA) to analyze the data according to the characteristics of the single-cell genotyping data. Finally, we were able to determine the chromosomal location (97.1%) and orientation (64.4%) of the contigs in the 456 Mb genome of G. nipponicus, sequenced using nanopores. This method promises to be a useful tool for determining the genomes of non-model organisms for which breeding systems have not yet been established by linkage analysis.

Variation in the laterality often exists within species and can be maintained by frequency-dependent selection. Although the molecular developmental mechanisms underlying the left-right axis formation have been investigated, the genomic mechanisms underlying variation in laterality remain largely unknown. The scale-eating cichlid Perissodus microlepis in Lake Tanganyika exhibit lateralized predation; lefty individuals with the mouth opening toward the right preferentially attack on the prey's left trunk, while righty individuals with the opposite opening attacks on the right trunk. Here, we performed RNA-sequencing and subsequent confirmation with quantitative-PCR in the telencephalon, optic tectum, and hindbrain of the cichlid and identified five genes (pkd1b, ntn1b, ansn, pde6g, and rbp4l1) that were differentially expressed between the hemispheres regardless of the laterality. Surprisingly, pkd1b and ntn1b are involved in nodal and netrin signalling, respectively, which are important for left–right asymmetry formation during early embryogenesis. This result indicates that nodal- and netrin-related signals may also play important roles in the maintenance of asymmetry in adult brain. By contrast, no genes showed reversal of lateral differences between lefty and righty individuals in any brain regions examined, suggesting that laterality in the scale-eating cichlid does not simply result from inversion of the left–right asymmetry of gene expression.

Differential gene expression can play an important role in phenotypic evolution and divergent adaptation. Although differential gene expression can be caused by both local- and distant-regulatory changes, we know little about their relative contribution to transcriptome evolution in natural populations. Here, we conducted expression quantitative trait loci (eQTL) analysis to investigate the genetic architecture underlying transcriptome divergence between marine and stream ecotypes of threespine sticklebacks (Gasterosteus aculeatus). We identified both local and distant eQTLs, some of which constitute hotspots, regions with a disproportionate number of significant eQTLs relative to the genomic background. The majority of local eQTLs including those in the hotspots caused expression changes consistent with the direction of transcriptomic divergence between ecotypes. Genome scan analysis showed that many local eQTLs overlapped with genomic regions of high differentiation. In contrast, nearly half of the distant eQTLs including those in the hotspots caused opposite expression changes, and few overlapped with regions of high differentiation, indicating that distant eQTLs may act as a constraint of transcriptome evolution. Finally, a comparison between two salinity conditions revealed that nearly half of eQTL hotspots were environment specific, suggesting that analysis of genetic architecture in multiple conditions is essential for predicting response to selection.

Intragenomic conflict, the conflict of interest between different genomic regions within an individual, is proposed as a mechanism driving both the rapid evolution of heterochromatin-related proteins and the establishment of intrinsic genomic incompatibility between species. Although molecular studies of laboratory model organisms have demonstrated the link between heterochromatin evolution and hybrid abnormalities, we know little about their link in natural systems. Previously, we showed that F1 hybrids between the Japan Sea stickleback and the Pacific Ocean stickleback show hybrid male sterility and found a region responsible for hybrid male sterility on the X chromosome, but did not identify any candidate genes. In this study, we first screened for genes rapidly evolving under positive selection during the speciation of Japanese sticklebacks to find genes possibly involved in intragenomic conflict. We found that the region responsible for hybrid male sterility contains a rapidly evolving gene encoding a heterochromatin-binding protein TRIM24B. We conducted biochemical experiments and showed that the binding affinity of TRIM24B to a heterochromatin mark found at centromeres and transposons, histone H4 lysine 20 trimethylation (H4K20me3), is reduced in the Japan Sea stickleback. In addition, mRNA expression levels of Trim24b were different between the Japan Sea and the Pacific Ocean testes. Further expression analysis of genes possibly in the TRIM24B-regulated pathway showed that some gypsy retrotransposons are overexpressed in the F1 hybrid testes. We, therefore, demonstrate that a heterochromatin-binding protein can evolve rapidly under positive selection and functionally diverge during stickleback speciation.

The Oryzias woworae species group, composed of O. asinua, O. wolasi, and O. woworae, is widely distributed in southeastern Sulawesi, an island in the Indo-Australian Archipelago. Deep-elongated body shape divergence is evident among these three species to the extent that it is used as a species-diagnostic character. These fishes inhabit a variety of habitats, ranging from upper streams to ponds, suggesting that the body shape divergence among the three species may reflect adaptation to local environments. First, our geometric morphometrics among eight local populations of this species group revealed that the three species cannot be separated by body shape and that riverine populations had more elongated bodies and longer caudal parts than lacustrine populations. Second, their phylogenetic relationships did not support the presence of three species; phylogenies using mitochondrial DNA and genomic data obtained from RNA-Seq revealed that the eight populations could not be sorted into three different clades representing three described species. Third, phylogenetic corrections of body shape variations and ancestral state reconstruction of body shapes demonstrated that body shape divergence between riverine and lacustrine populations persisted even if the phylogenies were considered and that body shape evolved rapidly irrespective of phylogeny. Sexual dimorphism in body shape was also evident, but the degree of dimorphism did not significantly differ between riverine and lacustrine populations after phylogenetic corrections, suggesting that sexual selection may not substantially contribute to geographical variations in body shape. Overall, these results indicate that the deep-elongated body shape divergence of the O. woworae species group evolved locally in response to habitat environments, such as water currents, and that a thorough taxonomic reexamination of the O. woworae species group may be necessary.

Adaptation to different salinities can drive and maintain divergence between populations of aquatic organisms. Anadromous and stream ecotypes of threespine stickleback (Gasterosteus aculeatus) are an excellent model to explore the genetic mechanisms underlying osmoregulation divergence. Using a parapatric pair of anadromous and stream stickleback ecotypes, we employed an integrated genomic approach to identify candidate genes important for adaptation to different salinity environments. Quantitative trait loci (QTL) mapping of plasma sodium concentrations under a seawater challenge experiment identified a significant QTL on chromosome 16. To identify candidate genes within this QTL, we first conducted RNA-seq and microarray analysis on gill tissue to find ecotypic differences in gene expression that were associated with plasma Na+ levels. This resulted in the identification of ten candidate genes. Quantitative PCR analysis on gill tissue of additional Japanese stickleback populations revealed that the majority of the candidate genes showed parallel divergence in expression levels. Second, we conducted whole-genome sequencing and found five genes that are predicted to have functionally important amino acid substitutions. Finally, we conducted genome scan analysis and found that eight of these candidate genes were located in genomic islands of high differentiation, suggesting that they may be under divergent selection. The candidate genes included those involved in ATP synthesis and hormonal signalling, whose expression or amino acid changes may underlie the variation in salinity tolerance. Further functional molecular analysis of these genes will reveal the causative genetic and genomic changes underlying divergent adaptation.

Aphids exhibit several polyphenisms in which discontinuous, alternative phenotypes are produced depending on environmental conditions. One representative example is the wing polyphenism, where winged and wingless females are produced through parthenogenesis. Previous work has shown that, in some aphid species, the density condition sensed by the mother aphid determines the developmental fate of embryos in her ovary, with high densities leading to winged progeny and low densities to wingless progeny. However, little is known about the molecular and physiological mechanisms underlying the wing polyphenism. To identify genes involved in the wing-morph determination in the vetch aphid, Megoura crassicauda, we compared maternal and embryonic transcripts between high- and low-density conditions using differential display, followed by quantitative real-time PCR (qRT-PCR). Under the high-density condition, two genes (Uba1 and Naca) were found to be upregulated in maternal tissues without ovaries, while one gene (ClpP) was upregulated in ovaries containing embryos. Uba1 and Naca encode factors that function in protein modification or transcriptional/translational regulation, respectively. In addition to differential display, candidate gene approaches focusing on morphogenetic and endocrine genes, i.e., wg, dpp, ap, hh, InR, IRS, Foxo, EcR, and USP, were also carried out. We found that wg was upregulated in maternal tissues under the high-density condition. The identified genes from both approaches are candidates for further study of their involvement in the transduction of density signals in mother aphids and/or the initial process of wing differentiation in embryos. J. Exp. Zool. 317:194–203, 2012. © 2012 Wiley Periodicals, Inc.

Most aphids show reproductive polyphenism, i.e. they alternate their reproductive modes from parthenogenesis to sexual reproduction in response to short photoperiods. Although juvenile hormone (JH) has been considered a likely candidate for regulating the transition from asexual to sexual reproduction after photoperiod sensing, there are few studies investigating the direct relationship between JH titres and the reproductive-mode change. In addition, the sequencing of the pea aphid genome has allowed identification of the genes involved in the JH pathway, which in turn allows us to examine their expression levels in relation to the reproductive-mode change. Using liquid chromatography-mass spectrometry in the pea aphid, JHIII titre was shown to be lower in aphids producing sexual morphs under short-day conditions than in aphids producing parthenogenetic morphs under long-day conditions. The expression levels of genes upstream and downstream of JH action were quantified by real-time quantitative reverse-transcription-PCR across the reproductive-mode change. The expression level of JH esterase, which is responsible for JH degradation, was significantly higher in aphids reared under short-day conditions. This suggests that the upregulation of the JH degradation pathway may be responsible for the lower JHIII titre in aphids exposed to short-days, leading to the production of sexual morphs.

Little is known about when, how or even if the wing development gene network elucidated in Drosophila is deployed in direct-developing insects. Here we identify the wing development genes (as determined in Drosophila) of the pea aphid (Acyrthosiphon pisum), which produces winged or unwinged adults in response to environmental cues. We find that the principal wing development genes studied in Drosophila are present in the aphid genome and that apterous and decapentaplegic exhibit duplications. We followed expression levels of 11 of these developmental genes at embryogenesis and across the nymphal instars. Six showed significant stage-specific expression level effects and apterous1 exhibited significantly different expression levels between winged and unwinged morphs, suggesting this gene acts proximately to realize polyphenic development.

The nutritional diversity of resources can affect the adaptive evolution of consumer metabolism and consumer diversification. The omega-3 long-chain polyunsaturated fatty acids eicosapentaenoic acid (EPA; 20:5n-3) and docosahexaenoic acid (DHA; 22:6n-3) have a high potential to affect consumer fitness, through their widespread effects on reproduction, growth and survival. However, few studies consider the evolution of fatty acid metabolism within an ecological context. In this review, we first document the extensive diversity in both primary producer and consumer fatty acid distributions amongst major ecosystems, between habitats and amongst species within habitats. We highlight some of the key nutritional contrasts that can shape behavioural and/or metabolic adaptation in consumers, discussing how consumers can evolve in response to the spatial, seasonal and community-level variation of resource quality. We propose a hierarchical trait-based approach for studying the evolution of consumers’ metabolic networks and review the evolutionary genetic mechanisms underpinning consumer adaptation to EPA and DHA distributions. In doing so, we consider how the metabolic traits of consumers are hierarchically structured, from cell membrane function to maternal investment, and have strongly environment-dependent expression. Finally, we conclude with an outlook on how studying the metabolic adaptation of consumers within the context of nutritional landscapes can open up new opportunities for understanding evolutionary diversification., Ecology Letters, 24 (8), ISSN:1461-023X, ISSN:1461-0248

Background One of the key characters of social insects is the division of labor, in which different tasks are allocated to various castes. In termites, one of the representative groups of social insects, morphological differences as well as behavioral differences can be recognized among castes. However, very little is known about the neuronal and molecular bases of caste differentiation and caste-specific behavior. In almost all termite species, soldiers play defensive roles in their colonies, and their morphology and behavior are largely different from workers (or pseudergates). Therefore, we predicted that some genes linked to defensive behavior and/or those required for neuronal changes are differentially expressed between workers and soldiers, or during the soldier differentiation, respectively. Results Using the brain and suboesophageal ganglion (SOG) of the damp-wood termite Hodotermopsis sjostedti, we first screened genes specifically expressed in soldiers or during soldier differentiation by the differential display method, followed by quantitative real-time polymerase chain reaction. No distinctive differences in expression patterns were detected between pseudergates and soldiers. In the course of soldier differentiation, however, five genes were found to be up-regulated in brain and/or SOG: 14-3-3epsilon, fibrillin2, beta-tubulin, ciboulot, and a hypothetical protein containing a SAP motif. Some of these genes are thought to be associated with cytoskeletal structure or motor-associated proteins in neuronal tissues. Conclusion The identified five genes could be involved in soldier-specific neuronal modifications, resulting in defensive behaviors in termite soldiers. The temporal expression patterns of these genes were consistent with the neuronal changes during soldier differentiation, suggesting that molecular machineries, in which the identified factors would participate, play important roles in behavioral differentiation of termite soldiers.