Magic angle spinning NMR of proteins: high-frequency dynamic nuclear polarization and (1)H detection.

Magic angle spinning (MAS) NMR studies of amyloid and membrane proteins and large macromolecular complexes are an important new approach to structural biology. However, the applicability of these experiments, which are based on (13)C- and (15)N-detected spectra, would be enhanced if the sensitivity were improved. Here we discuss two advances that address this problem: high-frequency dynamic nuclear polarization (DNP) and (1)H-detected MAS techniques. DNP is a sensitivity enhancement technique that transfers the high polarization of exogenous unpaired electrons to nuclear spins via microwave irradiation of electron-nuclear transitions. DNP boosts NMR signal intensities by factors of 10(2) to 10(3), thereby overcoming NMR's inherent low sensitivity. Alternatively, it permits structural investigations at the nanomolar scale. In addition, (1)H detection is feasible primarily because of the development of MAS rotors that spin at frequencies of 40 to 60 kHz or higher and the preparation of extensively (2)H-labeled proteins.

Ultrafast Magic Angle Spinning Solid-State NMR Spectroscopy: Advances in Methodology and Applications.

Solid-state NMR spectroscopy is one of the most commonly used techniques to study the atomic-resolution structure and dynamics of various chemical, biological, material, and pharmaceutical systems spanning multiple forms, including crystalline, liquid crystalline, fibrous, and amorphous states. Despite the unique advantages of solid-state NMR spectroscopy, its poor spectral resolution and sensitivity have severely limited the scope of this technique. Fortunately, the recent developments in probe technology that mechanically rotate the sample fast (100 kHz and above) to obtain "solution-like" NMR spectra of solids with higher resolution and sensitivity have opened numerous avenues for the development of novel NMR techniques and their applications to study a plethora of solids including globular and membrane-associated proteins, self-assembled protein aggregates such as amyloid fibers, RNA, viral assemblies, polymorphic pharmaceuticals, metal-organic framework, bone materials, and inorganic materials. While the ultrafast-MAS continues to be developed, the minute sample quantity and radio frequency requirements, shorter recycle delays enabling fast data acquisition, the feasibility of employing proton detection, enhancement in proton spectral resolution and polarization transfer efficiency, and high sensitivity per unit sample are some of the remarkable benefits of the ultrafast-MAS technology as demonstrated by the reported studies in the literature. Although the very low sample volume and very high RF power could be limitations for some of the systems, the advantages have spurred solid-state NMR investigation into increasingly complex biological and material systems. As ultrafast-MAS NMR techniques are increasingly used in multidisciplinary research areas, further development of instrumentation, probes, and advanced methods are pursued in parallel to overcome the limitations and challenges for widespread applications. This review article is focused on providing timely comprehensive coverage of the major developments on instrumentation, theory, techniques, applications, limitations, and future scope of ultrafast-MAS technology.

Molecular Mechanism of P53 Peptide Permeation through Lipid Membranes from Solid-State NMR Spectroscopy and Molecular Dynamics Simulations.

Macrocyclic peptides show promise in targeting high-value therapeutically relevant binding sites due to their high affinity and specificity. However, their clinical application is often hindered by low membrane permeability, which limits their effectiveness against intracellular targets. Previous studies focused on peptide conformations in various solvents, leaving a gap in understanding their interactions with and translocation through lipid bilayers. Addressing this, our study explores the membrane interactions of stapled peptides, a subclass of macrocyclic peptides, using solid-state nuclear magnetic resonance (ssNMR) spectroscopy and molecular dynamics (MD) simulations. We conducted ssNMR measurements on ATSP-7041M, a prototypical stapled peptide, to understand its interaction with lipid membranes, leading to an MD-informed model for peptide membrane permeation. Our findings reveal that ATSP-7041M adopts a stable α-helical structure upon membrane binding, facilitated by a cation-π interaction between its phenylalanine side chain and the lipid headgroup. This interaction makes the membrane-bound state energetically favorable, facilitating membrane affinity and insertion. The bound peptide displayed asymmetric insertion depths, with the C-terminus penetrating deeper (approximately 9 Å) than the N-terminus (approximately 4.3 Å) relative to the lipid headgroups. Contrary to expectations, peptide dynamics was not hindered by membrane binding and exhibited rapid motions similar to cell-penetrating peptides. These dynamic interactions and peptide-lipid affinity appear to be crucial for membrane permeation. MD simulations indicated a thermodynamically stable transmembrane conformation of ATSP-7041M, reducing the energy barrier for translocation. Our study offers an in silico view of ATSP-7041M's translocation from the extracellular to the intracellular region, highlighting the significance of peptide-lipid interactions and dynamics in enabling peptide transit through membranes.

Probing Molecular Packing of Lipid Nanoparticles from 31P Solution and Solid-State NMR.

Lipid nanoparticles (LNPs) are intricate multicomponent systems widely recognized for their efficient delivery of oligonucleotide cargo to host cells. Gaining insights into the molecular properties of LNPs is crucial for their effective design and characterization. However, analysis of their internal structure at the molecular level presents a significant challenge. This study introduces 31P nuclear magnetic resonance (NMR) methods to acquire structural and dynamic information about the phospholipid envelope of LNPs. Specifically, we demonstrate that the 31P chemical shift anisotropy (CSA) parameters serve as a sensitive indicator of the molecular assembly of distearoylphosphatidylcholine (DSPC) lipids within the particles. An analytical protocol for measuring 31P CSA is developed, which can be implemented using either solution NMR or solid-state NMR, offering wide accessibility and adaptability. The capability of this method is demonstrated using both model DSPC liposomes and real-world pharmaceutical LNP formulations. Furthermore, our method can be employed to investigate the impact of formulation processes and composition on the assembly of specifically LNP particles or, more generally, phospholipid-based delivery systems. This makes it an indispensable tool for evaluating critical pharmaceutical properties such as structural homogeneity, batch-to-batch reproducibility, and the stability of the particles.

Predicting the Stability of Formulations Containing Lyophilized Human Serum Albumin and Sucrose/Trehalose Using Solid-State NMR Spectroscopy.

Stabilization of proteins by disaccharides in lyophilized formulations depends on the interactions between the protein and the disaccharide (system homogeneity) and the sufficiently low mobility of the system. Human serum albumin (HSA) was lyophilized with disaccharides (sucrose and/or trehalose) in different relative concentrations. Solid-state nuclear magnetic resonance (ssNMR) spectroscopy 1H T1 and 1H T1ρ relaxation times were measured to determine the homogeneity of the lyophilized systems on 20-50 and 1-3 nm domains, respectively, with 1H T1 relaxation times also being used to determine the ß-relaxation rate. HSA/sucrose systems had longer 1H T1 relaxation times and were slightly more stable than HSA/trehalose systems in almost all cases shown. HSA/sucrose/trehalose systems have 1H T1 relaxation times between the HSA/sucrose and HSA/trehalose systems and did not result in a more stable system compared with binary systems. Inhomogeneity was evident in a sample containing relative concentrations of 10% HSA and 90% trehalose, suggesting trehalose crystallization during lyophilization. Under these stability conditions and with these ssNMR acquisition parameters, a 1H T1 relaxation time below 1.5 s correlated with an unstable sample, regardless of the disaccharide(s) used.

Insights into Factors Affecting Ethylene-Vinyl Acetate Copolymer Crystallinity in Islatravir Implant.

Islatravir, a highly potent nucleoside reverse transcriptase translocation inhibitor (NRTTI) for the treatment of HIV, has great potential to be formulated as ethylene-vinyl acetate (EVA) copolymer-based implants via hot melt extrusion. The crystallinity of EVA determines its physical and rheological properties and may impact the drug-eluting implant performance. Herein, we describe the systematic analysis of factors affecting the EVA crystallinity in islatravir implants. Differential scanning calorimetry (DSC) on EVA and solid-state NMR revealed drug loading promoted EVA crystallization, whereas BaSO4 loading had negligible impact on EVA crystallinity. The sterilization through γ-irradiation appeared to significantly impact the EVA crystallinity and surface characteristics of the implants. Furthermore, DSC analysis of thin implant slices prepared with an ultramicrotome indicated that the surface layer of the implant was more crystalline than the core. These findings provide critical insights into factors affecting the crystallinity, mechanical properties, and physicochemical properties of the EVA polymer matrix of extruded islatravir implants.

Small-Angle X-ray Scattering as a Powerful Tool for Phase and Crystallinity Assessment of Monoclonal Antibody Crystallites in Support of Batch Crystallization.

Crystalline suspensions of monoclonal antibodies (mAbs) have great potential to improve drug substance isolation and purification on a large scale and to be used for drug delivery via high-concentration formulations. Crystalline mAb suspensions are expected to have enhanced chemical and physical properties relative to mAb solutions delivered intravenously, making them attractive candidates for subcutaneous delivery. In contrast to small molecules, the development of protein crystalline suspensions is not a widely used approach in the pharmaceutical industry. This is mainly due to the challenges in finding crystalline hits and the suboptimal physical properties of the resulting crystallites when hits are found. Modern advances in instrumentation and increased knowledge of mAb crystallization have, however, resulted in higher probabilities of discovering crystal forms and improving their particle properties and characterization. In this regard, physical, analytical characterization plays a central role in the initial steps of understanding and later optimizing the crystallization of mAbs and requires careful selection of the appropriate tools. This contribution describes a novel crystal structure of the antibody pembrolizumab and demonstrates the usefulness of small-angle X-ray scattering (SAXS) for characterizing its crystalline suspensions. It illustrates the advantages of SAXS when used to (i) confirm crystallinity and crystal phase of crystallites produced in batch mode; (ii) confirm crystallinity under various conditions and detect variations in crystal phases, enabling fine-tuning of the crystallizations for phase control across multiple batches; (iii) monitor the physical response and stability of the crystallites in suspension with regard to filtration and washing; and (iv) monitor the physical stability of the crystallites upon drying. Overall, this work highlights how SAXS is an essential tool for mAb crystallization characterization.

Unraveling Pre-filled Syringe Needle Clogging: Exploring a Fresh Outlook Through Innovative Techniques.

OBJECTIVE: This study aimed to investigate the movement of liquid in the needle region of staked-in-needle pre-filled syringes using neutron imaging and synchrotron X-ray tomography. The objective was to gain insights into the dynamics of liquid presence and understand the factors contributing to needle clogging. METHODS: Staked-in-needle pre-filled syringes were examined using neutron radiography and synchrotron X-ray phase-contrast computed tomography. Neutron radiography provided a 2D visualization of liquid presence in the needle, while synchrotron X-ray tomography offered high-resolution 3D imaging to study detailed morphological features of the liquid. RESULTS: Neutron radiography revealed liquid presence in the needle region for as-received samples and after temperature and pressure cycling. Pressure cycling had a more pronounced effect on liquid formation. Synchrotron X-ray tomography confirmed the presence of liquid and revealed various morphologies, including droplets of different sizes, liquid segments blocking sections of the needle, and a thin layer covering the needle wall. Liquid presence was also observed between the steel needle and the glass barrel. CONCLUSIONS: The combination of neutron imaging and synchrotron X-ray tomography provided valuable insights into the dynamics of liquid movement in staked-in-needle pre-filled syringes. Temperature and pressure cycling were found to contribute to additional liquid formation, with pressure changes playing a significant role. The detailed morphological analysis enhanced the understanding of microstructural arrangements within the needle. This research contributes to addressing the issue of needle clogging and can guide the development of strategies to improve pre-filled syringe performance.

Deciphering molecular basis of pesticide-induced recurrent pregnancy loss: insights from transcriptomics analysis.

Recent studies have revealed a notable connection between pesticide exposure and Recurrent Pregnancy Loss (RPL), yet the precise molecular underpinning of this toxicity remains elusive. Through the alignment of Differentially Expressed Genes (DEGs) of healthy and RPL patients with the target genes of 9 pesticide components, we identified a set of 12 genes responsible for RPL etiology. Interestingly, biological process showed that besides RPL, those 12 genes also associated with preeclampsia and cardiovascular disease. Enrichment analysis showed the engagement of these genes associated with essential roles in the molecular transport of small molecules, as well as the aldosterone-regulated sodium reabsorption, endocrine and other factor-regulated calcium reabsorption, mineral absorption, ion homeostasis, and ion transport by P-type ATPases. Notably, the crosstalk targets between pesticide components played crucial roles in influencing RPL results, suggesting a role in attenuating pesticide agents that contribute to RPL. It is important to note that non-significant concentration of the pesticide components observed in both control and RPL samples should not prematurely undermine the potential for pesticides to induce RPL in humans. This study emphasizes the complexity of pesticide induced RPL and highlights avenues for further research and precautionary measures.

Solid-State NMR Spectroscopy to Probe State and Phase Transitions in Frozen Solutions.

Freezing is commonly encountered during the processing and storage of biomacromolecule products. Therefore, understanding the phase and state transitions in pharmaceutical frozen solutions is crucial for the rational development of biopharmaceuticals. Solid-state nuclear magnetic resonance spectroscopy (ssNMR) was used to analyze solutions containing sodium phosphate buffer, histidine, and trehalose. Upon freezing, crystallization of disodium phosphate hydrogen dodecahydrate (Na2HPO4·12H2O, DPDH) and histidine was identified using 31P and 13C ssNMR, respectively, and confirmed by synchrotron X-ray diffractometry (SXRD). Using histidine as a molecular probe and based on the chemical shifts of atoms of interest, the pH of the freeze concentrate was measured. The unfrozen water content in freeze concentrates was quantified by 1H single pulse experiments. 13C-insensitive nuclei enhancement by polarization transfer (INEPT) and cross-polarization (CP) experiments were used as orthogonal tools to characterize the solutes in a "mobile" and a more "solid-like" state in the freeze-concentrated solutions, respectively. The above analyses were applied to a commercial monoclonal antibody (mAb) formulation of dupilumab. This work further establishes ssNMR spectroscopy as a highly capable biophysical tool to investigate the attributes of biopharmaceuticals and thereby provide insights into process optimization and formulation development.

Characterizing Silicone Oil-Induced Protein Aggregation with Stimulated Raman Scattering Imaging.

Particles in biopharmaceutical products present high risks due to their detrimental impacts on product quality and safety. Identification and quantification of particles in drug products are important to understand particle formation mechanisms, which can help develop control strategies for particle formation during the formulation development and manufacturing process. However, existing analytical techniques such as microflow imaging and light obscuration measurement lack the sensitivity and resolution to detect particles with sizes smaller than 2 µm. More importantly, these techniques are not able to provide chemical information to determine particle composition. In this work, we overcome these challenges by applying the stimulated Raman scattering (SRS) microscopy technique to monitor the C-H Raman stretching modes of the proteinaceous particles and silicone oil droplets formed in the prefilled syringe barrel. By comparing the relative signal intensity and spectral features of each component, most particles can be classified as protein-silicone oil aggregates. We further show that morphological features are poor indicators of particle composition. Our method has the capability to quantify aggregation in protein therapeutics with chemical and spatial information in a label-free manner, potentially allowing high throughput screening or investigation of aggregation mechanisms.

Probing Molecular Packing of Amorphous Pharmaceutical Solids Using X-ray Atomic Pair Distribution Function and Solid-State NMR.

The structural investigation of amorphous pharmaceuticals is of paramount importance in comprehending their physicochemical stability. However, it has remained a relatively underexplored realm primarily due to the limited availability of high-resolution analytical tools. In this study, we utilized the combined power of X-ray pair distribution functions (PDFs) and solid-state nuclear magnetic resonance (ssNMR) techniques to probe the molecular packing of amorphous posaconazole and its amorphous solid dispersion at the molecular level. Leveraging synchrotron X-ray PDF data and employing the empirical potential structure refinement (EPSR) methodology, we unraveled the existence of a rigid conformation and discerned short-range intermolecular C-F contacts within amorphous posaconazole. Encouragingly, our ssNMR 19F-13C distance measurements offered corroborative evidence supporting these findings. Furthermore, employing principal component analysis on the X-ray PDF and ssNMR data sets enabled us to gain invaluable insights into the chemical nature of the intermolecular interactions governing the drug-polymer interplay. These outcomes not only furnish crucial structural insights facilitating the comprehension of the underlying mechanisms governing the physicochemical stability but also underscore the efficacy of synergistically harnessing X-ray PDF and ssNMR techniques, complemented by robust modeling strategies, to achieve a high-resolution exploration of amorphous structures.

Convolutional Neural Networks Enable Highly Accurate and Automated Subvisible Particulate Classification of Biopharmaceuticals.

Quantification of subvisible particles, which are generally defined as those ranging in size from 2 to 100 µm, is important as critical characteristics for biopharmaceutical formulation development. Micro Flow Imaging (MFI) provides quantifiable morphological parameters to study both the size and type of subvisible particles, including proteinaceous particles as well as non-proteinaceous features incl. silicone oil droplets, air bubble droplets, etc., thus enabling quantitative and categorical particle attribute reporting for quality control. However, limitations in routine MFI image analysis can hinder accurate subvisible particle classification. In this work, we custom-built a subvisible particle-aware Convolutional Neural Network, SVNet, which has a very small computational footprint, and achieves comparable performance to prior state-of-art image classification models. SVNet significantly improves upon current standard operating procedures for subvisible particulate assessments as confirmed by thorough real-world validation studies.

Effect of Buffer Salts on Physical Stability of Lyophilized and Spray-Dried Protein Formulations Containing Bovine Serum Albumin and Trehalose.

This study examined the effect of buffer salts on the physical stability of spray-dried and lyophilized formulations of a model protein, bovine serum albumin (BSA). BSA formulations with various buffers were dried by either lyophilization or spray drying. The protein powders were then characterized using solid-state Fourier transform infrared spectroscopy (ssFTIR), powder X-ray diffraction (PXRD), size exclusion chromatography (SEC), solid-state hydrogen/deuterium exchange with mass spectrometry (ssHDX-MS), and solid-state nuclear magnetic resonance spectroscopy (ssNMR). Particle characterizations such as Brunauer-Emmett-Teller (BET) surface area, particle size distribution, and particle morphology were also performed. Results from conventional techniques such as ssFTIR did not exhibit correlations with the physical stability of studied formulations. Deconvoluted peak areas of deuterated samples from the ssHDX-MS study showed a satisfactory correlation with the loss of the monomeric peak area measured by SEC (R2 of 0.8722 for spray-dried formulations and 0.8428 for lyophilized formulations) in the 90-day accelerated stability study conducted at 40°C. mDSC and PXRD was unable to measure phase separation in the samples right after drying. In contrast, ssNMR successfully detected the occurrence of phase separation between the succinic buffer component and protein in the lyophilized formulation, which results in a distribution of microenvironmental acidity and the subsequent loss of long-term stability. Moreover, our results suggested that buffer salts have less impact on physical stability for the spray-dried formulations than the lyophilized solids.

Assessing Antigen-Adjuvant Complex Stability Against Physical Stresses By wNMR.

This study applies an emerging analytical technology, wNMR (water proton nuclear magnetic resonance), to assess the stability of aluminum adjuvants and antigen-adjuvant complexes against physical stresses, including gravitation, flow and freeze/thaw. Results from wNMR are verified by conventional analytical technologies, including static light scattering and microfluidic imaging. The results show that wNMR can quickly and noninvasively determine whether an aluminum adjuvant or antigen-adjuvant complex sample has been altered by physical stresses.

Efficient analysis of pharmaceutical drug substances and products using a solid-state NMR CryoProbe.

Solid-state nuclear magnetic resonance (ssNMR) is a high-resolution and versatile spectroscopic tool for characterizing pharmaceutical solids. However, the inherent low sensitivity of NMR remains a significant challenge in the analysis of natural abundance drug substances and products. Here, we report, for the first time, the application of a CPMAS CryoProbe™ to improve the sensitivity of 13C and 15N detection by approximately 5 to 6 times for solid-state analysis of a commercial pharmaceutical drug posaconazole (POSA). The sensitivity enhancement enables two-dimensional (2D) 13C-13C and 1H-15N correlation experiments, which are otherwise time-prohibitive using regular MAS probes, for resonance assignment and structural elucidation. These polarization transfer and correlation experiments reveal drug-drug and drug-polymer interactions in amorphous POSA and its amorphous solid dispersion formulation. Our results demonstrated that the CPMAS CryoProbe™ can be widely applied for routine pharmaceutical analysis and advanced structural investigations with significantly enhanced efficiency and throughput.

Corrected solid-state 13 C nuclear magnetic resonance peak assignment and side-group quantification of hydroxypropyl methylcellulose acetyl succinate pharmaceutical excipients.

Hydroxypropyl methylcellulose acetyl succinate (HPMCAS) is widely used as a pharmaceutical excipient, making a detailed understanding of its tunable structure important for formulation design. Several recently reported peak assignments in the solid-state 13 C NMR spectrum of HPMCAS have been corrected here using peak integrals in quantitative spectra, spectral editing, empirical chemical-shift predictions based on solution NMR, and full spectrum simulation analogous to deconvolution. Unlike in cellulose, the strong peak at 84 ppm must be assigned to C2 and C3 methyl ethers, instead of regular C4 of cellulose. The proposed assignment of signals at <65 ppm to OCH sites, including C5 of cellulose, could not be confirmed. CH2 spectral editing showed two resolved OCH2 bands, a more intense one from O-CH2 ethers of C6 at >69 ppm and a smaller one from its esters and possibly residual CH2 -OH groups, near 63 ppm. The strong intensities of resolved signals of acetyl, succinoyl, and oxypropyl substituents indicated the substitution of >85% of the OH groups in HPMCAS. The side-group concentrations in three different grades of HPMCAS were quantified.

Quantifying Micromolar Crystallinity in Pharmaceutical Materials Utilizing 19F Solid-State NMR.

Quantifying low-level components in solid-state analysis presents a significant challenge for most thermal, diffractometric, vibrational, and spectroscopic techniques. In pharmaceutical analysis, identifying and quantifying the physical form of the drug substance in solid dosages is a critical task to ensure the quality of drug products. For example, recrystallization of active pharmaceutical ingredients in amorphous solid dispersions can compromise the stability and bioavailability of drug products. Herein, we have developed and demonstrated fluorine-19 solid-state nuclear magnetic resonance (19F ssNMR) methods and pushed the boundary to quantify minor crystalline contents in amorphous pharmaceuticals. Calibration curves suggest that 19F direct polarization and 1H-19F cross-polarization ssNMR can readily quantify 0.1% w/w crystalline compound I, a commercial fluorinated drug molecule developed by Merck & Co., Inc., Rahway, NJ, U.S.A., in its amorphous formulation. 1H-19F multiple cross-polarization (MultiCP) has been implemented, for the first time, and compared with conventional cross-polarization methods. Most importantly, a relaxation-filtered 19F ssNMR method was utilized to unambiguously identify and quantify as low as 0.04% w/w crystalline components, that is, 6 µmol in a 100 mg tablet at 25% drug loading, by suppressing the signal from the amorphous counterpart. Such a low level of detection offers high confidence and sensitivity to quantify trace amounts of phase change in pharmaceutical amorphous materials in the solid state, which can facilitate formulation development as well as quality control.

Gene Flow Increases Phylogenetic Structure and Inflates Cryptic Species Estimations: A Case Study on Widespread Philippine Puddle Frogs (Occidozyga laevis).

In cryptic amphibian complexes, there is a growing trend to equate high levels of genetic structure with hidden cryptic species diversity. Typically, phylogenetic structure and distance-based approaches are used to demonstrate the distinctness of clades and justify the recognition of new cryptic species. However, this approach does not account for gene flow, spatial, and environmental processes that can obfuscate phylogenetic inference and bias species delimitation. As a case study, we sequenced genome-wide exons and introns to evince the processes that underlie the diversification of Philippine Puddle Frogs-a group that is widespread, phenotypically conserved, and exhibits high levels of geographically based genetic structure. We showed that widely adopted tree- and distance-based approaches inferred up to 20 species, compared to genomic analyses that inferred an optimal number of five distinct genetic groups. Using a suite of clustering, admixture, and phylogenetic network analyses, we demonstrate extensive admixture among the five groups and elucidate two specific ways in which gene flow can cause overestimations of species diversity: 1) admixed populations can be inferred as distinct lineages characterized by long branches in phylograms; and 2) admixed lineages can appear to be genetically divergent, even from their parental populations when simple measures of genetic distance are used. We demonstrate that the relationship between mitochondrial and genome-wide nuclear $p$-distances is decoupled in admixed clades, leading to erroneous estimates of genetic distances and, consequently, species diversity. Additionally, genetic distance was also biased by spatial and environmental processes. Overall, we showed that high levels of genetic diversity in Philippine Puddle Frogs predominantly comprise metapopulation lineages that arose through complex patterns of admixture, isolation-by-distance, and isolation-by-environment as opposed to species divergence. Our findings suggest that speciation may not be the major process underlying the high levels of hidden diversity observed in many taxonomic groups and that widely adopted tree- and distance-based methods overestimate species diversity in the presence of gene flow. [Cryptic species; gene flow; introgression; isolation-by-distance; isolation-by-environment; phylogenetic network; species delimitation.].

Phylogenomic Analysis of Ultraconserved Elements Resolves the Evolutionary and Biogeographic History of Segmented Trapdoor Spiders.

The segmented trapdoor spiders (Liphistiidae) are the sole surviving family of the suborder Mesothelae, which forms the sister lineage to all other living spiders. Liphistiids have retained a number of plesiomorphic traits and their present-day distribution is limited to East and Southeast Asia. Studying this group has the potential to shed light on the deep evolutionary history of spiders, but the phylogeny and divergence times of the family have not been resolved with confidence. We performed phylogenomic and molecular dating analyses of 2765 ultraconserved element loci from 185 liphistiid taxa. Our analyses show that the crown group of Liphistiidae appeared in the mid-Cretaceous at 102 Ma (95% credibility interval 92-113 Ma), but it was not until the Neogene that much of the diversification within the family occurred in mainland Southeast and East Asia. This diversification was coincident with tectonic events such as the extension of the East Asian continental margin, as well as geological upheavals in Indochina induced by the collision between India and Asia. Our study highlights the important role of major tectonic events in shaping the evolutionary history, present-day diversity, and geographical distribution of mesothele and liphistiid spiders. [biogeography; concatenation; Liphistiidae; molecular dating; summary coalescent; UCEs.].