Methods for Estimating Personal Disease Risk and Phylogenetic Diversity of Hematopoietic Stem Cells.

An individual's chronological age does not always correspond to the health of different tissues in their body, especially in cases of disease. Therefore, estimating and contrasting the physiological age of tissues with an individual's chronological age may be a useful tool to diagnose disease and its progression. In this study, we present novel metrics to quantify the loss of phylogenetic diversity in hematopoietic stem cells (HSCs), which are precursors to most blood cell types and are associated with many blood-related diseases. These metrics showed an excellent correspondence with an age-related increase in blood cancer incidence, enabling a model to estimate the phylogeny-derived age (phyloAge) of HSCs present in an individual. The HSC phyloAge was generally older than the chronological age of patients suffering from myeloproliferative neoplasms (MPNs). We present a model that relates excess HSC aging with increased MPN risk. It predicted an over 200 times greater risk based on the HSC phylogenies of the youngest MPN patients analyzed. Our new metrics are designed to be robust to sampling biases and do not rely on prior knowledge of driver mutations or physiological assessments. Consequently, they complement conventional biomarker-based methods to estimate physiological age and disease risk.

TopHap: rapid inference of key phylogenetic structures from common haplotypes in large genome collections with limited diversity.

MOTIVATION: Building reliable phylogenies from very large collections of sequences with a limited number of phylogenetically informative sites is challenging because sequencing errors and recurrent/backward mutations interfere with the phylogenetic signal, confounding true evolutionary relationships. Massive global efforts of sequencing genomes and reconstructing the phylogeny of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strains exemplify these difficulties since there are only hundreds of phylogenetically informative sites but millions of genomes. For such datasets, we set out to develop a method for building the phylogenetic tree of genomic haplotypes consisting of positions harboring common variants to improve the signal-to-noise ratio for more accurate and fast phylogenetic inference of resolvable phylogenetic features. RESULTS: We present the TopHap approach that determines spatiotemporally common haplotypes of common variants and builds their phylogeny at a fraction of the computational time of traditional methods. We develop a bootstrap strategy that resamples genomes spatiotemporally to assess topological robustness. The application of TopHap to build a phylogeny of 68 057 SARS-CoV-2 genomes (68KG) from the first year of the pandemic produced an evolutionary tree of major SARS-CoV-2 haplotypes. This phylogeny is concordant with the mutation tree inferred using the co-occurrence pattern of mutations and recovers key phylogenetic relationships from more traditional analyses. We also evaluated alternative roots of the SARS-CoV-2 phylogeny and found that the earliest sampled genomes in 2019 likely evolved by four mutations of the most recent common ancestor of all SARS-CoV-2 genomes. An application of TopHap to more than 1 million SARS-CoV-2 genomes reconstructed the most comprehensive evolutionary relationships of major variants, which confirmed the 68KG phylogeny and provided evolutionary origins of major and recent variants of concern. AVAILABILITY AND IMPLEMENTATION: TopHap is available at https://github.com/SayakaMiura/TopHap. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

An Evolutionary Portrait of the Progenitor SARS-CoV-2 and Its Dominant Offshoots in COVID-19 Pandemic.

Global sequencing of genomes of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has continued to reveal new genetic variants that are the key to unraveling its early evolutionary history and tracking its global spread over time. Here we present the heretofore cryptic mutational history and spatiotemporal dynamics of SARS-CoV-2 from an analysis of thousands of high-quality genomes. We report the likely most recent common ancestor of SARS-CoV-2, reconstructed through a novel application and advancement of computational methods initially developed to infer the mutational history of tumor cells in a patient. This progenitor genome differs from genomes of the first coronaviruses sampled in China by three variants, implying that none of the earliest patients represent the index case or gave rise to all the human infections. However, multiple coronavirus infections in China and the United States harbored the progenitor genetic fingerprint in January 2020 and later, suggesting that the progenitor was spreading worldwide months before and after the first reported cases of COVID-19 in China. Mutations of the progenitor and its offshoots have produced many dominant coronavirus strains that have spread episodically over time. Fingerprinting based on common mutations reveals that the same coronavirus lineage has dominated North America for most of the pandemic in 2020. There have been multiple replacements of predominant coronavirus strains in Europe and Asia as well as continued presence of multiple high-frequency strains in Asia and North America. We have developed a continually updating dashboard of global evolution and spatiotemporal trends of SARS-CoV-2 spread (http://sars2evo.datamonkey.org/).

Molecular evidence for SARS-CoV-2 in samples collected from patients with morbilliform eruptions since late 2019 in Lombardy, northern Italy.

As a reference laboratory for measles and rubella surveillance in Lombardy, we evaluated the association between SARS-CoV-2 infection and measles-like syndromes, providing preliminary evidence for undetected early circulation of SARS-CoV-2. Overall, 435 samples from 156 cases were investigated. RNA from oropharyngeal swabs (N = 148) and urine (N = 141) was screened with four hemi-nested PCRs and molecular evidence for SARS-CoV-2 infection was found in 13 subjects. Two of the positive patients were from the pandemic period (2/12, 16.7%, March 2020-March 2021) and 11 were from the pre-pandemic period (11/44, 25%, August 2019-February 2020). Sera (N = 146) were tested for anti-SARS-CoV-2 IgG, IgM, and IgA antibodies. Five of the RNA-positive individuals also had detectable anti-SARS-CoV-2 antibodies. No strong evidence of infection was found in samples collected between August 2018 and July 2019 from 100 patients. The earliest sample with evidence of SARS-CoV-2 RNA was from September 12, 2019, and the positive patient was also positive for anti-SARS-CoV-2 antibodies (IgG and IgM). Mutations typical of B.1 strains previously reported to have emerged in January 2020 (C3037T, C14408T, and A23403G), were identified in samples collected as early as October 2019 in Lombardy. One of these mutations (C14408T) was also identified among sequences downloaded from public databases that were obtained by others from samples collected in Brazil in November 2019. We conclude that a SARS-CoV-2 progenitor capable of producing a measles-like syndrome may have emerged in late June-late July 2019 and that viruses with mutations characterizing B.1 strain may have been spreading globally before the first Wuhan outbreak. Our findings should be complemented by high-throughput sequencing to obtain additional sequence information. We highlight the importance of retrospective surveillance studies in understanding the early dynamics of COVID-19 spread and we encourage other groups to perform retrospective investigations to seek confirmatory proofs of early SARS-CoV-2 circulation.

PathFinder: Bayesian inference of clone migration histories in cancer.

SUMMARY: Metastases cause a vast majority of cancer morbidity and mortality. Metastatic clones are formed by dispersal of cancer cells to secondary tissues, and are not medically detected or visible until later stages of cancer development. Clone phylogenies within patients provide a means of tracing the otherwise inaccessible dynamic history of migrations of cancer cells. Here, we present a new Bayesian approach, PathFinder, for reconstructing the routes of cancer cell migrations. PathFinder uses the clone phylogeny, the number of mutational differences among clones, and the information on the presence and absence of observed clones in primary and metastatic tumors. By analyzing simulated datasets, we found that PathFinder performes well in reconstructing clone migrations from the primary tumor to new metastases as well as between metastases. It was more challenging to trace migrations from metastases back to primary tumors. We found that a vast majority of errors can be corrected by sampling more clones per tumor, and by increasing the number of genetic variants assayed per clone. We also identified situations in which phylogenetic approaches alone are not sufficient to reconstruct migration routes.In conclusion, we anticipate that the use of PathFinder will enable a more reliable inference of migration histories and their posterior probabilities, which is required to assess the relative preponderance of seeding of new metastasis by clones from primary tumors and/or existing metastases. AVAILABILITY AND IMPLEMENTATION: PathFinder is available on the web at https://github.com/SayakaMiura/PathFinder.

A new method for inferring timetrees from temporally sampled molecular sequences.

Pathogen timetrees are phylogenies scaled to time. They reveal the temporal history of a pathogen spread through the populations as captured in the evolutionary history of strains. These timetrees are inferred by using molecular sequences of pathogenic strains sampled at different times. That is, temporally sampled sequences enable the inference of sequence divergence times. Here, we present a new approach (RelTime with Dated Tips [RTDT]) to estimating pathogen timetrees based on a relative rate framework underlying the RelTime approach that is algebraic in nature and distinct from all other current methods. RTDT does not require many of the priors demanded by Bayesian approaches, and it has light computing requirements. In analyses of an extensive collection of computer-simulated datasets, we found the accuracy of RTDT time estimates and the coverage probabilities of their confidence intervals (CIs) to be excellent. In analyses of empirical datasets, RTDT produced dates that were similar to those reported in the literature. In comparative benchmarking with Bayesian and non-Bayesian methods (LSD, TreeTime, and treedater), we found that no method performed the best in every scenario. So, we provide a brief guideline for users to select the most appropriate method in empirical data analysis. RTDT is implemented for use via a graphical user interface and in high-throughput settings in the newest release of cross-platform MEGA X software, freely available from http://www.megasoftware.net.

Computational enhancement of single-cell sequences for inferring tumor evolution.

Motivation: Tumor sequencing has entered an exciting phase with the advent of single-cell techniques that are revolutionizing the assessment of single nucleotide variation (SNV) at the highest cellular resolution. However, state-of-the-art single-cell sequencing technologies produce data with many missing bases (MBs) and incorrect base designations that lead to false-positive (FP) and false-negative (FN) detection of somatic mutations. While computational methods are available to make biological inferences in the presence of these errors, the accuracy of the imputed MBs and corrected FPs and FNs remains unknown. Results: Using computer simulated datasets, we assessed the robustness performance of four existing methods (OncoNEM, SCG, SCITE and SiFit) and one new method (BEAM). BEAM is a Bayesian evolution-aware method that improves the quality of single-cell sequences by using the intrinsic evolutionary information in the single-cell data in a molecular phylogenetic framework. Overall, BEAM and SCITE performed the best. Most of the methods imputed MBs with high accuracy, but effective detection and correction of FPs and FNs is a challenge, especially for small datasets. Analysis of an empirical dataset shows that computational methods can improve both the quality of tumor single-cell sequences and their utility for biological inference. In conclusion, tumor cells descend from pre-existing cells, which creates evolutionary continuity in single-cell sequencing datasets. This information enables BEAM and other methods to correctly impute missing data and incorrect base assignments, but correction of FPs and FNs remains challenging when the number of SNVs sampled is small relative to the number of cells sequenced. Availability and implementation: BEAM is available on the web at https://github.com/SayakaMiura/BEAM.

Predicting clone genotypes from tumor bulk sequencing of multiple samples.

Motivation: Analyses of data generated from bulk sequencing of tumors have revealed extensive genomic heterogeneity within patients. Many computational methods have been developed to enable the inference of genotypes of tumor cell populations (clones) from bulk sequencing data. However, the relative and absolute accuracy of available computational methods in estimating clone counts and clone genotypes is not yet known. Results: We have assessed the performance of nine methods, including eight previously-published and one new method (CloneFinder), by analyzing computer simulated datasets. CloneFinder, LICHeE, CITUP and cloneHD inferred clone genotypes with low error (<5% per clone) for a majority of datasets in which the tumor samples contained evolutionarily-related clones. Computational methods did not perform well for datasets in which tumor samples contained mixtures of clones from different clonal lineages. Generally, the number of clones was underestimated by cloneHD and overestimated by PhyloWGS, and BayClone2, Canopy and Clomial required prior information regarding the number of clones. AncesTree and Canopy did not produce results for a large number of datasets. Overall, the deconvolution of clone genotypes from single nucleotide variant (SNV) frequency differences among tumor samples remains challenging, so there is a need to develop more accurate computational methods and robust software for clone genotype inference. Availability and implementation: CloneFinder is implemented in Python and is available from https://github.com/gstecher/CloneFinderAPI. Supplementary information: Supplementary data are available at Bioinformatics online.

Somatic evolutionary timings of driver mutations.

BACKGROUND: A unified analysis of DNA sequences from hundreds of tumors concluded that the driver mutations primarily occur in the earliest stages of cancer formation, with relatively few driver mutation events detected in the late-arising subclones. However, emerging evidence from the sequencing of multiple tumors and tumor regions per individual suggests that late-arising subclones with additional driver mutations are underestimated in single-sample analyses. METHODS: To test whether driver mutations generally map to early tumor development, we examined multi-regional tumor sequencing data from 101 individuals reported in 11 published studies. Following previous studies, we annotated mutations as early-arising when all tumors/regions had those mutations (ubiquitous). We then inferred the fraction of mutations occurring early and compared it with late-arising mutations that were found in only single tumors/regions. RESULTS: While a large fraction of driver mutations in tumors occurred relatively early in cancers, later driver mutations occurred at least as frequently as the early drivers in a substantial number of patients. This result was robust to many different approaches to annotate driver mutations. The relative frequency of early and late driver mutations varied among patients of the same cancer type and in different cancer types. We found that previous reports of the preponderance of early driver mutations were primarily informed by analysis of single tumor variant allele profiles, with which it is challenging to clearly distinguish between early and late drivers. CONCLUSIONS: The origin and preponderance of new driver mutations are not limited to early stages of tumor evolution, with different tumors and regions showing distinct driver mutations and, consequently, distinct characteristics. Therefore, tumors with extensive intratumor heterogeneity appear to have many newly acquired drivers.

Long-Term Outcomes of Suture vs. Ring Tricuspid Annuloplasty for Functional Tricuspid Regurgitation.

BACKGROUND: Whether there is a significant difference in the long-term surgical outcomes between suture annuloplasty and ring annuloplasty for the treatment of functional tricuspid regurgitation (TR) is still controversial. We compared the long-term outcomes of tricuspid annuloplasty (TAP) with and without an annuloplasty ring.Methods and Results:From January 1996 to December 2015, we consecutively enrolled 684 patients (mean age, 65.5 years; 60% women) undergoing TAP for functional TR: 312 underwent conventional suture annuloplasty (Group S) and 372 underwent ring annuloplasty (Group R). Baseline characteristics were comparable between the 2 groups, except for age and the prevalence of prior cardiac surgery. The mean follow-up period was 7.3 years (range, 0.3-20.3 years).There was no significant difference in overall survival, freedom from major adverse cardiac and cerebrovascular events, and freedom from recurrent moderate to severe TR between Groups S and R, although postoperative mean TR grade and sPAP were significantly lower in Group R. Multivariate analysis revealed that age, prior cardiac surgery, and preoperative severe TR were independent predictors of recurrent moderate or greater TR. CONCLUSIONS: There was no significant difference in the long-term surgical outcomes between ring and suture TAP, although postoperative mean TR grade was lower in the ring annuloplasty group.

Hypolipidemic effects of crude green tea polysaccharide on rats, and structural features of tea polysaccharides isolated from the crude polysaccharide.

Crude tea polysaccharide (crude TPS) was prepared from instant green tea by ethanol precipitation followed by ultrafiltration membrane treatment and its effects on blood lipid, liver lipid, and fecal lipid levels were examined with Sprague-Dawley rats fed a high-fat diet. Although crude TPS showed no effects on the serum lipid levels, it suppressed the liver lipid accumulation and increased the fecal excretion of dietary fat. Then, the structural features of crude TPS were investigated. After separation of crude TPS by DEAE-cellulose and gel-filtration column chromatography, two kinds of neutral tea polysaccharides (NTPS-LP and NTPS-HH) and an acidic polysaccharide (ATPS-MH) were obtained. According to monosaccharide composition, methylation, and NMR analyses, NTPS-LP, NPTS-HH, and ATPS-MH were presumed to be starch, arabinogalactan with ß-1,3-linked galactosyl backbone blanched at position 6 and with 1,5-linked arabinofuranosyl residues, and α-1,4-linked galacturonic acid backbone with arabinogalactan region, respectively.

Random fluctuation of selection coefficients and the extent of nucleotide variation in human populations.

It is well known that the selection coefficient of a mutant allele varies from generation to generation, and the effect of this factor on genetic variation has been studied by many theoreticians. However, no consensus has been reached. One group of investigators believes that fluctuating selection has an effect of enhancing genetic variation, whereas the other group contends that it has a diversity-reducing effect. In recent years, it has become possible to study this problem by using single nucleotide polymorphisms (SNPs) as well as exome sequence data. In this article we present the theoretical distributions of mutant nucleotide frequencies for the two models of fluctuating selection and then compare the distributions with the empirical distributions obtained from SNP and exome sequence data in human populations. Interestingly, both SNP and exome sequence data showed that the neutral mutation model fits the empirical distribution quite well. Furthermore, the mathematical models with diversity-enhancing and diversity-reducing effects also fit the empirical distribution reasonably well. This result implies that there is no need of distinguishing between the diversity-enhancing and diversity-reducing models of fluctuating selection and the nucleotide polymorphism in human populations can be explained largely by neutral mutations when long-term evolution is considered.

Fabrication Technology of Self-Dissolving Sodium Hyaluronate Gels Ultrafine Microneedles for Medical Applications with UV-Curing Gas-Permeable Mold.

Microneedles are of great interest in diverse fields, including cosmetics, drug delivery systems, chromatography, and biological sensing for disease diagnosis. Self-dissolving ultrafine microneedles of pure sodium hyaluronate hydrogels were fabricated using a UV-curing TiO2-SiO2 gas-permeable mold polymerized by sol-gel hydrolysis reactions in nanoimprint lithography processes under refrigeration at 5 °C, where thermal decomposition of microneedle components can be avoided. The moldability, strength, and dissolution behavior of sodium hyaluronate hydrogels with different molecular weights were compared to evaluate the suitability of ultrafine microneedles with a bottom diameter of 40 µm and a height of 80 µm. The appropriate molecular weight range and formulation of pure sodium hyaluronate hydrogels were found to control the dissolution behavior of self-dissolving ultrafine microneedles while maintaining the moldability and strength of the microneedles. This fabrication technology of ultrafine microneedles expands their possibilities as a next-generation technique for bioactive gels for controlling the blood levels of drugs and avoiding pain during administration.

Improving cellular phylogenies through the integrated use of mutation order and optimality principles.

The study of tumor evolution is being revolutionalized by single-cell sequencing technologies that survey the somatic variation of cancer cells. In these endeavors, reliable inference of the evolutionary relationship of single cells is a key step. However, single-cell sequences contain many errors and missing bases, which necessitate advancing standard molecular phylogenetics approaches for applications in analyzing these datasets. We have developed a computational approach that integratively applies standard phylogenetic optimality principles and patterns of co-occurrence of sequence variations to produce more expansive and accurate cellular phylogenies from single-cell sequence datasets. We found the new approach to also perform well for CRISPR/Cas9 genome editing datasets, suggesting that it can be useful for various applications. We apply the new approach to some empirical datasets to showcase its use for reconstructing recurrent mutations and mutational reversals as well as for phylodynamics analysis to infer metastatic cell migrations between tumors.

Bootstrap confidence for molecular evolutionary estimates from tumor bulk sequencing data.

Bulk sequencing is commonly used to characterize the genetic diversity of cancer cell populations in tumors and the evolutionary relationships of cancer clones. However, bulk sequencing produces aggregate information on nucleotide variants and their sample frequencies, necessitating computational methods to predict distinct clone sequences and their frequencies within a sample. Interestingly, no methods are available to measure the statistical confidence in the variants assigned to inferred clones. We introduce a bootstrap resampling approach that combines clone prediction and statistical confidence calculation for every variant assignment. Analysis of computer-simulated datasets showed the bootstrap approach to work well in assessing the reliability of predicted clones as well downstream inferences using the predicted clones (e.g., mapping metastatic migration paths). We found that only a fraction of inferences have good bootstrap support, which means that many inferences are tentative for real data. Using the bootstrap approach, we analyzed empirical datasets from metastatic cancers and placed bootstrap confidence on the estimated number of mutations involved in cell migration events. We found that the numbers of driver mutations involved in metastatic cell migration events sourced from primary tumors are similar to those where metastatic tumors are the source of new metastases. So, mutations with driver potential seem to keep arising during metastasis. The bootstrap approach developed in this study is implemented in software available at https://github.com/SayakaMiura/CloneFinderPlus.

Validation of the severe COVID-19 prognostic value of serum IL-6, IFN-λ3, CCL17, and calprotectin considering the timing of clinical need for prediction.

Although biomarkers to predict coronavirus disease 2019 (COVID-19) severity have been studied since the early pandemic, no clear guidelines on using them in clinical practice are available. Here, we examined the ability of four biomarkers to predict disease severity using conserved sera from COVID-19 patients who received inpatient care between January 1, 2020 and September 21, 2021 at the National Center for Global Health and Medicine, collected at the appropriate time for prediction. We predicted illness severity in two situations: 1) prediction of future oxygen administration for patients without oxygen support within 8 days of onset (Study 1) and 2) prediction of future mechanical ventilation support (excluding non-invasive positive pressure ventilation) or death of patients within 4 days of the start of oxygen administration (Study 2). Interleukin-6, IFN-λ3, thymus and activation-regulated chemokine, and calprotectin were measured retrospectively. Other laboratory and clinical information were collected from medical records. AUCs were calculated from ROC curves and compared for the predictive ability of the four biomarkers. Study 1 included 18 patients, five of whom had developed oxygen needs. Study 2 included 45 patients, 13 of whom required ventilator management or died. In Study 1, IFN-λ3 showed a good predictive ability with an AUC of 0.92 (95% CI 0.76-1.00). In Study 2, the AUC of each biomarker was 0.70-0.74. The number of biomarkers above the cutoff showed the possibility of good prediction with an AUC of 0.86 (95% CI 0.75-0.97). When two or more biomarkers were positive, sensitivity and specificity were 0.92 and 0.63, respectively. In terms of biomarker testing at times when prognostication may be clinically useful, IFN-λ3 was predictive of oxygenation demand and a combination of the four biomarkers was predictive of mechanical ventilator requirement.

MOCA for Integrated Analysis of Gene Expression and Genetic Variation in Single Cells.

In cancer, somatic mutations occur continuously, causing cell populations to evolve. These somatic mutations result in the evolution of cellular gene expression patterns that can also change due to epigenetic modifications and environmental changes. By exploring the concordance of gene expression changes with molecular evolutionary trajectories of cells, we can examine the role of somatic variation on the evolution of gene expression patterns. We present Multi-Omics Concordance Analysis (MOCA) software to jointly analyze gene expressions and genetic variations from single-cell RNA sequencing profiles. MOCA outputs cells and genes showing convergent and divergent gene expression patterns in functional genomics.

A phylogenetic approach to study the evolution of somatic mutational processes in cancer.

Cancer cell genomes change continuously due to mutations, and mutational processes change over time in patients, leaving dynamic signatures in the accumulated genomic variation in tumors. Many computational methods detect the relative activities of known mutation signatures. However, these methods may produce erroneous signatures when applied to individual branches in cancer cell phylogenies. Here, we show that the inference of branch-specific mutational signatures can be improved through a joint analysis of the collections of mutations mapped on proximal branches of the cancer cell phylogeny. This approach reduces the false-positive discovery rate of branch-specific signatures and can sometimes detect faint signatures. An analysis of empirical data from 61 lung cancer patients supports trends based on computer-simulated datasets for which the correct signatures are known. In lung cancer somatic variation, we detect a decreasing trend of smoking-related mutational processes over time and an increasing influence of APOBEC mutational processes as the tumor evolution progresses. These analyses also reveal patterns of conservation and divergence of mutational processes in cell lineages within patients.

Clone Phylogenetics Reveals Metastatic Tumor Migrations, Maps, and Models.

Dispersal routes of metastatic cells are not medically detected or even visible. A molecular evolutionary analysis of tumor variation provides a way to retrospectively infer metastatic migration histories and answer questions such as whether the majority of metastases are seeded from clones within primary tumors or seeded from clones within pre-existing metastases, as well as whether the evolution of metastases is generally consistent with any proposed models. We seek answers to these fundamental questions through a systematic patient-centric retrospective analysis that maps the dynamic evolutionary history of tumor cell migrations in many cancers. We analyzed tumor genetic heterogeneity in 51 cancer patients and found that most metastatic migration histories were best described by a hybrid of models of metastatic tumor evolution. Synthesizing across metastatic migration histories, we found new tumor seedings arising from clones of pre-existing metastases as often as they arose from clones from primary tumors. There were also many clone exchanges between the source and recipient tumors. Therefore, a molecular phylogenetic analysis of tumor variation provides a retrospective glimpse into general patterns of metastatic migration histories in cancer patients.

Fabrication of High-Resolution Fine Microneedles Derived from Hydrolyzed Hyaluronic Acid Gels in Vacuum Environment Imprinting Using Water Permeable Mold.

Hydrolyzed hyaluronic acid high-resolution fine microneedles of 13 µm in diameter and 24 µm in height were fabricated from hydrolyzed hyaluronic acid gels made in mixtures of water using vacuum environment imprint lithography processes with a water permeable mold. The gas traps of water and volatile solvents in the imprint materials cause transfer failure in the conventional water impermeable molds of quartz and metal. However, the water permeable mold allows the use of 67 wt% dilution water with high solubility to increase the fluidity of the hydrolyzed hyaluronic acid during the patterning of high-resolution fine microneedles for cosmetics and pharmaceuticals. This demonstration sets a new paradigm of functional pure gels for high-resolution nano-patterning applications with various cosmetic and pharmaceutical materials containing dilution water using a water permeable mold.