Effect of homeopathic medicines on growth, survival and gastrointestinal microbiota of juvenile scallop Argopecten ventricosus

Objective. To study the effect of homeopathic medicines on growth, survival and gastrointestinal (GIT) microbiota of Catarina scallop Argopecten ventricosus. Materials and methods. Five homeopathic (HOM) treatments derived from bacteria [(T1) ViP-ViA 1D, (T2) ViP-ViA 7C], minerals [(T3) AcF-MsS 1D, (T4) PhASiT 7C] or venoms [T5) ViT 31C] and three controls: [(C1) diluted ethanol 1:99, (C2) diluted/succussed ethanol 1C and (C3) distilled water] were evaluated (21 days) in triplicate. Microbiota was analysed by sequencing the V3-V5 region of the 16S rRNA genes. Results. The best growth in shell-length corresponded to T1 (117 μm d-1) and T2 (108 μm d-1) and the highest survival to T3 and T5, stating T3 as the best HOMtreatment. A clear separation was found in rarefaction curves of HOM-treated against un-treated control scallops. Significant differences (p≤0.05) were found for Phyla (Proteobacteria> Actinobacteria> Firmicutes> Bacterloidetes>Chloroflexi and for Genera: Symbiobacterium> Microbacterium> Methylobacillus> Bacillus> Paenibacillus> Burkholderia> Nostoc> Methylobacterium> Leucobacter). The genus Symbiobacterium was dominant in T5, finding significant differences (p≤0.05) with respect to all treatments. At species level, Microbacterium maritypicum (Actinobacteria) showed a greater relative abundance (p≤0.05) in T1 and T3 and Symbiobacterium toebii (Firmicutes) was also significantly higher (p≤0.05) in abundance in T5 and T2, both against initial T0. Conclusions. This study showed for a first time, the composition of GIT microbiota in A. ventricosus and focused on the potential applicability of homeopathy to improve overall performance and modulate the GIT microbiota of the species.


INTRODUCTION
Commercially important scallops are produced along the Pacific coast and the Gulf of California, Mexico, but scallop fishing is considered unsustainable (1). Catarina scallop Argopecten ventricosus (Sowerby II, 1842) is a highly valued and rapid growth cultivable species (2,3) but hatchery spat production (4) is affected by massive mortalities during first growth phases (5,6). To prevent/ reduce disease and avoid massive larvae and spat mortalities in hatchery, antibiotics and other chemotherapeutics are traditionally used, which have different side effects on the environment (7). They induce development and spread of antibiotic resistant strains (8,9), accumulate in tissues of treated organisms (10) and affect their gastrointestinal tract (GIT) microbiota (11).
The structure and activity of GIT has effect on the host and helps to keep equilibrium with the surrounding environment. Different strategies have been developed exploiting this interaction to modulate the GIT microbial community, in order to promote growth, health and enhanced aquaculture productivity (6, 11,12). This is insufficient, so more efficacious and ecosustainable approaches such as aquacultural homeopathy are urgently required (13-17), which are mineral, plant and animal derived products that exert their action at ultra-low doses according to the "similarity principle", increasing host natural immunity and resistance (18). Homeopathic treatments may also reduce or substitute the use of antibiotics, hormones, disinfectants, and non-ecofriendly substances; they also decrease production costs and increase quality and innocuity of farmed molluscs, fishes and crustaceans (13-16).
This study focused on evaluate the effect of homeopathic medicines on growth, survival, and modulation of the GIT microbiota in A. ventricosus.

MATERIALS AND METHODS
Organisms and maintenance. Juvenile A. ventricosus (1.98 ± 0.1 cm of average length) from Bahía de La Paz, Baja California Sur, México, were cleaned of epifauna with filtered and UV-sterilized seawater after arriving to laboratory and were acclimatised one week (22°C, 37 PSU salinity). During the bioassay, scallops were allocated in nursery upwellers with filtered and UV sterilised seawater and fed on microalgae Isochrysis galbana and Chaetoceros calcitrans (1:1) at 10.8% of dry weight of soft tissues/day (aproximatelly 2x10 8 cells org -1 d -1 ).
Homeopathic components of T1 and T2 consisted in decimal or centesimal dynamisations (dilution/ succussion), developed at CIBNOR from nosodic mother tinctures (MTs) of pathogenic and virulent strains of Vibrio parahaemolyticus (CAIM-170, ViP) and Vibrio alginolyticus (CAIM-57 ViA), related with high mortalities in bivalves and shrimp. Each MT was obtained from bacterial cultures at 1 x 10 8 CFU mL -1 . Briefly, and in accordance to Mazón-Suástegui et al (16), the following procedures were applied: a) Bacterial cells were centrifuged (8000 g, 4°C, 20 min) and harvested from 15 mL of bacterial culture and washed twice; b) The corresponding pellets were diluted in 7.5 mL of MilliQ water and then inactivated through three freezing-thawing cycles (-80°C and 24°C, respectively) and sonicated eight times at 30 s, each time, to disrupt the cell wall and intracellular organelles; c) Unbroken cells were removed by centrifugation (3000 g, 4°C, 20 min) and the supernatant was diluted (1:1 v/v) in ethanol 87º GL (Similia ® Mexico) and vortexed at 3200 rpm (BenchMixer ® , Edison, NJ, U.S.A.) for two min to obtain 15 mL of MT from each bacterial strain.
The "stock" dynamisations (ViP 6C and ViPA 6C) were obtained by serial dilution and succussion in ethanol (1:99) of the respective TM. The "experimental" dynamisations (ViP 7C and ViA 7C) were obtained from respective 6C but substituting ethanol with distilled water.
The homeopathic treatment T3 (AcF 1D + MsS 1D) was the mixture of two decimal (1:9) dynamisations of 86% analytical grade phosphoric acid (Faga Lab ® , Guamuchil, Sinaloa, Mexico) and a saturated solution of sodium metasilicate (Faga Lab ® , Guamuchil, Sinaloa, Mexico). AcF´s TM was defined as the commercial product without dilution. MsS´s TM was prepared by diluting 9.4 g in 45 mL of distilled water, at 25°C. The first "experimental" decimal dynamisation (1D) of MsS was obtained by diluting 5 mL of MsS in 45 mL of distilled water (1:9) and vortexing two minutes at 3200 rpm (Benchmark mixer ® , Benchmark Scientific Inc.). The "experimental" AcF 1D was obtained by the same procedure.
Experimental design. We used 24 experimental units with 36 L of filtered and sterilized seawater and a continuous recirculating upflow pattern (1.68 mL s-1) during 21 h d -1 . Each unit hatd 4 PVC cylinders with plastic mesh bottom, allocating 13 scallops cylinder -1 (52 scallops unit -1 ). Temperature and salinity were maintained at 23.5 ± 0.5°C and 37.5±0.5UPS, respectively. Organic waste was removed daily and siphoning up to 60% of the water to reduce environmental variations. The scallops were fed as previously described. The treatments were added directly to seawater for 21 days (100 μL L -1 ) stopping water and food flow for three hours to favour their absorption in epithelial and branchial tissues. At the end of the trial, growth and survival were recorded for each treatment and its three replicates, taking samples for GIT metagenomic studies.
Organic waste was removed daily by siphoning up to 60% of water volume to reduce environmental variations Microbiota analysis. The bacterial DNA was extracted a two-time point: 24 samples at the beginning (T0) and 192 at the end (T21) of bioassay (8 scallops replicate -1 ; 24 scallops treatment -1 ). The GIT of each scallop was fixed individually (100 mg tissue in 500 µL RNA-Later ® ) and preserved at -20°C.

D N A E x t r a c t i o n a n d I l l u m i n a M i S e q sequencing.
Tissue from the GIT of each scallop (216 samples) were homogenised using the homogeniser FastPrep, (MP Biomedicals ® ) at high speed (6 m s -1 ) for 30 s in 500 μL of buffer TE (50:50 Tris: EDTA, pH 8), followed by the addition of 20 μL of 10% SDS and incubation at 56°C for 30 min. Then, 250 μL of 7.5 M potassium acetate was added to lyse cells, and the mixture was incubated on ice for 15 min.
After centrifuging cell lysate at 8000 g at 4°C for 10 min, the water phase was extracted with phenol: chloroform: isoamyl alcohol (25:24:1) and centrifuged. The resultant water phase was extracted with chloroform: isoamyl alcohol (24:1) and centrifuged. Afterward, 350 μL of cold isopropanol was added to the supernatant, followed by overnight incubation at -20°C and DNA collection through centrifugation. The DNA pellet was washed twice with cold 70 o GL ethanol, centrifuged, dried on air and resuspended in 20 μL of nuclease free water (22). DNA pools of each replica were made after DNA extraction for a total of 27 samples. DNA concentration and quality were determined at A260 nm and A280 nm through NanoDrop spectrophotometer (NanoDrop 8000, Thermo Scientific ® ) and electrophoresis in 1.5% agarose gel using the Sybr Safe ® (Invitrogen) indicator. The total DNA of each individual sample, which was diluted in nuclease free water to a final concentration of 100 ng µL-1, was used to amplify ±700 bp of the hypervariable region V3-V5 of the 16S rRNA gen (positions 339 and 939 from Escherichia coli).

Mazón-Suástegui et al -Effect of homeopathy in scallop microbiota
Phylogenetic analysis. Sequences were processed using the classifier software RDP (Ribosomal Database Project) with an average sequencing quality (Phred quality score) > 30 and a similarity threshold (Confidence Cutoff) of 80%, able to classify sequences at genus level. The Greengenes 16S rRNA Database (http:// greengenes.lbl.gov "Greengenes 16S") was used for taxonomic classification at phylum, class, order, family, and genus, and only in some cases even at species level. Bioinformatics analysis was performed by LANGEBIO-CINVESTAV.
The software Estimates for Statistical analysis V.9.1.0 (25) was used to calculate the Shannon-Weaver (H) and Simpson (J) indices of diversity as well as Chao richness estimates and to generate rarefaction curves. The software XLSTAT version 2016.05.33324 (Addinsoft 1995-2016) was used to compare microbial communities and treatments through the Principal Component Analysis (PCA) (26) and to detect differences on phylotype abundance and treatment using oneway analysis of variance (ANOVA) and Tukey´s multiple comparison test. Only were those principal components with high values (> 1.0) considered statistically significant.
Data accessibility. Sequences generated were uploaded to the NCBI Sequence Read Archive (SRA) with the key: Bioproject PRJNA341370, and the SRA access number: SRP089926.
Analysis of alfa diversity. PCR-amplified 16S-rRNA gene sequencing of microbiota associated to A. ventricosus GIT generated 21,717,584 reads with 80% classification from 27 amplicons. The number of OTUs was shown to be homogenous and had high richness with a number of OTUs/ Chao-1 around 699/779 and 703/768 for treated and control groups, respectively, whereas the initial control group at the beginning of the experiment (T0), exhibited the highest value of 689/985 ( Table  1). The Shannon-Weaver (H′) and Simpson (1-λ) indices also showed a high and homogeneous diversity for all experimental groups. ventricosus treated with homeopathic (HOM) medicines for 21 days. Values are mean ± confidence intervals at 95%. Different letter denotes significant differences (p<0.05). Although rarefaction curves did not reach the asymptotic phase, a clear separation and differentiation was detected between the rarefaction curves corresponding to HOM-treated groups and un-treated groups. As an exception, the rarefaction curve of the low (decimal) dynamisation HOM-treatment T1 (ViP 1D + ViA 1D) grouped together with T5 and control treatments ( Figure  2). The groups C1, C2, C3, T1 and T5 reached the asymptotic phase whereas T0, T2, T3 and T4 did not clearly reach it; thus, a number of OTUs remained undetected for these groups. Proteobacteria was dominated by the classes α-(≈11%), β-(≈8%), γ-(≈6%, p ≤ 0.001) and δ-(≈4%, p ≤ 0.05). Significant differences between T0 and all other groups were recorded for the γ-Proteobacteria class. Significant differences were also detected between T 0 and the treatments T2, T3 and T4 for the δ-Proteobacteria class (Figure 4).
The variation on microbiota composition through PCA analysis of A. ventricosus treated with different homeopathy drugs may be explained by two principal components (PC1 = 51.64%; PC2 = 14.31%) with a cumulative variance of 65.95%. Treatments T2, T3, T4 and control C3 were found in the negative PC1 associated to the phyla Spirochaetes, Fibrobacteres, Tenericutes and Chlamydiae, and other less abundant phyla whereas T0, T1, T5 and the control C2 were found in the positive PC1 associated with Proteobacteria, Actinobacteria, Bacteroidetes, Firmicutes, Cyanobacteria and other less abundant phyla. These results agreed with the relative abundance and composition at Phyla level described above (Figure 6).

DISCUSSION
Diseases related to bacteria and viruses cause large mortalities and economic losses in global aquaculture industry and may even affect natural populations of molluscs, crustaceans, and fish. Bacteria are naturally present in the GIT microbiota of animals, including bivalve mollusc. Some species or strains are facultative pathogens that can attack susceptible animals. Pectinid molluscs, such as the Catarina scallop A. ventricosus feed by filtration and can become reservoir of potential pathogens transmissible to other organisms of the same or different species, including man (27).
In this study, the juvenile scallop treated with homeopathic nosodes T1 and T2 showed significant differences in shell-length-growth rate (μm d-1) with respect to control groups (p <0.05). The nosodes are widely used in homeopathy; they are broad-spectrum preparations prepared from biologic materials, such as cultures or clinical samples of microorganisms (e.g. bacteria, fungi and viruses), parasites, diseased tissues (cancerous tissues), or decomposition products from animals an human being (28).
Kiarazm et al (29)  The greatest survival rate corresponded to the treatments T3 and T5 with respect to the control group C3. HOM-treatment T3 consisted in decimal dynamisations from sodium metasilicate and phosphoric acid TMs. In rats, sodium metasilicate has a beneficial action, because silica influences the complex network of cytokine interaction that regulates the immune response (31).
Nutrition is of most importance environmental factor which interact with other biological processes of life, and has a profound influence on several immune defense mechanisms. Deficient or excess levels of nutrients alter the immune system, and this reflects in measurable effects on several immune functions during in vivo or in vitro assays, but also it can modify infectious disease processes (32).
Homeopathic Phosphoric acid has the ability to treat gastrointestinal conditions related to malnutrition and poor assimilation of food; is used to promote the force of exhaustion and general vigor, and it has been related to the improvement of overall performance in A. ventricosus (14). The sinergical interaction between sodium metasilicate and phosphoric acid in decimal (1D) homeopathic dilutions demonstrated to increase the immune response in juveniles of A. ventricosus; when this treatment was applied in a centesimal dynamisation (7C), it did not show a clear effect on immune response in Catarina scallop, but it increased the energetic reserves in mantle, muscle and digestive gland of the species, which are important energy storage tissues in mollusc bivalves (29).
T5 treatment is a commercial homeopathic medicine for human use (Vidatox ® ), whose active principle is the venom of the red scorpion Rhopalurus junceus, an endemic species of Cuba. Preclinical studies attribute antitumor properties (33) and according to Mazón-Suástegui et al (15), increased the survival and immunity of juvenile white shrimp Litopenaeus vannamei challenged with a pathogenic strain of Vibrio parahaemolyticus.
The microbiota associated with the GIT is a key factor to the development, immunity, and nutrition of the host (21). Most of those bacteria inhabits the surrounding environment (22) and are relevant to reduce host vulnerability against pathogens and their related infectious disease (22)(23)(24). In this study, the values of the Shannon index showed similar values in all the experimental groups from 2.75 to 2.83, which indicated that the microbiota associated with these organisms was diverse in all the groups.
A study by Trabal-Fernández et al (11) showed that microbiota associated to oyster spat (juvenile stage) was more diverse than those associated with adults of the species. In these sense, maturity of the GIT can influence the microbiota composition in marine organisms (34).
Proteobacteria was the dominant phylum in the GIT of A. ventricosus, both in the HOM-treated groups and in Un-treated controls, followed by Actinobacteria and Firmicutes. This Phylum was found predominant in early-and advanced juvenile and adult stages of the oysters C. corteziensis, C. gigas, and C. sikamea (11); in larval stages of C. gigas (35), C. virginica (36), Haliotis discus (37), Patella pellucida (38) and the scallop Patinopecten yessoensisis (39); as well as in fish and crustaceans such as Oncorhynchus mykiss and L. vannamei (40).
Most classes within Proteobacteria plays important roles in bivalve molluscs, being able to degrade cellulose and agar which are major components of their food, and to fix nitrogen in bivalve GIT (41). Our results agree with Lasa et al (42) who reported Proteobacteria as the most abundant phylum in oyster. In this study, Actinobacteria was the second phylum in relative abundance. Some Actinobacteria produces secondary metabolites which can protect the host against infections (43) and play an important role in mineralization of organic matter, immobilization of mineral nutrients, nitrogen fixation and protection of the environment, in addition to their cellulolytic- (44) and chitinolytic activity (45).
Firmicutes was the third most abundant phylum (≈12%) in the GIT microbiota of A. ventricosus in all the experimental groups, with a lower relative abundance than the two previous phyla. Firmicutes is a highly relevant group in aquatic environments and found in the oyster microbiota (11), being involved in complex enzymatic processes, such as the degradation and fermentation of polysaccharides (46). It is dominant in the GIT of herbivorous organisms and some Lactobacillus species can stimulate the immune system, protecting the host from invasion and establishment of pathogens (47).
The dominant classes in this study were Alphaand Betaproteobacteria; Alphaproteobacteria is dominant in marine environments (48) ( Figure  4). Our results agreed with Trabal-Fernández et al (11), who reported Proteobacteria as the most abundant in oysters; however, Alpha-and Gammaproteobacteria were more abundant in juvenile oysters, while the most abundant class in adults was Gammaproteobacteria, followed by Beta-and Alphaproteobacteria.
The predominant genera in the GIT of juvenile A. ventricosus were Symbiobacterium and Microbacterium, followed by Methylobacillus, Bacillus, Paenibacillus, Burkholderia, Nostoc, Methylobacterium and Leucobacter, common in the marine environment (42). Symbiobacterium was dominant (p≤0.05) in the HOM T5 treatment (Vidatox ® ), which acts stimulating the immune system. This a very widespread genus in the natural environment, but its limited knowledge hinders its isolation (49). Oyster shells are common sources of Symbiobacterium sp. and other bacteria related to its growth (50). Ishii et al (51) inferred that this group of bacteria could play an important role in the nitrogen cycle in anaerobic environments, so they are beneficial because participates in nitrogen fixation, favoring the nutrition of marine organisms.
The species Microbacterium maritypicum was more abundant (p≤0.05) in T1 and T3 with respect to the controls. This bacterium produces siderophores and may inhibit the growth of opportunistic pathogens lacking the ability to produce iron (52). Iron is a limiting bioactive metal in seawater but essential for growth of marine bacteria (53), so these HOMtreatments may contribute to reduce pathogens.
HOM treatments T2 and T5 increased the abundance of Symbiobacterium toebii with respect to T0 (p≤0.05). Symbiobacterium species are symbiotic (54) but show a mono-marked growth if CO 2 or bicarbonate is available (55). Joong-Jae et al (54) showed that Geobacillus kaustophilus, Escherichia coli and Bacillus subtilis had effects on the growth of S. toebii, indicating that there are growth promoting factors, widely present in various bacterial strains.
In conclusions the composition and diversity of the gastrointestinal microbiota of the marine mollusk A. ventricosus has been described for the first time, with dominance of the phylum Proteobacteria, Actinobacteria and Bacteroidetes. Homeopathic medication promoted growth and survival, but also the proliferation of Microbacterium maritypicum (T1, T3) and Symbiobacterium toebii. (T2, T5). This allows us to assume that aquaculture homeopathy has the potential to modulate that microbiota and improve seed production in the laboratory.

Conflict of interests
The authors declare no conflict of interests.