The Homeobox Genes of Caenorhabditis elegans and Insights into Their Spatio-Temporal Expression Dynamics during Embryogenesis

Sex-specific developmental gene expression atlas unveils dimorphic gene networks in C. elegans

Sex-specific traits and behaviors emerge during development by the acquisition of unique properties in the nervous system of each sex. However, the genetic events responsible for introducing these sex-specific features remain poorly understood. In this study, we create a comprehensive gene expression atlas of pure populations of hermaphrodites and males of the nematode Caenorhabditis elegans across development. We discover numerous differentially expressed genes, including neuronal gene families like transcription factors, neuropeptides, and G protein-coupled receptors. We identify INS-39, an insulin-like peptide, as a prominent male-biased gene expressed specifically in ciliated sensory neurons. We show that INS-39 serves as an early-stage male marker, facilitating the effective isolation of males in high-throughput experiments. Through complex and sex-specific regulation, ins-39 plays pleiotropic sexually dimorphic roles in various behaviors, while also playing a shared, dimorphic role in early life stress. This study offers a comparative sexual and developmental gene expression database for C. elegans. Furthermore, it highlights conserved genes that may underlie the sexually dimorphic manifestation of different human diseases.

Single-cell transcriptomics of the human parasite Schistosoma mansoni first intra-molluscan stage reveals tentative tegumental and stem-cell regulators

Schistosomiasis is a major Neglected Tropical Disease, caused by the infection with blood flukes in the genus Schistosoma. To complete the life cycle, the parasite undergoes asexual and sexual reproduction within an intermediate snail host and a definitive mammalian host, respectively. The intra-molluscan phase provides a critical amplification step that ensures a successful transmission. However, the cellular and molecular mechanisms underlying the development of the intra-molluscan stages remain poorly understood. Here, single cell suspensions from S. mansoni mother sporocysts were produced and sequenced using the droplet-based 10X Genomics Chromium platform. Six cell clusters comprising two tegument, muscle, neuron, parenchyma and stem/germinal cell clusters were identified and validated by in situ hybridisation. Gene Ontology term analysis predicted key biological processes for each of the clusters, including three stem/germinal sub-clusters. Furthermore, putative transcription factors predicted for stem/germinal and tegument clusters may play key roles during parasite development and interaction with the intermediate host.

The homeodomain transcription factor CEH-37 regulates PMK-1/p38 MAPK pathway to protect against intestinal infection via the phosphatase VHP-1

Increasing evidence indicate that the expression of defense genes at the right place and the right time are regulated by host-defense transcription factors. However, the precise mechanisms of this regulation are not well understood. Homeodomain transcription factors, encoded by homeobox genes, play crucial role for the development of multicellular eukaryotes. In this study, we demonstrated that homeodomain transcription factor CEH-37 (known as OTX2 in mammals) was a key transcription factor for host defense in Caenorhabditis elegans. Meanwhile, CEH-37 acted in the intestine to protect C. elegans against pathogen infection. We further showed that the homeodomain transcription factor CEH-37 positively regulated PMK-1/ p38 MAPK activity to promote the intestinal immunity via suppression phosphatase VHP-1. Furthermore, we demonstrated that this function was conserved, because the human homeodomain transcription factor OTX2 also exhibited protective function in lung epithelial cells during Pseudomonas aeruginosa infection. Thus, our work reveal that CEH-37/OTX2 is a evolutionarily conserved transcription factor for defense against pathogen infection. The finding provides a model in which CEH-37 decreases VHP-1 phosphatase activity, allowing increased stimulation of PMK-1/p38 MAPK phosphorylation cascade in the intestine for pathogen resistance.

An insight into the functional genomics and species classification of Eudiplozoon nipponicum (Monogenea, Diplozoidae), a haematophagous parasite of the common carp Cyprinus carpio

BACKGROUND

Monogenea (Platyhelminthes, Neodermata) are the most species-rich class within the Neodermata superclass of primarily fish parasites. Despite their economic and ecological importance, monogenean research tends to focus on their morphological, phylogenetic, and population characteristics, while comprehensive omics analyses aimed at describing functionally important molecules are few and far between. We present a molecular characterisation of monogenean representative Eudiplozoon nipponicum, an obligate haematophagous parasite infecting the gills of the common carp. We report its nuclear and mitochondrial genomes, present a functional annotation of protein molecules relevant to the molecular and biochemical aspect of physiological processes involved in interactions with the fish hosts, and re-examinate the taxonomic position of Eudiplozoon species within the Diplozoidae family.

RESULTS

We have generated 50.81 Gbp of raw sequencing data (Illumina and Oxford Nanopore reads), bioinformatically processed, and de novo assembled them into a genome draft 0.94 Gbp long, consisting of 21,044 contigs (N50 = 87 kbp). The final assembly represents 57% of the estimated total genome size (~ 1.64 Gbp), whereby repetitive and low-complexity regions account for ~ 64% of the assembled length. In total, 36,626 predicted genes encode 33,031 proteins and homology-based annotation of protein-coding genes (PCGs) and proteins characterises 14,785 (44.76%) molecules. We have detected significant representation of functional proteins and known molecular functions. The numbers of peptidases and inhibitors (579 proteins), characterised GO terms (16,016 unique assigned GO terms), and identified KEGG Orthology (4,315 proteins) acting in 378 KEGG pathways demonstrate the variety of mechanisms by which the parasite interacts with hosts on a macromolecular level (immunomodulation, feeding, and development). Comparison between the newly assembled E. nipponicum mitochondrial genome (length of 17,038 bp) and other diplozoid monogeneans confirms the existence of two distinct Eudiplozoon species infecting different fish hosts: Cyprinus carpio and Carassius spp.

CONCLUSIONS

Although the amount of sequencing data and characterised molecules of monogenean parasites has recently increased, a better insight into their molecular biology is needed. The E. nipponicum nuclear genome presented here, currently the largest described genome of any monogenean parasite, represents a milestone in the study of monogeneans and their molecules but further omics research is needed to understand these parasites' biological nature.

A Compilation of the Diverse miRNA Functions in Caenorhabditis elegans and Drosophila melanogaster Development

MicroRNAs are critical regulators of post-transcriptional gene expression in a wide range of taxa, including invertebrates, mammals, and plants. Since their discovery in the nematode, Caenorhabditis elegans, miRNA research has exploded, and they are being identified in almost every facet of development. Invertebrate model organisms, particularly C. elegans, and Drosophila melanogaster, are ideal systems for studying miRNA function, and the roles of many miRNAs are known in these animals. In this review, we compiled the functions of many of the miRNAs that are involved in the development of these invertebrate model species. We examine how gene regulation by miRNAs shapes both embryonic and larval development and show that, although many different aspects of development are regulated, several trends are apparent in the nature of their regulation.

Genome-Wide Analysis of WOX Multigene Family in Sunflower (Helianthus annuus L.)

The WUSCHEL-related homeobox (WOX) is a family of specific transcription factors involved in plant development and response to stress, characterized by the presence of a homeodomain. This study represents the first comprehensive characterization of the WOX family in a member of the Asteraceae family, the sunflower (H. annuus L.). Overall, we identified 18 putative HaWOX genes divided by phylogenetic analysis in three major clades (i.e., ancient, intermediate, and WUS). These genes showed conserved structural and functional motifs. Moreover, HaWOX has homogeneously distributed on H. annuus chromosomes. In particular, 10 genes originated after whole segment duplication events, underpinning a possible evolution of this family along with the sunflower genome. In addition, gene expression analysis evidenced a specific pattern of regulation of the putative 18 HaWOX during embryo growth and in ovule and inflorescence meristem differentiation, suggesting a pivotal role for this multigenic family in sunflower development. The results obtained in this work improved the understanding of the WOX multigenic family, providing a resource for future study on functional analysis in an economically valuable species such as sunflower.

Delineating mechanisms and design principles of Caenorhabditis elegans embryogenesis using in toto high-resolution imaging data and computational modeling

The nematode (roundworm) Caenorhabditis elegans is one of the most popular animal models for the study of developmental biology, as its invariant development and transparent body enable in toto cellular-resolution fluorescence microscopy imaging of developmental processes at 1-min intervals. This has led to the development of various computational tools for the systematic and automated analysis of imaging data to delineate the molecular and cellular processes throughout the embryogenesis of C. elegans, such as those associated with cell lineage, cell migration, cell morphology, and gene activity. In this review, we first introduce C. elegans embryogenesis and the development of techniques for tracking cell lineage and reconstructing cell morphology during this process. We then contrast the developmental modes of C. elegans and the customized technologies used for studying them with the ones of other animal models, highlighting its advantage for studying embryogenesis with exceptional spatial and temporal resolution. This is followed by an examination of the physical models that have been devised—based on accurate determinations of developmental processes afforded by analyses of imaging data—to interpret the early embryonic development of C. elegans from subcellular to intercellular levels of multiple cells, which focus on two key processes: cell polarization and morphogenesis. We subsequently discuss how quantitative data-based theoretical modeling has improved our understanding of the mechanisms of C. elegans embryogenesis. We conclude by summarizing the challenges associated with the acquisition of C. elegans embryogenesis data, the construction of algorithms to analyze them, and the theoretical interpretation.

Mutations in Growth-Related Genes Induced by EMS Treatment in Scallops

Background: The goal of genetic breeding is to select variants with mutations that are related to expected traits, such as fast growth. Artificial induction has been widely used to obtain strains with more mutations for further selection. Ethylmethylsulfone (EMS) is one of the most commonly used chemical mutagens in plant and microorganism breeding. However, the application of EMS mutagenesis in shellfish has not been reported. The aim of this study is to evaluate the potential use of EMS as a mutagen in scallop breeding, especially in characterization of mutations in growth-related genes. Results: Our results indicated that hatching of about 50% of fertilized eggs was blocked by treatment with 20 mM EMS for 3 h and the resulted larvae developed normally into adult stages. We then evaluated the mutagenic effects of EMS by sequencing the genomes of 4 adult scallops from the control group and 12 from the treatment group at 8 months after fertilization. On average, after removing shared types of mutations, there were 1,151,380 ± 258,188 SNPs (Single Nucleotide Polymorphisms) and 229,256 ± 51,714 InDels (insertion-deletion) in each animal in the EMS treatment group, while there were only134841 ± 10,115 SNPs and 42,605 ± 5,136 InDels in the control group. The average mutation rate in the genome of the EMS treatment group (0.0137 ± 0.0013%) was about 9 times that of the control group (0.0015 ± 0.0002%). GO (Gene Ontology) annotation and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analyses revealed that mutations induced by EMS occurred evenly in most biological processes, cellular components and functions, as well in most pathways. However, significant lower percentage of mutations were found in the exonic region, in non-synonymous or Stopgain/Stoploss SNPs and in coding domains, suggesting apparent DNA repair or selection during grow-out stage. Analyses of the growth-related genes with mutations indicated that mutations in MFS (Major Facilitator Superfamily) and Tubulin were only found in the large-sized group (Five largest scallops: Treated-1, Treated-2, Treated-3, Treated-4, and Treated-5) and Homeobox and Socs (Suppressor of cytokine signaling) only in the small group (Two smallest scallops: Treated-11 and Treated-12). These results suggested that these genes may be involved in the regulation of growth in these animals, although further verification is certainly warranted. Conclusion: Treatment of fertilized eggs with 20 mM EMS for 3 h induced 9 times more mutations in scallop genomes. We found that mutations in MFS and Tubulin may be related to fast growth in the large-sized group and those mutations in Homeobox and SOCs may be involved in the slow growth in the small-sized scallops. EMS can be used to accelerate selection of economically important traits in molluscs.

Emergence of distinct syntenic density regimes is associated with early metazoan genomic transitions

Background

Animal genomes are strikingly conserved in terms of local gene order (microsynteny). While some of these microsyntenies have been shown to be coregulated or to form gene regulatory blocks, the diversity of their genomic and regulatory properties across the metazoan tree of life remains largely unknown.

Results

Our comparative analyses of 49 animal genomes reveal that the largest gains of synteny occurred in the last common ancestor of bilaterians and cnidarians and in that of bilaterians. Depending on their node of emergence, we further show that novel syntenic blocks are characterized by distinct functional compositions (Gene Ontology terms enrichment) and gene density properties, such as high, average and low gene density regimes. This is particularly pronounced among bilaterian novel microsyntenies, most of which fall into high gene density regime associated with higher gene coexpression levels. Conversely, a majority of vertebrate novel microsyntenies display a low gene density regime associated with lower gene coexpression levels.

Conclusions

Our study provides first evidence for evolutionary transitions between different modes of microsyntenic block regulation that coincide with key events of metazoan evolution. Moreover, the microsyntenic profiling strategy and interactive online application (Syntenic Density Browser, available at: http://synteny.csb.univie.ac.at/) we present here can be used to explore regulatory properties of microsyntenic blocks and predict their coexpression in a wide-range of animal genomes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08304-2.

Gene cluster from plant to microbes: Their role in genome architecture, organism's development, specialized metabolism and drug discovery

Plants and microbes fulfil our daily requirements through different high-value chemicals, e.g., nutraceuticals, pharmaceuticals, cosmetics, and through varieties of fruits, crops, vegetables, and many more. Utmost care would therefore be taken for growth, development and sustainability of these important crops and medicinal plants and microbes. Homeobox genes and HOX clusters and their recently characterized expanded family members, including newly discovered homeobox, WOX gene from medicinal herb, Panax ginseng, significantly contributes in the growth and development of these organisms. On the other hand, secondary metabolites produced through secondary metabolism of plants and microbes are used as organisms defense as well as drugs/drug-like molecules for humans. Both the developmental HOX cluster and the biosynthetic gene-cluster (BGC) for secondary metabolites are organised in organisms genome. Genome mining and genomewide analysis of these clusters will definitely identify and characterize many more important molecules from unexplored plants and microbes and underexplored human microbiota and the evolution studies of these clusters will indicate their source of origin. Although genomics revolution now continues at a pace, till date only few hundred plant genome sequences are available. However, next-generation sequencing (NGS) technology now in market and may be applied even for plants with recalcitrant genomes, eventually may discover genomic potential towards production of secondary metabolites of diverse plants and micro-organisms present in the environment and microbiota. Additionally, the development of tools for genome mining e.g., antiSMASH, plantiSMASH, and more and more computational approaches that predicts hundreds of secondary metabolite BGCs will be discussed.

A genome-wide identification and analysis of the homeobox genes in the brown planthopper, Nilaparvata lugens (Hemiptera: Delphacidae)

The homeobox family is a large and diverse superclass of genes, many of which act as transcription factors that play important roles in tissue differentiation and embryogenesis in animals. The brown planthopper (BPH), Nilaparvata lugens, is the most destructive pest of rice in Asia, and high fecundity contributes significantly to its ecological success in natural and agricultural habits. Here, we identified 94 homeobox genes in BPH, which could be divided into 75 gene families and 9 classes. This number is comparable to the number of homeobox genes found in the honeybee Apis mellifera, but is slightly less than in Drosophila or the red flour beetle Tribolium castaneum. A spatio-temporal analysis indicated that most BPH homeobox genes were expressed in a development and tissue-specific manner, of which 21 genes were highly expressed in ovaries. RNA interference (RNAi)-mediated functional assay showed that 22 homeobox genes were important for nymph development and the nymph to adult transition, whereas 67 genes were dispensable during this process. Fecundity assay showed that knockdown of 13 ovary-biased genes (zfh1, schlank, abd-A, Lim3_2, Lmxb, Prop, ap_1, Not, lab, Hmx, vis, Pknox, and C15) led to the reproductive defect. This is the first comprehensive investigation into homeobox genes in a hemipteran insect and thus helps us to understand the functional significance of homeobox genes in insect reproduction.

Homeobox genes and the specification of neuronal identity

The enormous diversity of cell types that characterizes any animal nervous system is defined by neuron-type-specific gene batteries that endow cells with distinct anatomical and functional properties. To understand how such cellular diversity is genetically specified, one needs to understand the gene regulatory programmes that control the expression of cell-type-specific gene batteries. The small nervous system of the nematode Caenorhabditis elegans has been comprehensively mapped at the cellular and molecular levels, which has enabled extensive, nervous system-wide explorations into whether there are common underlying mechanisms that specify neuronal cell-type diversity. One principle that emerged from these studies is that transcription factors termed 'terminal selectors' coordinate the expression of individual members of neuron-type-specific gene batteries, thereby assigning unique identities to individual neuron types. Systematic mutant analyses and recent nervous system-wide expression analyses have revealed that one transcription factor family, the homeobox gene family, is broadly used throughout the entire C. elegans nervous system to specify neuronal identity as terminal selectors. I propose that the preponderance of homeobox genes in neuronal identity control is a reflection of an evolutionary trajectory in which an ancestral neuron type was specified by one or more ancestral homeobox genes, and that this functional linkage then duplicated and diversified to generate distinct cell types in an evolving nervous system.

Meta-analysis of Caenorhabditis elegans single-cell developmental data reveals multi-frequency oscillation in gene activation

MOTIVATION

The advent of in vivo automated techniques for single-cell lineaging, sequencing and analysis of gene expression has begun to dramatically increase our understanding of organismal development. We applied novel meta-analysis and visualization techniques to the EPIC single-cell-resolution developmental gene expression dataset for Caenorhabditis elegans from Bao, Murray, Waterston et al. to gain insights into regulatory mechanisms governing the timing of development.

RESULTS

Our meta-analysis of the EPIC dataset revealed that a simple linear combination of the expression levels of the developmental genes is strongly correlated with the developmental age of the organism, irrespective of the cell division rate of different cell lineages. We uncovered a pattern of collective sinusoidal oscillation in gene activation, in multiple dominant frequencies and in multiple orthogonal axes of gene expression, pointing to the existence of a coordinated, multi-frequency global timing mechanism. We developed a novel method based on Fisher's Discriminant Analysis to identify gene expression weightings that maximally separate traits of interest, and found that remarkably, simple linear gene expression weightings are capable of producing sinusoidal oscillations of any frequency and phase, adding to the growing body of evidence that oscillatory mechanisms likely play an important role in the timing of development. We cross-linked EPIC with gene ontology and anatomy ontology terms, employing Fisher's Discriminant Analysis methods to identify previously unknown positive and negative genetic contributions to developmental processes and cell phenotypes. This meta-analysis demonstrates new evidence for direct linear and/or sinusoidal mechanisms regulating the timing of development. We uncovered a number of previously unknown positive and negative correlations between developmental genes and developmental processes or cell phenotypes. Our results highlight both the continued relevance of the EPIC technique, and the value of meta-analysis of previously published results. The presented analysis and visualization techniques are broadly applicable across developmental and systems biology.

AVAILABILITY AND IMPLEMENTATION

Analysis software available upon request.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A Series of Tubes: The C. elegans Excretory Canal Cell as a Model for Tubule Development

Formation and regulation of properly sized epithelial tubes is essential for multicellular life. The excretory canal cell of C. elegans provides a powerful model for investigating the integration of the cytoskeleton, intracellular transport, and organismal physiology to regulate the developmental processes of tube extension, lumen formation, and lumen diameter regulation in a narrow single cell. Multiple studies have provided new understanding of actin and intermediate filament cytoskeletal elements, vesicle transport, and the role of vacuolar ATPase in determining tube size. Most of the genes discovered have clear homologues in humans, with implications for understanding these processes in mammalian tissues such as Schwann cells, renal tubules, and brain vasculature. The results of several new genetic screens are described that provide a host of new targets for future studies in this informative structure.

Unique homeobox codes delineate all the neuron classes of C. elegans

It is presently not known whether neuronal cell type diversity, defined by cell type-specific anatomical, biophysical, functional and molecular signatures, can be reduced to relatively simple molecular descriptors of neuronal identity ¹. Examination of the expression of all conserved homeodomain proteins encoded by the Caenorhabditis elegans genome ² reveals that the complete set of 118 C. elegans neuron classes can be described individually by unique combinations of homeodomain protein expression, thereby providing the simplest currently known descriptor of neuronal diversity. Computational as well as genetic loss of function analyses corroborate that homeodomain proteins not only provide unique descriptors of neuron type, but also play a critical role specifying neuronal identity. We speculate that the pervasive employment of homeobox genes in defining unique neuronal identities reflects the evolutionary history of neuronal cell-type specification.

LSM2-8 and XRN-2 contribute to the silencing of H3K27me3-marked genes through targeted RNA decay

In fission yeast and plants, RNA-processing and degradation contribute to heterochromatin silencing, alongside conserved pathways of transcriptional repression. It was unknown if similar pathways exist in metazoans. Here we describe a pathway of silencing in C. elegans somatic cells, in which the highly conserved RNA binding complex LSM2-8 contributes selectively to the repression of heterochromatic reporters and endogenous genes bearing the Polycomb mark, histone H3K27me3. It acts by degrading selected transcripts through the XRN-2 exoribonuclease. Disruption of the LSM2-8 pathway leads to mRNA stabilization. Unlike previously described pathways of heterochromatic RNA degradation, LSM2-8-mediated RNA degradation does not require nor deposit H3K9 methylation. Rather, loss of this pathway coincides with a localized reduction in H3K27me3 at lsm-8-sensitive loci. Thus, we have uncovered a mechanism of RNA degradation that selectively contributes to the silencing of a subset of H3K27me3-marked genes, revealing a previously unrecognized layer of post-transcriptional control in metazoan heterochromatin.

How Weird is The Worm? Evolution of the Developmental Gene Toolkit in Caenorhabditis elegans

Comparative developmental biology and comparative genomics are the cornerstones of evolutionary developmental biology. Decades of fruitful research using nematodes have produced detailed accounts of the developmental and genomic variation in the nematode phylum. Evolutionary developmental biologists are now utilising these data as a tool with which to interrogate the evolutionary basis for the similarities and differences observed in Nematoda. Nematodes have often seemed atypical compared to the rest of the animal kingdom-from their totally lineage-dependent mode of embryogenesis to their abandonment of key toolkit genes usually deployed by bilaterians for proper development-worms are notorious rule breakers of the bilaterian handbook. However, exploring the nature of these deviations is providing answers to some of the biggest questions about the evolution of animal development. For example, why is the evolvability of each embryonic stage not the same? Why can evolution sometimes tolerate the loss of genes involved in key developmental events? Lastly, why does natural selection act to radically diverge toolkit genes in number and sequence in certain taxa? In answering these questions, insight is not only being provided about the evolution of nematodes, but of all metazoans.

Novel exc Genes Involved in Formation of the Tubular Excretory Canals of Caenorhabditis elegans

Regulation of luminal diameter is critical to the function of small single-celled tubes, of which the seamless tubular excretory canals of Caenorhabditis elegans provide a tractable genetic model. Mutations in several sets of genes exhibit the Exc phenotype, in which canal luminal growth is visibly altered. Here, a focused reverse genomic screen of genes highly expressed in the canals found 18 genes that significantly affect luminal outgrowth or diameter. These genes encode novel proteins as well as highly conserved proteins involved in processes including gene expression, cytoskeletal regulation, and vesicular and transmembrane transport. In addition, two genes act as suppressors on a pathway of conserved genes whose products mediate vesicle movement from early to recycling endosomes. The results provide new tools for understanding the integration of cytoplasmic structure and physiology in forming and maintaining the narrow diameter of single-cell tubules.

Phylogenetic and mutational analyses of human LEUTX, a homeobox gene implicated in embryogenesis

Recently, human PAIRED-LIKE homeobox transcription factor (TF) genes were discovered whose expression is limited to the period of embryo genome activation up to the 8-cell stage. One of these TFs is LEUTX, but its importance for human embryogenesis is still subject to debate. We confirmed that human LEUTX acts as a TAATCC-targeting transcriptional activator, like other K50-type PAIRED-LIKE TFs. Phylogenetic comparisons revealed that Leutx proteins are conserved across Placentalia and comprise two conserved domains, the homeodomain, and a Leutx-specific domain containing putative transcriptional activation motifs (9aaTAD). Examination of human genotype resources revealed 116 allelic variants in LEUTX. Twenty-four variants potentially affect function, but they occur only heterozygously at low frequency. One variant affects a DNA-specificity determining residue, mutationally reachable by a one-base transition. In vitro and in silico experiments showed that this LEUTX mutation (alanine to valine at position 54 in the homeodomain) results in a transactivational loss-of-function to a minimal TAATCC-containing promoter and a 36 bp motif enriched in genes involved in embryo genome activation. A compensatory change in residue 47 restores function. The results support the notion that human LEUTX functions as a transcriptional activator important for human embryogenesis.

Xenacoelomorpha Survey Reveals That All 11 Animal Homeobox Gene Classes Were Present in the First Bilaterians

Homeodomain transcription factors are involved in many developmental processes across animals and have been linked to body plan evolution. Detailed classifications of these proteins identified 11 distinct classes of homeodomain proteins in animal genomes, each harboring specific sequence composition and protein domains. Although humans contain the full set of classes, Drosophila melanogaster and Caenorhabditis elegans each lack one specific class. Furthermore, representative previous analyses in sponges, ctenophores, and cnidarians could not identify several classes in those nonbilaterian metazoan taxa. Consequently, it is currently unknown when certain homeodomain protein classes first evolved during animal evolution. Here, we investigate representatives of the sister group to all remaining bilaterians, the Xenacoelomorpha. We analyzed three acoel, one nemertodermatid, and one Xenoturbella transcriptomes and identified their expressed homeodomain protein content. We report the identification of representatives of all 11 classes of animal homeodomain transcription factors in Xenacoelomorpha and we describe and classify their homeobox genes relative to the established animal homeodomain protein families. Our findings suggest that the genome of the last common ancestor of bilateria contained the full set of these gene classes, supporting the subsequent diversification of bilaterians.

OFF-responses of interneurons optimize avoidance behaviors depending on stimulus strength via electrical synapses

Optimization of the types and timing of avoidance behaviors depending on the intensity of a noxious stimulus is essential for survival; however, processing in the central nervous system and its developmental basis are largely unknown. Here, we report that Caenorhabditis elegans preferentially selects one of three different types of avoidance behaviors depending on the strength of the noxious stimulus. We screened 210 neuronal transcription factors using a combination of optogenetics and RNA interference methods and identified 19 candidates required for avoidance behaviors. One candidate, gene lin-32 (abnormal cellLINeage32), which encodes an atonal homolog, is required for the neural fate determination of AIB interneurons and functions by regulating the expression of electrical and chemical synapse genes, namely, inx-1 (innexin 1) and AMPA-type ionotropic glutamate receptor glr-1. When examined by Ca imaging, AIB showed an OFF calcium increase to the noxious stimulus. The OFF calcium increase was provoked only by strong stimulation, suggesting a role for optimization of the avoidance behavior. However, lin-32 mutants showed an impaired AIB OFF calcium increase, concomitant with a reduced occurrence of the dynamic avoidance behavior called the "omega turn". The AIB neural responses may be transferred to downstream inter/motor neurons projecting to the neck muscles via electrical synapses comprising inx-1. Finally, we found a correlation between powerful contractions of the neck muscles and omega turns. Thus, the central regulation of the magnitude and timing of activation of the AIB interneurons optimizes the probability of omega turn depending on the stimulus context.

For survival, animals exhibit appropriate behaviors depending on the type and strength of a stimulus; however, processing in the central nervous system and its developmental basis are largely unknown. Here, we performed RNA interference screening of 210 transcription factors by observing the avoidance behaviors evoked by a harmful stimulus and identified 19 candidates with behavioral defects. One candidate gene, lin-32, which encodes a homologue of a well-conserved transcription factor, atonal, that plays a key role in cell type specification, is required for the optimization of avoidance behaviors. lin-32 regulates the expression of various synaptic genes in the key neurons required for the optimization of avoidance behaviors. The key neurons respond at the time of the removal of the strong stimulation, and an appropriate magnitude and timing of the calcium increase are required to determine the type and timing of dynamic avoidance behaviors. However, the key neuronal response does not occur in the lin-32 mutants. Activation of the key neurons may be transferred to the downstream neurons projecting to the neck muscles via electrical synapses. Finally, we show a correlation between powerful neck muscle contractions and dynamic avoidance behaviors.

Expression and clinical significance of HOXA5, E-cadherin, and β-catenin in cervical squamous cell carcinoma

OBJECTIVES

HOXA5 has been identified as a biomarker in pathogenesis of several cancers, such as non-small cell lung cancer (NSCLC) and breast cancer cells. The role has not been explored in cervical squamous cell carcinoma (CSCC).

METHODS

Tissues of 120 cases with CSCC and 30 controls with chronic cervicitis were constructed from our archived surgical pathology files and staining with HOXA5. Additional antibodies to E-cadherin and β-catenin were stained for comparison. For each marker, low expression was defined as staining score 0 to 3 points, whereas high expression referred to 4 points and above. Fifty-four patients in this research with cervical cancer were followed up for prognostic assessment.

RESULT

HOXA5 had high expression in chronic cervicitis and low in CSCC (P=0.004). The positivity rates of HOXA5 in patients without muscular layer invasion (MLI) and lymphatic invasion (LI) was higher than that in metastasis (113 vs. 17; 117 vs. 3). Consistently, low expression of HOXA5 was more common in poorly differentiated carcinoma, CSCC subjects without MLI and LI. Expression of E-cadherin and β-catenin was parallel with the expression of HOXA5. Additionally, patients with higher expression of HOXA5 had much more favorable prognosis than those with lower expression among follow up of the 54 patients.

CONCLUSION

In parallel with E-cadherin and β-catenin, low expression of HOXA5 was more common in CSCC patients with poor differentiation, without MLI and LI, among those which showed poor prognosis.

Homeodomain proteins: an update

Here, we provide an update of our review on homeobox genes that we wrote together with Walter Gehring in 1994. Since then, comprehensive surveys of homeobox genes have become possible due to genome sequencing projects. Using the 103 Drosophila homeobox genes as example, we present an updated classification. In animals, there are 16 major classes, ANTP, PRD, PRD-LIKE, POU, HNF, CUT (with four subclasses: ONECUT, CUX, SATB, and CMP), LIM, ZF, CERS, PROS, SIX/SO, plus the TALE superclass with the classes IRO, MKX, TGIF, PBC, and MEIS. In plants, there are 11 major classes, i.e., HD-ZIP (with four subclasses: I to IV), WOX, NDX, PHD, PLINC, LD, DDT, SAWADEE, PINTOX, and the two TALE classes KNOX and BEL. Most of these classes encode additional domains apart from the homeodomain. Numerous insights have been obtained in the last two decades into how homeodomain proteins bind to DNA and increase their specificity by interacting with other proteins to regulate cell- and tissue-specific gene expression. Not only protein-DNA base pair contacts are important for proper target selection; recent experiments also reveal that the shape of the DNA plays a role in specificity. Using selected examples, we highlight different mechanisms of homeodomain protein-DNA interaction. The PRD class of homeobox genes was of special interest to Walter Gehring in the last two decades. The PRD class comprises six families in Bilateria, and tinkers with four different motifs, i.e., the PAIRED domain, the Groucho-interacting motif EH1 (aka Octapeptide or TN), the homeodomain, and the OAR motif. Homologs of the co-repressor protein Groucho are also present in plants (TOPLESS), where they have been shown to interact with small amphipathic motives (EAR), and in yeast (TUP1), where we find an EH1-like motif in MATα2.