Artificial selection on cis-element of Abl contributes cocoon yield increase in domestic silkworm
Creators: name: Xiang, Hui
Publisher: name: Zenodo
Year: 2022
Subjects: Abl
artifical selection
domestication
silk yield
silkworm

description: The silkworm (<em>Bombyx mori</em>) is an important silk-producing domestic insect. The quality and yield of the silk produced by <em>B. mori </em>exceeds that of its ancestor, the wild silkworm <em>B. mandarina</em>. However, to date little is known about the molecular mechanisms underlying domestication-related increases in silk yield. Here, we identified a gene associated with both domestication and silk-related quantitative trait locus (QTL): Abelson tyrosine protein kinase (<em>Abl</em>). Population genomic data for <em>B. mori</em> and <em>B. mandarina</em> identified obvious signatures of artificial selection in the genomic region bearing the <em>Abl</em><em> </em>gene. There were two fixed nucleotide substitutions (−244 and −1311) in the transcription factor binding motif of the upstream regulatory region of <em>B. mori Abl</em>. Compared to <em>Abl </em>in the wild silkworm <em>B. mandarina</em>, <em>B. mori Abl</em> exhibited significantly greater promoter activity and was upregulated in the silk gland from the last day of the 5<sup>th</sup> larval instar to pupal stage (P0). Compared to the wild type, CRISPR/Cas9-generated <em>Abl </em>loss-of-function mutants exhibited a higher sensitivity to diseases, shorter developmental duration, and reductions in economically important silk traits, including cocoon weight, pupal weight, and cocoon layer thickness. Comparison of silk-gland transcriptomes between the wild type and the mutant indicated that genes enriched in ribosome biosynthesis, splicing, RNA transport, and the carbon metabolism were significantly downregulated in the mutants. Weighted gene co-expression network analysis (WGCNA) confirmed that genes related to ribosome biosynthesis were pivotal, driving the significant differentiation in silk yield between <em>B. mori</em> and <em>B. mandarina</em>. Here, we demonstrated that artificial selection acts on the cis-elements of silk-trait QTL gene <em>Abl </em>in a novel case (the domestic silkworm), and that this selection pressure increased <em>Abl </em>promoter activity and expression level. By promoting protein translation and synthesis, artificial selection further enhanced the robustness of silkworm larvae and improved cocoon silk synthesis.
description: Funding provided by: National Natural Science Foundation of China<br>Crossref Funder Registry ID: http://dx.doi.org/10.13039/501100001809<br>Award Number: 32070411Funding provided by: National Natural Science Foundation of China<br>Crossref Funder Registry ID: http://dx.doi.org/10.13039/501100001809<br>Award Number: 32000383Funding provided by: Natural Science Foundation of Guangdong Province<br>Crossref Funder Registry ID: http://dx.doi.org/10.13039/501100003453<br>Award Number: 2019A1515011012Funding provided by: the Laboratory of Lingnan Modern Agriculture Project*<br>Crossref Funder Registry ID: <br>Award Number: NZ2021019
description: Three duplicate samples were set for RNA-seq of slik glands of the wild type and mutant at the 6<sup>th</sup> day of the 5<sup>th</sup> instar larval stage (5L6D), wandering stage (W) and prepupae (PP), respectively. Silk glands were dissected, stored in dry ice and then sent to Novogene Company (Tianjin, China) for RNA extraction and RNA-seq. Sequencing libraries were generated using NEBNext® UltraTM RNA Library Prep Kit for Illumina® (NEB, USA) following manufacturer's instructions and index codes were added to attribute sequences to each sample. Clustering of the index-coded samples was performed on a cBot Cluster Generation System using HiSeq 4000 PE Cluster Kit (Illumia) according to the manufacturer's instructions. After cluster generation, the library preparations were sequenced on an Illumina Hiseq 4000 platform and 150bp paired-end reads were generated. Sequence raw data were qualified, filtered, built, and mapped to the silkworm reference genome (http://www.silkdb.org/silkdb/) by HISAT2 . HTSeq v0.6.0 was used to count the read numbers mapped to each gene and further normalized by DEseq2 in the R package 1.16.1 (TNLIST, Beijing, China) .Principal component analysis (PCA) was performed using genes expressed in the slikglands (Normalized reads &gt;5) with an online software. The expression value of each gene was calculated and normalized using fragments per kilobasesof exon per million reads mapped (FPKM).