Ethics approval and consent to participate
The study was conducted in accordance with the Helsinki Declaration on this subject. On recruitment, participants signed a document in which they consented to the confidential treatment of their personal information and consented to the publication of anonymized data. The study protocol and the informed consent form were reviewed and approved by the Research Ethics Board of the Tokyo University of Agriculture (Authorized No. 2211).
Samples and collection
Scalp hair shaft samples were collected from 12 healthy Japanese and Taiwanese adults (six females and six males), aged 21 to 45 years, who consented to participate in this study. No volunteer was taking any medication during the experimental period. Hair shafts, hair swabs, and scalp swabs were collected using nitrile gloves, and hair shaft samples were cut using sterile scissors. After sampling, the hairs were chopped into lengths of ~5 mm with scissors and placed in plastic microtubes. Hair swab samples were directly taken from the hair using cotton swabs (Mentip for hospitals, Nihon Menbou Corporation, Saitama, Japan) and pre-moistened with 50 μL of sterile distilled water. Scalp swab samples were collected directly from the head crown using cotton swabs by rubbed them onto the scalp surface (between hair strands) to cover a total surface area of 2.5 cm2. The head of each swab was removed from the handle and placed in an Eppendorf tube.
Extraction of bacterial DNA
Bacterial DNA was extracted using the NucleoSpin® Tissue kit (MACHEREY-NAGEL, Düren, Germany) according to the manufacturer’s instructions, with a slight modification. First, hair shaft, hair swab, and scalp swab samples were immersed in 100 μL lysozyme solution (20 mg/ml lysozyme derived from egg white [Wako Pure Chemical Industries, Osaka, Japan] in 20 mM Tris-HCl and 0.2 mM EDTA, pH 8.0) for 30 min at 37 °C, as previously reported (34), and the DNA extracts obtained (100 μL) were stored at -20 °C until use.
Analysis of bacterial community structures using 16S rRNA gene sequencing
To analyze the bacterial community structures using the MiSeq™ platform (Illumina Inc., CA, USA), a three-step PCR method was performed using the extracted DNA samples. In the first-step PCR amplification, a universal primer set for the V4 region of the bacterial 16S rRNA gene (515F, 5ʹ-GTG CCA GCM GCC GCG GTA A-3ʹ and 806R, 5ʹ-GGA CTA CHV GGG TWT CTA AT-3ʹ) (48) was used. The 25 μL reaction mixture consisted of 12.5 μL of Kapa HiFi HotStart Ready Mix (Kapa Biosystems Inc., Wilmington, MA, USA), 0.5 μL of each primer (10 pM), and 11.5 μL of extracted bacterial DNA. The amplification program included an initial denaturation step at 95 °C for 3 min, followed by 40 cycles of denaturation at 98 °C for 30 s, annealing at 56 °C for 30 s, and elongation at 72 °C for 30 s. After electrophoresis through a 1.5 % (w/v) agarose gel, the targeted bands were excised from the gel with sterilized cutters, and the DNA extracted using the FastGene® Gel/PCR Extraction Kit (NIPPON Genetics Co., Ltd., Tokyo, Japan), according to the manufacturer’s instructions. The DNA concentration was measured using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). To prepare for 16S rRNA amplicon sequencing using MiSeq sequencing, templates are given tail, adapter, and index sequences in a two-step PCR; therefore, long-tailed primers are required, making amplification difficult. We were unable to perform direct amplification using a two-step PCR, probably because the amount of bacterial DNA obtained from the 3 cm hair shafts was very low. Therefore, we first performed a PCR using a universal primer set without any additional sequences. Consequently, we successfully obtained sufficient template fragments with a minimum number of reaction cycles.
For the second-step PCR, a universal primer set for the V4 region of the bacterial 16S rRNA gene and tailed sequences for MiSeq sequencing were used (1-515F, 5ʹ- TCG TCG GCA GCG TCA GAT GTG TAT AAG AGA CAG GTG CCA GCM GCC GCG GTA A-3ʹ and 1-806R, 5ʹ-GTC TCG TGG GCT CGG AGA TGT GTA TAA GAG ACA GGG ACT ACH VGG GTW TCT AAT-3ʹ) (49). Although it has been reported that this primer set poorly amplifies Propionibacterium in human skin (50), the results here showed good amplification of the predominant species in hair, Cutibacterium acnes (previously named Propionibacterium acnes). The 25 μL reaction mixture consisted of 1.0 µL of each primer (5 µM), which was heat-shocked at 95 °C for 5 min, 12.5 µL of Kapa HiFi HotStart Ready Mix, 12.5 ng of DNA obtained from the first-step PCR amplicon, and sterilized ultrapure water. The amplification program included an initial denaturation step at 95 °C for 3 min, followed by 20 cycles of denaturation at 98 °C for 30 s, annealing at 55 °C for 30 s, and elongation at 72 °C for 30 s. The PCR products were purified using the FastGene® Gel/PCR extraction kit according to the manufacturer’s instructions.
For the third-step PCR, a primer set with flow cell adapter sequences, index sequences, and tailed sequences was used (forward primer, 5ʹ-AAT GAT ACG GCG ACC ACC GAG ATC TAC AC-Index sequence-TCG TCG GCA GCG TC-3ʹ and reverse primer, 5ʹ-CAA GCA GAA GAC GGC ATA CGA GAT-Index sequence-GTC TCG TGG GCT CGG-3ʹ). The third-step PCR mixture (25 µL) comprised 12.5 µL of Kapa HiFi HotStart Ready Mix, 0.5 µL of each primer (10 pM), and 11.5 µL of the second-step PCR amplicon. The amplification program included an initial denaturation step at 95 °C for 3 min, followed by 8 cycles of denaturation at 98 °C for 30 s, annealing at 55 °C for 30 s, and elongation at 72 °C for 30 s. After electrophoresis in a 1.5 % (w/v) agarose gel, the target bands were excised with sterilized cutters, and the DNA extracted using the FastGene® Gel/PCR extraction kit, as described above. The DNA concentrations of the third-step PCR amplicons were quantified using the Qubit dsDNA HS assay kit (Thermo Fisher Scientific Inc., USA) according to the manufacturer’s instructions. Purified PCR products from each sample were mixed, denatured, and sequenced with a MiSeq System (Illumina Inc., USA) using the MiSeq reagent kit v3 (300 bp ×2 cycles with paired end reads; Illumina, Inc., USA) according to the manufacturer’s instructions. Good’s coverage values (>95 %) were obtained for all hair samples using the DNA extraction kit and PCR conditions described above. Good’s coverage values were estimated using QIIME™ 1.9.1 software (51).
Bioinformatics and statistical analysis
The index and universal sequences of each read were checked and reads with complete index sequences were selected as valid. USEARCH V8.1.1861 software (52) was used to merge paired-end reads and remove chimeric sequences. After the chimeric check, the reads were grouped into operational taxonomic units (OTUs) with >97 % similarity. Clostridium and Escherichia were filtered from the OTU tables for further analysis because they are frequently used in our laboratory and have rarely been detected in previous reports (34–36). Alpha diversity (observed OTUs and Shannon’s index) was evaluated at a 1 % OTU distance using the QIIME™ software package (51). For taxonomy-based analyses, the representative sequences of each OTU were analyzed using the EzBioCloud platform (53). Statistical analysis of the bacterial community structure was performed using the Kruskal–Wallis test in XLSTAT software ver. 2014 (http://www.xlstat.com/en/).
Data availability
Raw sequence data were deposited under the bioproject number PRJDB18517 and the biosample number SAMD00802426-00802533 in the DDBJ Sequence Read Archive. The top 15 most abundant OTU sequences were deposited in the DDBJ/ENA/GenBank databases under the accession numbers LC822820-LC822834.