Search

Brown planthopper (BPH) is the most destructive insect pest of rice. Drought is the most detrimental environmental stress. BPH infestation causes adaxial leaf-rolling and bulliform cells (BCs) shrinkage similar to drought. The BC-related abaxially curled leaf1 (ACL1) gene negatively regulates BPH resistance and drought tolerance, with decreased cuticular wax in the gain-of-function mutant ACL1-D. ACL1 shows an epidermis-specific expression. The TurboID system and multiple biochemical assays reveal that ACL1 interacts with the epidermal-characteristic rice outermost cell-specific (ROC) proteins. ROC4 and ROC5 positively regulate BPH resistance and drought tolerance through modulating cuticular wax and BCs, respectively. Overexpression of ROC4 and ROC5 both rescue ACL1-D mutant in various related phenotypes. ACL1 competes with ROC4/ROC5 in homo-dimer and hetero-dimer formation, and interacts with the repressive TOPLESS-related proteins. Altogether, we illustrate that ACL1-ROC4/5 complexes synergistically mediate drought tolerance and BPH resistance through regulating cuticular wax content and BC development in rice, a mechanism that might facilitate BPH-resistant breeding., (© 2024. The Author(s).)

Manipulating grain size demonstrates great potential for yield promotion in cereals since it is tightly associated with grain weight. Several pathways modulating grain size have been elaborated in rice, but possible crosstalk between the ingredients is rarely studied. OsmiR396 negatively regulates grain size through targeting OsGRF4 (GS2) and OsGRF8, and proves to be multi-functioning. Here we showed that expression of GS3 gene, a Gγ-protein encoding gene, that negatively regulates grain size, was greatly down-regulated in the young embryos of MIM396, GRF8OE and GS2OE plants, indicating possible regulation of GS3 gene by OsmiR396/GRF module. Meanwhile, multiple biochemical assays proved possible transcriptional regulation of OsGRF4 and OsGRF8 proteins on GS3 gene. Further genetic relation analysis revealed tight genetic association between not only OsmiR396 and GS3 gene, but also GS2 and GS3 gene. Moreover, we revealed possible regulation of GS2 on four other grain size-regulating G protein encoding genes. Thus, the OsmiR396 pathway and the G protein pathway cross talks to regulate grain size. Therefore, we established a bridge linking the miRNA-transcription factors pathway and the G-protein signaling pathway that regulates grain size in rice., (© 2024. The Author(s).)

Leaf morphology is a crucial agronomic characteristic of rice that influences crop yield directly. One primary cause of rice leaf rolling can be attributed to alterations in bulliform cells. Several HD-ZIP IV genes have been identified to be epidemical characterized and function in leaf rolling in rice. Still others need to be studied to fully understand the overall function of HD-ZIP IV family. Among the nine ROC genes encoding HD-ZIP IV family transcription factors in rice, ROC1 exhibits the highest expression in the leaves. Overexpression of ROC1 decreased the size of bulliform cells, and thus resulted in adaxially rolled leaves. To the contrary, knockout of ROC1 (ROC1KO) through Crispr-cas9 system enlarged bulliform cells, and thus led to abaxially rolled leaves. Moreover, ROC1KO plants were sensitive to drought. ROC1 could form homodimers on its own, and heterodimers with ROC5 and ROC8 respectively. Compared to ROC1KO plants, leaves of the ROC1 and ROC8 double knocked out plants (ROC1/8DKO) were more severely rolled abaxially due to enlarged bulliform cells, and ROC1/8DKO plants were more drought sensitive. However, overexpression of ROC8 could not restore the abaxial leaf phenotype of ROC1KO plants. Therefore, we proved that ROC1, a member of the HD-ZIP IV family, regulated leaf rolling and drought stress response through tight association with ROC5 and ROC8., (© 2024. The Author(s).)

Background: The brown planthopper (BPH) is a kind of piercing-sucking insect specific to rice, with the damage tops the list of pathogens and insects in recent years. microRNAs (miRNAs) are pivotal regulators of plant-environment interactions, while the mechanism underlying their function against insects is largely unknown., Results: Here, we confirmed that OsmiR319, an ancient and conserved miRNA, negatively regulated resistance to BPHs, with overexpression of OsmiR319 susceptible to BPH, while suppression of OsmiR319 resistant to BPH in comparison with wild type. Meanwhile, we identified several targets of OsmiR319 that may mediate BPH resistance. Among them, OsPCF5 was the most obviously induced by BPH feeding, and over expression of OsPCF5 was resistance to BPH. In addition, various biochemical assays verified that OsPCF5 interacted with several MYB proteins, such as OsMYB22, OsMYB30, and OsMYB30C.Genetically, we revealed that both OsMYB22 and OsMYB30C positively regulated BPH resistance. Genetic interaction analyses confirmed that OsMYB22 and OsMYB30C both function in the same genetic pathway with OsmiR319b to mediate BPH resistance., Conclusions: Altogether, we revealed that OsPCF5 regulates BPH resistance via association with several MYB proteins downstream of OsmiR319, these MYB proteins might function as regulators of BPH resistance through regulating the phenylpropane synthesis., (© 2024. The Author(s).)

Background: With the development of single-cell technology, many cell traits can be measured. Furthermore, the multi-omics profiling technology could jointly measure two or more traits in a single cell simultaneously. In order to process the various data accumulated rapidly, computational methods for multimodal data integration are needed., Results: Here, we present inClust+, a deep generative framework for the multi-omics. It's built on previous inClust that is specific for transcriptome data, and augmented with two mask modules designed for multimodal data processing: an input-mask module in front of the encoder and an output-mask module behind the decoder. InClust+ was first used to integrate scRNA-seq and MERFISH data from similar cell populations, and to impute MERFISH data based on scRNA-seq data. Then, inClust+ was shown to have the capability to integrate the multimodal data (e.g. tri-modal data with gene expression, chromatin accessibility and protein abundance) with batch effect. Finally, inClust+ was used to integrate an unlabeled monomodal scRNA-seq dataset and two labeled multimodal CITE-seq datasets, transfer labels from CITE-seq datasets to scRNA-seq dataset, and generate the missing modality of protein abundance in monomodal scRNA-seq data. In the above examples, the performance of inClust+ is better than or comparable to the most recent tools in the corresponding task., Conclusions: The inClust+ is a suitable framework for handling multimodal data. Meanwhile, the successful implementation of mask in inClust+ means that it can be applied to other deep learning methods with similar encoder-decoder architecture to broaden the application scope of these models., (© 2024. The Author(s).)

Since the discovery of the RNA interference (RNAi) mechanism, it has been widely used in the fields of gene function, biomedicine, and crop pest control. In the direction of agricultural application, this technology is highly expected, especially in the field of pest control, which is called "the third revolution in the history of pesticides". Currently, RNA biopesticides are developing rapidly all over the world. A genetically modified product (MON87411) has been approved for marketing, and a large number of agricultural companies are developing products based on direct spraying RNA biopesticides and submitting them for regulatory approval. The biggest problem that has emerged for spray RNA biopesticides is the technology for large-scale and low-cost production of dsRNA. At present, the bacterial fermentation production technology can realize large-scale dsRNA production with a yield of 4.23~182 mg/L bacterial solution. This article describes the experimental protocol for dsRNA production technology based on bacterial fermentation., (© 2024. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

Brown planthopper (BPH) is the most destructive insect pest to rice that causes tremendous yield loss each year in rice planting Asia and South-East Asia areas. Compared with traditional chemical-based treatment, utilization of plant endogenous resistance is a more effective and environmental-friendly way for BPH control. Accordingly, quite a few quantitative trait loci (QTLs) for BPH resistance were cloned using forward genetics. However, BPH is apt to change quickly into new biotypes to overcome plant resistance, therefore, new resistance resources and genes are continuously needed. miRNAs are important regulators in both plant development and physiological regulation including immunity, and might be used as effective supplements for BPH resistance QTLs. miR159 is an ancient and conserved miRNA. In this study, we found that each OsMIR159 gene in rice responded to BPH feeding very obviously, and genetic function assay proved them to negatively regulate BPH resistance, with STTM159 showing resistance to BPH, and over expression of OsmiR159d susceptible to BPH. One target genes of OsmiR159, OsGAMYBL2, positively regulated BPH resistance. Further biochemical studies revealed that OsGAMYBL2 could directly bind to the promoter of G-protein γ subunit encoding GS3 gene and repress its expression. And genetically, GS3 responded to BPH feeding promptly and negatively regulated BPH resistance, GS3 over expression plants were susceptible to BPH, while GS3 knock-out plants were resistant to BPH. Thus, we identified new function of OsmiR159-OsGAMYBL2 in mediating BPH response, and revealed a new OsmiR159-G protein pathway that mediates BPH resistance in rice., (© 2023. The Author(s).)

The brown planthopper (BPH) is the most destructive pest of rice. The MYB transcription factors are vital for rice immunity, but most are activators. Although MYB22 positively regulates rice resistance to BPH and has an EAR motif associated with active repression, it remains unclear whether it is a transcriptional repressor affecting rice-BPH interaction. Genetic analyses revealed that MYB22 regulates rice resistance to BPH via its EAR motif. Several biochemical experiments (e.g. transient transcription assay, Y2H, LCA, and BiFC) indicated that MYB22 is a transcriptional repressor that interacts with the corepressor TOPLESS via its EAR motif and recruits HDAC1 to form a tripartite complex. Flavonoid-3'-hydroxylase (F3'H) is a flavonoid biosynthesis pathway-related gene that negatively regulates rice resistance to BPH. Based on a bioinformatics analysis and the results of EMSA and transient transcription assays, MYB22 can bind directly to the F3'H promoter and repress gene expression along with TOPLESS and HDAC1. We revealed a transcriptional regulatory mechanism influencing the rice-BPH interaction that differs from previously reported mechanisms. Specifically, MYB22-TOPLESS-HDAC1 is a novel transcriptional repressor complex with components that synergistically and positively regulate rice resistance to BPH through the transcriptional repression of F3'H., (© 2023 The Authors New Phytologist © 2023 New Phytologist Foundation.)

Brown planthopper (BPH) and blast disease jointly or individually cause big yield losses every year. To identify genes and metabolites with potential contributions to the dual resistance against both biotic-stress factors, we carried out a transcriptome and metabolome analysis for susceptible and resistant rice varieties after BPH and rice blast infestations. Coexpression network analysis identified a modular pattern that had the highest correlation coefficients (0.81) after the BPH and rice blast (-0.81) treatments. In total, 134 phenylpropanoid biosynthesis pathway-related genes were detected in this group. We found that the flavanone 3-hydroxylase gene (OsF3H) had opposite expression trends in response to BPH and rice blast infestations whereas the OsF3'H had similar expression patterns. Genetics analysis confirmed that the OsF3H gene knockdown lines demonstrated the opposite resistance phenotypes against BPH and rice blast, whereas the OsF3'H knockout lines enhanced rice resistance against both pests. Consistently, our metabolomics analysis identified the metabolite eriodictyol, one putative essential product of these two genes, that was more highly accumulated in the resistant rice variety of RHT than in the susceptible variety MDJ. This study highlights a useful strategy for identifying more genes and metabolites that have potential synergistic effects on rice against to multiple biotic stresses., (© 2022 John Wiley & Sons Ltd.)

To elucidate the application value of insect endogenous protease and its inhibitor genes in pest control, we analyzed in detail the transcriptome sequence of the Asian corn borer, Ostrinia furnacalis . We obtained 12 protease genes and 11 protease inhibitor genes, and comprehensively analyzed of their spatiotemporal expression by qRT-PCR. In which, a previous unstudied serine protease inhibitor gene attracted our attention. It belongs to the canonical serine proteinase inhibitor family, a trypsin inhibitor-like cysteine-rich domain (TIL)-type protease inhibitor, but its TIL domain lacks two cysteine residues, and it was named as ACB-TIL . Its expression level is relatively very low in the absence of pathogen stimulation, and can be up-regulated expression induced by Gram-negative bacteria ( Escherichia coli ), virus (BmNPV), and dsRNA (dsEGFP), but cannot be induced by fungus spores ( Metarrhizium anisopliae ). Prokaryotic expressed ACB-TIL protein can significantly inhibit the melanization in vitro . Injecting this protein into insect body can inhibit the production of antimicrobial peptides of attacin, lebocin and gloverin. Inhibition of ACB-TIL by RNAi can cause the responses of other immune-, protease- and inhibitor-related genes. ACB-TIL is primarily involved in Asian corn borer humoral immunity in responses to Gram-negative bacteria and viruses. This gene can be a potential target for pest control since this will mainly affect insect immune response., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Guan, Hu, Li, An, Miao and Li.)

The latest progresses of experimental biology have generated a large number of data with different formats and lengths. Deep learning is an ideal tool to deal with complex datasets, but its inherent "black box" nature needs more interpretability. At the same time, traditional interpretable machine learning methods, such as linear regression or random forest, could only deal with numerical features instead of modular features often encountered in the biological field. Here, we present MultiCapsNet (https://github.com/wanglf19/MultiCapsNet), a new deep learning model built on CapsNet and scCapsNet, which possesses the merits such as easy data integration and high model interpretability. To demonstrate the ability of this model as an interpretable classifier to deal with modular inputs, we test MultiCapsNet on three datasets with different data type and application scenarios. Firstly, on the labeled variant call dataset, MultiCapsNet shows a similar classification performance with neural network model, and provides importance scores for data sources directly without an extra importance determination step required by the neural network model. The importance scores generated by these two models are highly correlated. Secondly, on single cell RNA sequence (scRNA-seq) dataset, MultiCapsNet integrates information about protein-protein interaction (PPI), and protein-DNA interaction (PDI). The classification accuracy of MultiCapsNet is comparable to the neural network and random forest model. Meanwhile, MultiCapsNet reveals how each transcription factor (TF) or PPI cluster node contributes to classification of cell type. Thirdly, we made a comparison between MultiCapsNet and SCENIC. The results show several cell type relevant TFs identified by both methods, further proving the validity and interpretability of the MultiCapsNet., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Wang, Miao, Nie, Zhang, Zhang and Cai.)

In the present study, we identified a novel, positive-sense single-stranded RNA virus in the Chinese black cutworm, Agrotis ipsilon . It has a genome length of 11,312 nucleotides, excluding the poly(A) tails, and contains five open reading frames. The ORF2 encodes the conserved domains of RNA helicase and RNA-dependent RNA polymerase, while ORF4 and 5 encode three viral proteins. Herein, the A. ipsilon virus was clustered with a Helicoverpa armigera Nora virus and was thus provisionally named " Agrotis ipsilon Nora virus" (AINV). AINV was successfully transmitted into a novel host, Spodoptera frugiperda , through injection, causing a stable infection. This found the possibility of horizontal AINV transmission among moths belonging to the same taxonomic family. Nonetheless, AINV infection was deleterious to S. frugiperda and mainly mediated by antiviral and amino acid metabolism-related pathways. Furthermore, the infection significantly increased the S. frugiperda larval period but significantly reduced its moth eclosion rate. It suggests that AINV is probably to be a parasitic virus of S. frugiperda ., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Li, Guan, Wu, Chen, Yuan, Miao and Li.)

RNA interference (RNAi) is a valuable and revolutionary technology that has been widely applied in medicine and agriculture. The application of RNAi in various industries requires large amounts of low-cost double-stranded RNA (dsRNA). Chemical synthesis can only produce short dsRNAs; long dsRNAs need to be synthesized biologically. Several microbial chassis cells, such as Escherichia coli , Saccharomyces cerevisiae , and Bacillus species, have been used for dsRNA synthesis. However, the titer, rate of production, and yield of dsRNA obtained by these microorganism-based strategies is still low. In this review, we summarize advances in microbial dsRNA production, and analyze the merits and faults of different microbial dsRNA production systems. This review provides a guide for dsRNA production system selection. Future development of efficient microbial dsRNA production systems is also discussed., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Guan, Chu, Han, Miao and Li.)

microRNAs (miRNAs) are promising targets for crop improvement of complex agricultural traits. Coordinated activity between/among different miRNAs may fine-tune specific developmental processes in diverse organisms. Grain size is a main factor determining rice (Oryza sativa L.) crop yield, but the network of miRNAs influencing this trait remains uncharacterized. Here we show that sequestering OsmiR396 through target mimicry (MIM396) can substantially increase grain size in several japonica and indica rice subspecies and in plants with excessive tillers and a high panicle density. Thus, OsmiR396 has a major role related to the regulation of rice grain size. The grain shape of Growth Regulating Factor8 (OsGRF8)-overexpressing transgenic plants was most similar to that of MIM396 plants, suggesting OsGRF8 is a major mediator of OsmiR396 in grain size regulation. A miRNA microarray analysis revealed changes to the expression of many miRNAs, including OsmiR408, in the MIM396 plants. Analyses of gene expression patterns and functions indicated OsmiR408 is an embryo-specific miRNA that positively regulates grain size. Silencing OsmiR408 expression (miR408KO) using CRISPR technology resulted in small grains. Moreover, we revealed the direct regulatory effects of OsGRF8 on OsMIR408 expression. A genetic analysis further showed that the large-grain phenotype of MIM396 plants could be complemented by miR408KO. Also, several hormone signaling pathways might be involved in the OsmiR396/GRF-meditated grain size regulation. Our findings suggest that genetic regulatory networks comprising various miRNAs, such as OsmiR396 and OsmiR408, may be crucial for controlling rice grain size. Furthermore, the OsmiR396/GRF module may be important for breeding new high-yielding rice varieties., (© The Author(s) 2021. Published by Oxford University Press on behalf of American Society of Plant Biologists.)

Single-cell genomics has advanced rapidly as trace-DNA amplification technologies evolved. However, current technologies are subject to a variety of pitfalls such as contamination, uneven genomic coverage, and amplification errors. Even for the "golden" strategy of single stem cell-derived clonal formation, high-fidelity amplification is applicable merely to single stem cells. It's still challenging to accurately define somatic mutations of a single cell in various cell types. Herein, we provided evidence, for the first time, to prove that induced pluripotent stem cells (iPS cells or iPSC), being a single somatic cell-derived clone, are recording almost identical (>90%) mutational profile of the initial cell progenitor. This finding demonstrates iPS technique, applicable to any cell type, can be utilized as a cell cloning strategy favorable for single-cell genomic amplification. This novel strategy is not limited by cell-type constraints or amplification artifacts, and thus enables our detailed investigation on the characteristics of somatic mutations in heterogeneous normal cells., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2020 The Author(s).)

Drosophila development is governed by distinct ecdysone steroid pulses that initiate spatially and temporally defined gene expression programs. The translation of these signals into tissue-specific responses is crucial for metamorphosis, but the mechanisms that confer specificity to systemic ecdysone pulses are far from understood. Here, we identify Bric-à-brac 2 (Bab2) as an ecdysone-responsive transcriptional repressor that controls temporal gene expression during larval to pupal transition. Bab2 is necessary to terminate Salivary gland secretion (Sgs) gene expression, while premature Bab2 expression blocks Sgs genes and causes precocious salivary gland histolysis. The timely expression of bab2 is controlled by the ecdysone-responsive transcription factor Broad, and manipulation of EcR/USP/Broad signaling induces inappropriate Bab2 expression and termination of Sgs gene expression. Bab2 directly binds to Sgs loci in vitro and represses all Sgs genes in vivo. Our work characterizes Bab2 as a temporal regulator of somatic gene expression in response to systemic ecdysone signaling., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.)

Ethylene (ET) and jasmonic acid (JA) play important roles in plant defenses against biotic stresses. Crosstalk between JA and ET has been well studied in mediating pathogen resistance, but its roles in piercing-sucking insect resistance are unclear. The brown planthopper (BPH; Nilaparvata lugens) is the most notorious piercing-sucking insect specific to rice (Oryza sativa) that severely affects yield. A genetic analysis revealed that OsEBF1 and OsEIL1, which are in the ET signaling pathway, positively and negatively regulated BPH resistance, respectively. Molecular and biochemical analyses revealed direct interactions between OsEBF1 and OsEIL1. OsEBF1, an E3 ligase, mediated the degradation of OsEIL1 through the ubiquitination pathway, indicating the negative regulation of the ET-signaling pathway in response to BPH infestation. An RNA sequencing analysis revealed that a JA biosynthetic pathway-related gene, OsLOX9, was downregulated significantly in the oseil1 mutant. Biochemical analyses, including yeast one-hybrid, dual luciferase, and electrophoretic mobility shift assay, confirmed the direct regulation of OsLOX9 by OsEIL1. This study revealed the synergistic and negative regulation of JA and ET pathways in response to piercing-sucking insect attack. The synergistic mechanism was realized by transcriptional regulation of OsEIL1 on OsLOX9. OsEIL1-OsLOX9 is a novel crosstalk site in these two phytohormone signaling pathways., (© 2019 The Authors. New Phytologist © 2019 New Phytologist Trust.)

Double-stranded RNA (dsRNA)-induced genes are usually related to RNA interference (RNAi) mechanisms and are involved in immune-related pathways. In a previous study, we found a lepidopteran-specific nuclease gene REase that was up-regulated by dsRNA and that affected RNAi efficiency in Asian corn borer ( Ostrinia furnacalis ). In this study, to verify the function of REase , the homologous gene HaREase in cotton bollworm ( Helicoverpa armigera ) was knocked out using CRISPR/Cas9 system. We found that the midgut epithelium structure was apparently not affected in the Δ HaREase mutant [Knock out (KO)]. Transcript sequencing results showed that most of the known insect immune-related genes were up-regulated in KO. When second instar larvae were fed artificial diet with Cry1Ac, a protoxin from Bacillus thuringiensis ( Bt ), in sublethal doses (2.5 or 4 μg/g), the growth rate of KO was repressed significantly. The dsRNA stability was also enhanced in midgut extraction of KO; however, RNAi efficiency was not obviously improved compared with the wild type (WT). The KO and WT were injected with dsEGFP (Enhanced green fluorescent protein) and subjected to transcriptome sequencing. The results showed that the expression levels of 14 nuclease genes were enhanced in KO after the dsRNA treatment. These findings revealed that HaREase expression level was not only related with dsRNA stability, but also with Bt resistance in cotton bollworm. When HaREase was knocked out, other immune- or nuclease-related genes were enhanced significantly. These results remind us that insect immune system is complex and pest control for cotton bollworm is an arduous task., (Copyright © 2019 Guan, Chen, Li, Hu, Miao, Wang and Yang.)

Multi-functional microRNAs (miRNAs) are emerging as key modulators of plant-pathogen interactions. Although the involvement of some miRNAs in plant-insect interactions has been revealed, the underlying mechanisms are still elusive. The brown planthopper (BPH) is the most notorious rice (Oryza sativa)-specific insect that causes severe yield losses each year and requires urgent biological control. To reveal the miRNAs involved in rice-BPH interactions, we performed miRNA sequencing and identified BPH-responsive OsmiR396. Sequestering OsmiR396 by overexpressing target mimicry (MIM396) in three genetic backgrounds indicated that OsmiR396 negatively regulated BPH resistance. Overexpression of one BPH-responsive target gene of OsmiR396, growth regulating factor 8 (OsGRF8), showed resistance to BPH. Furthermore, the flavonoid contents increased in both the OsmiR396-sequestered and the OsGRF8 overexpressing plants. By analysing 39 natural rice varieties, the elevated flavonoid contents were found to correlate with enhanced BPH resistance. Artificial applications of flavonoids to wild type (WT) plants also increased resistance to BPH. A BPH-responsive flavanone 3-hydroxylase (OsF3H) gene in the flavonoid biosynthetic pathway was proved to be directly regulated by OsGRF8. A genetic functional analysis of OsF3H revealed its positive role in mediating both the flavonoid contents and BPH resistance. And analysis of the genetic correlation between OsmiR396 and OsF3H showed that down-regulation of OsF3H complemented the BPH resistance characteristic and simultaneously decreased the flavonoid contents of the MIM396 plants. Thus, we revealed a new BPH resistance mechanism mediated by the OsmiR396-OsGRF8-OsF3H-flavonoid pathway. Our study suggests potential applications of miRNAs in BPH resistance breeding., (© 2019 The Authors. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.)

Background: Rice (Oryza sativa) panicle architecture is the major determinant of the ideal plant architecture that directly influence yield potential. Many genes influencing development of primary branches, secondary branches, spikelet and pedicel would also influence panicle architecture, which is thus a complex trait regulated by genes from various aspects. miR156, an extensively studied miRNA, has recently emerged as promising target for crop improvement because of its role in plant architecture regulation, such as the number of tillers, plant height and the panicle architecture. Increasing evidence suggests that miR156 might play an important role in panicle architecture regulation., Main Body: To study the detailed function of miR156 in rice panicle architecture regulation, we examined the genetic interaction or transcriptional regulation of miR156/OsSPL to other panicle regulating genes. Our results revealed that expression of many panicle related genes were influenced by miR156. Through biochemical analysis, we further proved that miR156 directly regulated the axillary meristem regulating gene, LAX1, at the transcription level. And the intimate relations between miR156 and LAX1, and miR156 and LAX2 were also uncovered by genetic analysis. On the other hand, a tight genetic linkage between miR156 and RCN2, the panicle branch promoting gene, was also detected, which suggested a buffering mechanism for the miR156 mediated panicle architecture regulation. Furthermore, genetic analysis also demonstrated that miR156 functioned in the same pathway with OsRA2 to regulate pedicel length., Short Conclusion: Altogether, miR156 integrates several genetic pathways mediated by genes such as LAX1, LAX2, RCN2 and OsRA2, and comprehensively regulates panicle development in rice. Based on these analysis, we concluded that miR156 acts as an important regulator for panicle architecture through influencing various aspects of panicle development.

Background: In plants, basic helix-loop-helix (bHLH) proteins form the largest transcription factor (TF) family. Among them, HLH proteins are a small group of atypical members that lack the basic domain, and form dimers with bHLH proteins. Although bHLH proteins have been proved to play important roles in plant development and physiology, the function of HLH proteins is rarely studied, not to mention in plant biotic resistance. Brown planthopper (BPH) is a kind of rice-specific insect that causes devastating yield losses each year., Results: In this study, we identified OsHLH61 gene that encodes HLH protein. OsHLH61 gene could be highly induced by BPH infestation. Furthermore, Methyl Jasmonic acid (Me-JA) and cis-12-oxo- phytodienoic acid (OPDA) induced expression of OsHLH61, while SA repressed it. We knocked down expression of OsHLH61 by RNA interference (RNAi), the transgenic plants were susceptible to BPH infestation. RNA-seq analysis revealed that some pathogen-related (PR) genes in the Salicylic acid (SA) signaling pathway that mediate plant immunity were obviously down-regulated in the OsHLH61 RNAi plants. Meanwhile, yeast two-hybrid assay and bimolecular luciferase complementation (BiLC) analysis identified bHLH096 to be an interacting factor of OsHLH61. Also, some PR genes were down-regulated in the OsbHLH96 over expressing lines. Expression of OsbHLH96 was inhibited. Besides, OsbHLH96 might interact with Jasmonate Zim-Domain3 (OsJAZ3)., Conclusion: Altogether, we identified an OsHLH61-OsbHLH96 complex that might mediate defense to BPH through regulating PR genes. And OsHLH61-OsbHLH96 might be important in mediating SA and JA signaling crosstalk.

Numerous non-inherited somatic mutations, distinct from those of germ-line origin, occur in somatic cells during DNA replication per cell-division. The somatic mutations, recording the unique genetic cell-lineage 'history' of each proliferating normal cell, are important but remain to be investigated because of their ultra-low frequency hidden in the genetic background of heterogeneous cells. Luckily, the recent development of single-cell genomics biotechnologies enables the screening and collection of the somatic mutations, especial single nucleotide variations (SNVs), occurring in normal cells. Here, we established DSMNC: a database of somatic mutations in normal cells (http://dsmnc.big.ac.cn/), which provides most comprehensive catalogue of somatic SNVs in single cells from various normal tissues. In the current version, the database collected ∼0.8 million SNVs accumulated in ∼600 single normal cells (579 human cells and 39 mouse cells). The database interface supports the user-friendly capability of browsing and searching the SNVs and their annotation information. DSMNC, which serves as a timely and valuable collection of somatic mutations in individual normal cells, has made it possible to analyze the burdens and signatures of somatic mutations in various types of heterogeneous normal cells. Therefore, DSMNC will significantly improve our understanding of the characteristics of somatic mutations in normal cells.

Exogenous dsRNA enters the insect body and can induce the RNAi effect only when it is cleaved into siRNA. However, what kinds of base composition are easier to cut and what kinds of siRNA will be produced in vivo is largely unknown. In this study, we found that dsRNA processing into siRNA has sequence preference and regularity in insects. We injected 0.04 mg/g dsRNA into Asian corn borers or cotton bollworms according to their body weight, and then the siRNAs produced in vivo were analyzed by RNA-Seq. We discovered that a large number of siRNAs were produced with GGU nucleotide residues at the 5'- and 3'-ends and produced a siRNA peak on the sequence. Once the GGU site is mutated, the number of siRNAs will decrease significantly and the siRNA peak will also lost. However, in the red flour beetle, a member of Coleoptera, dsRNA was cut at more diverse sites, such as AAG, GUG, and GUU; more importantly, these enzyme restriction sites have a high conservation base of A/U. Our discovery regarding dsRNA in vivo cleavage preference and regularity will help us understand the RNAi mechanism and its application.

Tiller number and plant height are two of the main features of plant architecture that directly influence rice yield. Auxin and miR156, an extensively studied small RNA (smRNA), are both broadly involved in plant development and physiology, suggesting a possible relationship between the two. In this study, we identified a rice T-DNA insertion cluster and dwarf (cd) mutant that has an increased tiller number and reduced plant height. The T-DNA insertion was in close proximity to the miR156f gene and was associated with its up-regulation. Plants overexpressing miR156f resembled the cd mutant. In contrast, plants overexpressing an miR156f target mimic (MIM156fOE) had a reduced tiller number and increased height. Genetic analysis showed that OsSPL7 is a target of miR156f that regulates plant architecture. Plants overexpressing OsSPL7 had a reduced tiller number, while OsSPL7 RNAi plants had an increased tiller number and a reduced height. We also found that OsSPL7 binds directly to the OsGH3.8 promoter to regulate its transcription. Overexpression of OsGH3.8 and OsGH3.8 RNAi partially complemented the MIM156fOE and cd mutant phenotypes, respectively. Our combined data show that the miR156f-OsSPL7-OsGH3.8 pathway regulates tiller number and plant height in rice, and this pathway may allow crosstalk between miR156 and auxin.

Key Message: OsEXPA10 gene coordinates the balance between rice development and biotic resistance. Expansins are proteins that can loosen the cell wall. Previous studies have indicated that expansin-encoding genes were involved in defense against abiotic stress, but little is known about the involvement of expansins in biotic stress. Brown planthopper (BPH) is one of the worst insect pests of rice in the Asia-Pacific planting area, and many efforts have been made to identify and clone BPH-resistance genes for use in breeding resistant cultivars. At the same time, rice blast caused by Magnaporthe grisea is one of the three major diseases that severely affect rice production worldwide. Here, we demonstrated that one rice expansin-encoding gene, OsEXPA10, functions in both rice growth and biotic resistance. Over expression of OsEXPA10 improved rice growth but also increased susceptibility to BPH infestation and blast attack, while knock-down OsEXPA10 gene expression resulted in reduced plant height and grain size, but also increased resistance to BPH and the blast pathogen. These results imply that OsEXPA10 mediates the balance between rice development and biotic resistance.

The panicle architecture of rice is an important characteristic that influences reproductive success and yield. It is largely determined by the number and length of the primary and secondary branches. The number of panicle branches is defined by the inflorescence meristem state between determinacy and indeterminacy; for example, the maize ramosa2 ( ra2 ) mutant has more branches in its tassel through loss of spikelet determinacy. Some genes and factors influencing the number of primary and secondary branches have been studied, but little is known about the molecular mechanism underlying pedicel development, which also influences panicle architecture. We report here that rice OsRAMOSA2 ( OsRA2 ) gene modifies panicle architecture through regulating pedicel length. Ectopic expression of OsRA2 resulted in a shortened pedicel while inhibition of OsRA2 through RNA interference produced elongated pedicel. In addition, OsRA2 influenced seed morphology. The OsRA2 protein localized to the nucleus and showed transcriptional activation in yeast; in accordance with its function in pedicel development, OsRA2 mRNA was enriched in the anlagen of axillary meristems, such as primary and secondary branch meristems and the spikelet meristems of young panicles. This indicates a conserved role of OsRA2 for shaping the initial steps of inflorescence architecture. Genetic analysis revealed that OsRA2 may control panicle architecture using the same pathway as that of the axillary meristem gene LAX1 ( LAX PANICLE1 ). Moreover, OsRA2 acted downstream of RCN2 in regulating pedicel and branch lengths, but upstream of RCN2 for control of the number of secondary branches, indicating that branch number and length development in the panicle were respectively regulated using parallel pathway. Functional conservation between OsRA2 and AtLOB , and the conservation and diversification of RA2 in maize and rice are also discussed.

Specification of floral organ identity is critical for the establishment of floral morphology and inflorescence architecture. Although multiple genes are involved in the regulation of floral organogenesis, our understanding of the underlying regulating network is still fragmentary. MADs-box genes are principle members in the ABCDE model that characterized floral organs. OsMADS1 specifies the determinacy of spikelet meristem and lemma/palea identity in rice. However, the pathway through which OsMADS1 regulates floral organs remains elusive; here, we identified the microRNA172 (miR172) family as possible regulators downstream of OsMADS1 . Genetic study revealed that overexpression of each miR172 gene resulted in elongated lemma/palea and indeterminacy of the floret, which resemble the phenotype of osmads1 mutant. On the contrary, overexpression of each target APETALA2 ( AP2 ) genes resulted in shortened palea/lemma. Expression level and specificity of miR172 was greatly influenced by OsMADS1 , as revealed by Northern blot analysis and In situ hybridization. Genetically, AP2-3 and AP2-2 over expression rescued the elongation and inconsistent development of the lemma/palea in OsMADS1RNAi transgenic plants. Our results suggested that in rice, OsMADS1 and miR172s/ AP2 s formed a regulatory network involved in floral organ development, particularly the elongation of the lemma and the palea.

Acetoacetate (AA) is a ketone body and acts as a fuel to supply energy for cellular activity of various tissues. Here, we uncovered a novel function of AA in promoting muscle cell proliferation. Notably, the functional role of AA in regulating muscle cell function is further evidenced by its capability to accelerate muscle regeneration in normal mice, and it ameliorates muscular dystrophy in mdx mice. Mechanistically, our data from multiparameter analyses consistently support the notion that AA plays a non-metabolic role in regulating muscle cell function. Finally, we show that AA exerts its function through activation of the MEK1-ERK1/2-cyclin D1 pathway, revealing a novel mechanism in which AA serves as a signaling metabolite in mediating muscle cell function. Our findings highlight the profound functions of a small metabolite as signaling molecule in mammalian cells., (© 2016 by The American Society for Biochemistry and Molecular Biology, Inc.)

Insect double-stranded (ds)RNA expression in transgenic crops can increase plant resistance to biotic stress; however, creating transgenic crops to defend against every insect pest is impractical. Arabidopsis Mob1A is required for organ growth and reproduction. When Arabidopsis roots were soaked in dsMob1A, the root lengths and numbers were significantly suppressed and plants could not bolt or flower. Twenty-four hours after rice roots were immersed in fluorescent-labelled dsEYFP (enhanced yellow fluorescent protein), fluorescence was observed in the rice sheath and stem and in planthoppers feeding on the rice. The expression levels of Ago and Dicer in rice and planthoppers were induced by dsEYFP. When rice roots were soaked in dsActin, their growth was also significantly suppressed. When planthoppers or Asian corn borers fed on rice or maize that had been irrigated with a solution containing the dsRNA of an insect target gene, the insect's mortality rate increased significantly. Our results demonstrate that dsRNAs can be absorbed by crop roots, trigger plant and insect RNAi and enhance piercing-sucking and stem-borer insect mortality rates. We also confirmed that dsRNA was stable under outdoor conditions. These results indicate that the root dsRNA soaking can be used as a bioinsecticide strategy during crop irrigation., (© 2015 John Wiley & Sons Ltd.)

The silkworm Bombyx mori is an oligophagous insect that feeds mainly on mulberry leaves. The olfactory system of silkworm is a good model to study olfaction in Lepidoptera. Here, we carried out shotgun proteomic analysis and MS sequencing of the silkmoth antennae. A total of 364 proteins were detected, 77 were female specific, 143 were male specific, and 144 were expressed in both male and female antennae. Five odorant-binding proteins, two chemosensory proteins, and one olfactory receptor were identified. They may play a major role in the perception of odorants. An esterase and an aldehyde dehydrogenase were found only in male antennae. Glutathione S-transferases (GSTs) and cytochrome P450s, also found in silkworm antennae, may be involved in the degradation of xenobiotics. Additionally, antioxidation proteins and immunity proteins were identified. Juvenile hormone binding proteins (JHBP), juvenile hormone resistance protein II, and juvenile hormone episode hydrolase (JHEH) were found in the proteomic analysis, which suggests that the antennae are a target for juvenile hormone in the silkworm. Our results provide insight into the expression of proteins in the antennae of silkworm and will facilitate the future functional analysis of silkworm antennae.

Culture-expanded human mesenchymal stem cells (MSCs) are increasingly used in clinics, yet full characterization of the genomic compositions of these cells is lacking. We present a whole-genome investigation on the genetic dynamics of cultured MSCs under ex vivo establishment (passage 1 [p1]) and serial expansion (p8 and p13). We detected no significant changes in copy-number alterations (CNAs) and low levels of single-nucleotide changes (SNCs) until p8. Strikingly, a significant number (677) of SNCs were found in p13 MSCs. Using a sensitive Droplet Digital PCR assay, we tested the nonsynonymous SNCs detected by whole-genome sequencing and found that they were preexisting low-frequency mutations in uncultured mononuclear cells (∼0.01%) and early-passage MSCs (0.1%-1% at p1 and p8) but reached 17%-36% in p13. Our data demonstrate that human MSCs maintain a stable genomic composition in the early stages of ex vivo culture but are subject to clonal growth upon extended expansion., (Copyright © 2014 The Authors. Published by Elsevier Inc. All rights reserved.)

Motivation: Copy number abbreviation (CNA) is one type of genomic aberration that is often induced by genome instability and is associated with diseases such as cancer. Determination of the genome-wide CNA profile is an important step in identifying the underlying mutation mechanisms. Genomic data based on next-generation sequencing technology are particularly suitable for determination of high-quality CNA profile. Now is an important time to reevaluate the use of sequencing techniques for CNA analysis, especially with the rapid growth of the different targeted genome and whole-genome sequencing strategies., Results: In this study, we provide a comparison of resequencing strategies, with regard to their utility, applied to the same hepatocellular carcinoma sample for copy number determination. These strategies include whole-genome, exome and restriction site-associated DNA (RAD) sequencing. The last of these strategies is a targeted sequencing technique that involves cutting the genome with a restriction enzyme and isolating the targeted sequences. Our data demonstrate that RAD sequencing is an efficient and comprehensive strategy that allows the cost-effective determination of CNAs. Further investigation of RAD sequencing data led to the finding that a precise measurement of the allele frequency would be a helpful complement to the read depth for CNA analysis for two reasons. First, knowledge of the allele frequency helps to resolve refined calculations of allele-specific copy numbers, which, in turn, identify the functionally important CNAs that are under natural selection on the parental alleles. Second, this knowledge enables deconvolution of CNA patterns in complex genomic regions.

Background: Cotton bollworm (Helicoverpa armigera) and oriental tobacco budworm (Helicoverpa assulta) are noctuid sibling species. Under artificial manipulation, they can mate and produce fertile offspring. As serious agricultural insect pests, cotton bollworms are euryphagous insects, but oriental tobacco budworms are oligophagous insects. To identify the differentially expressed genes that affect host recognition and host adaptation between the two species, we constructed digital gene expression tag profiles for four developmental stages of the two species. High-throughput sequencing results indicated that we have got more than 23 million 17nt clean tags from both species, respectively. The number of unique clean tags was nearly same in both species (approximately 357,000)., Results: According to the gene annotation results, we identified 83 and 68 olfaction related transcripts from H. armigera and H. assulta, respectively. At the same time, 1137 and 1138 transcripts of digestion enzymes were identified from the two species. Among the olfaction related transcripts, more odorant binding protein and G protein-coupled receptor were identified in H. armigera than in H. assulta. Among the digestion enzymes, there are more detoxification enzyme, e.g. P450, carboxypeptidase and ATPase in H. assulta than in H. armigera. These differences partially explain that because of the narrow host plant range of H. assulta, more detoxification enzymes would help them increase the food detoxification and utilization efficiency., Conclusions: This study supplied some differentially expressed genes affecting host selection and adaptation between the two sibling species. These genes will be useful information for studying on the evolution of host plant selection. It also provides some important target genes for insect species-specific control by RNAi technology.

Background: Brown planthopper (BPH), Nilaparvata lugens Stål, is one of the most destructive insect pests of rice. The molecular responses of plants to sucking insects resemble responses to pathogen infection. However, the molecular mechanism of BPH-resistance in rice remains unclear. Transcription factors (TF) are up-stream regulators of various genes that bind to specific DNA sequences, thereby controlling the transcription from DNA to mRNA. They are key regulators for transcriptional expression in biological processes, and are probably involved in the BPH-induced pathways in resistant rice varieties., Results: We conducted a microarray experiment to analyze TF genes related to BPH resistance in a Sri Lankan rice cultivar, Rathu Heenati (RHT). We compared the expression profiles of TF genes in RHT with those of the susceptible rice cultivar Taichun Native 1 (TN1). We detected 2038 TF genes showing differential expression signals between the two rice varieties. Of these, 442 TF genes were probably related to BPH-induced resistance in RHT and TN1, and 229 may be related to constitutive resistance only in RHT. These genes showed a fold change (FC) of more than 2.0 (P<0.05). Among the 442 TF genes related to BPH-induced resistance, most of them were readily induced in TN1 than in RHT by BPH feeding, for instance, 154 TF genes were up-regulated in TN1, but only 31 TF genes were up-regulated in RHT at 24 hours after BPH infestation; 2-4 times more TF genes were induced in TN1 than in RHT by BPH. At an FC threshold of >10, there were 37 induced TF genes and 26 constitutive resistance TF genes. Of these, 13 were probably involved in BPH-induced resistance, and 8 in constitutive resistance to BPH in RHT., Conclusions: We explored the molecular mechanism of resistance to BPH in rice by comparing expressions of TF genes between RHT and TN1. We speculate that the level of gene repression, especially for early TF genes, plays an important role in the defense response. The fundamental point of the resistance strategy is that plants protect themselves by reducing their metabolic level to inhibit feeding by BPH and prevent damage from water and nutrient loss. We have selected 21 TF genes related to BPH resistance for further analyses to understand the molecular responses to BPH feeding in rice.

MucK is a member of the major facilitator super family (MFS) and is within a subgroup involved in the uptake of organic acids. Here, we provide the first evidence that mucK is required for normal body color pattern in insect larvae. In the cts mutant strain of Bombyx mori (Linnaeus), the larval anal plates, prothoracic tergum, and the head are reddish brown, as opposed to the white color found in the wild type. Positional cloning of the cts gene was carried out and it was mapped within 90kb on the 16th chromosome. Two molecular markers, 182A4 and 242GA, were found co-segregated with the cts phenotype. Only one candidate gene in the region, BGIBMGA013242, which codes for a major facilitator super family protein named BmMucK, located between the two loci. Sequence analysis revealed that a ∼3kb deletion in the genome resulting in a 126bp deletion in the open reading frame (ORF) of BmmucK in the cts mutant. Fluorescence in situ hybridization (FISH) and quantitative reverse transcription PCR showed that BmmucK expression coincided with reddish brown pigmented regions in the cts mutant strain. BmmucK expression in the head and anal plates were markedly elevated at the molting stages. We performed RNA interference of the mucK homolog in Tribolium castaneum and observed brown coloration accompanied by lethality in ∼1/3 of the insects. The results suggest that MucK plays an important role in the pigmentation process of insect larvae., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

The green alga Haematococcus pluvialis can produce large amounts of pink carotenoid astaxanthin which is a high value ketocarotenoid. In our study, transcriptional expression patterns of eight carotenoid genes in H. pluvialis in response to SA were measured using qRT-PCR. Results indicated that both 25 and 50 mg/L salicylic acid (SA) could increase astaxanthin productivity and enhance transcriptional expression of eight carotenoid genes in H. pluvialis. But these genes exhibited different expression profiles. Moreover, SA25 (25 mg/L SA) induction had a greater effect on the transcriptional expression of ipi-1, psy, pds, crtR-B and lyc (more than 6-fold up-regulation) than on ipi-2, bkt and crtO, but SA50 (50 mg/L SA) treatment had a greater impact on the transcriptional expression of ipi-1, ipi-2, pds, crtR-B and lyc than on psy, bkt and crtO. Furthermore, astaxanthin biosynthesis under SA was up-regulated mainly by ipi-1, ipi-2, psy, crtR-B, bkt and crtO at transcriptional level, lyc at post-transcriptional level and pds at both levels. Summarily, these results suggest that SA constitute molecular signals in the network of astaxanthin biosynthesis. Induction of astaxanthin accumulation by SA without any other stimuli presents an attractive application potential in astaxanthin production with H. pluvialis., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Rathu Heenati (RHT) is a Sri Lankan rice cultivar that carries a brown planthopper (BPH) resistance gene, Bph3, and shows broad-spectrum resistance to all four biotypes of BPH. The BPH-resistance loci in RHT has been studied extensively and assigned to four different rice chromosomes (3, 4, 6, and 10) by different research groups, but the gene has not been cloned previously. An Affymetrix rice genome array containing 48,564 japonica and 1,260 indica sequences was used to analyze the potential resistance-related genes on the four chromosomes by comparative analysis of the differentially expressed genes between resistant and susceptible rice cultivars exposed to BPH attack. The microarray results showed that at least 17 genes related to induced resistance and at least 193 genes related to constitutive resistance in RHT. On chromosome 3, the AOC4 was hypothesized to be the most important candidate gene. On chromosome 6, no valuable candidate resistance gene was identified in the Bph3 localization region. In the three Quantitative trait locus regions of chromosomes 3, 4, and 10, the numbers of constitutive and induced resistance-related genes found were 17, 26, and 12, respectively. The major probe on chromosome 10 represents a constitutive expression gene with a very high absolute fold-change of 2,588.82. The microarray analysis indicated that BPH resistance in RHT is probably controlled by a series of resistance-related genes. This study provides valuable information for cloning, functional analysis and marker-assisted breeding of these BPH resistance genes.

The key of RNAi approach success for potential insect pest control is mainly dependent on careful target selection and a convenient delivery system. We adopted second-generation sequencing technology to screen RNAi targets. Illumina's RNA-seq and digital gene expression tag profile (DGE-tag) technologies were used to screen optimal RNAi targets from Ostrinia furnalalis. Total 14690 stage specific genes were obtained which can be considered as potential targets, and 47 were confirmed by qRT-PCR. Ten larval stage specific expression genes were selected for RNAi test. When 50 ng/µl dsRNAs of the genes DS10 and DS28 were directly sprayed on the newly hatched larvae which placed on the filter paper, the larval mortalities were around 40∼50%, while the dsRNAs of ten genes were sprayed on the larvae along with artificial diet, the mortalities reached 73% to 100% at 5 d after treatment. The qRT-PCR analysis verified the correlation between larval mortality and the down-regulation of the target gene expression. Topically applied fluorescent dsRNA confirmed that dsRNA did penetrate the body wall and circulate in the body cavity. It seems likely that the combination of DGE-tag with RNA-seq is a rapid, high-throughput, cost less and an easy way to select the candidate target genes for RNAi. More importantly, it demonstrated that dsRNAs are able to penetrate the integument and cause larval developmental stunt and/or death in a lepidopteron insect. This finding largely broadens the target selection for RNAi from just gut-specific genes to the targets in whole insects and may lead to new strategies for designing RNAi-based technology against insect damage.

The pigmentation of insects has served as an excellent model for the study of morphological trait evolution and developmental biology. The melanism (mln) mutant of the silkworm Bombyx mori is notable for its strong black coloration, phenotypic differences between larval and adult stages, and its widespread use in strain selection. Here, we report the genetic and molecular bases for the formation of the mln morphological trait. Fine mapping revealed that an arylalkylamine N-acetyltransferase (AANAT) gene co-segregates with the black coloration patterns. Coding sequence variations and expression profiles of AANAT are also associated with the melanic phenotypes. A 126 bp deletion in the mln genome causes two alternatively spliced transcripts with premature terminations. An enzymatic assay demonstrated the absolute loss of AANAT activity in the mutant proteins. We also performed RNA interference of AANAT in wild-type pupae and observed a significant proportion of adults with ectopic black coloration. These findings indicate that functional deletion of this AANAT gene accounts for the mln mutation in silkworm. AANAT is also involved in a parallel melanin synthesis pathway in which ebony plays a role, whereas no pigmentation defect has been reported in the Drosophila model or in other insects to date. To the best of our knowledge, the mln mutation is the first characterized mutant phenotype of insects with AANAT, and this result contributes to our understanding of dopamine metabolism and melanin pattern polymorphisms.

Background: Bombyx mori, the domesticated silkworm, is a well-studied model insect with great economic and scientific significance. Although more than 400 mutations have been described in silkworms, most have not been identified, especially those affecting economically-important traits. Simple sequence repeats (SSRs) are effective and economical tools for mapping traits and genetic improvement. The current SSR linkage map is of low density and contains few polymorphisms. The purpose of this work was to develop a dense and informative linkage map that would assist in the preliminary mapping and dissection of quantitative trait loci (QTL) in a variety of silkworm strains., Results: Through an analysis of > 50,000 genotypes across new mapping populations, we constructed two new linkage maps covering 27 assigned chromosomes and merged the data with previously reported data sets. The integrated consensus map contains 692 unique SSR sites, improving the density from 6.3 cM in the previous map to 4.8 cM. We also developed 497 confirmed neighboring markers for corresponding low-polymorphism sites, with 244 having polymorphisms. Large-scale statistics on the SSR type were suggestive of highly efficient markers, based upon which we searched 16,462 available genomic scaffolds for SSR loci. With the newly constructed map, we mapped single-gene traits, the QTL of filaments, and a number of ribosomal protein genes., Conclusion: The integrated map produced in this study is a highly efficient genetic tool for the high-throughput mapping of single genes and QTL. Compared to previous maps, the current map offers a greater number of markers and polymorphisms; thus, it may be used as a resource for marker-assisted breeding.

MicroRNAs (miRNAs) are endogenous non-coding genes that participate in post-transcription regulation by either degrading mRNA or blocking its translation. It is considered to be very important in regulating insect development and metamorphosis. We conducted a large-scale screening for miRNA genes in the silkworm Bombyx mori using sequence-by-synthesis (SBS) deep sequencing of mixed RNAs from egg, larval, pupal, and adult stages. Of 2,227,930 SBS tags, 1,144,485 ranged from 17 to 25 nt, corresponding to 256,604 unique tags. Among these non-redundant tags, 95,184 were matched to the silkworm genome. We identified 3,750 miRNA candidate genes using a computational pipeline combining RNAfold and TripletSVM algorithms. We confirmed 354 miRNA genes using miRNA microarrays and then performed expression profile analysis on these miRNAs for all developmental stages. While 106 miRNAs were expressed in all stages, 248 miRNAs were egg- and pupa-specific, suggesting that insect miRNAs play a significant role in embryogenesis and metamorphosis. We selected eight miRNAs for quantitative RT-PCR analysis; six of these were consistent with our microarray results. In addition, we searched for orthologous miRNA genes in mammals, a nematode, and other insects and found that most silkworm miRNAs are conserved in insects, whereas only a small number of silkworm miRNAs has orthologs in mammals and the nematode. These results suggest that there are many miRNAs unique to insects.

Lysozymes act as crucial bacteriolytic enzymes in insect immune system by hydrolyzing the beta (1-->4) bonds between N-acetylglucosamine and N-acetylmuramic acid in the peptidoglycan of prokaryotic cell walls. We have isolated and characterized a Helicoverpa armigera cDNA encoding an insect lysozyme named HaLyz. We amplified a fragment by PCR, using degenerate primers derived from the conservative amino acid sequences for performing 5' and 3' RACE. The full-length cDNA was 661 base pairs. The theoretical pI and molecular weight of the protein were computed to be 9.08 and 15.6 kDa, respectively. Prokaryotic expression of the HaLyz ORF by Escherichia coli confirmed the calculated molecular weight of the protein. The deduced 135 amino acids showed high homology with known lysozymes from other insects, ranging from 47% to 89% by BLASTp search in NCBI. Analyses revealed that this protein has a typical lysozyme C signature among amino acids 93-111, CNVTCAEMLLDDITKASTC. An interesting relation between immunity and larva to pupa metamorphosis in insects was discovered. Real time-PCR showed that HaLyz gene expression was transiently enhanced at the onset of metamorphosis of the cotton bollworm, Helicoverpa armigera. The gene expression was up-regulated after the injection of E. coli or entomopathogenic fungi, Beauveria bassiana, but showed different expression patterns.

Pair-rule genes (genes that are expressed only in alternate segments, odd or even) play an important role in translating the broad gradients of upstream genes into dual segment periodicity for body plan patterning in Drosophila. However, homologues of pair-rule genes show a remarkable diversity of expression patterns and functions in other insects. We cloned the homologue of runt in the silkworm Bombyx mori, an intermediate germband-type insect. Whole-mount in situ hybridization revealed three stripes arose one by one before gastrulation at the blastoderm stage. Five additional stripes were then generated sequentially as the growth zone elongated. Eight stripes appeared in a pair-rule manner with two-segment periodicity, each of which was confined to the posterior of an odd-numbered parasegment. The weaker segmental secondary stripes emerged de novo in even-numbered parasegments. The Bmrunt transcript vanished before blastokinesis and was then expressed again in the whole embryo. RNA interference for Bmrunt caused severely truncated, almost completely asegmental defects. This cadual-like phenotype suggests that Bmrunt does not function as a pair-rule gene in silkworm segmentation. Bmrunt is required for formation of most body segments and axis elongation in B. mori., (Copyright 2008 Wiley-Liss, Inc.)

Many genes act together during the complex process of insect larval and pupal development. 20-hydroxyecdysone interacts with juvenile hormone to control insect growth and development and then activates several transcription factors, i.e., Broad, E74, and E75, and, subsequently, the late target genes. To investigate this phenomenon, we used serial analysis of gene expression (SAGE) tag-based cDNA microarray analysis to monitor the global gene expression profile during larval development and larva-pupa metamorphosis of the silkworm Bombyx mori. Of the 330 clones that were dotted to the chip, 267 were obtained by generating longer cDNA fragments from SAGE tags for gene identification, and the others were obtained from SAGE tag-matched genes or expressed sequence tags from public databases. According to the gene expression profile, the genes were classified into 12 clusters using a self-organizing map analysis. The results were partially confirmed using real-time reverse transcription-polymerase chain reaction. We obtained 22 full-length cDNAs using rapid amplification of 5' cDNA ends, of which eight genes were novel in the silkworm. Our results indicated that use of a cDNA microarray based on SAGE tags is effective for identifying and examining some low-expression genes associated with insect development.

Necrophagous fly is the most common insect evidence collected during a death investigation. We analyzed the DNA polymorphism among the five forensically important fly species, namely, Lucilia sericata, Aldrichina grahami, Chrysomya megacephala, Parasarcophaga crassipalpis and Musca domestica using inter simple sequence repeat (ISSR) method. Nine ISSR primers selected from 18 primers could amplify 105 clear and stable bands, of which 95 bands were polymorphic. Some primers produced completely different band pattern in different species, indicating that they can be used to identify these species. Aiming at obtaining more reliable markers that might be universally used, we started an effort to convert species-specific ISSR fragments into the sequence-characterized amplified region (SCAR) markers that can be used for the molecular diagnosis of the five species.

Intersimple sequence repeat (ISSR) amplification was used to analyze genetic relationships among silkworm, Bombyx mori L., strains. Nineteen primers containing simple sequence repeat (SSR) motifs were tested for amplification on a panel of 42 strains, representative of the diversity of silkworm germplasm; 12 of the primers amplified distinct, reproducible bands. The primers amplified a total of 108 bands, of which 85 (78.7%) were polymorphic. The ISSR results suggested that within the dinucleotide class, the poly(CA) motif was more common than the poly(CT) motif. The ISSR amplification pattern was used to group the silkworm strains into seven subclusters based on their origin in an unweighted pair-group method with arithmetic average cluster analysis by using Nei's genetic distance. Seven major ecotypic silkworm groups were analyzed. Principal component analysis of the ISSR data supported the unweighted pair-group method with arithmetic average clustering. Therefore, ISSR amplification is a valuable method for determining genetic variability among silkworm varieties. This efficient genetic fingerprinting technique should be useful for characterizing the large numbers of silkworm strains held in national and international germplasm centers.