Spatial Expression Patterns Research Articles

Genetic dissection of quantitative trait loci for spikelets compactness in two Yanzhan1‐derived recombinant inbred line wheat populations

AbstractSpikelets compactness (SC) plays an important role in the yield improvement and Fusarium head blight (FHB) resistance of wheat. In this study, 151 recombinant inbred lines (RIL) derived from a cross between Yangmai 4 (YM4) and Yanzhan 1 (YZ1) and 143 RIL derived from a cross between Yangmai 5 (YM5) and YZ1 were used to construct two wheat hexaploidy genetic linkage maps, evaluate SC across three years and identify quantitative trait loci (QTL) for SC. A total of seven QTL were identified in two populations and explained the phenotypic variation ranging from 6.27% to 22.49%. Four QTL were co‐located in two genomic regions on chromosomes 2D and 6B, namely, QSc‐2D and QSc‐6B, respectively. Comparison analysis showed that QSc‐6B is a novel QTL, whereas QSc‐2D was located in the interval of Rht8. Significant additive effects of QSc‐2D and QSc‐6B on SC were observed in two RIL populations. Moreover, the positive alleles of these two QTL extremely significantly increased grain number per spike (GNS) but did not have a significantly negative effect on thousand‐grain weight (TGW), kernel length (KL) and kernel weight (KW). Whether in spray or point inoculation, QSc‐6B did not significantly influence FHB resistance in two populations. However, QSc‐2D significantly affected FHB resistance in spray inoculation in YM4/YZ1 population. A breeder‐friendly Kompetitive Allele‐Specific PCR (KASP) marker KASP_6B for QSc‐6B has been developed and verified in a set of 181 wheat cultivars/lines. Furthermore, we preliminary analysed the candidate genes within the physical interval of QSc‐6B by sequence similarity, spatial expression patterns and collinearity analysis. TraesCS6B03G1206900, TraesCS6B03G1214800 and TraesCS6B03G1218300 show specific expression in spikelet. These results provide valuable information for map‐based cloning of QSc‐6B and a useful tool for wheat breeding.

MicroRNA Biology on Brain Development and Neuroimaging Approach.

Proper brain development requires the precise coordination and orchestration of various molecular and cellular processes and dysregulation of these processes can lead to neurological diseases. In the past decades, post-transcriptional regulation of gene expression has been shown to contribute to various aspects of brain development and function in the central nervous system. MicroRNAs (miRNAs), short non-coding RNAs, are emerging as crucial players in post-transcriptional gene regulation in a variety of tissues, such as the nervous system. In recent years, miRNAs have been implicated in multiple aspects of brain development, including neurogenesis, migration, axon and dendrite formation, and synaptogenesis. Moreover, altered expression and dysregulation of miRNAs have been linked to neurodevelopmental and psychiatric disorders. Magnetic resonance imaging (MRI) is a powerful imaging technology to obtain high-quality, detailed structural and functional information from the brains of human and animal models in a non-invasive manner. Because the spatial expression patterns of miRNAs in the brain, unlike those of DNA and RNA, remain largely unknown, a whole-brain imaging approach using MRI may be useful in revealing biological and pathological information about the brain affected by miRNAs. In this review, we highlight recent advancements in the research of miRNA-mediated modulation of neuronal processes that are important for brain development and their involvement in disease pathogenesis. Also, we overview each MRI technique, and its technological considerations, and discuss the applications of MRI techniques in miRNA research. This review aims to link miRNA biological study with MRI analytical technology and deepen our understanding of how miRNAs impact brain development and pathology of neurological diseases.

Comparing the development of cortex-wide gene expression patterns between two species in a common reference frame

Advances in sequencing techniques have made comparative studies of gene expression a current focus for understanding evolutionary and developmental processes. However, insights into the spatial expression of genes have been limited by a lack of robust methodology. To overcome this obstacle, we developed methods and software tools for quantifying and comparing tissue-wide spatial patterns of gene expression within and between species. Here, we compare cortex-wide expression of RZRβ and Id2 mRNA across early postnatal development in mice and voles. We show that patterns of RZRβ expression in neocortical layer 4 are highly conserved between species but develop rapidly in voles and much more gradually in mice, who show a marked expansion in the relative size of the putative primary visual area across the first postnatal week. Patterns of Id2 expression, by contrast, emerge in a dynamic and layer-specific sequence that is consistent between the two species. We suggest that these differences in the development of neocortical patterning reflect the independent evolution of brains, bodies, and sensory systems in the 35 million years since their last common ancestor.

A model-based constrained deep learning clustering approach for spatially resolved single-cell data.

Spatially resolved scRNA-seq (sp-scRNA-seq) technologies provide the potential to comprehensively profile gene expression patterns in tissue context. However, the development of computational methods lags behind the advances in these technologies, which limits the fulfillment of their potential. In this study, we develop a deep learning approach for clustering sp-scRNA-seq data, named Deep Spatially constrained Single-cell Clustering (DSSC). In this model, we integrate the spatial information of cells into the clustering process in two steps: (1) the spatial information is encoded by using a graphical neural network model, and (2) cell-to-cell constraints are built based on the spatial expression pattern of the marker genes and added in the model to guide the clustering process. Then, a deep embedding clustering is performed on the bottleneck layer of autoencoder by Kullback-Leibler (KL) divergence along with the learning of feature representation. DSSC is the first model that can use information from both spatial coordinates and marker genes to guide cell/spot clustering. Extensive experiments on both simulated and real data sets show that DSSC boosts clustering performance significantly compared with the state-of-the-art methods. It has robust performance across different data sets with various cell type/tissue organization and/or cell type/tissue spatial dependency. We conclude that DSSC is a promising tool for clustering sp-scRNA-seq data.

TIST: Transcriptome and Histopathological Image Integrative Analysis for Spatial Transcriptomics

Sequencing-based spatial transcriptomics (ST) is an emerging technology to study in situ gene expression patterns at the whole-genome scale. Currently, ST data analysis is still complicated by high technical noises and low resolution. In addition to the transcriptomic data, matched histopathological images are usually generated for the same tissue sample along the ST experiment. The matched high-resolution histopathological images provide complementary cellular phenotypical information, providing an opportunity to mitigate the noises in ST data. We present a novel ST data analysis method called transcriptome and histopathological image integrative analysis for ST (TIST), which enables the identification of spatial clusters (SCs) and the enhancement of spatial gene expression patterns by integrative analysis of matched transcriptomic data and images. TIST devises a histopathological feature extraction method based on Markov random field (MRF) to learn the cellular features from histopathological images, and integrates them with the transcriptomic data and location information as a network, termed TIST-net. Based on TIST-net, SCs are identified by a random walk-based strategy, and gene expression patterns are enhanced by neighborhood smoothing. We benchmark TIST on both simulated datasets and 32 real samples against several state-of-the-art methods. Results show that TIST is robust to technical noises on multiple analysis tasks for sequencing-based ST data and can find interesting microstructures in different biological scenarios. TIST is available at http://lifeome.net/software/tist/ and https://ngdc.cncb.ac.cn/biocode/tools/BT007317.

Zebrafish Slit2 and Slit3 Act Together to Regulate Retinal Axon Crossing at the Midline

Slit-Robo signaling regulates midline crossing of commissural axons in different systems. In zebrafish, all retinofugal axons cross at the optic chiasm to innervate the contralateral tectum. Here, the mutant for the Robo2 receptor presents severe axon guidance defects, which were not completely reproduced in a Slit2 ligand null mutant. Since slit3 is also expressed around this area at the stage of axon crossing, we decided to analyze the possibility that it collaborates with Slit2 in this process. We found that the disruption of slit3 expression by sgRNA-Cas9 injection caused similar, albeit slightly milder, defects than those of the slit2 mutant, while the same treatment in the slit2−/−mz background caused much more severe defects, comparable to those observed in robo2 mutants. Tracking analysis of in vivo time-lapse experiments indicated differential but complementary functions of these secreted factors in the correction of axon turn errors around the optic chiasm. Interestingly, RT-qPCR analysis showed a mild increase in slit2 expression in slit3-deficient embryos, but not the opposite. Our observations support the previously proposed “repulsive channel” model for Slit-Robo action at the optic chiasm, with both Slits acting in different manners, most probably relating to their different spatial expression patterns.

No functional contribution of the gustatory receptor, Gr64b, co-expressed in olfactory sensory neurons of Drosophila melanogaster

Chemosensation is essential for the survival of insects. Activities like searching for food, mating, and oviposition in the fruit fly, Drosophila melanogaster are to a great extent governed by chemical cues detected via olfaction and gustation. This chemical information is conveyed to higher brain centers via populations of diverse olfactory sensory neurons (OSNs) and gustatory sensory neurons (GSNs) expressing olfactory receptors (ORs) and gustatory receptors (GRs), respectively. ORs are exclusively expressed in the antenna and in the maxillary palps, while GRs are widely expressed in the labellum, tarsi, genitalia etc. Interestingly, 14 GRs were previously reported to be expressed in the antenna of D. melanogaster. However, the spatial expression pattern for all GRs and their functional role are still unclear. Recent data challenge the dogma that single OSNs express a single OR. In the present study, we studied the expression of 12 previously reported GRs among sensory structures on the fly antenna using the Gal4-UAS binary expression system. We observed antennal expression of nine out of the 12 reported. Out of these nine, consistent expression was only apparent for Gr64b, and we reconfirmed its presence in OSNs innervating three glomeruli in the antennal lobe. These glomeruli are known to be innervated by ab5A, ab5B and ab8A OSNs, respectively. Next, we generated double labeling crosses with Gr64b and observed co-expression of Gr64b with Or47a, which is expressed in the ab5B neuron. To elucidate the functional role of Gr64b co-expressed with Or47a, we challenged Or47a-expressing OSNs in wild type and Gr64b–/– mutant flies with odor stimulation using the single sensillum recording technique in two satiation states (fed and starved). Notably, we did not observe any significant odor sensitivity or specificity changes in Gr64b mutants as compared to wild type flies. Taken together, our results reveal co-expression of GRs with ORs in olfactory sensory neurons, while the functional contribution of the GR in this context remains obscure.

Scalable and model-free detection of spatial patterns and colocalization.

The expeditious growth in spatial omics technologies enables the profiling of genome-wide molecular events at molecular and single-cell resolution, highlighting a need for fast and reliable methods to characterize spatial patterns. We developed SpaGene, a model-free method to discover spatial patterns rapidly in large-scale spatial omics studies. Analyzing simulation and a variety of spatially resolved transcriptomics data showed that SpaGene is more powerful and scalable than existing methods. Spatial expression patterns identified by SpaGene reconstruct unobserved tissue structures. SpaGene also successfully discovers ligand–receptor interactions through their colocalization.

Exploring the potential for an evolutionarily conserved role of the taste 1 receptor gene family in gut sensing mechanisms of fish.

In this study, we investigated the transcriptional spatio-temporal dynamics of the taste 1 receptor (T1R) gene family repertoire in seabream (Sparus aurata [sa]), during larval ontogeny and in adult tissues. In early larval development, saT1R expression arises heterochronously, i.e. the extraoral taste-related perception in the gastrointestinal tract (GIT) anticipates first exogenous feeding (at 9 days post hatching [dph]), followed by the buccal/intraoral perception from 14 dph onwards, supporting the hypothesis that the early onset of the molecular machinery underlying saT1R expression in the GIT is not induced by food but rather genetically hardwired. During adulthood, we characterized the expression patterns of saT1R within specific tissues (n = 4) distributed in oropharingeal, GIT and brain regions substantiating their functional versatility as chemosensory signaling players to a variety of biological functions beyond oral taste sensation. Further, we provided for the first time direct evidences in fish for mRNA co-expression of a subset of saT1R genes (mostly saT1R3, i.e. the common subunit of the heterodimeric T1R complexes for the detection of “sweet” and “umami” substances), with the selected gut peptides ghrelin (ghr), cholecystokinin (cck), hormone peptide yy (pyy) and proglucagon (pg). Each peptide defines the enteroendocrine cells (ECCs) identity, and establishes on morphological basis, a direct link for T1R chemosensing in the regulation of fish digestive processes. Finally, we analyzed the spatial gene expression patterns of 2 taste signaling components functionally homologous to the mammalian G(i)α subunit gustducin, namely saG(i)α1 and saG(i)α2, and demonstrated their co-localization with the saT1R3 in EECs, thus validating their direct involvement in taste-like transduction mechanisms of the fish GIT. In conclusion, data provide new insights in the evolutionary conservation of gut sensing in fish suggesting a conserved role for nutrient sensors modulating entero-endocrine secretion.

Spatially patterned hydrogen peroxide orchestrates stomatal development in Arabidopsis

Stomatal pores allow gas exchange between plant and atmosphere. Stomatal development is regulated by multiple intrinsic developmental and environmental signals. Here, we show that spatially patterned hydrogen peroxide (H2O2) plays an essential role in stomatal development. H2O2 is remarkably enriched in meristemoids, which is established by spatial expression patterns of H2O2-scavenging enzyme CAT2 and APX1. SPEECHLESS (SPCH), a master regulator of stomatal development, directly binds to the promoters of CAT2 and APX1 to repress their expression in meristemoid cells. Mutations in CAT2 or APX1 result in an increased stomatal index. Ectopic expression of CAT2 driven by SPCH promoter significantly inhibits the stomatal development. Furthermore, H2O2 activates the energy sensor SnRK1 by inducing the nuclear localization of the catalytic α-subunit KIN10, which stabilizes SPCH to promote stomatal development. Overall, these results demonstrate that the spatial pattern of H2O2 in epidermal leaves is critical for the optimal stomatal development in Arabidopsis.

Duplication of a manganese/cadmium transporter gene reduces cadmium accumulation in rice grain.

Global contamination of soils with toxic cadmium (Cd) is a serious health threat. Here we found that a tandem duplication of a gene encoding a manganese/Cd transporter, OsNramp5, was responsible for low-Cd accumulation in Pokkali, an old rice cultivar. This duplication doubled the expression of OsNramp5 gene but did not alter its spatial expression pattern and cellular localization. Higher expression of OsNramp5 increased uptake of Cd and Mn into the root cells but decreased Cd release to the xylem. Introgression of this allele into Koshihikari, an elite rice cultivar, through backcrossing significantly reduced Cd accumulation in the grain when cultivated in soil heavily contaminated with Cd but did not affect both grain yield and eating quality. This study not only reveals the molecular mechanism underlying low-Cd accumulation but also provides a useful target for breeding rice cultivars with low-Cd accumulation.

Empty-spiracles is maternally expressed and essential for neurodevelopment and early embryo determination in Rhodnius prolixus

Since empty-spiracles (ems) was identified and characterized in Drosophila melanogaster as a head-gap gene, several studies have been carried out in other insect orders to confirm its evolutionary conserved function. Using the blood-sucking bug Rhodnius prolixus as biological model, we found an ems transcript with three highly conserved regions: Box-A, Box-B, and the homeodomain. R. prolixus embryos silenced by parental RNAi for two of these ems conserved regions showed both maternal and zygotic defects. Rp-emsB fragment results in early lethal embryogenesis, with eggs without any embryonic structure inside. Rp-emsB expression pattern is only maternally expressed and localized in the ovary tropharium, follicular cells, and in the unfertilized female pronucleus. Rp-emsA fragment is zygotically expressed during early blastoderm formation until late developmental stages in two main patterns: anterior in the antennal segment, and in a segmentary in the neuroblast and tracheal pits. R. prolixus knockdown embryos for Rp-emsA showed an incomplete larval hatching, reduced heads, and severe neuromotor defects. Furthermore, in situ hybridization revealed a spatial and temporal expression pattern that highly correlates with Rp-ems observed function. Here,Rp-ems function in R. prolixus development was validated, showing that empty-spiracles does not act as a true head-gap gene, but it is necessary for proper head development and crucial for early embryo determination and neurodevelopment.

Assessment of a zinc finger protein gene (MPZC3H10) as potential RNAi target for green peach aphid Myzus persicae control.

RNA interference (RNAi) has potential application in pest control, and selection of the specific target gene is one of the key steps in RNAi. As an important effector, the zinc finger protein (ZFP) gene has high similarity among aphid species, and may have potential use in an RNAi-based pest control strategy. This study assessed the control efficiency of an RNAi target, MPZC3H10, a CCCH-type ZFP gene, against green peach aphid. ZC3H10 amino acid sequence similarity is more than 97.71% among the five tested aphid species: Myzus persicae, Aphis citricidus, Acyrthosiphon pisum, Diuraphis noxia and Rhopalosiphum maidis. However, no homologous sequence was found in the transcriptome of their ladybeetle predator, Propylaea japonica. Spatial expression patterns revealed that MPZC3H10 showed high expression in the muscle and fat body of M. persicae. The RNAi bioassay revealed that silencing of MPZC3H10 resulted in high mortality (53.33%) in M. persicae. By contrast, there were no observed negative effects on the growth and development of P. japonica when fed on aphids treated with double-stranded RNA (dsRNA) or injected with a "high dose" of dsRNA. Targeting MPZC3H10 showed promising efficiency for green peach aphid control via artificially designed dsRNA, and was safe for the predatory ladybeetle. © 2022 Society of Chemical Industry.

Synthetic gene circuits take root.

Complex spatial patterns of gene expression are engineered in plants to modulate root morphology.

Expression pattern of NLRC5 in the postnatal mouse brain

Nucleotide oligomerization domain-like receptors (NLRs), belonging to a large family of pattern recognition receptors, participate in the host’s first line of defense against invading pathogens. Caspase recruitment domain containing 5 (NLRC5), the largest member in the NLR family, is demonstrated to be involved in the innate immune response and inflammatory diseases far and wide. Recent studies report that NLRC5 is associated with some central nervous system (CNS) diseases. Besides, NLRC5 is a mastery regulator for the expression of MHC class I both in the immune system and the CNS, while MHC class I is expressed and exerts its function in the brain. Therefore, it is necessary to investigate the expression pattern of NLRC5 in the developing and adult CNS. In our study, postnatal brain sections of C57BL/6 J mice are analyzed for the expression of NLRC5 protein by immunofluorescence. In the postnatal stages of developing telencephalon, NLRC5 exhibits a spatial and temporal expression pattern. NLRC5 is time-specifically expressed in subfields of hippocampus and different layers of prefrontal cortex. Moreover, it is shown that NLRC5 is highly cell type specific. It can be expressed in large quantities by neurons and microglia, but rarely expressed by astrocytes. Taken together, our research is important for further understanding the biological characteristics of NLRC5 and its function in the CNS.

Characterization and analysis of the promoter region of monodehydroascorbate reductase 4 (CpMDAR4) in papaya.

Differential spatial and temporal expression patterns due to regulatory cis-elements and two different isoforms are detected among CpMDAR4 alleles in papaya. The aim of this research was to study the effects of cis-element differences between the X, Y and Yh alleles on the expression of CpMDAR4, a potential candidate gene for sex differentiation in papaya, using a transcriptional reporter system in a model species Arabidopsis thaliana. Possible effects of a retrotransposon insertion in the Y and Yh alleles on the transcription and expression of CpMDAR4 alleles in papaya flowers were also examined. When comparing promoters and cis-regulatory elements among genes in the non-recombining region of the sex chromosomes, paired genes exhibited differences. Our results showed that differences in the promoter sequences of the CpMDAR4 alleles drove the expression of a reporter gene to different flower tissues in Arabidopsis. β-glucuronidase staining analysis of T2 and T3 lines for constructs containing 5' deletions of native Y and Yh allele promoters showed the loss of specific expression of the reporter gene in the anthers, confirming the existence and location of cis-regulatory element POLLEN1LELAT52. The expression analysis of CpMDAR4 alleles in papaya flowers also showed that all alleles are actively expressed in different flower tissues, with the existence of a shorter truncated isoform, with unknown function, for the Y and Yh alleles due to an LTR-RT insertion in the Y and Yh chromosomes. The observed expression patterns in Arabidopsis thaliana flowers and the expression patterns of CpMDAR4 alleles in papaya flowers suggest that MDAR4 might have a role on development of reproductive organs in papaya, and that it constitutes an important candidate for sex differentiation.

A Bayesian modified Ising model for identifying spatially variable genes from spatial transcriptomics data

A recent technology breakthrough in spatial molecular profiling (SMP) has enabled the comprehensive molecular characterizations of single cells while preserving spatial information. It provides new opportunities to delineate how cells from different origins form tissues with distinctive structures and functions. One immediate question in SMP data analysis is to identify genes whose expressions exhibit spatially correlated patterns, called spatially variable (SV) genes. Most current methods to identify SV genes are built upon the geostatistical model with Gaussian process to capture the spatial patterns. However, the Gaussian process models rely on ad hoc kernels that could limit the models' ability to identify complex spatial patterns. In order to overcome this challenge and capture more types of spatial patterns, we introduce a Bayesian approach to identify SV genes via a modified Ising model. The key idea is to use the energy interaction parameter of the Ising model to characterize spatial expression patterns. We use auxiliary variable Markov chain Monte Carlo algorithms to sample from the posterior distribution with an intractable normalizing constant in the model. Simulation studies using both simulated and synthetic data showed that the energy-based modeling approach led to higher accuracy in detecting SV genes than those kernel-based methods. When applied to two real spatial transcriptomics (ST) datasets, the proposed method discovered novel spatial patterns that shed light on the biological mechanisms. In summary, the proposed method presents a new perspective for analyzing ST data.

630 Cdk2: A marker for invasive melanoma

Background: Although cyclin-dependent kinase 2 (CDK2) is known to be upregulated in certain cancers, its spatial expression pattern in melanocytic tumors has remained unknown. Objective: To characterize CDK2 expression in benign and malignant melanocytic proliferations. Methods: We performed CDK2 immunohistochemistry (IHC) on normal skin and a tumor panel comprising over 200 congenital nevi, dysplastic nevi, melanomas in situ, and invasive melanomas. We recorded their staining patterns, intensities, and compared these to expression of melanoma marker PRAME (Preferentially expressed Antigen in Melanoma).

Spatial Transcriptomics as a Novel Approach to Redefine Electrical Stimulation Safety.

Current standards for safe delivery of electrical stimulation to the central nervous system are based on foundational studies which examined post-mortem tissue for histological signs of damage. This set of observations and the subsequently proposed limits to safe stimulation, termed the “Shannon limits,” allow for a simple calculation (using charge per phase and charge density) to determine the intensity of electrical stimulation that can be delivered safely to brain tissue. In the three decades since the Shannon limits were reported, advances in molecular biology have allowed for more nuanced and detailed approaches to be used to expand current understanding of the physiological effects of stimulation. Here, we demonstrate the use of spatial transcriptomics (ST) in an exploratory investigation to assess the biological response to electrical stimulation in the brain. Electrical stimulation was delivered to the rat visual cortex with either acute or chronic electrode implantation procedures. To explore the influence of device type and stimulation parameters, we used carbon fiber ultramicroelectrode arrays (7 μm diameter) and microwire electrode arrays (50 μm diameter) delivering charge and charge density levels selected above and below reported tissue damage thresholds (range: 2–20 nC, 0.1–1 mC/cm2). Spatial transcriptomics was performed using Visium Spatial Gene Expression Slides (10x Genomics, Pleasanton, CA, United States), which enabled simultaneous immunohistochemistry and ST to directly compare traditional histological metrics to transcriptional profiles within each tissue sample. Our data give a first look at unique spatial patterns of gene expression that are related to cellular processes including inflammation, cell cycle progression, and neuronal plasticity. At the acute timepoint, an increase in inflammatory and plasticity related genes was observed surrounding a stimulating electrode compared to a craniotomy control. At the chronic timepoint, an increase in inflammatory and cell cycle progression related genes was observed both in the stimulating vs. non-stimulating microwire electrode comparison and in the stimulating microwire vs. carbon fiber comparison. Using the spatial aspect of this method as well as the within-sample link to traditional metrics of tissue damage, we demonstrate how these data may be analyzed and used to generate new hypotheses and inform safety standards for stimulation in cortex.

Corazonin signaling modulates the synthetic activity of male accessory gland in Grapholita molesta

Although neuropeptide corazonin (Crz) has been identified in numerous insect species, the research about its function in regulation of reproduction is still in its infancy. Herein, we characterized the Crz (GmolCrz) and its receptor (GmolCrzR) to investigate their reproductive function in Grapholita molesta. Both molecular docking result and cell-based receptor activity assay showed that GmolCrz could interact with GmolCrzR. Additionally, spatial expression patterns of GmolCrz and GmolCrzR in males were evaluated. Knockdown of GmolCrz or GmolCrzR significantly lengthened copulation duration and decreased fertility in males. In these males, we found that the production of sperm was normal, while the content of accessory gland proteins (Acps) in the accessory gland (AG) was strongly diminished. Furthermore, knockdown of GmolCrz or GmolCrzR in males had no effect on sperm and Acps transfer to females. RNA-seq and gene expression analyses further confirmed that genes involved in serine-type endopeptidase activity were significantly downregulated in the AG upon GmolCrzR knockdown. Finally, sperm activation assays demonstrated that this process was disrupted in the spermatophore of females mated with GmolCrz or GmolCrzR knockdown males, which may cause the decreased fertility in males. Our findings provide new insights into the functions of Crz signaling in a Lepidopteran insect.