High-resolution silkworm pan-genome provides genetic insights into artificial selection and ecological adaptation

Safe utilization and remediation potential of the mulberry-silkworm system in heavy metal-contaminated lands: A review

Mulberry cultivation and silkworm rearing hold a prominent position in the agricultural industries of many Asian countries, contributing to economic growth, sustainable development, and cultural heritage preservation. Applying the soil-mulberry-silkworm system (SMSS) to heavy metal (HM)-contaminated areas is significant economically, environmentally, and socially. The ultimate goal of this paper is to review the main research progress of SMSS under HM stress, examining factors affecting its safe utilization and remediation potential for HM-contaminated soils. HM tolerance of mulberry and silkworms relates to their growth stages. Based on the standards for HM contaminants in various mulberry and silkworm products and the bioconcentration factor of HMs at different parts of SMSS, we calculated maximum safe Cd and Pb levels for SMSS application on contaminated lands. Several remediation practices demonstrated mulberry's ability to grow on barren lands, absorb various HMs, while silkworm excreta can adsorb HMs and improve soil fertility. Considering multiple factors influencing HM tolerance and accumulation, we propose a decision model to guide SMSS application in polluted areas. Finally, we discussed the potential of using molecular breeding techniques to screen or develop varieties better suited for HM-contaminated regions. However, actual pollution scenarios are often complex, requiring consideration of multiple factors. More large-scale applications are crucial to enhance the theoretical foundation for applying SMSS in HM pollution risk areas.

Environmental concentrations of microplastic-induced gut microbiota and metabolite disruption in silkworm, Bombyx mori

Microplastics (MPs) existing extensively in various ecosystems can be ingested by marine organisms and enter the food chain, resulting the health risks from the presence of MPs in aquatic and terrestrial ecosystems. In the present study, an ideal model for Lepidoptera, the silkworm, Bombyx mori, was exposed to environmental concentrations (0.125 μg, 0.25 μg or 0.5 μg/diet) of MPs for 5 days, and the global changes in gut microbes and metabolites were subsequently examined via 16S rDNA sequencing and GC‒MS-based metabolomics. The results showed that MPs exposure did not seriously threaten survival but may regulate signaling pathways involved in development and cocoon production. MPs exposure induced gut microbiota perturbation according to the indices of α-diversity and β-diversity, and the functional prediction of the altered microbiome and associated metabolites demonstrated the potential roles of the altered microbiome following MPs exposure in the metabolic and physiological states of silkworm. The metabolites markedly altered following MPs exposure may play vital biological roles in energy metabolism, lipid metabolism, xenobiotic detoxification and the immune system by directly or indirectly affecting the physiological state of silkworms. These findings contribute to assessing the health risks of MPs exposure in model insects and provide novel insight into the toxicity mechanism of MPs.

A pan-genome of 69 Arabidopsis thaliana accessions reveals a conserved genome structure throughout the global species range

Although originally primarily a system for functional biology, Arabidopsis thaliana has, owing to its broad geographical distribution and adaptation to diverse environments, developed into a powerful model in population genomics. Here we present chromosome-level genome assemblies of 69 accessions from a global species range. We found that genomic colinearity is very conserved, even among geographically and genetically distant accessions. Along chromosome arms, megabase-scale rearrangements are rare and typically present only in a single accession. This indicates that the karyotype is quasi-fixed and that rearrangements in chromosome arms are counter-selected. Centromeric regions display higher structural dynamics, and divergences in core centromeres account for most of the genome size variations. Pan-genome analyses uncovered 32,986 distinct gene families, 60% being present in all accessions and 40% appearing to be dispensable, including 18% private to a single accession, indicating unexplored genic diversity. These 69 new Arabidopsis thaliana genome assemblies will empower future genetic research.

NextDenovo: an efficient error correction and accurate assembly tool for noisy long reads

Long-read sequencing data, particularly those derived from the Oxford Nanopore sequencing platform, tend to exhibit high error rates. Here, we present NextDenovo, an efficient error correction and assembly tool for noisy long reads, which achieves a high level of accuracy in genome assembly. We apply NextDenovo to assemble 35 diverse human genomes from around the world using Nanopore long-read data. These genomes allow us to identify the landscape of segmental duplication and gene copy number variation in modern human populations. The use of NextDenovo should pave the way for population-scale long-read assembly using Nanopore long-read data.

Gut Microbiota Analysis in Silkworms (Bombyx mori) Provides Insights into Identifying Key Bacterials for Inclusion in Artificial Diet Formulations

The gut microbiome significantly influences the health and productivity of silkworms (Bombyx mori), the cornerstone of sericulture. With the increasing use of cost-effective artificial diets in sericulture, it is crucial to understand how these diets impact the silkworm gut microbiomes. Here we employed 16S rRNA sequencing to delineate the impact of three distinct dietary regimens on the silkworm gut microbiomes: exclusive mulberry leaf diet (SY), exclusive artificial feed diet (SL), and a sequential transition from artificial feed to mulberry leaves (ZS). Our results unveiled stark differences in microbial diversity across the groups, with the ZS group displaying an intermediary complexity. LefSe and random forest analyses identified Methylobacteriaceae, Microbacterium, and Rhodococcus as significantly enriched in the ZS group, suggesting their potential to facilitate silkworms' adaptation to dietary transitions. Functional profiling revealed differential pathway regulation, indicating a metabolic reconfiguration in response to dietary modulations. Notably, the enrichment of Lactobacillus and Weissella in both the SL and ZS groups highlights their potential as probiotics in artificial diets. Our findings provide insights into the diet adaptation mechanisms of silkworm gut microbiota, paving the way for harnessing the intestinal bacteria to enhance silkworm health and silk production through targeted microbial interventions in sericulture practices.

Disruption of Zfh3 abolishes mulberry-specific monophagy in silkworm larvae

Feeding behavior is critical for insect survival and fitness. Most researchers have explored the molecular basis of feeding behaviors by identifying and elucidating the function of olfactory receptors (ORs) and gustatory receptors (GRs). Other types of genes, such as transcription factors, have rarely been investigated, and little is known about their potential roles. The silkworm (Bombyx mori) is a well-studied monophagic insect which primarily feeds on mulberry leaves, but the genetic basis of its monophagy is still not understood. In this report, we focused on a transcription factor encoded by the Zfh3 gene, which is highly expressed in the silkworm central and peripheral nervous systems, including brain, antenna, and maxilla. To investigate its function, Zfh3 was abrogated using clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR associated protein 9 (Cas9) mutagenesis. Since Zfh3 knockout homozygotes are not viable, we studied feeding behavior in heterozygotes, and found that disruption of Zfh3 affects both gustation and olfaction. Mutant larvae lose preference for mulberry leaves, acquire the ability to consume an expanded range of diets, and exhibit improved adaptation to the M0 artificial diet, which contains no mulberry leaves. These results provide the first demonstration that a transcription factor modulates feeding behaviors in an insect.

The BTB-ZF gene Bm-mamo regulates pigmentation in silkworm caterpillars

The color pattern of insects is one of the most diverse adaptive evolutionary phenotypes. However, the molecular regulation of this color pattern is not fully understood. In this study, we found that the transcription factor Bm-mamo is responsible for black dilute (bd) allele mutations in the silkworm. Bm-mamo belongs to the BTB zinc finger family and is orthologous to mamo in Drosophila melanogaster. This gene has a conserved function in gamete production in Drosophila and silkworms and has evolved a pleiotropic function in the regulation of color patterns in caterpillars. Using RNAi and clustered regularly interspaced short palindromic repeats (CRISPR) technology, we showed that Bm-mamo is a repressor of dark melanin patterns in the larval epidermis. Using in vitro binding assays and gene expression profiling in wild-type and mutant larvae, we also showed that Bm-mamo likely regulates the expression of related pigment synthesis and cuticular protein genes in a coordinated manner to mediate its role in color pattern formation. This mechanism is consistent with the dual role of this transcription factor in regulating both the structure and shape of the cuticle and the pigments that are embedded within it. This study provides new insight into the regulation of color patterns as well as into the construction of more complex epidermal features in some insects.

A review on complete silk gene sequencing and de novo assembly of artificial silk

Silk, especially spider and insect silk, is a highly versatile biomaterial with potential applications in biomedicine, materials science, and biomimetic engineering. The primary structure of silk proteins is the basis for the mechanical properties of silk fibers. Biotechnologies such as single-molecule sequencing have facilitated an increasing number of reports on new silk genes and assembled silk proteins. Therefore, this review aims to provide a comprehensive overview of the recent advances in representative spider and insect silk proteins, focusing on identification methods, sequence characteristics, and de novo design and assembly. The review discusses three identification methods for silk genes: polymerase chain reaction (PCR)-based sequencing, PCR-free cloning and sequencing, and whole-genome sequencing. Moreover, it reveals the main spider and insect silk proteins and their sequences. Subsequent de novo assembly of artificial silk is covered and future research directions in the field of silk proteins, including new silk genes, customizable artificial silk, and the expansion of silk production and applications are discussed. This review provides a basis for the genetic aspects of silk production and the potential applications of artificial silk in material science and biomedical engineering.

Comparative analysis of Presence-Absence gene Variations in five hard tick species: impact and functional considerations

Tick species are vectors of harmful human and animal diseases, and their expansion is raising concerns under the global environmental changes' scenario. Ticks host and transmit bacteria, protozoa and viruses, making the understanding of host-pathogen molecular pathways critical to development of effective disease control strategies. Despite the considerable sizes and repeat contents of tick genomes, individual tick genomics is perhaps the most effective approach to reveal genotypic traits of interest. Presence-Absence gene Variations (PAVs) can contribute to individual differences within species, with dispensable genes carried by subsets of individuals possibly underpinning functional significance at individual or population-levels. We exploited 350 resequencing datasets of Dermacentor silvarum, Haemaphysalis longicornis, Ixodes persulcatus, Rhipicephalus microplus and Rhipicephalus sanguineus hard tick specimens to reveal the extension of PAV and the conservation of dispensable genes among individuals and, comparatively, between species. Overall, we traced 550-3,346 dispensable genes per species and were able to reconstruct 5.3-7 Mb of genomic regions not included in the respective reference genomes, as part of the tick pangenomes. Both dispensable genes and de novo predicted genes indicated that PAVs preferentially impacted mobile genetic elements in these tick species.

From phyllosphere to insect cuticles: silkworms gather antifungal bacteria from mulberry leaves to battle fungal parasite attacks

BACKGROUND

Bacterial transfers from plants to insect herbivore guts have been well investigated. However, bacterial exchanges between plant phyllospheres and insect cuticles remain unclear, as does their related biological function.

RESULTS

Here, we report that the cuticular bacterial loads of silkworm larvae quickly increased after molting and feeding on the white mulberry (Morus alba) leaves. The isolation and examination of silkworm cuticular bacteria identified one bacterium Mammaliicoccus sciuri that could completely inhibit the spore germination of fungal entomopathogens Metarhizium robertsii and Beauveria bassiana. Interestingly, Ma. sciuri was evident originally from mulberry leaves, which could produce a secreted chitinolytic lysozyme (termed Msp1) to damage fungal cell walls. In consistency, the deletion of Msp1 substantially impaired bacterial antifungal activity. Pretreating silkworm larvae with Ma. sciuri cells followed by fungal topical infections revealed that this bacterium could help defend silkworms against fungal infections. Unsurprisingly, the protective efficacy of ΔMsp1 was considerably reduced when compared with that of wild-type bacterium. Administration of bacterium-treated diets had no negative effect on silkworm development; instead, bacterial supplementation could protect the artificial diet from Aspergillus contamination.

CONCLUSIONS

The results of this study evidence that the cross-kingdom transfer of bacteria from plant phyllospheres to insect herbivore cuticles can help protect insects against fungal parasite attacks. Video Abstract.

High-throughput and genome-scale targeted mutagenesis using CRISPR in a nonmodel multicellular organism, Bombyx mori

Large-scale genetic mutant libraries are powerful approaches to interrogating genotype-phenotype correlations and identifying genes responsible for certain environmental stimuli, both of which are the central goal of life science study. We produced the first large-scale CRISPR-Cas9-induced library in a nonmodel multicellular organism, Bombyx mori We developed a piggyBac-delivered binary genome editing strategy, which can simultaneously meet the requirements of mixed microinjection, efficient multipurpose genetic operation, and preservation of growth-defect lines. We constructed a single-guide RNA (sgRNA) plasmid library containing 92,917 sgRNAs targeting promoters and exons of 14,645 protein-coding genes, established 1726 transgenic sgRNA lines following microinjection of 66,650 embryos, and generated 300 mutant lines with diverse phenotypic changes. Phenomic characterization of mutant lines identified a large set of genes responsible for visual phenotypic or economically valuable trait changes. Next, we performed pooled context-specific positive screens for tolerance to environmental pollutant cadmium exposure, and identified KWMTBOMO12902 as a strong candidate gene for breeding applications in sericulture industry. Collectively, our results provide a novel and versatile approach for functional B. mori genomics, as well as a powerful resource for identifying the potential of key candidate genes for improving various economic traits. This study also shows the effectiveness, practicality, and convenience of large-scale mutant libraries in other nonmodel organisms.

An efficient and safe strategy for germ cell-specific automatic excision of foreign DNA in F1 hybrid transgenic silkworms

The safety of transgenic technology is a major obstacle in the popularization and use of transgenic silkworms and their products. In sericulture, only the first filial generation (F1 ) hybrid eggs produced by cross-breeding Japanese and Chinese original strains are usually used for the large-scale breeding of silkworms, but this may result in uncontrolled transgene dispersal during the popularization and application of the F1 hybrid transgenic eggs. To address this issue, we developed a safe and efficient strategy using the GAL4/Upstream activating sequence (UAS) system, the FLP/flippase recognition target (FRT) system, and the gonad-specific expression gene promoters (RSHP1p and Nanosp) for the germ cell-specific automatic excision of foreign DNA in the F1 hybrid transgenic silkworms. We established 2 types of activator strains, R1p::GAL4-Gr and Nsp::GAL4-Gr, containing the testis-specific GAL4 gene expression cassettes driven by RSHP1p or Nanosp, respectively, and 1 type of effector strain, UAS::FLP-Rg, containing the UAS-linked FLP gene expression cassette. The FLP recombinase-mediated sperm-specific complete excision of FRT-flanked target DNA in the F1 double-transgenic silkworms resulting from the hybridization of R1p::GAL4-Gr and UAS::FLP-Rg was 100%, whereas the complete excision efficiency resulting from the hybridization of Nsp::GAL4-Gr and UAS::FLP-Rg ranged from 13.73% to 80.3%. Additionally, we identified a gene, sw11114, that is expressed in both testis and ovary of Bombyx mori, and can be used to establish novel gonad-specific expression systems in transgenic silkworms. This strategy has the potential to fundamentally solve the safety issue in the production of F1 transgenic silkworm eggs and provides an important reference for the safety of transgenic technology in other insect species.

Cyclin B3 plays pleiotropic roles in female reproductive organogenesis and early embryogenesis in the silkworm, Bombyx mori

BACKGROUND

The reproductive system plays a crucial role in insect survival, reproduction and species specificity. Understanding the molecular mechanisms underlying reproductive organogenesis contributes to improving the efficiency of sterile insect technique marked by an eco-friendly pest management strategy. Lepidoptera is one of the largest orders of insects, most of which are major pests in agriculture and forestry. Our study aimed to screen the genes responsible for reproductive organogenesis and unravel the mechanism underlying female reproductive organ defects.

RESULTS

Morphological investigation of female reproductive organs showed a defective connection between oviductus geminus and oviductus communis on the second day of pupa (P2) in Speckled mutant silkworm. RNA_Seq identified a total of 18 049 transcripts that were expressed in the P2 female internal reproductive organs without ovary in Spc/+ compared to +Spc /+Spc . Differential expression analysis identified 312 up-regulated genes and 221 down-regulated genes in Spc/+. KEGG analysis identified 44 significantly enriched pathways. The results of qRT-PCR performed on 33 genes significantly matched the outcomes of the RNA_Seq. Dysfunction of Cyclin B3 resulted in a defective connection of the oviductus communis with the ovariole, dysfunction of oogenesis, and a petite body. Moreover, homozygous recessive lethality of Cyclin B3/Cyclin B3 occurred during early embryogenesis.

CONCLUSION

Our results suggest that Cyclin B3 is a pleiotropic functional gene that regulates early embryogenesis, oogenesis, development, and female reproductive organogenesis. These results showed that Cyclin B3 has significant effects on lepidopteran mortality, growth, and reproductive physiology, which might be considered a novel and potentially eco-friendly target for lepidopteran pest management. © 2023 Society of Chemical Industry.

Effects of Habitual Dietary Change on the Gut Microbiota and Health of Silkworms

Diet plays a crucial role in shaping the gut microbiota and overall health of animals. Traditionally, silkworms are fed fresh mulberry leaves, and artificial diets do not support good health. The aim of this study was to explore the relationship between the dietary transition from artificial diets to mulberry leaves and the effects on the gut microbiota and physiological changes in silkworms as a model organism. With the transition from artificial diets to mulberry leaves, the diversity of the silkworm gut microbiota increased, and the proportion of Enterococcus and Weissella, the dominant gut bacterial species in silkworms reared on artificial diets, decreased, whereas the abundance of Achromobacter and Rhodococcus increased. Dietary transition at different times, including the third or fifth instar larval stages, resulted in significant differences in the growth and development, immune resistance, and silk production capacity of silkworms. These changes might have been associated with the rapid adaptation of the intestinal microbiota of silkworms to dietary transition. This study preliminarily established a dietary transition-gut microbial model in silkworms based on the conversion from artificial diets to mulberry leaves, thus providing an important reference for future studies on the mechanisms through which habitual dietary changes affect host physiology through the gut microbiome.

Silk Bioconjugates: From Chemistry and Concept to Application

Medical silks have captured global interest. While silk sutures have a long track record in humans, silk bioconjugates are still in preclinical development. This perspective examines key advances in silk bioconjugation, including the fabrication of silk-protein conjugates, bioconjugated silk particles, and bioconjugated substrates to enhance cell-material interactions in two and three dimensions. Many of these systems rely on chemical modification of the silk biopolymer, often using carbodiimide and reactive ester chemistries. However, recent progress in enzyme-mediated and click chemistries has expanded the molecular toolbox to enable biorthogonal, site-specific conjugation in a single step when combined with recombinant silk fibroin tagged with noncanonical amino acids. This perspective outlines key strategies available for chemical modification, compares the resulting silk conjugates to clinical benchmarks, and outlines open questions and areas that require more work. Overall, this assessment highlights a domain of new sunrise capabilities and development opportunities for silk bioconjugates that may ultimately offer new ways of delivering improved healthcare.

SilkMeta: a comprehensive platform for sharing and exploiting pan-genomic and multi-omic silkworm data

The silkworm Bombyx mori is a domesticated insect that serves as an animal model for research and agriculture. The silkworm super-pan-genome dataset, which we published last year, is a unique resource for the study of global genomic diversity and phenotype-genotype association. Here we present SilkMeta (http://silkmeta.org.cn), a comprehensive database covering the available silkworm pan-genome and multi-omics data. The database contains 1082 short-read genomes, 546 long-read assembled genomes, 1168 transcriptomes, 294 phenotype characterizations (phenome), tens of millions of variations (variome), 7253 long non-coding RNAs (lncRNAs), 18 717 full length transcripts and a set of population statistics. We have compiled publications on functional genomics research and genetic stock deciphering (mutant map). A range of bioinformatics tools is also provided for data visualization and retrieval. The large batch of omics data and tools were integrated in twelve functional modules that provide useful strategies and data for comparative and functional genomics research. The interactive bioinformatics platform SilkMeta will benefit not only the silkworm but also the insect biology communities.

The loci of insect phenotypic evolution

Insects are important elements of terrestrial ecosystems because they pollinate plants, destroy crops, transmit diseases to livestock and humans, and are important components of food chains. Here, I used Gephebase, a manually curated database of genetic variants associated with natural and domesticated trait variation, to explore current knowledge about the genes and the mutations known to contribute to natural phenotypic variation in insects. Analysis of over 600 mutations reveals that data are concentrated toward certain species and traits and that experimental approaches have changed over time. The distribution of coding and cis-regulatory changes varies with traits, experimental approaches, and identified gene loci. Recent studies highlight the important role of standing variation, repeated mutations in hotspot genes, recombination, inversions, and introgression.

Gene presence/absence variation in Mytilus galloprovincialis and its implications in gene expression and adaptation

Presence/absence variation (PAV) is a well-known phenomenon in prokaryotes that was described for the first time in bivalves in 2020 in Mytilus galloprovincialis. The objective of the present study was to further our understanding of the PAV phenomenon in mussel biology. The distribution of PAV was studied in a mussel chromosome-level genome assembly, revealing a widespread distribution but with hotspots of dispensability. Special attention was given to the effect of PAV in gene expression, since dispensable genes were found to be inherently subject to distortions due to their sparse distribution among individuals. Furthermore, the high expression and strong tissue specificity of some dispensable genes, such as myticins, strongly supported their biological relevance. The significant differences in the repertoire of dispensable genes associated with two geographically distinct populations suggest that PAV is involved in local adaptation. Overall, the PAV phenomenon would provide a key selective advantage at the population level.

PanKmer: k-mer-based and reference-free pangenome analysis

SUMMARY

Pangenomes are replacing single reference genomes as the definitive representation of DNA sequence within a species or clade. Pangenome analysis predominantly leverages graph-based methods that require computationally intensive multiple genome alignments, do not scale to highly complex eukaryotic genomes, limit their scope to identifying structural variants (SVs), or incur bias by relying on a reference genome. Here, we present PanKmer, a toolkit designed for reference-free analysis of pangenome datasets consisting of dozens to thousands of individual genomes. PanKmer decomposes a set of input genomes into a table of observed k-mers and their presence-absence values in each genome. These are stored in an efficient k-mer index data format that encodes SNPs, INDELs, and SVs. It also includes functions for downstream analysis of the k-mer index, such as calculating sequence similarity statistics between individuals at whole-genome or local scales. For example, k-mers can be "anchored" in any individual genome to quantify sequence variability or conservation at a specific locus. This facilitates workflows with various biological applications, e.g. identifying cases of hybridization between plant species. PanKmer provides researchers with a valuable and convenient means to explore the full scope of genetic variation in a population, without reference bias.

AVAILABILITY AND IMPLEMENTATION

PanKmer is implemented as a Python package with components written in Rust, released under a BSD license. The source code is available from the Python Package Index (PyPI) at https://pypi.org/project/pankmer/ as well as Gitlab at https://gitlab.com/salk-tm/pankmer. Full documentation is available at https://salk-tm.gitlab.io/pankmer/.

iNaturalist is an open science resource for ecological genomics by enabling rapid and tractable records of initial observations of sequenced biological samples

The rapidly growing body of publicly available sequencing data for rare species and/or wild-caught samples is accelerating the need for detailed records of the samples used to generate datasets. Many already published datasets are unlikely to ever be reused, not due to problems with the data themselves, but due to their questionable or unverifiable origins. In this paper, I present iNaturalist-a pre-existing citizen science platform that allows people to post photo observations of organisms in nature-as a tool that allows genomics researchers to rapidly publish observations of samples used to generate sequencing datasets. This practice aligns with the values of the open science movement, and I also discuss how iNaturalist, along with other online resources, can be used to create an open genomics pipeline that enables future replication studies and ensures the value of genomics datasets to future research.

Characterization of the primary structure of the major silk gene, h-fibroin, across caddisfly (Trichoptera) suborders

Larvae of caddisflies (Trichoptera) produce silk to build various underwater structures allowing them to exploit a wide range of aquatic environments. The silk adheres to various substrates underwater and has high tensile strength, extensibility, and toughness and is of interest as a model for biomimetic adhesives. As a step toward understanding how the properties of underwater silk evolved in Trichoptera, we used genomic data to identify full-length sequences and characterize the primary structure of the major silk protein, h-fibroin, across the order. The h-fibroins have conserved termini and basic motif structure with high variation in repeating modules and variation in the percentage of amino acids, mainly proline. This finding might be linked to differences in mechanical properties related to the different silk usage and sets a starting point for future studies to screen and correlate amino acid motifs and other sequence features with quantifiable silk properties.

Telomere-to-telomere assemblies of 142 strains characterize the genome structural landscape in Saccharomyces cerevisiae

Pangenomes provide access to an accurate representation of the genetic diversity of species, both in terms of sequence polymorphisms and structural variants (SVs). Here we generated the Saccharomyces cerevisiae Reference Assembly Panel (ScRAP) comprising reference-quality genomes for 142 strains representing the species' phylogenetic and ecological diversity. The ScRAP includes phased haplotype assemblies for several heterozygous diploid and polyploid isolates. We identified circa (ca.) 4,800 nonredundant SVs that provide a broad view of the genomic diversity, including the dynamics of telomere length and transposable elements. We uncovered frequent cases of complex aneuploidies where large chromosomes underwent large deletions and translocations. We found that SVs can impact gene expression near the breakpoints and substantially contribute to gene repertoire evolution. We also discovered that horizontally acquired regions insert at chromosome ends and can generate new telomeres. Overall, the ScRAP demonstrates the benefit of a pangenome in understanding genome evolution at population scale.

Disruption of a BTB-ZF transcription factor causes female sterility and melanization in the larval body of the silkworm, Bombyx mori

The dilute black (bd) of the silkworm Bombyx mori is a recessive mutant that produces a grayish-black color in the larval integument, instead of the characteristic white color found in wild-type larvae. In addition, eggs produced by bd females are sterile due to a deficiency in the micropylar apparatus. We identified candidate genes responsible for the bd phenotype using publicly available RNA-seq data. One of these candidate genes was homologous to the maternal gene required for meiosis (mamo) of Drosophila melanogaster, which encodes a broad-complex, tramtrack, and bric-à-brac-zinc finger (BTB-ZF) transcription factor essential for female fertility. In three independent bd strains, the expression of the B. mori mamo (Bmmamo) was downregulated in the larval integument. Using a CRISPR/Cas9-mediated knockout strategy, we found that Bmmamo knockout mutants exhibit a grayish-black color in the larval integument and female infertility. Moreover, larvae obtained from the complementation cross between bd/+ mutants and heterozygous knockouts for the Bmmamo also exhibited a grayish-black color, indicating that Bmmamo is responsible for the bd phenotype. Gene expression analysis using Bmmamo knockout mutants suggested that the BmMamo protein suppresses the expression of melanin synthesis genes. Previous comparative genome analysis revealed that the Bmmamo was selected during silkworm domestication, and we found that Bmmamo expression in the larval integument is higher in B. mori than in the wild silkworm B. mandarina, suggesting that the Bmmamo is involved in domestication-associated pigmentation changes of the silkworm.

A review of the pangenome: how it affects our understanding of genomic variation, selection and breeding in domestic animals?

As large-scale genomic studies have progressed, it has been revealed that a single reference genome pattern cannot represent genetic diversity at the species level. While domestic animals tend to have complex routes of origin and migration, suggesting a possible omission of some population-specific sequences in the current reference genome. Conversely, the pangenome is a collection of all DNA sequences of a species that contains sequences shared by all individuals (core genome) and is also able to display sequence information unique to each individual (variable genome). The progress of pangenome research in humans, plants and domestic animals has proved that the missing genetic components and the identification of large structural variants (SVs) can be explored through pangenomic studies. Many individual specific sequences have been shown to be related to biological adaptability, phenotype and important economic traits. The maturity of technologies and methods such as third-generation sequencing, Telomere-to-telomere genomes, graphic genomes, and reference-free assembly will further promote the development of pangenome. In the future, pangenome combined with long-read data and multi-omics will help to resolve large SVs and their relationship with the main economic traits of interest in domesticated animals, providing better insights into animal domestication, evolution and breeding. In this review, we mainly discuss how pangenome analysis reveals genetic variations in domestic animals (sheep, cattle, pigs, chickens) and their impacts on phenotypes and how this can contribute to the understanding of species diversity. Additionally, we also go through potential issues and the future perspectives of pangenome research in livestock and poultry.

Toxicity Tests for Chemical Pesticide Registration: Requirement Differences among the United States, the European Union, Japan, and China?

Pesticide registration is a scientific, legal, and administrative process that checks if a pesticide is safe and effective for its intended use before it can be used. The toxicity test is a key part of pesticide registration, which includes human health and ecological effect testing. Different countries adopt their own toxicity test criteria for pesticide registration guidelines. However, these differences, which may help accelerate the progress of pesticide registration and reduce the number of animals used, are yet to be explored and compared. Herein, we outlined the details and compared the differences between the toxicity tests in the United States, the European Union, Japan, and China. Some differences lie in the types and waiver policy, while others are in new approach methodologies (NAMs). On the basis of these differences, there is great potential for the optimization of NAMs during the toxicity tests. It is expected that this perspective can contribute to developing and adopting NAMs.

Highly accurate genome assembly of an improved high-yielding silkworm strain, Nichi01

The silkworm (Bombyx mori) is an important lepidopteran model insect and an industrial domestic animal traditionally used for silk production. Here, we report the genome assembly of an improved Japanese strain Nichi01, in which the cocoon yield is comparable to that of commercial silkworm strains. The integration of PacBio Sequel II long-read and ddRAD-seq-based high-density genetic linkage map achieved the highest quality genome assembly of silkworms to date; 22 of the 28 pseudomolecules contained telomeric repeats at both ends, and only four gaps were present in the assembly. A total of 452 Mbp of the assembly with an N50 of 16.614 Mbp covered 99.3% of the complete orthologs of the lepidopteran core genes. Although the genome sequence of Nichi01 and that of the previously reported low-yielding tropical strain p50T assured their accuracy in most regions, we corrected several regions, misassembled in p50T, in our assembly. A total of 18,397 proteins were predicted using over 95 Gb of mRNA-seq derived from 10 different organs, covering 96.9% of the complete orthologs of the lepidopteran core genes. The final assembly and annotation files are available in KAIKObase (https://kaikobase.dna.affrc.go.jp/index.html) along with a genome browser and BLAST searching service, which would facilitate further studies and the breeding of silkworms and other insects.

Recessive embryonic lethal mutations uncovered in heterozygous condition in silkworm semiconsomic strains

In this study, we found two embryonic lethal mutations, t04 lethal (l-t04) and m04 lethal (l-m04), in semiconsomic strains T04 and M04, respectively. In these semiconsomic strains, the entire diploid genome, except for one chromosome 4 of the wild silkworm Bombyx mandarina, is substituted with chromosomes of the domesticated silkworm B. mori, and l-t04 and l-m04 mutations are located on B. mandarina-derived chromosome 4. To clarify the cause of the lethalities and the genes responsible for these mutations, positional cloning and CRISPR/Cas9 mediated knockout screening were performed. Finally, genetic complementation tests identified the mutations responsible for the l-t04 and l-m04 as the Bombyx homolog of imaginal discs arrested (Bmida) and TATA box binding protein-associated factor 5 (BmTaf5), respectively. Lethal stages of each knockout mutant indicated the importance of these genes in B. mori late embryogenesis. The lethal mutations responsible for l-t04 and l-m04 were not found in parental strains or wild B. mandarina collected from 39 distinct locations in Japan, indicating that both mutations were independently introduced during or after the development of the semiconsomic strains. We conclude that the recessive embryonic lethality in the T04 and M04 strains is due to deleterious mutations produced in B. mandarina-derived chromosome 4.

Combined CRISPR toolkits reveal the domestication landscape and function of the ultra-long and highly repetitive silk genes

Highly repetitive sequences play a major structural and function role in the genome. In the present study, we developed Cas9-assisted cloning and SMRT sequencing of long repetitive sequences (CACS) to sequence and manipulate highly repetitive genes from eukaryotic genomes. CACS combined Cas9-mediated cleavage of a target segment from an intact genome, Gibson assembly cloning, and PacBio SMRT sequencing. Applying CACS, we directly cloned and sequenced the complete sequences of fibroin heavy chain (FibH) genes from 17 domesticated (Bombyx mori) and 7 wild (Bombyx mandarina) silkworms. Our analysis revealed the unique fine structure organization, genetic variations, and domestication dynamics of FibH. We also demonstrated that the length of the repetitive regions determined the mechanical properties of silk fiber, which was further confirmed by Cas9 editing of FibH. CACS is a simple, robust, and efficient approach, providing affordable accessibility to highly repetitive regions of a genome. STATEMENT OF SIGNIFICANCE: Silkworm silk is the earliest and most widely used animal fiber, and its excellent performance mainly depends on the fibroin heavy chain (FibH) protein. The FibH gene is the main breakthrough in understanding the formation mechanism and improvement of silk fiber. In the study, we developed a CACS method for characterizing the fine structure and domestication landscape of 24 silkworm FibH genes. We used CRISPR/Cas9 to edit the repetitive sequence of FibH genes, revealing the relationship between FibH genes and mechanical properties of silkworm silk. Our study is helpful in modifying silk genes to manipulate other valuable highly repetitive sequences, and provides insight for silkworm breeding.

FibH Gene Complete Sequences (FibHome) Revealed Silkworm Pedigree

The highly repetitive and variable fibroin heavy chain (FibH) gene can be used as a silkworm identification; however, only a few complete FibH sequences are known. In this study, we extracted and examined 264 FibH gene complete sequences (FibHome) from a high-resolution silkworm pan-genome. The average FibH lengths of the wild silkworm, local, and improved strains were 19,698 bp, 16,427 bp, and 15,795 bp, respectively. All FibH sequences had a conserved 5' and 3' terminal non-repetitive (5' and 3' TNR, 99.74% and 99.99% identity, respectively) sequence and a variable repetitive core (RC). The RCs differed greatly, but they all shared the same motif. During domestication or breeding, the FibH gene mutated with hexanucleotide (GGTGCT) as the core unit. Numerous variations existed that were not unique to wild and domesticated silkworms. However, the transcriptional factor binding sites, such as fibroin modulator-binding protein, were highly conserved and had 100% identity in the FibH gene's intron and upstream sequences. The local and improved strains with the same FibH gene were divided into four families using this gene as a marker. Family I contained a maximum of 62 strains with the optional FibH (Opti-FibH, 15,960 bp) gene. This study provides new insights into FibH variations and silkworm breeding.

Deciphering the Genetic Basis of Silkworm Cocoon Colors Provides New Insights into Biological Coloration and Phenotypic Diversification

The genetic basis of phenotypic variation is a long-standing concern of evolutionary biology. Coloration has proven to be a visual, easily quantifiable, and highly tractable system for genetic analysis and is an ever-evolving focus of biological research. Compared with the homogenized brown-yellow cocoons of wild silkworms, the cocoons of domestic silkworms are spectacularly diverse in color, such as white, green, and yellow-red; this provides an outstanding model for exploring the phenotypic diversification and biological coloration. Herein, the molecular mechanism underlying silkworm green cocoon formation was investigated, which was not fully understood. We demonstrated that five of the seven members of a sugar transporter gene cluster were specifically duplicated in the Bombycidae and evolved new spatial expression patterns predominantly expressed in silk glands, accompanying complementary temporal expression; they synergistically facilitate the uptake of flavonoids, thus determining the green cocoon. Subsequently, polymorphic cocoon coloring landscape involving multiple loci and the evolution of cocoon color from wild to domestic silkworms were analyzed based on the pan-genome sequencing data. It was found that cocoon coloration involved epistatic interaction between loci; all the identified cocoon color-related loci existed in wild silkworms; the genetic segregation, recombination, and variation of these loci shaped the multicolored cocoons of domestic silkworms. This study revealed a new mechanism for flavonoids-based biological coloration that highlights the crucial role of gene duplication followed by functional diversification in acquiring new genetic functions; furthermore, the results in this work provide insight into phenotypic innovation during domestication.

Comparative Proteomic Analysis Provides New Insights into the Molecular Basis of Thermal-Induced Parthenogenesis in Silkworm (Bombyx mori)

Artificial parthenogenetic induction via thermal stimuli in silkworm is an important technique that has been used in sericultural production. However, the molecular mechanism underlying it remains largely unknown. We have created a fully parthenogenetic line (PL) with more than 85% occurrence and 80% hatching rate via hot water treatment and genetic selection, while the parent amphigenetic line (AL) has less than 30% pigmentation rate and less than 1% hatching rate when undergoing the same treatment. Here, isobaric tags for relative and absolute quantitation (iTRAQ)-based analysis were used to investigate the key proteins and pathways associated with silkworm parthenogenesis. We uncovered the unique proteomic features of unfertilized eggs in PL. In total, 274 increased abundance proteins and 211 decreased abundance proteins were identified relative to AL before thermal induction. Function analysis displayed an increased level of translation and metabolism in PL. After thermal induction, 97 increased abundance proteins and 187 decreased abundance proteins were identified. An increase in stress response-related proteins and decrease in energy metabolism suggested that PL has a more effective response to buffer the thermal stress than AL. Cell cycle-related proteins, including histones, and spindle-related proteins were decreased in PL, indicating an important role of this decrease in the process of ameiotic parthenogenesis.

Systematic identification of smORFs in domestic silkworm (Bombyx mori)

The silkworm (Bombyx mori) is not only an excellent model species, but also an important agricultural economic insect. Taking it as the research object, its advantages of low maintenance cost and no biohazard risks are considered. Small open reading frames (smORFs) are an important class of genomic elements that can produce bioactive peptides. However, the smORFs in silkworm had been poorly identified and studied. To further study the smORFs in silkworm, systematic genome-wide identification is essential. Here, we identified and analyzed smORFs in the silkworm using comprehensive methods. Our results showed that at least 738 highly reliable smORFs were found in B. mori and that 34,401 possible smORFs were partially supported. We also identified some differentially expressed and tissue-specific-expressed smORFs, which may be closely related to the characteristics and functions of the tissues. This article provides a basis for subsequent research on smORFs in silkworm, and also hopes to provide a reference point for future research methods for smORFs in other species.

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.