Patent Grant
11385215
This disclosure generally relates to soil health indicators, gene functions, and agronomy.
The soil microbiome includes thousands of organisms, including bacteria, fungi, nematodes, and insects, among other microbes. These organisms collectively contain tens to hundreds of thousands of unique genes, which encode proteins that form enzymes, which in turn catalyze reactions of microbial metabolic pathways. Metagenomics (also referred to as environmental genomics or community genomics) is the science of describing the profile of microbiomes detected in a biological sample such as soil, here including the functional genes underlying microbial metabolism. As one application, it is desirable to predict whether a farmer's field will produce a high or low crop yield, and whether the crops will develop disease.
Further, it is challenging to determine the impact of microbiome in soil on crop yield and disease pressure. For example, crop yield may be altered by microbial control over nutrient cycling such that added nutrients are lost due to microbial activity. This yield loss may not be a direct function of the abundance of the underlying microbes, but could also arise due to their interactions with other factors in soils. Moreover, in many cases these gene abundances may not directly correspond to the abundances of known, named organisms.
The disclosed embodiments have advantages and features which will be more readily apparent from the detailed description, the appended claims, and the accompanying figures (or drawings). A brief introduction of the figures is below.
An analytics system uses soil health indicators to determine metrics for soil samples, for example, indicating performance of crops grown in geographical locations having the soil samples.
In various embodiments, a method includes determining nucleic acid sequence reads of genetic material in a soil sample from a geographical location. The method further includes determining, by processing the nucleic acid sequence reads, a first set of measures of a plurality of gene functions represented in the nucleic acid sequence reads. The method further includes determining, by processing the plurality of gene functions, a second set of measures of a plurality of metabolic pathways of microorganisms present in soil at the geographical location. The method further includes determining, by processing the plurality of metabolic pathways, a third set of measures of a plurality of soil health indicators of the soil at the geographical location, the plurality of soil health indicators including a plurality of levels of granularity. The method further includes determining a measure of an agronomic attribute of the soil at the geographical location as a function of the first, second, and third sets of measures. In some embodiments, the measure of an agronomic attribute of the soil is transmitted to a client device for display.
In an embodiment, determining the second set of measures of the plurality of metabolic pathways of microorganisms comprises determining a function of measures of constituent genes from the plurality of gene functions. In an embodiment, determining the third set of measures of the plurality of soil health indicators comprises determining a function of the plurality of metabolic pathways corresponding to the nucleic acid sequence reads.
In an embodiment, the plurality of gene functions includes one or more of: nitrogen cycling, phosphorous cycling, carbon cycling, or oxygen availability. In an embodiment, the agronomic attribute is nitrogen mineralization and the plurality of metabolic pathways includes at least secreted proteases and urea mineralization.
In an embodiment, the method further includes determining a recommendation for treatment of the soil at the geographical location using at least the measure of the agronomic attribute. In an embodiment, the method further includes filtering out at least one metabolic pathway of the plurality of metabolic pathways.
In an embodiment, the plurality of levels of granularity includes a level closer to biological attributes and another level closer to agronomic attributes.
In various embodiments, a method includes receiving metadata describing a soil sample, where the metadata indicates one or more types of crops grown in a geographical location having the soil sample. The method further includes determining nucleic acid sequence reads of the soil sample. The method further includes determining, for each nucleic acid sequence read of at least a subset of the nucleic acid sequence reads, functional descriptors of the genes represented in nucleic acid sequence reads. The method further includes determining reference metrics of soil samples from geographical locations in which the one or more types of crops were grown. The method further includes determining a metric of the soil sample using the functional descriptors and the reference metrics. The method further includes transmitting the metric to a client device for display on a user interface.
In an embodiment, determining the metric of the soil sample comprises determining a value of a soil health indicator of the soil sample using the functional descriptors. The method further includes determining a distribution of values of the soil health indicator for the soil samples using the reference metrics. The method further includes determining a percentile of the value with respect to the distribution of values.
In an embodiment, determining the functional descriptors comprises determining a plurality of microbial functional genes (also referred to herein as functional genes or microbial genes) in the soil sample. The method further includes determining, for each of the plurality of genes, a count of the genes in the soil sample. The method further includes normalizing these counts using a total count of genes in the soil sample or by other methods such as single copy marker genes, sequence read depth, etc.
In an embodiment, determining the nucleic acid sequence reads of the soil sample comprises extracting microbial material from the soil sample. The method further includes generating nucleic acid sequence reads of the microbial material. The method further includes filtering the nucleic acid sequence reads by read quality scores.
In various embodiments, a method includes obtaining a soil sample from a geographical location. The method further includes receiving metadata indicating the geographical location. The method further includes determining a plurality of functional genes in the soil sample. The method further includes determining, for each of the plurality of functional genes, a measure of the functional gene in the soil sample. The method further includes determining functional descriptors of genes represented in the soil sample using the measures of the functional genes. The method further includes determining reference metrics of soil samples from geographical locations within a threshold distance of the geographical location. The method further includes determining a metric of the soil sample using the functional descriptors and the reference metrics. The method further includes transmitting the metric to a client device for display on a user interface.
In various embodiments, one or more processors may execute instructions stored by a non-transitory computer-readable storage medium to control a computer system to perform steps of any of the above methods. In various embodiments, a system includes a sampling tube for obtaining a soil sample from a geographical location. The system further includes one or more processors and a memory, the memory storing computer program instructions that when executed by the one or more processors cause the one or more processors to perform steps of any of the above methods.
The analytics system 100 determines metrics of soil samples using soil health indicators. A soil health indicator is defined as a value of microbial driven function pertinent to agricultural production. A soil health indicator may reflect soil mineral and organic element availability, plant growth promoting factors, interaction with plant pathogens, crop performance, or other indicators of soil function or health. A soil health indicator may be derived by processing nucleic acids of a soil sample, for example, by sequencing deoxyribonucleic acid (DNA) or ribonucleic acid (RNA) to determine functional descriptors associated with microbes (also referred to herein as microorganisms or organisms) present in the soil sample. Soil health indicators may be used to predict physical attributes of crops (e.g., stem size, plant height, or fruit size), crop yield, or cycling of mineral and organic elements in soils including nutrients, carbon, and oxygen.
The analytics system 100 obtains or receives soil samples from users (e.g., of the analytics system 100) such as farmers or other third parties (e.g., agriculture-related companies). In some embodiments, the analytics system 100 or an entity associated with the analytics system 100 provides a sampling tube 280 to a user, e.g., as part of a kit for collection of soil sample or related information. The user collects a soil sample using the sampling tube 280 and returns the sampling tube 280 (e.g., via mail or other delivery methods) to the analytics system 100 for processing. An interior of the sampling tube 280 may be sterilized and may include a preservative solution, for example, to help maintain conditions of the soil sample or microbes present in the soil sample. The analytics system 100 may determine or indicate to a user a target volume, mass, or weight, of the soil sample to be collected using the sampling tube 280. The analytics system 100 may also provide sampling recommendations or protocols to users. For example, the sampling recommendations indicate a range of depth for soil collection (e.g., 0-6 inches below ground level), which may vary based on type of crop, geographical location, or other factors. This range of depth instructs the user how to conduct the sampling of the soil, including to what depth the sample should be extracted in a particular region in California for a specific type of crop (e.g., corn or a particular varietal of corn).
In some embodiments, the sampling tubes 280 are associated with a label (e.g., barcode or QR code) for tracking or identification. The analytics system 100 associates information describing users or soil samples with identification keys of obtained soil samples. For example, a soil sample is identified with an agricultural field, geographical location, or one or more types of crops. The information may include metadata, which is further described below with reference to
The analytics system 100 determines a metric of a soil sample in view of a “crop community,” that is, reference information associated with the soil sample. For example, the reference information includes data of other soil samples having similar conditions, in which same types of one or more crops were grown, treated with similar management or agricultural practices, or having other traits in common with the soil sample. In an example use case, a crop community corresponds to soil from a geographical region such as a city, state, or country. The crop community can be further differentiated by crops such as one or more vegetables, fruits, or nuts, among others. The analytics system 100 stores metrics in the soil health indicators database 102.
The analytics system 100 stores reference information in the reference database 104. The analytics system 100 may receive reference information from one or more data sources 120 or client devices 110. For instance, users of the analytics system 100 provide soil samples and information (e.g., metadata) describing the soil samples that the analytics system 100 may use as reference information for other soil samples. The analytics system 100 stores information derived using the soil samples or metadata as reference information in the reference database 104. Moreover, the analytics system 100 may associate reference information with associated metadata such as geographic location, one or more crop types, crop season, or other information describing an agricultural field. Accordingly, the analytics system 100 may perform lookup for reference information by querying the reference database 104 using metadata. For example, the analytics system 100 can lookup reference information for agricultural fields within a certain region or radius of geographical locations.
The analytics system 100 provides metrics to users, e.g., for presentation on a client device 110 of a user. The analytics system 100 may also derive recommendations from metrics regarding agricultural techniques. Based on metrics or recommendations, farmers or other users are informed as to a variety of actions that determine inputs or practices to use on fields, when to plant, where to plant, which crops to plant, or which varietals of those crops to plant, among other insights that may improve crop or soil health or performance.
A client device 110 comprises one or more computing devices capable of processing data as well as transmitting and receiving data over the network 130. For example, a client device 110 may be a desktop computer, a laptop computer, a mobile phone, a tablet computing device, an Internet of Things (IoT) device, or any other device having computing and data communication capabilities. The analytics system 100 may provide information to the client device 110 for presentation to a farmer or another user. The information may include metrics or recommendations determined by the analytics system 100 regarding soil samples or crops.
Though not shown in
For purposes of explanation, this disclosure uses soil samples generally as example use cases, though the embodiments described herein may be adapted for systems and methods using other types of biological samples or physical samples. For instance, the biological sample may be at least in part a liquid or aqueous sample used for growing plants in a hydroponics system. As a different example, the biological sample may be a sample of a gut microbiome of a subject (e.g., a human or another type of organism), and the analytics system 100 may determine metrics associated with physiology or other attributes of the subject.
In various embodiments, the analytics system 100 determines functional gene based soil health indicators (SHI) for one or more of Nitrogen cycling, Phosphorus cycling, Soil Carbon availability, and Oxygen availability. Individual indicators for each of these element cycles represent simplified versions of more complex biological transformations, derived from even more detailed lists of metabolic pathways and their constituent genes, present across a large number of microbial organisms with different gene sequence variations. The analytics system 100 determines soil health indicators using curation of multiple data sources for gene sequences, pathways, and element transformations, which may not necessarily be designed for application to soils or agriculture. The curation process may include one or more techniques known to one skilled in the art of genomics, metabolism, biogeochemistry, or agronomy. Techniques may include processes to determine where to obtain data, how to group the data, and how to analyze or distill the data, and how to experimentally validate the data, among others. Moreover, complex biological and chemical cycles may be further distilled into simplified indicators which are interpretable and relevant to agricultural operators using expert judgment, benchmarking, and validation.
In some embodiments, simplified, composite SHI are derived from genes which provide a representation of biological reactions transforming molecules containing each element in soils. Lists of genes involved in these transformations may be obtained and cross validated from multiple sources, many of which may be incomplete sources compared to the collected information. Scientific literature may be used in conjunction with public databases on metabolic pathways or schemes of gene categorization including the MetaCyc database, the Kyoto Encyclopedia of Genes and Genomes (KEGG) and SEED gene ontologies. Further data sources may be used to obtain additional gene annotation sources or models, including the UniProt, Pfam, and InterPro databases. These databases generally represent known molecular biology across organisms as organized for varying purposes which are not commonly organized to represent element cycling, soils, or agriculture.
As a general note, several of the most relevant genes are known to be horizontally transferred among microbes (such that organism names or name hierarchies are not necessarily deterministic of gene count). The number of these gene copies may vary within a given organism name or group (taxonomy). Further, microbial genes affecting plant growth and nutrition may have complex interactions with the soil environment, such that relationships vary by soil types or geographical locations. Methods to derive indicators of soil function (e.g., including gene groupings) may need to account for such interactions and localities during the validation of different indicator formulations or computed derivations.
The below descriptions provide various working examples of soil health indicators and their generation.
A. Processing Gene Lists from Databases
Existing genetic and metabolic databases may have incomplete representation of metabolism in soil, which makes it challenging to map information from these databases to soil functions. Some metabolic databases are organized around cellular metabolism (e.g., human or microbial), and may not map clearly to soil metabolism, including the omission of one or more metabolic functions.
The SEED protein database includes a list of genes that code for enzymes involved in phosphorus cycling. Yet this database is not organized around metabolic pathways and includes several genes without known environmental functions. Examples of such functions may include genes involved in intracellular phosphorus cycling, which may not have bearing on extracellular phosphorus cycling. These extraneous genes may present a false signal in agricultural soils.
The analytics system 100 filters out these extraneous genes from the data underlying soil health indicator calculations. Conversely, once relevant genes for phosphorus cycling are identified, it may become apparent that these genes are indeed present in metabolic databases (e.g. KEGG or Metacyc) and other protein databases (e.g., Pfam or Uniprot), even though the genes may not necessarily be explicitly organized or represented in metabolic maps.
B. Resolving Gene Ambiguity in Genomic Databases
In some embodiments, resulting element cycling gene lists are evaluated for consistency and inclusion among sources. The analytics system 100 may apply further curation to consider genes that might provide ambiguous information as the result of genes being involved in multiple pathways, with reversible reactions, having high sequence identity with similar genes with different functions (homologs), etc.
In some embodiments, from these gene lists, individual SHIs are calculated using hierarchical categorization schemes including both commonly recognized metabolic pathways, and expert curated classification (grouping) schemes, which may include excluding some genes with ambiguous functions in soils. Gene based indicators of soil health may be derived using expert judgment, benchmarking against agricultural data, and other validation processes to deliver simplified but accurate representations of underlying biology and chemistry which are relevant in an agricultural context.
Genomic or metabolic databases present genes as a collection of sequences with similar function or protein structure. These collections of genes may include sequences from organisms in which function differs from an agronomic attribute. For example, the KEGG database organizes genes by KEGG orthology (KO), which is a collection of genes with a similar function. Gene sequences from individual organisms associated with a given KO may be classified with a single identifier. In some embodiments, the analytics system 100 matches unknown DNA sequences from agricultural soils to a gene category or KO.
A given KO may include genes from organisms that perform different functions in the environment. For example, the KO identifier K00370 (and genomic neighbor K00371) are a collection of three similar genes involved in nitrogen cycling with different functions: narG, narZ, and nxrA (or for K00371, narH, narY, and nxrB). Notably, these genes are also represented with a single identifier in the Pfam database, which does not explicitly identify underlying differences in genes or the chemical reactions carried out by them. Reduction of nitrite (NO2—) to nitrate (NO3—) by nitrite oxidizing bacteria is carried out by the gene nxrA, while the opposite reaction (NO3— to NO2—) is carried out by narG. Yet the sequences from these genes are sufficiently similar to be grouped together in the KEGG database. Furthermore, narG is involved in two different metabolic pathways: dissimilatory nitrate reduction (NO3— to NH4) and denitrification (NO3— to N2O or N2). The analytics system 100 can use one or more different processes as described below to resolve the ambiguity of soil sequences.
In some embodiments, the analytics system 100 labels genes ambiguous, and several other genes may have similar ambiguity (e.g., pmo/amo genes, K10944, K10945, K10946, which may oxidize methane or ammonia). For example, as shown in Table 1 and Table 5 (below), there is a category for ambiguous nitrate reductases. The analytics system 100 processes genes such as these with ambiguous functions in one or more different ways (described below) with respect to derivation of gene counts and subsequent counts of metabolic pathways and soil health indicators.
In some embodiments, the analytics system 100 removes the counts of these ambiguous genes, which may be non-informative for metabolic pathway abundances. In other embodiments, the analytics system 100 examines, in soils, the correlation of these ambiguous genes with other non-ambiguous genes in a pathway, and assigns or splits the counts of the genes in a pathway accordingly.
In another embodiment, the analytics system 100 may process ambiguous gene categories by splitting the underlying sequences in the reference data to create new unambiguous categories. This may be accomplished, for example, through expert knowledge of the organisms performing these functions (and pruning of reference DNA sequences accordingly), or by construction of a phylogenetic tree of gene sequences, where coherent branches of related sequences are known to carry out given sub-functions.
C. Mixed Metabolic Signals Underlying Single Soil Functions
A different challenge with the formulation of soil health indicators occurs where a given function of soils relevant to agriculture is associated with multiple metabolic pathways, which are related by their net output but can present different signals in the environment. For purposes of explanation, the following example references data on functions representing nitrogen mineralization (i.e., conversion of organic nitrogen forms in soil organic matter to NH4). The analytics system 100, in this example, uses nitrogen mineralization as an agronomic attribute (e.g., function) to determine an amount of nitrogen obtained by a crop from soil organic matter, which may reduce the need for additional fertilizer. The analytics system 100 uses a combination of agronomic, biogeochemical (e.g., nutrient cycling), and metabolic knowledge to curate relevant genes. For example, nitrogen mineralization may involve the breakdown of organic nitrogen containing compounds like proteins and amino acids.
In an embodiment, a list of the most directly relevant genes for proteases (protein degradation) includes those that are released extracellularly into soil by bacteria or neutral and alkaline proteases (npr and apr). These genes release urea from proteins, such that ureases (ure) may also be an indicator of protease activity. Moreover, the analytics system 100 considers urea transporter genes as an additional index of this activity. However, urea is also a form of nitrogen that may be commonly applied to soil in agriculture. An additional form of organic nitrogen in soils is chitin (degraded by the chi gene), a non-polymer protein found in insects and fungi. A common enzyme assay for nitrogen mineralization in biogeochemical studies of soil is NAG (N-aceyl-β-D-glucosaminidase), which measures potential degradation of chitin.
Each of the above described gene functions for microbial utilization of soil protein, urea, and chitin fall within the biogeochemical and agronomic descriptions of soil nitrogen mineralization. The analytics system 100 considers such consistency or prior (gene and metabolic) inconsistencies to determine a soil health indicator. The effectiveness of a soil health indicator in reflecting these processes may vary based on the mixture of one or more underlying processes, their individual responses to different soil chemical factors, or potential signal driven by other factors (e.g., shaping gene abundances) and environmental noise (e.g., unknown variation). Resolution of soil health indicator effectiveness and determination of the ideal gene and metabolic factors underlying an indicator of soil health may require data obtained from real world soil samples with characterized variation or data matching expected environmental controls (e.g., nitrogen mineralization rate data or chemical proxies like soil carbon-to-nitrogen ratios, and amino sugars). Further, the analytics system 100 determines the weights of underlying genes and metabolic pathways in an indicator by inspection of relationships or machine learning in data comparing soil genetic potential to expected environmental relationships.
D. Calculation of Soil Health Indicators
In various embodiments, the analytics system 100 determines soil health indicators by assigning soil DNA sequence data to the corresponding genes in reference databases, in order to determine counts of each gene within each soil sample. For example, the analytics system 100 may use bioinformatics processes to assign DNA sequences match to nitrogen mineralization genes (nprS, aprA, cynS, URE, ureA, ureB, ureC, urtA, urtB, urtC, urtD, urtE, HEX, HEXA B, csn, nagZ, sacC, cynS, chiA) as shown in Table 5 (below). The analytics system 100 may normalize gene counts by total reads or gene hits, rarifaction, normalization by single copy marker genes, or other transformations. The analytics system 100 may combine reads or normalized read counts within a protein where distinct subunits of a protein are represented. For example, in the urease gene, counts are represented as the mean counts for the three subunits (ureA, ureB, ureC). As another example, the nif gene for nitrogen fixation proteins is also composed of several subunits, nifH, nifD, nifK, nifW. The nifH subunit codes for the metalloprotein containing the active site of nitrogen fixation. The analytics system 100 uses these codes to index biogeochemical activity of nif. A different perspective on these genes is held from a biochemical perspective, where several organisms are known to harbor nifH genes, but are incapable of nitrogen fixation due to missing helper proteins from nifD, such that nifD might be considered a more suitable indicator of nitrogen fixation activity.
Once gene counts are determined, the analytics system 100 combines the counts into a function of the counts (e.g., an aggregate count) for a given metabolic pathway (e.g., nitrogen pathways in Table 1 and Table 5) of which these genes are members. For nitrogen mineralization, these pathways may include urea mineralization (URE, ureA, ureB, ureC), urea transporters (urtA, urtB, urtC, urtD, urtE), chitinases (HEX, HEXA B, csn, nagZ, sacC, cynS, chiA), and secreted proteases (nprS, aprA, cynS), as shown in Table 5. The analytics system 100 may weight, add, or eliminate from indicators these genes within each of the pathways based on relationships in agronomic data, for example, relative fit with rates of nitrogen mineralization.
The analytics system 100 represents or combines counts of metabolic pathways in one or more ways in an agronomic indicator of soil health, where nitrogen mineralization can be viewed as a source of nitrogen from soils (reducing the need for external fertilization of crops). Examples of different combinations of these pathways are shown in Table 1, which includes a list of pathways and corresponding metabolic class (e.g., nitrogen mineralization). Several other potential combinations of pathways are shown in Table 1, where soil health indicator levels (SHI levels 1, 2, and 3) represent increasing granularity, which may be more representative of the underlying biology. Useful indicators balance granularity to be representative of agronomic outcomes with clear interpretability.
In the example of nitrogen mineralization, at SHI level 1 (shown in Table 1), the analytics system 100 considers each pathway as a factor of soil nitrogen supply. SHI level 2 includes more granularity, which separates nitrogen mineralization into organic nitrogen (Proteases and Chitinases pathways) and urea mineralization (urea mineralization and transporters). SHI level 3 includes even further granularity, retaining separate metabolism for soil protein (Proteases), chitin (Chitinases), and urea (urea mineralization, here disregarding urea transporters). In each of these instances (SHI levels 1, 2, 3), the analytics system 100 aggregates gene hit counts within the metabolic pathways belonging to the indicator (SHI level category), after being aggregated from gene counts to pathway counts. This aggregation from pathways to indicators may be accomplished by calculation of a function (e.g., including ratios, averages, weighted sums, etc.), with weights obtained by elimination, manual curation, statistical fitting, or machine learning. The analytics system 100 uses the aggregated counts to determine soil health indicators (e.g., including crop or regional specific percentiles). The analytics system 100 may provide for display (on a client device) soil health indicators for a given farm or series of fields.
E. Example Soil Health Indicator for Nitrogen Mineralization
Panel (a) shows that adding (e.g., normalized) gene count in pathways belonging to the nitrogen mineralization class does not provide a soil health indicator that tracks an expected increasing trend with soil C:N ratio. This ratio is thought to reflect soil and microbial nitrogen availability, and hence unmet demand for N driving energetically expensive nitrogen mineralization. Inspection of the underlying metabolic pathways reveals signals are mixed under this soil health indicator. For example, the class nitrogen mineralization includes (b) urea mineralization and (c) urea transporters. Panel (b) shows that the pathway for urea mineralization has an increasing relationship with soil C:N. On the other hand, panel (c) shows that urea transporters present no clear signal. As shown in panel (d), chitinase metabolism also has an inversely proportional relationship with soil C:N ratios. Panel (e) shows that secreted proteases have a proportional relationship with soil C:N ratios, as expected by first principles, with increasing C:N reflecting greater nitrogen scarcity.
The data in
Having eliminated urea transporters and chitin as less relevant, the analytics system 100 determines the soil health indicator shown in panel (f) as the sum of urea mineralization (ureases) and secreted proteases. These metabolic pathways included in the resulting soil health indicator are internally consistent because urea results as the product of nitrogen mineralization in proteins. In addition, these metabolic pathways have consistent signals within agronomic soils, unlike other pathways of nitrogen mineralization (as shown in
F. Nitrogen Cycling
Nitrogen cycling Soil Health Indicators calculated from functional gene data may represent various facets of the microbial nitrogen cycle in soil. These indicators may be broadly constructed to represent pathways of nitrogen loss, transformation, retention, or gain in agricultural soils. Nitrogen is a critical element for photosynthesis and plant biomass growth, and may be carefully managed in agricultural soils. Nitrogen cycling SHIs provide information on how added N fertilizer may be retained or lost through soil microbial processes, and may further indicate deficiencies or excess. This information can be used in agriculture to target fertilization rates to increase or maximize crop productivity while reducing or minimizing excess fertilizer use, which may increase costs and negatively impact water quality, for example.
Forms of nitrogen considered may include at least a part or the entire range of N oxidation states, including N2, N2O, NO, NO2−, NO3−, NH3+, NH4, and organic (carbon-bound) N compounds (e.g. amino acids/proteins, chitin, nitrile, etc). Metabolic processes driving transformation among these forms broadly include denitrification, nitrification, uptake of nitrate and ammonium, utilization of urea, chitin, and proteins, and atmospheric nitrogen fixation.
Simplified representations of these cycles for communicating agricultural SHIs may include soil N loss, retention, and supply; mobilization, immobilization, nitrate retention, urea and organic N mineralization, and N-fixation, etc. Examples of simplified representations of agricultural Soil Health Indicators from major N cycling pathways are given in Table 1. A detailed list of genes and their identifiers linking to sequence databases are given for the major pathway categories in Table 5 in the Appendix. For nitrogen cycling genes, these tables further list genes with ambiguous functions which may be tracked by the analytics system 100.
Calculation of these example indicator values may be made by summing total normalized gene counts within each biological or hierarchical SHI grouping. For example, representation of broad biological classes may be derived from the sum of gene counts for each gene in each biological class in Table 5. However, this level of detail may be too fine for some agricultural applications, and more condensed SHIs shown in Table 1 can be calculated, for example, as the sum of counts of biological groupings within each coarser level SHI class shown in Table 1. In this embodiment, counts of genes are aggregated by biological classes using Table 5, then further aggregated into simplified SHIs using Table 1. In various embodiments, the SHI levels are organized by granularity. In the example shown in Table 5, Level 3 is providing more biological detail, and Level 1 is more agriculturally interpretable (e.g., to a farmer). In addition, Level 2 is intermediate complexity grouping between Level 1 and Level 3. Multiple levels are shown to illustrate the process of deriving which level or combination of biology and agronomy is potentially most suitable. The analytics system 100 can use data driven support to combine related biological cycles without canceling or reducing signals carried by the underlying biology.
However, data transformations underlying SHI calculations are not limited to the hierarchies in these tables, and may further include calculation of ratios, averages, weighted sums, etc., including statistical fitting or machine learning. For example, indication of a process rate for denitrification may involve several factors represented by individual SHIs. This process may depend on genes SHIs for one or more of denitrification, nitrification, carbon availability, and oxygen availability. Machine learning (i.e., LASSO regression) may be used to fit the best combination of weights for genes or SHIs using a validation dataset.
G. Phosphorus Cycling
Phosphorus cycling Soil Health Indicators calculated from functional gene data may represent various processes of microbial phosphorus cycling in soil. These indicators may be broadly constructed to represent pathways of phosphorus transformations including microbial uptake, mineral P solubilization, or mineralization of various organic P forms in agricultural soils. Phosphorus is a critical element for plant biomass growth, and may be carefully managed in agricultural soils. Phosphorus cycling SHIs may indicate soil P deficiencies or P available to be mineralized by microbes for plant use. This information can be used agriculture to target fertilization rates to increase or maximize crop productivity while reducing or minimizing excess fertilizer use, which may increase costs and negatively impact water quality.
Forms of phosphorus considered may include free inorganic orthophosphate (PO42) or PO42 bound to minerals containing cations e.g., K+, Ca2+, Mn2+, zn2+, Fe3+, Al3+. Potential utilization of organic P forms like phytate, phosphonates, etc. may also be represented in gene based soil health indicators, as are various forms of high and low affinity phosphorus transporters (active and passive transport, respectively).
Simplified representations of these metabolic process for agricultural management may include as examples: counts or ratios of P uptake vs. P scavenging; P transporters vs. mineral P solubilization vs. organic P hydrolysis, etc. (Table 2). Each of these broad biologically based categories represents a collection of metabolic pathways, genes, and enzymes shown in detail in Table 6.
Calculation of these example indicator values may be made by summing total normalized gene counts within each biological or hierarchical SHI grouping. For example, broad biological classes may be derived from the sum of gene counts for each gene in each biological class in Table 6, which may be further aggregated into simplified SHIs using mapping shown in Table 2. However, data transformations underlying SHI calculations are not limited to the simple hierarchies in these tables and may further include calculation of ratios or weighted sums, averages, or ratios, including statistical data fitting by regression or machine learning.
H. Carbon Cycling
Carbon cycling Soil Health Indicators calculated from functional gene data may represent microbial potential to utilize various carbon substrates, for example, that are the primary constituents of soil organic matter. These SHIs are biological indicators of the relative availability of carbon for microbial decomposition, or inversely soil carbon resistance to breakdown and potential for its retention as soil organic matter. Carbon availability or retention may be of considerable interest for management of agriculture soils, as available C may drive microbial mediated processes e.g., nitrogen mineralization or microbial immobilization. Retention of soil carbon by the accumulation of less available forms (e.g., often utilized by microbes with slower metabolism) is of also of considerable agricultural interest as soil carbon accumulation can promote soil water retention along with microbial biomass that can help retain added nutrients and support plant beneficial microbes.
Carbon cycling SHIs represent pathways responsible for the breakdown of major soil carbon classes, including simple and complex polymer carbohydrates and polysaccharides, cellulose, hemicellulose, lignin, and various alkane and aromatic compounds. For the purposes of informing agricultural management practices, simplified indices like carbohydrate to aromatic gene ratios, or normalized counts of these groups or finer scale categories may be of interest.
Example formulations of simplified categories for communication of agricultural soil health indicators are shown in (Table 3). Each of the soil carbon compound classes represented therein may be degraded by enzymes which are encoded by numerous microbial genes (Table 7). These finer scale categories may also be of interest for informing predictions of other element cycles like nitrogen or phosphorus, given the potential overlap of groups of microbes specializing in certain transformations of both carbon and nutrients.
Calculation of these example indicator values may be made by summing total normalized gene counts within each biological or hierarchical SHI grouping. For example, broad biological classes may be derived from the sum of gene counts for each gene in each biological class in Table 7, which may be further aggregated into simplified SHIs using mapping shown in Table 3. However, data transformations underlying SHI calculations are not limited to the simple hierarchies in these tables, and may further include calculation of ratios or weighted sums, averages, or ratios, including statistical data fitting by regression or machine learning.
I. Oxygen Availability
Soil oxygen availability may be of agricultural interest given the role of higher or lower oxygen availability in cycling of nitrogen, phosphorus, or carbon, or waterlogging and anaerobic effects on plant root growth and rot. Microbial indicators of oxygen availability, for example, might reveal soils are poorly drained, with potential agricultural interventions including altered tillage, addition of mineral or organic materials to improve soil structure and drainage, or engineering solutions like drainage ditches, tile drainage, or altered to irrigation rates or timing.
Oxygen availability Soil Health Indicators calculated from functional gene data may represent the prevalence of oxygen in soil, as reflected by the presence of genetic elements involved in various microbial respiration pathways. These pathways may include aerobic respiration (including low and high affinity systems) and pathways of anaerobic metabolism including denitrification, iron reduction, sulfate reduction, and methanogenesis from organic substrates or H2/CO2 reduction.
Simplified representations of oxygen availability from the relative abundance of these pathways may include aerobic and anaerobic gene counts or ratios, further consideration of microaerobic or individual pathways, or even weighted sums of pathways, given the ability to multiply gene counts by quantitative scales based on the redox potential or free energy yield of these pathways, or scaled representations thereof (Table 4). Individual genes comprising each respiratory pathway in this SHI are shown in Table 8, compiled by translating genes into KEGG or Pfam identifiers used to retrieve counts from gene annotation tools.
Calculation of these example indicator values may be made by summing total normalized gene counts within each biological or hierarchical SHI grouping. However, data transformations underlying SHI calculations are not limited to the simple hierarchies in these tables, and may further include calculation of ratios or weighted sums, averages, or ratios, including data fitting machine learning.
The analytics system 100 receives 210 metadata describing a soil sample 275. In some embodiments, the metadata may indicate one or more crops grown in a geographical location having the soil sample 275. Example types of crop include corn, soybean, lettuce, strawberry, potato, among other types of fruits, nuts, vegetables, or plants. The cropping history of a geographical location may include a rotation of multiple types of crops, e.g., based on seasonality or the geographical location. In other embodiments, the metadata may indicate other information such as the geographical location, a current crop grown in the geographical location, or attributes describing treatment of the soil sample 275. The geographical location may be defined by global positioning system (GPS) coordinates or other suitable information, e.g., a neighborhood, city, state, country, or identification number. Example attributes describing treatment of the soil sample include agricultural techniques such as no-till farming, use of a cover crop to manage soil qualities (e.g., erosion, fertility, disease, or biodiversity), carbon farming, strip-till, and conservation agriculture. Attributes may also describe water or fertilizer usage, whether a crop is organic, temperature, precipitation, or climate, among other types of crop or soil related information. Metadata may also indicate a soil type of the soil sample 275. The analytics system 100 may process soil samples of different soil types, for example, sandy, silt, clay, loamy, and peat, among others. Additionally, the analytics system 100 may store soil chemical data including soil pH, carbon, and nutrients (e.g., nitrogen or nitrate, ammonium, phosphorus or phosphate, potassium, and calcium concentrations, among others).
The analytics system 100 determines 220 nucleic acid sequence reads of the soil sample 275. Referring now to
A soil sample 275 is obtained 222 using any of the methods previously described with reference to
In some embodiments, the nucleic acids of the microbial material may be processed prior to library preparation. For example, target genes or genome regions may be enriched for polymerase chain reaction (PCR) amplification or amplicon sequencing. Targeted DNA primers may be used to flank a region of interest. Alternatively, in shotgun sequencing, the microbial material may be prepared for sequencing of the entire content, e.g., microbes in a crop community of the processed soil sample. In some use cases, DNA fragment size may be controlled chemically using size selection gel beads, physically using ultrasonic shearing, or enzymatically using transposase fragmentation.
Sequencing library preparation is performed 224 on the extracted microbial material. Library preparation may include attaching sequencing adapters or tags to nucleic acids to facilitate reading of the nucleic acids. Sequencing tags may be unique to each sample (e.g., serving as a barcode) and enable identification of sequenced data associated with each sample in a multiplexed run with multiple samples. Libraries may also be prepared using other suitable methods such as ligation or transposase. In some use cases, library preparation includes protocols from sequencer original equipment manufacturers (OEMs), third party kit providers, or other resources. The analytics system 100 may store data from library preparation for future processing or analyses of other soil samples.
Once the sequencing library is prepared, the library or a portion of the library can be sequenced such that nucleic acid sequence reads of the microbial material are generated 225 using one or more techniques. In some embodiments, a sequencer 285 performs sequencing (e.g., of DNA or RNA) and outputs sequence reads of the microbial material. The sequencer 285 may provide the output sequence reads to the analytics system 100. The sequencer 285 can be communicatively coupled to the analytics system 100 through a wireless, wired, or a combination of wireless and wired communication technologies. In some embodiments, the nucleic acid sequence reads are generated using next generation sequencing (NGS) techniques including synthesis technology (ILLUMINA®), pyrosequencing (454 LIFE SCIENCES), ion semiconductor technology (Ion Torrent sequencing), single-molecule real-time sequencing (PACIFIC BIOSCIENCES®), or nanopore sequencing (OXFORD NANOPORE TECHNOLOGIES). DNA sequencing can also be performed as described in Sanger et al. (PNAS 74:5463, (1977)) and the Amersham International plc sequencing handbook, which methods are incorporated by reference herein. The analytics system 100 filters 226 the nucleic acid sequence reads, e.g., for quality control. In particular, the analytics system 100 may remove sequence reads having artificial multiplexing barcode or adapter sequences. In addition, the analytics system 100 may determine that a sequence read is low quality responsive to determining that a length of the sequence read is less than a threshold value, the sequence read includes at least a threshold number of ambiguous bases, or a read quality score (e.g., determined using a third-party tool) is less than a threshold score. The analytics system 100 may discard low quality sequence reads. The analytics system may also partition sequence reads using identification barcodes for demultiplexing batches of sequence reads generated from multiple samples.
Returning to
The analytics system 100 then determines 230 functional descriptors of the genes represented by the nucleic acid sequence reads. In various embodiments, a functional descriptor may be a functional gene measure in a soil sample. A functional descriptor may also be a metabolic pathway or SHI, or a collection of genes thereof. The analytics system 100 may store the functional descriptors in the gene database 106, e.g., in a table or another suitable type of data structure. In one embodiment, for each nucleic acid sequence read of at least a subset of the nucleic acid sequence reads, the analytics system 100 determines functional descriptors of the microbial gene associated with the nucleic acid sequence read.
In an embodiment, the analytics system 100 determines microbial functional genes present in the soil sample 275. For each of the present functional genes, the analytics system 100 determines a measure (e.g., count) of the functional gene in the soil sample 275. The analytics system 100 may normalize the measure using a total measure of functional genes in the soil sample, e.g., by dividing individual counts over a total count. In other embodiments, the analytics system 100 may normalize based on total sequence reads per sample, total annotations, total genome copies measured, etc. Normalization of the counts may account for sources of systematic variation or error in the counts, e.g., as stored in a table or another type of data structure. In other embodiments, instead of using nucleic acid sequencing, the analytics system 100 determines counts of organisms' functions using quantitative PCR (qPCR) or Droplet Digital PCR (ddPCR).
In some embodiments, data on the genetic functional potential of soil microbes is determined by comparing quality control filtered sequence reads from soil DNA to a reference database, e.g., using bioinformatics tools.
In some embodiments, the reference database contains genomic DNA sequence reads from known microorganisms, and information on functions (identity) of genes corresponding to these sequences. Examples of these databases may include NCBI's RefSeq database, the UniProt, Interpro, Pfam, SEED, COG, or KEGGs databases, etc. Databases may contain different systems of organizing genes (gene ontologies) into unique functional units (gene or protein identifiers). Genomic databases may also contain DNA sequence read data from multiple organisms for each gene or protein. Some databases may represent genes as consensus sequences or models (e.g., Pfam).
In some embodiments, identification of soil sample DNA sequence read queries by comparison to a reference database is made by a bioinformatics “pipeline.” These software pipelines are characterized by the ability execute multiple DNA read processing tasks in sequence, e.g., quality filtering and demultiplexing of raw sequencer data, determination of functional attributes, normalization and tabulation, etc. Such pipelines may contain public or proprietary software tools, with or without custom modifications including packaging to orchestrate handoffs between programs. Tools for identifying reads by matching the reference database (gene annotation) may use various methods including exact or fuzzy matching of DNA or more often translated protein sequences, or “k-mers” unique strings of sequences.
Examples of read annotation tools include BLAST, LAST, RapSearch2, DIAMOND, etc. Alternately, gene annotations may be determined using consensus probabilistic models (e.g., hidden Markov models) instead of alignment to genes of different taxa, e.g., using HMMer. Results from of these methods may be returned as a table representing the best matching gene annotation for each sequence read, along with a quality score. These read tables may be filtered to retain high confidence scoring gene hits, using metrics like bit score, e-value, match percent identity, etc.
In some embodiments, read tables are subsequently processed to obtain a table of counts for each gene present for each sample considered, and may be normalized. Gene count tables are compiled by aggregating read tables by their corresponding gene annotations. These gene count tables may be normalized by a variety of methods including total gene counts, sequence read depth or single copy maker genes in each sample, or other methods including centered log ratios.
In some embodiments, normalized gene count tables are used to obtain various Soil Health Indicators (SHIs). Genes counts used in the calculations for these SHIs are obtained for each sample from gene count tables, then raw SHIs values are calculated for each sample.
In some embodiments, raw SHI values are further scaled to their community percentile values, based on their crop communities. For example, each sample is assigned a crop community based on cropping, management, and geospatial (GPS) location. For each crop community, a database is compiled of raw SHI values for samples previously analyzed with correspondingly similar crop and location reported at the time of their collection. Each new sample may be given a community percentile value for each SHI or SHI category, based on the rank percentile of a given raw sample SHI among corresponding SHIs in the community database. Other data transformations may be applied prior to or following calculation of community percentile values.
The analytics system 100 determines 240 reference metrics of soil samples, e.g., from geographical locations or communities in which the one or more types of crop were grown. The reference metrics may include a distribution of values of soil health indicators retrieved from the soil health indicators database 102. Generally, the analytics system 100 may retrieve the reference metrics (or “crop community values”) from soil health indicators determined for soil samples of other users of the analytics system 100 or from other sources of reference information. For example, the analytics system 100 determines reference metrics from other soil samples within a threshold distance (e.g., 10, 50, 100, or 200 miles) from the soil sample. In a different example where the metadata indicates cropping history, the analytics system 100 determines reference metrics from other soil samples in which at least one common crop is currently or was previously grown. The analytics system 100 may use crop rotation systems as a basis for these community metrics, for example, corn-soybean rotations compared to continuous corn cultivation over consecutive years in a given geographical area. Furthermore, the analytics system 100 may determines reference metrics from other soil samples treated with similar or same agricultural techniques as those treated to the soil sample.
The analytics system 100 determines 250 a metric of the soil sample 275 using the functional descriptors and the reference metrics. The metric may be a measure of an agronomic attribute of the soil sample 275 or the soil at the geographical location generally. In an embodiment, the analytics system 100 determines a value of a soil health indicator using the functional descriptors. As previously described, the soil health indicator may be calculated as the function of measures of one or more types of microbial metabolic genes, e.g., associated with nitrogen or phosphorous cycling, utilization of different forms of carbon, or indicators of oxygen availability, plant growth promoting bacteria, root disease resistance, or post-harvest disease susceptibility. The analytics system 100 may determine the metric by performing one or more statistical transformations of the value of the soil health indicator. For example, the analytics system 100 determines a percentile of the value of the soil health indicator with respect to a distribution of soil health indicator values, as provided by the reference metrics. The percentiles may be scaled from 0 to 100%. In other embodiments, the analytics system 100 scales the value of soil health indicator to a different range such as 0.0 to 1.0 or 0 to 10, which may not necessarily be a percentile range.
In some embodiments, the analytics system 100 determines ranges of the reference metrics. The analytics system 100 may organize values of a soil health indicator for a set of fields (e.g., based on reference information of a community), within a threshold geographical location (or having another common characteristic or metadata), into buckets of a range of percentiles. For example, one bucket includes the top 10% of values of a soil health indicator associated with capacity of a given nutrient. Another bucket includes the next 10% of values of the soil health indicator, and so forth until a bucket including the bottom 10% of values of the soil health indicator. In other embodiments, the buckets may be associated with different intervals such as 20%, 25%, or 50%. When determining the metric for the soil sample, the analytics system 100 may identify a bucket to which the value of the soil health indicator of the soil sample belongs. For instance, the analytics system 100 determines that the value, of the soil health indicator of the soil sample collected from a geographical location, falls within the top 10% of values for nitrogen capacity of farms in the geographical location. Accordingly, the analytics system 100 may determine “0-10%” or “10%” as the metric.
In a different embodiment, the analytics system 100 may determine the metric according to standard deviations of the value of the soil health indicator away from an average value of the soil health indicator based on reference metrics. In some embodiments, the analytics system 100 may normalize the reference metrics to a logarithmic scale.
The analytics system 100 transmits 260 the metric to a client device 110 for display on a user interface, e.g., as shown in
In an optional step in some embodiments, soil at the geographical location (from which the soil sample 275 is obtained) is treated 270 according to the metric. For example, the metric may indicate that a crop has lower levels of nutrients in comparison to average levels for crops of the same or similar type, or crops grown in similar conditions or geographical locations. In response, a farmer may provide additional fertilizer, fumigation, water, cover crop, or other types of substances to the crop or soil to modify levels of oxygen, nitrogen, phosphorous, potassium, or carbon of the soil. In some embodiments, the analytics system 100 may receive new soil samples from a field after a treatment is applied to the field, e.g., according to metrics or recommendations provided by the analytics system 100. The analytics system 100 determines updated metrics (or recommendations) by processing the new soil samples and transmits the updated metrics to the client device 110 for presentation. Thus, the farmer may evaluate effect of the treatment by comparing the metrics before and after applying the treatment. The analytics system 100 may also receive additional soil samples from a field continuously over a period of time (e.g., weekly, monthly, or at arbitrary sample collection times) and track performance or health of the field by identifying trends in the determined metrics. The analytics system 100 may determine trends in context of crop community data.
In one embodiment, the analytics system 100 may provide a command to a client device 110 or another type of device to automatically treat the soil with a treatment loaded onto the device. For instance, the device is a manned or autonomous tractor for applying fertilizer, water, or other substance to soil or crops.
In an embodiment, the analytics system 100 determines 810 nucleic acid sequence reads of genetic material in a soil sample (e.g., soil sample 275 as shown in
The analytics system 100 determines 820, by processing the nucleic acid sequence reads (e.g., using a sequencer 285 as shown in
The analytics system 100 determines 850 a measure of an agronomic attribute of the soil at the geographical location as a function of the first, second, and third sets of measures. The agronomic attribute describes an aspect of the soil related to crop production. In one example, the agronomic attribute is nitrogen mineralization and the plurality of metabolic pathways includes at least secreted proteases and urea mineralization.
In other embodiments, the analytics system 100 does not necessarily use all three sets of the measures. For example, the analytics system 100 may determine the measure of the agronomic attribute using measures of gene functions and soil health indicators, without having to use measures of metabolic pathways. In some embodiments, the analytics system 100 may also filter out at least one metabolic pathway of the metabolic pathways. For instance, if the analytics system 100 determines that a given metabolic pathway is not indicative of (e.g., strongly correlated to) an expected agronomic attribute, the analytics system 100 determines not to use that given metabolic pathway in the determination of the measure of the agronomic attribute.
In an optional step 860, the soil at the geographical location may be treated according to the measure of the agronomic attribute. For example, an amount of fertilizer supplied to the soil may be increased or decreased based on the nitrogen mineralization metric for the soil. In some embodiments, the analytics system 100 can determine a recommendation for treatment of the soil using at least the measure of the agronomic attribute. The analytics system 100 may also transmit the measure of the agronomic attribute to a client device for display to a user. The analytics system 100 may calculate a community percentile metric of the measure of the agronomic attribute using reference metrics.
A. Nitrogen Cycling
Ammonification is the process of converting inorganic nitrate (NO3−) to ammonium (NH4+), which is free outside the cell. This process is also known as dissimilatory nitrate reduction or DNRA. In agricultural soils, this microbial process may help stem the loss of inorganic nitrogen by solubilization and leaching, as NO3− is soluble in soil solution whereas NH4+ is not (due to differences in polarity and ionic adsorption). Thus, ammonification may help stop the loss of nitrate by converting it to ammonium. Greater abundance of ammonification genes in these clay soils may reflect greater N retention.
Nitrification is the process of converting ammonium or ammonia (NH4+ or NH3) to nitrate (NO3−), which may be further divided by different bacterial groups carrying out conversions between NH4+/NH3→NO2− and then NO2−→NO3−. In agricultural settings, many N fertilizers may be applied as forms of NH4+ to reduce their mobility and leaching. Nitrifying bacteria can convert this NH4+/NH3 into the readily soluble NO3−. While both NH4+/NH3 and NO3— are available for uptake by plant roots, NO3— may be lost either to leaching or to the atmosphere as N2 as the product of denitrifying bacterial respiration. Thus, nitrification represents a conversion to forms of N that are potentially lost in agricultural soils, and the potential for this process may be greater in soils with less clay content.
B. Phosphorous Cycling
Soil P immobilization (uptake into microbial biomass, here inorganic P transporter genes) and Soil P scavenging (SHI level 2 includes mineral P solubilization, phosphatase and phosphonate lysis genes, described in Table 2) are two opposing processes. The first represents incorporation of inorganic P into microbial cells, while the second represents the potential to produce inorganic P from soil organic matter or mineral forms in the case of P solubilization.
However, the abundance of each of these gene categories may commonly reflect microbial effort allocated to acquire P, either by incorporating it or liberating it into solution. Thus, while metabolic functions differ between P solubilization and P scavenging, they may display common patterns along soil gradients.
Increasing microbial allocation to P acquisition via both P immobilization and P solubilization with soil pH (panel a) likely reflects pH driven changes in soil P availability. Soil pH affects the solubility and oxidation state of soil metals, which adsorb inorganic P (phosphate or PO4−). More acidic soils have more reduced metals including iron and aluminum, which in turn have less ability to adsorb inorganic P. Therefore, P availability is expected to be higher in acidic soils (lower pH), and lower in more basic soils (higher pH). The panel a) shows P acquisition gene SHIs (including P immobilization and P solubilization) increase with soil pH, and presumably therefore with P limitation.
P limitation may be assessed by comparing its relative availability, often when compared to N as indicated by N:P ratios. As N:P ratios increase, the relative availability of P decreases, reflecting declining P availability. Panel b) shows that microbial P acquisition gene abundances increase with increasing soil extractable N:P, a proxy for P limitation. Critical N:P ratios for plants and soils commonly are commonly divided as: N limited <N:P=16:1>P limited. However, soil extractable N:P (more direct measure of inorganic nutrient availability) tends to align with a different scale described by “ecoenzymatic stoichiometry” where the ratio of enzymes for transforming these elements is in balance is at N:P=1:1.
C. Carbon Cycling
As shown in panel a), community percentiles of SHI level 2 category “L2.0 Polymers” increase with soil clay content. The level 2 Polymers SHI category includes genes for the degradation of most plant biomass constituents (see Tables 1, 5). These polymers are known to bind to soil clay minerals, and higher gene abundances for their degradation reflects greater microbial allocation to access them in clay soils where they are more tightly bound.
In an opposite trend, genes grouped under the Nitroaromatics class (“C_arom_nitro”) decreased with soil clay content. These genes degrade nitroaromatic compounds, which could reflect scavenging of organic matter for nitrogen under limiting conditions. One interpretation of the decrease nitroaromatic degrading genes with increasing clay might be that low clay soils are more likely to be sandy, and sandy soils are known to be nitrogen poor. Higher abundances of nitroaromatic degrading genes in these soils would reflect increased nitrogen scavenging from aromatics in these lower nitrogen soils.
As shown in panel b), the community percentile abundance of genes degrading nitroaromatics decreased with increasing ratios of soil extractable nitrogen:phosphorus (Extractable N:P). As this ratio increases, nitrogen becomes increasingly available relative to phosphorus, such that higher numbers reflect decreasing nitrogen limitation. Decreasing nitroaromatic degrading genes with extractable N:P likely reflect lesser allocation to organic N scavenging from C compounds as N is more readily available at higher N:P. Extractable N:P ratios may better reflect the availability of inorganic forms of these elements when compared to total soil N:P.
Higher abundance of polymer degrading genes with increasing soil N:P may be in part attributable to the association of accumulated soil organic matter with both nitrogen content and ionic adsorption of ammonium.
D. Oxygen Availability
The foregoing description of the embodiments of the invention has been presented for the purpose of illustration; it is not intended to be exhaustive or to limit the invention to the precise forms disclosed. Persons skilled in the relevant art can appreciate that many modifications and variations are possible in light of the above disclosure.
Some portions of this description describe the embodiments of the invention in terms of algorithms and symbolic representations of operations on information. These algorithmic descriptions and representations are commonly used by those skilled in the data processing arts to convey the substance of their work effectively to others skilled in the art. These operations, while described functionally, computationally, or logically, are understood to be implemented by computer programs or equivalent electrical circuits, microcode, or the like. Furthermore, it has also proven convenient at times, to refer to these arrangements of operations as modules, without loss of generality. The described operations and their associated modules may be embodied in software, firmware, hardware, or any combinations thereof.
Any of the steps, operations, or processes described herein may be performed or implemented with one or more hardware or software modules, alone or in combination with other devices. In one embodiment, a software module is implemented with a computer program product including a computer-readable non-transitory medium containing computer program code, which can be executed by a computer processor for performing any or all of the steps, operations, or processes described.
Embodiments of the invention may also relate to a product that is produced by a computing process described herein. Such a product may include information resulting from a computing process, where the information is stored on a non-transitory, tangible computer readable storage medium and may include any embodiment of a computer program product or other data combination described herein.
As used herein any reference to “one embodiment” or “an embodiment” means that a particular element, feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment. The appearances of the phrase “in one embodiment” in various places in the specification are not necessarily all referring to the same embodiment.
As used herein, the terms “comprises,” “comprising,” “includes,” “including,” “has,” “having” or any other variation thereof, are intended to cover a non-exclusive inclusion. For example, a process, method, article, or apparatus that comprises a list of elements is not necessarily limited to only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Further, unless expressly stated to the contrary, “or” refers to an inclusive or and not to an exclusive or. For example, a condition A or B is satisfied by any one of the following: A is true (or present) and B is false (or not present), A is false (or not present) and B is true (or present), and both A and B are true (or present).
Finally, the language used in the specification has been principally selected for readability and instructional purposes, and it may not have been selected to delineate or circumscribe the inventive subject matter. It is therefore intended that the scope of the invention be limited not by this detailed description, but rather by any claims that issue on an application based hereon. Accordingly, the disclosure of the embodiments of the invention is intended to be illustrative, but not limiting, of the scope of the invention, which is set forth in the following claims.
A. Nitrogen Cycling
B. Phosphorus Cycling
C. Carbon Cycling
D. Oxygen Availability
This application claims the benefit of priority to U.S. Provisional Application No. 62/810,338, filed on Feb. 25, 2019, which is incorporated herein by reference in its entirety for all purposes.
| Number | Name | Date | Kind |
|---|---|---|---|
| 7255989 | Jeannin et al. | Aug 2007 | B1 |
| 20160232621 | Ethington | Aug 2016 | A1 |
| 20180148778 | Mercer | May 2018 | A1 |
| 20200024636 | Claypool et al. | Jan 2020 | A1 |
| Entry |
|---|
| Nielsen et al., Microorganisms as indicators of soil health, 2002, NERI Technical Report, 388, p. 1-85 (Year: 2002). |
| Vestergaard et al., Making big data smart—how to use metagenomics to understand soil quality, 2017, Biol Fertil Soils, 53, p. 479-A84 (Year: 2017). |
| Meneghine et al., Metagenomic analysis of soil and freshwater from zoo agricultural area with organic fertilization, 2017, PLOS one, p. 1-20 (Year: 2017). |
| Shokralla et al., 2012, Next-generation sequencing technologies for environmental DNA research, 2012, 21, p. 1794-1805 (Year: 2012). |
| Jeffries et al., Metagenomic Functional Potential Predicts Degradation Rates of a Model Organophosphorus Xenobiotic in Pesticide Contaminated Soils, 2018, frontiers in Microbiology, 9(147), p. 1-12 (Year: 2018). |
| An et al., Statistical Approach of Functional Profiling for a Microbial Community, 2014, PLOS ONE, 9(9), p. 1-11 (Year: 2014). |
| Srinvasulu, Enzymes and Pesticides, 2014, OMICS Group, p. 1-10 (Year: 2014). |
| Young et al. 2014, Limitations and recommendations for successful DNA extraction from forensic soil samples: A review, 2014, Science and Justice, 2014, p. 238-244 (Year: 2014). |
| Rinot et al., Soil health assessment: A critical review of current methodologies and a proposed new approach, 2019, Science of the Total Environment, p. 1484-1491 (Year: 2019). |
| PCT Invitation to Pay Additional Fees and, Where Applicable, Protest Fee, PCT Application No. PCT/US2020/19634, dated Apr. 14, 2020, 2 pages. |
| Kuske, C. R et al., “Diverse Uncultivated Bacterial Groups from Soils of the Arid Southwestern United States That are Present in Many Geographic Regions,” Applied and Environmental Microbiology, Sep. 1997, pp. 3614-3621, vol. 63, No. 9, [Online] [Retrieved on Aug. 5, 2020], Retrieved from the Internet<URL: https://aem.asm.org/content/aem/63/9/3614.full.pdf>. |
| PCT International Search Report and Written Opinion, PCT Application No. PCT/US2020/019634, dated Jun. 15, 2020, 14 pages. |
| Sanger, F. et al., “DNA sequencing with chain-terminating inhibitors,” PNAS, 1977, vol. 72, No. 12, pp. 5463-5467. |
| Number | Date | Country | |
|---|---|---|---|
| 20200271636 A1 | Aug 2020 | US |
| Number | Date | Country | |
|---|---|---|---|
| 62810338 | Feb 2019 | US |
