Analysis of the Soil Microbiome of a Los Angeles Urban Farm

&is research aimed to address health and food security challenges affecting Los Angeles and beyond in our own backyard. One possible way to combat these challenges is through biotechnology.&e purpose of this research was to characterize the arable soils at Los Angeles Pierce College and their potential uses. &e methods employed were collection of soil physical and chemical properties, cultivation of soil microorganisms, and 16S rDNA sequencing with a citizen science platform. Statistical analysis in R revealed associations between field conditions and life in the soil. &e testing of plant-based antimicrobial agents on antibioticresistant bacteria revealed that plant oils could be more efficacious for control than common antibiotics. We found evidence of plant growth-promoting, antibiotic-producing, and bioremediating bacteria in different fields with distinct microclimates based on the 16S rDNA taxonomy results and inferred functional profiles. &e major findings included an abundance of Streptomyces, Bacillus, Methylobacteria, and others with desirable functions. Numerous Streptomyces isolates cultivated on selective media demonstrated the ability to reduce nitrate. &ese results will empower future trials and product development by identifying target soil microorganisms and the most viable fields.


Introduction
On a broad scale, it is urgent to address threats to public health and food security. ese issues include increased food demand with a growing population, a limited number of arable acres, and threats to human health such as antibiotic-resistant bacteria and chemically contaminated land and aquifers [1][2][3]. One possible way to combat these challenges is through biotechnology. e Pierce Farm in Los Angeles has the potential to promote human health by harnessing natural biological pathways responsible for detoxification arising from the soil food web. In addition, the Farm can cycle nutrients from the landscape and provide a local source of produce to decrease the city's foodprint. e urban foodprint is the burden created by unsustainable consumption and waste generation [4].
Microorganisms can be unified with social and economic growth in a sustainable manner [5,6]. For example, soil microorganisms can be biopharmed which perform the desired functions such as antibiotic production and bioremediation. e first step is to mine genomic data for desirable functions [7,8].
e purpose of this research is to characterize the arable soils at Los Angeles Pierce College and their potential uses. We are interested in the interaction between humans, plants, livestock, and soil microorganisms in urban fields which may help us respond to health challenges. Recently, soil microbiomes have received attention by researchers who searched for novel taxa and evaluated effects of environmental factors such as urbanization, fire, and different plant hosts [9][10][11][12][13].
us, soil physical and chemical properties data and 13 DNA samples have been collected during 2018 and 2019. e data include bacterial colony forming unit (CFU) counts, pH, EC (electrical conductivity), TDS (total dissolved solids), texture, N, P, and K concentrations, and organic matter.
We hypothesized that the amount of bacterial growth in the soil can be predicted by soil physical and chemical properties. What are the most important factors for determining the number of colony forming units that may be cultivated from a field? A transformed linear statistical model was evaluated for hypothesis testing.
Field 19 at Pierce Farm is bordered by a newly constructed apartment development and can help us study the effect of urbanization on microbiomes. We hypothesized that there would be a higher proportion of antibiotic resistance genes in Field 19 due to population density [14].
Inferred functional profiles provided evidence to test this hypothesis.
We hypothesized that plant growth-promoting bacteria will be present in the Pierce College Arboretum with the ability to reduce nitrogen and complete the nitrogen cycle. To evaluate these claims, inferred functional profiles generated in Piphillin were considered.
Streptomyces were isolated and tested for nitrate reduction activity. e effects of antibiotic-resistant bacteria are an important issue for human and livestock health and a high priority for the World Health Organization [15]. Lavendula stoeches, Ocimum basilicum, and Eucalyptus pulverulenta oils were tested using a standard microbiological assay to investigate effects on antibiotic-resistant bacteria isolated from several fields on the Pierce Farm. Our hypothesis was that the bacterial isolates would be better controlled by plant extracts than antibiotics. We hypothesized that Pierce Farm can also provide antibiotic-producing bacteria and possibly new antibiotics for biopharming. We looked to sequencing taxonomy results and inferred functional data to find out which fields may be most fruitful in this area.
We searched for endophytic bioremediating bacterial sequences in Field 21 within 16S rDNA sequencing results. We hypothesized that these bacteria can clean up soils contaminated by TCE (trichloroethylene) and hydrocarbons. Functional analysis data and taxonomy sequencing results were used to support or disclaim this hypothesis.

Description of Areas of Interest.
Each region of the Pierce Farm tested had a distinct microclimate. A brief description of each field is given in Table 1.

Antibiotic-Resistant Bacteria.
We tested Lavendula stoeches, Ocimum basilicum, and Eucalyptus pulverulenta oils to investigate effects on antibiotic-resistant bacteria. e plants are known to contain terpenes with antibacterial properties such as eucalyptol and beta-carophyllene [16]. e first step was isolation of antibiotic-resistant bacteria from the rhizosphere of the LAPC Farm. e sources of the antibiotic-resistant bacteria are summarized in Table 3. Penicillin (10 units)-resistant isolates and tetracycline (30 mcg)-resistant isolates were used for the plant extract assays.
Plant essential oils were procured by Soxhlet extraction which ran for three cycles. Standard microbiological assays outlined by Redfern et al. [17] were performed to determine the efficacy of the oils in controlling bacterial growth. Bacterial isolates were cultured overnight in nutrient broth and pipetted onto nutrient agar plates. 20 microliters of the crude ethanolic extracts from dry plant biomass were pipetted onto paper discs atop the nutrient agar plates. A summary of the essential oil treatments is given in Table 2.
e diameter of the area of inhibition (AOI) of bacterial growth was measured. e standard errors were graphed about the mean AOIs in Excel, allowing visual inspection of the dataset for a significant difference between oil treatments [18]. In spring 2019, basil oil was tested alongside eucalyptus oil. is was to determine if there were any antibiotic-resistant isolates that reacted better to treatment with basil or lavender oil, in spite of the trend toward eucalyptus oil having the best performance.

Metagenomic Analysis.
For the 16S rDNA sequences, ranacapa was used for diversity analysis [19]. Piphillin was used for inference of bacterial functions [20]. 16S ribosomal DNA sequences are the typical sequences used to identify bacteria.
Ranacapa diversity analysis, functional analysis with Piphillin, and STAMP were carried out according to the instructions in the Pipeline for Undergraduate Microbiome Analysis (PUMA) Manual [19].
STAMP was used for analysis and visualization of metagenomic profiles [21]. Soil classification data were collected using the sampling and testing methods published in the Introduction to Soil Science lab manual by St. Clair et al. [22]. e metagenomic soil DNA samples were collected according to the CalEDNA citizen science instructional videos and sent to UCLA for sequencing using Illumina technology. e total DNA was extracted using the Pow-erSoil DNA Extraction kit, followed by amplification using 16S metabarcoding primers and purification [23]. Quality control was carried out in QIIME.
e sampling locations at the Pierce Farm and the Arboretum, shown in Figure 1, were spread out in order to capture the most diverse picture of the soil microbiome.

Statistical Analysis of Soil Physical, Chemical, and Biological Properties Data.
e results of an exploratory data analysis in the dlookr package for R 3.6.0 revealed that pH and moisture percentage were the only normally distributed values [24]. Most of the variables were skewed and exhibited a better distribution with a log transformation [25,26]. Also, since the environmental data collected here had various units of measure, a transformation was needed [27,28].
Results for multivariate regression of transformed data for CFU on plate D as a function of field, pH, organic matter content, N, P, K, and moisture were significant (R 2 � 0.89).
is supports the initial hypothesis. However, none of the individual variables had a significant p value. It is possible that using the presence of antibiotic-producing bacteria as a factor would help to improve the goodness of fit, if that data had been available.
Simple linear regression models showed that Field was significant for CFU D, organic matter, moisture, and P. Contrasts revealed no significant differences between fields for the moisture parameter. e two that were the most different were Marquis and Field 12 (p � 0.19). e results of Tukey's HSD on moisture did not show honest significant differences between groups. is determination is supported by Figure 3. It should be noted that 2018-2019 was an atypically heavy rainfall year.
e Arboretum showed a trend toward higher levels of organic matter by percentage, although there was a wide range of values for the areas sampled. e range was similar to all fields except Field 12 which is cultivated and had significantly lower organic matter by percentage. Phosphorus levels were highest in Fields 15, 28, Arboretum, and Marquis, although not significantly higher. All fields had similar levels of phosphorus. e highest number of bacterial colonies grew on D plates (10 − 4 dilution) in Marquis field, which is significantly higher than those cultivated Field 12, and similar to other fields on the farm (Figure 4).

Principal Components Analysis (PCA) of Soil Classification Data.
e PCA was run using the prcomp() function in R with imputed data from dlookr package and plotted using ggplot2.
e principal component analysis showed some separation in the data by Term between Spring 2019 and Fall 2019. is is because in Spring 2019, samples tended to have lower values for TDS and P as a group, possibly owing to leaching related to the amount of rainfall received during fall and winter [29,30].
As shown in Figure 5, however, there was no clear separation or clustering of samples by texture, term, or field which emphasizes similarity between samples overall. e variables that contributed the most to the PCA were the imputed values for the TDS and K concentrations in parts per million. Samples that separated out high on PC1 and positive on PC2 tended to have imputed missing values for K and lower values for the TDS experiment. e highest sample on PC2 was a Spring 2019 sandy sample from Marquis C which had a low TDS (32 ppm) and low K-measured value (25 ppm). is makes sense due to the low cation exchange capacity of sand. In the imputed data set, pH was negatively correlated with N, and EC was positively correlated with P. Figure 6 summarizes the strength of the correlation between the soil classification factors. e correlation matrix revealed moderate negative associations (0.6-0.8) between the CFU counts for plates C (10 − 3 ) and D (10 − 4 ) and percentage of silt by     Figure 2. Table 4 shows selected metadata for the 16S isolates.  hydrometer method; this was mirrored by moderate positive associations between the CFU counts on plates C and D with percentage of sand by the hydrometer method. is suggests the bacteria that were cultivated on these plates were mostly aerobes. In addition, there was a moderate positive association between the CFU counts on plate C with pH. ere were moderate positive associations between P and TDS and slight negative associations between pH and both N and K. EC and TDS showed a strong correlation as expected; there is a similar association between clay and sand percentages by hydrometer by nature of the relationship between the factors.

Antibiotic-Resistant Bacteria.
During Fall 2018, essential oil trials on antibiotic-resistant bacteria, both eucalyptus and lavender oils, were effective in controlling the bacterial growth better than the antibiotics tested. e antibiotics tetracycline and penicillin sometimes had no effect, and the isolates displayed an AOI of 0. e mean of AOI for the eucalyptus treatment was 17.35 mm; the mean of AOI for the lavender treatment was 13.37 mm. A t-test confirmed that the bacteria were significantly more susceptible to the eucalyptus extract versus the lavender product (p � 0.05), as shown in Figure 7.
A visual inspection of the graphs of the standard errors about the mean for Spring 2019 trials revealed that basil and lavender were better at controlling some isolates, as shown in Figure 8.
Basil was the best at controlling antibiotic-resistant isolate 3 from Marquis C, and lavender was best at controlling isolate 7 from Field 28; the other isolates were better controlled by eucalyptus. Overall, this supported the hypothesis that plant extracts could control resistant bacteria better than common antibiotics.

Diversity
Analysis. Next, the 16S rDNA sequencing results will be considered, beginning with diversity. e alpha diversity observed describes the number of taxa present, but does not describe the abundance [19]. It is interesting that sample S3 shows the lowest number of bacterial taxa in Figure 9, although the highest amount of plant diversity resides in the Arboretum.

Statistical Analysis of Inferred Metagenomic Profiles.
When a principal coordinate analysis is carried out in STAMP and visualized with respect to the type of plant rhizosphere where soil was collected and metagenomic profiles, the samples that separate out the most are dead brome in Field 19 and cactus from the Arboretum. is separation is visible in Figure 10. e sycamore-associated sample from the Arboretum separates out PC1 vs. PC2 and on PC2 vs. PC3. is suggests that most of the variation in the functional hierarchy of the samples is coming from Field 19 and two different samples from the Arboretum. e results from the STAMP analysis are statistically significant for comparisons between Field 21 and the Arboretum only due to the need for additional DNA samples from other fields. However, information from all of the samples can help to guide researchers in the future about where to sample to find bacteria with the indicated functions. In the L2 pathway, Field 19 did not have a trend toward a higher proportion of antimicrobial drug resistance. e Arboretum had a higher proportion of drug resistance antimicrobial genes in the L2 pathway, as shown in Figure 11. is did not support the hypothesis for Field 19 about drug resistance. e sample from Field 19 in fact had a trend toward a lower proportion of human disease in L1 and bacterial disease pathway genes in L2 than the other fields studied. However, Field 19 had a trend toward the highest proportion   Applied and Environmental Soil Science     Observed diversity S10 S12 S13 S11 S14 S2 S3 S5 S6 S7 S8 S9 S1 Sample ID  Applied and Environmental Soil Science of vancomycin-resistant genes in the L3 pathway which partially supports the hypothesis. Several fields exhibited functional profiles that suggested the bacteria could be isolated and cultured to make a product. Neomycin, kanamycin, and gentamycin biosynthesis (L3) had a trend toward being the highest in terms of the proportion of sequences related to that function. In the L3 pathway, Field 19 had a high proportion of hits in potential bioremediation genes for chlorocyclohexane, chloroalkene, styrene, atrazine, and toluene degradation. Field 21 had lower numbers of reads sequenced for this function.
is partially supports our hypothesis and lets us know that in the future we should try to culture bacteria from Field 19 which is downgrade from Field 21.
Field 19 also had a high proportion of genes related to remediation of plastic byproducts, including BPA, styrene, and caprolactum degradation. ere was a trend toward a high proportion of sequences attributed to penicillin production occurring in Field 19. In addition, the highest proportion of streptomycin and tetracycline production genes was attributed by Field 19.    Applied and Environmental Soil Science 9 Field 19 evidenced a high proportion of the functions related to terpene degradation, including geraniol, limonene, and pinene degradation. Taking samples from there in the future could help researchers explain why some bacteria are more tolerant of certain plant extracts than other extracts. Understanding the underlying genetic basis of antibiotic resistance can aid the development of treatments. e Arboretum had a higher proportion of mineral absorption genes in the L3 pathway and a higher proportion of zeatin synthesis genes in the L3 pathway. Zeatin is a cytokine which promotes plant growth, bushy growth habit, and vegetative shooting [30]. Furthermore, the Arboretum had a higher proportion of hits in the plant-pathogen interaction pathway in L3.  Sum of S10 Sum of S9 Sum of S8 Sum of S7 Sum of S11 Sum of S12 Sum of S13 Sum of S2  e Arboretum also had a higher proportion of antibiotic resistance genes in L2 and the highest proportion of Staph infection pathway genes in L3. Furthermore, the Arboretum exhibited a trend toward having the highest proportion of sequences related to E. coli biofilm formation.
is could be a source of genes for radiation cleanup. e Arboretum similarly had a high proportion of hits related to Pseudomonas aeruginosa biofilm formation which could be useful for plant growth-promoting compost teas.
Field 15 showed a trend toward a high proportion of sequences involved in polycyclic aromatic hydrocarbon (PAH) degradation. PAHs in the urban environment are associated with carcinogenicity and have been shown to be remediated by microbes such as Methylobacteria sp. [31][32][33]. Field 15 exhibited a high proportion of sequences related to nitrogen metabolism as well. is makes sense due to the horse manure composting in the area.
Field 21 had a trend toward the highest proportion of hits in the L3 thermogenesis pathway, as well as the highest proportion of hits in the L3 methane metabolism pathway. e proportion of sequences for methane metabolism was significantly higher in the Arboretum than it was in Field 21  ( Figure 12). Methanogenic bacteria can be associated with sulfate reduction as part of their metabolism [36]. Field 21 also had the highest proportion of plant hormone signal transduction pathway functions in L3. In L2, Field 21 had a trend toward a higher proportion of sequences related to environmental adaptation. In L3, Field 21, and Field 19 had similarly high proportions of dioxin degradation pathway hits.
e highest proportion of siderophore production functions occurred at low moisture. is data point corresponded to Field 19.
Marquis C had the highest proportion of Salmonella infection genes in L3. is makes sense due to the presence of chicken waste and bedding in the compost. Marquis C and D ranked high in the proportion of sequences related to xylene degradation which could also be a source of genetic material for bioremediation.

Relative Abundance of Methylobacteria Taxa.
Evidence of Methylobacteria populii growing in Field 21 was found at the location shown in Figure 13. is endophytic bacterium, which is usually symbiotic with poplar trees, was growing in the free soil near the plant rootzone of a California native mature poplar tree [35]. e sum of Methylobacteria sequences by species is shown in Figure 14. A possible application of these bacteria would be to provide inoculated California native black poplar trees for habitat restoration and bioremediation near lands contaminated by hydrocarbons [36]. For example, nearby Pierce College is the former nuclear engine testing site, the Santa Susana Field Laboratory.
e Field Lab contains soils which are contaminated by TCE [37]. TCE is a potent carcinogen and a risk to developing babies. e most effective method for TCE removal from water is aeration combined with activated carbon treatment [38]. A challenge with TCE degradation is that it resists oxidation (Ibid). Recently, Boeing has had the Santa Susana Field Lab Reclassified as an Open Space Preserve [39].

Relative Abundance of Plant Growth Promoters and
Antibiotic Producers. We have found evidence of Bacillus and Streptomyces species which are known to produce antibiotics growing in the Pierce College Arboretum, Field 15, and the Nature Canyon (Field 21). It is possible that there are other bacteria or subspecies present that are capable of producing novel antibiotics. For example, recently, researchers in the EU and the US have discovered new antibiotics from soil. is is partially owing to the iChip, which is a new device for growing bacteria which were previously recalcitrant to common bacterial cultivation methods [40]. e distribution of reads for Streptomyces spp. spans all of the fields at Pierce College. is determination is supported by Figure 15. Each field has also provided isolated strains of putative Streptomyces (Table 5).
ese results support our hypothesis wherein we predicted the presence of plant growth-promoting bacteria on the farm. e distribution of Bacillus-related genuses for all fields shows Bacillus was the most common of those queried, as illustrated by Figure 16. Each field has been confirmed to provide sequence reads and isolates for plant growth promotion and antibiotics, which supports our hypothesis. In addition, the natural presence of Bacillus subtilis in these fields may be an organic competitive advantage for

Clostridium Pseudomonas
Sum of S1 Sum of S4 Sum of S3 Sum of S6 Sum of S5 Sum of S10 Sum of S9 Sum of S8 Sum of S7 Sum of S11 Sum of S12 Sum of S13 Sum of S2 Sum of S2 Sum of S13 Sum of S12 Sum of S11 Sum of S7 Sum of S8 Sum of S9 Sum of S10 Sum of S5 Sum of S6 Sum of S3 Sum of S4 Sum of S1 production of commonly fumigated crops such as strawberries and carrots [41,42], which are likely to grow well in the area. Our 16S sequencing taxonomy data indicates the presence of plant growth-promoting bacteria such as the biopesticide Bacillus subtilis [43], nitrogen-fixing Oxalobacteraceae spp., and auxin-producing Streptomyces spp. on the Pierce College Farm. e Oxalobacteraceae spp. also help break down organic matter and mineralizing plant nutrients, and many Streptomyces spp. are capable of chelating iron and reducing nitrogen. Both are capable of solubilizing phosphate .
ere was a high relative abundance of the archaea Candidatus nitrocosmicus which oxidizes ammonia and is part of the nitrogen cycle [44]. Pseudomonas spp. are known to enhance yield in red pepper and were also present in Pierce fields [45]. Clostridium spp. are sulfite-reducers, and compost tea containing Clostridia has been shown to increase greenhouse yield of peppers [46].
We isolated strains of putative Streptomyces spp. grown on selective media and tested them for nitrate reductase activity. e outcomes of the nitrate reductase test are summarized in Table 5. Most of the samples exhibited nitrogen reduction capabilities except for two isolates from Marquis D, putative Streptomyces isolates 14St and 15St. e samples which had the highest number of reads from zeatin-producing bacterial genuses were S6 from Field 15 and S1 from Field 19. Agrobacterium sequences were not highly abundant in the samples. Rhizobium had the most hits from a variety of fields and was represented everywhere on the farm, including the Arboretum. e samples with the highest number of reads for Clostridium spp. were S10 from Marquis D. Figure 17 shows the relative abundance of these bacteria.
We should grow plants for beauty, food, and varietal trials using sustainable local inputs to offset the high cost associated with the high productivity of urban farms [47]. Currently, avocado and citrus orchards are proposed. Avocados, strawberries, clover, maintaining secluded preserves a sustainable urban food garden to feed hungry students, and maintaining the urban tree park are popular ideas. Since several of the fields have not been cultivated in over five years, obtaining a USDA Organic certification is a viable option which would add value to the land. e social benefit of bioremediation is the amount of revenue saved by not having to pay costs of remediation or medical bills for humans who experience ill effects from contamination.
For example, recently NASA announced that it is not feasible to clean up its portion of the Santa Susana Field Laboratory located near Pierce College because the cost exceeds $500 million [48].

Conclusion
ese results will empower future trials and product development by identifying target soil microorganisms and the most viable fields. Cultivation of these soil microorganisms along with crops can give our community a healthier and more sustainable future.
Terpene extracts from eucalyptus, basil, and lavender may have therapeutic value since our hypothesis was supported ( Table 6). Use of plant material or extracts supplemented in animal feed may lead to decreased need for antibiotics. An extension of this research would be to test growing the plants which produce the phytoactive chemicals near livestock grazing areas; the crops may help to control antibiotic-resistant bacteria or at least keep them from spreading.
Future work to quantify the phytoactive compounds in the plant extracts is necessary. e possibility of synergistic effects of antibiotics and medicinal herbs as discussed by Budimir [49] should be investigated.
Different field conditions recorded such as field location, pH, organic matter content, N, P, K, and moisture were shown to be related to the amount of bacterial growth in the soil. It is clear that fertility is a collection of these properties within the top 6 inches of the plant rhizosphere. Further research should record the presence of antibiotic-producing Low amount of human disease and antimicrobial functions found in Field 19 bacteria on agar plates as a factor in number of CFUs present in order to improve the statistical model. Soil is a renewable resource over a long period of time. e land at Pierce College is classified as prime farmland (USDA NRCS). In the meantime, the value of the land for other enterprises increases over the years as buildable land becomes more scarce and this is one of a few remaining open spaces. Preserving this unique, biodiverse resource in the urban landscape is crucial for education and for productivity in its own right.
Preserving these plots of land for agriculture and education could yield tangible economic benefits from products such as food, antibiotics, and biofertilizers. Education also gives social benefits, as well as ornamental gardens. Imagine a world without beauty! e bacteria which can be isolated from the Pierce College urban farm can potentially aid in pollution remediation. Trees, turf, and other plants perform carbon fixation and have a cooling effect on the neighborhood. Further research should consider other soil microflora and macrofauna, such as the nematological composition of the fields.

Data Availability
Data from the 16S survey are available at https://www.ncbi. nlm.nih.gov/sra/PRJNA556308. Disclosure e content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. Daila Melendez is currently a horticulture BS student at Oregon State University, Corvallis, OR, USA 97331. Serina Reitz is currently a landscape architecture BS student at California State Polytechnic University, Pomona, CA, USA 91768. Destini Kananipour is currently a biology BA student at California State University, Northridge, CA, USA 91330.

Conflicts of Interest
e authors declare that there are no conflicts of interest.