Nucleotide diversity and selection signature in the domesticated silkworm, Bombyx mori, and wild silkworm, Bombyx mandarina

Bridging Nature and Engineering: Protein-Derived Materials for Bio-Inspired Applications

The sophisticated, elegant protein-polymers designed by nature can serve as inspiration to redesign and biomanufacture protein-based materials using synthetic biology. Historically, petro-based polymeric materials have dominated industrial activities, consequently transforming our way of living. While this benefits humans, the fabrication and disposal of these materials causes environmental sustainability challenges. Fortunately, protein-based biopolymers can compete with and potentially surpass the performance of petro-based polymers because they can be biologically produced and degraded in an environmentally friendly fashion. This paper reviews four groups of protein-based polymers, including fibrous proteins (collagen, silk fibroin, fibrillin, and keratin), elastomeric proteins (elastin, resilin, and wheat glutenin), adhesive/matrix proteins (spongin and conchiolin), and cyanophycin. We discuss the connection between protein sequence, structure, function, and biomimetic applications. Protein engineering techniques, such as directed evolution and rational design, can be used to improve the functionality of natural protein-based materials. For example, the inclusion of specific protein domains, particularly those observed in structural proteins, such as silk and collagen, enables the creation of novel biomimetic materials with exceptional mechanical properties and adaptability. This review also discusses recent advancements in the production and application of new protein-based materials through the approach of synthetic biology combined biomimetics, providing insight for future research and development of cutting-edge bio-inspired products. Protein-based polymers that utilize nature's designs as a base, then modified by advancements at the intersection of biology and engineering, may provide mankind with more sustainable products.

Midgut serine proteinases participate in dietary adaptations of the castor (Eri) silkworm Samia ricini Anderson transferred from Ricinus communis to an ancestral host, Ailanthus excelsa Roxb

Dietary change influenced the life-history traits, nutritional utilization, and midgut serine proteinases in the larvae of the domesticated polyphagous S. ricini, transferred from R. communis (common name: castor; family Euphorbiaceae; the host plant implicated in its domestication) to A. excelsa (common name: Indian tree of heaven; family Simaroubaceae; an ancestral host of wild Samia species). Significantly higher values for fecundity and body weight were observed in larvae feeding on R. communis (Scr diet), and they took less time to reach pupation than insects feeding on A. excelsa (Scai diet). Nevertheless, the nutritional index for efficiency of conversion of digested matter (ECD) was similar for larvae feeding on the two plant species, suggesting the physiological adaptation of S. ricini (especially older instars) to an A. excelsa diet. In vitro protease assays and gelatinolytic zymograms using diagnostic substrates and protease inhibitors revealed significantly elevated levels (p ≤ 0.05) of digestive trypsins, which may be associated with the metabolic costs influencing slow growth in larvae feeding on A. excelsa. RT-PCR with semidegenerate serine proteinase gene-specific primers, and cloning and sequencing of 3' cDNA ends identified a large gene family comprising at least two groups of putative chymotrypsins (i.e., Sr I and Sr II) resembling invertebrate brachyurins/collagenases with wide substrate specificities, and five groups of putative trypsins (i.e., Sr III, Sr IV, Sr V, Sr VII, and Sr VIII). Quantitative RT-PCR indicated that transcripts belonging to the Sr I, Sr III, Sr IV, and Sr V groups, especially the Sr IV group (resembling achelase I from Lonomia achelous), were expressed differentially in the midguts of fourth instars reared on the two plant species. Sequence similarity indicated shared lineages with lepidopteran orthologs associated with expression in the gut, protein digestion, and phytophagy. The results obtained are discussed in the context of larval serine proteinases in dietary adaptations, domestication, and exploration of new host plant species for commercial rearing of S. ricini.

Comparative Silk Transcriptomics Illuminates Distinctive Impact of Artificial Selection in Silkworm Modern Breeding

Early domestication and the following improvement are two important processes in the cocoon silk evolution of silkworms. In contrast to early domestication, understanding of the improvement process is still fuzzy. By systematically comparing the larval silk gland transcriptomes of the wild, early domestic, and improved silkworms, we highlighted a novel landscape of transcriptome in the silk glands of improved ones. We first clarified that silk cocoon protein genes were up-regulated in modern breeding but not in early domestication. Furthermore, we found that differentially expressed genes (DEGs) between improved and early domestic silkworms (2711), as well as between improved and wild silkworms (2264), were obviously more than those between the early domestic and wild silkworms (158), with 1671 DEGs specific in the improved silkworm (IS-DEGs). Hierarchical clustering of all the DEGs consistently indicated that improved silkworms were significantly diverged from the early domestic and wild silkworms, suggesting that modern breeding might cause prompt and drastic dynamic changes of gene expression in the silk gland. We further paid attention to these 1671 IS-DEGs and were surprised to find that down-regulated genes were enriched in basic organonitrogen compound biosynthesis, RNA biosynthesis, and ribosome biogenesis processes, which are generally universally expressed, whereas those up-regulated genes were enriched in organonitrogen compound catabolic processes and functions involving in the dynamic regulation of protein post-translation of modification. We finally highlighted one candidate improvement gene among these up-regulated IS-DEGs, i.e., GDAP2, which may play roles in silk behavior and the overall robustness of the improved silkworm. The findings strongly suggest that modern breeding may facilitate effective control of the basic consumption of nitrogen and a stronger switch of nitrogen resources from other tissues to the silk glands, for an efficient supply for silk production, and implies the importance of brain behavior and robustness in silk yield improvement of modern breeding.

Phylogenetic relationships among Bombycinae (Lepidoptera, Bombycoidea, and Bombycidae) based on mitochondrial genomes

The subfamily Bombycinae Latreille, [1802] is an important silk-producing group, including well-known economical insects. Although there are many studies on the development of these economic insects, the relationships between genera/species of this subfamily are still unclear. Two data sets of mitochondrial genomes, 13 protein-coding genes (13PCGs) and 13PCGs-AA, were used to estimate phylogenetic relationships based on the maximum likelihood and Bayesian inference methods. The results strongly support the subfamily Bombycinae as a monophyletic group divided into two clades.

Population Evolution, Genetic Diversity and Structure of the Medicinal Legume, Glycyrrhiza uralensis and the Effects of Geographical Distribution on Leaves Nutrient Elements and Photosynthesis

Glycyrrhiza uralensis is a valuable medicinal legume, which occurs widely in arid and semi-arid regions. G. uralensis demand has risen steeply due to its high medical and commercial value. Interpret genome-wide information can stimulate the G. uralensis development as far as its increased bioactive compound levels, and plant yield are concerned. In this study, leaf nutrient concentration and photosynthetic chlorophyll index of G. uralensis were evaluated to determine the G. uralensis growth physiology in three habitats. We observed that G. uralensis nutrient levels and photosynthesis differed significantly in three regions (p < 0.05). Whole-genome re-sequencing of the sixty G. uralensis populations samples from different habitats was performed using an Illumina HiSeq sequencing platform to elucidate the distribution patterns, population evolution, and genetic diversity of G. uralensis. 150.06 Gb high-quality clean data was obtained after strict filtering. The 895237686 reads were mapped against the reference genome, with an average 89.7% mapping rate and 87.02% average sample coverage rate. A total of 6985987 SNPs were identified, and 117970 high-quality SNPs were obtained after filtering, which were subjected to subsequent analysis. Principal component analysis (PCA) based on interindividual SNPs and phylogenetic analysis based on interindividual SNPs showed that the G. uralensis samples could be categorized into central, southern, and eastern populations, which reflected strong genetic differentiation due to long periods of geographic isolation. In this study, a total of 131 candidate regions were screened, and 145 candidate genes (such as Glyur001802s00036258, Glyur003702s00044485, Glyur001802s00036257, Glyur007364s00047495, Glyur000028s00003476, and Glyur000398s00034457) were identified by selective clearance analysis based on Fst and θπ values. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed significant enrichment of 110 GO terms including carbohydrate metabolic process, carbohydrate biosynthetic process, carbohydrate derivative biosynthetic process, and glucose catabolic process (p < 0.05). Alpha-linolenic acid metabolism, biosynthesis of unsaturated fatty acids, and fatty acid degradation pathways were significantly enriched (p < 0.05). This study provides information on the genetic diversity, genetic structure, and population adaptability of the medicinal legumes, G. uralensis. The data obtained in this study provide valuable information for plant development and future optimization of breeding programs for functional genes.

Phylogeographic Relationships among Bombyx mandarina (Lepidoptera: Bombycidae) Populations and Their Relationships to B. mori Inferred from Mitochondrial Genomes

We report 37 mitochondrial genome (mitogenome) sequences of Bombyx mori strains (Lepidoptera: Bombycidae) and four of B. mandarina individuals, each preserved and collected, respectively, in South Korea. These mitogenome sequences combined with 45 public data showed a substantial genetic reduction in B. mori strains compared to the presumed ancestor B. mandarina, with the highest diversity detected in the Chinese origin B. mori. Chinese B. mandarina were divided into northern and southern groups, concordant to the Qinling-Huaihe line, and the northern group was placed as an immediate progenitor of monophyletic B. mori strains in phylogenetic analyses, as has previously been detected. However, one individual that was in close proximity to the south Qinling-Huaihe line was exceptional, belonging to the northern group. The enigmatic South Korean population of B. mandarina, which has often been regarded as a closer genetic group to Japan, was most similar to the northern Chinese group, evidencing substantial gene flow between the two regions. Although a substantial genetic divergence is present between B. mandarina in southern China and Japan, a highly supported sister relationship between the two regional populations may suggest the potential origin of Japanese B. mandarina from southern China instead of the Korean peninsula.

Release of Mediator Enzyme β-Hexosaminidase and Modulated Gene Expression Accompany Hemocyte Degranulation in Response to Parasitism in the Silkworm Bombyx mori

In insects infections trigger hemocyte-mediated immune reactions including degranulation by exocytosis; however, involvement of mediator enzymes in degranulation process is unknown in insects. We report here that in silkworm Bombyx mori, infection by endoparasitoid Exorista bombycis and microsporidian Nosema bombycis activated granulation in granulocytes and promoted degranulation of accumulated structured granules. During degranulation the mediator lysosomal enzyme β-hexosaminidase showed increased activity and expression of β-hexosaminidase gene was enhanced. The events were confirmed in vitro after incubation of uninfected hemocytes with E. bombycis larval tissue protein. On infection, cytotoxicity marker enzyme lactate dehydrogenase (LDH) was released from the hemocytes illustrating cell toxicity. Strong positive correlation (R² = 0.71) between LDH activity and β-hexosaminidase released after the infection showed parasitic-protein-induced hemocyte damage and accompanied release of the enzymes. Expression of β-hexosaminidase gene was enhanced in early stages after infection followed by down regulation. The expression showed positive correlation (R² = 0.705) with hexosaminidase activity pattern. B. mori hexosaminidase showed 98% amino acid similarity with that of B. mandarina showing origin from same ancestral gene; however, 45-60% varied from other lepidopterans showing diversity. The observation signifies the less known association of hexosaminidase in degranulation of hemocytes induced by parasitic infection in B. mori and its divergence in different species.

Silkworm and spider silk electrospinning: a review

Issues of fossil fuel and plastic pollution are shifting public demand toward biopolymer-based textiles. For instance, silk, which has been traditionally used during at least 5 milleniums in China, is re-emerging in research and industry with the development of high-tech spinning methods. Various arthropods, e.g. insects and arachnids, produce silky proteinic fiber of unique properties such as resistance, elasticity, stickiness and toughness, that show huge potential for biomaterial applications. Compared to synthetic analogs, silk presents advantages of low density, degradability and versatility. Electrospinning allows the creation of nonwoven mats whose pore size and structure show unprecedented characteristics at the nanometric scale, versus classical weaving methods or modern techniques such as melt blowing. Electrospinning has recently allowed to produce silk scaffolds, with applications in regenerative medicine, drug delivery, depollution and filtration. Here we review silk production by the spinning apparatus of the silkworm Bombyx mori and the spiders Aranea diadematus and Nephila Clavipes. We present the biotechnological procedures to get silk proteins, and the preparation of a spinning dope for electrospinning. We discuss silk's mechanical properties in mats obtained from pure polymer dope and multi-composites. This review highlights the similarity between two very different yarn spinning techniques: biological and electrospinning processes.

Alternative moth-eye nanostructures: antireflective properties and composition of dimpled corneal nanocoatings in silk-moth ancestors

Moth-eye nanostructures are a well-known example of biological antireflective surfaces formed by pseudoregular arrays of nipples and are often used as a template for biomimetic materials. Here, we provide morphological characterization of corneal nanostructures of moths from the Bombycidae family, including strains of domesticated Bombyx mori silk-moth, its wild ancestor Bombyx mandarina, and a more distantly related Apatelodes torrefacta. We find high diversification of the nanostructures and strong antireflective properties they provide. Curiously, the nano-dimple pattern of B. mandarina is found to reduce reflectance as efficiently as the nanopillars of A. torrefacta. Access to genome sequence of Bombyx further permitted us to pinpoint corneal proteins, likely contributing to formation of the antireflective nanocoatings. These findings open the door to bioengineering of nanostructures with novel properties, as well as invite industry to expand traditional moth-eye nanocoatings with the alternative ones described here.

The potential for adaptive maintenance of diversity in insect antimicrobial peptides

Genes involved in immune defence are among the fastest evolving in the genomes of many species. Interestingly, however, genes encoding antimicrobial peptides (AMPs) have shown little evidence for adaptive divergence in arthropods, despite the centrality of these peptides in direct killing of microbial pathogens. This observation, coupled with a failure to detect phenotypic consequence of genetic variation in AMPs, has led to the hypothesis that individual AMPs make minor contributions to overall immune defence and that AMPs instead act as a collective cocktail. Recent data, however, have suggested an alternative explanation for the apparent lack of adaptive divergence in AMP genes. Molecular evolutionary and phenotypic data have begun to suggest that variant AMP alleles may be maintained through balancing selection in invertebrates, a pattern similar to that observed in several vertebrate AMPs. Signatures of balancing selection include high rates of non-synonymous polymorphism, trans-species amino acid polymorphisms, and convergence of amino acid states across the phylogeny. In this review, we revisit published literature on insect AMP genes and analyse newly available population genomic datasets in Drosophila, finding enrichment for patterns consistent with adaptive maintenance of polymorphism.This article is part of the themed issue 'Evolutionary ecology of arthropod antimicrobial peptides'.

Temperature- and sex-related effects of serine protease alleles on larval development in the Glanville fritillary butterfly

The body reserves of adult Lepidoptera are accumulated during larval development. In the Glanville fritillary butterfly, larger body size increases female fecundity, but in males fast larval development and early eclosion, rather than large body size, increase mating success and hence fitness. Larval growth rate is highly heritable, but genetic variation associated with larval development is largely unknown. By comparing the Glanville fritillary population living in the Åland Islands in northern Europe with a population in Nantaizi in China, within the source of the post-glacial range expansion, we identified candidate genes with reduced variation in Åland, potentially affected by selection under cooler climatic conditions than in Nantaizi. We conducted an association study of larval growth traits by genotyping the extremes of phenotypic trait distributions for 23 SNPs in 10 genes. Three genes in clip-domain serine protease family were associated with larval growth rate, development time and pupal weight. Additive effects of two SNPs in the prophenoloxidase-activating proteinase-3 (ProPO3) gene, related to melanization, showed elevated growth rate in high temperature but reduced growth rate in moderate temperature. The allelic effects of the vitellin-degrading protease precursor gene on development time were opposite in the two sexes, one genotype being associated with long development time and heavy larvae in females but short development time in males. Sexually antagonistic selection is here evident in spite of sexual size dimorphism.

Signatures of selection in tilapia revealed by whole genome resequencing

Natural selection and selective breeding for genetic improvement have left detectable signatures within the genome of a species. Identification of selection signatures is important in evolutionary biology and for detecting genes that facilitate to accelerate genetic improvement. However, selection signatures, including artificial selection and natural selection, have only been identified at the whole genome level in several genetically improved fish species. Tilapia is one of the most important genetically improved fish species in the world. Using next-generation sequencing, we sequenced the genomes of 47 tilapia individuals. We identified a total of 1.43 million high-quality SNPs and found that the LD block sizes ranged from 10–100 kb in tilapia. We detected over a hundred putative selective sweep regions in each line of tilapia. Most selection signatures were located in non-coding regions of the tilapia genome. The Wnt signaling, gonadotropin-releasing hormone receptor and integrin signaling pathways were under positive selection in all improved tilapia lines. Our study provides a genome-wide map of genetic variation and selection footprints in tilapia, which could be important for genetic studies and accelerating genetic improvement of tilapia.

Demographic history and gene flow during silkworm domestication

BACKGROUND

Gene flow plays an important role in domestication history of domesticated species. However, little is known about the demographic history of domesticated silkworm involving gene flow with its wild relative.

RESULTS

In this study, four model-based evolutionary scenarios to describe the demographic history of B. mori were hypothesized. Using Approximate Bayesian Computation method and DNA sequence data from 29 nuclear loci, we found that the gene flow at bottleneck model is the most likely scenario for silkworm domestication. The starting time of silkworm domestication was estimated to be approximate 7,500 years ago; the time of domestication termination was 3,984 years ago. Using coalescent simulation analysis, we also found that bi-directional gene flow occurred during silkworm domestication.

CONCLUSIONS

Estimates of silkworm domestication time are nearly consistent with the archeological evidence and our previous results. Importantly, we found that the bi-directional gene flow might occur during silkworm domestication. Our findings add a dimension to highlight the important role of gene flow in domestication of crops and animals.

Analysis of the genomic sequence of Philosamia cynthia nucleopolyhedrin virus and comparison with Antheraea pernyi nucleopolyhedrin virus

Background

Two species of wild silkworms, the Chinese oak silkworm (Antheraea pernyi) and the castor silkworm Philosamia cynthia ricini, can acquire a serious disease caused by Nucleopolyhedrin Viruses (NPVs) (known as AnpeNPV and PhcyNPV, respectively). The two viruses have similar polyhedral morphologies and their viral fragments share high sequence similarity. However, the physical maps of the viral genomes and cross-infectivity of the viruses are different. The genome sequences of two AnpeNPV isolates have been published.

Results

We sequenced and analyzed the full-length genome of PhcyNPV to compare the gene contents of the two viruses. The genome of PhcyNPV is 125, 376 bp, with a G + C content of 53.65%, and encodes 138 open reading frames (ORFs) of at least 50 amino acids (aa) (GenBank accession number: JX404026). Between PhcyNPV and AnpeMNPV-L and -Z isolates, 126 ORFs are identical, including 30 baculovirus core genes. Nine ORFs were only found in PhcyNPV. Four genes, cath, v-chi, lef 10 and lef 11, were not found in PhcyNPV. However, most of the six genes required for infectivity via the oral route were found in PhcyNPV and in the two AnpeNPV isolates, with high sequence similarities. The pif-3 gene of PhcyNPV contained 59 aa extra amino acids at the N-terminus compared with AnpeNPV.

Conclusions

Most of the genes in PhcyNPV are similar to the two AnpeNPV isolates, including the direction of expression of the ORFs. Only a few genes were missing from PhcyNPV. These data suggest that PhcyNPV and AnpeNPV might be variants of each other, and that the differences in cross-infection might be caused by gene mutations.

Little gene flow between domestic silkmoth Bombyx mori and its wild relative Bombyx mandarina in Japan, and possible artificial selection on the CAD gene of B. mori

We analyzed PCR-amplified carbamoyl-phosphate synthetase 2, aspartate transcarbamylase, and dihydroorotase (CAD) gene fragments from 146 Bombyx mori native strains and found extremely low levels of DNA polymorphism. Two haplotypes were identified, one of which was predominant. CAD haplotype analysis of 42 samples of Japanese B. mandarina revealed four haplotypes. No common haplotype was shared between the two species and at least five base substitutions were detected. This result was suggestive of low levels of gene flow between the two species. The nucleotide diversity (π) scores of the two samples differed markedly: lower π values were estimated for B. mori native strains than Japanese B. mandarina. We further analyzed 12 Chinese B. mandarina derived from seven areas of China, including Taiwan. The results clearly indicated that the π score was ~80-fold greater in Chinese B. mandarina than in B. mori. The extremely low level of DNA polymorphism in B. mori compared to its wild relatives suggested that the CAD gene itself or its tightly linked regions are possible targets for silkworm domestication.

Management increases genetic diversity of honey bees via admixture

The process of domestication often brings about profound changes in levels of genetic variation in animals and plants. The honey bee, Apis mellifera, has been managed by humans for centuries for both honey and wax production and crop pollination. Human management and selective breeding are believed to have caused reductions in genetic diversity in honey bee populations, thereby contributing to the global declines threatening this ecologically and economically important insect. However, previous studies supporting this claim mostly relied on population genetic comparisons of European and African (or Africanized) honey bee races; such conclusions require reassessment given recent evidence demonstrating that the honey bee originated in Africa and colonized Europe via two independent expansions. We sampled honey bee workers from two managed populations in North America and Europe as well as several old-world progenitor populations in Africa, East and West Europe. Managed bees had highly introgressed genomes representing admixture between East and West European progenitor populations. We found that managed honey bees actually have higher levels of genetic diversity compared with their progenitors in East and West Europe, providing an unusual example whereby human management increases genetic diversity by promoting admixture. The relationship between genetic diversity and honey bee declines is tenuous given that managed bees have more genetic diversity than their progenitors and many viable domesticated animals.