Transcriptional landscape of small non-coding RNAs reveals diversity of categories and functions in molluscs.

Small non-coding RNAs (sncRNAs) are non-coding RNA molecules that play various roles in metazoans. Among the sncRNAs, microRNAs (miRNAs) guide post-translational gene regulation during cellular development, proliferation, apoptosis, and differentiation, while PIWI-interacting RNAs (piRNAs) suppress transposon activity to safeguard the genome from detrimental insertion mutagenesis. While an increasing number of piRNAs are being identified in the soma and germlines of various organisms, they are scarcely reported in molluscs. To unravel the small RNA (sRNA) expression patterns and genomic function in molluscs, we generated a comprehensive sRNA dataset by sRNA sequencing (sRNA-seq) of eight mollusc species. Abundant miRNAs were identified and characterized in all investigated molluscs, and ubiquitous piRNAs were discovered in both somatic and gonadal tissues in six of the investigated molluscs, which are more closely associated with transposon silencing. Tens of piRNA clusters were also identified based on the genomic mapping results, which varied among different tissues and species. Our dataset serves as important reference data for future genomic and genetic studies on sRNAs in these molluscs and related species, especially in elucidating the ancestral state of piRNAs in bilaterians.

Pilot study of a comprehensive resource estimation method from environmental DNA using universal D-loop amplification primers.

Many studies have investigated the ability of environmental DNA (eDNA) to identify the species. However, when individual species are to be identified, accurate estimation of their abundance using traditional eDNA analyses is still difficult. We previously developed a novel analytical method called HaCeD-Seq (haplotype count from eDNA by sequencing), which focuses on the mitochondrial D-loop sequence for eels and tuna. In this study, universal D-loop primers were designed to enable the comprehensive detection of multiple fish species by a single sequence. To sequence the full-length D-loop with high accuracy, we performed nanopore sequencing with unique molecular identifiers (UMI). In addition, to determine the D-loop reference sequence, whole genome sequencing was performed with thin coverage, and complete mitochondrial genomes were determined. We developed a UMI-based Nanopore D-loop sequencing analysis pipeline and released it as open-source software. We detected 5 out of 15 species (33%) and 10 haplotypes out of 35 individuals (29%) among the detected species. This study demonstrates the possibility of comprehensively obtaining information related to population size from eDNA. In the future, this method can be used to improve the accuracy of fish resource estimation, which is currently highly dependent on fishing catches.

Structure Determination of Kahalalide Analogues Based on Metagenomic Analysis of a Bryopsis sp. Marine Green Alga.

Two kahalalide analogues were isolated from a Bryopsis sp. marine green alga. Even though our initial structure determination of the peptides by NMR and MS identified them as kahalalide Z1 (KZ1; 3) and Z2 (KZ2; 4), the absolute configuration of the Thr residues by Marfey's analysis was different from those found in kahalalide F (KF), 3, and 4. To ascertain the absolute configuration of the amino acid residues genetically, we conducted a metagenomic analysis for symbiotic bacteria in the alga, leading to the biosynthetic gene cluster (BGC) responsible for producing the kahalalides named kahalalides Z3 (KZ3; 1) and Z4 (KZ4; 2). The identification of amino acid residues based on the A-domain suggested these peptides possess the amino acid sequence d-allo-Thr-l-Val-l-Val-d-Val residues at the N-terminus, instead of the d-Val-l-Thr-l-Val-d-Val residues found in KF, 3, and 4. The N-terminal amino acid sequence including absolute configuration was unambiguously determined by a comparison of LCMS data of synthetic tetrapeptides and the hydrolysates derived from 1 and 2. This structural difference is caused by swapping the substrate specificities of the first two A-domains.

Zebrafish Danio rerio myotomal muscle structure and growth from a spatial transcriptomics perspective.

Fish exhibit different muscle structures and growth characteristics compared with mammals. We used a spatial transcriptomics approach and examined myotomal muscle sections from zebrafish. Adult muscles were divided into eight regions according to spatial gene expression characteristics. Slow muscle was located in the wedge-shaped region near the lateral line and at the base of the dorsal fin, intermediate muscle was located in a ribbon-shaped region adjacent to slow muscle, and fast muscle was located in the deep region of the trunk, surrounded by intermediate muscle; the interior of fast muscle was further divided into 6 parts by their transcriptomic features. Combined analysis of adult and larval data revealed that adult muscles contain specific regions similar to larval muscles. These regions showed active myogenesis and a high expression of genes associated with muscle hyperplasia. This is the first study to apply spatial transcriptomics to fish myotomal muscle structure and growth.

Genetic characterization of 1210 Japanese pedigrees with inherited retinal diseases by whole-exome sequencing.

Inherited retinal diseases (IRDs) comprise a phenotypically and genetically heterogeneous group of ocular disorders that cause visual loss via progressive retinal degeneration. Here, we report the genetic characterization of 1210 IRD pedigrees enrolled through the Japan Eye Genetic Consortium and analyzed by whole exome sequencing. The most common phenotype was retinitis pigmentosa (RP, 43%), followed by macular dystrophy/cone- or cone-rod dystrophy (MD/CORD, 13%). In total, 67 causal genes were identified in 37% (448/1210) of the pedigrees. The first and second most frequently mutated genes were EYS and RP1, associated primarily with autosomal recessive (ar) RP, and RP and arMD/CORD, respectively. Examinations of variant frequency in total and by phenotype showed high accountability of a frequent EYS missense variant (c.2528G>A). In addition to the two known EYS founder mutations (c.4957dupA and c.8805C>G) of arRP, we observed a frequent RP1 variant (c.5797C>T) in patients with arMD/CORD.

Induction of endoplasmic reticulum stress markers in an acromegaly model.

Acromegaly is a growth hormone (GH) excess pathological condition in humans. Acromegaly is associated with somatic disfigurement and a wide range of systemic manifestations such as arthritis, neuropathy, carpal tunnel syndrome, reproductive disorders, metabolic disorders, and gastrointestinal complications. The influence of excess GH on the cellular level could aid in understanding the root causes of acromegaly-related health complications. Previously, we found that GH excess induces DNA damage to somatic cells and reduces the stem cells number and causes premature aging. In this study, an in-depth analysis of the acromegaly RNAseq data revealed the disruption of important biological cellular processes. Gene set enrichment analysis, heatmap, and enrichment analysis of acromegaly RNAseq data revealed induction of endoplasmic reticulum (ER) stress markers in various organs. Interestingly, the induction of ER stress was even more apparent than in aged zebrafish. Splicing of box-binding protein-1 (XBP1) mRNA is a hallmark of ER stress. Therefore, we quantified spliced XBP1 mRNA in different organs of our acromegaly model. Thus, our study emphasizes the importance of ER stress in GH oversecretion, which is important for understanding the health complications of acromegaly.

Novel mutations in malonyl-CoA-acyl carrier protein transacylase provoke autosomal recessive optic neuropathy.

Inherited optic neuropathies are rare eye diseases of optic nerve dysfunction that present in various genetic forms. Previously, mutation in three genes encoding mitochondrial proteins has been implicated in autosomal recessive forms of optic atrophy that involve progressive degeneration of optic nerve and retinal ganglion cells (RGC). Using whole exome analysis, a novel double homozygous mutation p.L81R and pR212W in malonyl CoA-acyl carrier protein transacylase (MCAT), a mitochondrial protein involved in fatty acid biosynthesis, has now been identified as responsible for an autosomal recessive optic neuropathy from a Chinese consanguineous family. MCAT is expressed in RGC that are rich in mitochondria. The disease variants lead to structurally unstable MCAT protein with significantly reduced intracellular expression. RGC-specific knockdown of Mcat in mice, lead to an attenuated retinal neurofiber layer, that resembles the phenotype of optic neuropathy. These results indicated that MCAT plays an essential role in mitochondrial function and maintenance of RGC axons, while novel MCAT p.L81R and p.R212W mutations can lead to optic neuropathy.

Environmental DNA study on aquatic ecosystem monitoring and management: Recent advances and prospects.

Environmental DNA (eDNA) is organismal DNA that can be detected in the environment and is derived from cellular material of organisms shed into aquatic or terrestrial environments. It can be sampled and monitored using molecular methods, which is important for the early detection of invasive and native species as well as the discovery of rare and cryptic species. While few reviews have summarized the latest findings on eDNA for most aquatic animal categories in the aquatic ecosystem, especially for aquatic eDNA processing and application. In the present review, we first performed a bibliometric network analysis of eDNA studies on aquatic animals. Subsequently, we summarized the abiotic and biotic factors affecting aquatic eDNA occurrence. We also systematically discussed the relevant experiments and analyses of aquatic eDNA from various aquatic organisms, including fish, molluscans, crustaceans, amphibians, and reptiles. Subsequently, we discussed the major achievements of eDNA application in studies on the aquatic ecosystem and environment. The application of eDNA will provide an entirely new paradigm for biodiversity conservation, environment monitoring, and aquatic species management at a global scale.

Discovery and functional understanding of MiRNAs in molluscs: a genome-wide profiling approach.

Small non-coding RNAs play a pivotal role in gene regulation, repression of transposable element and viral activity in various organisms. Among the various categories of these small non-coding RNAs, microRNAs (miRNAs) guide post-translational gene regulation in cellular development, proliferation, apoptosis, oncogenesis, and differentiation. Here, we performed a genome-wide computational prediction of miRNAs to improve the understanding of miRNA observation and function in molluscs. As an initial step, hundreds of conserved miRNAs were predicted in 35 species of molluscs through genome scanning. Afterwards, the miRNAs' population, isoforms, organization, and function were characterized in detail. Furthermore, the key miRNA biogenesis factors, including AGO2, DGCR8, DICER, DROSHA, TRABP2, RAN, and XPO5, were elucidated based on homologue sequence searching. We also summarized the miRNAs' function in biomineralization, immune and stress response, as well as growth and development in molluscs. Because miRNAs play a vital role in various lifeforms, this study will provide insight into miRNA biogenesis and function in molluscs, as well as other invertebrates.

A new PDE6A missense variant p.Arg544Gln in rod-cone dystrophy.

PURPOSE: Thus far, only one Japanese patient with autosomal recessive rod-cone dystrophy (AR-RCD) associated with the phosphodiesterase 6A gene (PDE6A) has been reported. The purpose of this study was to analyze the clinical features of a Japanese female patient with AR-RCD with a novel missense variant in PDE6A. METHODS: We performed whole-exome sequencing (WES) to identify the disease-causing variant and a comprehensive ophthalmic examination including full-field electroretinography (ERG). RESULTS: WES analysis revealed that the patient carried a novel homozygous missense variant (c.1631G > A; p.Arg544Gln) in PDE6A. Her unaffected parents carried the heterozygous variant. The patient reported night blindness in her early 20 s. At the age of 25 years, she underwent a comprehensive ophthalmic examination. Her corrected visual acuity was 20/13 in the right and 20/10 in the left eyes. Fundus images showed degenerative changes with bone spicule pigmentation in the mid-peripheral retina, and peripheral retinal vessels were not attenuated. Ultra-wide-field fundus autofluorescence images demonstrated large hypoautofluorescent regions corresponding to the degenerative changes, surrounded by hyperautofluorescence. Cross-sectional optical coherence tomography demonstrated a preserved ellipsoid zone and retinal thickness in the center of the macula, with perifoveal atrophy. ERG responses were subnormal, revealing that rod-mediated responses were more affected than cone-mediated responses, consistent with findings observed in RCD. CONCLUSIONS: This is the second case of a patient with AR-RCD associated with PDE6A in the Japanese population. These findings will contribute to a better clinical understanding of PDE6A-associated RCD and valuable insights for gene therapy trials.

Francisellosis of Yesso scallops Mizuhopecten yessoensis in Japan is caused by a novel type of Francisella halioticida.

Francisella halioticida, the causative agent of francisellosis of the giant abalone Haliotis gigantea, has also been isolated from Yesso scallops Mizuhopecten yessoensis, which presented with orange/pinkish lesions in the adductor muscle and experienced high mortality. However, it is not clear whether the F. halioticida isolated from the giant abalone and Yesso scallops are phenotypically and genetically identical to each other. The present study revealed that isolates from the giant abalone and Yesso scallops were phenotypically different, with slower growth in modified eugon broth and a lack of prolyl aminopeptidase and phenylalanine aminopeptidase in Yesso scallop isolates. Additionally, we found that 3 of 8 housekeeping genes were different between them. Based on these phenotypic and genetic differences, we propose that F. halioticida isolated from Yesso scallops in Japan be designated as the 'J-scallop type' to distinguish it from strains from abalone ('abalone type'). Whole-genome sequencing analysis of a strain belonging to the J-scallop type showed that the overall similarity between the J-scallop and abalone type strains was estimated to be 99.84%. In accordance with a lack of prolyl aminopeptidase activity, in general, all of the J-scallop type strains examined have a 1 bp deletion in the responsible gene encoding prolyl aminopeptidase. This deletion was confirmed in all F. halioticida in diseased Yesso scallops examined, suggesting that in Japan, francisellosis of Yesso scallops is caused by a novel type of F. halioticida and not by the abalone type.

A Review of Discovery Profiling of PIWI-Interacting RNAs and Their Diverse Functions in Metazoans.

PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs (sncRNAs) that perform crucial biological functions in metazoans and defend against transposable elements (TEs) in germ lines. Recently, ubiquitously expressed piRNAs were discovered in soma and germ lines using small RNA sequencing (sRNA-seq) in humans and animals, providing new insights into the diverse functions of piRNAs. However, the role of piRNAs has not yet been fully elucidated, and sRNA-seq studies continue to reveal different piRNA activities in the genome. In this review, we summarize a set of simplified processes for piRNA analysis in order to provide a useful guide for researchers to perform piRNA research suitable for their study objectives. These processes can help expand the functional research on piRNAs from previously reported sRNA-seq results in metazoans. Ubiquitously expressed piRNAs have been discovered in the soma and germ lines in Annelida, Cnidaria, Echinodermata, Crustacea, Arthropoda, and Mollusca, but they are limited to germ lines in Chordata. The roles of piRNAs in TE silencing, gene expression regulation, epigenetic regulation, embryonic development, immune response, and associated diseases will continue to be discovered via sRNA-seq.

Clinical and genetic characteristics of 10 Japanese patients with PROM1-associated retinal disorder: A report of the phenotype spectrum and a literature review in the Japanese population.

Variants in the PROM1 gene are associated with cone (-rod) dystrophy, macular dystrophy, and other phenotypes. We describe the clinical and genetic characteristics of 10 patients from eight Japanese families with PROM1-associated retinal disorder (PROM1-RD) in a nationwide cohort. A literature review of PROM1-RD in the Japanese population was also performed. The median age at onset/examination of 10 patients was 31.0 (range, 10-45)/44.5 (22-73) years. All 10 patients showed atrophic macular changes. Seven patients (70.0%) had spared fovea to various degrees, approximately half of whom had maintained visual acuity. Generalized cone (-rod) dysfunction was demonstrated in all nine subjects with available electrophysiological data. Three PROM1 variants were identified in this study: one recurrent disease-causing variant (p.Arg373Cys), one novel putative disease-causing variant (p.Cys112Arg), and one novel variant of uncertain significance (VUS; p.Gly53Asp). Characteristic features of macular atrophy with generalized cone-dominated retinal dysfunction were shared among all 10 subjects with PROM1-RD, and the presence of foveal sparing was crucial in maintaining visual acuity. Together with the three previously reported variants [p.R373C, c.1551+1G>A (pathogenic), p.Asn580His (likely benign)] in the literature of Japanese patients, one prevalent missense variant (p.Arg373Cys, 6/9 families, 66.7%) detected in multiple studies was determined in the Japanese population, which was also frequently detected in the European population.

Clinical Course and Electron Microscopic Findings in Lymphocytes of Patients with DRAM2-Associated Retinopathy.

DRAM2-associated retinopathy is a rare inherited retinal dystrophy, and its outcome has not been determined. A single retinal involvement by a mutation of the DRAM2 gene is unexplained. We found three unrelated patients with a disease-causing DRAM2 variant in a biallelic state from 1555 Japanese individuals of 1314 families with inherited retinal dystrophy. We reviewed their medical records and examined their peripheral lymphocytes by transmission electron microscopy (TEM). Patient 1 was a 38-year-old woman who complained of night blindness and reduced vision. She developed macular degeneration at age 43 years. Patients 2 and 3 were a man and a woman both of whom noticed night blindness in their 30s. Both had a degeneration in the macula and midperiphery in their 40s, which progressed to a diffuse retinal degeneration in their 60s when their vision was reduced to hand motions. Three novel DRAM2 variants were identified. TEM of the lymphocytes of Patients 1 and 2 showed abnormal structures in 40.6% and 0.3% of the peripheral lymphocytes, respectively. We concluded that the DRAM2-associated retinopathy of our patients was a progressive rod-cone dystrophy, and the visual outcome was poor. The systemic effect of DRAM2 mutations may be compensable and have variations.

Gene expression profiles at different stages for formation of pearl sac and pearl in the pearl oyster Pinctada fucata.

BACKGROUND: The most critical step in the pearl formation during aquaculture is issued to the proliferation and differentiation of outer epithelial cells of mantle graft into pearl sac. This pearl sac secretes various matrix proteins to produce pearls by a complex physiological process which has not been well-understood yet. Here, we aimed to unravel the genes involved in the development of pearl sac and pearl, and the sequential expression patterns of different shell matrix proteins secreted from the pearl sac during pearl formation by pearl oyster Pinctada fucata using high-throughput transcriptome profiling. RESULTS: Principal component analysis (PCA) showed clearly different gene expression profiles between earlier (before 1 week) and later stages (1 week to 3 months) of grafting. Immune-related genes were highly expressed between 0 h - 24 h (donor dependent) and 48 h - 1 w (host dependent), and in the course of wound healing process pearl sac was developed by two weeks of graft transplantation. Moreover, for the first time, we identified some stem cell marker genes including ABCG2, SOX2, MEF2A, HES1, MET, NRP1, ESR1, STAT6, PAX2, FZD1 and PROM1 that were expressed differentially during the formation of pearl sac. The expression profiling of 192 biomineralization-related genes demonstrated that most of the shell matrix proteins (SMPs) involved in prismatic layer formation were first up-regulated and then gradually down-regulated indicating their involvement in the development of pearl sac and the onset of pearl mineralization. Most of the nacreous layer forming SMPs were up-regulated at 2 weeks after the maturation of pearl sac. Nacrein, MSI7 and shematrin involved in both layer formation were highly expressed during 0 h - 24 h, down-regulated up to 1 week and then up-regulated again after accomplishment of pearl sac formation. CONCLUSIONS: Using an RNA-seq approach we unraveled the expression pattern of the key genes involved in the development of pearl sac and pearl as a result of host immune response after grafting. These findings provide valuable information in understanding the molecular mechanism of pearl formation and immune response in P. fucata.

Late-onset night blindness with peripheral flecks accompanied by progressive trickle-like macular degeneration.

PURPOSE: To report the clinical and genetic characteristics of 6 cases with late-onset night blindness with peripheral flecks accompanied by progressive trickle-like macular degeneration. METHODS: Clinical and genetic data were collected from 6 independent patients who complained of night blindness in their fifth to eighth decade of life. The ophthalmological examinations included ophthalmoscopy, fundus autofluorescence (FAF), and full-field electroretinography (ERG). Whole exome sequencing with target gene analysis was performed to determine the causative genes and variants. RESULTS: All of the patients first complained of night blindness at the ages of 40-71 years. Funduscopic examinations demonstrated white or atrophic flecks scattered in the posterior pole and peripheral retina bilaterally. FAF showed patchy hypo-autofluorescence spots in the posterior pole similar to that of the trickling type of age-related macular degeneration (AMD). The region of abnormal FAF rapidly expanded with age, and one eye developed a choroidal neovascularization. The full-field scotopic ERGs with 20 min of dark adaptation were severely reduced or extinguished in all cases. There was partial recovery of the ERGs after 180 min of dark adaptation. The cone ERGs were reduced in all cases. Whole exome sequencing revealed no pathogenic variants of 301 retinal disease-associated genes. CONCLUSIONS: The six cases had some common features with the flecked retina syndrome, familial drusen, and late-onset retinal degeneration although none had pathogenic variants causative for these disorders. These cases may represent a subset of severe trickling AMD or a new clinical entity of acquired pan-retinal visual cycle deficiency of unknown etiology.

Novel biallelic loss-of-function KCNV2 variants in cone dystrophy with supernormal rod responses.

PURPOSE: To report clinical and genetic features including long-term full-field electroretinography (FF-ERG) findings of a patient with cone dystrophy with supernormal rod responses (CDSRR). METHODS: Ophthalmological medical records including FF-ERG were retrospectively reviewed. Genetic analysis using whole-exome sequencing (WES) was performed. Identified KCNV2 variants were confirmed by Sanger sequencing. RESULTS: A 30-year-old female patient was referred to our hospital for assessment of decreased vision from childhood. Funduscopy showed macular atrophy in both eyes. FF-ERG showed decreased amplitudes and delayed peak time of b-waves for dark-adapted (DA) 0.01 ERG, increased b/a-wave ratio with a slightly diminished a-wave for DA 3.0 and DA 25.7 ERG, residual a-waves and almost extinguished b-waves for light-adapted (LA) 3.0 ERG, and extremely diminished amplitudes in LA 30-Hz flicker responses. At 45 years of age, funduscopy showed progressive macular atrophy, whereas the responses for her FF-ERG remained unchanged compared to those observed at 30 years of age. WES identified the compound heterozygous KCNV2 variants (p.W67X and p.D174GfsX198) in the patient. These variants have previously been unreported as pathogenic variants. Each parent had one of the variants. Subsequently, the patient was finally diagnosed with CDSRR with the novel compound heterozygous KCNV2 variants. CONCLUSIONS: Biallelic loss-of-function KCNV2 variants (p.W67X and p.D174GfsX198) were identified as the cause of CDSRR. Long-term FF-ERG findings demonstrated there were no ERG changes during 15 years of observation, indicating that there was no evidence of progressive peripheral retinal dysfunction, in spite of worsening macular atrophy.

LRRTM4-C538Y novel gene mutation is associated with hereditary macular degeneration with novel dysfunction of ON-type bipolar cells.

The macula is a unique structure in higher primates, where cone and rod photoreceptors show highest density in the fovea and the surrounding area, respectively. The hereditary macular dystrophies represent a heterozygous group of rare disorders characterized by central visual loss and atrophy of the macula and surrounding retina. Here we report an atypical absence of ON-type bipolar cell response in a Japanese patient with autosomal dominant macular dystrophy (adMD). To identify a causal genetic mutation for the adMD, we performed whole-exome sequencing (WES) on four affected and four-non affected members of the family for three generations, and identified a novel p.C538Y mutation in a post-synaptic gene, LRRTM4. WES analysis revealed seven rare genetic variations in patients. We further referred to our in-house WES data from 1360 families with inherited retinal diseases, and found that only p.C538Y mutation in LRRTM4 was associated with adMD-affected patients. Combinatorial filtration using public database of single-nucleotide polymorphism frequency and genotype-phenotype annotated database identified novel mutation in atypical adMD.

Early onset flecked retinal dystrophy associated with new compound heterozygous RPE65 variants.

Purpose: To report genetic and clinical features of two unrelated Japanese patients with early onset flecked retinal dystrophy. Methods: Patients underwent comprehensive ophthalmic examinations that included electroretinography (ERG) after 30 min and 24 h of dark adaptation (DA). Disease-causing gene variants were identified with whole exome sequencing (WES), with identified candidates confirmed with direct sequencing. Results: WES identified compound heterozygous RPE65 variants in both patients. Variants in patient 1 included c.1543C>T (p.R515W) and c.683A>C (p.Q228P), while patient 2 exhibited c.1028T>A (p.L343*) and c.683A>C (p.Q228P). Although variants p.R515W and p.L343* have been previously reported as pathogenic, variant p.Q228P was reported as uncertain significance. Each unaffected parent carried the variant heterozygously. Both patients had similar ophthalmic findings, including decreased visual acuity with early onset night blindness, numerous dense white dots/flecks occurring mainly outside the vascular arcades, a diffuse and/or disrupted ellipsoid line as shown with optical coherence tomography, and non-recordable rod and combined responses along with decreased cone responses after 30 min of DA. After 24 h of DA, both patients exhibited marked or partial recovery of the combined responses. Conclusions: The results indicate that the recovery of combined or residual cone responses might be associated with a mild form of RPE65-related early onset flecked retinal dystrophy with new compound heterozygous variants.

Collaborative environmental DNA sampling from petal surfaces of flowering cherry Cerasus × yedoensis 'Somei-yoshino' across the Japanese archipelago.

Recent studies have shown that environmental DNA is found almost everywhere. Flower petal surfaces are an attractive tissue to use for investigation of the dispersal of environmental DNA in nature as they are isolated from the external environment until the bud opens and only then can the petal surface accumulate environmental DNA. Here, we performed a crowdsourced experiment, the "Ohanami Project", to obtain environmental DNA samples from petal surfaces of Cerasus × yedoensis 'Somei-yoshino' across the Japanese archipelago during spring 2015. C. × yedoensis is the most popular garden cherry species in Japan and clones of this cultivar bloom simultaneously every spring. Data collection spanned almost every prefecture and totaled 577 DNA samples from 149 collaborators. Preliminary amplicon-sequencing analysis showed the rapid attachment of environmental DNA onto the petal surfaces. Notably, we found DNA of other common plant species in samples obtained from a wide distribution; this DNA likely originated from the pollen of the Japanese cedar. Our analysis supports our belief that petal surfaces after blossoming are a promising target to reveal the dynamics of environmental DNA in nature. The success of our experiment also shows that crowdsourced environmental DNA analyses have considerable value in ecological studies.