Genetic characterization of cucumber genetic resources in the NARO Genebank indicates their multiple dispersal trajectories to the East.

KEY MESSAGE: Genotyping-by-sequencing of 723 worldwide cucumber genetic resources revealed that cucumbers were dispersed eastward via at least three distinct routes, one to Southeast Asia and two from different directions to East Asia. The cucumber (Cucumis sativus) is an economically important vegetable crop cultivated and consumed worldwide. Despite its popularity, the manner in which cucumbers were dispersed from their origin in South Asia to the rest of the world, particularly to the east, remains a mystery due to the lack of written records. In this study, we performed genotyping-by-sequencing (GBS) on 723 worldwide cucumber accessions, mainly deposited in the Japanese National Agriculture and Food Research Organization (NARO) Genebank, to characterize their genetic diversity, relationships, and population structure. Analyses based on over 60,000 genome-wide single-nucleotide polymorphisms identified by GBS revealed clear genetic differentiation between Southeast and East Asian populations, suggesting that they reached their respective region independently, not progressively. A deeper investigation of the East Asian population identified two subpopulations with different fruit characteristics, supporting the traditional classification of East Asian cucumbers into two types thought to have been introduced by independent routes. Finally, we developed a core collection of 100 accessions representing at least 93.2% of the genetic diversity present in the entire collection. The genetic relationships and population structure, their associations with geographic distribution and phenotypic traits, and the core collection presented in this study are valuable resources for elucidating the dispersal history and promoting the efficient use and management of genetic resources for research and breeding in cucumber.

Chromatographic printed array strip (C-PAS) method for cultivar-specific identification of sweetpotato cultivars 'Beniharuka' and 'Fukumurasaki'.

Sweetpotato (Ipomoea batatas) cultivars grown in Japan are highly valued for their excellent sweetness, high quality, and good texture. The export volume of sweetpotato from Japan has been rising rapidly, with a 10-fold increase on a weight basis over the last 10 years. However, since sweetpotato is propagated vegetatively from storage roots, it is easy to cultivate and propagate this crop, prompting concerns that Japanese sweetpotato cultivars/lines are being exported overseas, cultivated without permission, or reimported. Therefore, a rapid and accurate cultivar identification methodology is needed. In this study, we comprehensively analyzed the insertion sites of Cl8 retrotransposon to develop a cultivar identification technique for the Japanese cultivars 'Beniharuka' and 'Fukumurasaki'. These two cultivars were successfully distinguished from other cultivars using a minimum of two marker sets. Using the chromatographic printed array strip (C-PAS) method for DNA signal detection, 'Beniharuka' and 'Fukumurasaki' can be precisely identified using a single strip of chromatographic paper based on multiplex DNA signals derived from the amplicons of the Cl8 insertion sites. Since this method can detect DNA signals in only ~15 minutes, we expect that our method will facilitate rapid, reliable, and convenient cultivar discrimination for on-site inspection of sweetpotato.

A target cultivar-specific identification system based on the chromatographic printed array strip method for eight prominent Japanese citrus cultivars.

Citrus is a major cultivated crop in Japan, and new cultivars are of great interest in the Japanese and global market. Recently, the infringement of breeders' rights to citrus cultivars bred in Japan has become a problem related to the agricultural product export strategy promoted by the Japanese government. Cultivar identification systems using DNA markers are an effective tool for protecting breeders' rights. Here, a novel target cultivar-specific identification system using the chromatographic printed array strip method was developed for eight prominent Japanese citrus cultivars. A polymorphic InDel fragment specific to each cultivar was explored through the screening of published citrus InDel markers and next-generation sequencing of retrotransposon libraries. The cultivar-specific DNA marker set for each cultivar comprised 1-3 polymorphic InDel fragments in combination with a PCR-positive DNA marker for the ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit gene. The DNA markers were detected within 3 hours from DNA extraction to the detection by the C-PAS4 membrane stick following multiplex PCR. The developed system is superior as a convenient, rapid, and cost-effective DNA diagnostic method during inspection. The proposed target cultivar-specific identification system is expected to serve as an efficient tool for the injunction of suspicious registered cultivars, contributing to the protection of breeders' rights.

Elucidation of genetic variation and population structure of melon genetic resources in the NARO Genebank, and construction of the World Melon Core Collection.

Numerous genetic resources of major crops have been introduced from around the world and deposited in Japanese National Agriculture and Food Research Organization (NARO) Genebank. Understanding their genetic variation and selecting a representative subset ("core collection") are essential for optimal management and efficient use of genetic resources. In this study, we conducted genotyping-by-sequencing (GBS) to characterize the genetic relationships and population structure in 755 accessions of melon genetic resources. The GBS identified 39,324 single-nucleotide polymorphisms (SNPs) that are distributed throughout the melon genome with high density (one SNP/10.6 kb). The phylogenetic relationships and population structure inferred using this SNP dataset are highly associated with the cytoplasm type and geographical origin. Our results strongly support the recent hypothesis that cultivated melon was established in Africa and India through multiple independent domestication events. Finally, we constructed a World Melon Core Collection that covers at least 82% of the genetic diversity and has a wide range of geographical origins and fruit morphology. The genome-wide SNP dataset, phylogenetic relationships, population structure, and the core collection provided in this study should largely contribute to genetic research, breeding, and genetic resource preservation in melon.

QTL analysis and GWAS of agronomic traits in sweetpotato (Ipomoea batatas L.) using genome wide SNPs.

While sweetpotato (Ipomoea batatas L.) improvement has generally been done by field-based selection, molecular genetic studies on traits of interest, i.e., molecular markers are needed for enhancing the breeding program of this world's 7th most important crop, as such markers facilitate marker-assisted selection. Here, we performed a combined approach of QTLs analyses and GWAS of storage root ß-carotene content (BC), dry-matter (DM) and starch content (SC) using the genetic linkage maps constructed with 5,952 and 5,640 SNPs obtained from F1 progenies between cultivars 'J-Red' and 'Choshu'. BC was negatively correlated with DM (r = -0.45) and SC (r = -0.51), while DM was positively correlated with SC (r = 0.94). In both parental maps, a total of five, two and five QTL regions on linkage groups 7 and 8 were associated with BC, DM and SC, respectively. In GWAS of BC, one strong signal (P = 1.04 × 10-9) was observed on linkage group 8, which co-located with one of the above QTL regions. The SNPs markers found here, particularly for ß-carotene, would be useful base resources for future marker-assisted selection program with this trait.

DNA markers based on retrotransposon insertion polymorphisms can detect short DNA fragments for strawberry cultivar identification.

In this study, DNA markers were developed for discrimination of strawberry (Fragaria × ananassa L.) cultivars based on retrotransposon insertion polymorphisms. We performed a comprehensive genomic search to identify retrotransposon insertion sites and subsequently selected one retrotransposon family, designated CL3, which provided reliable discrimination among strawberry cultivars. Through analyses of 75 strawberry cultivars, we developed eight cultivar-specific markers based on CL3 retrotransposon insertion sites. Used in combination with 10 additional polymorphic markers, we differentiated 35 strawberry cultivars commonly cultivated in Japan. In addition, we demonstrated that the retrotransposon-based markers were effective for PCR detection of DNA extracted from processed food materials, whereas a SSR marker was ineffective. These results indicated that the retrotransposon-based markers are useful for cultivar discrimination for processed food products, such as jams, in which DNA may be fragmented or degraded.

Genome-Wide Association Studies (GWAS) for Yield and Weevil Resistance in Sweet potato (Ipomoea batatas (L.) Lam).

KEY MESSAGE: We apply the GWAS to sweet potato genome, and identified the SNPs associated with yield and weevil resistance. The sweet potato (Ipomoea batatas (L.) Lam) is a highly heterozygous, outcrossing, polyploid species, which presents challenges for genetic analysis. Therefore, we considered that genome-wide association studies (GWAS) may be applied to the study of the sweet potato genome. The yield of two sweet potato varieties [Purple Sweet Lord (PSL) and 90IDN-47] was assessed at two locations (Kumamoto and Okinawa prefectures) in Japan in 2013 and the yield scores were used for GWAS. The results showed that there were several single nucleotide polymorphisms (SNP) above the significance thresholds in PSL; two peaks were detected in Kumamoto and Okinawa on the Ib03-3 and Ib01-4 linkage groups of PSL, respectively. As for 90IDN-47, one relatively high peak was detected in Kumamoto on the Ib13-8 linkage group. Interestingly, although high peaks above significance thresholds were detected in Kumamoto and Okinawa in PSL, the peaks were located in different linkage groups. This result suggests that the genetic regions controlling yield may change in response to environmental conditions. Additionally, we investigated the degree of weevil damage to the plants, which is the greatest problem in sweet potato cultivation in Okinawa. In this experiment, no SNPs were identified above the significance thresholds. However, one relatively high peak was found in the 90IDN-47 genotype, which showed resistance to weevils. On the other hand, one relatively high peak was also detected in the PSL genotype, which showed susceptibility to weevils. These results suggest that two regions could affect weevil resistance and may contain the gene(s) controlling weevil resistance.

Impairment of Lhca4, a subunit of LHCI, causes high accumulation of chlorophyll and the stay-green phenotype in rice.

Chlorophyll is an essential molecule for acquiring light energy during photosynthesis. Mutations that result in chlorophyll retention during leaf senescence are called 'stay-green' mutants. One of the several types of stay-green mutants, Type E, accumulates high levels of chlorophyll in the pre-senescent leaves, resulting in delayed yellowing. We isolated delayed yellowing1-1 (dye1-1), a rice mutant whose yellowing is delayed in the field. dye1-1 accumulated more chlorophyll than the wild-type in the pre-senescent and senescent leaves, but did not retain leaf functionality in the 'senescent green leaves', suggesting that dye1-1 is a Type E stay-green mutant. Positional cloning revealed that DYE1 encodes Lhca4, a subunit of the light-harvesting complex I (LHCI). In dye1-1, amino acid substitution occurs at the location of a highly conserved amino acid residue involved in pigment binding; indeed, a severely impaired structure of the PSI-LHCI super-complex in dye1-1 was observed in a blue native PAGE analysis. Nevertheless, the biomass and carbon assimilation rate of dye1-1 were comparable to those in the wild-type. Interestingly, Lhcb1, a trimeric LHCII protein, was highly accumulated in dye1-1, in the chlorophyll-protein complexes. The high accumulation of LHCII in the LHCI mutant dye1 suggests a novel functional interaction between LHCI and LHCII.

Genetic linkage analysis using DNA markers in sweetpotato.

Sweetpotato is one of the most important food crop species in the world, with more than 104,000,000 tons produced each year, and the breeding of superior cultivars with agronomically important traits, such as improved disease resistance, yield, and nutrient richness, is necessary, especially in developing countries. However, as a result of the crop's complex genomic architecture, which results from its hexaploidy (2n = 6× = 90), high heterozygosity, huge genome, and outcrossing nature, the analysis of genetic linkage in sweetpotato has been challenging. In addition, the lack of whole genome sequences or gene annotations for cultivated hexaploids has interrupted the validation of mapped QTL regions and gene cloning. In spite of these technical difficulties, linkage map construction and QTL mapping analysis have been reported. This review summarizes the results of these linkage analyses, which used SSR, AFLP, and retrotransposon-based molecular markers, and describes future directions for the genetic analysis and marker-assisted breeding of this important but genetically complex crop species.

Construction of a linkage map based on retrotransposon insertion polymorphisms in sweetpotato via high-throughput sequencing.

Sweetpotato (Ipomoea batatas L.) is an outcrossing hexaploid species with a large number of chromosomes (2n = 6x = 90). Although sweetpotato is one of the world's most important crops, genetic analysis of the species has been hindered by its genetic complexity combined with the lack of a whole genome sequence. In the present study, we constructed a genetic linkage map based on retrotransposon insertion polymorphisms using a mapping population derived from a cross between 'Purple Sweet Lord' (PSL) and '90IDN-47' cultivars. High-throughput sequencing and subsequent data analyses identified many Rtsp-1 retrotransposon insertion sites, and their allele dosages (simplex, duplex, triplex, or double-simplex) were determined based on segregation ratios in the mapping population. Using a pseudo-testcross strategy, 43 and 47 linkage groups were generated for PSL and 90IDN-47, respectively. Interestingly, most of these insertions (~90%) were present in a simplex manner, indicating their utility for linkage map construction in polyploid species. Additionally, our approach led to savings of time and labor for genotyping. Although the number of markers herein was insufficient for map-based cloning, our trial analysis exhibited the utility of retrotransposon-based markers for linkage map construction in sweetpotato.

Efficient screening of long terminal repeat retrotransposons that show high insertion polymorphism via high-throughput sequencing of the primer binding site.

Retrotransposons have been used frequently for the development of molecular markers by using their insertion polymorphisms among cultivars, because multiple copies of these elements are dispersed throughout the genome and inserted copies are inherited genetically. Although a large number of long terminal repeat (LTR) retrotransposon families exist in the higher eukaryotic genomes, the identification of families that show high insertion polymorphism has been challenging. Here, we performed an efficient screening of these retrotransposon families using an Illumina HiSeq2000 sequencing platform with comprehensive LTR library construction based on the primer binding site (PBS), which is located adjacent to the 5' LTR and has a motif that is universal and conserved among LTR retrotransposon families. The paired-end sequencing library of the fragments containing a large number of LTR sequences and their insertion sites was sequenced for seven strawberry (Fragaria × ananassa Duchesne) cultivars and one diploid wild species (Fragaria vesca L.). Among them, we screened 24 families with a "unique" insertion site that appeared only in one cultivar and not in any others, assuming that this type of insertion should have occurred quite recently. Finally, we confirmed experimentally the selected LTR families showed high insertion polymorphisms among closely related cultivars.

mPing: The bursting transposon.

Though transposable elements (TEs) have been considered as an efficient source of evolution, it has never been possible to test this hypothesis because most of TE insertions had occurred millions of years ago, or because currently active TEs have very few copies in a host genome. However, mPing, the first active DNA transposon in rice, was revealed to hold a key to answer this question. mPing has attained high copy numbers and still retained very high activity in a traditional rice strain, which enabled direct observation of behavior and impact of a bursting TE. A comprehensive analysis of mPing insertion sites has revealed it avoids exons but prefers promoter regions and thus moderately affects transcription of neighboring genes. Some of the mPing insertions have introduced possibly useful expression profile to adjacent genes that indicated TE's potential in de novo formation of gene regulatory network.

Surgical resection of an intraluminal tumor in the azygos vein with an unknown primary site causing superior vena cava syndrome.

Intraluminal tumor in the azygos vein is a rare disease that can cause superior vena cava (SVC) syndrome. Radiotherapy and endovascular stenting with or without chemotherapy are reported to have a high clinical success rate for the management of SVC syndrome with malignancy, but a poor survival rate. Here, we report a 69-year-old man who presented with swelling of the face and upper extremities, who was diagnosed with SVC syndrome caused by an intraluminal tumor in the azygos vein. Enhanced chest computed tomography revealed an intraluminal mass with a filling defect from the azygos vein to the SVC, with no extravascular extension or dissemination of the primary tumor. Surgical resection of the mass en bloc with the azygos vein and SVC reconstruction was performed. A poorly differentiated carcinoma was diagnosed on postoperative pathological evaluation. Twelve months after resection, the patient was well with no signs of recurrent disease. This case highlights that surgical resection should be considered as a treatment of choice for the management of SVC syndrome caused by an intraluminal malignancy in the azygos vein.

Pulmonary endarterectomy for chronic thromboembolic pulmonary hypertension with bronchial obstruction by a carcinoid tumor.

Pulmonary endarterectomy (PEA) is a standard treatment for chronic thromboembolic pulmonary hypertension (CTEPH). CTEPH combined with bronchial obstruction by a tumor is rare but should be assessed carefully because PEA for obstructed segments can be less therapeutic and make the subsequent surgical resection challenging. This report describes a case of CTEPH with bronchial obstruction by a typical carcinoid tumor in a 75-year-old man. On-site evaluation and removal of the obstructive tumor were performed bronchoscopically, increasing the effectiveness of subsequent PEA for all affected pulmonary segments. This report illustrates a PEA strategy to treat CTEPH with bronchial tumor obstruction.

Mapping of Nematode Resistance in Hexaploid Sweetpotato Using an Next-Generation Sequencing-Based Association Study.

The southern root-knot nematode (SRKN; Meloidogyne incognita) is a typical parasitic nematode that affects sweetpotato [Ipomoea batatas (L.) Lam.], causing a significant decrease in crop yield and commercial value. In Japan, the SRKN is classified into 10 races: SP1-SP5, SP6-1, SP6-2, and SP7-SP9, with the dominant race differing according to the cultivation area. Soil insecticides have previously been used to reduce the soil density of SRKNs; however, this practice is both costly and labor intensive. Therefore, the development of SRKN-resistant sweetpotato lines and cultivars is necessary. However, due to the complexity of polyploid inheritance and the highly heterogeneous genomic composition of sweetpotato, genetic information and research for this species are significantly lacking compared to those for other major diploid crop species. In this study, we utilized the recently developed genome-wide association approach, which uses multiple-dose markers to assess autopolyploid species. We performed an association analysis to investigate resistance toward SRKN-SP2, which is the major race in areas with high sweetpotato production in Japan. The segregation ratio of resistant and susceptible lines in the F1 mapping population derived from the resistant "J-Red" and susceptible "Choshu" cultivars was fitted to 1: 3, suggesting that resistance to SP2 may be regulated by two loci present in the simplex. By aligning the double digest restriction-site associated DNA sequencing reads to the published Ipomoea trifida reference sequence, 46,982 single nucleotide polymorphisms (SNPs) were identified (sequencing depth > 200). The association study yielded its highest peak on chromosome 7 (Chr07) and second highest peak on chromosome 3 (Chr03), presenting as a single-dose in both loci. Selective DNA markers were developed to screen for resistant plants using the SNPs identified on Chr03 and Chr07. Our results showed that SRKN-SP2-resistant plants were selected with a probability of approximately 70% when combining the two selective DNA markers. This study serves as a model for the identification of genomic regions that control agricultural traits and the elucidation of their effects, and is expected to greatly advance marker-assisted breeding and association studies in polyploid crop species.

Exploiting the miniature inverted-repeat transposable elements insertion polymorphisms as an efficient DNA marker system for genome analysis and evolutionary studies in wheat and related species.

Transposable elements (TEs) constitute ~80% of the complex bread wheat genome and contribute significantly to wheat evolution and environmental adaptation. We studied 52 TE insertion polymorphism markers to ascertain their efficiency as a robust DNA marker system for genetic studies in wheat and related species. Significant variation was found in miniature inverted-repeat transposable element (MITE) insertions in relation to ploidy with the highest number of "full site" insertions occurring in the hexaploids (32.6 ± 3.8), while the tetraploid and diploid progenitors had 22.3 ± 0.6 and 15.0 ± 3.5 "full sites," respectively, which suggested a recent rapid activation of these transposons after the formation of wheat. Constructed phylogenetic trees were consistent with the evolutionary history of these species which clustered mainly according to ploidy and genome types (SS, AA, DD, AABB, and AABBDD). The synthetic hexaploids sub-clustered near the tetraploid species from which they were re-synthesized. Preliminary genotyping in 104 recombinant inbred lines (RILs) showed predominantly 1:1 segregation for simplex markers, with four of these markers already integrated into our current DArT-and SNP-based linkage map. The MITE insertions also showed stability with no single excision observed. The MITE insertion site polymorphisms uncovered in this study are very promising as high-potential evolutionary markers for genomic studies in wheat.

Comprehensive survey of transposon mPing insertion sites and transcriptome analysis for identifying candidate genes controlling high protein content of rice.

Rice is the most important crop species in the world, being staple food of more than 80% of people in Asia. About 80% of rice grain is composed of carbohydrates (starch), with its protein content as low as 7-8%. Therefore, increasing the protein content of rice offers way to create a stable protein source that contributes to improving malnutrition and health problems worldwide. We detected two rice lines harboring a significantly higher protein content (namely, HP5-7 and HP7-5) in the EG4 population. The EG4 strain of rice is a unique material in that the transposon mPing has high transpositional activity and high copy numbers under natural conditions. Other research indicated that mPing is abundant in the gene-rich euchromatic regions, suggesting that mPing amplification should create new allelic variants, novel regulatory networks, and phenotypic changes in the EG4 population. Here, we aimed to identify the candidate genes and/or mPing insertion sites causing high protein content by comprehensively identifying the mPing insertion sites and carrying out an RNA-seq-based transcriptome analysis. By utilizing the next-generation sequencing (NGS)-based methods, ca. 570 mPing insertion sites were identified per line in the EG4 population. Our results also indicated that mPing apparently has a preference for inserting itself in the region near a gene, with 38 genes in total found to contain the mPing insertion in the HP lines, of which 21 and 17 genes were specific to HP5-7 and HP7-5, respectively. Transcriptome analysis revealed that most of the genes related to protein synthesis (encoding glutelin, prolamin, and globulin) were up-regulated in HP lines relative to the control line. Interestingly, the differentially expressed gene (DEG) analysis revealed that the expression levels of many genes related to photosynthesis decreased in both HP lines; this suggests the amount of starch may have decreased, indirectly contributing to the increased protein content. The high-protein lines studied here are expected to contribute to the development of high protein-content rice by introducing valuable phenotypic traits such as high and stable yield, disease resistance, and abundant nutrients.

Transcriptome Analysis Reveals Key Genes Involved in Weevil Resistance in the Hexaploid Sweetpotato.

Because weevils are the most damaging pests of sweetpotato, the development of cultivars resistant to weevil species is considered the most important aspect in sweetpotato breeding. However, the genes and the underlying molecular mechanisms related to weevil resistance are yet to be elucidated. In this study, we performed an RNA sequencing-based transcriptome analysis using the resistant Kyushu No. 166 (K166) and susceptible Tamayutaka cultivars. The weevil resistance test showed a significant difference between the two cultivars at 30 days after the inoculation, specifically in the weevil growth stage and the suppressed weevil pupation that was only observed in K166. Differential expression and gene ontology analyses revealed that the genes upregulated after inoculation in K166 were related to phosphorylation, metabolic, and cellular processes. Because the weevil resistance was considered to be related to the suppression of larval pupation, we investigated the juvenile hormone (JH)-related genes involved in the inhibition of insect metamorphosis. We found that the expression of some terpenoid-related genes, which are classified as plant-derived JHs, was significantly increased in K166. This is the first study involving a comprehensive gene expression analysis that provides new insights about the genes and mechanisms associated with weevil resistance in sweetpotato.

Genetic Mapping in Autohexaploid Sweet Potato with Low-Coverage NGS-Based Genotyping Data.

Next-generation sequencing (NGS)-based genotyping methods can generate numerous genetic markers in a single experiment and have contributed to plant genetic mapping. However, for high precision genetic analysis, the complicated genetic segregation mode in polyploid organisms requires high-coverage NGS data and elaborate analytical algorithms. In the present study, we propose a simple strategy for the genetic mapping of polyploids using low-coverage NGS data. The validity of the strategy was investigated using simulated data. Previous studies indicated that accurate allele dosage estimation from low-coverage NGS data (read depth < 40) is difficult. Therefore, we used allele dosage probabilities calculated from read counts in association analyses to detect loci associated with phenotypic variations. The allele dosage probabilities showed significant detection power, although higher allele dosage estimation accuracy resulted in higher detection power. On the contrary, differences in the segregation patterns between the marker and causal genes resulted in a drastic decrease in detection power even if the marker and casual genes were in complete linkage and the allele dosage estimation was accurate. These results indicated that the use of a larger number of markers is advantageous, even if the accuracy of allele dosage estimation is low. Finally, we applied the strategy for the genetic mapping of autohexaploid sweet potato (Ipomoea batatas) populations to detect loci associated with agronomic traits. Our strategy could constitute a cost-effective approach for preliminary experiments done performed to large-scale studies.

An in vitro model of temporal enhancement of epithelium barrier permeability by low-energy shock waves without contrast agents.

One of the commonly used techniques for drug delivery is to temporarily increase the permeability of tissue barriers. Acoustic energies such as ultrasound and shock waves are known to modulate tissue permeability. Recently, it was found that shock waves modulate the blood-brain barrier in the rat brain without injection of contrast agents such as microbubbles. This finding implies that modulation of other tissue barriers by shock wave exposure without contrast agents may be possible. To examine whether the modulation is also possible with other tissue barriers, we here investigated whether shock waves would modulate an in vitro tissue barrier model consisting of epithelial cells cultured on culture inserts. The permeability of the epithelium sheets evaluated by trans-epithelial electrical resistance (TEER) was increased following shock waves at a peak pressure of 11 MPa. The increased permeability recovered within 2 h. This enhancement was realized with one-shot low-energy shock waves having an acoustic energy of 0.013 mJ/mm2. Monitoring the peak pressures in every exposure revealed that the minimum peak pressure required for the enhancement is 2.9 MPa. These results indicate that shock wave exposure has the potential to temporarily increase the permeability of epithelium barriers to enhance drug delivery without contrast agents. Graphical abstract Enhancements of epithelial barrier permeability were evaluated by trans-epithelial electrical resistance (TEER) before and after shock wave exposures.