1.Astaurov BL, Rovinskaya IS. Chromosome complex of Ussuri geographical race of Bombyx mandarina M. with special reference to the problem of the origin of the domesticated silkworm, Bombyx mori. Cytology.1959;1:327–332.
2.Shimada T, Kurimoto Y, Kobayashi M. Phylogenetic relationship of silkmoths inferred from sequence data of the arylphorin gene. Mol Phylogenet Evol.1995;4(3):223–234.
3.Sun W, Yu H, Shen Y, Banno Y, Xiang Z, Zhang Z. Phylogeny and evolutionary history of the silkworm. Sci China Life Sci.2012;55(6):483–496.
4.Goldsmith MR, Shimada T, Abe H. The genetics and genomics of the silkworm, Bombyx mori. Annu Rev Entomol.2005;50:71–100.
5.Xia Q, Li S, Feng Q. Advances in silkworm studies accelerated by the genome sequencing of Bombyx mori. Annu Rev Entomol.2014;59:513–536.
6.Inoue S, Tanaka K, Arisaka F, Kimura S, Ohtomo K, Mizuno S. Silk fibroin of Bombyx mori is secreted, assembling a high molecular mass elementary unit consisting of H-chain, L-chain, and P25, with a 6:6:1 molar ratio. J Biol Chem.2000;275(51):40517–40528.
7.Gamo T. Genetic variants of the Bombyx mori silkworn encoding sericin proteins of different lengths. Biochem Genet.1982;20(1–2):165–177.
8.Chevillard M, Couble P, Prudhomme JC. Complete nucleotide sequence of the gene encoding the Bombyx mori silk protein P25 and predicted amino acid sequence of the protein. Nucleic Acids Res.1986;14(15):6341–6342.
9.Michaille JJ, Couble P, Prudhomme JC, Garel A. A single gene produces multiple sericin messenger RNAs in the silk gland of Bombyx mori. Biochimie.1986;68(10–11):1165–1173.
10.Yamaguchi K, Kikuchi Y, Takagi T, Kikuchi A, Oyama F, Shimura K, Mizuno S. Primary structure of the silk fibroin light chain determined by cDNA sequencing and peptide analysis. J Mol Biol.1989;210(1):127–139.
11.Zhou CZ, Confalonieri F, Medina N, Zivanovic Y, Esnault C, Yang T, Jacquet M, Janin J, Duguet M, Perasso R, et al. Fine organization of Bombyx mori fibroin heavy chain gene. Nucleic Acids Res.2000;28(12):2413–2419.
12.Li B, Wang X, Hou C, Xu A, Li M. Genetic analysis of quantitative trait loci for cocoon and silk production quantity in Bombyx mori (Lepidoptera: Bombycidae). Eur J Entomol.2013;110(2):205–213.
13.Zhan S, Huang J, Guo Q, Zhao Y, Li W, Miao X, Goldsmith MR, Li M, Huang Y. An integrated genetic linkage map for silkworms with three parental combinations and its application to the mapping of single genes and QTL. BMC Genomics.2009;10:389.
14.Li C, Zuo W, Tong X, Hu H, Qiao L, Song J, Xiong G, Gao R, Dai F, Lu C. A composite method for mapping quantitative trait loci without interference of female achiasmatic and gender effects in silkworm, Bombyx mori. Anim Genet.2015;46(4):426–432.
15.Li C, Tong X, Zuo W, Luan Y, Gao R, Han M, Xiong G, Gai T, Hu H, Dai F, et al. QTL analysis of cocoon shell weight identifies BmRPL18 associated with silk protein synthesis in silkworm by pooling sequencing. Scientific Reports.2017;7(1):17985.
16.Li J, Qin S, Yu H, Zhang J, Liu N, Yu Y, Hou C, Li M. Comparative Transcriptome Analysis Reveals Different Silk Yields of Two Silkworm Strains. PLoS One.2016;11(5):e0155329.
17.Luan Y, Zuo W, Li C, Gao R, Zhang H, Tong X, Han M, Hu H, Lu C, Dai F. Identification of Genes that Control Silk Yield by RNA Sequencing Analysis of Silkworm (Bombyx mori) Strains of Variable Silk Yield. Int J of Mol Sci.2018;19(12):3718.
18.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.
19.Li JY, Cai F, Ye XG, Liang JS, Li JK, Wu MY, Zhao D, Jiang ZD, You ZY, Zhong BX. Comparative Proteomic Analysis of Posterior Silk Glands of Wild and Domesticated Silkworms Reveals Functional Evolution during Domestication. J of Proteome Res.2017;16(7):2495–2507.
20.Dai ZJ, Sun W, Zhang Z. Comparative analysis of iTRAQ-based proteomes for cocoons between the domestic silkworm (Bombyx mori) and wild silkworm (Bombyx mandarina).. J. Proteomics.2019;192:366–373.
21.Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nature Protocols.2012;7(3):562–578.
22.Ma L, Xu H, Zhu J, Ma S, Liu Y, Jiang RJ, Xia Q, Li S. Ras1CA overexpression in the posterior silk gland improves silk yield. Cell Res.2011;21(6):934–943.
23.Zhang CD, Li FF, Chen XY, Huang MH, Zhang J, Cui H, Pan MH, Lu C. DNA replication events during larval silk gland development in the silkworm, Bombyx mori. J Insect Physiol.2012;58(7):974–978.
24.Oh H, Irvine KD. In vivo analysis of Yorkie phosphorylation sites. Oncogene.2009;28(17):1916–1927.
25.Oh H, Irvine KD. In vivo regulation of Yorkie phosphorylation and localization. Development.2008;135(6):1081–1088.
26.Zhou QZ, Fang SM, Zhang Q, Yu QY, Zhang Z. Identification and comparison of long non-coding RNAs in the silk gland between domestic and wild silkworms. Insect Sci.2018;25(4):604–616.
27.Xiang H, Liu X, Li M, Zhu Y, Wang L, Cui Y, Liu L, Fang G, Qian H, Xu A, et al. The evolutionary road from wild moth to domestic silkworm. Nature Ecology & Evolution.2018(8):1268–1279.
28.Mackay TF. The genetic architecture of quantitative traits. Annu Rev Genet.2001;35:303–339.
29.Lynch M, Walsh B: Genetics and Analysis of Quantitative Traits. Oxford: Oxford University Press; 1998.
30.Georges M, Charlier C, Hayes B. Harnessing genomic information for livestock improvement. Nature Reviews Genetics.2019;20(3):135–156.
31.Zhang P, Liu S, Song HS, Zhang G, Jia Q, Li S. YorkieCA overexpression in the posterior silk gland improves silk yield in Bombyx mori. J Insect Physiol.2017;100:93–99.
32.Chen F, Porter D, Vollrath F. Structure and physical properties of silkworm cocoons. J R Soc Interface.2012;9(74):2299–2308.
33.Guan J, Zhu W, Liu B, Yang K, Vollrath F, Xu J. Comparing the microstructure and mechanical properties of Bombyx mori and Antheraea pernyi cocoon composites. Acta Biomater.2017;47:60–70.
34.Zhang Y-Y, Wu Z-Y, Zhou Y, Zhao F, Li J. An analysis on morphological structure of cocoon filaments from diffrerent layer of mulberry silkworm and eri-silkworm cocoons. Science of Sericulture.2013;39(6):1126–1130.
35.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 Biology.2013;14(4):R36.
36.Patro R, Mount SM, Kingsford C. Sailfish enables alignment-free isoform quantification from RNA-seq reads using lightweight algorithms. Nat Biotechnol.2014;32(5):462–464.
37.Anders S, Pyl PT, Huber W. HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics.2014;31:166–169.
38.Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics.2010;26(1):139–140.
39.Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biology.2010;11(3):R25.
40.Wu Y, Cheng T, Liu C, Liu D, Zhang Q, Long R, Zhao P, Xia Q. Systematic Identification and Characterization of Long Non-Coding RNAs in the Silkworm, Bombyx mori. PLoS One.2016;11(1):e0147147.
41.Benjamini Y, Hochberg Y. Controlling the False Discovery Rate - a Practical And Powerful Approach To Multiple Testing. J Roy Stat Soc B Met.1995;57(1):289–300.
42.Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics.2008;9:559.
43.Miller JA, Horvath S, Geschwind DH. Divergence of human and mouse brain transcriptome highlights Alzheimer disease pathways. Proc Natl Acad Sci U S A.2010;107(28):12698–12703.
44.Zhang B, Horvath S. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol.2005;4:Article17.
45.Ravasz E, Somera AL, Mongru DA, Oltvai ZN, Barabasi AL. Hierarchical organization of modularity in metabolic networks. Science.2002;297(5586):1551–1555.
46.Yang L, Li Y, Wei Z, Chang X. Coexpression network analysis identifies transcriptional modules associated with genomic alterations in neuroblastoma. Biochimica et Biophysica Acta.2018;1864(6):2341–2348.
47.Lu X, Li QT, Xiong Q, Li W, Bi YD, Lai YC, Liu XL, Man WQ, Zhang WK, Ma B, et al. The transcriptomic signature of developing soybean seeds reveals the genetic basis of seed trait adaptation during domestication. The Plant Journal.2016;86(6):530–544.
48.Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol.2011;29(7):644–652.
49.Finn RD, Attwood TK, Babbitt PC, Bateman A, Bork P, Bridge AJ, Chang HY, Dosztanyi Z, El-Gebali S, Fraser M, et al. InterPro in 2017-beyond protein family and domain annotations. Nucleic Acids Res.2017;45(D1):D190-D199.
50.Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biology.2010;11(2):R14.
51.Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, Kong L, Gao G, Li CY, Wei L. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res.2011;39:W316–322.
52.Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics.2010;26(5):589–595.
53.McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, et al. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res.2010;20(9):1297–1303.
54.Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Genome Project Data Processing S. The Sequence Alignment/Map format and SAMtools. Bioinformatics.2009;25(16):2078–2079.
55.Garrison E, Marth G. Haplotype-based variant detection from short-read sequencing. https://arxiv.org/abs/1207.3907 (2012). Accessed 15 Oct 2019.
56.Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, Handsaker RE, Lunter G, Marth GT, Sherry ST, et al. The variant call format and VCFtools. Bioinformatics.2011;27(15):2156–2158.
57.Qiu Q, Wang L, Wang K, Yang Y, Ma T, Wang Z, Zhang X, Ni Z, Hou F, Long R, et al. Yak whole-genome resequencing reveals domestication signatures and prehistoric population expansions. Nature Communications.2015;6:10283.
58.Wang MS, Zhang RW, Su LY, Li Y, Peng MS, Liu HQ, Zeng L, Irwin DM, Du JL, Yao YG, et al. Positive selection rather than relaxation of functional constraint drives the evolution of vision during chicken domestication. Cell Res.2016;26(5):556–573.