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.

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.

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.

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.

Chromatin remodelling beyond transcription: the INO80 and SWR1 complexes.

Chromatin-modifying factors have essential roles in DNA processing pathways that dictate cellular functions. The ability of chromatin modifiers, including the INO80 and SWR1 chromatin-remodelling complexes, to regulate transcriptional processes is well established. However, recent studies reveal that the INO80 and SWR1 complexes have crucial functions in many other essential processes, including DNA repair, checkpoint regulation, DNA replication, telomere maintenance and chromosome segregation. During these diverse nuclear processes, the INO80 and SWR1 complexes function cooperatively with their histone substrates, gamma-H2AX and H2AZ. This research reveals that INO80 and SWR1 ATP-dependent chromatin remodelling is an integral component of pathways that maintain genomic integrity.

The nonhistone, N-terminal tail of an essential, chimeric H2A variant regulates mitotic H3-S10 dephosphorylation.

H2A.Y is an essential, divergent Tetrahymena thermophila histone variant. It has a long nonhistone N terminus that contains leucine-rich repeats (LRR) and an LRR cap domain with similarity to Sds22p, a regulator of yeast protein phosphatase 1 (PP1) activity in the nucleus. In growing cells, H2A.Y is incorporated into micronuclei only during S phase, which occurs immediately after micronuclear mitosis. Depletion of H2A.Y causes prolonged retention of mitosis-associated histone H3-S10 phosphorylation and mitotic abnormalities that mimic S10E mutation. In cells where H2A.Y is depleted, an inducible chimeric gene, in which the H2A.Y N terminus is attached to H2A.X, is shown to regulate micronuclear H3-S10 phosphorylation. H2A.Y can also be specifically coimmunoprecipitated with a Tetrahymena PP1 ortholog (Ppo1p). Taken together, these results argue that the N terminus of H2A.Y functions to regulate H3-S10 dephosphorylation. This striking in vivo case of "cross-talk" between a H2A variant and a specific post-translational modification of another histone demonstrates a novel function for a histone variant.

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.

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).

Phylogeny-Based Systematization of Arabidopsis Proteins with Histone H1 Globular Domain.

H1 (or linker) histones are basic nuclear proteins that possess an evolutionarily conserved nucleosome-binding globular domain, GH1. They perform critical functions in determining the accessibility of chromatin DNA to trans-acting factors. In most metazoan species studied so far, linker histones are highly heterogenous, with numerous nonallelic variants cooccurring in the same cells. The phylogenetic relationships among these variants as well as their structural and functional properties have been relatively well established. This contrasts markedly with the rather limited knowledge concerning the phylogeny and structural and functional roles of an unusually diverse group of GH1-containing proteins in plants. The dearth of information and the lack of a coherent phylogeny-based nomenclature of these proteins can lead to misunderstandings regarding their identity and possible relationships, thereby hampering plant chromatin research. Based on published data and our in silico and high-throughput analyses, we propose a systematization and coherent nomenclature of GH1-containing proteins of Arabidopsis (Arabidopsis thaliana [L.] Heynh) that will be useful for both the identification and structural and functional characterization of homologous proteins from other plant species.

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.

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.

Linker histones: History and current perspectives.

Although the overall structure of the fifth histone (linker histone, H1) is understood, its location on the nucleosome is only partially defined. Whilst it is clear that H1 helps condense the chromatin fibre, precisely how this is achieved remains to be determined. H1 is not a general gene repressor in that although it must be displaced from transcription start sites for activity to occur, there is only partial loss along the body of genes. How the deposition and removal of H1 occurs in particular need of further study. Linker histones are highly abundant nuclear proteins about which we know too little.

Specificities and genomic distribution of somatic mammalian histone H1 subtypes.

Histone H1 is a structural component of chromatin that may have a role in the regulation of chromatin dynamics. Unlike core histones, the linker histone H1 family is evolutionarily diverse and many organisms have multiple H1 variants or subtypes, distinguishable between germ-line and somatic cells. In mammals, the H1 family includes seven somatic H1 variants with a prevalence that varies between cell types and over the course of differentiation, H1.1 to H1.5 being expressed in a replication-dependent manner, whereas H1.0 and H1X are replication-independent. Until recently, it has not been known whether the different variants had specific roles in the regulation of nuclear processes or were differentially distributed across the genome. To address this, an increasing effort has been made to investigate divergent features among H1 variants, regarding their structure, expression patterns, chromatin dynamics, post-translational modifications and genome-wide distribution. Although H1 subtypes seem to have redundant functions, several reports point to the idea that they are also differently involved in specific cellular processes. Initial studies investigating the genomic distribution of H1 variants have started to suggest that despite a wide overlap, different variants may be enriched or preferentially located at different chromatin types, but this may depend on the cell type, the relative abundance of the variants, the differentiation state of the cell, or whether cells are derived from a neoplastic process. Understanding the heterogeneity of the histone H1 family is crucial to elucidate their role in chromatin organization, gene expression regulation and other cellular processes.

Proteomic analysis reveals the differential histone programs between male germline cells and vegetative cells in Lilium davidii.

In flowering plants, male germline fate is determined after asymmetric division of the haploid microspore. Daughter cells have distinct fates: the generative cell (GC) undergoes further mitosis to generate sperm cells (SCs), and the vegetative cell (VC) terminally differentiates. However, our understanding of the mechanisms underlying germline development remains limited. Histone variants and modifications define chromatin states, and contribute to establishing and maintaining cell identities by affecting gene expression. Here, we constructed a lily protein database, then extracted and detailed histone entries into a comprehensive lily histone database. We isolated large amounts of nuclei from VCs, GCs and SCs from lily, and profiled histone variants of all five histone families in all three cell types using proteomics approaches. We revealed 92 identities representing 32 histone variants: six for H1, 11 for H2A, eight for H2B, five for H3 and two for H4. Nine variants, including five H1, two H2B, one H3 and one H4 variant, specifically accumulated in GCs and SCs. We also detected H3 modification patterns in the three cell types. GCs and SCs had almost identical histone profiles and similar H3 modification patterns, which were significantly different from those of VCs. Our study also revealed the presence of multiple isoforms, and differential expression patterns between isoforms of a variant. The results suggest that differential histone programs between the germline and companion VCs may be established following the asymmetric division, and are important for identity establishment and differentiation of the male germline as well as the VC.

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.

Cryptic speciation and limited hybridization within Lumbricus earthworms (Clitellata: Lumbricidae).

Cryptic mitochondrial (mt) lineages are known to exist in the earthworm morphospecies Lumbricus rubellus and L. terrestris. The latter was recently split into two species, L. terrestris and L. herculeus, based on large genetic distances and a statistical difference in body size. There is support for the separation of some lineages in L. rubellus into species, whereas other lineages, separated by similar mt genetic distances, have been found to be part of the same species. However, no study has evaluated the status of the cryptic mt lineages in L. terrestris-L. herculeus and L. rubellus using nuclear genes. We use a combination of methods to reveal extensive cryptic speciation and limited hybridization in Lumbricus, based on one nuclear (H3) and one mitochondrial (COI) marker. Using a Bayesian multi-locus species delimitation method, as well as single gene haplotype networks and gene trees, we delimit seven well supported cryptic species within the morphospecies L. rubellus, and confirm the split within the species-pair L. terrestris-L. herculeus. Limited hybridization was found between the most common species of L. rubellus (A) in northern Europe and two other species (B and H) in this complex, as well as between L. terrestris and L. herculeus. Deep mt divergence was found within L. terrestris s.str. but no support for further splitting of this taxon was found. Both L. rubellus and L. terrestris are well studied model organisms, and considering that cryptic species and hybridization were found within them, it is important that they are properly identified in future studies.

Phylogenetic relationships of Mediterranean and North-East Atlantic Cantharidinae and notes on Stomatellinae (Vetigastropoda: Trochidae).

The subfamily Cantharidinae Gray, 1857 (Trochoidea: Trochidae) includes 23 recognized genera and over 200 known living species. These marine top shell snails are microphagous grazers that generally live in shallow rocky shores and in macroalgae and seagrass beds of sub-tropical and temperate waters from the Central and Western Indo-Pacific biogeographic regions to the Mediterranean Sea and the Eastern Atlantic Ocean. Recent molecular phylogenetic studies revising the family Trochidae supported the monophyly of the subfamily Cantharidinae and its sister group relationship to the subfamily Stomatellinae. These studies and others has thus far mostly focused on Indo-Pacific members of the subfamily Cantharidinae whereas here, we investigated phylogenetic relationships among their counterparts from the Mediterranean Sea and the North-eastern (NE) Atlantic Ocean including 33 species of genera Gibbula, Jujubinus, Phorcus, Clelandella, and Callumbonella. The Mediterranean and NE Atlantic taxa were supplemented with 30 Indo-Pacific Cantharidinae species plus 19 members of the sister group subfamily Stomatellinae. Phylogenetic trees were constructed using Bayesian inference and maximum likelihood with two datasets comprised of partial sequences of four or six mitochondrial (cox1, rrnL, rrnS, and cob) and nuclear (28S rRNA and histone H3) genes. A clade comprised of all Mediterranean and NE Atlantic taxa was recovered with high support, but its sister group among the Indo-Pacific lineages could not be determined with confidence (although the assignment of "Trochus" kotschyi to Priotrochus could be rejected). Within the Mediterranean and NE Atlantic clade, genera Phorcus and Jujubinus were recovered as reciprocally monophyletic, and the deep-sea genera Clelandella and Callumbonella were placed with high support as sister to Jujubinus. However, the genus Gibbula as currently defined was not monophyletic and constituent species were divided into three major clades and two independent lineages. Phylogenetic relationships among Phorcus, Jujubinus (plus Clelandella and Callumbonella), and the different clades of Gibbula were not fully resolved but received higher support in the phylogenetic analyses based on six genes. A first approach to resolve phylogenetic relationships within Stomatellinae was conducted showing that the diversity of the subfamily is highly underestimated at present, and that Calliotrochus is possibly a member of this subfamily. A chronogram was reconstructed using an uncorrelated relaxed lognormal molecular clock and the origin of the Mediterranean and NE Atlantic clade was dated right after the Azolla phase in the Middle Eocene about 48 million years ago whereas diversification of major clades (genera) followed the eastern closure of the Tethys Ocean in the Middle Miocene about 14 million years ago.