T-Box Genes in Drosophila Limb Development

Protocol for dissecting Drosophila pupae and visualizing RNA expression using hybridization chain reaction

Visualizing RNA expression in the Drosophila epidermis during pupal development is challenging because the tissue is fragile during early pupal stages and increasingly impermeable at later stages. Here, we present a protocol for tissue dissection and detection of RNA in situ. We describe steps for using the hybridization chain reaction (HCR) in early and late Drosophila pupal stages (and larval imaginal discs). This protocol facilitates the study of dynamically changing patterns in gene expression during pupal development.

TBX3 is essential for establishment of the posterior boundary of anterior genes and upregulation of posterior genes together with HAND2 during the onset of limb bud development

During limb bud formation, axis polarities are established as evidenced by the spatially restricted expression of key regulator genes. In particular, the mutually antagonistic interaction between the GLI3 repressor and HAND2 results in distinct and non-overlapping anterior-distal Gli3 and posterior Hand2 expression domains. This is a hallmark of the establishment of antero-posterior limb axis polarity, together with spatially restricted expression of homeodomain and other transcriptional regulators. Here, we show that TBX3 is required for establishment of the posterior expression boundary of anterior genes in mouse limb buds. ChIP-seq and differential gene expression analysis of wild-type and mutant limb buds identifies TBX3-specific and shared TBX3-HAND2 target genes. High sensitivity fluorescent whole-mount in situ hybridisation shows that the posterior expression boundaries of anterior genes are positioned by TBX3-mediated repression, which excludes anterior genes such as Gli3, Alx4, Hand1 and Irx3/5 from the posterior limb bud mesenchyme. This exclusion delineates the posterior mesenchymal territory competent to establish the Shh-expressing limb bud organiser. In turn, HAND2 is required for Shh activation and cooperates with TBX3 to upregulate shared posterior identity target genes in early limb buds.

The transcription factor optomotor-blind restricts apterous expression through TrxG and PcG genes

The establishment of body pattern is a fundamental process in developmental biology. In Drosophila, the wing disc is subdivided into dorsal (D) and ventral (V) compartments by the D/V boundary. The dorsal fate is adopted by expressing the selector gene apterous (ap). ap expression is regulated by three combinational cis-regulatory modules which are activated by EGFR pathway, Ap-Vg auto-regulatory and epigenetic mechanisms. Here, we found that the Tbx family transcription factor Optomotor-blind (Omb) restricted ap expression in the ventral compartment. Loss of omb induced autonomous initiation of ap expression in the middle third instar larvae in the ventral compartment. Oppositely, over-activation of omb inhibited ap in the medial pouch. All three enhancers apE, apDV and apP were upregulated in omb null mutants, indicating a combinational regulation of ap modulators. However, Omb affected ap expression neither by directly regulating EGFR signaling, nor via Vg regulation. Therefore, a genetic screen of epigenetic regulators, including the Trithorax group (TrxG) and Polycomb group (PcG) genes was performed. We found that knocking down the TrxG gene kohtalo (kto), domino (dom) or expressing the PcG gene grainy head (grh), the ectopic ap in omb mutants was repressed. The inhibition of apDV by kto knockdown and grh activation could contribute to ap repression. Moreover, Omb and the EGFR pathway are genetically parallel in ap regulation in the ventral compartment. Collectively, Omb is a repressive signal for ap expression in the ventral compartment, which requires TrxG and PcG genes.

Transgenerational inheritance of wing development defects in Drosophila melanogaster induced by cadmium

The transgenerational inheritance of phenotype induced by environmental factors is a new focus in epigenetic research. In this study, Drosophila melanogaster (F₀) was cultured in the medium containing cadmium (Cd, 4.5 mg/kg) from eggs to adults, and offspring (F₁-F₄) were continuously kept in standard medium (without cadmium). The phenotype analysis showed that cadmium induced developmental defects on wings and apoptosis in the wing disc cells of Drosophila (F₀). The wing defects were transmitted for at least four generations even without Cd afterwards. And the effect on the mRNA expression of wing development related genes (shg, omb, F-actin, Mekk1) can be maintained for at least two or three generations. More importantly, under cadmium stress, the post-translational modification (PTM) on the histones H3K4me3 in the third instar larvae and ovaries or testes of adult flies increased significantly, while the levels of H3K9me3 and H3K27me3 decreased significantly. The expression of histone methylation related genes (dSet-1, ash1, Lsd1) increased significantly and these changes can be transmitted to offspring from one or two generations in ovaries or testes. These results suggest that the phenotypic defects of wings caused by cadmium can be inherited to the offspring, and this transgenerational inheritance effect may be related to the epigenetic regulation of histone methylation. Therefore, the adaptability of offspring should be considered when evaluating the toxicity and environmental risk of cadmium.

Non-additive gene interactions underpin molecular and phenotypic responses in honey bee larvae exposed to imidacloprid and thymol

Understanding the cumulative risk of chemical mixtures at environmentally realistic concentrations is a key challenge in honey bee ecotoxicology. Ecotoxicogenomics, including transcriptomics, measures responses in individual organisms at the molecular level which can provide insights into the mechanisms underlying phenotypic responses induced by one or more stressors and link impacts on individuals to populations. Here, fifth instar honey bee larvae were sampled from a previously reported field experiment exploring the phenotypic impacts of environmentally realistic chronic exposures of the pesticide imidacloprid (5 μg.kg^-1 for six weeks) and the acaricide thymol (250 g.kg^-1 applied via Apiguard gel in-hive for four weeks), both separately and in combination. RNA-seq was used to discover individual and interactive chemical effects on larval gene expression and to uncover molecular mechanisms linked to reported adult and colony phenotypes. The separate and combined treatments had distinct gene expression profiles which represented differentially affected signaling and metabolic pathways. The molecular signature of the mixture was characterised by additive interactions in canonical stress responses associated with oxidative stress and detoxification, and non-additive interactions in secondary responses including developmental, neurological, and immune pathways. Novel emergent impacts on eye development genes correlated with long-term defects in visual learning performance as adults. This is consistent with these chemicals working through independent modes of action that combine to impact common downstream pathways, and highlights the importance of establishing mechanistic links between molecular and phenotypic responses when predicting effects of chemical mixtures on ecologically relevant population outcomes.

The T-box gene optomotor-blind organizes proximodistal leg patterning in the beetle Tribolium castaneum by repressing dorsal Dpp pathway activity

Leg axis formation in Drosophila is organized by Wingless (Wg) and Decapentaplegic (Dpp) that control a number of downstream factors to pattern the dorsoventral (DV) and proximodistal (PD) axis. The T-box genes are important downstream factors mainly involved in dorsoventral leg axis formation. The ventral side is specified by H15 and midline, whereas optomotor-blind (omb) and Dorsocross (Doc1) are factors to specify dorsal cell fates. We show here that omb also organizes PD leg axis patterning in the beetle Tribolium castaneum. In the legs, Tc-omb is expressed along the dorsal side and represses ventral factors like wg and H15. Intriguingly, removing Tc-omb function leads to the activation of the Dpp pathway along the dorsal side of the legs, thus mimicking normal dpp expression in Drosophila. Dpp activity along the dorsal side leads to altered expression of proximal-distal patterning genes such as Distal-less (Dll) and dachshund (dac). Our results indicate a cell-autonomous activation of Dll and repression of dac by dpp. These findings are compatible with the cross-regulatory "cascade model" of proximal-distal leg imaginal disc patterning of Drosophila.

A chelicerate Wnt gene expression atlas: novel insights into the complexity of arthropod Wnt-patterning

The Wnt genes represent a large family of secreted glycoprotein ligands that date back to early animal evolution. Multiple duplication events generated a set of 13 Wnt families of which 12 are preserved in protostomes. Embryonic Wnt expression patterns (Wnt-patterning) are complex, representing the plentitude of functions these genes play during development. Here, we comprehensively investigated the embryonic expression patterns of Wnt genes from three species of spiders covering both main groups of true spiders, Haplogynae and Entelegynae, a mygalomorph species (tarantula), as well as a distantly related chelicerate outgroup species, the harvestman Phalangium opilio. All spiders possess the same ten classes of Wnt genes, but retained partially different sets of duplicated Wnt genes after whole genome duplication, some of which representing impressive examples of sub- and neo-functionalization. The harvestman, however, possesses a more complete set of 11 Wnt genes but with no duplicates. Our comprehensive data-analysis suggests a high degree of complexity and evolutionary flexibility of Wnt-patterning likely providing a firm network of mutational protection. We discuss the new data on Wnt gene expression in terms of their potential function in segmentation, posterior elongation, and appendage development and critically review previous research on these topics. We conclude that earlier research may have suffered from the absence of comprehensive gene expression data leading to partial misconceptions about the roles of Wnt genes in development and evolution.

Characterization of Long Non-coding RNAs Modified by m6A RNA Methylation in Skeletal Myogenesis

Proper development of mammalian skeletal muscle relies on precise gene expression regulation. Our previous studies revealed that muscle development is regulated by both mRNA and long non-coding RNAs (lncRNAs). Accumulating evidence has demonstrated that N⁶-methyladenosine (m⁶A) plays important roles in various biological processes, making it essential to profile m⁶A modification on a transcriptome-wide scale in developing muscle. Patterns of m⁶A methylation in lncRNAs in developing muscle have not been uncovered. Here, we reveal differentially expressed lncRNAs and report temporal m⁶A methylation patterns in lncRNAs expressed in mouse myoblasts and myotubes by RNA-seq and methylated RNA immunoprecipitation (MeRIP) sequencing. Many lncRNAs exhibit temporal differential expression, and m⁶A-lncRNAs harbor the consensus m⁶A motif “DRACH” along lncRNA transcripts. Interestingly, we found that m⁶A methylation levels of lncRNAs are positively correlated with the transcript abundance of lncRNAs. Overexpression or knockdown of m⁶A methyltransferase METTL3 alters the expression levels of these lncRNAs. Furthermore, we highlight that the function of m⁶A genic lncRNAs might correlate to their nearby mRNAs. Our work reveals a fundamental expression reference of m⁶A-mediated epitranscriptomic modifications in lncRNAs that are temporally expressed in developing muscle.

TBX2/3 is required for regeneration of dorsal-ventral and medial-lateral polarity in planarians

The molecular mechanisms responsible for axis establishment during non-embryonic processes remain elusive. The planarian flatworm is an ideal model organism to study body axis polarization and patterning in vivo. Here, we identified a homolog of the TBX2/3 in the planarian Dugesia japonica. RNA interference (RNAi) knockdown of TBX2/3 results in the ectopic formation of protrusions in the midline of the dorsal surface which shows an abnormal expression of midline and ventral cell markers. Additionally, the TBX2/3 RNAi animals also show the duplication of expression of the boundary marker at the lateral edge. Furthermore, TBX2/3 is expressed in muscle cells and co-expressed with bmp4. Inhibition of bone morphogenetic protein (BMP) signaling reduces the expression of TBX2/3 at the midline. These results suggest that TBX2/3 RNAi results in phenotypic characters caused by inhibition of the BMP signal, indicating that TBX2/3 is required for DV and ML patterning, and might be a downstream gene of BMP signaling.

T-box transcription factors Dorsocross and optomotor-blind control Drosophila leg patterning in a functionally redundant manner

The T-box genes are essential transcription factors during limb development. In Drosophila, Dorsocross (Doc) and optomotor-blind (omb), members of the Tbx2 and Tbx6 families, are best studied in the Drosophila wing development. Despite prominently expressed in leg discs, the specific function of these genes in leg growth is still not revealed. Here we demonstrated that Doc and omb regulated the morphogenesis of leg intermediate regions in a functionally redundant manner. Loss of Doc or omb individually did not result in any developmental defects of the legs, but loss of both genes induced significant defects in femur and proximal tibia of the adult legs. These genes located in the dorsal domain, where the Doc region expanded and cross-overlapped with the omb region corresponding to the presumptive leg intermediate region. We detected that the normal epithelial folds in the leg discs were disrupted along with dorsal repression of cell proliferation and activation of cell apoptosis when Doc and omb were both reduced. Furthermore, the dorsal expression of dachshund (dac), a canonical leg developmental gene specifying the leg intermediate region, was maintained by Doc and omb. Meanwhile, the Notch pathway was compromised in the dorsal domain when these genes were reduced, which might contribute to the joint defect of the adult leg intermediate regions. Our study provides cytological and genetic evidence for understanding the redundant function of Doc and omb in leg morphogenesis.

Genome streamlining in a minute herbivore that manipulates its host plant

The tomato russet mite, Aculops lycopersici, is among the smallest animals on earth. It is a worldwide pest on tomato and can potently suppress the host's natural resistance. We sequenced its genome, the first of an eriophyoid, and explored whether there are genomic features associated with the mite's minute size and lifestyle. At only 32.5 Mb, the genome is the smallest yet reported for any arthropod and, reminiscent of microbial eukaryotes, exceptionally streamlined. It has few transposable elements, tiny intergenic regions, and is remarkably intron-poor, as more than 80% of coding genes are intronless. Furthermore, in accordance with ecological specialization theory, this defense-suppressing herbivore has extremely reduced environmental response gene families such as those involved in chemoreception and detoxification. Other losses associate with this species' highly derived body plan. Our findings accelerate the understanding of evolutionary forces underpinning metazoan life at the limits of small physical and genome size.

Transcriptome analysis of Aedes aegypti Aag2 cells in response to dengue virus-2 infection

Background

Dengue virus (DENV) is a flavivirus transmitted by mosquitoes that is prevalent in tropical and subtropical countries and has four serotypes (DENV1-4). Aedes aegypti, as the main transmission vector of DENV, exhibits strong infectivity and transmission. With the aim of obtaining a better understanding of the Ae. aegypti-DENV interaction, the transcriptome changes in DENV-2-infected Aag2 cells were studied to describe the immune responses of mosquitoes using the Ae. aegypti Aag2 cell line as a model.

Methods

RNAseq technology was used to sequence the transcripts of the Ae. aegypti Aag2 cell line before and after infection with DENV-2. A bioinformatics analysis was then performed to assess the biological functions of the differentially expressed genes, and the sequencing data were verified by quantitative reverse transcription-polymerase chain reaction (qRT-PCR).

Results

The transcriptome analysis generated 8866 unigenes that were found in both groups, 225 unigenes that were only found in the infection group, and 683 unigenes that only existed in the control group. A total of 1199 differentially expressed genes, including 1014 upregulated and 185 downregulated genes, were identified. The bioinformatics analysis showed that the differentially expressed genes were mainly involved in the longevity regulating pathway, circadian rhythm, DNA replication, and peroxisome, purine, pyrimidine, and drug metabolism. The qRT-PCR verification results showed the same trend, which confirmed that the expression of the differentially expressed genes had changed, and that the transcriptome sequencing data were reliable.

Conclusions

This study investigated the changes in the transcriptome levels in the DENV-2-infected Ae. aegypti Aag2 cell line, which provides a faster and effective method for discovering genes related to Ae. aegypti pathogen susceptibility. The findings provide basic data and directions for further research on the complex mechanism underlying host-pathogen interactions.

Two sets of candidate crustacean wing homologues and their implication for the origin of insect wings

The origin of insect wings is a biological mystery that has fascinated scientists for centuries. Identification of tissues homologous to insect wings from lineages outside of Insecta will provide pivotal information to resolve this conundrum. Here, through expression and clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9) functional analyses in Parhyale, we show that a gene network similar to the insect wing gene network (preWGN) operates both in the crustacean terga and in the proximal leg segments, suggesting that the evolution of a preWGN precedes the emergence of insect wings, and that from an evo-devo perspective, both of these tissues qualify as potential crustacean wing homologues. Combining these results with recent wing origin studies in insects, we discuss the possibility that both tissues are crustacean wing homologues, which supports a dual evolutionary origin of insect wings (that is, novelty through a merger of two distinct tissues). These outcomes have a crucial impact on the course of the intellectual battle between the two historically competing wing origin hypotheses.

BmBlimp-1 gene encoding a C2H2 zinc finger protein is required for wing development in the silkworm Bombyx mori

Cys2-His2 zinc finger (C2H2-ZF) proteins represent the most common class of transcription factors. These factors have great potential for the management of developmental progression by regulating the specific spatiotemporal expression of genes. In this study, we cloned one C2H2-ZF protein gene of Bombyx mori, BGIBMGA000319, that is orthologous to B-lymphocyte-induced maturation protein-1 (Blimp-1); we thus named it as Bombyx mori Blimp-1 (BmBlimp-1). In the silkworm, the BmBlimp-1 gene is specifically upregulated during day 2 of the pupal to adult stage and is highly expressed in wing discs on day 3 of the pupa. Knockdown of its expression level in the pupal stage results in a crumpled-winged silkworm moth. Using the predicted DNA-binding sequences of BmBlimp-1 to search the silkworm genome to screen target genes of BmBlimp-1, 7049 genes were identified to have at least one binding site of BmBlimp-1 on their 1 kb upstream and downstream genome regions. Comparisons of those genes with a reported pupal wing disc transcriptome data resulted in 4065 overlapping genes being retrieved. GO enrichment analysis of the overlapping genes showed that most of the genes were enriched in the binding term. Combining functional annotation and real-time quantitative PCR, 15 genes were identified as the candidate target genes of BmBlimp-1, including several wing cuticular protein genes, chitin synthase A, and wing disc development genes, such as Wnt1, cubitus interruptus (ci) and engrailed (en). Moreover, the amino acid sequence of the zinc finger motif of Blimp-1 gene was highly conserved among the 15 insect species. We propose that BmBlimp-1 is an important regulatory factor in silkworm wing development.

FoxB, a new and highly conserved key factor in arthropod dorsal-ventral (DV) limb patterning

Forkhead box (Fox) transcription factors evolved early in animal evolution and represent important components of conserved gene regulatory networks (GRNs) during animal development. Most of the researches concerning Fox genes, however, are on vertebrates and only a relatively low number of studies investigate Fox gene function in invertebrates. In addition to this shortcoming, the focus of attention is often restricted to a few well-characterized Fox genes such as FoxA (forkhead), FoxC (crocodile) and FoxQ2. Although arthropods represent the largest and most diverse animal group, most other Fox genes have not been investigated in detail, not even in the arthropod model species Drosophila melanogaster. In a general gene expression pattern screen for panarthropod Fox genes including the red flour beetle Tribolium castaneum, the pill millipede Glomeris marginata, the common house spider Parasteatoda tepidariorum, and the velvet worm Euperipatoides kanangrensis, we identified a Fox gene with a highly conserved expression pattern along the ventral ectoderm of arthropod and onychophoran limbs. Functional investigation of FoxB in Parasteatoda reveals a hitherto unrecognized important function of FoxB upstream of wingless (wg) and decapentaplegic (dpp) in the GRN orchestrating dorsal–ventral limb patterning.

Transcriptional repression of the ectodomain sheddase ADAM10 by TBX2 and potential implication for Alzheimer's disease

BACKGROUND

The ADAM10-mediated cleavage of transmembrane proteins regulates cellular processes such as proliferation or migration. Substrate cleavage by ADAM10 has also been implicated in pathological situations such as cancer or Morbus Alzheimer. Therefore, identifying endogenous molecules, which modulate the amount and consequently the activity of ADAM10, might contribute to a deeper understanding of the enzyme's role in both, physiology and pathology.

METHOD

To elucidate the underlying cellular mechanism of the TBX2-mediated repression of ADAM10 gene expression, we performed overexpression, RNAi-mediated knockdown and pharmacological inhibition studies in the human neuroblastoma cell line SH-SY5Y. Expression analysis was conducted by e.g. real-time RT-PCR or western blot techniques. To identify the binding region of TBX2 within the ADAM10 promoter, we used luciferase reporter assay on deletion constructs and EMSA/WEMSA experiments. In addition, we analyzed a TBX2 loss-of-function Drosophila model regarding the expression of ADAM10 orthologs by qPCR. Furthermore, we quantified the mRNA level of TBX2 in post-mortem brain tissue of AD patients.

RESULTS

Here, we report TBX2 as a transcriptional repressor of ADAM10 gene expression: both, the DNA-binding domain and the repression domain of TBX2 were necessary to effect transcriptional repression of ADAM10 in neuronal SH-SY5Y cells. This regulatory mechanism required HDAC1 as a co-factor of TBX2. Transcriptional repression was mediated by two functional TBX2 binding sites within the core promoter sequence (- 315 to - 286 bp). Analysis of a TBX2 loss-of-function Drosophila model revealed that kuzbanian and kuzbanian-like, orthologs of ADAM10, were derepressed compared to wild type. Vice versa, analysis of cortical brain samples of AD-patients, which showed reduced ADAM10 mRNA levels, revealed a 2.5-fold elevation of TBX2, while TBX3 and TBX21 levels were not affected.

CONCLUSION

Our results characterize TBX2 as a repressor of ADAM10 gene expression and suggest that this regulatory interaction is conserved across tissues and species.

bantam microRNA is a negative regulator of the Drosophila decapentaplegic pathway

Decapentaplegic (Dpp), the Drosophila homolog of the vertebrate bone morphogenetic protein (BMP2/4), is crucial for patterning and growth in many developmental contexts. The Dpp pathway is regulated at many different levels to exquisitely control its activity. We show that bantam (ban), a microRNA, modulates Dpp signaling activity. Over expression of ban decreases phosphorylated Mothers against decapentaplegic (Mad) levels and negatively affects Dpp pathway transcriptional target genes, while null mutant clones of ban upregulate the pathway. We provide evidence that dpp upregulates ban in the wing imaginal disc, and attenuation of Dpp signaling results in a reduction of ban expression, showing that they function in a feedback loop. Furthermore, we show that this feedback loop is important for maintaining anterior-posterior compartment boundary stability in the wing disc through regulation of optomotor blind (omb), a known target of the pathway. Our results support a model that ban functions with dpp in a negative feedback loop.

Genetics of clubfoot; recent progress and future perspectives

Clubfoot or talipes equinovarus (TEV) is an inborn three-dimensional deformity of leg, ankle and foot. It results from structural defects of several tissues of foot and lower leg leading to abnormal positioning of foot and ankle joints. TEV can lead to long-lasting functional disability, malformation and discomfort if left untreated. Substantial progress has been achieved in the management and diagnosis of limb defects; however, not much is known about the molecular players and signalling pathways underlying TEV disorder. The homeostasis and development of the limb depends on the complex interactions between the lateral plate mesoderm cells and outer ectoderm. These complex interactions include HOX signalling and PITX1-TBX4 pathways. The susceptibility to develop TEV is determined by a number of environmental and genetic factors, although the nature and level of interplay between them remains unclear. Familial occurrence and inter and intra phenotypic variability of TEV is well documented. Variants in genes that code for contractile proteins of skeletal myofibers might play a role in the aetiology of TEV but, to date, no strong candidate genes conferring increased risk have emerged, although variants in TBX4, PITX1, HOXA, HOXC and HOXD clusters genes, NAT2 and others have been shown to be associated with TEV. The mechanisms by which variants in these genes confer risk and the nature of the physical and genetic interaction between them remains to be determined. Elucidation of genetic players and cellular pathways underlying TEV will certainly increase our understanding of the pathophysiology of this deformity.