Cis-regulatory polymorphism at fiz ecdysone oxidase contributes to polygenic evolutionary response to malnutrition in Drosophila

We investigate the contribution of a candidate gene, fiz (fezzik), to complex polygenic adaptation to juvenile malnutrition in Drosophila melanogaster. Experimental populations maintained for >250 generations of experimental evolution to a nutritionally poor larval diet (Selected populations) evolved several-fold lower fiz expression compared to unselected Control populations. Here we show that this divergence in fiz expression is mediated by a cis-regulatory polymorphism. This polymorphism, originally sampled from a natural population in Switzerland, is distinct from a second cis-regulatory SNP previously identified in non-African D. melanogaster populations, implying that two independent cis-regulatory variants promoting high fiz expression segregate in non-African populations. Enzymatic analyses of Fiz protein expressed in E. coli demonstrate that it has ecdysone oxidase activity acting on both ecdysone and 20-hydroxyecdysone. Four of five fiz paralogs annotated to ecdysteroid metabolism also show reduced expression in Selected larvae, implying that malnutrition-driven selection favored general downregulation of ecdysone oxidases. Finally, as an independent test of the role of fiz in poor diet adaptation, we show that fiz knockdown by RNAi results in faster larval growth on the poor diet, but at the cost of greatly reduced survival. These results imply that downregulation of fiz in Selected populations was favored by selection on the nutritionally poor diet because of its role in suppressing growth in response to nutrient shortage. However, they suggest that fiz downregulation is only adaptive in combination with other changes evolved by Selected populations, which ensure that the organism can sustain the faster growth promoted by fiz downregulation.

Transposable Elements: Distribution, Polymorphism, and Climate Adaptation in Populus

Transposable elements (TEs) are a class of mobile genetic elements that make effects on shaping rapid phenotypic traits of adaptive significance. TE insertions are usually related to transcription changes of nearby genes, and thus may be subjected to purifying selection. Based on the available genome resources of Populus, we found that the composition of Helitron DNA family were highly variable and could directly influence the transcription of nearby gene expression, which are involving in stress-responsive, programmed cell death, and apoptosis pathway. Next, we indicated TEs are highly enriched in Populus trichocarpa compared with three other congeneric poplar species, especially located at untranslated regions (3'UTRs and 5'UTRs) and Helitron transposons, particularly 24-nt siRNA-targeted, are significantly associated with reduced gene expression. Additionally, we scanned a representative resequenced Populus tomentosa population, and identified 9,680 polymorphic TEs loci. More importantly, we identified a Helitron transposon located at the 3'UTR, which could reduce WRKY18 expression level. Our results highlight the importance of TE insertion events, which could regulate gene expression and drive adaptive phenotypic variation in Populus.

The Landscapes of Full-Length Transcripts and Splice Isoforms as Well as Transposons Exonization in the Lepidopteran Model System, Bombyx mori

The domesticated silkworm, Bombyx mori, is an important model system for the order Lepidoptera. Currently, based on third-generation sequencing, the chromosome-level genome of Bombyx mori has been released. However, its transcripts were mainly assembled by using short reads of second-generation sequencing and expressed sequence tags which cannot explain the transcript profile accurately. Here, we used PacBio Iso-Seq technology to investigate the transcripts from 45 developmental stages of Bombyx mori. We obtained 25,970 non-redundant high-quality consensus isoforms capturing ∼60% of previous reported RNAs, 15,431 (∼47%) novel transcripts, and identified 7,253 long non-coding RNA (lncRNA) with a large proportion of novel lncRNA (∼56%). In addition, we found that transposable elements (TEs) exonization account for 11,671 (∼45%) transcripts including 5,980 protein-coding transcripts (∼32%) and 5,691 lncRNAs (∼79%). Overall, our results expand the silkworm transcripts and have general implications to understand the interaction between TEs and their host genes. These transcripts resource will promote functional studies of genes and lncRNAs as well as TEs in the silkworm.

Identification and effect of Zf-AD-containing C2H2 zinc finger genes on BmNPV replication in the silkworm (Bombyx mori)

Zf-AD-containing C2H2 zinc -finger genes (ZAD) are uniquely present and have lineage-specific expansion in arthropods. Arthropods are also the hosts of Baculoviruses. We studied the possible relationship between the lineage-specific expansion of ZAD genes and arthropod-Baculovirus co-evolution. We used the silkworm (Bombyx mori) as a model. We identified 73 ZAD genes (BmZAD) in the silkworm. Sequence-based similarity analysis showed that nine clusters involving 28 BmZADs may have undergone species-specific expansion in the silkworm. Expression pattern analysis showed that the BmZADs were divided into five groups. Group I comprised 10 genes with high expression in multiple tissues, suggesting that BmZADs may play roles in the development of various tissues. We identified six BmZADs that could be induced by the Nucleopolyhedrovirus (BmNPV). Among them, BmZAD69 expression is capable of responding to BmNPV infection, and the ZAD domain is indispensable for the function of BmZAD69 in BmNPV replication. We also detected a 3 bp deletion at 1.7 kb upstream of BmZAD69, which may make it more sensitive to BmNPV infection, and thus elevate the BmNPV resistance in Qiufeng_N, a strain with strong virus resistance. These data suggest that BmZADs may be involved in BmNPV infection and that ZAD genes may play a role in arthropod-Baculovirus co-evolution.

Whole-genome resequencing reveals loci under selection during silkworm improvement

Breeding or genetic improvement refers to the process of artificial selection following domestication; as such, it has had a major influence on modern agriculture and animal production. Improvement generally focuses on traits that greatly affect the economic performance. Therefore, understanding the genetic basis underlying improvement will contribute to the identification of genes controlling economic traits and will facilitate future crop and animal breeding. However, genome-wide study of the molecular basis underlying improvement remains rare. The silkworm is a unique, entirely domesticated economically important invertebrate; genetic improvement has had a huge effect on the silkworm regarding silk-related traits. Herein, we performed whole-genomic sequencing on local and genetically improved silkworm lines to identify the genomic regions under strong selection in silkworm breeding/improvement. By genomic-wide selective sweeping analysis, we identified 24 genomic regions with strong selection signals, eight of which contained 13 candidate genes underlying silkworm breeding. Interestingly, six of these genes were annotated with functions related to neural signal response. Among the six genes, BGIBMGA004050 encodes silkworm CREB-regulated_transcription_coactivator_1 (BmCRTC1), which was reported to be involved in energy-sensing pathways. These results suggested that improvement may have affected the nervous system of the silkworm. This research will provide new insights into the genetic basis underlying the genetic improvement of silkworms and possibly of other species.

Overexpression of a Malus baccata NAC Transcription Factor Gene MbNAC25 Increases Cold and Salinity Tolerance in Arabidopsis

NAC (no apical meristem (NAM), Arabidopsis thaliana transcription activation factor (ATAF1/2) and cup shaped cotyledon (CUC2)) transcription factors play crucial roles in plant development and stress responses. Nevertheless, to date, only a few reports regarding stress-related NAC genes are available in Malus baccata (L.) Borkh. In this study, the transcription factor MbNAC25 in M. baccata was isolated as a member of the plant-specific NAC family that regulates stress responses. Expression of MbNAC25 was induced by abiotic stresses such as drought, cold, high salinity and heat. The ORF of MbNAC25 is 1122 bp, encodes 373 amino acids and subcellular localization showed that MbNAC25 protein was localized in the nucleus. In addition, MbNAC25 was highly expressed in new leaves and stems using real-time PCR. To analyze the function of MbNAC25 in plants, we generated transgenic Arabidopsis plants that overexpressed MbNAC25. Under low-temperature stress (4 °C) and high-salt stress (200 mM NaCl), plants overexpressing MbNAC25 enhanced tolerance against cold and drought salinity conferring a higher survival rate than that of wild-type (WT). Correspondingly, the chlorophyll content, proline content, the activities of antioxidant enzymes superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT) were significantly increased, while malondialdehyde (MDA) content was lower. These results indicated that the overexpression of MbNAC25 in Arabidopsis plants improved the tolerance to cold and salinity stress via enhanced scavenging capability of reactive oxygen species (ROS).

Long Noncoding RNAs and Repetitive Elements: Junk or Intimate Evolutionary Partners?

Our recent ability to sequence entire genomes, along with all of their transcribed RNAs, has led to the surprising finding that only ∼1% of the human genome is used to encode proteins. This finding has led to vigorous debate over the functional importance of the transcribed but untranslated portions of the genome. Currently, scientists tend to assume coding genes are functional until proven not to be, while the opposite is true for noncoding genes. This review takes a new look at the evidence for and against widespread noncoding gene functionality. We focus in particular on long noncoding RNA (noncoding RNAs longer than 200 nucleotides) genes and their 'junk' associates, transposable elements, and satellite repeats. Taken together, the suggestion put forward is that more of this junk DNA may be functional than nonfunctional and that noncoding RNAs and transposable elements act symbiotically to drive evolution.

A non-LTR retrotransposon activates anthocyanin biosynthesis by regulating a MYB transcription factor in Capsicum annuum

The flavonoid compound anthocyanin is an important plant metabolite with nutritional and aesthetic value as well as anti-oxidative capacity. MYB transcription factors are key regulators of anthocyanin biosynthesis in plants. In pepper (Capsicum annuum), the CaAn2 gene, encoding an R2R3 MYB transcription factor, regulates anthocyanin biosynthesis. However, no functional study or structural analysis of functional and dysfunctional CaAn2 alleles has been performed. Here, to elucidate the function of CaAn2, we generated transgenic Nicotiana benthamiana and Arabidopsis thaliana plants expressing CaAn2. All of the tissues in these plants were purple. Promoter analysis of CaAn2 in purple C. annuum 'KC00134' plants revealed the insertion of a non-long terminal repeat (LTR) retrotransposon designated Ca-nLTR-A. To determine the promoter activity and functional domain of Ca-nLTR-A, various constructs carrying different domains of Ca-nLTR-A fused with GUS were transformed into N. benthamiana. Promoter analysis showed that the 3' untranslated region (UTR) of the second open reading frame of Ca-nLTR-A is responsible for CaAn2 expression in 'KC00134'. Sequence analysis of Ca-nLTR-A identified transcription factor binding sites known to regulate anthocyanin biosynthesis. This study indicates that insertion of a non-LTR retrotransposon in the promoter may activate expression of CaAn2 by recruiting transcription factors at the 3' UTR and thus provides the first example of exaptation of a non-LTR retrotransposon into a new promoter in plants.

Transposable Elements Contribute to the Adaptation of Arabidopsis thaliana

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.

The evolutionary road from wild moth to domestic silkworm

The Silk Road, which derives its name from the trade of silk produced by the domestic silkworm Bombyx mori, was an important episode in the development and interaction of human civilizations. However, the detailed history behind silkworm domestication remains ambiguous, and little is known about the underlying genetics with respect to important aspects of its domestication. Here, we reconstruct the domestication processes and identify selective sweeps by sequencing 137 representative silkworm strains. The results present an evolutionary scenario in which silkworms may have been initially domesticated in China as trimoulting lines, then subjected to independent spreads along the Silk Road that gave rise to the development of most local strains, and further improved for modern silk production in Japan and China, having descended from diverse ancestral sources. We find that genes with key roles in nitrogen and amino acid metabolism may have contributed to the promotion of silk production, and that circadian-related genes are generally selected for their adaptation. We additionally identify associations between several candidate genes and important breeding traits, thereby advancing the applicable value of our resources.

Drosophila relics hobo and hobo-MITEs transposons as raw material for new regulatory networks

Hypermutable strains of Drosophila simulans have been studied for 20 years. Several mutants were isolated and characterized, some of which had phenotypes associated with alteration in development; for example, showing ectopic legs with eyes being expressed in place of antennae. The causal agent of this hypermutability is a non-autonomous hobo-related sequence (hoboVA). Around 100 mobilizable copies of this element are present in the D. simulans genome, and these are likely mobilized by the autonomous and canonical hobo element. We have shown that hoboVA has transcription factor binding sites for the developmental genes, hunchback and even-skipped, and that this transposon is expressed in embryos, following the patterns of these genes. We suggest that hobo and hobo-related elements can be material for the emergence of new regulatory networks.

Segmental duplications: evolution and impact among the current Lepidoptera genomes

Background

Structural variation among genomes is now viewed to be as important as single nucleoid polymorphisms in influencing the phenotype and evolution of a species. Segmental duplication (SD) is defined as segments of DNA with homologous sequence.

Results

Here, we performed a systematic analysis of segmental duplications (SDs) among five lepidopteran reference genomes (Plutella xylostella, Danaus plexippus, Bombyx mori, Manduca sexta and Heliconius melpomene) to understand their potential impact on the evolution of these species. We find that the SDs content differed substantially among species, ranging from 1.2% of the genome in B. mori to 15.2% in H. melpomene. Most SDs formed very high identity (similarity higher than 90%) blocks but had very few large blocks. Comparative analysis showed that most of the SDs arose after the divergence of each linage and we found that P. xylostella and H. melpomene showed more duplications than other species, suggesting they might be able to tolerate extensive levels of variation in their genomes. Conserved ancestral and species specific SD events were assessed, revealing multiple examples of the gain, loss or maintenance of SDs over time. SDs content analysis showed that most of the genes embedded in SDs regions belonged to species-specific SDs (“Unique” SDs). Functional analysis of these genes suggested their potential roles in the lineage-specific evolution. SDs and flanking regions often contained transposable elements (TEs) and this association suggested some involvement in SDs formation. Further studies on comparison of gene expression level between SDs and non-SDs showed that the expression level of genes embedded in SDs was significantly lower, suggesting that structure changes in the genomes are involved in gene expression differences in species.

Conclusions

The results showed that most of the SDs were “unique SDs”, which originated after species formation. Functional analysis suggested that SDs might play different roles in different species. Our results provide a valuable resource beyond the genetic mutation to explore the genome structure for future Lepidoptera research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1007-y) contains supplementary material, which is available to authorized users.

Regulatory activities of transposable elements: from conflicts to benefits

Transposable elements (TEs) are a prolific source of tightly regulated, biochemically active non-coding elements, such as transcription factor-binding sites and non-coding RNAs. Many recent studies reinvigorate the idea that these elements are pervasively co-opted for the regulation of host genes. We argue that the inherent genetic properties of TEs and the conflicting relationships with their hosts facilitate their recruitment for regulatory functions in diverse genomes. We review recent findings supporting the long-standing hypothesis that the waves of TE invasions endured by organisms for eons have catalysed the evolution of gene-regulatory networks. We also discuss the challenges of dissecting and interpreting the phenotypic effect of regulatory activities encoded by TEs in health and disease.

Multiple cis-acting elements involved in up-regulation of a cytochrome P450 gene conferring resistance to deltamethrin in smal brown planthopper, Laodelphax striatellus (Fallén)

As well as arising from single point mutations in binding sites or detoxifying enzymes, it is likely that insecticide resistance mechanisms are frequently controlled by multiple genetic factors, resulting in resistance being inherited as a quantitative trait. However, empirical evidence for this is still rare. Here we analyse the causes of up-regulation of CYP6FU1, a monoxygenase implicated in resistance to deltamethrin in the rice pest Laodelphax striatellus. The 5'-flanking region of this gene was cloned and sequenced from individuals of a susceptible and a resistant strain. A luminescent reporter assay was used to evaluate different 5'-flanking regions and their fragments for promoter activity. Mutations enhancing promoter activity in various fragments were characterized, singly and in combination, by site mutation recovery. Nucleotide diversity in flanking sequences was greatly reduced in deltamethrin-resistant insects compared to susceptible ones. Phylogenetic sequence analysis found that CYP6FU1 had five different types of 5'-flanking region. All five types were present in a susceptible strain but only a single type showing the highest promoter activity was present in a resistant strain. Four cis-acting elements were identified whose influence on up-regulation was much more pronounced in combination than when present singly. Of these, two were new transcription factor (TF) binding sites produced by mutations, another one was also a new TF binding site alternated from an existing one, and the fourth was a unique transcription start site. These results demonstrate that multiple cis-acting elements are involved in up-regulating CYP6FU1 to generate a resistance phenotype.

Homeodomain Protein Scr Regulates the Transcription of Genes Involved in Juvenile Hormone Biosynthesis in the Silkworm

The silkworm Dominant trimolting (Moltinism, M³) mutant undergoes three larval molts and exhibits precocious metamorphosis. In this study, we found that compared with the wild-type (WT) that undergoes four larval molts, both the juvenile hormone (JH) concentration and the expression of the JH-responsive gene Krüppel homolog 1 (Kr-h1) began to be greater in the second instar of the M³ mutant. A positional cloning analysis revealed that only the homeodomain transcription factor gene Sex combs reduced (Scr) is located in the genomic region that is tightly linked to the M³ locus. The expression level of the Scr gene in the brain-corpora cardiaca-corpora allata (Br-CC-CA) complex, which controls the synthesis of JH, was very low in the final larval instar of both the M³ and WT larvae, and exhibited a positive correlation with JH titer changes. Importantly, luciferase reporter analysis and electrophoretic mobility shift assay (EMSA) demonstrated that the Scr protein could promote the transcription of genes involved in JH biosynthesis by directly binding to the cis-regulatory elements (CREs) of homeodomain protein on their promoters. These results conclude that the homeodomain protein Scr is transcriptionally involved in the regulation of JH biosynthesis in the silkworm.

Identification, Diversity and Evolution of MITEs in the Genomes of Microsporidian Nosema Parasites

Miniature inverted-repeat transposable elements (MITEs) are short, non-autonomous DNA transposons, which are widespread in most eukaryotic genomes. However, genome-wide identification, origin and evolution of MITEs remain largely obscure in microsporidia. In this study, we investigated structural features for de novo identification of MITEs in genomes of silkworm microsporidia Nosema bombycis and Nosema antheraeae, as well as a honeybee microsporidia Nosema ceranae. A total of 1490, 149 and 83 MITE-related sequences from 89, 17 and five families, respectively, were found in the genomes of the above-mentioned species. Species-specific MITEs are predominant in each genome of microsporidian Nosema, with the exception of three MITE families that were shared by N. bombycis and N. antheraeae. One or multiple rounds of amplification occurred for MITEs in N. bombycis after divergence between N. bombycis and the other two species, suggesting that the more abundant families in N. bombycis could be attributed to the recent amplification of new MITEs. Significantly, some MITEs that inserted into the homologous protein-coding region of N. bombycis were recruited as introns, indicating that gene expansion occurred during the evolution of microsporidia. NbS31 and NbS24 had polymorphisms in different geographical strains of N. bombycis, indicating that they could still be active. In addition, several small RNAs in the MITEs in N. bombycis are mainly produced from both ends of the MITEs sequence.

Exploring the Phenotypic Space and the Evolutionary History of a Natural Mutation in Drosophila melanogaster

A major challenge of modern Biology is elucidating the functional consequences of natural mutations. Although we have a good understanding of the effects of laboratory-induced mutations on the molecular- and organismal-level phenotypes, the study of natural mutations has lagged behind. In this work, we explore the phenotypic space and the evolutionary history of a previously identified adaptive transposable element insertion. We first combined several tests that capture different signatures of selection to show that there is evidence of positive selection in the regions flanking FBti0019386 insertion. We then explored several phenotypes related to known phenotypic effects of nearby genes, and having plausible connections to fitness variation in nature. We found that flies with FBti0019386 insertion had a shorter developmental time and were more sensitive to stress, which are likely to be the adaptive effect and the cost of selection of this mutation, respectively. Interestingly, these phenotypic effects are not consistent with a role of FBti0019386 in temperate adaptation as has been previously suggested. Indeed, a global analysis of the population frequency of FBti0019386 showed that climatic variables explain well the FBti0019386 frequency patterns only in Australia. Finally, although FBti0019386 insertion could be inducing the formation of heterochromatin by recruiting HP1a (Heterochromatin Protein 1a) protein, the insertion is associated with upregulation of sra in adult females. Overall, our integrative approach allowed us to shed light on the evolutionary history, the relevant fitness effects, and the likely molecular mechanisms of an adaptive mutation and highlights the complexity of natural genetic variants.

Transposable elements as agents of rapid adaptation may explain the genetic paradox of invasive species

Rapid adaptation of invasive species to novel habitats has puzzled evolutionary biologists for decades, especially as this often occurs in the face of limited genetic variability. Although some ecological traits common to invasive species have been identified, little is known about the possible genomic/genetic mechanisms that may underlie their success. A common scenario in many introductions is that small founder population sizes will often lead to reduced genetic diversity, but that invading populations experience large environmental perturbations, such as changes in habitat and environmental stress. Although sudden and intense stress is usually considered in a negative context, these perturbations may actually facilitate rapid adaptation by affecting genome structure, organization and function via interactions with transposable elements (TEs), especially in populations with low genetic diversity. Stress-induced changes in TE activity can alter gene action and can promote structural variation that may facilitate the rapid adaptation observed in new environments. We focus here on the adaptive potential of TEs in relation to invasive species and highlight their role as powerful mutational forces that can rapidly create genetic diversity. We hypothesize that activity of transposable elements can explain rapid adaptation despite low genetic variation (the genetic paradox of invasive species), and provide a framework under which this hypothesis can be tested using recently developed and emerging genomic technologies.