Histone H2B.8 compacts flowering plant sperm through chromatin phase separation.

Sperm chromatin is typically transformed by protamines into a compact and transcriptionally inactive state1,2. Sperm cells of flowering plants lack protamines, yet they have small, transcriptionally active nuclei with chromatin condensed through an unknown mechanism3,4. Here we show that a histone variant, H2B.8, mediates sperm chromatin and nuclear condensation in Arabidopsis thaliana. Loss of H2B.8 causes enlarged sperm nuclei with dispersed chromatin, whereas ectopic expression in somatic cells produces smaller nuclei with aggregated chromatin. This result demonstrates that H2B.8 is sufficient for chromatin condensation. H2B.8 aggregates transcriptionally inactive AT-rich chromatin into phase-separated condensates, which facilitates nuclear compaction without reducing transcription. Reciprocal crosses show that mutation of h2b.8 reduces male transmission, which suggests that H2B.8-mediated sperm compaction is important for fertility. Altogether, our results reveal a new mechanism of nuclear compaction through global aggregation of unexpressed chromatin. We propose that H2B.8 is an evolutionary innovation of flowering plants that achieves nuclear condensation compatible with active transcription.

Identification of cellular proteins associated with human cytomegalovirus (HCMV) DNA replication suggests novel cellular and viral interactions.

Upon entry of Human cytomegalovirus (HCMV) into the host cell, the viral genome is transported to the nucleus where it serves as a template for transcription and genome replication. Production of new viral genomes is a coordinated effort between viral and cellular proteins. While the core replication proteins are encoded by the virus, additional cellular proteins support the process of genome synthesis. We used accelerated native isolation of proteins on nascent DNA (aniPOND) to study protein dynamics on nascent viral DNA during HCMV infection. Using this method, we identified specific viral and cellular proteins that are associated with nascent viral DNA. These included transcription factors, transcriptional regulators, DNA damage and repair factors, and chromatin remodeling complexes. The association of these identified proteins with viral DNA was confirmed by immunofluorescent imaging, chromatin-immunoprecipitation analyses, and shRNA knockdown experiments. These data provide evidence for the requirement of cellular factors involved in HCMV replication.

Histone variant H2A.Z regulates zygotic genome activation.

During embryogenesis, the genome shifts from transcriptionally quiescent to extensively active in a process known as Zygotic Genome Activation (ZGA). In Drosophila, the pioneer factor Zelda is known to be essential for the progression of development; still, it regulates the activation of only a small subset of genes at ZGA. However, thousands of genes do not require Zelda, suggesting that other mechanisms exist. By conducting GRO-seq, HiC and ChIP-seq in Drosophila embryos, we demonstrate that up to 65% of zygotically activated genes are enriched for the histone variant H2A.Z. H2A.Z enrichment precedes ZGA and RNA Polymerase II loading onto chromatin. In vivo knockdown of maternally contributed Domino, a histone chaperone and ATPase, reduces H2A.Z deposition at transcription start sites, causes global downregulation of housekeeping genes at ZGA, and compromises the establishment of the 3D chromatin structure. We infer that H2A.Z is essential for the de novo establishment of transcriptional programs during ZGA via chromatin reorganization.

Systemic depletion of histone macroH2A1.1 boosts hippocampal synaptic plasticity and social behavior in mice.

Gene expression and epigenetic processes in several brain regions regulate physiological processes such as cognitive functions and social behavior. MacroH2A1.1 is a ubiquitous variant of histone H2A that regulates cell stemness and differentiation in various organs. Whether macroH2A1.1 has a modulatory role in emotional behavior is unknown. Here, we employed macroH2A1.1 knock-out (-/- ) mice to perform a comprehensive battery of behavioral tests, and an assessment of hippocampal synaptic plasticity (long-term potentiation) accompanied by whole hippocampus RNA sequencing. MacroH2A1.1-/- mice exhibit a stunningly enhancement both of sociability and of active stress-coping behavior, reflected by the increased social behavior in social activity tests and higher mobility time in the forced swim test, respectively. They also display an increased hippocampal synaptic plasticity, accompanied by significant neurotransmission transcriptional networks changes. These results suggest that systemic depletion of histone macroH2A1.1 supports an epigenetic control necessary for hippocampal function and social behavior.

A survey of human histone H1 subtypes interaction networks: Implications for histones H1 functioning.

To show a spectrum of histone H1 subtypes (H1.1-H1.5) activity realized through the protein-protein interactions, data selected from APID resources were processed with sequence-based bioinformatics approaches. Histone H1 subtypes participate in over half a thousand interactions with nuclear and cytosolic proteins (ComPPI database) engaged in the enzymatic activity and binding of nucleic acids and proteins (SIFTER tool). Small-scale networks of H1 subtypes (STRING network) have similar topological parameters (P > .05) which are, however, different for networks hubs between subtype H1.1 and H1.4 and subtype H1.3 and H1.5 (P < .05) (Cytoscape software). Based on enriched GO terms (g:Profiler toolset) of interacting proteins, molecular function and biological process of networks hubs is related to RNA binding and ribosome biogenesis (subtype H1.1 and H1.4), cell cycle and cell division (subtype H1.3 and H1.5) and protein ubiquitination and degradation (subtype H1.2). The residue propensity (BIPSPI predictor) and secondary structures of interacting surfaces (GOR algorithm) as well as a value of equilibrium dissociation constant (ISLAND predictor) indicate that a type of H1 subtypes interactions is transient in term of the stability and medium-strong in relation to the strength of binding. Histone H1 subtypes bind interacting partners in the intrinsic disorder-dependent mode (FoldIndex, PrDOS predictor), according to the coupled folding and binding and mutual synergistic folding mechanism. These results evidence that multifunctional H1 subtypes operate via protein interactions in the networks of crucial cellular processes and, therefore, confirm a new histone H1 paradigm relating to its functioning in the protein-protein interaction networks.

Primate-specific histone variants.

Canonical histones (H2A, H2B, H3, and H4) are present in all eukaryotes where they package genomic DNA and participate in numerous cellular processes, such as transcription regulation and DNA repair. In addition to the canonical histones, there are many histone variants, which have different amino acid sequences, possess tissue-specific expression profiles, and function distinctly from the canonical counterparts. A number of histone variants, including both core histones (H2A/H2B/H3/H4) and linker histones (H1/H5), have been identified to date. Htz1 (H2A.Z) and CENP-A (CenH3) are present from yeasts to mammals, and H3.3 is present from Tetrahymena to humans. In addition to the prevalent variants, others like H3.4 (H3t), H2A.Bbd, and TH2B, as well as several H1 variants, are found to be specific to mammals. Among them, H2BFWT, H3.5, H3.X, H3.Y, and H4G are unique to primates (or Hominidae). In this review, we focus on localization and function of primate- or hominidae-specific histone variants.

The evolution and functional divergence of the histone H2B family in plants.

Chromatin regulation of eukaryotic genomes depends on the formation of nucleosome complexes between histone proteins and DNA. Histone variants, which are diversified by sequence or expression pattern, can profoundly alter chromatin properties. While variants in histone H2A and H3 families are well characterized, the extent of diversification of histone H2B proteins is less understood. Here, we report a systematic analysis of the histone H2B family in plants, which have undergone substantial divergence during the evolution of each major group in the plant kingdom. By characterising Arabidopsis H2Bs, we substantiate this diversification and reveal potential functional specialization that parallels the phylogenetic structure of emergent clades in eudicots. In addition, we identify a new class of highly divergent H2B variants, H2B.S, that specifically accumulate during chromatin compaction of dry seed embryos in multiple species of flowering plants. Our findings thus identify unsuspected diverse properties among histone H2B proteins in plants that has manifested into potentially novel groups of histone variants.

The atypical histone variant H3.15 promotes callus formation in Arabidopsis thaliana.

Plants are capable of regenerating new organs after mechanical injury. The regeneration process involves genome-wide reprogramming of transcription, which usually requires dynamic changes in the chromatin landscape. We show that the histone 3 variant HISTONE THREE RELATED 15 (H3.15) plays an important role in cell fate reprogramming during plant regeneration in Arabidopsis H3.15 expression is rapidly induced upon wounding. Ectopic overexpression of H3.15 promotes cell proliferation to form a larger callus at the wound site, whereas htr15 mutation compromises callus formation. H3.15 is distinguished from other Arabidopsis histones by the absence of the lysine residue 27 that is trimethylated by the POLYCOMB REPRESSIVE COMPLEX 2 (PRC2) in constitutively expressed H3 variants. Overexpression of H3.15 promotes the removal of the transcriptional repressive mark H3K27me3 from chromatin, which results in transcriptional de-repression of downstream genes, such as WUSCHEL RELATED HOMEOBOX 11 (WOX11). Our results reveal a new mechanism for a release from PRC2-mediated gene repression through H3.15 deposition into chromatin, which is involved in reprogramming cell fate to produce pluripotent callus cells.

Molecular phylogeny of Caudofoveata (Mollusca) challenges traditional views.

The shell-less, worm-shaped Caudofoveata (=Chaetodermomorpha) is one of the least known groups of molluscs. The taxon consists of 141 recognized species found from intertidal environments to the deep-sea where they live burrowing in sediment. Evolutionary relationships of the group have been debated, but few studies based on morphological or molecular data have investigated the phylogeny of the group. Here we use molecular phylogenetics to resolve relationships among and within families of Caudofoveata. Phylogenetic analyses were performed using selected mitochondrial and nuclear genes from species from all recognized families of Caudofoveata. In resulting trees and contrary to traditional views, Prochaetodermatidae forms the sister clade to a clade containing the other two currently recognized families, Chaetodermatidae and Limifossoridae. The monophyly of Prochaetodermatidae is highly supported, but Limifossoridae and Chaetodermatidae are not recovered as monophyletic. Most of the caudofoveate genera are also not recovered as monophyletic in our analyses. Thus results from our molecular data suggest that the current classification of Caudofoveata is in need of revision, and indicate evolutionary scenarios that differ from previously proposed hypotheses based on morphology.

Immunocytochemical analysis of valproic acid induced histone H3 and H4 acetylation during differentiation of rat adipose derived stem cells into neuron-like cells.

Valproic acid (VPA) is an inhibitor of histone deacetylases (HDACs) that can regulate differentiation and proliferation of stem cells by epigenetic mechanisms. We investigated VPA induced histone H3 and H4 acetylation in adipose derived stem cells (ADSCs) transdifferentiated into neuron-like cells (NLCs). Rat ADSCs were transdifferentiated into neural stem cells (NSCs) that had been generated from neurospheres. The NSCs were differentiated into NLCs by induction with different concentrations of VPA at 24, 48 and 72 h. The NLCs were evaluated using anti-H3 and -H4 antibodies, and ADSCs, NSCs and NLCs were evaluated using immunofluorescence. The ADSCs were immunoreactive to CD90 and CD49d, but not to CD45 and CD31. Both the neurospheres and NSCs were immunostained with nestin and neurofilament 68. The neurospheres expressed Musashi1, Sox2 and Neu N genes as determined by RT-PCR. Our dose-response study indicated that the optimal concentration of VPA was 1 mM at 72 h. Histone acetylation levels of H3 and H4 immunostaining intensities in NLCs were significantly greater than for ADSCs and NSCs. VPA alters H4 and H3 acetylation immunoreactivities of ADSCs transdifferentiated into NLCs.

Increasing the Separation Capacity of Intact Histone Proteoforms Chromatography Coupling Online Weak Cation Exchange-HILIC to Reversed Phase LC UVPD-HRMS.

Top-down proteomics is an emerging analytical strategy to characterize combinatorial protein post-translational modifications (PTMs). However, sample complexity and small mass differences between chemically closely related proteoforms often limit the resolution attainable by separations employing a single liquid chromatographic (LC) principle. In particular, for ultramodified proteins like histones, extensive and time-consuming fractionation is needed to achieve deep proteoform coverage. Herein, we present the first online nanoflow comprehensive two-dimensional liquid chromatography (nLC×LC) platform top-down mass spectrometry analysis of histone proteoforms. The described two-dimensional LC system combines weak cation exchange chromatography under hydrophilic interaction LC conditions (i.e., charge- and hydrophilicity-based separation) with reversed phase liquid chromatography (i.e., hydrophobicity-based separation). The two independent chemical selectivities were run at nanoflows (300 nL/min) and coupled online with high-resolution mass spectrometry employing ultraviolet photodissociation (UVPD-HRMS). The nLC×LC workflow increased the number of intact protein masses observable relative to one-dimensional approaches and allowed characterization of hundreds of proteoforms starting from limited sample quantities (∼1.5 µg).

Exploring intrinsically disordered proteins in Chlamydomonas reinhardtii.

The content of intrinsically disordered protein (IDP) is related to organism complexity, evolution, and regulation. In the Plantae, despite their high complexity, experimental investigation of IDP content is lacking. We identified by mass spectrometry 682 heat-resistant proteins from the green alga, Chlamydomonas reinhardtii. Using a phosphoproteome database, we found that 331 of these proteins are targets of phosphorylation. We analyzed the flexibility propensity of the heat-resistant proteins and their specific features as well as those of predicted IDPs from the same organism. Their mean percentage of disorder was about 20%. Most of the IDPs (~70%) were addressed to other compartments than mitochondrion and chloroplast. Their amino acid composition was biased compared to other classic IDPs. Their molecular functions were diverse; the predominant ones were nucleic acid binding and unfolded protein binding and the less abundant one was catalytic activity. The most represented proteins were ribosomal proteins, proteins associated to flagella, chaperones and histones. We also found CP12, the only experimental IDP from C. reinhardtii that is referenced in disordered protein database. This is the first experimental investigation of IDPs in C. reinhardtii that also combines in silico analysis.

Modification of the histone tetramer at the H3-H3 interface impacts tetrasome conformations and dynamics.

Nucleosomes consisting of a short piece of deoxyribonucleic acid (DNA) wrapped around an octamer of histone proteins form the fundamental unit of chromatin in eukaryotes. Their role in DNA compaction comes with regulatory functions that impact essential genomic processes such as replication, transcription, and repair. The assembly of nucleosomes obeys a precise pathway in which tetramers of histones H3 and H4 bind to the DNA first to form tetrasomes, and two dimers of histones H2A and H2B are subsequently incorporated to complete the complex. As viable intermediates, we previously showed that tetrasomes can spontaneously flip between a left-handed and right-handed conformation of DNA-wrapping. To pinpoint the underlying mechanism, here we investigated the role of the H3-H3 interface for tetramer flexibility in the flipping process at the single-molecule level. Using freely orbiting magnetic tweezers, we studied the assembly and structural dynamics of individual tetrasomes modified at the cysteines close to this interaction interface by iodoacetamide (IA) in real time. While such modification did not affect the structural properties of the tetrasomes, it caused a 3-fold change in their flipping kinetics. The results indicate that the IA-modification enhances the conformational plasticity of tetrasomes. Our findings suggest that subnucleosomal dynamics may be employed by chromatin as an intrinsic and adjustable mechanism to regulate DNA supercoiling.

New World spittlebugs (Hemiptera: Cercopidae: Ischnorhininae): Dated molecular phylogeny, classification, and evolution of aposematic coloration.

The spittlebug family Cercopidae (Hemiptera: Auchenorrhyncha: Cicadomorpha: Cercopoidea) is distributed worldwide, with highest species diversity in the tropics. Several included species are economically important pests of major agricultural crops and cultivated pasture grasses. Taxonomically, Cercopidae is divided into two subfamilies: the paraphyletic Old World Cercopinae and the monophyletic New World Ischnorhininae. Results are here presented from an investigation of phylogenetic relationships within Ischnorhininae based on DNA sequences from seven loci (18S rDNA, 28S rDNA, Histone 2A, Histone 3, Wingless, Cytochrome Oxidase I, and Cytochrome Oxidase II) generated from exemplars of 119 spittlebug species. The resulting topology is used to test alternative higher-level classification hypotheses of Ischnorhininae and, with fossil-calibration, dates were estimated for major events in the evolutionary history of Cercopidae, including a much earlier divergence date (around 68-50 Mya) than previously reported in the literature. In addition, for the first time in Cercopidae, ancestral states of some predation avoidances strategies were reconstructed, with results suggesting an origin of aposematic coloration in the Cercopidae ancestor, with subsequent independent losses of aposematic coloration in multiple lineages.

A preliminary molecular phylogeny of shield-bearer moths (Lepidoptera: Adeloidea: Heliozelidae) highlights rich undescribed diversity.

Heliozelidae are a widespread, evolutionarily early diverging family of small, day-flying monotrysian moths, for which a comprehensive phylogeny is lacking. We generated the first molecular phylogeny of the family using DNA sequences of two mitochondrial genes (COI and COII) and two nuclear genes (H3 and 28S) from 130 Heliozelidae specimens, including eight of the twelve known genera: Antispila, Antispilina, Coptodisca, Heliozela, Holocacista, Hoplophanes, Pseliastis, and Tyriozela. Our results provide strong support for five major Heliozelidae clades: (i) a large widespread clade containing the leaf-mining genera Antispilina, Coptodisca and Holocacista and some species of Antispila, (ii) a clade containing most of the described Antispila, (iii) a clade containing the leaf-mining genus Heliozela and the monotypic genus Tyriozela, (iv) an Australian clade containing Pseliastis and (v) an Australian clade containing Hoplophanes. Each clade includes several new species and potentially new genera. Collectively, our data uncover a rich and undescribed diversity that appears to be especially prevalent in Australia. Our work highlights the need for a major taxonomic revision of the family and for generating a robust molecular phylogeny using multi-gene approaches in order to resolve the relationships among clades.

Molecular phylogeny of Pasiphaeidae (Crustacea, Decapoda, Caridea) reveals systematic incongruence of the current classification.

Caridean shrimps constitute one of the most diverse groups of decapod crustaceans, notwithstanding their poorly resolved infraordinal relationships. One of the systematically controversial families in Caridea is the predominantly pelagic Pasiphaeidae, comprises 101 species in seven genera. Pasiphaeidae species exhibit high morphological disparity, as well as ecological niche width, inhabiting shallow to very deep waters (>4000m). The present work presents the first molecular phylogeny of the family, based on a combined dataset of six mitochondrial and nuclear gene markers (12S rDNA, 16S rDNA, histone 3, sodium-potassium ATPase α-subunit, enolase and ATP synthase ß-subunit) from 33 species belonged to six genera of Pasiphaeidae with 19 species from 12 other caridean families as outgroup taxa. Maximum likelihood and Bayesian inference analyses conducted on the concatenated dataset of 2265bp suggest the family Pasiphaeidae is not monophyletic, with Psathyrocaris more closely related to other carideans than to the other five pasiphaeid genera included in this analysis. Leptochela occupies a sister position to the remaining genera and is genetically quite distant from them. At the generic level, the analysis supports the monophyly of Pasiphaea, Leptochela and Psathyrocaris, while Eupasiphae is shown to be paraphyletic, closely related to Parapasiphae and Glyphus. The present molecular result strongly implies that certain morphological characters used in the present systematic delineation within Pasiphaeidae may not be synapomorphies and the classification within the family needs to be urgently revised.

Muscle-relevant genes marked by stable H3K4me2/3 profiles and enriched MyoD binding during myogenic differentiation.

Post-translational modifications of histones play a key role in the regulation of gene expression during development and differentiation. Numerous studies have shown the dynamics of combinatorial regulation by transcription factors and histone modifications, in the sense that different combinations lead to distinct expression outcomes. Here, we investigated gene regulation by stable enrichment patterns of histone marks H3K4me2 and H3K4me3 in combination with the chromatin binding of the muscle tissue-specific transcription factor MyoD during myogenic differentiation of C2C12 cells. Using k-means clustering, we found that specific combinations of H3K4me2/3 profiles over and towards the gene body impact on gene expression and marks a subset of genes important for muscle development and differentiation. By further analysis, we found that the muscle key regulator MyoD was significantly enriched on this subset of genes and played a repressive role during myogenic differentiation. Among these genes, we identified the pluripotency gene Patz1, which is repressed during myogenic differentiation through direct binding of MyoD to promoter elements. These results point to the importance of integrating histone modifications and MyoD chromatin binding for coordinated gene activation and repression during myogenic differentiation.

Effects of tectonics and large scale climatic changes on the evolutionary history of Hyalomma ticks.

Hyalomma Koch, 1844 are ixodid ticks that infest mammals, birds and reptiles, to which 27 recognized species occur across the Afrotropical, Palearctic and Oriental regions. Despite their medical and veterinary importance, the evolutionary history of the group is enigmatic. To investigate various taxonomic hypotheses based on morphology, and also some of the mechanisms involved in the diversification of the genus, we sequenced and analysed data derived from two mtDNA fragments, three nuclear DNA genes and 47 morphological characters. Bayesian and Parsimony analyses based on the combined data (2242 characters for 84 taxa) provided maximum resolution and strongly supported the monophyly of Hyalomma and the subgenus Euhyalomma Filippova, 1984 (including H. punt Hoogstraal, Kaiser and Pedersen, 1969). A predicted close evolutionary association was found between morphologically similar H. dromedarii Koch, 1844, H. somalicum Tonelli Rondelli, 1935, H. impeltatum Schulze and Schlottke, 1929 and H. punt, and together they form a sister lineage to H. asiaticum Schulze and Schlottke, 1929, H. schulzei Olenev, 1931 and H. scupense Schulze, 1919. Congruent with morphological suggestions, H. anatolicum Koch, 1844, H. excavatum Koch, 1844 and H. lusitanicum Koch, 1844 form a clade and so also H. glabrum Delpy, 1949, H. marginatum Koch, 1844, H. turanicum Pomerantzev, 1946 and H. rufipes Koch, 1844. Wide scale continental sampling revealed cryptic divergences within African H. truncatum Koch, 1844 and H. rufipes and suggested that the taxonomy of these lineages is in need of a revision. The most basal lineages in Hyalomma represent taxa currently confined to Eurasia and molecular clock estimates suggest that members of the genus started to diverge approximately 36.25 million years ago (Mya). The early diversification event coincides well with the collision of the Indian and Eurasian Plates, an event that was also characterized by large scale faunal turnover in the region. Using S-Diva, we also propose that the closure of the Tethyan seaway allowed for the genus to first enter Africa approximately 17.73Mya. In concert, our data supports the notion that tectonic events and large scale global changes in the environment contributed significantly to produce the rich species diversity currently found in the genus Hyalomma.

Pleistocene range shifts, refugia and the origin of widespread species in western Palaearctic water beetles.

Quaternary glacial cycles drove major shifts in both the extent and location of the geographical ranges of many organisms. During glacial maxima, large areas of central and northern Europe were inhospitable to temperate species, and these areas are generally assumed to have been recolonized during interglacials by range expansions from Mediterranean refugia. An alternative is that this recolonization was from non-Mediterranean refugia, in central Europe or western Asia, but data on the origin of widespread central and north European species remain fragmentary, especially for insects. We studied three widely distributed lineages of freshwater beetles (the Platambus maculatus complex, the Hydraena gracilis complex, and the genus Oreodytes), all restricted to running waters and including both narrowly distributed southern endemics and widespread European species, some with distributions spanning the Palearctic. Our main goal was to determine the role of the Pleistocene glaciations in shaping the diversification and current distribution of these lineages. We sequenced four mitochondrial and two nuclear genes in populations drawn from across the ranges of these taxa, and used Bayesian probabilities and Maximum Likelihood to reconstruct their phylogenetic relationships, age and geographical origin. Our results suggest that all extant species in these groups are of Pleistocene origin. In the H. gracilis complex, the widespread European H. gracilis has experienced a rapid, recent range expansion from northern Anatolia, to occupy almost the whole of Europe. However, in the other two groups widespread central and northern European taxa appear to originate from central Asia, rather than the Mediterranean. These widespread species of eastern origin typically have peripherally isolated forms in the southern Mediterranean peninsulas, which may be remnants of earlier expansion-diversification cycles or result from incipient isolation of populations during the most recent Holocene expansion. The accumulation of narrow endemics of such lineages in the Mediterranean may result from successive cycles of range expansion, with subsequent speciation (and local extinction in glaciated areas) through multiple Pleistocene climatic cycles.

Molecular phylogeny of the Achatinoidea (Mollusca: Gastropoda).

This study presents a multi-gene phylogenetic analysis of the Achatinoidea and provides an initial basis for a taxonomic re-evaluation of family level groups within the superfamily. A total of 5028 nucleotides from the nuclear rRNA, actin and histone 3 genes and the 1st and 2nd codon positions of the mitochondrial cytochrome c oxidase subunit I gene were sequenced from 24 species, representing six currently recognised families. Results from maximum likelihood, neighbour joining, maximum parsimony and Bayesian inference trees revealed that, of currently recognised families, only the Achatinidae are monophyletic. For the Ferussaciidae, Ferussacia folliculus fell separately to Cecilioides gokweanus and formed a sister taxon to the rest of the Achatinoidea. For the Coeliaxidae, Coeliaxis blandii and Pyrgina umbilicata did not group together. The Subulinidae was not resolved, with some subulinids clustering with the Coeliaxidae and Thyrophorellidae. Three subfamilies currently included within the Subulinidae based on current taxonomy likewise did not form monophyletic groups.