BackgroundSmall heat shock proteins (sHSPs) are products of heat shock response and of other stress responses, and ubiquitous in all three domains of life, archaea, bacteria, and eukarya. They mainly function as molecular chaperones to protect proteins from being denatured in extreme conditions. Study on insect sHSPs could provide some insights into evolution of insects that have adapted to diverse niches in the world.ResultsTaking advantage of the newly assembled genome sequence, we performed a genome-wide analysis of the candidate sHSP genes in the silkworm, Bombyx mori. Based on known silkworm sHSP sequences, we identified 16 silkworm sHSP genes. Most of them are distributed on two silkworm chromosomes 5 and 27, respectively. 15 of 16 silkworm sHSPs have expression evidence. The comparative analysis of insect sHSPs from B. mori, Drosophila melanogaster, Apis mellifera, Tribolium castaneum, and Anopheles gambiae revealed that there is only one orthologous cluster whereas remaining clusters are species-specific on the phylogenetic tree. This suggested that most of sHSPs might have diverged in function across insects investigated. In addition, the data presented in this study also revealed that sHSPs in the insect orthologous cluster are highly conserved in both sequence and expression pattern. In sum, insect sHSPs show a completely different evolutionary pattern from that found in vertebrate sHSPs.ConclusionB. mori has the largest number of insect sHSP genes characterized to date, including 16 genes. The inference that most species-specific sHSPs might have diverged in function across insects investigated will help us understand the adaptability of these insects to diverse environments.

BackgroundOrganisms need to adapt to keep pace with a changing environment. Examining recent range expansion aids our understanding of how organisms evolve to overcome environmental constraints. However, how organisms adapt to climate changes is a crucial biological question that is still largely unanswered. The plant Arabidopsis thaliana is an excellent system to study this fundamental question. Its origin is in the Iberian Peninsula and North Africa, but it has spread to the Far East, including the most south-eastern edge of its native habitats, the Yangtze River basin, where the climate is very different.ResultsWe sequenced 118 A. thaliana strains from the region surrounding the Yangtze River basin. We found that the Yangtze River basin population is a unique population and diverged about 61,409 years ago, with gene flows occurring at two different time points, followed by a population dispersion into the Yangtze River basin in the last few thousands of years. Positive selection analyses revealed that biological regulation processes, such as flowering time, immune and defense response processes could be correlated with the adaptation event. In particular, we found that the flowering time gene SVP has contributed to A. thaliana adaptation to the Yangtze River basin based on genetic mapping.Conclusions A. thaliana adapted to the Yangtze River basin habitat by promoting the onset of flowering, a finding that sheds light on how a species can adapt to locales with very different climates.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-017-1378-9) contains supplementary material, which is available to authorized users.

Rapid phenotypic changes in traits of adaptive significance are crucial for organisms to thrive in changing environments. How such phenotypic variation is achieved rapidly, despite limited genetic variation in species that experience a genetic bottleneck is unknown.Capsella rubella, an annual and inbreeding forb (Brassicaceae), is a great system for studying this basic question. Its distribution is wider than those of its congeneric species, despite an extreme genetic bottleneck event that severely diminished its genetic variation. Here, we demonstrate that transposable elements (TEs) are an important source of genetic variation that could account for its high phenotypic diversity. TEs are (i) highly enriched inC. rubellacompared with its outcrossing sister speciesCapsella grandiflora, and (ii) 4.2% of polymorphic TEs inC. rubellaare associated with variation in the expression levels of their adjacent genes. Furthermore, we show that frequent TE insertions atFLOWERING LOCUS C (FLC)in natural populations ofC. rubellacould explain 12.5% of the natural variation in flowering time, a key life history trait correlated with fitness and adaptation. In particular, we show that a recent TE insertion at the 3′ UTR ofFLCaffects mRNA stability, which results in reducing its steady-state expression levels, to promote the onset of flowering. Our results highlight that TE insertions can drive rapid phenotypic variation, which could potentially help with adaptation to changing environments in a species with limited standing genetic variation.

Transposable elements (TEs) are mobile genetic elements with very high mutation rates that play important roles in shaping genome architecture and regulating phenotypic variation. However, the extent to which TEs influence the adaptation of organisms in their natural habitats is largely unknown. Here, we scanned 201 representative resequenced genomes from the model plant Arabidopsis thaliana and identified 2,311 polymorphic TEs from noncentromeric regions. We found expansion and contraction of different types of TEs in different A. thaliana populations. More importantly, we identified two TE insertions that are likely candidates to play a role in adaptive evolution. Our results highlight the importance of variations in TEs for the adaptation of plants in general in the context of rapid global climate change.

BackgroundHorizontal gene transfer (HGT), a source of genetic variation, is generally considered to facilitate hosts' adaptability to environments. However, convincing evidence supporting the significant contribution of the transferred genes to the evolution of metazoan recipients is rare.ResultsIn this study, based on sequence data accumulated to date, we used a unified method consisting of similarity search and phylogenetic analysis to detect horizontally transferred genes (HTGs) between prokaryotes and five insect species including Drosophila melanogaster, Anopheles gambiae, Bombyx mori, Tribolium castaneum and Apis mellifera. Unexpectedly, the candidate HTGs were not detected in D. melanogaster, An. gambiae and T. castaneum, and 79 genes in Ap. mellifera sieved by the same method were considered as contamination based on other information. Consequently, 14 types of 22 HTGs were detected only in the silkworm. Additionally, 13 types of the detected silkworm HTGs share homologous sequences in species of other Lepidopteran superfamilies, suggesting that the majority of these HTGs were derived from ancient transfer events before the radiation of Ditrysia clade. On the basis of phylogenetic topologies and BLAST search results, donor bacteria of these genes were inferred, respectively. At least half of the predicted donor organisms may be entomopathogenic bacteria. The predicted biochemical functions of these genes include four categories: glycosyl hydrolase family, oxidoreductase family, amino acid metabolism, and others.ConclusionsThe products of HTGs detected in this study may take part in comprehensive physiological metabolism. These genes potentially contributed to functional innovation and adaptability of Lepidopteran hosts in their ancient lineages associated with the diversification of angiosperms. Importantly, our results imply that pathogens may be advantageous to the subsistence and prosperity of hosts through effective HGT events at a large evolutionary scale.

De novo genes, which originate from ancestral nongenic sequences, are one of the most important sources of protein-coding genes. This origination process is crucial for the adaptation of organisms. However, how de novo genes arise and become fixed in a population or species remains largely unknown. Here, we identified 782 de novo genes from the model plant Arabidopsis thaliana and divided them into three types based on the availability of translational evidence, transcriptional evidence, and neither transcriptional nor translational evidence for their origin. Importantly, by integrating multiple types of omics data, including data from genomes, epigenomes, transcriptomes, and translatomes, we found that epigenetic modifications (DNA methylation and histone modification) play an important role in the origination process of de novo genes. Intriguingly, using the transcriptomes and methylomes from the same population of 84 accessions, we found that de novo genes that are transcribed in approximately half of the total accessions within the population are highly methylated, with lower levels of transcription than those transcribed at other frequencies within the population. We hypothesized that, during the origin of de novo gene alleles, those neutralized to low expression states via DNA methylation have relatively high probabilities of spreading and becoming fixed in a population. Our results highlight the process underlying the origin of de novo genes at the population level, as well as the importance of DNA methylation in this process.

BackgroundIn contrast to positive selection, which reduces genetic variation by fixing beneficial alleles, balancing selection maintains genetic variation within a population or species and plays crucial roles in adaptation in diverse organisms. However, which genes, genome-wide, are under balancing selection and the extent to which these genes are involved in adaptation are largely unknown.ResultsWe performed a genome-wide scan for genes under balancing selection across two plant species, Arabidopsis thaliana and its relative Capsella rubella, which diverged about 8 million generations ago. Among hundreds of genes with shared coding-region polymorphisms, we find evidence for long-term balancing selection in five genes: AT1G35220, AT2G16570, AT4G29360, AT5G38460, and AT5G44000. These genes are involved in the response to biotic and abiotic stress and other fundamental biochemical processes. More intriguingly, for these genes, we detected significant ecological diversification between the two haplotype groups, suggesting that balancing selection has been very important for adaptation.ConclusionsOur results indicate that beyond the well-known S-locus genes and resistance genes, many loci are under balancing selection. These genes are mostly correlated with resistance to stress or other fundamental functions and likely play a more important role in adaptation to diverse habitats than previously thought.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-017-1342-8) contains supplementary material, which is available to authorized users.

The position at which the second gold(i)-phosphine group is attached was experimentally found to play a noticeable role in intersystem crossing rates of gold(i) naphthalene derivatives. However, the physical origin is ambiguous. Herein we have employed generalized trajectory-based surface-hopping dynamics simulations to simulate the excited-state relaxation dynamics of these gold(i) naphthalene compounds including both the intersystem crossing process from the initially populated first excited singlet states S to triplet manifolds and internal conversion processes within these triplet states. Our predicted intersystem crossing rates are consistent with experiments very well. On the basis of the present results, we have found that (1) ultrafast and subpicosecond intersystem crossing processes are mainly caused by small energy gaps and large spin-orbit couplings between S and T; (2) adding the second gold(i)-phosphine group does not increase spin-orbit couplings between S and T but decrease their values remarkably, which implies that heavy-atom effects are state-specific, not state-universal; (3) the position at which the second gold(i)-phosphine group is attached has a remarkable influence on the electronic structures of S and T and their relative energies, which affect energy gaps and spin-orbit couplings between S and T and eventually modulate intersystem crossing rates from S to T. These new insights are very useful for the design of gold-containing compounds with excellent photoluminescence properties. Finally, this work also exemplifies that different isomers of a compound could have distinct excited-state relaxation dynamics.