Transcriptome profile in Drosophila Kc and S2 embryonic cell lines

Drosophila melanogaster cell lines are an important resource for a range of studies spanning genomics, molecular genetics, and cell biology. Amongst these valuable lines are Kc167 (Kc) and Schneider 2 (S2) cells, which were originally isolated in the late 1960s from embryonic sources and have been used extensively to investigate a broad spectrum of biological activities including cell-cell signaling and immune system function. Whole-genome tiling microarray analysis of total RNA from these two cell types was performed as part of the modENCODE project over a decade ago and revealed that they share a number of gene expression features. Here, we expand on these earlier studies by using deep-coverage RNA-sequencing approaches to investigate the transcriptional profile in Kc and S2 cells in detail. Comparison of the transcriptomes reveals that ∼75% of the 13,919 annotated genes are expressed at a detectable level in at least one of the cell lines, with the majority of these genes expressed at high levels in both cell lines. Despite the overall similarity of the transcriptional landscape in the two cell types, 2,588 differentially expressed genes are identified. Many of the genes with the largest fold change are known only by their CG designations, indicating that the molecular control of Kc and S2 cell identity may be regulated in part by a cohort of relatively uncharacterized genes. Our data also indicate that both cell lines have distinct hemocyte-like identities, but share active signaling pathways and express a number of genes in the network responsible for dorsal-ventral patterning of the early embryo. © The Author(s) 2023. Published by Oxford University Press on behalf of the Genetics Society of America

Global sensitivity analysis of a dynamic model for gene expression in Drosophila embryos

It is well known that gene regulation is a tightly controlled process in early organismal development. However, the roles of key processes involved in this regulation, such as transcription and translation, are less well understood, and mathematical modeling approaches in this field are still in their infancy. In recent studies, biologists have taken precise measurements of protein and mRNA abundance to determine the relative contributions of key factors involved in regulating protein levels in mammalian cells. We now approach this question from a mathematical modeling perspective. In this study, we use a simple dynamic mathematical model that incorporates terms representing transcription, translation, mRNA and protein decay, and diffusion in an early Drosophila embryo. We perform global sensitivity analyses on this model using various different initial conditions and spatial and temporal outputs. Our results indicate that transcription and translation are often the key parameters to determine protein abundance. This observation is in close agreement with the experimental results from mammalian cells for various initial conditions at particular time points, suggesting that a simple dynamic model can capture the qualitative behavior of a gene. Additionally, we find that parameter sensitivites are temporally dynamic, illustrating the importance of conducting a thorough global sensitivity analysis across multiple time points when analyzing mathematical models of gene regulation

Characterization of the intergenic RNA profile at abdominal-A and Abdominal-B in the Drosophila bithorax complex

The correct spatial expression of two Drosophila bithorax complex (BX-C) genes, abdominal-A (abdA) and Abdominal-B (AbdB), is dependent on the 100-kb intergenic infraabdominal (iab) region. The iab region is known to contain a number of different domains (iab2 through iab8) that harbor cis-regulatory elements responsible for directing expression of abdA and AbdB in the second through eighth abdominal segments. Here, we use in situ hybridization to perform high-resolution mapping of the transcriptional activity in the iab control regions. We show that transcription of the control regions themselves is abundant and precedes activation of the abdA and AbdB genes. As with the homeotic genes of the BX-C, the transcription patterns of the RNAs from the iab control regions demonstrate colinearity with the sequence of the iab regions along the chromosome and the domains in the embryo under the control of the specific iab regions. These observations suggest that the intergenic RNAs may play a role in initiating cis regulation at the BX-C early in development