About the Authors:
Songzhen He
Affiliation: State Key Laboratory of Silkworm Genome Biology, Key Laboratory for Sericulture Functional Genomics and Biotechnology of Agricultural Ministry, Southwest University, Chongqing, 400715, China
Xiaoling Tong
Affiliation: State Key Laboratory of Silkworm Genome Biology, Key Laboratory for Sericulture Functional Genomics and Biotechnology of Agricultural Ministry, Southwest University, Chongqing, 400715, China
Kunpeng Lu
Affiliation: State Key Laboratory of Silkworm Genome Biology, Key Laboratory for Sericulture Functional Genomics and Biotechnology of Agricultural Ministry, Southwest University, Chongqing, 400715, China
Yaru Lu
Affiliation: State Key Laboratory of Silkworm Genome Biology, Key Laboratory for Sericulture Functional Genomics and Biotechnology of Agricultural Ministry, Southwest University, Chongqing, 400715, China
Jiangwen Luo
Affiliation: State Key Laboratory of Silkworm Genome Biology, Key Laboratory for Sericulture Functional Genomics and Biotechnology of Agricultural Ministry, Southwest University, Chongqing, 400715, China
Wenhao Yang
Affiliation: State Key Laboratory of Silkworm Genome Biology, Key Laboratory for Sericulture Functional Genomics and Biotechnology of Agricultural Ministry, Southwest University, Chongqing, 400715, China
Min Chen
Affiliation: State Key Laboratory of Silkworm Genome Biology, Key Laboratory for Sericulture Functional Genomics and Biotechnology of Agricultural Ministry, Southwest University, Chongqing, 400715, China
Min-jin Han
Affiliation: State Key Laboratory of Silkworm Genome Biology, Key Laboratory for Sericulture Functional Genomics and Biotechnology of Agricultural Ministry, Southwest University, Chongqing, 400715, China
Hai Hu
Affiliation: State Key Laboratory of Silkworm Genome Biology, Key Laboratory for Sericulture Functional Genomics and Biotechnology of Agricultural Ministry, Southwest University, Chongqing, 400715, China
Cheng Lu
Affiliation: State Key Laboratory of Silkworm Genome Biology, Key Laboratory for Sericulture Functional Genomics and Biotechnology of Agricultural Ministry, Southwest University, Chongqing, 400715, China
Fangyin Dai
* E-mail: [email protected]
Affiliation: State Key Laboratory of Silkworm Genome Biology, Key Laboratory for Sericulture Functional Genomics and Biotechnology of Agricultural Ministry, Southwest University, Chongqing, 400715, China
Introduction
Pigmentation plays an important role in the adaptive processes of insects and provides a model system for biology [1,2]. One of the most common primary pigmentation processes is melanism, which occurs throughout Insecta and is an important source of insect polymorphism [3–7]. It is caused by the accumulation of melanin [4,8], and the genes, proteins and metabolites pertaining to melanin metabolism pathway in insects have been extensively studied [8,9]. Functional defects or abnormal expression of some key genes in the melanin metabolism pathway can also lead to melanism [6,9,10]. In contrast, some studies found that some melanic morphs had no variations in the classic melanin metabolism pathway genes indicating the mechanistic diversity of melanism [7,11,12]. An example for this mechanism is the industrial melanism in the peppered moth (Biston betularia) [7,11]. These studies suggest a knowledge gap in the genetic basis and regulatory mechanism of insect melanism.
In addition to insect polymorphism, melanism also plays an important role in ecological adaption. For example, in numerous cases, melanism is associated with advantageous phenotypes, such as higher fecundity, increased pathogen resistance, higher vigor, higher resistance to UV, etc. [2,13–18]. Deciphering the molecular basis underlying the relationship between melanism and the related phenotypic-effects could complement our existing knowledge in ecological mechanisms related to melanism, and clarify potential mechanisms in adaptive evolution.
The silkworm, Bombyx mori, is the most advanced lepidopteran model with over 100 body color mutants described thus far [19,20]. Studies of body color mutants has enhanced our fundamental understanding of the silkworm pigmentation processes and the genetic mechanisms underlying insect melanism [6,10,12,21–25].
Recently, a new melanic mutant in the silk moths termed sex-controlled melanism (sml) was discovered by researchers in the Silkworm Gene Bank of Southwest University, China [26]. It is a spontaneous autosomal recessive mutant [26], and has a rare wing- and body-color variation in adult B. mori. The sml female moths exhibit the white appearance of typical wild-types and are without obvious abnormal phenotypes, while the sml male moths exhibit the melanic phenotype, and interestingly also show extended adult longevity as well as higher vigor than typical wild-types. Our previous hereditary analysis demonstrated that this melanic phenotype was independent of the genetic variations, which were reported to be associated with melanism [26]. This suggests that an unknown novel gene likely plays a role in creating this melanic morph. Thus, we considered sml as a superior model to study the molecular mechanism of melanism and the phenotypic-effects associated with this melanic phenotype.
To obtain a comprehensive view of gene expression related to melanism and the interesting phenotypic-effects in sml males, in the present study, we used a transcriptome approach and investigated the differentially expressed genes (DEGs) between sml males and sml females as well as a wild-type strain (Dazao) in the wings and integuments at the initial stage of melanin formation and deposition. We identified 59 unique DEGs in the melanic mutant, and detected important genes related to the melanic morph as well as the related phenotypes in the sml male moth. Our findings would facilitate the understanding of the molecular basis underlying melanism and the related phenotypic-effects.
Materials and Methods
Silkworm Strains
The sex-controlled melanism mutant (sml) and wild-type strain Dazao, were obtained from the Silkworm Gene Bank in Southwest University (Chongqing, China) and were reared on mulberry leaves under a 12 h light/12 h dark photoperiod at 25°C and 75% relative humidity. Under these conditions, the pupal stage lasted about 9 days.
Statistical Analysis of Survival Time
To assay survival time, silkworm strains sml, N4, Nis. and Dazao were reared under the conditions mentioned above. The survival of unmated moths was recorded every 6 h. Data were statistically analyzed by Log-rank (Mantel-Cox) test using GraphPad Prism 5 software. P-values less than 0.05 were considered to be significant.
RNA Preparation and Illumina Sequencing
Wings and integuments without pupal case were dissected within the first 5 hours on the 8th day of pupal stage, which represents the initial stage of melanin formation and deposition. For each sample, the mixture of wings and integuments from five individuals were pooled as one tissue sample. Total RNA was isolated from the four samples using TRIzol Reagent (Invitrogen, USA) according to manufacturer’s instructions. RNA contamination and degradation were verified by separating RNA on a 1% agarose gel. RNA purity and concentration were assessed using a Nanodrop ND-2000C spectrophotometer (Thermo, USA). RNA integrity was checked using the Agilent Bioanalyzer 2100 system (Agilent Technologies, USA) according to the manufacturer’s instructions. Briefly, mRNA was purified from 3 μg total RNA per sample using oligo (dT) magnetic beads. Transcriptome libraries were generated using standard protocols as previously reported [27,28]. The serial numbers were as follows: sml males (with the serial number as sml_M), sml females (with the serial number as sml_F), Dazao males (with the serial number as DZ_M) and Dazao females (with the serial number as DZ_F). The libraries were sequenced on an Illumina Hiseq2000 system (Illumina, USA) and 125-bp paired-end reads were generated. Transcriptome library construction and sequencing were performed by NovoGene Corporation (Beijing, China). Raw data are deposited in the NCBI Short Read Archive (http://www.ncbi.nlm.nih.gov/sra/) under the SRA accession number SRR3068622.
Data Analysis
Raw reads from the Illumina sequencing in FASTQ format were processed using in-house Perl scripts to remove reads containing poly-N, reads containing adapter and low quality reads. All successive analyses were carried out using clean data with high quality.
First, the GC content, Q20, Q30, and sequence duplication level of the clean sequences were calculated. Then, the clean data were mapped to the silkworm reference genome (Bombyx mori, ftp://ftp.ensemblgenomes.org/pub/release-20/metazoa/fasta/bombyx_mori/dna/) using TopHat v2.0.12 [29]. Cufflinks v2.1.1 [30] was used to identify and construct novel and known transcripts from the TopHat mapping results. HTSeq v0.6.1 [31] was used to count the reads mapped to each transcript and reads per kilo-base per million mapped reads (RPKM) was used to quantify transcript expression. The read counts were adjusted using edgeR 3.0.8 [32] through one scaling normalized factor. Differential expression analysis was performed using DEGSeq 1.12.0 [33]. The P-values were adjusted by the Benjamini-Hochberg method [34] for multiple testing. An adjusted P-value of 0.005 and log2 (Fold-change) of ±1 were set as the threshold to determine statistically significant differential expression.
GOseq [35] was used to perform Gene Ontology (GO) enrichment analysis of differentially expressed transcripts. GO terms with adjusted P-value less than 0.05 were considered significantly enriched transcripts. KOBAS v2.0 [36] was used to test the statistical enrichment of differentially expressed genes in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. KEGG terms with corrected P-value less than 0.05 were considered as significantly enriched genes.
Real-Time qRT-PCR
By random selection, 20 genes with differential expression were chosen for verification by quantitative RT-PCR (qRT-PCR), which was performed on a CFX96 Touch Real-Time PCR System (Bio-Rad, USA) with iTaq™ Universal SYBR® Green SuperMix (Bio-Rad, USA). Gene expression levels were normalized to the expression of B. mori eukaryotic translation initiation factor 4A (microarray ID: sw22934). Relative expression was calculated using the 2-ΔΔCt method [37]. Primer sequences used for qRT-PCR are described in S7 Table.
Results
sml Male Begins to Melanize in Late Pupal Stage Due to Melanic Body Color of the Adult
Homozygous sml male moths exhibit the melanic phenotype with grayish-black scales on the wings and body [26], while females exhibit the white appearance of typical wild-types (Fig 1).
[Figure omitted. See PDF.]
Fig 1. Phenotypes of the sml mutant and wild type (Dazao).
Visible differences were not observed before P8. (A) Phenotypes of sml mutation were investigated from pupation to the 1st day after eclosion. Pigmentation in sml males was initially observed in P8. Pupae in P9 and moths of sml males exhibited melanic phenotypes, while females were normal. (B) Phenotype of wild type (Dazao) from the 1st day of pupae to the 1st day after eclosion. P1-P9 represents the 1st day to the 9th day of pupae; M represents the 1st day after eclosion; ♂indicates male and ♀ indicates female.
https://doi.org/10.1371/journal.pone.0155061.g001
To compare the pigmentation in sml and wild-type, phenotypes were investigated every day from pupation to the 1st day after eclosion (Fig 1). Visible differences were not observed before P8 (the 8th day of pupae). Pigmentation in sml males were initially observed by the naked eye in P8, suggesting that melanin formation and accumulation is initiated as early as P8. Further investigation revealed that in sml males the colors of other sclerotic tissues such as pupal cases, moth's dorsal plate and tentacles were not visibly different from sml females or wild-type individuals (S1 Fig), indicating that the phenotype of sml male was due to melanic body color of the adult.
The sml Male Moth Exhibits Significantly Extended Adult Longevity
In most cases, melanism plays an important role in ecological adaption and endows selective advantage on insects [3]. The sml female moths with the typical white appearance show normal phenotypes (S2 Fig), including normal oviposition amount, adult longevity, and so on. However, as previous reports of melanic insects [16,38,39], we found that the sml male moths with a melanic phenotype were more active than wild-type male moths. Interestingly, the adult longevity of sml male moths was also significantly longer than wild-type males, and even longer than the females (Fig 2A). In contrast, the lifespan of wild-type female silkworms was conspicuously longer than the males [40] (Fig 2B, 2C and 2D) consistent with the longer lifespan of females throughout the animal kingdom. These findings revealed that the melanic individuals in sml strain were more viable suggesting that melanism could extend the adult longevity of insect under some conditions. Thus, sml mutants could be a superior model to study the molecular mechanisms of melanism and the associated phenotypic-effects.
[Figure omitted. See PDF.]
Fig 2. Survival time of unmated silk-moth.
The survival time of sml male moths is significantly longer than sml female moths, and is in contrast to the trend in typical silk-moths. (A) The survival time of sml male moths (n = 64) is significantly longer than females (n = 58). (B) Survival of Nis. male moths (n = 45) is significantly shorter than females (n = 34). (C) Survival of N4 male moths (n = 79) is significantly shorter than females (n = 101). (D) The survival time of Dazao male moths (n = 65) is significantly shorter than females (n = 58). P-values less than 0.05 were considered significant. ♂ indicates male and ♀ indicates female.
https://doi.org/10.1371/journal.pone.0155061.g002
Sequencing and Assembly
To investigate global transcriptional changes associated with this melanic phenotype and the associated phenotypic-effects in sml male, both male and female samples from sml and Dazao strains at the initial stage of melanin formation and deposition (P8) were selected for comparative analysis. Four transcriptome libraries were constructed, and 52.0 to 63.5 million raw reads were generated for each library (Table 1). A summary of sequencing data is outlined in Table 1. Analysis of quality scores demonstrated that the quality of data was robust. After quality control, 6.26 to 7.7 Gb clean bases (about 14 to 17-folds coverage of the B. mori genome) were obtained for each sample (Table 1), and 79.71 to 91.34% clean reads were matched to the silkworm genome (S1 Table). A total of 1,572 novel transcripts were identified according to the gene structure annotation file (ftp:/ftp.ensemblgenomes.org/pub/release-20/metazoa/gtf/bombyx_mori/) (S2 Table).
[Figure omitted. See PDF.]
Table 1. Summary of transcriptome sequencing data.
https://doi.org/10.1371/journal.pone.0155061.t001
Analysis of Differentially Expressed Genes
Differentially expressed transcripts related to the melanic morph in sml males could provide molecular insights into the mechanisms underlying melanism and the associated phenotypic-effects. After analysis, a total of 874 DEGs were identified among the four samples by comparing any two samples (S3 Table). A hierarchical clustering graph was generated to observe gene expression patterns from all the DEGs based on log10 (RPKM+1) values of the samples (Fig 3). The regular distribution of DEGs is shown in the hierarchical clustering graph with the four regions marked with boxes of different colors as follows: DEGs between different strains (light-blue boxes), unique DEGs in sml male (red box), unique DEGs in sml female (purple box) and DEGs between different genders (orange box) (Fig 3).
[Figure omitted. See PDF.]
Fig 3. Hierarchical clustering graph of DEGs.
Hierarchical clustering graph of the 874 DEGs identified among the four samples by comparing any two samples based on log10 (RPKM+1) values. Yellow bands represent high levels of gene expression; blue bands represent low levels of gene expression. Four different regions are marked by boxes with different colors. The light-blue boxes indicate DEGs between different strains; the red box indicates unique DEGs in sml male; the purple box indicates unique DEGs in sml female; the orange box indicates DEGs between different genders. DZ represents Dazao; M represents male; F represents female.
https://doi.org/10.1371/journal.pone.0155061.g003
To survey the molecular function of DEGs related to melanic morph, functional annotation was subsequently performed based on GO categories. DEGs obtained from the three pairs of melanic-typical comparisons (sml_M vs sml_F, sml_M vs DZ_M and sml_M vs DZ_F) were categorized into functional groups: biological process, cellular component and molecular function, respectively (S4 Table). GO terms that were significantly enriched (corrected P-value less than 0.05) are shown in Fig 4. The three pairwise melanic-typical comparisons shared nine common significantly enriched GO terms (Fig 4) including metabolic process (chitin metabolic process GO:0006030, amino sugar metabolic process GO:0006040, glucosamine-containing compound metabolic process GO:1901071, aminoglycan metabolic process GO:0006022), oxidation−reduction process (GO:0055114), extracellular region (GO:0005576), structural constituent of cuticle (GO:0042302) and binding (chitin binding GO:0008061 and carbohydrate derivative binding GO:0097367). Interestingly, apart from the oxidation−reduction process (GO:0055114) and the extracellular region (GO:0005576), most genes in the remaining seven common GO terms were up-regulated in the sml male (S4 Table). Such up-regulated gene expression of the seven GO terms in the melanic sml male sample suggested that these molecular functions could be closely related to the melanic morph and related phenotypes.
[Figure omitted. See PDF.]
Fig 4. GO classification of DEGs.
GO classification of DEGs obtained respectively from three pair of melanic-typical comparisons (sml_M vs sml_F, sml_M vs DZ_M and sml_M vs DZ_F). Significantly enriched GO terms are shown. Asterisks (*) indicate the nine common significantly enriched GO terms (corrected P-value less than 0.05) shared by the three pairwise melanic-typical comparisons. DZ represents Dazao; M represents male; F represents female.
https://doi.org/10.1371/journal.pone.0155061.g004
Meanwhile, the DEGs obtained from the three pairwise melanic-typical comparisons were classified into the reference pathways in the KEGG database, respectively [41] (S5 Table). The significantly enriched KEGG terms (corrected P-value less than 0.05) from this analysis are shown in S3 Fig. The carbon pool by folate pathway (bmor00670) and glycine, serine and threonine metabolism pathway (bmor00260) were significantly enriched for the DEGs between sml male vs sml female (S3 Fig, sml_M vs sml_F). Likewise, two significantly enriched pathways, fatty acid degradation (bmor00071) and citrate cycle (TCA cycle) (bmor00020) were identified in the comparisons between sml_M vs DZ_M (S3 Fig). In the comparisons between sml_M vs DZ_F, the pyruvate metabolism pathway (bmor00620) was significantly enriched (S3 Fig, sml_M vs DZ_F).
Since the four samples used in transcriptome analysis originated from different strains and different genders, the DEGs identified in the comparison between sml_M vs sml_F are likely linked to differences in gender and/or body color. Similarly, the DEGs detected in the comparison between sml_M vs DZ_M could be closely related to differences in the genetic background between strains and body color. On the other hand, the DEGs obtained in the comparisons between sml_M vs DZ_F could be likely related to the differences in genetic background among the strains, sex differences, and body color.
To eliminate the interference from genetic background and sex differences, a Venn diagram was constructed to detect the major DEGs associated with the melanic-morph. As shown in Fig 5A, there were 168 DEGs, 503 DEGs and 623 DEGs detected in the three pairwise melanic-typical comparisons; sml_M vs sml_F, sml_M vs DZ_M and sml_M vs DZ_F, respectively. Among these, 72 DEGs were common to the pairwise melanic-typical comparisons (Fig 5A). From these, 13 DEGs were excluded because they showed expression differences in the typical-typical comparison between Dazao male vs Dazao female (DZ_M vs DZ_F) and are likely related to sex differences rather than body color diversities (Fig 5B). Therefore, the remaining 59 DEGs were selected as candidates to study body color diversities in further analyses.
[Figure omitted. See PDF.]
Fig 5. Venn diagrams of DEGs.
(A) Venn diagram of DEGs in the three pairwise melanic-typical comparisons. The 72 overlapped genes were common to the three pairwise comparisons. (B) Venn diagram between the 72 overlapped genes and the DEGs in the typical-typical comparison between DZ_M vs DZ_F. DZ represents Dazao; M represents male; F represents female.
https://doi.org/10.1371/journal.pone.0155061.g005
The 59 DEGs were further analyzed manually based on functional annotation and transcriptome expression data (S3 Table). Among these DEGs, 49 were up-regulated and 10 were down-regulated in the melanic sml male sample compared to the other three typical samples (S6 Table). Intriguingly, up to 24 genes (about 41% of the candidate DEGs) were cuticular protein encoding genes (GO:0042302 and GO:0005576), including CPR (cuticular protein with the R&R Consensus), CPG (glycine-rich cuticular protein), CPT (cuticular protein with a Tweedle motif) and CPH (hypothetical cuticular protein) genes, and 22 genes were uniquely up-regulated in the sml male sample (Table 2), suggesting that these cuticular protein genes may be related to melanism. Two genes with key roles in insect melanism, i.e. Bm-laccase2 (BGIBMGA006745) and Bm-yellow (BGIBMGA001149) (GO:0055114) [42,43], were significantly up-regulated in the sml male (Table 2 Melanin metabolism). Meanwhile, the laccase2 B-type isoform transcript, BGIBMGA006746, was also significantly up-regulated in sml male (Table 2 Melanin metabolism) suggesting that this isoform could also play a role in the melanism process at late pupal stage. But, other genes involved in melanization pathway, such as pale and Ddc, didn't show expression differences in the melanic sample compared to the other three white samples (Table 2 and S6 Table). Moreover, BGIBMGA011430 (GO:0004096) and BGIBMGA000667 (GO:0005179), encoding a catalase and an insulin-like peptide (Bombyxin E-1) respectively, were also found among the candidate DEGs (Table 2). These two genes might be related to the extended life span of sml in adulthood.
[Figure omitted. See PDF.]
Table 2. Candidate DEGs of interest in the melanic sml male.
https://doi.org/10.1371/journal.pone.0155061.t002
Validation of Transcriptome Results by qRT-PCR
To validate the transcriptome data at the P8 stage, 20 DEGs were randomly selected for qRT-PCR analysis, which also showed similar down- or up-regulated trends consistent with the Illumina sequencing data (S4 Fig) indicating the reliability of the comparative analysis of our transcriptomes.
Discussion
Melanism is critical for insect polymorphism [3–7] and is also an important adaptive trait that usually imparts advantageous phenotypes, and plays important roles in natural selection and adaptive evolution [2,13–18]. However, due to the mechanistic diversity, there are still knowledge gaps in the molecular mechanisms underlying melanism and derived phenotypic-effects. In the present study, we used a transcriptome approach to perform comparative analysis between sml, a newly-discovered silkworm melanic mutant, and a wild-type strain (Dazao). Our findings revealed the important genes associated with this melanic morph as well as the derived phenotypes.
Candidate DEGs Underlying the Melanic Morph
In many insects, the substances and genes pertaining to melanin metabolism pathway has been scientifically studied so far. The synthesis of melanin begins with the hydroxylation of tyrosine [9]. Tyrosine hydroxylase (pale) converts tyrosine into dopa, and dopa decarboxylase (Ddc) catalyzes dopa into dopamine [8]. In this pathway, dopa and dopamine serve as key precursors for melanin. The final steps in the melanism process are the formation of melanin and its deposition in the cuticle. Yellow and laccase2 are thought to play major roles in these steps [42–45]. Our previous finding in sml revealed that an unknown novel gene plays a role in creating this melanic morph [26]. In the present study, we found that Bm-laccase2 and Bm-yellow were significantly up-regulated in sml male, but other genes involved in insect melanization pathway didn't show expression differences in the melanic sample compared to the other three white samples (Table 2 and S6 Table). Since the unique differentially expressed genes in the melanic mutant function toward the end of the melanin formation pathway, we presume that sml melanic alleles likely function via the melanin metabolism pathway by up-regulating the expression of yellow and laccase2 genes.
Moreover, laccase2 is thought to play important roles in insect cuticle pigmentation and sclerotization, and the laccase2 gene encodes two alternative splicing isoforms, i.e. A-type isoform and B-type isoform [46–50]. It was reported that both isoforms were able to oxidize the same endogenous substrates, including dopamine and dopa [51], suggesting that the isoforms may have similar (or the same) functions. However, there is as yet no evidence for the presence of the B-type isoform in cuticle [46–50,52]. It is noteworthy that all integument samples for analysis in previous studies of lepidopteran insects, including Bombyx mori and Manduca sexta, were gathered from larval stage, pre-pupal stage or newly-ecdysed pupae [48–50]. In this study, we found the laccase2 B-type isoform transcript, BGIBMGA006746, was not only expressed at P8 stage, but also significantly up-regulated in the melanic mutant (Table 2). According to these results, we presume that the laccase2 B-type isoform was likely involved in silkworm cuticle pigmentation and sclerotization, especially cuticular melanism process, at late pupal stage.
Previous studies found that interactions between cuticular proteins and chitin are essential for the exoskeletal structure and physical properties of the cuticle [53,54]. In addition, color patterns and exoskeletal structures, including structures of cuticle surface coat and scales, were found to be closely related [55–57]. Furthermore, the expression of cuticular protein genes varied in insects with different color patterns [57]. Futahashi et al. proposed that cuticular proteins were likely related to color patterns and that different kinds of cuticular proteins may play roles in transporting or maintaining different cuticular pigments [57]. In our present study, we found 59 candidate DEGs that are likely associated with melanism, since our screening strategy was strictly based on body color diversity (Fig 5). Among the 59 DEGs, 24 (about 41% of the candidate DEGs) were cuticular protein coding genes, and 22 of these were uniquely up-regulated in melanic sml male (Table 2 and S6 Table). Thus, we presume that these up-regulated cuticular protein genes in sml male could specifically play a role in transporting and/or maintaining melanin.
Additionally, Janssen et al. demonstrated the correlations between scale structure and scale color, and showed that black scales had denser ultra-structures [56]. Besides, the present studies showed that CPRs and CPTs, with chitin-binding domains, play important roles in the cuticle structure [53,54,58–60]. And it was also reported that CPGs could play a part in insect cuticle sclerotization by cross-linking of proteins [61]. Interestingly, more than half (about 59%) of the up-regulated cuticular protein genes in our melanic mutant encode CPRs, CPTs or CPGs (Table 2). According to these findings, we presume that the numerous up-regulated cuticular proteins which were related to exoskeleton structure and cuticle sclerotization in the black exoskeletons might result in the increased density of the black scales. Furthermore, we assume that interaction between melanin and cuticular proteins should have an important contribution for the structure and physical properties of the black scales.
Candidate DEGs Underlying the Derived Phenotypes
Previous studies showed that some melanic individuals with higher activity might have shorter adult longevity [16]. In contrast, the melanic sml male moths not only had higher activity but also had longer survival time (Fig 2). Among the candidate DEGs related to melanism, we found that a catalase gene (BGIBMGA011430) was significantly up-regulated in sml male, and an insulin-like peptide gene (Bombyxin E-1) was significantly down-regulated in sml male (Table 2 and S6 Table). It was reported that over-expression of catalase or mimetics extended the mean life-span of organisms [62,63]. And previous studies also showed aging in Drosophila could be slowed by reducing the insulin/IGF-1 signaling pathway (IIS) levels induced by decreasing insulin production [64,65]. According to the above information, we hypothesize that the longer survival of the melanic sml male moths might be associated with the enhanced antioxidant defense system and the reduction of IIS. Our future research will focus on this long survival phenotype and the underlying molecular mechanisms in the sml silkworm strain.
Moreover, previous reports have indicated that melanic insects have higher resistance to solar radiation than typical-colored insects [2]. Hu et al. demonstrated that the melanin granules deposited in the cuticle may help silkworms to resist UVA-induced damage [18]. Besides, dopamine is not only a major melanin precursor, but also an important neurotransmitter, influencing olfactory responses, memory retrieval, learning ability, locomotor activity, etc. in insects [66–68]. Furthermore, phenoloxidases serve as key enzymes in melanism and also participate in humoral immunity [69]. Finally, many developmental genes as well as their protein products are known to affect multiple traits [3]. Based on reports from various insects, we presume that the melanism associate phenotypic-effects could result from two probable causes: 1) direct influence of melanin granules or their precursors, and 2) additional functions of melanic alleles other than those participating in pigmentation.
In conclusion, our results revealed the DEGs related to the melanic morph and the derived phenotypic-effects in sml male moth. Based on our findings, we presume that the unknown sml melanic alleles might function by up-regulating the expression of yellow and laccase2. Additionally, we discovered the laccase2 B-type isoform was likely involved in the silkworm cuticular melanism process at late pupal stage. This is the first report to detect 22 cuticular protein family genes likely associated with melanin transport and/or maintenance in silk moths. Besides, we also hypothesize that the longer survival of the melanic sml male moths might be associated with the enhanced antioxidant defense system and the reduction of insulin/IGF-1 signaling pathway (IIS). Finally, we predict two probable causes for the derivation of phenotypic-effects from insect melanism. However, because only one developmental stage (P8, the initial stage of melanin formation and deposition) was selected for the study, there might be a risk of missing some information of candidate genes. Further studies are needed to address how cuticular proteins function in the melanism process, and to understand the regulatory mechanisms that cause melanism-induced phenotypic-effects at a wider range of developmental stages. The findings of this study provide a better understanding of the molecular mechanisms of melanization in the sml strain. They also facilitate the understanding of the molecular basis underlying melanism and the derived phenotypic-effects.
Supporting Information
[Figure omitted. See PDF.]
S1 Fig. Phenotypes of sclerotic tissues in the sml mutant and wild type (Dazao).
(A) Phenotypes of the sml strain in P9 (the 9th day of pupae). (B) The pupal case color of sml male was not different from the typical individuals. (C) The dorsal plate color of sml male was not visibly different from the male Dazao. (D) The tentacles color in Dazao and sml were not visibly different.
https://doi.org/10.1371/journal.pone.0155061.s001
(TIF)
S2 Fig. The sml female moths show normal phenotypes.
(A) There are no significant differences in fecundity between sml and wild-type females (Student’s t-test, two-tailed, n = 5). (B) There are no significant differences between sml and wild-type females in the adult longevity (Student’s t-test, two-tailed, n = 58). P-values less than 0.05 were considered significant. ♀ indicates female.
https://doi.org/10.1371/journal.pone.0155061.s002
(TIF)
S3 Fig. Scatterplot of significantly enriched KEGG pathway.
Rich factor is the ratio of the DEGs number to the total number of annotated genes in a certain pathway. Color and size of the dots represent the range of the q value and the gene number, respectively.
https://doi.org/10.1371/journal.pone.0155061.s003
(TIF)
S4 Fig. qRT-PCR validation of DEGs.
Fold changes of each gene in the four samples were calculated by dividing the expression level in the sml male sample. (A) RPKM fold changes of sequencing results. (B) qRT-PCR fold changes. All data are mean ± S.D (n = 3). DZ represents Dazao; M represents male; F represents female.
https://doi.org/10.1371/journal.pone.0155061.s004
(TIF)
S1 Table. Summary of RNA-seq data mapped to the silkworm reference genome.
https://doi.org/10.1371/journal.pone.0155061.s005
(DOCX)
S2 Table. Gene structural annotations and sequences of the novel transcripts.
https://doi.org/10.1371/journal.pone.0155061.s006
(XLS)
S3 Table. Expression levels of all the DEGs among the four samples by comparing any two samples.
https://doi.org/10.1371/journal.pone.0155061.s007
(XLSX)
S4 Table. Gene Ontology (GO) classification of the DEGs obtained respectively from the three melanic-typical comparisons (sml_M vs sml_F, sml_M vs DZ_M and sml_M vs DZ_F).
https://doi.org/10.1371/journal.pone.0155061.s008
(XLSX)
S5 Table. KEGG classification of the DEGs obtained respectively from the three melanic-typical comparisons.
https://doi.org/10.1371/journal.pone.0155061.s009
(XLSX)
S6 Table. Functional annotation of the 59 candidate DEGs associated with the melanic morph.
https://doi.org/10.1371/journal.pone.0155061.s010
(XLSX)
S7 Table. Primer sequences used in this work.
https://doi.org/10.1371/journal.pone.0155061.s011
(DOCX)
Acknowledgments
We wish to thank Dr. Niannian Li and Dr. Xin Ding for integument dissection, Dr. Liang Qiao and Dr. Jiangbo Song for valuable advices.
Author Contributions
Conceived and designed the experiments: SH XT FD. Performed the experiments: SH KL YL JL WY MC. Analyzed the data: SH YL JL MH. Contributed reagents/materials/analysis tools: FD HH CL. Wrote the paper: SH XT FD.
Citation: He S, Tong X, Lu K, Lu Y, Luo J, Yang W, et al. (2016) Comparative Analysis of Transcriptomes among Bombyx mori Strains and Sexes Reveals the Genes Regulating Melanic Morph and the Related Phenotypes. PLoS ONE 11(5): e0155061. https://doi.org/10.1371/journal.pone.0155061
1. Kettlewell B, Majerus M, Clarke C, Grant B, Ford EB, Wells J, et al. The Evolution of Melanism: A study of recurring necessity; with special reference to industrial melanism in the Lepidoptera. Oxford: Clarendon Press. 1973.
2. Majerus M. Melanism: Evolution in action. Oxford University Press. 1998.
3. True JR. Insect melanism: the molecules matter. TRENDS in Ecology and Evolution. 2003; 18: 640–647.
4. Wittkopp PJ, Carroll SB, Kopp A. Evolution in black and white: genetic control of pigment patterns in Drosophila. Trends Genet. 2003; 19: 495–504. pmid:12957543
5. Ninomiya Y, Tanaka K, Hayakawa Y. Mechanisms of black and white stripe pattern formation in the cuticles of insect larvae. J Insect Physiol. 2006; 52: 638–645. pmid:16618489
6. Dai FY, Qiao L, Tong XL, Cao C, Chen P, Chen J, et al. Mutations of an arylalkylamine-N-acetyltransferase, Bm-iAANAT, are responsible for silkworm melanism mutant. J Biol Chem. 2010; 285: 19553–19560. pmid:20332088
7. van't Hof AE, Edmonds N, Dalikova M, Marec F, Saccheri IJ. Industrial melanism in British peppered moths has a singular and recent mutational origin. Science. 2011; 332: 958–960. pmid:21493823
8. Wright TR. The genetics of biogenic amine metabolism, sclerotization, and melanization in Drosophila melanogaster. Adv Genet. 1987; 24: 127–222. pmid:3124532
9. Qiao L, Li Y, Xiong G, Liu X, He S, Tong X, et al. Effects of altered catecholamine metabolism on pigmentation and physical properties of sclerotized regions in the silkworm melanism mutant. PLoS One. 2012; 7: e42968. pmid:22937004
10. Dai F, Qiao L, Cao C, Liu X, Tong X, He S, et al. Aspartate Decarboxylase is Required for a Normal Pupa Pigmentation Pattern in the Silkworm, Bombyx mori. Sci Rep. 2015; 5: 10885. pmid:26077025
11. van't Hof AE, Saccheri IJ. Industrial melanism in the peppered moth is not associated with genetic variation in canonical melanisation gene candidates. PLoS One. 2010; 5: e10889. pmid:20526362
12. Ito K, Katsuma S, Kuwazaki S, Jouraku A, Fujimoto T, Sahara K, et al. Mapping and recombination analysis of two moth colour mutations, Black moth and Wild wing spot, in the silkworm Bombyx mori. Heredity (Edinb). 2015; 116: 52–59.
13. Brakefield PM. Industrial melanism: Do we have the answers? Trends Ecol Evol. 1987; 2: 117–122. pmid:21227832
14. Gershenson SM. A melanic form of the oak silkworm Antheraea pernyi (Lepidoptera, Attacidae). Vestnik Zool. 1994; 6: 46–51.
15. Verhoog MDea. Melanic moths and the ability to encapsulate parasitoid eggs and larvae. Proc Neth Entomol Soc. 1996; 7: 127–133.
16. Verhoog MDea. The influence of genes for melanism in the activity of the flour moth,Ephestia kuehniella. AnimBehav. 1998; 56: 683–688.
17. Kea Wilson. Melanism and disease resistance in insects. Ecol Lett. 2001; 4: 637–649.
18. Hu YG, Shen YH, Zhang Z, Shi GQ. Melanin and urate act to prevent ultraviolet damage in the integument of the silkworm, Bombyx mori. Arch Insect Biochem Physiol. 2013; 83: 41–55. pmid:23575996
19. Lu C, Dai F, Xiang Z. Studies on the Mutation Strains of the Bombyx mori Gene Bank. Sci Agric Sin. 2003; 36: 968–975.
20. Yea Banno. A Guide to the Silkworm Mutants—Gene Name and Gene (Silkwork Genetics Division, Institute of Genetic Resources, Kyushu University, 2005).
21. Futahashi R, Sato J, Meng Y, Okamoto S, Daimon T, Yamamoto K, et al. yellow and ebony are the responsible genes for the larval color mutants of the silkworm Bombyx mori. Genetics. 2008; 180: 1995–2005. pmid:18854583
22. Meng Y, Katsuma S, Daimon T, Banno Y, Uchino K, Sezutsu H, et al. The silkworm mutant lemon (lemon lethal) is a potential insect model for human sepiapterin reductase deficiency. J Biol Chem. 2009; 284: 11698–11705. pmid:19246455
23. Ito K, Katsuma S, Yamamoto K, Kadono-Okuda K, Mita K, Shimada T. Yellow-e determines the color pattern of larval head and tail spots of the silkworm Bombyx mori. J Biol Chem. 2010; 285: 5624–5629. pmid:19996320
24. Liu C, Yamamoto K, Cheng TC, Kadono-Okuda K, Narukawa J, Liu SP, et al. Repression of tyrosine hydroxylase is responsible for the sex-linked chocolate mutation of the silkworm, Bombyx mori. Proc Natl Acad Sci U S A. 2010; 107: 12980–12985. pmid:20615980
25. Fujii T, Abe H, Kawamoto M, Katsuma S, Banno Y, Shimada T. Albino (al) is a tetrahydrobiopterin (BH4)-deficient mutant of the silkworm Bombyx mori. Insect Biochem Mol Biol. 2013; 43: 594–600. pmid:23567588
26. Dai F, Hu H, He S, Lu C. Linkage Analysis of Bombyx mori Sex-controlled Melanism (sml) Mutation Using Morphological and Molecular Markers. Science of Sericulture. 2013; 39(3): 0442–0447.
27. Fang SM, Hu BL, Zhou QZ, Yu QY, Zhang Z. Comparative analysis of the silk gland transcriptomes between the domestic and wild silkworms. BMC Genomics. 2015; 16: 60. pmid:25887670
28. Wang X, Xiong M, Lei C, Zhu F. The developmental transcriptome of the synanthropic fly Chrysomya megacephala and insights into olfactory proteins. BMC Genomics. 2015; 16: 20. pmid:25612629
29. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013; 14: R36. pmid:23618408
30. Roberts A, Pimentel H, Trapnell C, Pachter L. Identification of novel transcripts in annotated genomes using RNA-Seq. Bioinformatics. 2011; 27: 2325–2329. pmid:21697122
31. Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015; 31: 166–169. pmid:25260700
32. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010; 26: 139–140. pmid:19910308
33. Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010; 11: R106. pmid:20979621
34. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B. 1995; 57: 289–300.
35. Young M, Wakefield M, Smyth G, Oshlack A. Method Gene ontology analysis for RNA-seq: accounting for selection bias. 2010.
36. Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011; 39(Web Server issue): W316–322. pmid:21715386
37. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2 (-Delta Delta C(T)) Method. Methods. 2001; 25(4): 402–408. pmid:11846609
38. Hall JC. Genetics of circadian rhythms. Annu Rev Genet. 1990; 24: 659–697. pmid:2088180
39. Kyriacou CP, Burnett B, Connolly K. The behavioural basis of overdominance in competitive mating success at the ebony locus of Drosophila melanogaster. Animal Behaviour. 1978; 26: 1195–1206.
40. Chen C, Song J, Chen M, Li Z, Tong X, Hu H, et al. Rhodiola rosea extends lifespan and improves stress tolerance in silkworm, Bombyx mori. Biogerontology. 2015.
41. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000; 28: 27–30. pmid:10592173
42. Futahashi R, Banno Y, Fujiwara H. Caterpillar color patterns are determined by a two-phase melanin gene prepatterning process: new evidence from tan and laccase2. Evol Dev. 2010; 12: 157–167. pmid:20433456
43. Walter MF, Black BC, Afshar G, Kermabon AY, Wright TR, Biessmann H. Temporal and spatial expression of the yellow gene in correlation with cuticle formation and dopa decarboxylase activity in Drosophila development. Dev Biol. 1991; 147: 32–45. pmid:1879614
44. Mitchell HK, Weber-Tracy UM, Schaar G. Aspects of Cuticle Formation in Drosophila melanogaster. J Exp Zool. 1971; 176: 429–437.
45. Futahashi R, Tanaka K, Matsuura Y, Tanahashi M, Kikuchi Y, Fukatsu T. Laccase2 is required for cuticular pigmentation in stinkbugs. Insect Biochem Mol Biol. 2011; 41: 191–196. pmid:21167282
46. Arakane Y, Muthukrishnan S, Beeman RW, Kanost MR, Kramer KJ. Laccase 2 is the phenoloxidase gene required for beetle cuticle tanning. Proc Natl Acad Sci U S A. 2005; 102: 11337–11342. pmid:16076951
47. He N, Botelho JM, McNall RJ, Belozerov V, Dunn WA, Mize T, et al. Proteomic analysis of cast cuticles from Anopheles gambiae by tandem mass spectrometry. Insect Biochem Mol Biol. 2007; 37: 135–146. pmid:17244542
48. Dittmer NT, Gorman MJ, Kanost MR. Characterization of endogenous and recombinant forms of laccase-2, a multicopper oxidase from the tobacco hornworm, Manduca sexta. Insect Biochem Mol Biol. 2009; 39: 596–606. pmid:19576986
49. Yatsu J, Asano T. Cuticle laccase of the silkworm, Bombyx mori: purification, gene identification and presence of its inactive precursor in the cuticle. Insect Biochem Mol Biol. 2009; 39: 254–262. pmid:19168135
50. Asano T, Taoka M, Yamauchi Y, Craig Everroad R, Seto Y, Isobe T, et al. Re-examination of a alpha-chymotrypsin-solubilized laccase in the pupal cuticle of the silkworm, Bombyx mori: Insights into the regulation system for laccase activation during the ecdysis process. Insect Biochem Mol Biol. 2014; 55C: 61–69. pmid:25460512
51. Gorman MJ, Sullivan LI, Nguyen TD, Dai H, Arakane Y, Dittmer NT, et al. Kinetic properties of alternatively spliced isoforms of laccase-2 from Tribolium castaneum and Anopheles gambiae. Insect Biochem Mol Biol. 2012; 42: 193–202. pmid:22198355
52. Gorman MJ, Dittmer NT, Marshall JL, Kanost MR. Characterization of the multicopper oxidase gene family in Anopheles gambiae. Insect Biochem Mol Biol. 2008; 38: 817–824. pmid:18675911
53. Moussian B. Recent advances in understanding mechanisms of insect cuticle differentiation. Insect Biochem Mol Biol. 2010; 40: 363–375. pmid:20347980
54. Qiao L, Xiong G, Wang RX, He SZ, Chen J, Tong XL, et al. Mutation of a cuticular protein, BmorCPR2, alters larval body shape and adaptability in silkworm, Bombyx mori. Genetics. 2014; 196: 1103–1115. pmid:24514903
55. G L.E., F H.S., S T.D., H D.J.. Correlations of ultrastructure and pigmentation suggest how genes control development of wing scales of Heliconius butterflies. J Res Lepid. 1988; 26:141–160.
56. Janssen JM, Monteiro A, Brakefield PM. Correlations between scale structure and pigmentation in butterfly wings. Evol Dev. 2001; 3: 415–423. pmid:11806637
57. Futahashi R, Shirataki H, Narita T, Mita K, Fujiwara H. Comprehensive microarray-based analysis for stage-specific larval camouflage pattern-associated genes in the swallowtail butterfly, Papilio xuthus. BMC Biol. 2012; 10: 46. pmid:22651552
58. Cornman RS, Togawa T, Dunn WA, He N, Emmons AC, Willis JH. Annotation and analysis of a large cuticular protein family with the R&R Consensus in Anopheles gambiae. BMC Genomics. 2008; 9: 22. pmid:18205929
59. Tang L, Liang J, Zhan Z, Xiang Z, He N. Identification of the chitin-binding proteins from the larval proteins of silkworm, Bombyx mori. Insect Biochem Mol Biol. 2010; 40: 228–234. pmid:20149871
60. Andersen SO. Are structural proteins in insect cuticles dominated by intrinsically disordered regions? Insect Biochem Mol Biol. 2011; 41: 620–627. pmid:21477652
61. Suzuki Y, Matsuoka T, Iimura Y, Fujiwara H. Ecdysteroid-dependent expression of a novel cuticle protein gene BMCPG1 in the silkworm, Bombyx mori. Insect Biochem Mol Biol. 2002; 32: 599–607. pmid:12020834
62. Fabrizio P, Liou LL, Moy VN, Diaspro A, Valentine JS, Gralla EB, et al. SOD2 functions downstream of Sch9 to extend longevity in yeast. Genetics. 2003; 163: 35–46. pmid:12586694
63. Melov S, Ravenscroft J, Malik S, Gill MS, Walker DW, Clayton PE, et al. Extension of life-span with superoxide dismutase/catalase mimetics. Science. 2000; 289: 1567–1569. pmid:10968795
64. Wessells RJ, Fitzgerald E, Cypser JR, Tatar M, Bodmer R. Insulin regulation of heart function in aging fruit flies. Nat Genet. 2004; 36: 1275–1281. pmid:15565107
65. Broughton SJ, Piper MD, Ikeya T, Bass TM, Jacobson J, Driege Y, et al. Longer lifespan, altered metabolism, and stress resistance in Drosophila from ablation of cells making insulin-like ligands. Proc Natl Acad Sci U S A. 2005; 102: 3105–3110. pmid:15708981
66. Bicker G, Menzel R. Chemical codes for the control of behaviour in arthropods. Nature. 1989; 337: 33–39. pmid:2562906
67. Wicker-Thomas C, Hamann M. Interaction of dopamine, female pheromones, locomotion and sex behavior in Drosophila melanogaster. J Insect Physiol. 2008; 54: 1423–1431. pmid:18762189
68. Tempel BL, Livingstone MS, Quinn WG. Mutations in the dopa decarboxylase gene affect learning in Drosophila. Proc Natl Acad Sci U S A. 1984; 81: 3577–3581. pmid:6427773
69. Lavine MD, Strand MR. Insect hemocytes and their role in immunity. Insect Biochem Mol Biol. 2002; 32: 1295–1309. pmid:12225920
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2016 He et al. This is an open access article distributed under the terms of the Creative Commons Attribution License: http://creativecommons.org/licenses/by/4.0/ (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
As a source of insect polymorphism, melanism plays an important role in ecological adaption and usually endows advantageous phenotypic-effects on insects. However, due to the mechanistic diversity, there are knowledge gaps in the molecular mechanisms underlying melanism and the related phenotypes. In silk moths, a recessive melanic mutant (sex-controlled melanism, sml) strain exhibits extended adult longevity. We took a transcriptome approach to perform a comparative analysis between this sml strain and a wild-type strain (Dazao). Our analysis resulted in the identification of 59 unique differentially expressed genes in the melanic mutant. Two key genes (laccase2 and yellow) involved in melanin formation were significantly up-regulated in melanic individuals. The laccase2 B-type isoform (BGIBMGA006746) was found to likely participate in the silkworm cuticular melanism process at late pupal stage. Moreover, we discovered 22 cuticular protein encoding genes with the possible function in melanin transport and/or maintenance. Based on our findings, we presume that the longer survival of the melanic sml male moths might be associated with the enhanced antioxidant defense systems and a reduction in the insulin/IGF-1 signaling pathway (IIS). These findings will facilitate the understanding of the molecular basis underlying melanism and the derived phenotypic-effects.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer